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.

Measuring Shared Electrons in Extended Molecular Systems: Covalent Bonds from Plane-Wave Representation of Wave Function.

In the study of materials and macromolecules by first-principle methods, the bond order is a useful tool to represent molecules, bulk materials and interfaces in terms of simple chemical concepts. Despite the availability of several methods to compute the bond order, most applications have been limited to small systems because a high spatial resolution of the wave function and an all-electron representation of the electron density are typically required. Both limitations are critical for large-scale atomistic calculations, even within approximate density-functional theory (DFT) approaches. In this work, we describe our methodology to quickly compute delocalization indices for all atomic pairs, while keeping the same representation of the wave function used in most compute-intensive DFT calculations on high-performance computing equipment. We describe our implementation into a post-processing tool, designed to work with Quantum ESPRESSO, a popular open-source DFT package. In this way, we recover a description in terms of covalent bonds from a representation of wave function containing no explicit information about atomic types and positions.

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.

Investigation of structural, electronic and magnetic properties of breathing metal-organic framework MIL-47(Mn): a first principles approach.

The structural, electronic and magnetic properties of the MIL-47(Mn) metal-organic framework are investigated using first principles calculations. We find that the large-pore structure is the ground state of this material. We show that upon transition from the large-pore to the narrow-pore structure, the magnetic ground-state configuration changes from antiferromagnetic to ferromagnetic, consistent with the computed values of the intra-chain coupling constant. Furthermore, the antiferromagnetic and ferromagnetic configuration phases have intrinsically different electronic behavior: the former is semiconducting, the latter is a metal or half-metal. The change of electronic properties during breathing posits MIL-47(Mn) as a good candidate for sensing and other applications. Our calculated electronic band structure for MIL-47(Mn) presents a combination of flat dispersionless and strongly dispersive regions in the valence and conduction bands, indicative of quasi-1D electronic behavior. The spin coupling constants are obtained by mapping the total energies onto a spin Hamiltonian. The inter-chain coupling is found to be at least one order of magnitude smaller than the intra-chain coupling for both large and narrow pores. Interestingly, the intra-chain coupling changes sign and becomes five times stronger going from the large pore to the narrow pore structure. As such MIL-47(Mn) could provide unique opportunities for tunable low-dimensional magnetism in transition metal oxide systems.

Quantum Crystallography: Current Developments and Future Perspectives.

Crystallography and quantum mechanics have always been tightly connected because reliable quantum mechanical models are needed to determine crystal structures. Due to this natural synergy, nowadays accurate distributions of electrons in space can be obtained from diffraction and scattering experiments. In the original definition of quantum crystallography (QCr) given by Massa, Karle and Huang, direct extraction of wavefunctions or density matrices from measured intensities of reflections or, conversely, ad hoc quantum mechanical calculations to enhance the accuracy of the crystallographic refinement are implicated. Nevertheless, many other active and emerging research areas involving quantum mechanics and scattering experiments are not covered by the original definition although they enable to observe and explain quantum phenomena as accurately and successfully as the original strategies. Therefore, we give an overview over current research that is related to a broader notion of QCr, and discuss options how QCr can evolve to become a complete and independent domain of natural sciences. The goal of this paper is to initiate discussions around QCr, but not to find a final definition of the field.

Reproducibility in density functional theory calculations of solids.

The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Reaction pathways for oxygen evolution promoted by cobalt catalyst.

The in-depth understanding of the molecular mechanisms regulating the water oxidation catalysis is of key relevance for the rationalization and the design of efficient oxygen evolution catalysts based on earth-abundant transition metals. Performing ab initio DFT+U molecular dynamics calculations of cluster models in explicit water solution, we provide insight into the pathways for oxygen evolution of a cobalt-based catalyst (CoCat). The fast motion of protons at the CoCat/water interface and the occurrence of cubane-like Co-oxo units at the catalyst boundaries are the keys to unlock the fast formation of O-O bonds. Along the resulting pathways, we identified the formation of Co(IV)-oxyl species as the driving ingredient for the activation of the catalytic mechanism, followed by their geminal coupling with O atoms coordinated by the same Co. Concurrent nucleophilic attack of water molecules coming directly from the water solution is discouraged by high activation barriers. The achieved results suggest also interesting similarities between the CoCat and the Mn4Ca-oxo oxygen evolving complex of photosystem II.

The interactions of nitrogen dioxide with graphene-stabilized Rh clusters: a DFT study.

We study the interactions of NO2 gas molecules with Rh nanoparticles supported on graphene, using first-principles molecular dynamics in the Car-Parrinello scheme. The stability, morphology, adsorption energies of various models of Rhx nanoparticles (x = 1, 3, 10, 20) supported on graphene, and the binding of NO2 molecules to the Rh clusters, together with its effect on the graphene properties, are reported. Metastable flat structures anchored to the substrate that can bind NO2 to Rh via both N and O atoms are identified, with adsorption energies in the range 60-70 kcal per mole per molecule.

Zn induced structural aggregation patterns of ß-amyloid peptides by first-principle simulations and XAS measurements.

We show in this paper that in the presence of Zn ions a peculiar structural aggregation pattern of ß-amyloid peptides in which metal ions are sequentially coordinated to either three or four histidines of nearby peptides is favored. To stabilize this configuration a deprotonated imidazole ring from one of the histidines forms a bridge connecting two adjacent Zn ions. Though present in zeolite imidazolate frameworks, remarkably in biological compounds this peculiar Zn-imidazolate-Zn topology is only found in enzymes belonging to the Cu,Zn-superoxide dismutase family in the form of an imidazolate bridging Cu and Zn. The results we present are obtained by combining X-ray absorption spectroscopy experimental data with detailed first-principle molecular dynamics simulations.

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials.

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

Calculation of near-edge x-ray-absorption fine structure at finite temperatures: spectral signatures of hydrogen bond breaking in liquid water.

We calculate the near-edge x-ray-absorption fine structure of H(2)O in the gas, hexagonal ice, and liquid phases using heuristic density-functional based methods. We present a detailed comparison of our results with experiment. The differences between the ice and water spectra can be rationalized in terms of the breaking of hydrogen bonds around the absorbing molecule. In particular the increase in the pre-edge absorption feature from ice to water is shown to be due to the breaking of a donor hydrogen bond. We also find that in water approximately 19% of hydrogen bonds are broken.