1
|
Rick SW, Thompson WH. Effects of polarizability and charge transfer on water dynamics and the underlying activation energies. J Chem Phys 2023; 158:2890774. [PMID: 37191215 DOI: 10.1063/5.0151253] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 05/02/2023] [Indexed: 05/17/2023] Open
Abstract
A large number of force fields have been proposed for describing the behavior of liquid water within classical atomistic simulations, particularly molecular dynamics. In the past two decades, models that incorporate molecular polarizability and even charge transfer have become more prevalent, in attempts to develop more accurate descriptions. These are frequently parameterized to reproduce the measured thermodynamics, phase behavior, and structure of water. On the other hand, the dynamics of water is rarely considered in the construction of these models, despite its importance in their ultimate applications. In this paper, we explore the structure and dynamics of polarizable and charge-transfer water models, with a focus on timescales that directly or indirectly relate to hydrogen bond (H-bond) making and breaking. Moreover, we use the recently developed fluctuation theory for dynamics to determine the temperature dependence of these properties to shed light on the driving forces. This approach provides key insight into the timescale activation energies through a rigorous decomposition into contributions from the different interactions, including polarization and charge transfer. The results show that charge transfer effects have a negligible effect on the activation energies. Furthermore, the same tension between electrostatic and van der Waals interactions that is found in fixed-charge water models also governs the behavior of polarizable models. The models are found to involve significant energy-entropy compensation, pointing to the importance of developing water models that accurately describe the temperature dependence of water structure and dynamics.
Collapse
Affiliation(s)
- Steven W Rick
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| |
Collapse
|
2
|
Han B, Isborn CM, Shi L. Incorporating Polarization and Charge Transfer into a Point-Charge Model for Water Using Machine Learning. J Phys Chem Lett 2023; 14:3869-3877. [PMID: 37067482 DOI: 10.1021/acs.jpclett.3c00036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Rigid nonpolarizable water models with fixed point charges have been widely employed in molecular dynamics simulations due to their efficiency and reasonable accuracy for the potential energy surface. However, the dipole moment surface of water is not necessarily well-described by the same fixed charges, leading to failure in reproducing dipole-related properties. Here, we developed a machine-learning model trained against electronic structure data to assign point charges for water, and the resulting dipole moment surface significantly improved the predictions of the dielectric constant and the low-frequency IR spectrum of liquid water. Our analysis reveals that within our atom-centered point-charge description of the dipole moment surface, the intermolecular charge transfer is the major source of the peak intensity at 200 cm-1, whereas the intramolecular polarization controls the enhancement of the dielectric constant. The effects of exact Hartree-Fock exchange in the hybrid density functional on these properties are also discussed.
Collapse
Affiliation(s)
- Bowen Han
- Chemistry and Biochemistry, University of California, Merced, California 95343, United States
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Biochemistry, University of California, Merced, California 95343, United States
| |
Collapse
|
3
|
Begušić T, Blake GA. Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality. Nat Commun 2023; 14:1950. [PMID: 37029146 PMCID: PMC10082090 DOI: 10.1038/s41467-023-37667-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
4
|
Töpfer K, Käser S, Meuwly M. Double proton transfer in hydrated formic acid dimer: Interplay of spatial symmetry and solvent-generated force on reactivity. Phys Chem Chem Phys 2022; 24:13869-13882. [PMID: 35620978 PMCID: PMC9176184 DOI: 10.1039/d2cp01583h] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The double proton transfer (DPT) reaction in the hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed – as in a conventional empirical force field – and conformationally fluctuating – as available from the ML-based PES – charge models for FAD show a significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% (∼1 kcal mol−1) between 300 K and 600 K. The rate for DPT is largest, ∼1 ns−1, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O–H⋯O hydrogen bond whereas the second proton transfer is significantly controlled by the O–O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors. Simulation of double proton transfer in formic acid dimer by reactive ML potential in explicit molecular mechanics water solvent.![]()
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
5
|
Begušić T, Tao X, Blake GA, Miller TF. Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy. J Chem Phys 2022; 156:131102. [PMID: 35395895 DOI: 10.1063/5.0087156] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
6
|
Torii H. Singular value decomposition analysis of the electron density changes occurring upon electrostatic polarization of water. RSC Adv 2022; 12:2564-2573. [PMID: 35425301 PMCID: PMC8979083 DOI: 10.1039/d1ra06649h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2021] [Accepted: 01/12/2022] [Indexed: 12/31/2022] Open
Abstract
In-depth elucidation of how molecules are electrically polarized would be one key factor for understanding the properties of those molecules under various thermodynamic and/or spatial conditions. Here this problem is tackled for the case of hydrogen-bonded water by conducting singular value decomposition of the electron density changes that occur upon electrostatic polarization. It is shown that all those electron density changes are approximately described as linear combinations of ten orthonormal basis “vectors”. One main component is the interatomic charge transfer through each OH bond, while some others are characterized as the atomic dipolar polarizations, meaning that both of these components are important for the electrostatic polarization of water. The interaction parameters that reasonably well reproduce the induced dipole moments are derived, which indicate the extent of mixing of the two components in electrostatic polarization. The main features of the electron density changes that occur upon electrostatic polarization of water are elucidated by conducting singular value decomposition analysis of those changes.![]()
Collapse
Affiliation(s)
- Hajime Torii
- Department of Applied Chemistry and Biochemical Engineering, Faculty of Engineering, Shizuoka University 3-5-1 Johoku, Naka-ku Hamamatsu 432-8561 Japan +81-53-478-1624 +81-53-478-1624.,Department of Optoelectronics and Nanostructure Science, Graduate School of Science and Technology, Shizuoka University 3-5-1 Johoku, Naka-ku Hamamatsu 432-8561 Japan
| |
Collapse
|
7
|
Banks PA, Burgess L, Ruggiero MT. The necessity of periodic boundary conditions for the accurate calculation of crystalline terahertz spectra. Phys Chem Chem Phys 2021; 23:20038-20051. [PMID: 34518858 DOI: 10.1039/d1cp02496e] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Terahertz vibrational spectroscopy has emerged as a powerful spectroscopic technique, providing valuable information regarding long-range interactions - and associated collective dynamics - occurring in solids. However, the terahertz sciences are relatively nascent, and there have been significant advances over the last several decades that have profoundly influenced the interpretation and assignment of experimental terahertz spectra. Specifically, because there do not exist any functional group or material-specific terahertz transitions, it is not possible to interpret experimental spectra without additional analysis, specifically, computational simulations. Over the years simulations utilizing periodic boundary conditions have proven to be most successful for reproducing experimental terahertz dynamics, due to the ability of the calculations to accurately take long-range forces into account. On the other hand, there are numerous reports in the literature that utilize gas phase cluster geometries, to varying levels of apparent success. This perspective will provide a concise introduction into the terahertz sciences, specifically terahertz spectroscopy, followed by an evaluation of gas phase and periodic simulations for the assignment of crystalline terahertz spectra, highlighting potential pitfalls and good practice for future endeavors.
Collapse
Affiliation(s)
- Peter A Banks
- Department of Chemistry, University of Vermont, 82 University Place, Burlington, Vermont 05405, USA.
| | - Luke Burgess
- Department of Chemistry, University of Vermont, 82 University Place, Burlington, Vermont 05405, USA.
| | - Michael T Ruggiero
- Department of Chemistry, University of Vermont, 82 University Place, Burlington, Vermont 05405, USA.
| |
Collapse
|
8
|
Shiraga K, Fujii Y, Koreeda A, Tanaka K, Arikawa T, Ogawa Y. Dynamical Collectivity and Nuclear Quantum Effects on the Intermolecular Stretching Mode of Liquid Water. J Phys Chem B 2021; 125:1632-1639. [PMID: 33393775 DOI: 10.1021/acs.jpcb.0c10154] [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/28/2022]
Abstract
This study investigated the broadband terahertz and low-frequency Raman spectroscopy of liquid water (H2O, D2O, and H218O) over 2 decades of frequency to address long-standing challenges regarding the interpretation of the intermolecular stretching mode at around 5 THz. We experimentally demonstrated that the intermolecular stretching mode of liquid water obtained via terahertz spectroscopy is significantly redshifted and broadened compared with that via Raman. This result was rationalized by the enhanced dynamical collectivity probed by terahertz spectroscopy, although both have a common origin in the kinetic motion. Their temperature and isotope dependences emphasize the significance of oscillation mass in determining the intermolecular stretching lineshape, while quantum effects cannot be overlooked in both terahertz and low-frequency Raman spectra.
Collapse
Affiliation(s)
- Keiichiro Shiraga
- Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan.,RIKEN Center for Integrative Medical Sciences, Yokohama 230-0045, Japan
| | - Yasuhiro Fujii
- Department of Physical Sciences, Ritsumeikan University, Kusatsu 525-8577, Japan
| | - Akitoshi Koreeda
- Department of Physical Sciences, Ritsumeikan University, Kusatsu 525-8577, Japan
| | - Koichiro Tanaka
- Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan.,Institute for Integrated Cell-Material Sciences (iCeMS), Kyoto University, Kyoto 606-8501, Japan
| | - Takashi Arikawa
- Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Yuichi Ogawa
- Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
9
|
Han B, Isborn CM, Shi L. Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models. J Chem Theory Comput 2021; 17:889-901. [PMID: 33405925 DOI: 10.1021/acs.jctc.0c01102] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.
Collapse
Affiliation(s)
- Bowen Han
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
10
|
Sidler D, Hamm P. A Feynman diagram description of the 2D-Raman-THz response of amorphous ice. J Chem Phys 2020; 153:044502. [PMID: 32752676 DOI: 10.1063/5.0018485] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The 2D-Raman-THz response in all possible time-orderings (Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman) of amorphous water ice is calculated in two ways: from atomistic molecular dynamics simulations and with the help of a Feynman diagram model, the latter of which power-expands the potential energy surface and the dipole and polarizability surfaces up to leading order. Comparing both results allows one to dissect the 2D-Raman-THz response into contributions from mechanical anharmonicity, as well as electrical dipole and polarizability anharmonicities. Mechanical anharmonicity dominates the 2D-Raman-THz response of the hydrogen-bond stretching and hydrogen-bond bending bands of water, and dipole anharmonicity dominates that of the librational band, while the contribution of polarizability anharmonicity is comparably weak. A distinct echo of the hydrogen-bond stretching band is observed for the THz-Raman-THz pulse sequence, again dominated by mechanical anharmonicity. A peculiar mechanism is discussed, which is based on the coupling between the many normal modes within the hydrogen-bond stretching band and which will inevitably generate such an echo for an amorphous structure.
Collapse
Affiliation(s)
- David Sidler
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
11
|
Lewis NHC, Iscen A, Felts A, Dereka B, Schatz GC, Tokmakoff A. Vibrational Probe of Aqueous Electrolytes: The Field Is Not Enough. J Phys Chem B 2020; 124:7013-7026. [DOI: 10.1021/acs.jpcb.0c05510] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Nicholas H. C. Lewis
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Aysenur Iscen
- Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Alanna Felts
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Bogdan Dereka
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - George C. Schatz
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
12
|
Koner D, Salehi SM, Mondal P, Meuwly M. Non-conventional force fields for applications in spectroscopy and chemical
reaction dynamics. J Chem Phys 2020; 153:010901. [DOI: 10.1063/5.0009628] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Seyedeh Maryam Salehi
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Padmabati Mondal
- Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra
Pradesh, India
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel,
Switzerland and Department of Chemistry, Brown University, Providence, Rhode Island, USA
| |
Collapse
|
13
|
Carlson S, Brünig FN, Loche P, Bonthuis DJ, Netz RR. Exploring the Absorption Spectrum of Simulated Water from MHz to Infrared. J Phys Chem A 2020; 124:5599-5605. [DOI: 10.1021/acs.jpca.0c04063] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Shane Carlson
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Florian N. Brünig
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Philip Loche
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Douwe Jan Bonthuis
- Institute of Theoretical and Computational Physics, Graz University of Technology, 8010 Graz, Austria
| | - Roland R. Netz
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| |
Collapse
|
14
|
Schröder C, Lyons A, Rick SW. Polarizable MD simulations of ionic liquids: How does additional charge transfer change the dynamics? Phys Chem Chem Phys 2020; 22:467-477. [DOI: 10.1039/c9cp05478b] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A new model for treating charge transfer in ionic liquids is developed and applied to 1-ethyl-3-methylimidazolium tetrafluoroborate. The model allows for us to examine the roles of charge transfer, polarizability, and charge scaling effects on the dynamics of ionic liquids.
Collapse
Affiliation(s)
- Christian Schröder
- University of Vienna
- Faculty of Chemistry
- Department of Computational Biological Chemistry
- A-1090 Vienna
- Austria
| | - Alex Lyons
- University of New Orleans
- Department of Chemistry
- New Orleans
- USA
| | - Steven W. Rick
- University of New Orleans
- Department of Chemistry
- New Orleans
- USA
| |
Collapse
|
15
|
Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk. J Phys Chem B 2019; 123:9846-9861. [PMID: 31647873 DOI: 10.1021/acs.jpcb.9b03258] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.
Collapse
|
16
|
Affiliation(s)
- Peter Hamm
- Department of Chemistry, University of Zurich, Winterthurerstr. 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
17
|
Sidler D, Hamm P. Feynman diagram description of 2D-Raman-THz spectroscopy applied to water. J Chem Phys 2019; 150:044202. [DOI: 10.1063/1.5079497] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- David Sidler
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | - Peter Hamm
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|