1
|
Piskulich ZA, Laage D, Thompson WH. A structure-dynamics relationship enables prediction of the water hydrogen bond exchange activation energy from experimental data. Chem Sci 2024; 15:2197-2204. [PMID: 38332825 PMCID: PMC10848719 DOI: 10.1039/d3sc04495e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Accepted: 12/23/2023] [Indexed: 02/10/2024] Open
Abstract
It has long been understood that the structural features of water are determined by hydrogen bonding (H-bonding) and that the exchange of, or "jumps" between, H-bond partners underlies many of the dynamical processes in water. Despite the importance of H-bond exchanges there is, as yet, no direct method for experimentally measuring the timescale of the process or its associated activation energy. Here, we identify and exploit relationships between water's structural and dynamical properties that provide an indirect route for determining the H-bond exchange activation energy from experimental data. Specifically, we show that the enthalpy and entropy determining the radial distribution function in liquid water are linearly correlated with the activation energies for H-bond jumps, OH reorientation, and diffusion. Using temperature-dependent measurements of the radial distribution function from the literature, we demonstrate how these correlations allow us to infer the value of the jump activation energy, Ea,0, from experimental results. This analysis gives Ea,0 = 3.43 kcal mol-1, which is in good agreement with that predicted by the TIP4P/2005 water model. We also illustrate other approaches for estimating this activation energy consistent with these estimates.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas Lawrence KS 66045 USA
| | - Damien Laage
- PASTEUR, Départment de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS Paris 75005 France
| | - Ward H Thompson
- Department of Chemistry, University of Kansas Lawrence KS 66045 USA
| |
Collapse
|
2
|
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
|
3
|
Malik R, Chandra A, Das B, Chandra A. Temperature Dependence of Non-Condon Effects in Two-Dimensional Vibrational Spectroscopy of Water. J Phys Chem B 2023; 127:2488-2498. [PMID: 36893383 DOI: 10.1021/acs.jpcb.2c06794] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Abstract
Non-Condon effects in vibrational spectroscopy refers to the dependence of a molecule's vibrational transition dipole and polarizability on the coordinates of the surrounding environment. Earlier studies have shown that such effects can be pronounced for hydrogen-bonded systems like liquid water. Here, we present a theoretical study of two-dimensional vibrational spectroscopy under the non-Condon and Condon approximations at varying temperatures. We have performed calculations of both two-dimensional infrared and two-dimensional vibrational Raman spectra to gain insights into the temperature dependence of non-Condon effects in nonlinear vibrational spectroscopy. The two-dimensional spectra are calculated for the OH vibration of interest in the isotopic dilution limit where the coupling between the oscillators is ignored. Generally, both the infrared and Raman line shapes undergo red shifts with decrease in temperature due to strengthening of hydrogen bonds and decrease in the fraction of OH modes with weaker or no hydrogen bonds. The infrared line shape is further red-shifted under the non-Condon effects at a given temperature, while the Raman line shape does not show any such red shift due to non-Condon effects. The spectral dynamics becomes slower on decrease of temperature due to slower hydrogen bond relaxation and, for a given temperature, the spectral diffusion occurs at a faster rate upon inclusion of non-Condon effects. The time scales of spectral diffusion extracted from different metrics agree well with each other and also with experiments. The changes in the spectrum due to non-Condon effects are found to be more significant at lower temperatures.
Collapse
Affiliation(s)
- Ravi Malik
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Abhilash Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Banshi Das
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
4
|
Hu K, Shirakashi R. Dynamic Electric Field Alignment Determines the Water Rotational Motion around Protein. J Phys Chem B 2023; 127:1376-1384. [PMID: 36749793 DOI: 10.1021/acs.jpcb.2c07405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Water rotational dynamics in biomolecular solution is crucial to evaluating and controlling biomolecule stability. In this molecular dynamics simulation (MD) study on lysozyme solutions, we present how the exerted internal electric field determines water rotational dynamics. We find that the relaxation time of water rotation is equivalent to that of the reorientation of the exerted overall electric field for every single water molecule, regardless of its translation mode. Namely, water molecular rotation synchronizes with the exerted field reorientation. We also map the reorientation process of the electric field at fixed points relative to protein in the solution, which displays the local hydration dynamics commensurate with the reported time-dependent fluorescence Stokes shift (TDFSS) measurements. Comparing the spatial distribution of local field reorientation relaxation time with that of rotational relaxation time, we further suggest that water rotation dynamics are subject to the reorientation of the local overall field within the hydration layer. While outside the hydration layer, the relaxation time of the local electric field reorientation is short enough (subpicosecond) to assume the δ function, showing the electric force with randomly changing orientation is applied to each water molecule.
Collapse
Affiliation(s)
- Kang Hu
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro City, Tokyo 153-8505, Japan.,Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Bunkyo-ku, Tokyo 113-0033, Japan
| | - Ryo Shirakashi
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro City, Tokyo 153-8505, Japan
| |
Collapse
|
5
|
Borkowski AK, Campbell NI, Thompson WH. Direct calculation of the temperature dependence of 2D-IR spectra: Urea in water. J Chem Phys 2023; 158:064507. [PMID: 36792517 DOI: 10.1063/5.0135627] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
A method for directly calculating the temperature derivative of two-dimensional infrared (2D-IR) spectra from simulations at a single temperature is presented. The approach is demonstrated by application to the OD stretching spectrum of isotopically dilute aqueous (HOD in H2O) solutions of urea as a function of concentration. Urea is an important osmolyte because of its ability to denature proteins, which has motivated significant interest in its effect on the structure and dynamics of water. The present results show that the temperature dependence of both the linear IR and 2D-IR spectra, which report on the underlying energetic driving forces, is more sensitive to urea concentration than the spectra themselves. Additional physical insight is provided by calculation of the contributions to the temperature derivative from different interactions, e.g., water-water, water-urea, and urea-urea, present in the system. Finally, it is demonstrated how 2D-IR spectra at other temperatures can be obtained from only room temperature simulations.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| | - N Ian Campbell
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| |
Collapse
|
6
|
Borkowski AK, Thompson WH. Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions. J Phys Chem B 2022; 126:6700-6712. [PMID: 36004804 DOI: 10.1021/acs.jpcb.2c03957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Hofmeister series is frequently used to rank ions based on their behavior from chaotropes ("structure breakers"), which weaken the surrounding hydrogen-bond network, to kosmotropes ("structure makers"), which enhance it. Here, we use fluctuation theory to investigate the energetic and entropic driving forces underlying the Hofmeister series for aqueous alkali-halide solutions. Specifically, we exploit the OH stretch infrared (IR) spectrum in isotopically dilute HOD/D2O solutions as a probe of the effect of the salt on the water properties for different concentrations and choice of halide anion. Fluctuation theory is used to calculate the temperature derivative of these IR spectra, including decomposition of the derivative into different energetic contributions. These contributions are used to determine the thermodynamic driving forces in terms of effective internal energy and entropic contributions. This analysis implicates entropic contributions as the key factor in the Hofmeister series behavior of the OH stretch IR spectra, while the effective internal energy is nearly ion-independent.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
7
|
Gomez A, Piskulich ZA, Thompson WH, Laage D. Water Diffusion Proceeds via a Hydrogen-Bond Jump Exchange Mechanism. J Phys Chem Lett 2022; 13:4660-4666. [PMID: 35604934 DOI: 10.1021/acs.jpclett.2c00825] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The self-diffusion of water molecules plays a key part in a broad range of essential processes in biochemistry, medical imaging, material science, and engineering. However, its molecular mechanism and the role played by the water hydrogen-bond network rearrangements are not known. Here we combine molecular dynamics simulations and analytic modeling to determine the molecular mechanism of water diffusion. We establish a quantitative connection between the water diffusion coefficient and hydrogen-bond jump exchanges, and identify the features that determine the underlying energetic barrier. We thus provide a unified framework to understand the coupling between translational, rotational, and hydrogen-bond dynamics in liquid water. It explains why these different dynamics do not necessarily exhibit identical temperature dependences although they all result from the same hydrogen-bond exchange events. The consequences for the understanding of water diffusion in supercooled conditions and for water transport in complex aqueous systems, including ionic, biological, and confined solutions, are discussed.
Collapse
Affiliation(s)
- Axel Gomez
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Damien Laage
- PASTEUR, Department of Chemistry, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
8
|
Piskulich ZA, Laage D, Thompson WH. Using Activation Energies to Elucidate Mechanisms of Water Dynamics. J Phys Chem A 2021; 125:9941-9952. [PMID: 34748353 DOI: 10.1021/acs.jpca.1c08020] [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/29/2022]
Abstract
Recent advances in the calculation of activation energies are shedding new light on the dynamical time scales of liquid water. In this Perspective, we examine how activation energies elucidate the central, but not singular, role of the exchange of hydrogen-bond (H-bond) partners that rearrange the H-bond network of water. The contributions of other motions to dynamical time scales and their associated activation energies are discussed along with one case, vibrational spectral diffusion, where H-bond exchanges are not mechanistically significant. Nascent progress on outstanding challenges, including descriptions of non-Arrhenius effects and activation volumes, are detailed along with some directions for future investigations.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Damien Laage
- PASTEUR, Department de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, Paris 75005, France
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
9
|
Piskulich ZA, Laird BB. Molecular Simulations of Phase Equilibria and Transport Properties in a Model CO 2-Expanded Lithium Perchlorate Electrolyte. J Phys Chem B 2021; 125:9341-9349. [PMID: 34351157 DOI: 10.1021/acs.jpcb.1c05369] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Carbon-dioxide (CO2)-expanded liquids, in which a significant mole fraction of CO2 is dissolved into an organic solvent, have been of significant interest, especially as catalytic support media. Because the CO2 mole fraction and density can be controlled over a significant range by changing the CO2 partial pressure, the transport properties of these solvents are highly tunable. Recently, these liquids have garnered interest as potential electrolyte solutions for catalytic electrochemistry; however, little is currently known about the influence of the electrolyte on CO2 expansion. In the present work, we use molecular-dynamics simulations to study diffusion and viscosity in a model lithium perchlorate electrolyte in CO2-expanded acetonitrile and demonstrate that these properties are highly dependent on the concentration of the electrolyte. Our present results indicate that the electrolyte slows down diffusion of both CO2 and MeCN, and that the slowed diffusion in the former is driven by changes in the activation entropy, whereas slowed diffusion in the latter is driven by changes in the activation energy.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Brian B Laird
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
10
|
Piskulich ZA, Thompson WH. Examining the Role of Different Molecular Interactions on Activation Energies and Activation Volumes in Liquid Water. J Chem Theory Comput 2021; 17:2659-2671. [PMID: 33819026 DOI: 10.1021/acs.jctc.0c01217] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
There are a large number of force fields available to model water in molecular dynamics simulations, which each have their own strengths and weaknesses in describing the behavior of the liquid. One particular weakness in many of these models is their description of dynamics away from ambient conditions, where their ability to reproduce measurements is mixed. To investigate this issue, we use the recently developed fluctuation theory for dynamics to directly evaluate measures of the local temperature and pressure dependence: the activation energy and the activation volume. We examine these activation parameters for hydrogen-bond jump exchange times, OH reorientation times, and diffusion coefficients calculated from the SPC/E, SPC/Fw, TIP3P-PME, TIP3P-PME/Fw, OPC3, TIP4P/2005, TIP4P/Ew, E3B2, and E3B3 water models. Activation energy decompositions available through the fluctuation theory approach provide mechanistic insight into the origins of different temperature dependences between the various models, as well as the influence of three-body effects and flexibility.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
11
|
Piskulich ZA, Laage D, Thompson WH. On the role of hydrogen-bond exchanges in the spectral diffusion of water. J Chem Phys 2021; 154:064501. [PMID: 33588543 DOI: 10.1063/5.0041270] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dynamics of a vibrational frequency in a condensed phase environment, i.e., the spectral diffusion, has attracted considerable interest over the last two decades. A significant impetus has been the development of two-dimensional infrared (2D-IR) photon-echo spectroscopy that represents a direct experimental probe of spectral diffusion, as measured by the frequency-frequency time correlation function (FFCF). In isotopically dilute water, which is perhaps the most thoroughly studied system, the standard interpretation of the longest timescale observed in the FFCF is that it is associated with hydrogen-bond exchange dynamics. Here, we investigate this connection by detailed analysis of both the spectral diffusion timescales and their associated activation energies. The latter are obtained from the recently developed fluctuation theory for the dynamics approach. The results show that the longest timescale of spectral diffusion obtained by the typical analysis used cannot be directly associated with hydrogen-bond exchanges. The hydrogen-bond exchange time does appear in the decay of the water FFCF, but only as an additional, small-amplitude (<3%) timescale. The dominant contribution to the long-time spectral diffusion dynamics is considerably shorter than the hydrogen-bond exchange time and exhibits a significantly smaller activation energy. It thus arises from hydrogen-bond rearrangements, which occur in between successful hydrogen-bond partner exchanges, and particularly from hydrogen bonds that transiently break before returning to the same acceptor.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| | - Damien Laage
- PASTEUR, Department of Chemistry, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| |
Collapse
|
12
|
Piskulich ZA, Thompson WH. Temperature Dependence of the Water Infrared Spectrum: Driving Forces, Isosbestic Points, and Predictions. J Phys Chem Lett 2020; 11:7762-7768. [PMID: 32852956 DOI: 10.1021/acs.jpclett.0c02301] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The temperature derivative of the infrared (IR) spectrum of HOD/D2O is directly calculated from simulations at a single temperature using a fluctuation theory approach. It is demonstrated, on the basis of an energetic decomposition of the derivative, that the blue shift with increasing temperature is associated with the competition between electrostatic and Lennard-Jones interactions. The same competition gives rise, where their contributions cancel, to a near isosbestic point. The derivative is further used to define an effective internal energy (and entropy) associated with the IR spectrum, and it is shown how a van't Hoff relation can be used to accurately predict the spectrum over a wide range of temperatures. These predictions also explain why a precise isosbestic point is not observed.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|