Quantum ESPRESSO: One Further Step toward the Exascale.

We review the status of the Quantum ESPRESSO software suite for electronic-structure calculations based on plane waves, pseudopotentials, and density-functional theory. We highlight the recent developments in the porting to GPUs of the main codes, using an approach based on OpenACC and CUDA Fortran offloading. We describe, in particular, the results achieved on linear-response codes, which are one of the distinctive features of the Quantum ESPRESSO suite. We also present extensive performance benchmarks on different GPU-accelerated architectures for the main codes of the suite.

Quantum ESPRESSO toward the exascale.

Quantum ESPRESSO is an open-source distribution of computer codes for quantum-mechanical materials modeling, based on density-functional theory, pseudopotentials, and plane waves, and renowned for its performance on a wide range of hardware architectures, from laptops to massively parallel computers, as well as for the breadth of its applications. In this paper, we present a motivation and brief review of the ongoing effort to port Quantum ESPRESSO onto heterogeneous architectures based on hardware accelerators, which will overcome the energy constraints that are currently hindering the way toward exascale computing.

Pd doping, conformational, and charge effects on the dichroic response of a monolayer protected Au38(SR)24 nanocluster.

TDDFT simulations of the absorption and CD spectra of a Pd2Au36(SC2H4Ph)24 monolayer-protected cluster (MPC) are carried out with the aim of investigating the effects of doping, conformational degrees of freedom of the thiolates' end-groups, and charge states on the optical and dichroic response of a prototypical MPC species. Clear signatures of Pd doping in both absorption and CD spectra are found to be a consequence of the participation of Pd (4d) states in the ligand-based d-band and on the unoccupied MOs of lower energy. Exploration of conformational space points to a much greater sensitivity of optical rotation to the conformation of the end-groups of the organic monolayer compared to absorption. Finally, the effect of charge is mainly seen as a decreased dependence of the dichroic response on conformation. The agreement between the TDDFT predictions and the available experimental data is good, and enables an assignment of absorption and CD bands to specific classes of one-particle excitations.

Individual Component Map of Rotatory Strength and Rotatory Strength Density Plots As Analysis Tools of Circular Dichroism Spectra of Complex Systems.

We introduce individual component maps of rotatory strength (ICM-RS) and rotatory strength density (RSD) plots as analysis tools of chiro-optical linear response spectra deriving from time-dependent density functional theory (TDDFT) simulations. ICM-RS and RSD allow one to visualize the origin of chiro-optical response in momentum or real space, including signed contributions and therefore highlighting cancellation terms that are ubiquitous in chirality phenomena, and should be especially useful in analyzing the spectra of complex systems. As test cases, we use ICM-RS and RSD to analyze circular dichroism spectra of selected (Ag-Au)30(SR)18 monolayer-protected metal nanoclusters, showing the potential of the proposed tools to derive insight and understanding, and eventually rational design, in chiro-optical studies of complex systems.

Au38(SPh)24: Au38 Protected with Aromatic Thiolate Ligands.

Au38(SR)24 is one of the most extensively investigated gold nanomolecules along with Au25(SR)18 and Au144(SR)60. However, so far it has only been prepared using aliphatic-like ligands, where R = -SC6H13, -SC12H25 and -SCH2CH2Ph. Au38(SCH2CH2Ph)24 when reacted with HSPh undergoes core-size conversion to Au36(SPh)24, and existing literature suggests that Au38(SPh)24 cannot be synthesized. Here, contrary to prevailing knowledge, we demonstrate that Au38(SPh)24 can be prepared if the ligand exchanged conditions are optimized, under delicate conditions, without any formation of Au36(SPh)24. Conclusive evidence is presented in the form of matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS), electrospray ionization mass spectra (ESI-MS) characterization, and optical spectra of Au38(SPh)24 in a solid glass form showing distinct differences from that of Au38(S-aliphatic)24. Theoretical analysis confirms experimental assignment of the optical spectrum and shows that the stability of Au38(SPh)24 is not negligible with respect to that of its aliphatic analogous, and contains a significant component of ligand-ligand attractive interactions. Thus, while Au38(SPh)24 is stable at RT, it converts to Au36(SPh)24 either on prolonged etching (longer than 2 hours) at RT or when etched at 80 °C.

A new time dependent density functional algorithm for large systems and plasmons in metal clusters.

A new algorithm to solve the Time Dependent Density Functional Theory (TDDFT) equations in the space of the density fitting auxiliary basis set has been developed and implemented. The method extracts the spectrum from the imaginary part of the polarizability at any given photon energy, avoiding the bottleneck of Davidson diagonalization. The original idea which made the present scheme very efficient consists in the simplification of the double sum over occupied-virtual pairs in the definition of the dielectric susceptibility, allowing an easy calculation of such matrix as a linear combination of constant matrices with photon energy dependent coefficients. The method has been applied to very different systems in nature and size (from H2 to [Au147](-)). In all cases, the maximum deviations found for the excitation energies with respect to the Amsterdam density functional code are below 0.2 eV. The new algorithm has the merit not only to calculate the spectrum at whichever photon energy but also to allow a deep analysis of the results, in terms of transition contribution maps, Jacob plasmon scaling factor, and induced density analysis, which have been all implemented.

Vibrationally resolved high-resolution NEXAFS and XPS spectra of phenanthrene and coronene.

We performed a combined experimental and theoretical study of the C1s Near-Edge X-ray Absorption Fine-Structure (NEXAFS) spectroscopy and X-ray Photoelectron Spectroscopy in the gas phase of two polycyclic aromatic hydrocarbons (phenanthrene and coronene), typically formed in combustion reactions. In the NEXAFS of both molecules, a double-peak structure appears in the C1s → LUMO region, which differ by less than 1 eV in transition energies. The vibronic coupling is found to play an important role in such systems. It leads to weakening of the lower-energy peak and strengthening of the higher-energy one because the 0 - n (n > 0) vibrational progressions of the lower-energy peak appear in nearly the same region of the higher-energy peak. Vibrationally resolved theoretical spectra computed within the Frank-Condon (FC) approximation and linear coupling model agree well with the high-resolution experimental results. We find that FC-active normal modes all correspond to in-plane vibrations.