Structure and Dynamics of Liquid Water from ab Initio Molecular Dynamics—Comparison of BLYP, PBE, and revPBE Density Functionals with and without van der Waals Corrections

Grid-Based Projector Augmented Wave (GPAW) Implementation of Quantum Mechanics/Molecular Mechanics (QM/MM) Electrostatic Embedding and Application to a Solvated Diplatinum Complex

A multiscale density functional theory-quantum mechanics/molecular mechanics (DFT-QM/MM) scheme is presented, based on an efficient electrostatic coupling between the electronic density obtained from a grid-based projector augmented wave (GPAW) implementation of density functional theory and a classical potential energy function. The scheme is implemented in a general fashion and can be used with various choices for the descriptions of the QM or MM regions. Tests on H₂O clusters, ranging from dimer to decamer show that no systematic energy errors are introduced by the coupling that exceeds the differences in the QM and MM descriptions. Over 1 ns of liquid water, Born-Oppenheimer QM/MM molecular dynamics (MD) are sampled combining 10 parallel simulations, showing consistent liquid water structure over the QM/MM border. The method is applied in extensive parallel MD simulations of an aqueous solution of the diplatinum [Pt₂(P₂O₅H₂)₄]^4- complex (PtPOP), spanning a total time period of roughly half a nanosecond. An average Pt-Pt distance deviating only 0.01 Å from experimental results, and a ground-state Pt-Pt oscillation frequency deviating by <2% from experimental results were obtained. The simulations highlight a remarkable harmonicity of the Pt-Pt oscillation, while also showing clear signs of Pt-H hydrogen bonding and directional coordination of water molecules along the Pt-Pt axis of the complex.

Protonation-modulated localization of excess electrons in histidine aqueous solutions revealed by ab initio molecular dynamics simulations: anion-centered versus cation-centered localization

In this work, we present an ab initio molecular dynamics simulation study on the interaction of an excess electron (EE) with histidine in its aqueous solution. Two different configurations of histidine (imidazole group protonated or not) are considered to reflect its different existing forms in neutral or slightly acidic surroundings. The simulation results indicate that localizations of EEs in different aqueous histidine solutions are quite different and are strongly affected by protonation of the side chain imidazole group and are thus pH-controlled. In neutral aqueous histidine solution, an EE localizes onto the carboxyl anionic group of the amino acid backbone after a relatively lengthy diffuse state, performing just like in an aliphatic amino acid solution. But in weakly acidic solution in which the side chain imidazole group is protonated, an EE undergoes a short lifetime diffuse state and finally localizes on the protonated imidazole group. We carefully examine these two different localization dynamics processes and analyze the competition between different dominating groups in their corresponding electron localization mechanisms. To explain the difference, we investigate the frontier molecular orbitals of these two systems and find that their energy levels and compositions are important to determine these differences. These findings can provide helpful information to understand the interaction mechanisms of low energy EEs with amino acids and even oligopeptides, especially with aromatic rings.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis

The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH₃, NH₄⁺, CN^-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

Ionic diffusion and proton transfer in aqueous solutions of alkali metal salts

We report on a series of ab initio molecular dynamics investigations on LiCl, NaCl, and KCl aqueous solutions under the effect of static electric fields.

Ab initio Electronic Structure of Liquid Water

Self-consistent GW calculations with efficient vertex corrections are employed to determine the electronic structure of liquid water. Nuclear quantum effects are taken into account through ab initio path-integral molecular dynamics simulations. We reveal a sizable band-gap renormalization of up to 0.7 eV due to hydrogen-bond quantum fluctuations. Our calculations lead to a band gap of 8.9 eV, in accord with the experimental estimate. We further resolve the ambiguities in the band-edge positions of liquid water. The valence-band maximum and the conduction-band minimum are found at -9.4 and -0.5  eV with respect to the vacuum level, respectively.

Water adsorption on bimetallic PtRu/Pt(111) surface alloys

The adsorption of water on bimetallic PtRu/Pt(111) surface alloys has been studied based on periodic density functional theory calculations including dispersion corrections. The Ru atoms of the PtRu surface alloy interact more strongly with water than Pt atoms, as far as both single water molecules and ice-like hexagonal structures are concerned. Within the surface alloy layer, the lateral ligand effect reducing the local reactivity of the surface atoms with increasing Ru content is more dominant than the opposing geometric effect due to the tensile strain. The structural preference for the Ru atoms also prevails at room temperature, as ab initio molecular dynamics simulations show.

Ab initio molecular dynamics study of Se(iv) species in aqueous environment

An ab initio molecular dynamics investigation is carried out on various water-borne Se(iv) species, H₂SeO₃, HSeO₃^- and SeO₃^2-, in aqueous environment. Consistent with the reported acid dissociation constants, in neutral solution H₂SeO₃ exchanges protons with the surrounding water molecules establishing a dynamic equilibrium with HSeO₃^-. The SeO₃^2- species is found to be stable only in basic environment, which is emulated in the present simulation through introducing a hydroxide ion, OH^-, in the system. The hydration structure, hydrogen bonding and spectroscopic signatures of the species are comprehensively analyzed. The influence of the solute's hydration structure on the structural and dynamic response of the solvent is discussed. The correlation between the strength as well as the number of hydrogen bonds accepted by the solute on its vibrational properties are analyzed.

Reactive wetting properties of TiO2 nanoparticles predicted by ab initio molecular dynamics simulations

Small-sized wet TiO2 nanoparticles have been investigated by ab initio molecular dynamics simulations. Chemical and physical adsorption of water on the TiO2-water interface was studied as a function of water content, ranging from dry nanoparticles to wet nanoparticles with monolayer coverage of water. The surface reactivity was shown to be a concave function of water content and driven by surface defects. The local coordination number at the defect was identified as the key factor to decide whether water adsorption proceeds through dissociation or physisorption on the surface. A consistent picture of TiO2 nanoparticle wetting at the microscopic level emerges, which corroborates existing experimental data and gives further insight into the molecular mechanisms behind nanoparticle wetting. These calculations will facilitate the engineering of metal oxide nanoparticles with a controlled catalytic water activity.

Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration

A challenge in density functional theory is developing functionals that simultaneously describe intermolecular electron correlation and electron delocalization. Recent exchange-correlation functionals address those two issues by adding corrections important at long ranges: an atom-centered pairwise dispersion term to account for correlation and a modified long-range component of the electron exchange term to correct for delocalization. Here we investigate how those corrections influence the accuracy of binding free energy predictions for sodium-water clusters. We find that the dual-corrected ωB97X-D functional gives cluster binding energies closest to high-level ab initio methods (CCSD(T)). Binding energy decomposition shows that the ωB97X-D functional predicts the smallest ion-water (pairwise) interaction energy and larger multibody contributions for a four-water cluster than most other functionals - a trend consistent with CCSD(T) results. Also, ωB97X-D produces the smallest amounts of charge transfer and the least polarizable waters of the density functionals studied, which mimics the lower polarizability of CCSD. When compared with experimental binding free energies, however, the exchange-corrected CAM-B3LYP functional performs best (error <1 kcal/mol), possibly because of its parametrization to experimental formation enthalpies. For clusters containing more than four waters, "split-shell" coordination must be considered to obtain accurate free energies in comparison with experiment.

Adsorption-Induced Surface Stresses of the Water/Quartz Interface: Ab Initio Molecular Dynamics Study

Adsorption-induced deformation is expansion or contraction of a solid due to adsorption on its surface. This phenomenon is important for a wide range of applications, from chemomechanical sensors to methane recovery from geological formations. The strain of the solid is driven by the change of the surface stress due to adsorption. Using ab initio molecular dynamics, we calculate the surface stresses for the dry α-quartz surfaces, and investigate how these stresses change when the surfaces are exposed to water. We find that the nonhydroxylated surface shows small and approximately isotropic changes in stress, while the hydroxylated surface, which interacts more strongly with the polar water molecules, shows larger and qualitatively anisotropic (opposite sign in xx and yy) surface stress changes. All of these changes are several times larger than the surface tension of water itself. The anisotropy and possibility of positive surface stress change can explain experimentally observed surface area contraction due to adsorption.

First-principles study of the structure of water layers on flat and stepped Pb electrodes

On the basis of perodic density functional theory (DFT) calculations, we have addressed the geometric structures and electronic properties of water layers on flat and stepped Pb surfaces. In contrast to late d-band metals, on Pb(111) the energy minimum structure does not correspond to an ice-like hexagonal arrangement at a coverage of 2/3, but rather to a distorted structure at a coverage of 1 due to the larger lattice constant of Pb. At stepped Pb surfaces, the water layers are pinned at the step edge and form a complex network consisting of rectangles, pentagons and hexagons. The thermal stability of the water layers has been studied by using ab initio molecular dynamics simulations (AIMD) at a temperature of 140 K. Whereas the water layer on Pb(111) is already unstable at this temperature, the water layers on Pb(100), Pb(311), Pb(511) and Pb(711) exhibit a higher stability because of stronger water-water interactions. The vibrational spectra of the water layers at the stepped surfaces show a characteristic splitting into three modes in the O-H stretch region.

Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics

We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.

Reaction Mechanism for Direct Proton Transfer from Carbonic Acid to a Strong Base in Aqueous Solution I: Acid and Base Coordinate and Charge Dynamics

Protonation by carbonic acid H2CO3 of the strong base methylamine CH3NH2 in a neutral contact pair in aqueous solution is followed via Car-Parrinello molecular dynamics simulations. Proton transfer (PT) occurs to form an aqueous solvent-stabilized contact ion pair within 100 fs, a fast time scale associated with the compression of the acid-base hydrogen-bond (H-bond), a key reaction coordinate. This rapid barrierless PT is consistent with the carbonic acid-protonated base pKa difference that considerably favors the PT, and supports the view of intact carbonic acid as potentially important proton donor in assorted biological and environmental contexts. The charge redistribution within the H-bonded complex during PT supports a Mulliken picture of charge transfer from the nitrogen base to carbonic acid without altering the transferring hydrogen's charge from approximately midway between that of a hydrogen atom and that of a proton.

Influence of Solvation and Dynamics on the Mechanism and Kinetics of Nucleophilic Aromatic Substitution Reactions in Liquid Ammonia

The role of the solvent and the influence of dynamics on the kinetics and mechanism of the SNAr reaction of several halonitrobenzenes in liquid ammonia, using both static calculations and dynamic ab initio molecular dynamics simulations, are investigated. A combination of metadynamics and committor analysis methods reveals how this reaction can change from a concerted, one-step mechanism in gas phase to a stepwise pathway, involving a metastable Meisenheimer complex, in liquid ammonia. This clearly establishes, among others, the important role of the solvent and highlights the fact that accurately treating solvation is of crucial importance to correctly unravel the reaction mechanism. It is indeed shown that H-bond formation of the reacting NH3 with the solvent drastically reduces the barrier of NH3 addition. The halide elimination step, however, is greatly facilitated by proton transfer from the reacting NH3 to the solvent. Furthermore, the free energy surface strongly depends on the halide substituent and the number of electron-withdrawing nitro substituents.

Anisotropic structure and dynamics of the solvation shell of a benzene solute in liquid water from ab initio molecular dynamics simulations

The anisotropic structure and dynamics of the hydration shell of a benzene solute in liquid water have been investigated by means of ab initio molecular dynamics simulations using the BLYP (Becke–Lee–Yang–Parr) and dispersion corrected BLYP-D functionals.

Demulsification of heavy oil-in-water emulsions by reduced graphene oxide nanosheets

The purpose of this work is to investigate the relationship between the chemical structure of the GO materials and the demulsification performance.

Ab initio molecular dynamics study of an aqueous NaCl solution under an electric field

Ab initio molecular dynamics simulations of salty water under an electric field reveal two regimes of the relative mobilities of chlorine and sodium ions. When water dissociation and proton transfer are actived at strong field intensities, the presence of the ions hinders the efficiency of the proton transfer mechanism.

Neural network molecular dynamics simulations of solid–liquid interfaces: water at low-index copper surfaces

Molecular dynamics simulation of the water–copper interface have been carried out using high-dimensional neural network potential based on density functional theory.

Probing the structural and dynamical properties of liquid water with models including non-local electron correlation

Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.

Hydrophobicity: effect of density and order on water's rotational slowing down

Ab initio molecular dynamics (AIMD) simulations of over 4.5 ns were performed in the temperature range of T = 260-350 K with van der Waals corrections to investigate the relationship between local water density and tetrahedral order in bulk water and in the presence of a hydrophobe, tetramethylurea (TMU). We demonstrate that in bulk water, defects consisting of 5- and higher coordinated water are a major contributor to dynamics. Close to a hydrophobe, 3-coordinated defects take over. The co-existence of these defects gives rise to very different local densities. We propose that the slowing down of rotational motion close to a hydrophobe is induced by an interplay between density and order with the slowing down decreasing in the following order: (i) low-density ordered-water, (ii) normal-density ordered-water, (iii) high-density ordered-water and (iv) disordered-water. The proportions of these water environments vary with temperature. These local environments also support the idea of water's polymorphism, i.e., the existence of the high- and low-density states in supercooled water.

Ultrafast Vibrational Echo Spectroscopy of Liquid Water from First-Principles Simulations

Vibrational echo spectroscopy has become a powerful technique to study vibrational spectral diffusion in water and aqueous solutions. The dynamics of vibrational spectral diffusion is intimately related to the hydrogen bond fluctuations in liquid water and other hydrogen bonded liquids. Earlier theoretical calculations of vibrational echo spectroscopy of aqueous systems were based on classical molecular dynamics simulations involving empirical force fields of water. In the current work, we have employed the method of ab initio molecular dynamics simulation to calculate the spectral observables of vibrational echo and two-dimensional infrared (2D-IR) spectroscopy of liquid water at room temperature under Condon and cumulant approximations. The time scales extracted from the temporal decay of the frequency-time correlation function (FTCF), short-time slope of three pulse photon echo (SP3E), dynamic line width (DLW), and the slope of nodal line of 2D-IR spectra are found to be in reasonably close agreement with each other which reinforces the assertion that signatures of FTCF can be captured using three pulse photon echo and 2D-IR spectroscopy.

Wetting of water on hexagonal boron nitride@Rh(111): a QM/MM model based on atomic charges derived for nano-structured substrates

The wetting of water on corrugated and flat hexagonal boron nitride (h-BN) monolayers on Rh(111) is studied within a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. Water is treated by QM methods, whereas the interactions between liquid and substrate are described at the MM level. The electrostatic properties of the substrate are reproduced by assigning specifically generated partial charges to each MM atom. We propose a method to determine restrained electrostatic potential (RESP) charges that can be applied to periodic systems. Our approach is based on the Gaussian and plane waves algorithm and allows an easy tuning of charges for nano-structured substrates. We have successfully applied it to reproduce the electrostatic potential of the corrugated and flat h-BN/Rh(111) known as nanomesh. Molecular dynamics simulations of water films in contact with these substrates are performed and structural and dynamic properties of the interfaces are analyzed. Based on this analysis and on the interaction energies between water film and substrate, we found that the corrugated nanomesh is slightly more hydrophilic. On a macroscopic scale, we expect a smaller contact angle for a droplet on the corrugated surface.