1
|
Li M, Wang P, Yu X, Su Y, Zhao J. Impact of Nuclear Quantum Effects on the Structural Properties of Protonated Water Clusters. J Phys Chem A 2024; 128:5954-5962. [PMID: 39007820 DOI: 10.1021/acs.jpca.4c03340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
Nuclear quantum effects (NQEs) play a crucial role in hydrogen-bonded systems due to quantum tunneling and proton fluctuation. Our understanding of how NQEs affect microstructures mainly focuses on bulk phases of liquids and solids but remains deficient for water clusters, including their hydrogen nuclei, hydrogen-bonded configurations, and temperature dependence. Here, we conducted ab initio molecular dynamics (MD) and path integral MD simulations to investigate the influence of NQEs on the structural properties of protonated water clusters H+(H2O)n (n = 3, 6, 9, 12). The results reveal that the NQEs become less evident as the cluster size increases due to the competition between NQEs and electrostatic interactions. Simulations of several H+(H2O)6 isomers at different temperatures indicate that the effect of elevated temperature on proton transfer is related to the initial structure. Interestingly, the process of proton transfer also involves the interconversion between Zundel-type and Eigen-type isomers. These findings significantly deepen our understanding of ion-water and water-water interactions, opening new avenues for the study of hydrated ion clusters and related systems.
Collapse
Affiliation(s)
- Mengxu Li
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | | | - Xueke Yu
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Yan Su
- Key Laboratory of Materials Modification by Laser, Ion and Electron Beams (Dalian University of Technology), Ministry of Education, Dalian 116024, China
| | - Jijun Zhao
- Guangdong Basic Research Center of Excellence for Structure and Fundamental Interactions of Matter, Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, School of Physics, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
2
|
Ojha D, Henao A, Zysk F, Kühne TD. Nuclear quantum effects on the vibrational dynamics of the water-air interface. J Chem Phys 2024; 160:204114. [PMID: 38804494 DOI: 10.1063/5.0204071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 04/11/2024] [Indexed: 05/29/2024] Open
Abstract
We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.
Collapse
Affiliation(s)
- Deepak Ojha
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Andrés Henao
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Frederik Zysk
- Dynamics of Condensed Matter and Center for Sustainable Systems Design, Department of Chemistry, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
| | - Thomas D Kühne
- Center for Advanced Systems Understanding (CASUS), Untermarkt 20, D-02826 Görlitz, Germany, Helmholtz Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, D-01328 Dresden, Germany, and TU Dresden, Institute of Artificial Intelligence, Chair of Computational System Sciences, Nöthnitzer Straße 46, D-01187 Dresden, Germany
| |
Collapse
|
3
|
Thakur AC, Remsing RC. Nuclear quantum effects in the acetylene:ammonia plastic co-crystal. J Chem Phys 2024; 160:024502. [PMID: 38189604 DOI: 10.1063/5.0179161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 12/20/2023] [Indexed: 01/09/2024] Open
Abstract
Organic molecular solids can exhibit rich phase diagrams. In addition to structurally unique phases, translational and rotational degrees of freedom can melt at different state points, giving rise to partially disordered solid phases. The structural and dynamic disorder in these materials can have a significant impact on the physical properties of the organic solid, necessitating a thorough understanding of disorder at the atomic scale. When these disordered phases form at low temperatures, especially in crystals with light nuclei, the prediction of material properties can be complicated by the importance of nuclear quantum effects. As an example, we investigate nuclear quantum effects on the structure and dynamics of the orientationally disordered, translationally ordered plastic phase of the acetylene:ammonia (1:1) co-crystal that is expected to exist on the surface of Saturn's moon Titan. Titan's low surface temperature (∼90 K) suggests that the quantum mechanical behavior of nuclei may be important in this and other molecular solids in these environments. By using neural network potentials combined with ring polymer molecular dynamics simulations, we show that nuclear quantum effects increase orientational disorder and rotational dynamics within the acetylene:ammonia (1:1) co-crystal by weakening hydrogen bonds. Our results suggest that nuclear quantum effects are important to accurately model molecular solids and their physical properties in low-temperature environments.
Collapse
Affiliation(s)
- Atul C Thakur
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Richard C Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
4
|
Sharma A, Kumar V, Chakraborty S. Micro-Solvation of Propofol in Propylene Glycol-Water Binary Mixtures: Molecular Dynamics Simulation Studies. J Phys Chem B 2023; 127:11011-11022. [PMID: 37972382 DOI: 10.1021/acs.jpcb.3c04932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
The water microstructure around propofol plays a crucial role in controlling their solubility in the binary mixture. The unusual nature of such a water microstructure can influence both translational and reorientational dynamics, as well as the water hydrogen bond network near propofol. We have carried out all-atom molecular dynamics simulations of five different compositions of the propylene glycol (PG)/water binary mixture containing propofol (PFL) molecules to investigate the differential behavior of water microsolvation shells around propofol, which is likely to control the propofol solubility. It is evident from the simulation snapshots for various compositions that the PG at high molecular ratio favors the water cluster and extended chainlike network that percolates within the PG matrix, where the propofol is in the dispersed state. We estimated that the radial distribution function indicates higher ordered water microstructure around propofol for high PG content, as compared to the lower PG content in the PG/water mixture. So, the hydrophilic PG regulates the stability of the water micronetwork around propofol and its solubility in the binary mixture. We observed that the translational and rotational mobility of water belonging to the propofol microsolvation shell is hindered for high PG content and relaxed toward the low PG molecular ratio in the PG/water mixture. It has been noticed that the structural relaxation of the hydrogen bond formed between the propofol and the water molecules present in the propofol microsolvation shell for all five compositions is found to be slower for high PG content and becomes faster on the way to low PG content in the mixture. Simultaneously, we calculated the intermittent residence time correlation function of the water molecules belonging to the microsolvation shell around the propofol for five different compositions and found a faster short time decay followed up with long time components. Again, the origin of such long time decay is primarily from the structural relaxation of the microsolvation shell around the propofol, where the high PG content shows the slower structural relaxation that turns faster as the PG content approaches to the other end of the compositions. So, our studies showed that the slower structural relaxation of the microsolvation shell around propofol for a high PG molecular ratio in the PG/water mixture correlate well with the extensive ordering of the water microstructure and restricted water mobility and facilitates the dissolution process of propofol in the binary mixture.
Collapse
Affiliation(s)
- Anupama Sharma
- Department of Computational Sciences, School of Basic Sciences, Central University of Punjab, Bathinda 151401, India
| | - Vishal Kumar
- Department of Computational Sciences, School of Basic Sciences, Central University of Punjab, Bathinda 151401, India
| | - Sudip Chakraborty
- Department of Computational Sciences, School of Basic Sciences, Central University of Punjab, Bathinda 151401, India
| |
Collapse
|
5
|
Heindel JP, Herman KM, Xantheas SS. Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development. Annu Rev Phys Chem 2023; 74:337-360. [PMID: 37093659 DOI: 10.1146/annurev-physchem-062422-023532] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.
Collapse
Affiliation(s)
- Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington, USA; ,
| |
Collapse
|
6
|
Plé T, Mauger N, Adjoua O, Inizan TJ, Lagardère L, Huppert S, Piquemal JP. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J Chem Theory Comput 2023; 19:1432-1445. [PMID: 36856658 DOI: 10.1021/acs.jctc.2c01233] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, and parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostatting. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will help to shed new light on the quantum nature of important interactions in biological matter.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | | | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Simon Huppert
- Institut des Nanosciences de Paris (INSP), CNRS UMR 7588, and Sorbonne Université, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France.,Institut Universitaire de France, 75005 Paris, France.,Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
7
|
Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. Improving Condensed-Phase Water Dynamics with Explicit Nuclear Quantum Effects: The Polarizable Q-AMOEBA Force Field. J Phys Chem B 2022; 126:8813-8826. [PMID: 36270033 DOI: 10.1021/acs.jpcb.2c04454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We introduce a new parametrization of the AMOEBA polarizable force field for water denoted Q-AMOEBA, for use in simulations that explicitly account for nuclear quantum effects (NQEs). This study is made possible thanks to the recently introduced adaptive Quantum Thermal Bath (adQTB) simulation technique which computational cost is comparable to classical molecular dynamics. The flexible Q-AMOEBA model conserves the initial AMOEBA functional form, with an intermolecular potential including an atomic multipole description of electrostatic interactions (up to quadrupole), a polarization contribution based on the Thole interaction model and a buffered 14-7 potential to model van der Waals interactions. It has been obtained by using a ForceBalance fitting strategy including high-level quantum chemistry reference energies and selected condensed-phase properties targets. The final Q-AMOEBA model is shown to accurately reproduce both gas-phase and condensed-phase properties, notably improving the original AMOEBA water model. This development allows the fine study of NQEs on water liquid phase properties such as the average H-O-H angle compared to its gas-phase equilibrium value, isotope effects, and so on. Q-AMOEBA also provides improved infrared spectroscopy prediction capabilities compared to AMOEBA03. Overall, we show that the impact of NQEs depends on the underlying model functional form and on the associated strength of hydrogen bonds. Since adQTB simulations can be performed at near classical computational cost using the Tinker-HP package, Q-AMOEBA can be extended to organic molecules, proteins, and nucleic acids opening the possibility for the large-scale study of the importance of NQEs in biophysics.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 75005 Paris, France
| |
Collapse
|
8
|
Shaimardanov AR, Shulga DA, Palyulin VA. Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields? J Phys Chem A 2022; 126:6278-6294. [PMID: 36054931 DOI: 10.1021/acs.jpca.2c02722] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.
Collapse
Affiliation(s)
- Arslan R Shaimardanov
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russian Federation
| | - Dmitry A Shulga
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russian Federation
| | - Vladimir A Palyulin
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russian Federation
| |
Collapse
|
9
|
Naseem-Khan S, Lagardère L, Narth C, Cisneros GA, Ren P, Gresh N, Piquemal JP. Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:3607-3621. [PMID: 35575306 PMCID: PMC10851344 DOI: 10.1021/acs.jctc.2c00029] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We present the extension of the Sum of Interactions Between Fragments Ab initio Computed (SIBFA) many-body polarizable force field to condensed-phase molecular dynamics (MD) simulations. The quantum-inspired SIBFA procedure is grounded on simplified integrals obtained from localized molecular orbital theory and achieves full separability of its intermolecular potential. It embodies long-range multipolar electrostatics (up to quadrupole) coupled to a short-range penetration correction (up to charge-quadrupole), exchange repulsion, many-body polarization, many-body charge transfer/delocalization, exchange dispersion, and dispersion (up to C10). This enables the reproduction of all energy contributions of ab initio symmetry-adapted perturbation theory (SAPT(DFT)) gas-phase reference computations. The SIBFA approach has been integrated within the Tinker-HP massively parallel MD package. To do so, all SIBFA energy gradients have been derived and the approach has been extended to enable periodic boundary conditions simulations using smooth particle mesh Ewald. This novel implementation also notably includes a computationally tractable simplification of the many-body charge transfer/delocalization contribution. As a proof of concept, we perform a first computational experiment defining a water model fitted on a limited set of SAPT(DFT) data. SIBFA is shown to enable a satisfactory reproduction of both gas-phase energetic contributions and condensed-phase properties highlighting the importance of its physically motivated functional form.
Collapse
Affiliation(s)
- Sehr Naseem-Khan
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Louis Lagardère
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- IP2CT, FR 2622, CNRS, Sorbonne Université, 75005 Paris, France
| | | | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Nohad Gresh
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
| | - Jean-Philip Piquemal
- LCT, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
10
|
Li C, Paesani F, Voth GA. Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature. J Chem Theory Comput 2022; 18:2124-2131. [PMID: 35263110 PMCID: PMC9059465 DOI: 10.1021/acs.jctc.1c01223] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. We find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.
Collapse
Affiliation(s)
- Chenghan Li
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, 60637
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, Materials Science and Engineering, and San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, 60637
| |
Collapse
|
11
|
Naseem-Khan S, Piquemal JP, Cisneros GA. Improvement of the Gaussian Electrostatic Model by separate fitting of Coulomb and exchange-repulsion densities and implementation of a new dispersion term. J Chem Phys 2021; 155:194103. [PMID: 34800949 PMCID: PMC8598263 DOI: 10.1063/5.0072380] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Accepted: 10/28/2021] [Indexed: 11/15/2022] Open
Abstract
The description of each separable contribution of the intermolecular interaction is a useful approach to develop polarizable force fields (polFFs). The Gaussian Electrostatic Model (GEM) is based on this approach, coupled with the use of density fitting techniques. In this work, we present the implementation and testing of two improvements of GEM: the Coulomb and exchange-repulsion energies are now computed with separate frozen molecular densities and a new dispersion formulation inspired by the Sum of Interactions Between Fragments Ab initio Computed polFF, which has been implemented to describe the dispersion and charge-transfer interactions. Thanks to the combination of GEM characteristics and these new features, we demonstrate a better agreement of the computed structural and condensed properties for water with experimental results, as well as binding energies in the gas phase with the ab initio reference compared with the previous GEM* potential. This work provides further improvements to GEM and the items that remain to be improved and the importance of the accurate reproduction for each separate contribution.
Collapse
Affiliation(s)
- Sehr Naseem-Khan
- Department of Chemistry, University of North Texas, Denton, Texas 76201, USA
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR 7616 CNRS, 75005 Paris, France
| | - G. Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, USA
| |
Collapse
|
12
|
Dubey V, Dueby S, Daschakraborty S. Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective. Phys Chem Chem Phys 2021; 23:19964-19986. [PMID: 34515269 DOI: 10.1039/d1cp02202d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | - Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
13
|
Kutus B, Shalit A, Hamm P, Hunger J. Dielectric response of light, heavy and heavy-oxygen water: isotope effects on the hydrogen-bonding network's collective relaxation dynamics. Phys Chem Chem Phys 2021; 23:5467-5473. [PMID: 33650591 DOI: 10.1039/d0cp06460b] [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
Isotopic substitutions largely affect the dielectric relaxation dynamics of hydrogen-bonded liquid water; yet, the role of the altered molecular masses and nuclear quantum effects has not been fully established. To disentangle these two effects we study the dielectric relaxation of light (H216O), heavy (D216O) and heavy-oxygen (H218O) water at temperatures ranging from 278 to 338 K. Upon 16O/18O exchange, we find that the relaxation time of the collective orientational relaxation mode of water increases by 4-5%, in quantitative agreement with the enhancement of viscosity. Despite the rotational character of dielectric relaxation, the increase is consistent with a translational mass factor. For H/D substitution, the slow-down of the relaxation time is more pronounced and also shows a strong temperature dependence. In addition to the classical mass factor, the enhancement of the relaxation time for D216O can be described by an apparent temperature shift of 7.2 K relative to H216O, which is higher than the 6.5 K shift reported for viscosity. As this shift accounts for altered zero-point energies, the comparison suggests that the underlying thermally populated states relevant to the activation of viscous flow and dielectric relaxation differ.
Collapse
Affiliation(s)
- Bence Kutus
- Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | - Andrey Shalit
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Johannes Hunger
- Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| |
Collapse
|
14
|
Agosta L, Dzugutov M, Hermansson K. Supercooled liquid-like dynamics in water near a fully hydrated titania surface: Decoupling of rotational and translational diffusion. J Chem Phys 2021; 154:094708. [PMID: 33685161 DOI: 10.1063/5.0039693] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO2 on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.
Collapse
Affiliation(s)
- Lorenzo Agosta
- Department of Chemistry-Ångström, Uppsala University, S-75121 Uppsala, Sweden
| | - Mikhail Dzugutov
- Department of Chemistry-Ångström, Uppsala University, S-75121 Uppsala, Sweden
| | - Kersti Hermansson
- Department of Chemistry-Ångström, Uppsala University, S-75121 Uppsala, Sweden
| |
Collapse
|
15
|
Mao Y, Loipersberger M, Horn PR, Das A, Demerdash O, Levine DS, Prasad Veccham S, Head-Gordon T, Head-Gordon M. From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis. Annu Rev Phys Chem 2021; 72:641-666. [PMID: 33636998 DOI: 10.1146/annurev-physchem-090419-115149] [Citation(s) in RCA: 47] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N2, and BF bound to a [Ru(II)(NH3)5]2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.
Collapse
Affiliation(s)
- Yuezhi Mao
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Matthias Loipersberger
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Paul R Horn
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Akshaya Das
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Omar Demerdash
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Daniel S Levine
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Srimukh Prasad Veccham
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Teresa Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| |
Collapse
|
16
|
Liquid water contains the building blocks of diverse ice phases. Nat Commun 2020; 11:5757. [PMID: 33188195 PMCID: PMC7666157 DOI: 10.1038/s41467-020-19606-y] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 10/23/2020] [Indexed: 11/24/2022] Open
Abstract
Water molecules can arrange into a liquid with complex hydrogen-bond networks and at least 17 experimentally confirmed ice phases with enormous structural diversity. It remains a puzzle how or whether this multitude of arrangements in different phases of water are related. Here we investigate the structural similarities between liquid water and a comprehensive set of 54 ice phases in simulations, by directly comparing their local environments using general atomic descriptors, and also by demonstrating that a machine-learning potential trained on liquid water alone can predict the densities, lattice energies, and vibrational properties of the ices. The finding that the local environments characterising the different ice phases are found in water sheds light on the phase behavior of water, and rationalizes the transferability of water models between different phases. Molecular understanding of water is challenging due to the structural complexity of liquid water and the large number of ice phases. Here the authors use a machine-learning potential trained on liquid water to demonstrate the structural similarity of liquid water and that of 54 real and hypothetical ice phases.
Collapse
|
17
|
Lambros E, Paesani F. How good are polarizable and flexible models for water: Insights from a many-body perspective. J Chem Phys 2020; 153:060901. [DOI: 10.1063/5.0017590] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
18
|
Yao Y, Kanai Y. Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional. J Chem Phys 2020; 153:044114. [DOI: 10.1063/5.0012815] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
19
|
Bilbrey JA, Heindel JP, Schram M, Bandyopadhyay P, Xantheas SS, Choudhury S. A look inside the black box: Using graph-theoretical descriptors to interpret a Continuous-Filter Convolutional Neural Network (CF-CNN) trained on the global and local minimum energy structures of neutral water clusters. J Chem Phys 2020; 153:024302. [DOI: 10.1063/5.0009933] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jenna A. Bilbrey
- Computing and Analytics Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Joseph P. Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Malachi Schram
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Sotiris S. Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| | - Sutanay Choudhury
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, USA
| |
Collapse
|
20
|
Berta D, Ferenc D, Bakó I, Madarász Á. Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water. J Chem Theory Comput 2020; 16:3316-3334. [PMID: 32268067 PMCID: PMC7304866 DOI: 10.1021/acs.jctc.9b00703] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
![]()
Nuclear quantum effects
have significant contributions to thermodynamic
quantities and structural properties; furthermore, very expensive
methods are necessary for their accurate computation. In most calculations,
these effects, for instance, zero-point energies, are simply neglected
or only taken into account within the quantum harmonic oscillator
approximation. Herein, we present a new method, Generalized Smoothed
Trajectory Analysis, to determine nuclear quantum effects from molecular
dynamics simulations. The broad applicability is demonstrated with
the examples of a harmonic oscillator and different states of water.
Ab initio molecular dynamics simulations have been performed for ideal
gas up to the temperature of 5000 K. Classical molecular dynamics
have been carried out for hexagonal ice, liquid water, and vapor at
atmospheric pressure. With respect to the experimental heat capacity,
our method outperforms previous calculations in the literature in
a wide temperature range at lower computational cost than other alternatives.
Dynamic and structural nuclear quantum effects of water are also discussed.
Collapse
Affiliation(s)
- Dénes Berta
- Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar Tudósok Körútja 2, H-1117 Budapest, Hungary.,Department of Chemistry, Kings College London, 7 Trinity Street, SE1 1DB London, United Kingdom
| | - Dávid Ferenc
- Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar Tudósok Körútja 2, H-1117 Budapest, Hungary
| | - Imre Bakó
- Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar Tudósok Körútja 2, H-1117 Budapest, Hungary
| | - Ádám Madarász
- Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar Tudósok Körútja 2, H-1117 Budapest, Hungary
| |
Collapse
|
21
|
Rakshit A, Bandyopadhyay P, Heindel JP, Xantheas SS. Atlas of putative minima and low-lying energy networks of water clusters n = 3-25. J Chem Phys 2019; 151:214307. [PMID: 31822087 DOI: 10.1063/1.5128378] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a database consisting of the putative minima and ∼3.2 × 106 local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.
Collapse
Affiliation(s)
- Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
22
|
Yong X, Burnham CJ, English NJ, Tse JS. Classical and path-integral molecular-dynamics study on liquid water and ice melting using non-empirical TTM2.1-F model. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652774] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Xue Yong
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| | - Christian J. Burnham
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Dublin, Ireland
| | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Canada
| |
Collapse
|
23
|
Carnevale D, Pelupessy P, Bodenhausen G. Cross-term Splittings Due to the Orientational Inequivalence of Proton Magnetic Shielding Tensors: Do Water Molecules Trapped in Crystals Hop or Tunnel? J Phys Chem Lett 2019; 10:3224-3231. [PMID: 31059264 DOI: 10.1021/acs.jpclett.9b00914] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Water molecules trapped in crystals of barium chlorate monohydrate have been investigated by magic-angle spinning (MAS) proton NMR spectroscopy in the temperature range 110-300 K. At high temperatures, a single spinning sideband pattern is observed. Below 150 K, however, two interleaved spinning sideband manifolds appear, with distinct centerbands that do not coincide with the average isotropic chemical shift seen at high temperatures. This hitherto unknown "cross-term splitting" results from the interplay of the homonuclear dipole-dipole coupling and two anisotropic proton shielding tensors that have identical principal components but nonequivalent orientations. The resulting cross terms cannot be averaged out by rotation about the magic angle. The analysis of the exchange-induced broadening, coalescence, and narrowing of the cross-term splitting in MAS spectra allows one to estimate the rate of exchange of the two protons between 140 and 190 K. The experimental data is compared with 2H and 1H NMR studies of the same sample reported in the literature. Density functional theory methods are utilized to estimate the thermal activation energy for a 2-fold hopping process of proton exchange about the H-O-H bisector. The Bell-Limbach model allows one to take into account contributions due to incoherent quantum tunneling in the low-temperature regime.
Collapse
Affiliation(s)
- Diego Carnevale
- Laboratoire des Biomolécules, LBM, Département de Chimie , École Normale Supérieure , PSL University, Sorbonne Université, CNRS, 75005 Paris , France
| | - Philippe Pelupessy
- Laboratoire des Biomolécules, LBM, Département de Chimie , École Normale Supérieure , PSL University, Sorbonne Université, CNRS, 75005 Paris , France
| | - Geoffrey Bodenhausen
- Laboratoire des Biomolécules, LBM, Département de Chimie , École Normale Supérieure , PSL University, Sorbonne Université, CNRS, 75005 Paris , France
| |
Collapse
|
24
|
Berger A, Ciardi G, Sidler D, Hamm P, Shalit A. Impact of nuclear quantum effects on the structural inhomogeneity of liquid water. Proc Natl Acad Sci U S A 2019; 116:2458-2463. [PMID: 30692247 PMCID: PMC6377483 DOI: 10.1073/pnas.1818182116] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The 2D Raman-terahertz (THz) response of liquid water is studied in dependence of temperature and isotope substitution ([Formula: see text]O, [Formula: see text]O, and [Formula: see text]O). In either case, a very short-lived (i.e., between 75 and 95 fs) echo is observed that reports on the inhomogeneity of the low-frequency intermolecular modes and hence, on the heterogeneity of the hydrogen bond networks of water. The echo lifetime slows down by about 20% when cooling the liquid from room temperature to the freezing point. Furthermore, the echo lifetime of [Formula: see text]O is [Formula: see text] slower than that of [Formula: see text]O, and both can be mapped on each other by introducing an effective temperature shift of [Formula: see text] K. In contrast, the temperature-dependent echo lifetimes of [Formula: see text]O and [Formula: see text]O are the same within error. [Formula: see text]O and [Formula: see text]O have identical masses, yet [Formula: see text]O is much closer to [Formula: see text]O in terms of nuclear quantum effects. It is, therefore, concluded that the echo is a measure of the structural inhomogeneity of liquid water induced by nuclear quantum effects.
Collapse
Affiliation(s)
- Arian Berger
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Gustavo Ciardi
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - David Sidler
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Andrey Shalit
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| |
Collapse
|
25
|
Silletta EV, Tuckerman ME, Jerschow A. Unusual Proton Transfer Kinetics in Water at the Temperature of Maximum Density. PHYSICAL REVIEW LETTERS 2018; 121:076001. [PMID: 30169046 DOI: 10.1103/physrevlett.121.076001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Revised: 06/26/2018] [Indexed: 06/08/2023]
Abstract
Water exhibits numerous anomalous properties, many of which remain poorly understood. One of its intriguing behaviors is that it exhibits a temperature of maximum density (TMD) at 4 °C. We provide here new experimental evidence for hitherto unknown abrupt changes in proton transfer kinetics at the TMD. In particular, we show that the lifetime of OH^{-} ions has a maximum at this temperature, in contrast to hydronium ions. Furthermore, base-catalyzed proton transfer shows a sharp local minimum at this temperature, and activation energies change abruptly as well. The measured lifetimes agree with earlier theoretical predictions as the temperature approaches the TMD. Similar results are also found for heavy water at its own TMD. These findings point to a high propensity of forming fourfold coordinated OH^{-} solvation complexes at the TMD, underlining the asymmetry between hydroxide and hydronium transport. These results could help to further elucidate the unusual properties of water and related liquids.
Collapse
Affiliation(s)
- Emilia V Silletta
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Mark E Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, USA
- Courant Institute of Mathematical Science, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| | - Alexej Jerschow
- Department of Chemistry, New York University, New York, New York 10003, USA
| |
Collapse
|
26
|
Han Y, Jin J, Wagner JW, Voth GA. Quantum theory of multiscale coarse-graining. J Chem Phys 2018; 148:102335. [PMID: 29544317 DOI: 10.1063/1.5010270] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Coarse-grained (CG) models serve as a powerful tool to simulate molecular systems at much longer temporal and spatial scales. Previously, CG models and methods have been built upon classical statistical mechanics. The present paper develops a theory and numerical methodology for coarse-graining in quantum statistical mechanics, by generalizing the multiscale coarse-graining (MS-CG) method to quantum Boltzmann statistics. A rigorous derivation of the sufficient thermodynamic consistency condition is first presented via imaginary time Feynman path integrals. It identifies the optimal choice of CG action functional and effective quantum CG (qCG) force field to generate a quantum MS-CG (qMS-CG) description of the equilibrium system that is consistent with the quantum fine-grained model projected onto the CG variables. A variational principle then provides a class of algorithms for optimally approximating the qMS-CG force fields. Specifically, a variational method based on force matching, which was also adopted in the classical MS-CG theory, is generalized to quantum Boltzmann statistics. The qMS-CG numerical algorithms and practical issues in implementing this variational minimization procedure are also discussed. Then, two numerical examples are presented to demonstrate the method. Finally, as an alternative strategy, a quasi-classical approximation for the thermal density matrix expressed in the CG variables is derived. This approach provides an interesting physical picture for coarse-graining in quantum Boltzmann statistical mechanics in which the consistency with the quantum particle delocalization is obviously manifest, and it opens up an avenue for using path integral centroid-based effective classical force fields in a coarse-graining methodology.
Collapse
Affiliation(s)
- Yining Han
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jaehyeok Jin
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jacob W Wagner
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
27
|
|
28
|
Poltavsky I, DiStasio RA, Tkatchenko A. Perturbed path integrals in imaginary time: Efficiently modeling nuclear quantum effects in molecules and materials. J Chem Phys 2018; 148:102325. [DOI: 10.1063/1.5006596] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- Igor Poltavsky
- Physics and Materials Science Research Unit, University of Luxembourg, Luxembourg L-1511, Luxembourg
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, Luxembourg L-1511, Luxembourg
| |
Collapse
|
29
|
Feng Y, Wang Z, Guo J, Chen J, Wang EG, Jiang Y, Li XZ. The collective and quantum nature of proton transfer in the cyclic water tetramer on NaCl(001). J Chem Phys 2018; 148:102329. [PMID: 29544296 DOI: 10.1063/1.5004737] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Proton tunneling is an elementary process in the dynamics of hydrogen-bonded systems. Collective tunneling is known to exist for a long time. Atomistic investigations of this mechanism in realistic systems, however, are scarce. Using a combination of ab initio theoretical and high-resolution experimental methods, we investigate the role played by the protons on the chirality switching of a water tetramer on NaCl(001). Our scanning tunneling spectroscopies show that partial deuteration of the H2O tetramer with only one D2O leads to a significant suppression of the chirality switching rate at a cryogenic temperature (T), indicating that the chirality switches by tunneling in a concerted manner. Theoretical simulations, in the meantime, support this picture by presenting a much smaller free-energy barrier for the translational collective proton tunneling mode than other chirality switching modes at low T. During this analysis, the virial energy provides a reasonable estimator for the description of the nuclear quantum effects when a traditional thermodynamic integration method cannot be used, which could be employed in future studies of similar problems. Given the high-dimensional nature of realistic systems and the topology of the hydrogen-bonded network, collective proton tunneling may exist more ubiquitously than expected. Systems of this kind can serve as ideal platforms for studies of this mechanism, easily accessible to high-resolution experimental measurements.
Collapse
Affiliation(s)
- Yexin Feng
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Zhichang Wang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Jing Guo
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ji Chen
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - En-Ge Wang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ying Jiang
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Xin-Zheng Li
- International Center for Quantum Materials and School of Physics, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
30
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
31
|
Reddy SK, Moberg DR, Straight SC, Paesani F. Temperature-dependent vibrational spectra and structure of liquid water from classical and quantum simulations with the MB-pol potential energy function. J Chem Phys 2017; 147:244504. [DOI: 10.1063/1.5006480] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Sandeep K. Reddy
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Daniel R. Moberg
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, Materials Science and Engineering, and San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
32
|
Demerdash O, Wang L, Head‐Gordon T. Advanced models for water simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1355] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Omar Demerdash
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California, Davis Davis CA USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry University of California Berkeley CA USA
- Department of Chemistry University of California Berkeley CA USA
- Department of Bioengineering University of California Berkeley CA USA
- Department of Chemical and Biomolecular Engineering University of California Berkeley CA USA
| |
Collapse
|
33
|
Yong X, Tse JS, English NJ. optPBE-vdW density functional theory study of liquid water and pressure-induced structural evolution in ice Ih. CAN J CHEM 2017. [DOI: 10.1139/cjc-2017-0201] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The accuracy of several local and non-local van der Waals (vdW) corrected exchange correlation functionals on the description of the effect of pressure on ice has been investigated. In a preliminary survey, the non-local vdW correction used in conjunction with the optPBE functional was shown to provide the best overall agreement on the structural parameters of ice Ih with experiments. More importantly, this combination reproduced correctly the recently observed crystal → crystal transformation in ice Ih at 80 K prior to amorphisation. The predicted transition pressure of 1.9 GPa is somewhat higher, showing that the current generation of vdW functionals are still not sufficiently accurate for the ice system. The existence of an intermediate crystalline state with a shear-hexagonal structure confirms the earlier prediction that the collapse of crystalline structure under compression originates from the softening of phonon modes in ice Ih’s basal plane.
Collapse
Affiliation(s)
- Xue Yong
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SK S7N 5E2, Canada
| | - John S. Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SK S7N 5E2, Canada
| | - Niall J. English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
34
|
Hamm P, Fanourgakis GS, Xantheas SS. A surprisingly simple correlation between the classical and quantum structural networks in liquid water. J Chem Phys 2017; 147:064506. [DOI: 10.1063/1.4993166] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Peter Hamm
- Department of Chemistry, University of Zurich, Zurich 8057, Switzerland
| | - George S. Fanourgakis
- Environmental Chemical Processes Laboratory, Department of Chemistry, University of Crete, P.O. Box 2208, 70013 Heraklion, Greece
| | - Sotiris S. Xantheas
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, USA
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
35
|
Wilkins DM, Manolopoulos DE, Pipolo S, Laage D, Hynes JT. Nuclear Quantum Effects in Water Reorientation and Hydrogen-Bond Dynamics. J Phys Chem Lett 2017; 8:2602-2607. [PMID: 28530836 DOI: 10.1021/acs.jpclett.7b00979] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We combine classical and ring polymer molecular dynamics simulations with the molecular jump model to provide a molecular description of the nuclear quantum effects (NQEs) on water reorientation and hydrogen-bond dynamics in liquid H2O and D2O. We show that while the net NQE is negligible in D2O, it leads to a ∼13% acceleration in H2O dynamics compared to a classical description. Large angular jumps-exchanging hydrogen-bond partners-are the dominant reorientation pathway (just as in a classical description); the faster reorientation dynamics arise from the increased jump rate constant. NQEs do not change the jump amplitude distribution, and no significant tunneling is found. The faster jump dynamics are quantitatively related to decreased structuring of the OO radial distribution function when NQEs are included. This is explained, via a jump model analysis, by competition between the effects of water's librational and OH stretch mode zero-point energies on the hydrogen-bond strength.
Collapse
Affiliation(s)
- David M Wilkins
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne , 1015 Lausanne, Switzerland
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, University of Oxford , South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Silvio Pipolo
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ Paris 06, CNRS, PSL Research University , 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR , 75005 Paris, France
| | - Damien Laage
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ Paris 06, CNRS, PSL Research University , 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR , 75005 Paris, France
| | - James T Hynes
- PASTEUR, Département de Chimie, École Normale Supérieure, UPMC Univ Paris 06, CNRS, PSL Research University , 75005 Paris, France
- Department of Chemistry and Biochemistry, University of Colorado , Boulder, Colorado 80309-0215, United States
| |
Collapse
|
36
|
Affiliation(s)
- Emilie B. Guidez
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| | - Mark S. Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, United States
| |
Collapse
|
37
|
Cisneros G, Wikfeldt KT, Ojamäe L, Lu J, Xu Y, Torabifard H, Bartók AP, Csányi G, Molinero V, Paesani F. Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions. Chem Rev 2016; 116:7501-28. [PMID: 27186804 PMCID: PMC5450669 DOI: 10.1021/acs.chemrev.5b00644] [Citation(s) in RCA: 278] [Impact Index Per Article: 34.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2015] [Indexed: 12/17/2022]
Abstract
Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.
Collapse
Affiliation(s)
| | - Kjartan Thor Wikfeldt
- Science
Institute, University of Iceland, VR-III, 107, Reykjavik, Iceland
- Department
of Physics, Albanova, Stockholm University, S-106 91 Stockholm, Sweden
| | - Lars Ojamäe
- Department
of Chemistry, Linköping University, SE-581 83 Linköping, Sweden
| | - Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yao Xu
- Lehrstuhl
Physikalische Chemie II, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Hedieh Torabifard
- Department
of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Albert P. Bartók
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
38
|
Abstract
Whereas the interactions between water molecules are dominated by strongly directional hydrogen bonds (HBs), it was recently proposed that relatively weak, isotropic van der Waals (vdW) forces are essential for understanding the properties of liquid water and ice. This insight was derived from ab initio computer simulations, which provide an unbiased description of water at the atomic level and yield information on the underlying molecular forces. However, the high computational cost of such simulations prevents the systematic investigation of the influence of vdW forces on the thermodynamic anomalies of water. Here, we develop efficient ab initio-quality neural network potentials and use them to demonstrate that vdW interactions are crucial for the formation of water's density maximum and its negative volume of melting. Both phenomena can be explained by the flexibility of the HB network, which is the result of a delicate balance of weak vdW forces, causing, e.g., a pronounced expansion of the second solvation shell upon cooling that induces the density maximum.
Collapse
|
39
|
Agmon N, Bakker HJ, Campen RK, Henchman RH, Pohl P, Roke S, Thämer M, Hassanali A. Protons and Hydroxide Ions in Aqueous Systems. Chem Rev 2016; 116:7642-72. [PMID: 27314430 DOI: 10.1021/acs.chemrev.5b00736] [Citation(s) in RCA: 287] [Impact Index Per Article: 35.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.
Collapse
Affiliation(s)
- Noam Agmon
- The Fritz Haber Research Center, Institute of Chemistry, The Hebrew University of Jerusalem , Jerusalem 91904, Israel
| | - Huib J Bakker
- FOM Institute AMOLF , Science Park 104, 1098 XG Amsterdam, The Netherlands
| | - R Kramer Campen
- Fritz Haber Institute of the Max Planck Society , Faradayweg 4-6, 14195 Berlin, Germany
| | - Richard H Henchman
- Manchester Institute of Biotechnology, School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, United Kingdom
| | - Peter Pohl
- Johannes Kepler University Linz , Institute of Biophysics, Gruberstrasse 40, 4020 Linz, Austria
| | - Sylvie Roke
- Laboratory for Fundamental BioPhotonics (LBP), Institute of Bioengineering (IBI), and Institute of Material Science (IMX), School of Engineering (STI), and Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL) , CH-1015, Lausanne, Switzerland
| | - Martin Thämer
- Fritz Haber Institute of the Max Planck Society , Faradayweg 4-6, 14195 Berlin, Germany.,Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago , Chicago, Illinois 60637, United States
| | - Ali Hassanali
- CMSP Section, The Abdus Salaam International Center for Theoretical Physics , I-34151 Trieste, Italy
| |
Collapse
|
40
|
Cheng B, Behler J, Ceriotti M. Nuclear Quantum Effects in Water at the Triple Point: Using Theory as a Link Between Experiments. J Phys Chem Lett 2016; 7:2210-2215. [PMID: 27203358 DOI: 10.1021/acs.jpclett.6b00729] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
One of the most prominent consequences of the quantum nature of light atomic nuclei is that their kinetic energy does not follow a Maxwell-Boltzmann distribution. Deep inelastic neutron scattering (DINS) experiments can measure this effect. Thus, the nuclear quantum kinetic energy can be probed directly in both ordered and disordered samples. However, the relation between the quantum kinetic energy and the atomic environment is a very indirect one, and cross-validation with theoretical modeling is therefore urgently needed. Here, we use state of the art path integral molecular dynamics techniques to compute the kinetic energy of hydrogen and oxygen nuclei in liquid, solid, and gas-phase water close to the triple point, comparing three different interatomic potentials and validating our results against equilibrium isotope fractionation measurements. We will then show how accurate simulations can draw a link between extremely precise fractionation experiments and DINS, therefore establishing a reliable benchmark for future measurements and providing key insights to increase further the accuracy of interatomic potentials for water.
Collapse
Affiliation(s)
- Bingqing Cheng
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne , 1015 Lausanne, Switzerland
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum , 44801 Bochum, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne , 1015 Lausanne, Switzerland
| |
Collapse
|
41
|
Giese TJ, Panteva MT, Chen H, York DM. Multipolar Ewald methods, 2: applications using a quantum mechanical force field. J Chem Theory Comput 2016; 11:451-61. [PMID: 25691830 PMCID: PMC4325604 DOI: 10.1021/ct500799g] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Indexed: 01/08/2023]
Abstract
A fully quantum mechanical force field (QMFF) based on a modified “divide-and-conquer” (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.
Collapse
Affiliation(s)
- Timothy J Giese
- Center for Integrative Proteomics Research, BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | | | | | | |
Collapse
|
42
|
Amann-Winkel K, Bellissent-Funel MC, Bove LE, Loerting T, Nilsson A, Paciaroni A, Schlesinger D, Skinner L. X-ray and Neutron Scattering of Water. Chem Rev 2016; 116:7570-89. [DOI: 10.1021/acs.chemrev.5b00663] [Citation(s) in RCA: 137] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Katrin Amann-Winkel
- Department
of Physics, AlbaNova University Center, Stockholm University, SE-106
91, Stockholm, Sweden
| | | | - Livia E. Bove
- IMPMC, CNRS-UMR 7590, Université P&M Curie, 75252 Paris, France
- Institute
of Condensed Matter Physics, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Thomas Loerting
- Institute
of Physical Chemistry, University of Innsbruck, 6020 Innsbruck, Austria
| | - Anders Nilsson
- Department
of Physics, AlbaNova University Center, Stockholm University, SE-106
91, Stockholm, Sweden
| | - Alessandro Paciaroni
- Dipartimento
di Fisica e Geologia, Università di Perugia, Via Alessandro
Pascoli, I-06123 Perugia, Italy
| | - Daniel Schlesinger
- Department
of Physics, AlbaNova University Center, Stockholm University, SE-106
91, Stockholm, Sweden
| | - Lawrie Skinner
- Mineral
Physics Institute, Stony Brook University, Stony Brook, New York 11794-2100, United States
| |
Collapse
|
43
|
Nagata Y, Ohto T, Backus EHG, Bonn M. Molecular Modeling of Water Interfaces: From Molecular Spectroscopy to Thermodynamics. J Phys Chem B 2016; 120:3785-96. [DOI: 10.1021/acs.jpcb.6b01012] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Tatsuhiko Ohto
- Graduate
School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Ellen H. G. Backus
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Mischa Bonn
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
44
|
Śmiechowski M. Molecular hydrogen solvated in water – A computational study. J Chem Phys 2015; 143:244505. [DOI: 10.1063/1.4938571] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Maciej Śmiechowski
- Department of Physical Chemistry, Chemical Faculty, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
| |
Collapse
|
45
|
Del Ben M, Hutter J, VandeVondele J. Probing the structural and dynamical properties of liquid water with models including non-local electron correlation. J Chem Phys 2015; 143:054506. [PMID: 26254660 DOI: 10.1063/1.4927325] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.
Collapse
Affiliation(s)
- Mauro Del Ben
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Joost VandeVondele
- Department of Materials, ETH Zurich, Wolfgang-Pauli-Strasse 27, CH-8093 Zurich, Switzerland
| |
Collapse
|
46
|
Brorsen KR, Willow SY, Xantheas SS, Gordon MS. The Melting Temperature of Liquid Water with the Effective Fragment Potential. J Phys Chem Lett 2015; 6:3555-3559. [PMID: 26722723 DOI: 10.1021/acs.jpclett.5b01702] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The direct simulation of the solid-liquid water interface with the effective fragment potential (EFP) via the constant enthalpy and pressure (NPH) ensemble was used to estimate the melting temperature (T(m)) of ice-I(h). Initial configurations and velocities, taken from equilibrated constant pressure and temperature (NPT) simulations at P = 1 atm and T = 305 K, 325 K and 399 K, respectively, yielded corresponding T(m) values of 378 ± 16 K, 382 ± 14 K and 384 ± 15 K. These estimates are consistently higher than experiment, albeit to the same degree as previously reported estimates using density functional theory (DFT)-based Born-Oppenheimer simulations with the Becke-Lee-Yang-Parr functional plus dispersion corrections (BLYP-D).
Collapse
Affiliation(s)
- Kurt R Brorsen
- Ames Laboratory, US-DOE and Department of Chemistry, Iowa State University , Ames, Iowa 50011, United States
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sotiris S Xantheas
- Physical Sciences Division, Pacific Northwest National Laboratory , 902 Battelle Boulevard, P.O. Box 999, MS K1-83, Richland, Washington 99352, United States
| | - Mark S Gordon
- Ames Laboratory, US-DOE and Department of Chemistry, Iowa State University , Ames, Iowa 50011, United States
| |
Collapse
|
47
|
Yin J, Fenley AT, Henriksen NM, Gilson MK. Toward Improved Force-Field Accuracy through Sensitivity Analysis of Host-Guest Binding Thermodynamics. J Phys Chem B 2015; 119:10145-55. [PMID: 26181208 DOI: 10.1021/acs.jpcb.5b04262] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Improving the capability of atomistic computer models to predict the thermodynamics of noncovalent binding is critical for successful structure-based drug design, and the accuracy of such calculations remains limited by nonoptimal force field parameters. Ideally, one would incorporate protein-ligand affinity data into force field parametrization, but this would be inefficient and costly. We now demonstrate that sensitivity analysis can be used to efficiently tune Lennard-Jones parameters of aqueous host-guest systems for increasingly accurate calculations of binding enthalpy. These results highlight the promise of a comprehensive use of calorimetric host-guest binding data, along with existing validation data sets, to improve force field parameters for the simulation of noncovalent binding, with the ultimate goal of making protein-ligand modeling more accurate and hence speeding drug discovery.
Collapse
Affiliation(s)
- Jian Yin
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California 92093-0736, United States
| | - Andrew T Fenley
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California 92093-0736, United States
| | - Niel M Henriksen
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California 92093-0736, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, California 92093-0736, United States
| |
Collapse
|
48
|
Bankura A, Santra B, DiStasio RA, Swartz CW, Klein ML, Wu X. A systematic study of chloride ion solvation in water using van der Waals inclusive hybrid density functional theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1059959] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Arindam Bankura
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Charles W. Swartz
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University , Philadelphia, PA, USA
| | - Xifan Wu
- Department of Physics, Temple University , Philadelphia, PA, USA
| |
Collapse
|
49
|
Affiliation(s)
- Biswajit Santra
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | | | - Fausto Martelli
- Department of Chemistry, Princeton University , Princeton, NJ, USA
| | - Roberto Car
- Department of Chemistry, Princeton University , Princeton, NJ, USA
- Department of Physics, Princeton University , Princeton, NJ, USA
- Princeton Institute for the Science and Technology of Materials, Princeton University , Princeton, NJ, USA
- Program in Applied and Computational Mathematics, Princeton University , Princeton, NJ, USA
| |
Collapse
|
50
|
DiStasio RA, Santra B, Li Z, Wu X, Car R. The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water. J Chem Phys 2015; 141:084502. [PMID: 25173016 DOI: 10.1063/1.4893377] [Citation(s) in RCA: 183] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.
Collapse
Affiliation(s)
- Robert A DiStasio
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Zhaofeng Li
- Department of Physics, Princeton University, Princeton, New Jersey 08544, USA
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|