The interplay between dynamic heterogeneities and structure of bulk liquid water: A molecular dynamics simulation study

Journal Article · Sun Jun 28 00:00:00 EDT 2015 · Journal of Chemical Physics

DOI:https://doi.org/10.1063/1.4922930· OSTI ID:22490849

Demontis, Pierfranco; Suffritti, Giuseppe B. ^[1]; Gulín-González, Jorge ^[2]; Masia, Marco ^[1]; Sant, Marco ^[1]

---

1. Dipartimento di Chimica e Farmacia, Università degli studi di Sassari, Sassari (Italy)
2. 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)

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{sup ∗} ∼ 315 ± 5 K, was spotted at T{sup ∗} ∼ 283 K and T{sup ∗} ∼ 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.

• OSTI ID: 22490849
• Journal Information: Journal of Chemical Physics, Vol. 142, Issue 24; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
• Country of Publication: United States
• Language: English

Related Subjects

37 INORGANIC
ORGANIC
PHYSICAL AND ANALYTICAL CHEMISTRY
97 MATHEMATICAL METHODS AND COMPUTING
ATOMS
COMPUTERIZED SIMULATION
CORRELATION FUNCTIONS
DISTRIBUTION
DISTRIBUTION FUNCTIONS
HYDROGEN
LIQUIDS
MOLECULAR DYNAMICS METHOD
MOLECULES
ORDER PARAMETERS
POTENTIALS
RELAXATION
TEMPERATURE RANGE 0400-1000 K
WATER