1
|
Montero de Hijes P, Dellago C, Jinnouchi R, Schmiedmayer B, Kresse G. Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks. J Chem Phys 2024; 160:114107. [PMID: 38506284 DOI: 10.1063/5.0197105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Accepted: 03/03/2024] [Indexed: 03/21/2024] Open
Abstract
In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.
Collapse
Affiliation(s)
- Pablo Montero de Hijes
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- University of Vienna, Faculty of Earth Sciences, Geography and Astronomy, Josef-Holaubuek-Platz 2, 1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
| | - Ryosuke Jinnouchi
- Toyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | | | - Georg Kresse
- University of Vienna, Faculty of Physics, Kolingasse 14, A-1090 Vienna, Austria
- VASP Software GmbH, Berggasse 21, A-1090 Vienna, Austria
| |
Collapse
|
2
|
Villard J, Bircher MP, Rothlisberger U. Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder. Chem Sci 2024; 15:4434-4451. [PMID: 38516095 PMCID: PMC10952088 DOI: 10.1039/d3sc05828j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 02/12/2024] [Indexed: 03/23/2024] Open
Abstract
The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.
Collapse
Affiliation(s)
- Justin Villard
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| | - Martin P Bircher
- Computational and Soft Matter Physics, Universität Wien Wien A-1090 Austria
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| |
Collapse
|
3
|
Borrelli W, Mei KJ, Park SJ, Schwartz BJ. Partial Molar Solvation Volume of the Hydrated Electron Simulated Via DFT. J Phys Chem B 2024; 128:2425-2431. [PMID: 38422045 PMCID: PMC10945486 DOI: 10.1021/acs.jpcb.3c05091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 02/03/2024] [Accepted: 02/14/2024] [Indexed: 03/02/2024]
Abstract
Different simulation models of the hydrated electron produce different solvation structures, but it has been challenging to determine which simulated solvation structure, if any, is the most comparable to experiment. In a recent work, Neupane et al. [J. Phys. Chem. B 2023, 127, 5941-5947] showed using Kirkwood-Buff theory that the partial molar volume of the hydrated electron, which is known experimentally, can be readily computed from an integral over the simulated electron-water radial distribution function. This provides a sensitive way to directly compare the hydration structure of different simulation models of the hydrated electron with experiment. Here, we compute the partial molar volume of an ab-initio-simulated hydrated electron model based on density-functional theory (DFT) with a hybrid functional at different simulated system sizes. We find that the partial molar volume of the DFT-simulated hydrated electron is not converged with respect to the system size for simulations with up to 128 waters. We show that even at the largest simulation sizes, the partial molar volume of DFT-simulated hydrated electrons is underestimated by a factor of 2 with respect to experiment, and at the standard 64-water size commonly used in the literature, DFT-based simulations underestimate the experimental solvation volume by a factor of ∼3.5. An extrapolation to larger box sizes does predict the experimental partial molar volume correctly; however, larger system sizes than those explored here are currently intractable without the use of machine-learned potentials. These results bring into question what aspects of the predicted hydrated electron radial distribution function, as calculated by DFT-based simulations with the PBEh-D3 functional, deviate from the true solvation structure.
Collapse
Affiliation(s)
| | | | - Sanghyun J. Park
- Department of Chemistry and
Biochemistry, University of California,
Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J. Schwartz
- Department of Chemistry and
Biochemistry, University of California,
Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|
4
|
Tal A, Bischoff T, Pasquarello A. Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections. Proc Natl Acad Sci U S A 2024; 121:e2311472121. [PMID: 38427604 PMCID: PMC10927489 DOI: 10.1073/pnas.2311472121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 01/31/2024] [Indexed: 03/03/2024] Open
Abstract
We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.
Collapse
Affiliation(s)
- Alexey Tal
- Chaire de Simulation à l’Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, LausanneCH-1015, Switzerland
| | - Thomas Bischoff
- Chaire de Simulation à l’Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, LausanneCH-1015, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l’Echelle Atomique, Ecole Polytechnique Fédérale de Lausanne, LausanneCH-1015, Switzerland
| |
Collapse
|
5
|
Kar R, Mandal S, Thakkur V, Meyer B, Nair NN. Speeding-up Hybrid Functional-Based Ab Initio Molecular Dynamics Using Multiple Time-stepping and Resonance-Free Thermostat. J Chem Theory Comput 2023; 19:8351-8364. [PMID: 37933121 DOI: 10.1021/acs.jctc.3c00964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2023]
Abstract
Ab initio molecular dynamics (AIMD) based on density functional theory (DFT) has become a workhorse for studying the structure, dynamics, and reactions in condensed matter systems. Currently, AIMD simulations are primarily carried out at the level of generalized gradient approximation (GGA), which is at the second rung of DFT functionals in terms of accuracy. Hybrid DFT functionals, which form the fourth rung in the accuracy ladder, are not commonly used in AIMD simulations as the computational cost involved is 100 times or higher. To facilitate AIMD simulations with hybrid functionals, we propose here an approach using multiple time stepping with adaptively compressed exchange operator and resonance-free thermostat, that could speed up the calculations by ∼30 times or more for systems with a few hundred of atoms. We demonstrate that by achieving this significant speed up and making the compute time of hybrid functional-based AIMD simulations at par with that of GGA functionals, we are able to study several complex condensed matter systems and model chemical reactions in solution with hybrid functionals that were earlier unthinkable to be performed.
Collapse
Affiliation(s)
- Ritama Kar
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India
| | - Sagarmoy Mandal
- Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Nägelsbachstr. 25, Erlangen 91052, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-Universität Erlangen-Nürnberg, Martensstr. 1, Erlangen 91058, Germany
| | - Vaishali Thakkur
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India
| | - Bernd Meyer
- Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Nägelsbachstr. 25, Erlangen 91052, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-Universität Erlangen-Nürnberg, Martensstr. 1, Erlangen 91058, Germany
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India
| |
Collapse
|
6
|
Di Liberto G, Giordano L. Role of solvation model on the stability of oxygenates on Pt(111): A comparison between microsolvation, extended bilayer, and extended metal/water interface. ELECTROCHEMICAL SCIENCE ADVANCES 2023. [DOI: 10.1002/elsa.202100204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Affiliation(s)
| | - Livia Giordano
- Department of Materials Science University of Milano‐Bicocca Milano Italy
| |
Collapse
|
7
|
Liu R, Zhang C, Liang X, Liu J, Wu X, Chen M. Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling. J Chem Phys 2022; 157:024503. [DOI: 10.1063/5.0094944] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | | | | | | | - Xifan Wu
- Physics, Temple University, United States of America
| | - Mohan Chen
- College of Engineering, Peking University, China
| |
Collapse
|
8
|
Zhang L, Wang H, Muniz MC, Panagiotopoulos AZ, Car R, E W. A deep potential model with long-range electrostatic interactions. J Chem Phys 2022; 156:124107. [DOI: 10.1063/5.0083669] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.
Collapse
Affiliation(s)
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, People’s Republic of China
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | - Weinan E
- School of Mathematical Sciences, Peking University, Beijing 100871, People’s Republic of China
- AI for Science Institute, Beijing, People’s Republic of China
- Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
9
|
Mandal S, Kar R, Klöffel T, Meyer B, Nair NN. Improving the scaling and performance of multiple time stepping-based molecular dynamics with hybrid density functionals. J Comput Chem 2022; 43:588-597. [PMID: 35147988 DOI: 10.1002/jcc.26816] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Revised: 01/07/2022] [Accepted: 01/18/2022] [Indexed: 12/18/2022]
Abstract
Density functionals at the level of the generalized gradient approximation (GGA) and a plane-wave basis set are widely used today to perform ab initio molecular dynamics (AIMD) simulations. Going up in the ladder of accuracy of density functionals from GGA (second rung) to hybrid density functionals (fourth rung) is much desired pertaining to the accuracy of the latter in describing structure, dynamics, and energetics of molecular and condensed matter systems. On the other hand, hybrid density functional based AIMD simulations are about two orders of magnitude slower than GGA based AIMD for systems containing ~100 atoms using ~100 compute cores. Two methods, namely MTACE and s-MTACE, based on a multiple time step integrator and adaptively compressed exchange operator formalism are able to provide a speed-up of about 7-9 in performing hybrid density functional based AIMD. In this work, we report an implementation of these methods using a task-group based parallelization within the CPMD program package, with the intention to take advantage of the large number of compute cores available on modern high-performance computing platforms. We present here the boost in performance achieved through this algorithm. This work also identifies the computational bottleneck in the s-MTACE method and proposes a way to overcome it.
Collapse
Affiliation(s)
- Sagarmoy Mandal
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India.,Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.,Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Ritama Kar
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| | - Tobias Klöffel
- Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.,Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Bernd Meyer
- Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.,Erlangen National High Performance Computing Center (NHR@FAU), Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur, India
| |
Collapse
|
10
|
Priyadarsini A, Mallik BS. Structure and rotational dynamics of water around hydrogen peroxide. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
11
|
Ringe S, Hörmann NG, Oberhofer H, Reuter K. Implicit Solvation Methods for Catalysis at Electrified Interfaces. Chem Rev 2021; 122:10777-10820. [PMID: 34928131 PMCID: PMC9227731 DOI: 10.1021/acs.chemrev.1c00675] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Implicit solvation
is an effective, highly coarse-grained approach
in atomic-scale simulations to account for a surrounding liquid electrolyte
on the level of a continuous polarizable medium. Originating in molecular
chemistry with finite solutes, implicit solvation techniques are now
increasingly used in the context of first-principles modeling of electrochemistry
and electrocatalysis at extended (often metallic) electrodes. The
prevalent ansatz to model the latter electrodes and the reactive surface
chemistry at them through slabs in periodic boundary condition supercells
brings its specific challenges. Foremost this concerns the difficulty
of describing the entire double layer forming at the electrified solid–liquid
interface (SLI) within supercell sizes tractable by commonly employed
density functional theory (DFT). We review liquid solvation methodology
from this specific application angle, highlighting in particular its
use in the widespread ab initio thermodynamics approach
to surface catalysis. Notably, implicit solvation can be employed
to mimic a polarization of the electrode’s electronic density
under the applied potential and the concomitant capacitive charging
of the entire double layer beyond the limitations of the employed
DFT supercell. Most critical for continuing advances of this effective
methodology for the SLI context is the lack of pertinent (experimental
or high-level theoretical) reference data needed for parametrization.
Collapse
Affiliation(s)
- Stefan Ringe
- Department of Energy Science and Engineering, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea.,Energy Science & Engineering Research Center, Daegu Institute of Science and Technology (DGIST), Daegu 42988, Republic of Korea
| | - Nicolas G Hörmann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany.,Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstraße 4, D-85747 Garching, Germany.,Chair for Theoretical Physics VII and Bavarian Center for Battery Technology (BayBatt), University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Karsten Reuter
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
12
|
Liu J, He X. Ab initio molecular dynamics simulation of liquid water with fragment-based quantum mechanical approach under periodic boundary conditions. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing 210009, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry at New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
13
|
Gittus OR, Bresme F. Thermophysical properties of water using reactive force fields. J Chem Phys 2021; 155:114501. [PMID: 34551553 DOI: 10.1063/5.0057868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.
Collapse
Affiliation(s)
- Oliver R Gittus
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| | - Fernando Bresme
- Department of Chemistry, Molecular Sciences Research Hub Imperial College, London W12 0BZ, United Kingdom
| |
Collapse
|
14
|
Krishnamoorthy A, Nomura KI, Baradwaj N, Shimamura K, Rajak P, Mishra A, Fukushima S, Shimojo F, Kalia R, Nakano A, Vashishta P. Dielectric Constant of Liquid Water Determined with Neural Network Quantum Molecular Dynamics. PHYSICAL REVIEW LETTERS 2021; 126:216403. [PMID: 34114857 DOI: 10.1103/physrevlett.126.216403] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Accepted: 03/30/2021] [Indexed: 06/12/2023]
Abstract
The static dielectric constant ϵ_{0} and its temperature dependence for liquid water is investigated using neural network quantum molecular dynamics (NNQMD). We compute the exact dielectric constant in canonical ensemble from NNQMD trajectories using fluctuations in macroscopic polarization computed from maximally localized Wannier functions (MLWF). Two deep neural networks are constructed. The first, NNQMD, is trained on QMD configurations for liquid water under a variety of temperature and density conditions to learn potential energy surface and forces and then perform molecular dynamics simulations. The second network, NNMLWF, is trained to predict locations of MLWF of individual molecules using the atomic configurations from NNQMD. Training data for both the neural networks is produced using a highly accurate quantum-mechanical method, DFT-SCAN that yields an excellent description of liquid water. We produce 280×10^{6} configurations of water at 7 temperatures using NNQMD and predict MLWF centers using NNMLWF to compute the polarization fluctuations. The length of trajectories needed for a converged value of the dielectric constant at 0°C is found to be 20 ns (40×10^{6} configurations with 0.5 fs time step). The computed dielectric constants for 0, 15, 30, 45, 60, 75, and 90°C are in good agreement with experiments. Our scalable scheme to compute dielectric constants with quantum accuracy is also applicable to other polar molecular liquids.
Collapse
Affiliation(s)
- Aravind Krishnamoorthy
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Ken-Ichi Nomura
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Nitish Baradwaj
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Kohei Shimamura
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Pankaj Rajak
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Ankit Mishra
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Shogo Fukushima
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Fuyuki Shimojo
- Department of Physics, Kumamoto University, Kumamoto 860-8555, Japan
| | - Rajiv Kalia
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Aiichiro Nakano
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| | - Priya Vashishta
- Collaboratory for Advanced Computing and Simulations, Department of Chemical Engineering and Materials Science, Department of Physics & Astronomy, and Department of Computer Science, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
15
|
Piaggi PM, Panagiotopoulos AZ, Debenedetti PG, Car R. Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional. J Chem Theory Comput 2021; 17:3065-3077. [PMID: 33835819 DOI: 10.1021/acs.jctc.1c00041] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.
Collapse
Affiliation(s)
- Pablo M Piaggi
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Athanassios Z Panagiotopoulos
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States.,Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States.,Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, United States
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Princeton University, Princeton, New Jersey 08544, United States.,Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
16
|
Drużbicki K, Gaboardi M, Fernandez-Alonso F. Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities. Polymers (Basel) 2021; 13:1440. [PMID: 33947108 PMCID: PMC8125526 DOI: 10.3390/polym13091440] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 04/27/2021] [Indexed: 12/19/2022] Open
Abstract
This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.
Collapse
Affiliation(s)
- Kacper Drużbicki
- Materials Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastian, Spain;
- Polish Academy of Sciences, Center of Molecular and Macromolecular Studies, Sienkiewicza 112, 90-363 Lodz, Poland
| | - Mattia Gaboardi
- Elettra—Sincrotrone Trieste S.C.p.A., S.S. 14 km 163.5 in Area Science Park, 34149 Trieste, Italy;
| | - Felix Fernandez-Alonso
- Materials Physics Center, CSIC-UPV/EHU, Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastian, Spain;
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastian, Spain
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
- IKERBASQUE, Basque Foundation for Science, Plaza Euskadi 5, 48009 Bilbao, Spain
| |
Collapse
|
17
|
Mandal S, Thakkur V, Nair NN. Achieving an Order of Magnitude Speedup in Hybrid-Functional- and Plane-Wave-Based Ab Initio Molecular Dynamics: Applications to Proton-Transfer Reactions in Enzymes and in Solution. J Chem Theory Comput 2021; 17:2244-2255. [PMID: 33740375 DOI: 10.1021/acs.jctc.1c00009] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ab initio molecular dynamics (MD) with hybrid density functionals and a plane wave basis is computationally expensive due to the high computational cost of exact exchange energy evaluation. Recently, we proposed a strategy to combine adaptively compressed exchange (ACE) operator formulation and a multiple time step integration scheme to reduce the computational cost significantly [J. Chem. Phys. 2019, 151, 151102 ]. However, it was found that the construction of the ACE operator, which has to be done at least once in every MD time step, is computationally expensive. In the present work, systematic improvements are introduced to further speed up by employing localized orbitals for the construction of the ACE operator. By this, we could achieve a computational speedup of an order of magnitude for a periodic system containing 32 water molecules. Benchmark calculations were carried out to show the accuracy and efficiency of the method in predicting the structural and dynamical properties of bulk water. To demonstrate the applicability, computationally intensive free-energy computations at the level of hybrid density functional theory were performed to investigate (a) methyl formate hydrolysis reaction in neutral aqueous media and (b) proton-transfer reaction within the active-site residues of the class C β-lactamase enzyme.
Collapse
Affiliation(s)
- Sagarmoy Mandal
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India.,Interdisciplinary Center for Molecular Materials and Computer Chemistry Center, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Nägelsbachstr. 25, Erlangen 91052, Germany
| | - Vaishali Thakkur
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur (IITK), Kanpur 208016, India
| |
Collapse
|
18
|
Contreras R, Lodeiro L, Rozas-Castro N, Ormazábal-Toledo R. On the role of water in the hydrogen bond network in DESs: an ab initio molecular dynamics and quantum mechanical study on the urea–betaine system. Phys Chem Chem Phys 2021; 23:1994-2004. [DOI: 10.1039/d0cp06078j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We herein report an ab initio molecular dynamics study on a natural DES composed of urea and betaine in a 3 : 2 ratio, as a test case for evaluating the water effect.
Collapse
Affiliation(s)
- Renato Contreras
- Departamento de Química, Facultad de Ciencias
- Universidad de Chile
- Las Palmeras 3425
- Casilla 653
- Chile
| | - Lucas Lodeiro
- Departamento de Química, Facultad de Ciencias
- Universidad de Chile
- Las Palmeras 3425
- Casilla 653
- Chile
| | - Nicolás Rozas-Castro
- Departamento de Química, Facultad de Ciencias
- Universidad de Chile
- Las Palmeras 3425
- Casilla 653
- Chile
| | - Rodrigo Ormazábal-Toledo
- Departamento de Química, Facultad de Ciencias
- Universidad de Chile
- Las Palmeras 3425
- Casilla 653
- Chile
| |
Collapse
|
19
|
Cassone G. Nuclear Quantum Effects Largely Influence Molecular Dissociation and Proton Transfer in Liquid Water under an Electric Field. J Phys Chem Lett 2020; 11:8983-8988. [PMID: 33035059 DOI: 10.1021/acs.jpclett.0c02581] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Proton transfer in liquid water controls acid-base chemistry, crucial enzyme reactions, and the functioning of fuel cells. Externally applied static electric fields in water are capable of dissociating molecules and transferring protons across the H-bond network. However, the impact of nuclear quantum effects (NQEs) on these fundamental field-induced phenomena has not yet been reported. By comparing state-of-the-art ab initio molecular dynamics (AIMD) and path integral AIMD simulations of water under electric fields, I show that quantum delocalization of the proton lowers the molecular ionization threshold to approximately one-third. Moreover, also the water behavior as a protonic semiconductor is considerably modified by the inclusion of NQEs. In fact, when the quantum nature of the nuclei is taken into account, the proton conductivity is ∼50% larger. This work proves that NQEs sizably affect the protolysis phenomenon and proton transfer in room-temperature liquid water.
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute for Chemical-Physical Processes, National Research Council, Viale F. Stagno d'Alcontres 37, 98158 Messina, Italy
| |
Collapse
|
20
|
Calio PB, Hocky GM, Voth GA. Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water. J Chem Theory Comput 2020; 16:5675-5684. [DOI: 10.1021/acs.jctc.0c00558] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Paul B. Calio
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Glen M. Hocky
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
21
|
Ab Initio Molecular Dynamics Study of Methanol-Water Mixtures under External Electric Fields. Molecules 2020; 25:molecules25153371. [PMID: 32722281 PMCID: PMC7435743 DOI: 10.3390/molecules25153371] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 07/17/2020] [Accepted: 07/20/2020] [Indexed: 11/25/2022] Open
Abstract
Intense electric fields applied on H-bonded systems are able to induce molecular dissociations, proton transfers, and complex chemical reactions. Nevertheless, the effects induced in heterogeneous molecular systems such as methanol-water mixtures are still elusive. Here we report on a series of state-of-the-art ab initio molecular dynamics simulations of liquid methanol-water mixtures at different molar ratios exposed to static electric fields. If, on the one hand, the presence of water increases the proton conductivity of methanol-water mixtures, on the other, it hinders the typical enhancement of the chemical reactivity induced by electric fields. In particular, a sudden increase of the protonic conductivity is recorded when the amount of water exceeds that of methanol in the mixtures, suggesting that important structural changes of the H-bond network occur. By contrast, the field-induced multifaceted chemistry leading to the synthesis of e.g., hydrogen, dimethyl ether, formaldehyde, and methane observed in neat methanol, in 75:25, and equimolar methanol-water mixtures, completely disappears in samples containing an excess of water and in pure water. The presence of water strongly inhibits the chemical reactivity of methanol.
Collapse
|
22
|
Mandal S, Nair NN. Efficient computation of free energy surfaces of chemical reactions using ab initio molecular dynamics with hybrid functionals and plane waves. J Comput Chem 2020; 41:1790-1797. [PMID: 32407582 DOI: 10.1002/jcc.26222] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Revised: 04/06/2020] [Accepted: 04/29/2020] [Indexed: 11/10/2022]
Abstract
Ab initio molecular dynamics (AIMD) simulations employing density functional theory (DFT) and plane waves are routinely carried out using density functionals at the level of generalized gradient approximation (GGA). AIMD simulations employing hybrid density functionals are of great interest as it offers a more accurate description of structural and dynamic properties than the GGA functionals. However, the computational cost for carrying out calculations using hybrid functionals and plane wave basis set is at least two orders of magnitude higher than that using GGA functionals. Recently, we proposed a strategy that combined the adaptively compressed exchange operator formulation and the multiple time step integration scheme to reduce the computational cost by an order of magnitude [J. Chem. Phys. 151, 151102 (2019)]. In this work, we demonstrate the application of this method to study chemical reactions, in particular, formamide hydrolysis in an alkaline aqueous medium. By actuating our implementation with the well-sliced metadynamics scheme, we can compute the two-dimensional free energy surface of this reaction at the level of hybrid-DFT. This work also investigates the accuracy of the PBE0 (hybrid) and the PBE (GGA) functionals in predicting the free energetics of this chemical reaction.
Collapse
Affiliation(s)
- Sagarmoy Mandal
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| |
Collapse
|
23
|
Priyadarsini A, Dasari S, Mallik BS. Thermophysical Properties and Angular Jump Dynamics of Water: A Comparative DFT and DFT-Dispersion-Based Molecular Dynamics Study. J Phys Chem A 2020; 124:6039-6049. [DOI: 10.1021/acs.jpca.0c02909] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Adyasa Priyadarsini
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi 502285, Sangareddy, Telangana, India
| | - Sathish Dasari
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi 502285, Sangareddy, Telangana, India
| | - Bhabani S. Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi 502285, Sangareddy, Telangana, India
| |
Collapse
|
24
|
Ko HY, Jia J, Santra B, Wu X, Car R, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance. J Chem Theory Comput 2020; 16:3757-3785. [PMID: 32045232 DOI: 10.1021/acs.jctc.9b01167] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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 description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Princeton University, Princeton, New Jersey 08544, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
25
|
Self-interaction error overbinds water clusters but cancels in structural energy differences. Proc Natl Acad Sci U S A 2020; 117:11283-11288. [PMID: 32393631 DOI: 10.1073/pnas.1921258117] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.
Collapse
|
26
|
Priyadarshini A, Biswas A, Chakraborty D, Mallik BS. Structural and Thermophysical Anomalies of Liquid Water: A Tale of Molecules in the Instantaneous Low- and High-Density Regions. J Phys Chem B 2020; 124:1071-1081. [DOI: 10.1021/acs.jpcb.9b11596] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Adyasa Priyadarshini
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285 Sangareddy, Telangana, India
| | - Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285 Sangareddy, Telangana, India
| | - Debashree Chakraborty
- Department of Chemistry, National Institute of Technology Karnataka, Surathkal, 575025 Mangalore, Karnataka, India
| | - Bhabani S. Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285 Sangareddy, Telangana, India
| |
Collapse
|
27
|
Sakti AW, Nishimura Y, Nakai H. Recent advances in quantum‐mechanical molecular dynamics simulations of proton transfer mechanism in various water‐based environments. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1419] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Aditya W. Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
| | - Hiromi Nakai
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering Waseda University Tokyo Japan
| |
Collapse
|
28
|
Zhao G, Shi S, Xie H, Xu Q, Ding M, Zhao X, Yan J, Wang D. Equation of state of water based on the SCAN meta-GGA density functional. Phys Chem Chem Phys 2020; 22:4626-4631. [DOI: 10.1039/c9cp06362e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
By ab initio molecular dynamics simulations, the newly developed SCAN meta-GGA functional is proved better than the widely used PBE-GGA functional in describing the equation of state of water.
Collapse
Affiliation(s)
- Gang Zhao
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Shuyi Shi
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Huijuan Xie
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Qiushuang Xu
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Mingcui Ding
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Xuguang Zhao
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Jinliang Yan
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| | - Dehua Wang
- School of Physics and Optoelectronic Engineering
- Ludong University
- Yantai 264025
- P. R. China
| |
Collapse
|
29
|
LaCount MD, Gygi F. Ensemble first-principles molecular dynamics simulations of water using the SCAN meta-GGA density functional. J Chem Phys 2019; 151:164101. [DOI: 10.1063/1.5124957] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Michael D. LaCount
- Department of Computer Science, University of California, Davis, California 95616, USA
| | - François Gygi
- Department of Computer Science, University of California, Davis, California 95616, USA
| |
Collapse
|
30
|
Ko HY, Zhang L, Santra B, Wang H, E W, DiStasio Jr RA, Car R. Isotope effects in liquid water via deep potential molecular dynamics. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652366] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ, USA
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, PA, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing, People's Republic of China
| | - Weinan E
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Mathematics, Princeton University, Princeton, NJ, USA
| | | | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ, USA
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Physics and Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ, USA
| |
Collapse
|
31
|
Biswas A, Priyadarsini A, Mallik BS. Dynamics and Spectral Response of Water Molecules around Tetramethylammonium Cation. J Phys Chem B 2019; 123:8753-8766. [DOI: 10.1021/acs.jpcb.9b05466] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285, Sangareddy, Telangana India
| | - Adyasa Priyadarsini
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285, Sangareddy, Telangana India
| | - Bhabani S. Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi-502285, Sangareddy, Telangana India
| |
Collapse
|
32
|
Nguyen NL, Ma H, Govoni M, Gygi F, Galli G. Finite-Field Approach to Solving the Bethe-Salpeter Equation. PHYSICAL REVIEW LETTERS 2019; 122:237402. [PMID: 31298883 DOI: 10.1103/physrevlett.122.237402] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Revised: 12/13/2018] [Indexed: 05/26/2023]
Abstract
We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation and the calculation of the screened Coulomb interaction in a finite field. The method does not require either the explicit evaluation of dielectric matrices or of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition, it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wave functions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.
Collapse
Affiliation(s)
- Ngoc Linh Nguyen
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - He Ma
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Marco Govoni
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
- Materials Science Division and Institute for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Francois Gygi
- Department of Computer Science, University of California Davis, Davis, California 95616, USA
| | - Giulia Galli
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
- Materials Science Division and Institute for Molecular Engineering, Argonne National Laboratory, Lemont, Illinois 60439, USA
| |
Collapse
|
33
|
Elton DC, Fritz M, Fernández-Serra M. Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost. Phys Chem Chem Phys 2018; 21:409-417. [PMID: 30534683 DOI: 10.1039/c8cp06077k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.
Collapse
Affiliation(s)
- Daniel C Elton
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA.
| | | | | |
Collapse
|
34
|
Mandal S, Debnath J, Meyer B, Nair NN. Enhanced sampling and free energy calculations with hybrid functionals and plane waves for chemical reactions. J Chem Phys 2018; 149:144113. [DOI: 10.1063/1.5049700] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sagarmoy Mandal
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Jayashrita Debnath
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Bernd Meyer
- Interdisciplinary Center of Molecular Materials (ICMM) and Computer-Chemistry-Center (CCC), Friedrich-Alexander-Universität Erlangen-Nürnberg, Nägelsbachstraße 25, 91052 Erlangen, Germany
| | - Nisanth N. Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
35
|
Ruiz Pestana L, Marsalek O, Markland TE, Head-Gordon T. The Quest for Accurate Liquid Water Properties from First Principles. J Phys Chem Lett 2018; 9:5009-5016. [PMID: 30118601 DOI: 10.1021/acs.jpclett.8b02400] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of density functional theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, display accuracy comparable to a computationally expensive dispersion-corrected hybrid functional, revPBE0-D3. We show that the meta-GGA DFT functional reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared spectra with comparable accuracy of a much more expensive hybrid functional. This work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei and that high-temperature simulations are a poor proxy for properly treating NQEs.
Collapse
Affiliation(s)
| | - Ondrej Marsalek
- Charles University , Faculty of Mathematics and Physics , Ke Karlovu 3 , 121 16 Prague 2 , Czech Republic
| | - Thomas E Markland
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | | |
Collapse
|
36
|
Lopes PA, Santaella GM, Lima CAS, Vasconcelos KDF, Groppo FC. Evaluation of soft tissues simulant materials in cone beam computed tomography. Dentomaxillofac Radiol 2018; 48:20180072. [PMID: 30004256 DOI: 10.1259/dmfr.20180072] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVES To evaluate different materials in simulating soft tissues and to analyze the influence of these materials on the mean (MPIV) and standard deviation of pixel intensity values comparing them to a gold-standard in CBCT images. METHODS Images of three piglet heads with their soft tissues intact (gold-standard) and different simulant materials were acquired: ice, modelling wax, and ballistic gelatin, with the same thickness of the original soft tissues. The pixel intensities were measured in dental, bone and soft tissues regions, in the mandible and maxilla, for all the groups. Analysis of variance, Dunnet's, Pearson's and linear regression tests were performed. RESULTS The simulators did not significantly change the MPIV of teeth in comparison with the gold-standard (p = 0.1017). Only ice (p = 0.0156) affected the MPIV of bone. Wax (p = 0.001) and ice (p = 0.0076), but not ballistic gelatin (p = 0.5814), altered the MPIV of soft tissue regions. When assessing the influence of the location (mandible or maxilla) among the simulants, the differences were significant only for the soft tissue regions. Standard deviation was not influenced by simulants (p > 0.05), but ballistic gelatin presented the lower variability. CONCLUSIONS The ballistic gelatin was the best soft tissue simulant since it had the lowest influence on the pixel intensity values for all regions.
Collapse
Affiliation(s)
- Priscila A Lopes
- Department of Oral Diagnosis, Area of Oral Radiology, Piracicaba Dental School, University of Campinas, Piracicaba, Sao Paulo, Brazil
| | - Gustavo M Santaella
- Department of Oral Diagnosis, Area of Oral Radiology, Piracicaba Dental School, University of Campinas, Piracicaba, Sao Paulo, Brazil
| | - Carlos Augusto S Lima
- Department of Oral Diagnosis, Area of Oral Radiology, Piracicaba Dental School, University of Campinas, Piracicaba, Sao Paulo, Brazil
| | - Karla de Faria Vasconcelos
- Department of Oral Diagnosis, Area of Oral Radiology, Piracicaba Dental School, University of Campinas, Piracicaba, Sao Paulo, Brazil
| | - Francisco C Groppo
- Department of Physiological Sciences, University of Campinas, Piracicaba, São Paulo, Brazil
| |
Collapse
|
37
|
Abstract
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.
Collapse
Affiliation(s)
- Hao Wang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry and Department of Physics, Duke University, Durham, North Carolina 27708, United States
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
38
|
Goel H, Ling S, Ellis BN, Taconi A, Slater B, Rai N. Predicting vapor liquid equilibria using density functional theory: A case study of argon. J Chem Phys 2018; 148:224501. [PMID: 29907054 DOI: 10.1063/1.5025726] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Sanliang Ling
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Breanna Nicole Ellis
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Anna Taconi
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| |
Collapse
|
39
|
Gaiduk AP, Gustafson J, Gygi F, Galli G. First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional. J Phys Chem Lett 2018; 9:3068-3073. [PMID: 29768015 DOI: 10.1021/acs.jpclett.8b01017] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.
Collapse
Affiliation(s)
- Alex P Gaiduk
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
| | - Jeffrey Gustafson
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| | - François Gygi
- Department of Computer Science , University of California , Davis , California 95616 , United States
| | - Giulia Galli
- Institute for Molecular Engineering , The University of Chicago , Chicago , Illinois 60637 , United States
- Materials Science Division , Argonne National Laboratory , Argonne , Illinois 60439 , United States
- Department of Chemistry , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
40
|
Zheng L, Chen M, Sun Z, Ko HY, Santra B, Dhuvad P, Wu X. Structural, electronic, and dynamical properties of liquid water by ab initio molecular dynamics based on SCAN functional within the canonical ensemble. J Chem Phys 2018; 148:164505. [DOI: 10.1063/1.5023611] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lixin Zheng
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Mohan Chen
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Zhaoru Sun
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Pratikkumar Dhuvad
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
41
|
Kenmoe S, Lisovski O, Piskunov S, Bocharov D, Zhukovskii YF, Spohr E. Water Adsorption on Clean and Defective Anatase TiO2 (001) Nanotube Surfaces: A Surface Science Approach. J Phys Chem B 2018; 122:5432-5440. [DOI: 10.1021/acs.jpcb.7b11697] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stephane Kenmoe
- Department of Theoretical Chemistry, University of Duisburg-Essen, Universitätsstr. 2, D-45141 Essen, Germany
| | - Oleg Lisovski
- Institute of Solid State Physics, University of Latvia, 8 Kengaraga Str., Riga LV-1063, Latvia
| | - Sergei Piskunov
- Institute of Solid State Physics, University of Latvia, 8 Kengaraga Str., Riga LV-1063, Latvia
| | - Dmitry Bocharov
- Institute of Solid State Physics, University of Latvia, 8 Kengaraga Str., Riga LV-1063, Latvia
| | - Yuri F. Zhukovskii
- Institute of Solid State Physics, University of Latvia, 8 Kengaraga Str., Riga LV-1063, Latvia
| | - Eckhard Spohr
- Department of Theoretical Chemistry, University of Duisburg-Essen, Universitätsstr. 2, D-45141 Essen, Germany
- Center of Computational Sciences and Simulation, University of Duisburg-Essen, D-45141 Essen, Germany
| |
Collapse
|
42
|
Machida M, Kato K, Shiga M. Nuclear quantum effects of light and heavy water studied by all-electron first principles path integral simulations. J Chem Phys 2018; 148:102324. [DOI: 10.1063/1.5000091] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Affiliation(s)
- Masahiko Machida
- CCSE, Japan Atomic Energy Agency (JAEA), 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Koichiro Kato
- Mizuho Information and Research Institute, Inc., 2-3, Kandanishiki-cho, Chiyoda-ku, Tokyo 101-8443, Japan
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency (JAEA), 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
43
|
Hydroxide diffuses slower than hydronium in water because its solvated structure inhibits correlated proton transfer. Nat Chem 2018. [DOI: 10.1038/s41557-018-0010-2] [Citation(s) in RCA: 133] [Impact Index Per Article: 22.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
44
|
|
45
|
Abstract
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.
Collapse
|
46
|
Galib M, Duignan TT, Misteli Y, Baer MD, Schenter GK, Hutter J, Mundy CJ. Mass density fluctuations in quantum and classical descriptions of liquid water. J Chem Phys 2017; 146:244501. [DOI: 10.1063/1.4986284] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Mirza Galib
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Timothy T. Duignan
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Yannick Misteli
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Jürg Hutter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, USA
| |
Collapse
|
47
|
Zhang L, Li W, Fang T, Li S. Accurate Relative Energies and Binding Energies of Large Ice–Liquid Water Clusters and Periodic Structures. J Phys Chem A 2017; 121:4030-4038. [DOI: 10.1021/acs.jpca.7b03376] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lei Zhang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Wei Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Tao Fang
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- Institute of Theoretical
and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry
of MOE, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
48
|
Mao Y, Shao Y, Dziedzic J, Skylaris CK, Head-Gordon T, Head-Gordon M. Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis. J Chem Theory Comput 2017; 13:1963-1979. [PMID: 28430427 DOI: 10.1021/acs.jctc.7b00089] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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 (NH3, NH4+, 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.
Collapse
Affiliation(s)
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K.,Faculty of Applied Physics and Mathematics, Gdańsk University of Technology , Gdańsk 80-233, Poland
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K
| | | | - Martin Head-Gordon
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
49
|
Ruiz Pestana L, Mardirossian N, Head-Gordon M, Head-Gordon T. Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals. Chem Sci 2017; 8:3554-3565. [PMID: 30155200 PMCID: PMC6092720 DOI: 10.1039/c6sc04711d] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2016] [Accepted: 02/24/2017] [Indexed: 01/23/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Luis Ruiz Pestana
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , USA .
| | - Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , USA .
- Kenneth S. Pitzer Center for Theoretical Chemistry , Department of Chemistry , University of California , Berkeley , USA
- Departments of Chemistry , Bioengineering , Chemical and Biomolecular Engineering , University of California , Berkeley , USA
| |
Collapse
|
50
|
Gaiduk AP, Galli G. Local and Global Effects of Dissolved Sodium Chloride on the Structure of Water. J Phys Chem Lett 2017; 8:1496-1502. [PMID: 28267335 DOI: 10.1021/acs.jpclett.7b00239] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Determining how the structure of water is modified by the presence of salts is instrumental to understanding the solvation of biomolecules and, in general, the role played by salts in biochemical processes. However, the extent of hydrogen bonding disruption induced by salts remains controversial. We performed extensive first-principles simulations of solutions of a simple salt (NaCl) and found that, while the cation does not significantly change the structure of water beyond the first solvation shell, the anion has a further reaching effect, modifying the hydrogen-bond network even outside its second solvation shell. We found that a distinctive fingerprint of hydrogen bonding modification is the change in polarizability of water molecules. Molecular dipole moments are instead insensitive probes of long-range modifications induced by Na+ and Cl- ions. Though noticeable, the long-range effect of Cl- is expected to be too weak to affect solubility of large biomolecules.
Collapse
Affiliation(s)
- Alex P Gaiduk
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States
| | - Giulia Galli
- Institute for Molecular Engineering, The University of Chicago , Chicago, Illinois 60637, United States
- Materials Science Division, Argonne National Laboratory , Argonne, Illinois 60439, United States
| |
Collapse
|