Search

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable d...

A comprehensive microscopic understanding of ambient liquid water is a major challenge for $ab$ $initio$ simulations as it simultaneously requires an accurate quantum mechanical description of the underlying potential energy surface (PES) as well as extensive sampling of configuration space. Due to the presence of light atoms (e.g., H or D), nuclear quantum fluctuations lead to observable changes in the structural properties of liquid water (e.g., isotope effects), and therefore provide yet another challenge for $ab$ $initio$ approaches. In this work, we demonstrate that the combination of dispersion-inclusive hybrid density functional theory (DFT), the Feynman discretized path-integral (PI) approach, and machine learning (ML) constitutes a versatile $ab$ $initio$ based framework that enables extensive sampling of both thermal and nuclear quantum fluctuations on a quite accurate underlying PES. In particular, we employ the recently developed deep potential molecular dynamics (DPMD) model---a neural-network representation of the $ab$ $initio$ PES---in conjunction with a PI approach based on the generalized Langevin equation (PIGLET) to investigate how isotope effects influence the structural properties of ambient liquid H$_2$O and D$_2$O. Through a detailed analysis of the interference differential cross sections as well as several radial and angular distribution functions, we demonstrate that this approach can furnish a semi-quantitative prediction of these subtle isotope effects., Comment: 19 pages, 5 figures, and 1 table

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553-823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

We carry out X-ray absorption spectroscopy experiment at oxygen K-edge in croconic acid (C5H2O5) crystal as a prototype of ferroelectric organic molecular solid, whose electric polarization is generated by proton transfer. The experimental spectrum is well reproduced by the electron-hole excitation theory simulations from configuration generated by ab initio molecular dynamics simulation. When inversion symmetry is broken in ferroelectric state, the hydrogen bonding environment on the two bonded molecules become inequivalent. Such a difference is sensitively probed by the bound excitation in the pre-edge, which are strongly localized on the excited molecules. Our analysis shows that a satellite peak in the pre-edge will emerge at higher excitation energy which serves as a clear signature of ferroelectricity in the material., Comment: 6 pages, 3 figures

We explore how anharmonicity, nuclear quantum effects (NQE), many-body dispersion interactions, and Pauli repulsion influence thermal properties of dispersion-bound molecular crystals. Accounting for anharmonicity with ab initio molecular dynamics yields cell parameters accurate to within $2%$ of experiment for a set of pyridinelike molecular crystals at finite temperatures and pressures. From the experimental thermal expansion curve, we find that pyridine-I has a Debye temperature just above its melting point, indicating sizable NQE across the entire crystalline range of stability. We find that NQE lead to a substantial volume increase in pyridine-I $(\ensuremath{\approx}40$% more than classical thermal expansion at 153 K) and attribute this to intermolecular Pauli repulsion promoted by intramolecular quantum fluctuations. When predicting delicate properties such as the thermal expansivity, we show that many-body dispersion interactions and more sophisticated density functional approximations improve the accuracy of the theoretical model.

We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the strongly constrained and appropriately normed (SCAN) meta-generalized-gradient approximation (GGA) functional approximation and carry out systematic comparisons with the results obtained from the GGA-level Perdew-Burke-Ernzerhof (PBE) functional and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.

Quantum 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 pseudo-potential and projector-augmented-wave approaches. Quantum 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 theirs 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., Comment: Psi-k highlight September 2017: psi-k.net/dowlnload/highlights/Highlight_137.pdf; J. Phys.: Condens. Matter, accepted

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.

Biological membranes are essential for cell life and hydration. Water provides the driving force for the assembly and stability of many cell components. Here, we study the structural properties of water in a phospholipid membrane. We characterize the local structures, inspecting the intermediate range order (IRO) and adopting a sensitive local order metric recently proposed by Martelli et al. that measures and grades the degree of overlap of the local environment with the structures of perfect ice. Close to the membrane, water acquires a high IRO and changes its dynamical properties; i.e., its translational and rotational degrees of freedom slow in a region that extends over ≃ 1 nm from the membrane interface. Surprisingly, we show that at distances as far as ≃ 2.5 nm from the interface, although the bulk-like dynamics are recovered, the IRO of water is still slightly higher than that in the bulk under the same thermodynamic conditions. Therefore, the water-membrane interface has a structural effect at ambient conditions that propagates further than the often-invoked 1-nm length scale. Consequently, this should be considered when analyzing experimental data of water confined by membranes and could help us to understand the role of water in biological systems.