Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Structural dynamics of supercooled water from quasielastic neutron scattering and molecular simulations.

Qvist, Johan LU ; Schober, Helmut and Halle, Bertil LU (2011) In Journal of Chemical Physics 134(14).
Abstract
One of the outstanding challenges presented by liquid water is to understand how molecules can move on a picosecond time scale despite being incorporated in a three-dimensional network of relatively strong H-bonds. This challenge is exacerbated in the supercooled state, where the dramatic slowing down of structural dynamics is reminiscent of the, equally poorly understood, generic behavior of liquids near the glass transition temperature. By probing single-molecule dynamics on a wide range of time and length scales, quasielastic neutron scattering (QENS) can potentially reveal the mechanistic details of water's structural dynamics, but because of interpretational ambiguities this potential has not been fully realized. To resolve these... (More)
One of the outstanding challenges presented by liquid water is to understand how molecules can move on a picosecond time scale despite being incorporated in a three-dimensional network of relatively strong H-bonds. This challenge is exacerbated in the supercooled state, where the dramatic slowing down of structural dynamics is reminiscent of the, equally poorly understood, generic behavior of liquids near the glass transition temperature. By probing single-molecule dynamics on a wide range of time and length scales, quasielastic neutron scattering (QENS) can potentially reveal the mechanistic details of water's structural dynamics, but because of interpretational ambiguities this potential has not been fully realized. To resolve these issues, we present here an extensive set of high-quality QENS data from water in the range 253-293 K and a corresponding set of molecular dynamics (MD) simulations to facilitate and validate the interpretation. Using a model-free approach, we analyze the QENS data in terms of two motional components. Based on the dynamical clustering observed in MD trajectories, we identify these components with two distinct types of structural dynamics: picosecond local (L) structural fluctuations within dynamical basins and slower interbasin jumps (J). The Q-dependence of the dominant QENS component, associated with J dynamics, can be quantitatively rationalized with a continuous-time random walk (CTRW) model with an apparent jump length that depends on low-order moments of the jump length and waiting time distributions. Using a simple coarse-graining algorithm to quantitatively identify dynamical basins, we map the Newtonian MD trajectory on a CTRW trajectory, from which the jump length and waiting time distributions are computed. The jump length distribution is Gaussian and the rms jump length increases from 1.5 to 1.9 Å as the temperature increases from 253 to 293 K. The rms basin radius increases from 0.71 to 0.75 Å over the same range. The waiting time distribution is exponential at all investigated temperatures, ruling out significant dynamical heterogeneity. However, a simulation at 238 K reveals a small but significant dynamical heterogeneity. The macroscopic diffusion coefficient deduced from the QENS data agrees quantitatively with NMR and tracer results. We compare our QENS analysis with existing approaches, arguing that the apparent dynamical heterogeneity implied by stretched exponential fitting functions results from the failure to distinguish intrabasin (L) from interbasin (J) structural dynamics. We propose that the apparent dynamical singularity at ∼220 K corresponds to freezing out of J dynamics, while the calorimetric glass transition corresponds to freezing out of L dynamics. (Less)
Please use this url to cite or link to this publication:
author
; and
organization
publishing date
type
Contribution to journal
publication status
published
subject
in
Journal of Chemical Physics
volume
134
issue
14
article number
144508
publisher
American Institute of Physics (AIP)
external identifiers
  • wos:000289582300031
  • pmid:21495765
  • scopus:79954617035
  • pmid:21495765
ISSN
0021-9606
DOI
10.1063/1.3578472
language
English
LU publication?
yes
id
48bdf626-b623-48a8-b099-f3aac0f4cebf (old id 1937073)
date added to LUP
2016-04-01 11:16:50
date last changed
2022-01-26 06:52:21
@article{48bdf626-b623-48a8-b099-f3aac0f4cebf,
  abstract     = {{One of the outstanding challenges presented by liquid water is to understand how molecules can move on a picosecond time scale despite being incorporated in a three-dimensional network of relatively strong H-bonds. This challenge is exacerbated in the supercooled state, where the dramatic slowing down of structural dynamics is reminiscent of the, equally poorly understood, generic behavior of liquids near the glass transition temperature. By probing single-molecule dynamics on a wide range of time and length scales, quasielastic neutron scattering (QENS) can potentially reveal the mechanistic details of water's structural dynamics, but because of interpretational ambiguities this potential has not been fully realized. To resolve these issues, we present here an extensive set of high-quality QENS data from water in the range 253-293 K and a corresponding set of molecular dynamics (MD) simulations to facilitate and validate the interpretation. Using a model-free approach, we analyze the QENS data in terms of two motional components. Based on the dynamical clustering observed in MD trajectories, we identify these components with two distinct types of structural dynamics: picosecond local (L) structural fluctuations within dynamical basins and slower interbasin jumps (J). The Q-dependence of the dominant QENS component, associated with J dynamics, can be quantitatively rationalized with a continuous-time random walk (CTRW) model with an apparent jump length that depends on low-order moments of the jump length and waiting time distributions. Using a simple coarse-graining algorithm to quantitatively identify dynamical basins, we map the Newtonian MD trajectory on a CTRW trajectory, from which the jump length and waiting time distributions are computed. The jump length distribution is Gaussian and the rms jump length increases from 1.5 to 1.9 Å as the temperature increases from 253 to 293 K. The rms basin radius increases from 0.71 to 0.75 Å over the same range. The waiting time distribution is exponential at all investigated temperatures, ruling out significant dynamical heterogeneity. However, a simulation at 238 K reveals a small but significant dynamical heterogeneity. The macroscopic diffusion coefficient deduced from the QENS data agrees quantitatively with NMR and tracer results. We compare our QENS analysis with existing approaches, arguing that the apparent dynamical heterogeneity implied by stretched exponential fitting functions results from the failure to distinguish intrabasin (L) from interbasin (J) structural dynamics. We propose that the apparent dynamical singularity at ∼220 K corresponds to freezing out of J dynamics, while the calorimetric glass transition corresponds to freezing out of L dynamics.}},
  author       = {{Qvist, Johan and Schober, Helmut and Halle, Bertil}},
  issn         = {{0021-9606}},
  language     = {{eng}},
  number       = {{14}},
  publisher    = {{American Institute of Physics (AIP)}},
  series       = {{Journal of Chemical Physics}},
  title        = {{Structural dynamics of supercooled water from quasielastic neutron scattering and molecular simulations.}},
  url          = {{http://dx.doi.org/10.1063/1.3578472}},
  doi          = {{10.1063/1.3578472}},
  volume       = {{134}},
  year         = {{2011}},
}