Ab initio molecular dynamics study of water at constant pressure using converged basis sets and empirical dispersion corrections

Characterization of the Hydrogen-Bond Network in High-Pressure Water by Deep Potential Molecular Dynamics

The hydrogen-bond (H-bond) network of high-pressure water is investigated by neural-network-based molecular dynamics (MD) simulations with first-principles accuracy. The static structure factors (SSFs) of water at three densities, i.e., 1, 1.115, and 1.24 g/cm3, are directly evaluated from 512 water MD trajectories, which are in quantitative agreement with the experiments. We propose a new method to decompose the computed SSF and identify the changes in the SSF with respect to the changes in H-bond structures. We find that a larger water density results in a higher probability for one or two non-H-bonded water molecules to be inserted into the inner shell, explaining the changes in the tetrahedrality of water under pressure. We predict that the structure of the accepting end of water molecules is more easily influenced by the pressure than by the donating end. Our work sheds new light on explaining the SSF and H-bond properties in related fields.

Hydroxide Diffusion in Functionalized Cylindrical Nanopores as Idealized Models of Anion Exchange Membrane Environments: An Ab Initio Molecular Dynamics Study

Anion exchange membranes (AEMs) have attracted significant interest for their applications in fuel cells and other electrochemical devices in recent years. Understanding water distributions and hydroxide transport mechanisms within AEMs is critical to improving their performance as concerns hydroxide conductivity. Recently, nanoconfined environments have been used to mimic AEM environments. Following this approach, we construct nanoconfined cylindrical pore structures using graphane nanotubes (GNs) functionalized with trimethylammonium cations as models of local AEM morphology. These structures were then used to investigate hydroxide transport using ab initio molecular dynamics (AIMD). The simulations showed that hydroxide transport is suppressed in these confined environments relative to the bulk solution although the mechanism is dominated by structural diffusion. One factor causing the suppressed hydroxide transport is the reduced proton transfer (PT) rates due to changes in hydroxide and water solvation patterns under confinement compared to bulk solution as well as strong interactions between hydroxide ions and the tethered cation groups.

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water has been a challenging task from ab initio molecular dynamics (AIMD). The difficulty mainly comes from a lack of accurate and efficient exchange-correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly-constrained and appropriately normed (SCAN) meta-GGA functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange-correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice in describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate hydroxide diffuses slower than hydronium in water, which is consistent with the experiments.

Assessing the Interplay between Functional-Driven and Density-Driven Errors in DFT Models of Water

We investigate the interplay between functional-driven and density-driven errors in different density functional approximations within density functional theory (DFT) and the implications of these errors for simulations of water with DFT-based data-driven potentials. Specifically, we quantify density-driven errors in two widely used dispersion-corrected functionals derived within the generalized gradient approximation (GGA), namely BLYP-D3 and revPBE-D3, and two modern meta-GGA functionals, namely strongly constrained and appropriately normed (SCAN) and B97M-rV. The effects of functional-driven and density-driven errors on the interaction energies are first assessed for the water clusters of the BEGDB dataset. Further insights into the nature of functional-driven errors are gained from applying the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) to the interaction energies, which demonstrates that functional-driven errors are strongly correlated with the nature of the interactions. We discuss cases where density-corrected DFT (DC-DFT) models display higher accuracy than the original DFT models and cases where reducing the density-driven errors leads to larger deviations from the reference energies due to the presence of large functional-driven errors. Finally, molecular dynamics simulations are performed with data-driven many-body potentials derived from DFT and DC-DFT data to determine the effect that minimizing density-driven errors has on the description of liquid water. Besides rationalizing the performance of widely used DFT models of water, we believe that our findings unveil fundamental relations between the shortcomings of some common DFT approximations and the requirements for accurate descriptions of molecular interactions, which will aid the development of a consistent, DFT-based framework for the development of data-driven and machine-learned potentials for simulations of condensed-phase systems.

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

Enhancing the Accuracy of Ab Initio Molecular Dynamics by Fine Tuning of Effective Two-Body Interactions: Acetonitrile as a Test Case

Grimme's dispersion-corrected density functional theory (DFT-D) methods have emerged among the most practical approaches to perform accurate quantum mechanical calculations on molecular systems ranging from small clusters to microscopic and mesoscopic samples, i.e., including hundreds or thousands of molecules. Moreover, DFT-D functionals can be easily integrated into popular ab initio molecular dynamics (MD) software packages to carry out first-principles condensed-phase simulations at an affordable computational cost. Here, starting from the well-established D3 version of the dispersion-correction term, we present a simple protocol to improve the accurate description of the intermolecular interactions of molecular clusters of growing size, considering acetonitrile as a test case. Optimization of the interaction energy was performed with reference to diffusion quantum Monte Carlo calculations, successfully reaching the same inherent accuracy of the latter (statistical error of ∼0.1 kcal/mol per molecule). The refined DFT-D3 model was then used to perform ab initio MD simulations of liquid acetonitrile, again showing significant improvements toward available experimental data with respect to the default correction.

Proton in the ring: spectroscopy and dynamics of proton bonding in macrocycle cavities

The proton bond is a paradigmatic quantum molecular interaction and a major driving force of supramolecular chemistry. The ring cavities of crown ethers provide an intriguing environment, promoting competitive proton sharing with multiple coordination anchors. This study shows that protons confined in crown ether cavities form dynamic bonds that migrate to varying pairs of coordinating atoms when allowed by the flexibility of the macrocycle backbone. Prototypic native crown ethers (12-crown-4, 15-crown-5 and 18-crown-6) and aza-crown ethers (cyclen, 1-aza-18-crown-6 and hexacyclen) are investigated. For each system, Infrared action spectroscopy experiments and ab initio Molecular Dynamics computations are employed to elucidate the structural effects associated with proton diffusion and its entanglement with the conformational and vibrational dynamics of the protonated host.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Effects of atomic species and interatomic distance on the interactions in one-dimensional nanomaterials

Noncovalent van der Waals (vdW) interactions are significant for the constitution of nanomaterials; however, they are not well understood in one-dimensional materials. Herein, we employ density functional theory (DFT) methods to address this issue and find that the many-body effects of vdW interactions within the one-dimensional wires composed of atoms chosen from the second period (B, C, N, O, F) vary with the interatomic distance of the wires. Furthermore, the atomic species effectively regulate the transition threshold of the many-body effects of vdW interactions. In the case of the adsorption of n-heptane (C7H16) on the wires, the atomic species alters the interactions between the wires and the molecule by modulating the coupling vibration between wires and C7H16 molecules. Correspondingly, replacing a portion of Pb with Tl atoms could contribute to the stability of the organic-inorganic hybrid halide perovskites with one-dimensional structures. Our findings not only contribute to the understanding of vdW interactions in one-dimensional structures with second-period atoms (B, C, N, O, F) but also provide clues for improving the stability of perovskites with one-dimensional structures.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide Salt on the Deprotonation of Formic Acid in Aqueous Solution

The deprotonation of formic acid is investigated using metadynamics in tandem with Born-Oppenheimer molecular dynamics simulations. We compare our findings for formic acid in pure water with previous studies before examining formic acid in aqueous solutions of lithium bromide. We carefully consider different definitions for the collective variable(s) used to drive the metadynamics, emphasizing that the variables used must include all of the possible reactive atoms in the system, in this case carboxylate oxygens and water hydrogens. This ensures that all the various possible proton exchange events can be accommodated and the collective variable(s) can distinguish the protonated and deprotonated states, even over rather long ab initio simulation runs (ca. 200-300 ps). Our findings show that the formic acid deprotonation barrier and the free energy of the deprotonated state are higher in concentrated lithium bromide, in agreement with the available experimental data for acids in salt solution. We show that the presence of Br^- in proximity to the formic acid hydroxyl group effectively inhibits deprotonation. Our study extends previous work on acid deprotonation in pure water and at air-water interfaces to more complex multicomponent systems of importance in atmospheric and marine chemistry.

Structure and Dynamics of Water at the Water-Air Interface Using First-Principles Molecular Dynamics Simulations. II. NonLocal vs Empirical van der Waals Corrections

van der Waals (vdW) correction schemes have been recognized to be essential for an accurate description of liquid water in first-principles molecular dynamics simulation. The description of the structure and dynamics of water is governed by the type of the vdW corrections. So far, two vdW correction schemes have been often used: empirical vdW corrections and nonlocal vdW corrections. In this paper, we assess the influence of the empirical vs nonlocal vdW correction schemes on the structure and dynamics of water at the water–air interface. Since the structure of water at the water–air interface is established by a delicate balance of hydrogen bond formation and breaking, the simulation at the water–air interface provides a unique platform to testify as to the heterogeneous interaction of water. We used the metrics [Ohto et al. J. Chem. Theory Comput., 2019, 15, 595−60230468702] which are directly connected with the sum-frequency generation spectroscopic measurement. We find that the overall performance of nonlocal vdW methods is either similar or worse compared to the empirical vdW methods. We also investigated the performance of the optB88-DRSLL functional, which showed slightly less accuracy than the revPBE-D3 method. We conclude that the revPBE-D3 method shows the best performance for describing the interfacial water.

Vibration Spectral Dynamics of Weakly Coordinating Water Molecules near an Anion: FPMD Simulations of an Aqueous Solution of Tetrafluoroborate

The extent to which the ions affect the nearby water molecules will decide the structure-making or breaking nature of those ions in aqueous solutions. The effects of a weakly coordinating anion on the structure, dynamics, and vibrational properties of water molecules are not so significant as compared to an anion capable of making strong ion-water hydrogen bonds. The present work deals with the first-principles molecular dynamics study of an aqueous solution of such a weakly coordinating anion, tetrafluoroborate (BF₄^-), using dispersion-corrected DFT-based first-principles molecular dynamics (FPMD) simulations. Various structural, dynamical, and spectral properties, such as radial distribution functions (RDFs), rotational dynamics, vibrational density of states (VDOS), hydrogen bond as well as dangling OH autocorrelation functions, and residence dynamics, were calculated to investigate the effects of the anion on nearby water molecules. The process of spectral diffusion was assessed through a time series wavelet transformation of trajectories obtained from FPMD simulations. The first ion-water solvation shell extends up to 5.5 Å, containing around 20 water molecules. The lifetime of the ion-water hydrogen bond is found to be 1.19 ps, whereas the water-water hydrogen bond lifetime is found to be 1.13 ps. Inside the solvation shell, the persistence time of dangling OH chromophores and the average frequency of OH modes inside the solvation shell are found to be more compared to bulk. Three time scales are found for solvation shell OH modes from the frequency-frequency correlation function. A very short time scale is found for the intact ion-water interaction; the short time scale is for the ion-water hydrogen bond, and the long time scale is for escape dynamics of water molecules from the ion solvation shell. From the mean squared displacement, it is found that solvation water molecules diffuse slower than the bulk. However, solvation shell water molecules show faster relaxation from the analysis of rotational anisotropy. Within the longer time scale of spectral diffusion, this process (which is related to various dynamics of the molecules) is not yet complete, as compared to fast anisotropic decay. This fact is similar to the experimental finding of spectral diffusion and anisotropy time scales in the aqueous solution of borohydride anion. The calculated results are also compared with available experimental data wherever possible.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide on the Structure of the Aqueous Solution-Air Interface

We present the results of ab initio molecular dynamics simulations of the solution–air interface of aqueous lithium bromide (LiBr). We find that, in agreement with the experimental data and previous simulation results with empirical polarizable force field models, Br^– anions prefer to accumulate just below the first molecular water layer near the interface, whereas Li⁺ cations remain deeply buried several molecular layers from the interface, even at very high concentration. The separation of ions has a profound effect on the average orientation of water molecules in the vicinity of the interface. We also find that the hydration number of Li⁺ cations in the center of the slab N_{c,Li⁺–H2O} ≈ 4.7 ± 0.3, regardless of the salt concentration. This estimate is consistent with the recent experimental neutron scattering data, confirming that results from nonpolarizable empirical models, which consistently predict tetrahedral coordination of Li⁺ to four solvent molecules, are incorrect. Consequently, disruption of the hydrogen bond network caused by Li⁺ may be overestimated in nonpolarizable empirical models. Overall, our results suggest that empirical models, in particular nonpolarizable models, may not capture all of the properties of the solution–air interface necessary to fully understand the interfacial chemistry.

Force Field for Water Based on Neural Network

We developed a novel neural network-based force field for water based on training with high-level ab initio theory. The force field was built based on an electrostatically embedded many-body expansion method truncated at binary interactions. The many-body expansion method is a common strategy to partition the total Hamiltonian of large systems into a hierarchy of few-body terms. Neural networks were trained to represent electrostatically embedded one-body and two-body interactions, which require as input only one and two water molecule calculations at the level of ab initio electronic structure method CCSD/aug-cc-pVDZ embedded in the molecular mechanics water environment, making it efficient as a general force field construction approach. Structural and dynamic properties of liquid water calculated with our force field show good agreement with experimental results. We constructed two sets of neural network based force fields: nonpolarizable and polarizable force fields. Simulation results show that the nonpolarizable force field using fixed TIP3P charges has already behaved well, since polarization effects and many-body effects are implicitly included due to the electrostatic embedding scheme. Our results demonstrate that the electrostatically embedded many-body expansion combined with neural network provides a promising and systematic way to build next-generation force fields at high accuracy and low computational costs, especially for large systems.

Advanced Potential Energy Surfaces for Molecular Simulation

Advanced potential energy surfaces are defined as theoretical models that explicitly include many-body effects that transcend the standard fixed-charge, pairwise-additive paradigm typically used in molecular simulation. However, several factors relating to their software implementation have precluded their widespread use in condensed-phase simulations: the computational cost of the theoretical models, a paucity of approximate models and algorithmic improvements that can ameliorate their cost, underdeveloped interfaces and limited dissemination in computational code bases that are widely used in the computational chemistry community, and software implementations that have not kept pace with modern high-performance computing (HPC) architectures, such as multicore CPUs and modern graphics processing units (GPUs). In this Feature Article we review recent progress made in these areas, including well-defined polarization approximations and new multipole electrostatic formulations, novel methods for solving the mutual polarization equations and increasing the MD time step, combining linear-scaling electronic structure methods with new QM/MM methods that account for mutual polarization between the two regions, and the greatly improved software deployment of these models and methods onto GPU and CPU hardware platforms. We have now approached an era where multipole-based polarizable force fields can be routinely used to obtain computational results comparable to state-of-the-art density functional theory while reaching sampling statistics that are acceptable when compared to that obtained from simpler fixed partial charge force fields.

Born-Oppenheimer Molecular Dynamics Simulations of a Bromate Ion in Water Reveal Its Dual Kosmotropic and Chaotropic Behavior

The solvation structure and dynamics of a bromate (BrO₃^-) ion in water are studied by means of Born-Oppenheimer molecular dynamics simulations at two different temperatures using the Becke-Lee-Yang-Parr functional with Grimme D3 dispersion corrections. The bromate ion possesses a pyramidal structure, and it has two types of solvation sites, namely, the bromine and oxygen atoms. We have looked at different radial and orientational distributions of water molecules around the bromate ion and also investigated their hydrogen bonding properties. The solvation structure of the bromate ion is also compared with that of the iodate (IO₃^-) ion, which is structurally rather similar to the bromate ion and was found to have some unusual solvation properties in water. It is found that the bromate ion follows a similar trend as that followed by the iodate ion as far as the solvation structure is concerned. However, the effect of the former on surrounding water is found to be much weaker than that of the latter. On the dynamical side, we have looked at diffusion, residence dynamics, and also the orientational and hydrogen bond relaxation of water molecules around the BrO₃^- ion and compared them with those of the bulk. Dynamical results are presented for both H₂O and D₂O around the BrO₃^- ion. Interpretation of the dynamical results in terms of structure-making (kosmotropic)/-breaking (chaotropic) properties of the BrO₃^- ion reveals that the bromine atom of this ion acts as a water structure breaker, whereas the three oxygens act as water structure makers. Thus, in spite of being a single ion, the bromate ion has dual characteristics and the experimentally observed kosmotropic ability of this ion is actually a trade-off between a chaotropic site (the bromine atom) and three kosmotropic sites (three oxygen atoms) that are present in the ion.

Validation of Molecular Simulation: An Overview of Issues

Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Q^exp of experimentally observable quantities Q. One must consider 1) the accuracy of Q^exp , 2) the accuracy of the function Q(r^N ) used to calculate a Q-value based on a molecular configuration r^N of N particles, 3) the sensitivity of the function Q(r^N ) to the configuration r^N , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.

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.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

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.

How van der Waals interactions determine the unique properties of water

Whereas the interactions between water molecules are dominated by strongly directional hydrogen bonds (HBs), it was recently proposed that relatively weak, isotropic van der Waals (vdW) forces are essential for understanding the properties of liquid water and ice. This insight was derived from ab initio computer simulations, which provide an unbiased description of water at the atomic level and yield information on the underlying molecular forces. However, the high computational cost of such simulations prevents the systematic investigation of the influence of vdW forces on the thermodynamic anomalies of water. Here, we develop efficient ab initio-quality neural network potentials and use them to demonstrate that vdW interactions are crucial for the formation of water's density maximum and its negative volume of melting. Both phenomena can be explained by the flexibility of the HB network, which is the result of a delicate balance of weak vdW forces, causing, e.g., a pronounced expansion of the second solvation shell upon cooling that induces the density maximum.

Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications

Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.

Water adsorption in SAPO-34: elucidating the role of local heterogeneities and defects using dispersion-corrected DFT calculations

The chabazite-type silicoaluminophosphate SAPO-34 is a promising adsorbent for applications in thermal energy storage using water adsorption-desorption cycles. In order to develop a microscopic understanding of the impact of local heterogeneities and defects on the water adsorption properties, the interaction of different models of SAPO-34 with water was studied using dispersion-corrected density-functional theory (DFT-D) calculations. In addition to SAPO-34 with isolated silicon atoms, the calculations considered models incorporating two types of heterogeneities (silicon islands, aluminosilicate domains), and two defect-containing (partially and fully desilicated) systems. DFT-D optimisations were performed for systems with small amounts of adsorbed water, in which all H2O molecules can interact with framework protons, and systems with large amounts of adsorbed water (30 H2O molecules per unit cell). At low loadings, the host-guest interaction energy calculated for SAPO-34 with isolated Si atoms amounts to approximately -90 kJ mol(-1). While the presence of local heterogeneities leads to the creation of some adsorption sites that are energetically slightly more favourable, the interaction strength is drastically reduced in systems with defects. At high water loadings, energies in the range of -70 kJ mol(-1) are obtained for all models. The DFT-D interaction energies are in good agreement with experimentally measured heats of water adsorption. A detailed analysis of the equilibrium structures was used to gain insights into the binding modes at low coverages, and to assess the extent of framework deprotonation and changes in the coordination environment of aluminium atoms at high water loadings.

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.

Improving intermolecular interactions in DFTB3 using extended polarization from chemical-potential equalization

Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMSD) interaction energy is improved from 6.07 kcal/mol to 1.49 kcal/mol for interactions with one charged species, whereas the RMSD is improved from 5.60 kcal/mol to 1.73 for a set of 9 salt bridges, compared to uncorrected DFTB3. For large water clusters and complexes that are dominated by dispersion interactions, the already satisfactory performance of the DFTB3-D3 model is retained; polarizabilities of neutral molecules are also notably improved. Overall, the CPE extension of DFTB3-D3 provides a more balanced description of different types of non-covalent interactions than Neglect of Diatomic Differential Overlap type of semi-empirical methods (e.g., PM6-D3H4) and PBE-D3 with modest basis sets.

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water

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 (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(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, POOO(θ), 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.

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

Molecular simulation of water and hydration effects in different environments: challenges and developments for DFTB based models

We discuss the description of water and hydration effects that employs an approximate density functional theory, DFTB3, in either a full QM or QM/MM framework. The goal is to explore, with the current formulation of DFTB3, the performance of this method for treating water in different chemical environments, the magnitude and nature of changes required to improve its performance, and factors that dictate its applicability to reactions in the condensed phase in a QM/MM framework. A relatively minor change (on the scale of kBT) in the O-H repulsive potential is observed to substantially improve the structural properties of bulk water under ambient conditions; modest improvements are also seen in dynamic properties of bulk water. This simple change also improves the description of protonated water clusters, a solvated proton, and to a more limited degree, a solvated hydroxide. By comparing results from DFTB3 models that differ in the description of water, we confirm that proton transfer energetics are adequately described by the standard DFTB3/3OB model for meaningful mechanistic analyses. For QM/MM applications, a robust parametrization of QM-MM interactions requires an explicit consideration of condensed phase properties, for which an efficient sampling technique was developed recently and is reviewed here. The discussions help make clear the value and limitations of DFTB3 based simulations, as well as the developments needed to further improve the accuracy and transferability of the methodology.