1
|
The collective burst mechanism of angular jumps in liquid water. Nat Commun 2023; 14:1345. [PMID: 36906703 PMCID: PMC10008639 DOI: 10.1038/s41467-023-37069-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Accepted: 02/24/2023] [Indexed: 03/13/2023] Open
Abstract
Understanding the microscopic origins of collective reorientational motions in aqueous systems requires techniques that allow us to reach beyond our chemical imagination. Herein, we elucidate a mechanism using a protocol that automatically detects abrupt motions in reorientational dynamics, showing that large angular jumps in liquid water involve highly cooperative orchestrated motions. Our automatized detection of angular fluctuations, unravels a heterogeneity in the type of angular jumps occurring concertedly in the system. We show that large orientational motions require a highly collective dynamical process involving correlated motion of many water molecules in the hydrogen-bond network that form spatially connected clusters going beyond the local angular jump mechanism. This phenomenon is rooted in the collective fluctuations of the network topology which results in the creation of defects in waves on the THz timescale. The mechanism we propose involves a cascade of hydrogen-bond fluctuations underlying angular jumps and provides new insights into the current localized picture of angular jumps, and its wide use in the interpretations of numerous spectroscopies as well in reorientational dynamics of water near biological and inorganic systems. The role of finite size effects, as well as of the chosen water model, on the collective reorientation is also elucidated.
Collapse
|
2
|
Yu H, Zhang Q, Zhuang W. Comparative analysis of hydration layer reorientation dynamics of antifreeze protein and protein cytochrome P450. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2203038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Antifreeze proteins (AFPs) inhibit ice re-crystallization by a mechanism remaining largely elusive. Dynamics of AFPs’ hydration water and its involvement in the antifreeze activity have not been identified conclusively. We herein, by simulation and theory, examined the water reorientation dynamics in the first hydration layer of an AFP from the spruce budworm, Choristoneura fumiferana, compared with a protein cytochrome P450 (CYP). The increase of potential acceptor water molecules around donor water molecules leads to the acceleration of hydrogen bond exchange between water molecules. Therefore, the jump reorientation of water molecules around the AFP active region is accelerated. Due to the mutual coupling and excitation of hydrogen bond exchange, with the acceleration of hydrogen bond exchange, the rearrangement of the hydrogen bond network and the frame reorientation of water are accelerated. Therefore, the water reorientation dynamics of AFP is faster than that of CYP. The results of this study provide a new physical image of antifreeze protein and a new understanding of the antifreeze mechanism of antifreeze proteins.
Collapse
Affiliation(s)
- Hongfeng Yu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Qiang Zhang
- College of Chemistry and Material Sciences, Inner Mongolia Minzu University, Tongliao 028043, China
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, China
- Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China
| |
Collapse
|
3
|
Liu C, Zhang Y, Zhang J, Wang J, Li W, Wang W. Interplay between translational diffusion and large-amplitude angular jumps of water molecules. J Chem Phys 2018; 148:184502. [DOI: 10.1063/1.5017935] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chao Liu
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Yangyang Zhang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Jian Zhang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Jun Wang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Wenfei Li
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Wei Wang
- National Laboratory of Solid State Microstructure, Department of Physics, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| |
Collapse
|
4
|
Zhang Q, Chen H, Wu T, Jin T, Pan Z, Zheng J, Gao Y, Zhuang W. The opposite effects of sodium and potassium cations on water dynamics. Chem Sci 2017; 8:1429-1435. [PMID: 28451283 PMCID: PMC5390786 DOI: 10.1039/c6sc03320b] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2016] [Accepted: 10/13/2016] [Indexed: 01/05/2023] Open
Abstract
Water rotational dynamics in NaSCN and KSCN solutions at a series of concentrations are investigated using femtosecond infrared spectroscopy and theory. Femtosecond infrared measurements, consistent with previous NMR observations, detect that sodium slows down while potassium accelerates the water O-H bond rotation. Results of reported neutron scattering measurements, on the other hand, suggested that these two cations have similar structure-breaking effects on water, and therefore should both accelerate water rotation through the presumably dominating large-amplitude angular jump component. To explain this discrepancy, theoretical studies with both classical and ab initio models were carried out, which indicate that both ions indeed accelerate the large-amplitude angular jump rotation of the water molecules, while the observed cation specific effect originates from the non-negligible opposite impact of the sodium and potassium cations on the diffusive rotation of water molecules.
Collapse
Affiliation(s)
- Qiang Zhang
- State Key Laboratory of Structural Chemistry , Fujian Institute of Research on the Structure of Matter , Chinese Academy of Sciences , Fuzhou , Fujian 350002 , China .
- Department of Chemistry , Bohai University , Jinzhou 121013 , China
| | - Hailong Chen
- Department of Chemistry , Rice University , Houston , TX 77005 , USA .
| | - Tianmin Wu
- Department of Chemical Physics , University of Science and Technology of China , Hefei , Anhui 230026 , P. R. China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry , Xiamen University , Xiamen , Fujian 361005 , China
| | - Tan Jin
- State Key Laboratory of Structural Chemistry , Fujian Institute of Research on the Structure of Matter , Chinese Academy of Sciences , Fuzhou , Fujian 350002 , China .
| | - Zhijun Pan
- Department of Chemical Physics , University of Science and Technology of China , Hefei , Anhui 230026 , P. R. China
| | - Junrong Zheng
- Department of Chemistry , Rice University , Houston , TX 77005 , USA .
- College of Chemistry and Molecular Engineering , Beijing National Laboratory for Molecular Sciences , Peking University , Beijing 100871 , China .
| | - Yiqin Gao
- College of Chemistry and Molecular Engineering , Beijing National Laboratory for Molecular Sciences , Peking University , Beijing 100871 , China .
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry , Fujian Institute of Research on the Structure of Matter , Chinese Academy of Sciences , Fuzhou , Fujian 350002 , China .
| |
Collapse
|
5
|
Henchman RH. Water's dual nature and its continuously changing hydrogen bonds. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:384001. [PMID: 27447299 DOI: 10.1088/0953-8984/28/38/384001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
A model is proposed for liquid water that is a continuum between the ordered state with predominantly tetrahedral coordination, linear hydrogen bonds and activated dynamics and a disordered state with a continuous distribution of multiple coordinations, multiple types of hydrogen bond, and diffusive dynamics, similar to that of normal liquids. Central to water's heterogeneous structure is the ability of hydrogen to donate to either one acceptor in a conventional linear hydrogen bond or to multiple acceptors as a furcated hydrogen. Linear hydrogen bonds are marked by slow, activated kinetics for hydrogen-bond switching to more crowded acceptors and sharp first peaks in the hydrogen-oxygen radial distribution function. Furcated hydrogens, equivalent to free, broken, dangling or distorted hydrogens, have barrierless, rapid kinetics and poorly defined first peaks in their hydrogen-oxygen radial distribution function. They involve the weakest donor in a local excess of donors, such that barrierless whole-molecule vibration rapidly swaps them between the linear and furcated forms. Despite the low number of furcated hydrogens and their transient existence, they are readily created in a single hydrogen-bond switch and free up the dynamics of numerous surrounding molecules, bringing about the disordered state. Hydrogens in the ordered state switch with activated dynamics to make the non-tetrahedral coordinations of the disordered state, which can also combine to make the ordered state. Consequently, the ordered and disordered states are both connected by diffusive dynamics and differentiated by activated dynamics, bringing about water's continuous heterogeneity.
Collapse
Affiliation(s)
- Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, UK. School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| |
Collapse
|
6
|
Demontis P, Gulín-González J, Masia M, Sant M, Suffritti GB. The interplay between dynamic heterogeneities and structure of bulk liquid water: A molecular dynamics simulation study. J Chem Phys 2016; 142:244507. [PMID: 26133441 DOI: 10.1063/1.4922930] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In order to study the interplay between dynamical heterogeneities and structural properties of bulk liquid water in the temperature range 130-350 K, thus including the supercooled regime, we use the explicit trend of the distribution functions of some molecular properties, namely, the rotational relaxation constants, the atomic mean-square displacements, the relaxation of the cross correlation functions between the linear and squared displacements of H and O atoms of each molecule, the tetrahedral order parameter q and, finally, the number of nearest neighbors (NNs) and of hydrogen bonds (HBs) per molecule. Two different potentials are considered: TIP4P-Ew and a model developed in this laboratory for the study of nanoconfined water. The results are similar for the dynamical properties, but are markedly different for the structural characteristics. In particular, for temperatures higher than that of the dynamic crossover between "fragile" (at higher temperatures) and "strong" (at lower temperatures) liquid behaviors detected around 207 K, the rotational relaxation of supercooled water appears to be remarkably homogeneous. However, the structural parameters (number of NNs and of HBs, as well as q) do not show homogeneous distributions, and these distributions are different for the two water models. Another dynamic crossover between "fragile" (at lower temperatures) and "strong" (at higher temperatures) liquid behaviors, corresponding to the one found experimentally at T(∗) ∼ 315 ± 5 K, was spotted at T(∗) ∼ 283 K and T(∗) ∼ 276 K for the TIP4P-Ew and the model developed in this laboratory, respectively. It was detected from the trend of Arrhenius plots of dynamic quantities and from the onset of a further heterogeneity in the rotational relaxation. To our best knowledge, it is the first time that this dynamical crossover is detected in computer simulations of bulk water. On the basis of the simulation results, the possible mechanisms of the two crossovers at molecular level are discussed.
Collapse
Affiliation(s)
- Pierfranco Demontis
- Dipartimento di Chimica e Farmacia, Università degli studi di Sassari, Sassari, Italy
| | - Jorge Gulín-González
- Grupo de Matemática y Física Computacionales, Universidad de las Ciencias Informáticas (UCI), Carretera a San Antonio de los Baños, Km 21/2, La Lisa, La Habana, Cuba
| | - Marco Masia
- Dipartimento di Chimica e Farmacia, Università degli studi di Sassari, Sassari, Italy
| | - Marco Sant
- Dipartimento di Chimica e Farmacia, Università degli studi di Sassari, Sassari, Italy
| | - Giuseppe B Suffritti
- Dipartimento di Chimica e Farmacia, Università degli studi di Sassari, Sassari, Italy
| |
Collapse
|
7
|
Mukherjee B. Microscopic origin of temporal heterogeneities in translational dynamics of liquid water. J Chem Phys 2015; 143:054503. [PMID: 26254657 DOI: 10.1063/1.4927709] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Liquid water is known to reorient via a combination of large angular jumps (due to exchange of hydrogen bonding (H-bond) partners) and diffusive orientations. Translation of the molecule undergoing the orientational jump and its initial and final H-bond acceptors plays a key role in the microscopic reorientation process. Here, we partition the translational dynamics into those occurring during intervals when rotating water molecules (and their initial and final H-bonding partners) undergo orientational jump and those arising when molecules wait between consecutive orientational jumps. These intervals are chosen in such a way that none of the four possible H-bonds involving the chosen water molecule undergo an exchange process within its duration. Translational dynamics is analysed in terms of the distribution of particle displacements, van Hove functions, and its moments. We observe that the translational dynamics, calculated from molecular dynamics simulations of liquid water, is fastest during the orientational jumps and slowest during periods of waiting. The translational dynamics during all temporal intervals shows an intermediate behaviour. This is the microscopic origin of temporal dynamic heterogeneity in liquid water, which is mild at 300 K and systematically increases with supercooling. Study of such partitioned dynamics in supercooled water shows increased disparity in dynamics occurring in the two different types of intervals. Nature of the distribution of particle displacements in supercooled water is investigated and it reveals signatures non-Gaussian behaviour.
Collapse
Affiliation(s)
- Biswaroop Mukherjee
- Thematic Unit of Excellence-Computational Materials Science, S. N. Bose National Centre for Basic Sciences, Block JD, Sector-III, Salt Lake, Kolkata 700098, India
| |
Collapse
|
8
|
Hartkamp R, Coasne B. Structure and transport of aqueous electrolytes: From simple halides to radionuclide ions. J Chem Phys 2014; 141:124508. [DOI: 10.1063/1.4896380] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Suffritti GB, Demontis P, Gulín-González J, Masia M. Distributions of single-molecule properties as tools for the study of dynamical heterogeneities in nanoconfined water. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2014; 26:155103. [PMID: 24675399 DOI: 10.1088/0953-8984/26/15/155103] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
The explicit trend of the distribution functions of single-molecule rotational relaxation constants and atomic mean-square displacement are used to study the dynamical heterogeneities in nanoconfined water. The trend of the single-molecule properties distributions is related to the dynamic heterogeneities, and to the dynamic crossovers found in water clusters of different shapes and sizes and confined in a variety of zeolites. This was true in all the cases that were considered, in spite of the various shapes and sizes of the clusters. It is confirmed that the high temperature dynamical crossover occurring in the temperature range 200-230 K can be interpreted at a molecular level as the formation of almost translationally rigid clusters, characterized by some rotational freedom, hydrogen bond exchange and translational jumps as cage-to-cage processes. We also suggest a mechanism for the low temperature dynamical crossover (LTDC), falling in the temperature range 150-185 K, through which the adsorbed water clusters are made of nearly rigid sub-clusters, slightly mismatched, and thus permitting a relatively free librational motion at their borders. It appears that the condition required for LTDC to occur is the presence of highly heterogeneous environments for the adsorbed molecules, with some dangling hydrogen bonds or weaker than water-water hydrogen bonds. Under these conditions some dynamics are permitted at very low temperature, although most rotational motion is frozen. Therefore, it is unlikely, though not entirely excluded, that LTDC will be found in supercooled bulk water where no heterogeneous interface is present.
Collapse
Affiliation(s)
- G B Suffritti
- Dipartimento di Chimica e Farmacia, Università di Sassari and INSTM, Unità di ricerca di Sassari, Via Vienna 2, I-07100 Sassari, Italy
| | | | | | | |
Collapse
|
10
|
Lepodise LM, Horvat J, Lewis RA. Collective librations of water molecules in the crystal lattice of rubidium bromide: experiment and simulation. Phys Chem Chem Phys 2013; 15:20252-61. [PMID: 24165585 DOI: 10.1039/c3cp53667j] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Terahertz spectroscopy of RbBr reveals four prominent absorption lines at room temperature and a further 15 lines at 10 K. Via density-functional-theory (DFT) numerical modelling using the PBE0 hybrid GGA functional, all the absorptions are identified as correlated librations of water molecules in the RbBr lattice. Each libration mode is a combination of rocking, wagging and twisting motions of the water molecules. The number of libration lines and numerical modelling show that the C2v symmetry of water in RbBr is broken. Our modelling shows that the distribution of libration amplitudes and phases for different water molecules in the RbBr unit cell varies greatly between the different modes. All librational lines red-shift with increasing temperature. The rate of change for most lines is in the range 60-90 MHz K(-1) (or (2-3) × 10(-3) cm(-1) K(-1)). Two lines shift more rapidly with temperature, at rates of 240 and 300 MHz K(-1) (or (8 and 10) × 10(-3) cm(-1) K(-1)), respectively. Furthermore, the temperature dependence of the linewidth distinguishes two groups of lines. For one group, with weak linear temperature dependence of linewidth, cubic anharmonic terms in the RbBr crystal field are significant. This group is mainly associated with fully symmetric correlated librations. For the second group, with strong non-linear temperature dependence of the linewidth, quartic anharmonic terms in the RbBr crystal field are significant. However, the distribution of libration amplitudes, as well as the type of libration modes, influence the temperature dependence of the red shift, the linewidth, and the intensity, as well. Our combined experimental and theoretical investigation confirms the necessity of obtaining low-temperature data to observe all the calculated modes; moreover, the richness of detail in the temperature dependence of the data invites further modelling spanning a range of temperatures.
Collapse
Affiliation(s)
- L M Lepodise
- Institute for Superconducting and Electronic Materials and School of Physics, University of Wollongong, Wollongong, New South Wales 2522, Australia.
| | | | | |
Collapse
|