Q uantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Q uantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice Ih at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner. Significance Water is vital to our everyday life, but its structure at a molecular level is still not fully understood from either experiment or theory. The latter is hampered by our inability to construct a purely predictive, first principles model. The difficulty in modeling water lies in capturing the delicate interplay among the many strong and weak forces that govern its behavior and phase diagram. Herein, molecular simulations with a recently proposed nonempirical quantum mechanical approach (the SCAN density functional) yield an excellent description of the structural, electronic, and dynamic properties of liquid water. SCAN (strongly constrained and appropriately normed)-based approaches, which describe diverse types of bonds in materials on an equal, accurate footing, will likely enable efficient and reliable modeling of aqueous phase chemistry.

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H₂O)₁₂₈ significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, S_OO(Q), and corresponding oxygen-oxygen radial distribution function, g_OO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, P_OOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.

The adsorption of aromatic molecules on metal surfaces plays a key role in condensed matter physics and functional materials. Depending on the strength of the interaction between the molecule and the surface, the binding is typically classified as either physisorption or chemisorption. Van der Waals (vdW) interactions contribute significantly to the binding in physisorbed systems, but the role of the vdW energy in chemisorbed systems remains unclear. Here we study the interaction of benzene with the (111) surface of transition metals, ranging from weak adsorption (Ag and Au) to strong adsorption (Pt, Pd, Ir, and Rh). When vdW interactions are accurately accounted for, the barrier to adsorption predicted by standard density-functional theory (DFT) calculations essentially vanishes, producing a metastable precursor state on Pt and Ir surfaces. Notably, vdW forces contribute more to the binding of covalently bonded benzene than they do when benzene is physisorbed. Comparison to experimental data demonstrates that some of the recently developed methods for including vdW interactions in DFT allow quantitative treatment of both weakly and strongly adsorbed aromatic molecules on metal surfaces, extending the already excellent performance found for molecules in the gas phase.

The first principles methods, density-functional theory and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen bonding in ambient and high-pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from hydrogen bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from density-functional theory exchange-correlation functionals which neglect vdW forces are greatly overestimated.