1
|
O’Neill N, Shi BX, Fong K, Michaelides A, Schran C. To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water. J Phys Chem Lett 2024; 15:6081-6091. [PMID: 38820256 PMCID: PMC11181334 DOI: 10.1021/acs.jpclett.4c01030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 05/23/2024] [Accepted: 05/24/2024] [Indexed: 06/02/2024]
Abstract
The extent of ion pairing in solution is an important phenomenon to rationalize transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However, direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wave function based machine learning potentials with the random phase approximation (RPA) and second order Møller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.
Collapse
Affiliation(s)
- Niamh O’Neill
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Benjamin X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Kara Fong
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| | - Christoph Schran
- Cavendish
Laboratory, Department of Physics, University
of Cambridge, Cambridge CB3 0HE, United
Kingdom
- Lennard-Jones
Centre, University of Cambridge, Trinity Ln, Cambridge CB2 1TN, United Kingdom
| |
Collapse
|
2
|
Matin S, Allen AEA, Smith J, Lubbers N, Jadrich RB, Messerly R, Nebgen B, Li YW, Tretiak S, Barros K. Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment. J Chem Theory Comput 2024. [PMID: 38307009 DOI: 10.1021/acs.jctc.3c01051] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2024]
Abstract
Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. We investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.
Collapse
Affiliation(s)
- Sakib Matin
- Department of Physics, Boston University, Boston, Massachusetts 02215, United States
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Alice E A Allen
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Justin Smith
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- NVIDIA Corp., Santa Clara, California 95051, United States
| | - Nicholas Lubbers
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ryan B Jadrich
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Richard Messerly
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Benjamin Nebgen
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Ying Wai Li
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| | - Kipton Barros
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87546, United States
| |
Collapse
|
3
|
Yu Q, Qu C, Houston PL, Nandi A, Pandey P, Conte R, Bowman JM. A Status Report on "Gold Standard" Machine-Learned Potentials for Water. J Phys Chem Lett 2023; 14:8077-8087. [PMID: 37656898 PMCID: PMC10510435 DOI: 10.1021/acs.jpclett.3c01791] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 08/28/2023] [Indexed: 09/03/2023]
Abstract
Owing to the central importance of water to life as well as its unusual properties, potentials for water have been the subject of extensive research over the past 50 years. Recently, five potentials based on different machine learning approaches have been reported that are at or near the "gold standard" CCSD(T) level of theory. The development of such high-level potentials enables efficient and accurate simulations of water systems using classical and quantum dynamical approaches. This Perspective serves as a status report of these potentials, focusing on their methodology and applications to water systems across different phases. Their performances on the energies of gas phase water clusters, as well as condensed phase structural and dynamical properties, are discussed.
Collapse
Affiliation(s)
- Qi Yu
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Independent
Researcher, Toronto, Ontario M9B 0E3, Canada
| | - Paul L. Houston
- Department
of Chemistry and Chemical Biology, Cornell
University, Ithaca, New York 14853, United States
- Department of Chemistry
and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Apurba Nandi
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511, Luxembourg City, Luxembourg
| | - Priyanka Pandey
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento
di Chimica, Università degli Studi
di Milano, via Golgi 19, 20133 Milano, Italy
| | - Joel M. Bowman
- Department
of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
4
|
Cinq N, Simon A, Louisnard F, Cuny J. Accurate SCC-DFTB Parametrization of Liquid Water with Improved Atomic Charges and Iterative Boltzmann Inversion. J Phys Chem B 2023; 127:7590-7601. [PMID: 37603798 DOI: 10.1021/acs.jpcb.3c03479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
This work presents improvements of the description of liquid water within the self-consistent-charge density-functional based tight-binding scheme combining the use of Weighted Mulliken (WMull) charges and optimized O-H repulsive potential through the iterative Boltzmann inversion (IBI) process. The quality of the newly developed models is validated considering pair radial distribution functions (RDFs), as well as other structural, energetic, thermodynamic, and dynamic properties. The use of WMull charges certainly improves the agreement with experimental data, however leading to over-structured RDFs at short distance, that can be further improved by considering an optimized O-H repulsive potential obtained by the IBI process. Three different schemes were used to optimize this potential: (i) optimization including short O-H distances. This led to accurate RDFs as well as improved self-diffusion coefficient and heat of vaporization, while the proton transfer energy barrier is severely deteriorated; (ii) optimization starting at long distance. The proton transfer energy barrier is recovered while the heat of vaporization is deteriorated and the O-H RDF is less accurate at short distance; (iii) optimization within the path-integral molecular dynamics scheme which allows us to exclude nuclear quantum effects from the repulsive potential. The latter potential, in conjunction with the WMull improved atomic charges, provides similar results as (i) for structural, dynamic, and thermodynamic properties while recovering a large part of the proton transfer energy barrier. It therefore offers a good compromise to study both dynamic properties and chemistry within liquid water at a quantum chemical level.
Collapse
Affiliation(s)
- Nicolas Cinq
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Aude Simon
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Fernand Louisnard
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| | - Jérôme Cuny
- Laboratoire de Chimie et Physique Quantiques (LCPQ), FeRMI Institute, Université de Toulouse [UT3] and CNRS, Toulouse F-31062, France
| |
Collapse
|
5
|
Chen MS, Lee J, Ye HZ, Berkelbach TC, Reichman DR, Markland TE. Data-Efficient Machine Learning Potentials from Transfer Learning of Periodic Correlated Electronic Structure Methods: Liquid Water at AFQMC, CCSD, and CCSD(T) Accuracy. J Chem Theory Comput 2023; 19:4510-4519. [PMID: 36730728 DOI: 10.1021/acs.jctc.2c01203] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Obtaining the atomistic structure and dynamics of disordered condensed-phase systems from first-principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure and provide a data-efficient approach to obtain machine-learned condensed-phase potential energy surfaces using AFQMC, CCSD, and CCSD(T) from a very small number (≤200) of energies by leveraging a transfer learning scheme starting from lower-tier electronic structure methods. We demonstrate the effectiveness of this approach for liquid water by performing both classical and path integral molecular dynamics simulations on these machine-learned potential energy surfaces. By doing this, we uncover the interplay of dynamical electron correlation and nuclear quantum effects across the entire liquid range of water while providing a general strategy for efficiently utilizing periodic correlated electronic structure methods to explore disordered condensed-phase systems.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York10027, United States
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York10010, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| |
Collapse
|
6
|
Huo J, Chen J, Liu P, Hong B, Zhang J, Dong H, Li S. Microscopic Mechanism of Proton Transfer in Pure Water under Ambient Conditions. J Chem Theory Comput 2023. [PMID: 37365994 DOI: 10.1021/acs.jctc.3c00244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
Water molecules and the associated proton transfer (PT) are prevalent in chemical and biological systems and have been a hot research topic. Spectroscopic characterization and ab initio molecular dynamics (AIMD) simulations have previously revealed insights into acidic and basic liquids. Presumably, the situation in the acidic/basic solution is not necessarily the same as in pure water; in addition, the autoionization constant for water is only 10-14 under ambient conditions, making the study of PT in pure water challenging. To overcome this issue, we modeled periodic water box systems containing 1000 molecules for tens of nanoseconds based on a neural network potential (NNP) with quantum mechanical accuracy. The NNP was generated by training a dataset containing the energies and atomic forces of 17 075 configurations of periodic water box systems, and these data points were calculated at the MP2 level that considers electron correlation effects. We found that the size of the system and the duration of the simulation have a significant impact on the convergence of the results. With these factors considered, our simulations showed that hydronium (H3O+) and hydroxide (OH-) ions in water have distinct hydration structures, thermodynamic and kinetic properties, e.g., the longer-lasting and more stable hydrated structure of OH- ions than that of H3O+, as well as a significantly higher free energy barrier for the OH--associated PT than that of H3O+, leading the two to exhibit completely different PT behaviors. Given these characteristics, we further found that PT via OH- ions tends not to occur multiple times or between many molecules. In contrast, PT via H3O+ can synergistically occur among multiple molecules and prefers to adopt a cyclic pattern among three water molecules, while it occurs mostly in a chain pattern when more water molecules are involved. Therefore, our studies provide a detailed and solid microscopic explanation for the PT process in pure water.
Collapse
Affiliation(s)
- Jun Huo
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
| | - Jianghao Chen
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Pei Liu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Benkun Hong
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| | - Jian Zhang
- School of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Hao Dong
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
- State Key Laboratory of Analytical Chemistry for Life Science, Nanjing University, Nanjing 210023, China
- Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, China
| |
Collapse
|
7
|
Wróbel P, Kubisiak P, Eilmes A. Hydrogen Bonding and Infrared Spectra of Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide/Water Mixtures: A View from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:10922-10932. [PMID: 36516319 PMCID: PMC9806834 DOI: 10.1021/acs.jpcb.2c06947] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Simulations of ab initio molecular dynamics have been performed for mixtures of ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) ionic liquid and water. Statistics of donors and acceptors of hydrogen bonds has revealed that with increasing water content, hydrogen bonds between EMIM cations and TFSI anions are replaced by bonds to water molecules. In the mixture of liquids, the total number of bonds (from EMIM cations or water molecules) formed by TFSI acceptors increases. IR spectra obtained from ab initio molecular dynamics trajectories are in good agreement with literature data for ionic liquid/water systems. Analysis of oscillations of individual C-H and O-H bonds has shown correlations between vibrational frequencies and hydrogen bonds formed by an EMIM cation or water molecule and has indicated that the changes in the IR spectrum result from the decreased number of water-water hydrogen bonds in the mixture. The tests of DFTB methodology with tailored parameterizations have yielded reasonably good description of the IR spectrum of bulk water, whereas available parameterizations have failed in satisfactory reproduction of the IR spectrum of EMIM-TFSI/water mixtures in the region above 3000 cm-1.
Collapse
|
8
|
Liu J, He X. Recent advances in quantum fragmentation approaches to complex molecular and condensed‐phase systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy China Pharmaceutical University Nanjing China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering East China Normal University Shanghai China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering East China Normal University Shanghai China
- New York University‐East China Normal University Center for Computational Chemistry New York University Shanghai Shanghai China
| |
Collapse
|