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
|
Piskulich ZA, Borkowski AK, Thompson WH. A Maxwell relation for dynamical timescales with application to the pressure and temperature dependence of water self-diffusion and shear viscosity. Phys Chem Chem Phys 2023; 25:12820-12832. [PMID: 37129891 DOI: 10.1039/d3cp01386c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
A Maxwell relation for a reaction rate constant (or other dynamical timescale) obtained under constant pressure, p, and temperature, T, is introduced and discussed. Examination of this relationship in the context of fluctuation theory provides insight into the p and T dependence of the timescale and the underlying molecular origins. This Maxwell relation motivates a suggestion for the general form of the timescale as a function of pressure and temperature. This is illustrated by accurately fitting simulation results and existing experimental data on the self-diffusion coefficient and shear viscosity of liquid water. A key advantage of this approach is that each fitting parameter is physically meaningful.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, KS 66045, USA.
| | | | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, KS 66045, USA.
| |
Collapse
|
3
|
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
|
4
|
Weigl P, Weißheit S, Pabst F, Kolmar H, Thiele CM, Walther T, Blochowicz T. Triplet States Reveal Slow Local Dynamics in the Solvation Shell of Biomolecules. J Phys Chem B 2022; 126:6324-6330. [PMID: 35973008 DOI: 10.1021/acs.jpcb.2c03784] [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
Protein hydration shell dynamics plays a pivotal role in biochemical processes such as protein folding, enzyme function, molecular recognition and interaction with biological membranes. Thus, it is crucial to understand the mobility of the solvation shell at the surface of biomolecules. Triplet state solvation dynamics can reveal the slow dynamics of the solvation shell. This is done in the present work without adding separate dye molecules but instead by using a phosphorescent subgroup of the biomolecule itself. In particular, we study a small heptapeptide in a glycerol-water mixture under cryoconservation conditions so that the system can be supercooled without crystallization. We find a significant slowing of molecules in the solvation shell in the millisecond range compared to the bulk. This opens up the possibility to unravel the nature of relaxation processes in the solvation shell usually overlapping at room temperature.
Collapse
Affiliation(s)
- Peter Weigl
- Institute for Condensed Matter Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany.,Institute for Applied Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany
| | - Susann Weißheit
- Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Florian Pabst
- Institute for Condensed Matter Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany
| | - Harald Kolmar
- Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Christina Marie Thiele
- Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, 64287 Darmstadt, Germany
| | - Thomas Walther
- Institute for Applied Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany
| | - Thomas Blochowicz
- Institute for Condensed Matter Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
Corsaro C, Fazio E. From Critical Point to Critical Point: The Two-States Model Describes Liquid Water Self-Diffusion from 623 to 126 K. Molecules 2021; 26:molecules26195899. [PMID: 34641442 PMCID: PMC8512083 DOI: 10.3390/molecules26195899] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Revised: 09/23/2021] [Accepted: 09/26/2021] [Indexed: 12/02/2022] Open
Abstract
Liquid’s behaviour, when close to critical points, is of extreme importance both for fundamental research and industrial applications. A detailed knowledge of the structural–dynamical correlations in their proximity is still today a target to reach. Liquid water anomalies are ascribed to the presence of a second liquid–liquid critical point, which seems to be located in the very deep supercooled regime, even below 200 K and at pressure around 2 kbar. In this work, the thermal behaviour of the self-diffusion coefficient for liquid water is analyzed, in terms of a two-states model, for the first time in a very wide thermal region (126 K < T < 623 K), including those of the two critical points. Further, the corresponding configurational entropy and isobaric-specific heat have been evaluated within the same interval. The two liquid states correspond to high and low-density water local structures that play a primary role on water dynamical behavior over 500 K.
Collapse
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
Borkowski AK, Piskulich ZA, Thompson WH. Examining the Hofmeister Series through Activation Energies: Water Diffusion in Aqueous Alkali-Halide Solutions. J Phys Chem B 2021; 125:350-359. [PMID: 33382267 DOI: 10.1021/acs.jpcb.0c09965] [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/30/2022]
Abstract
The effect of ions on the properties of aqueous solutions is often categorized in terms of the Hofmeister series that ranks them from chaotropes ("structure-breakers"), which weaken the surrounding hydrogen-bond network to kosmotropes ("structure-makers"), which enhance it. Here, we investigate the Hofmeister series in ∼1 M sodium-halide solutions using molecular dynamics simulations to calculate the effect of the identity and proximity of the halide anion on both the water diffusion coefficient and its activation energy. A recently developed method for calculating the activation energy from a single-temperature simulation is used, which also permits a rigorous decomposition into contributions from different interactions and motions. The mechanisms of the salt effects on the water dynamics are explored by separately considering water molecules based on their location relative to the ions. The results show that water diffusion is accelerated moving down the halide group from F- to I-. The behavior of the diffusion activation energy, Ea, is more complex, indicating a significant role for entropic effects. However, water molecules in the first or second solvation shell of an ion exhibit a decrease in Ea moving down the halide series and Na+ exhibits a larger effect than any of the anions. The Ea for water molecules within the second solvation shell of an ion are modest, indicating a short-ranged nature of the ion influence.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - 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
|
12
|
Monroe JI, Hatch HW, Mahynski NA, Shell MS, Shen VK. Extrapolation and interpolation strategies for efficiently estimating structural observables as a function of temperature and density. J Chem Phys 2020; 153:144101. [PMID: 33086808 DOI: 10.1063/5.0014282] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Thermodynamic extrapolation has previously been used to predict arbitrary structural observables in molecular simulations at temperatures (or relative chemical potentials in open-system mixtures) different from those at which the simulation was performed. This greatly reduces the computational cost in mapping out phase and structural transitions. In this work, we explore the limitations and accuracy of thermodynamic extrapolation applied to water, where qualitative shifts from anomalous to simple-fluid-like behavior are manifested through shifts in the liquid structure that occur as a function of both temperature and density. We present formulas for extrapolating in volume for canonical ensembles and demonstrate that linear extrapolations of water's structural properties are only accurate over a limited density range. On the other hand, linear extrapolation in temperature can be accurate across the entire liquid state. We contrast these extrapolations with classical perturbation theory techniques, which are more conservative and slowly converging. Indeed, we show that such behavior is expected by demonstrating exact relationships between extrapolation of free energies and well-known techniques to predict free energy differences. An ideal gas in an external field is also studied to more clearly explain these results for a toy system with fully analytical solutions. We also present a recursive interpolation strategy for predicting arbitrary structural properties of molecular fluids over a predefined range of state conditions, demonstrating its success in mapping qualitative shifts in water structure with density.
Collapse
Affiliation(s)
- Jacob I Monroe
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Harold W Hatch
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Nathan A Mahynski
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - M Scott Shell
- University of California - Santa Barbara, Santa Barbara, California 93106, USA
| | - Vincent K Shen
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
13
|
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
|
14
|
Piskulich ZA, Laage D, Thompson WH. Activation energies and the extended jump model: How temperature affects reorientation and hydrogen-bond exchange dynamics in water. J Chem Phys 2020; 153:074110. [DOI: 10.1063/5.0020015] [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)
- Zeke A. Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 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, Kansas 66045, USA
| |
Collapse
|