Evidence of a liquid liquid phase transition in H2O and D2O from path-integral molecular dynamics simulations

被引:0
|
作者
Eltareb, Ali [1 ,3 ]
Lopez, Gustavo E. [2 ,4 ]
Giovambattista, Nicolas [1 ,3 ,4 ]
机构
[1] CUNY Brooklyn Coll, Dept Phys, Brooklyn, NY 11210 USA
[2] CUNY, Lehman Coll, Dept Chem, Bronx, NY 10468 USA
[3] CUNY, Grad Ctr, PhD Program Phys, New York, NY 10016 USA
[4] CUNY, Grad Ctr, PhD Program Chem, New York, NY 10016 USA
来源
SCIENTIFIC REPORTS | 2022年 / 12卷 / 01期
基金
美国国家卫生研究院;
关键词
BULK SUPERCOOLED WATER; ISOTHERMAL COMPRESSIBILITY; THERMODYNAMIC PROPERTIES; 1ST-ORDER TRANSITION; CRITICAL-POINT; HEAT-CAPACITY; DENSITY; TEMPERATURE; PRESSURE; SINGULARITY;
D O I
10.1038/s41598-022-09525-x
中图分类号
O [数理科学和化学]; P [天文学、地球科学]; Q [生物科学]; N [自然科学总论];
学科分类号
07 ; 0710 ; 09 ;
摘要
We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H2O and D2O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density rho(T), isothermal compressibility kappa(T)(T), and self-diffusion coefficients D(T) of H2O and D2O are in excellent agreement with available experimental data; the isobaric heat capacity C-P(T) obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H2O and D2O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H2O and D2O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H2O, from PIMD simulations, is located at P-c = 167 +/- 9 MPa, T-c = 159 +/- 6 K, and rho(c) = 1.02 +/- 0.01 g/cm(3). Isotope substitution effects are important; the LLCP location in q-TIP4P/F D2O is estimated to be P-c = 176 +/- 4 MPa, T-c = 177 +/- 2 K, and rho(c) = 1.13 +/- 0.01 g/cm(3). Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, P-c = 203 +/- 4 MPa, T-c = 175 +/- 2 K, and rho(c) = 1.03 +/- 0.01 g/cm(3)). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of T-c for D2O and, particularly, H2O suggest that improved water models are needed for the study of supercooled water.
引用
收藏
页数:14
相关论文
共 50 条
  • [31] Path-integral molecular dynamics simulations of hydrated hydrogen chloride cluster HCl(H2O)4 on a semiempirical potential energy surface
    Takayanagi, Toshiyuki
    Takahashi, Kenta
    Kakizaki, Akira
    Shiga, Motoyuki
    Tachikawa, Masanori
    CHEMICAL PHYSICS, 2009, 358 (03) : 196 - 202
  • [32] Theoretical study on mechanisms of structural rearrangement and ionic dissociation in the HCl(H2O)4 cluster with path-integral molecular dynamics simulations
    Sugawara, Shuichi
    Yoshikawa, Takehiro
    Takayanagi, Toshiyuki
    Tachikawa, Masanori
    CHEMICAL PHYSICS LETTERS, 2011, 501 (4-6) : 238 - 244
  • [33] Vapor phase hydration of glasses in H2O and D2O
    Schatz, TR
    Buechele, AC
    Mooers, CF
    Wysoczanski, R
    Pegg, IL
    ENVIRONMENTAL ISSUES AND WASTE MANAGEMENT TECHNOLOGIES IN THE CERAMIC AND NUCLEAR INDUSTRIES VIII, 2003, 143 : 253 - 262
  • [34] PHASE DIAGRAMS OF H2O AND D2O AT HIGH PRESSURES
    PISTORIUS, CW
    RAPOPORT, E
    CLARK, JB
    JOURNAL OF CHEMICAL PHYSICS, 1968, 48 (12): : 5509 - +
  • [35] Stimulated Raman spectra from H2O, D2O, HDO, and solutions of NaClO4 in H2O and D2O
    Colles, M. J.
    Walrafen, G. E.
    Wecht, K. W.
    CHEMICAL PHYSICS LETTERS, 1970, 4 (10) : 621 - 624
  • [36] CORRESPONDENCE OF FUNDAMENTAL AND COMBINATION BANDS IN INFRARED-SPECTRA OF LIQUID H2O AND D2O
    BONNER, OD
    INFRARED PHYSICS, 1972, 12 (02): : 109 - &
  • [37] THE STRUCTURE OF ABSORPTION BANDS OF LIQUID H2O, D2O AND HDO, DETERMINED BY THE STRUCTURE OF THEIR MOLECULES
    BULANIN, MO
    OPTIKA I SPEKTROSKOPIYA, 1957, 2 (05): : 557 - 561
  • [38] X-ray emission from pure and dilute H2O and D2O in a liquid microjet: Hydrogen bonds and nuclear dynamics
    Lange, Kathrin M.
    Soldatov, Mikhail
    Golnak, Ronny
    Gotz, Malte
    Engel, Nicholas
    Koennecke, Rene
    Rubensson, Jan-Erik
    Aziz, Emad F.
    PHYSICAL REVIEW B, 2012, 85 (15)
  • [39] SOME VIBRATIONAL PROPERTIES OF SOLID AND LIQUID H2O AND D2O DERIVED FROM SLOW NEUTRON MEASUREMENTS
    LARSSON, KE
    DAHLBORG, U
    ARKIV FOR FYSIK, 1962, 22 (05): : 418 - 419
  • [40] Orthogonalyzed H2O and D2O species obtained from infrared spectra of liquid water at several temperatures
    Max, Jean-Joseph
    Larouche, Pascal
    Chapados, Camille
    JOURNAL OF MOLECULAR STRUCTURE, 2017, 1149 : 457 - 472