Computer simulation study of metastable ice VII and amorphous phases obtained by its melting

Jan Slovák, Hideki Tanaka

Research output: Contribution to journalArticle

11 Citations (Scopus)

Abstract

Molecular dynamics simulations of metastable ice VII and cubic ice Ic are carried out in order to examine (1) the ability of commonly used water interaction potentials to reproduce the properties of ices, and (2) the possibility of generating low-density amorphous (LDA) structures by heating ice VII, which is known to transform to LDA at ~135 K at normal pressure [S. Klotz, J. M. Besson, G. Hamel, R. J. Nelmes, J. S. Loveday, and W. G. Marshall, Nature (London) 398, 681 (1999)]. We test four simple empirical interaction potentials of water: TIP4P [W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983)], SPC/E [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. B 91, 6269 (1987)], TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000)], and ST2 [F. H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)]. We have found that TIP5P ice VII melts at 210 K, TIP4P at 90 K, and SPC/E at 70 K. Only TIP5P water after transition has a structure similar to that of LDA. TIP4P and SPC/E have almost identical structures, dissimilar to any known water or amorphous phases, but upon heating both slowly evolve towards LDA-like structure. ST2 ice VII is remarkably stable up to 430 K. TIP4P and SPC/E predict correctly the cubic ice collapse into a high-density amorphous ice (HDA) at ~1 GPa whereas TIP5P remains stable up to ~5 GPa. The densities of the simulated ice phases differ significantly, depending on the potential used, and are generally higher than experimental values. The importance of proper treatment of long-range electrostatic interactions is also discussed.

Original languageEnglish
Article number204512
JournalThe Journal of Chemical Physics
Volume122
Issue number20
DOIs
Publication statusPublished - May 22 2005

Fingerprint

Ice
ice
Melting
computerized simulation
melting
Computer simulation
Water
water
Heating
heating
interactions
Coulomb interactions
Molecular dynamics
electrostatics
molecular dynamics

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

Computer simulation study of metastable ice VII and amorphous phases obtained by its melting. / Slovák, Jan; Tanaka, Hideki.

In: The Journal of Chemical Physics, Vol. 122, No. 20, 204512, 22.05.2005.

Research output: Contribution to journalArticle

@article{f113727857cd4bb6a50b853a96d886e2,
title = "Computer simulation study of metastable ice VII and amorphous phases obtained by its melting",
abstract = "Molecular dynamics simulations of metastable ice VII and cubic ice Ic are carried out in order to examine (1) the ability of commonly used water interaction potentials to reproduce the properties of ices, and (2) the possibility of generating low-density amorphous (LDA) structures by heating ice VII, which is known to transform to LDA at ~135 K at normal pressure [S. Klotz, J. M. Besson, G. Hamel, R. J. Nelmes, J. S. Loveday, and W. G. Marshall, Nature (London) 398, 681 (1999)]. We test four simple empirical interaction potentials of water: TIP4P [W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983)], SPC/E [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. B 91, 6269 (1987)], TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000)], and ST2 [F. H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)]. We have found that TIP5P ice VII melts at 210 K, TIP4P at 90 K, and SPC/E at 70 K. Only TIP5P water after transition has a structure similar to that of LDA. TIP4P and SPC/E have almost identical structures, dissimilar to any known water or amorphous phases, but upon heating both slowly evolve towards LDA-like structure. ST2 ice VII is remarkably stable up to 430 K. TIP4P and SPC/E predict correctly the cubic ice collapse into a high-density amorphous ice (HDA) at ~1 GPa whereas TIP5P remains stable up to ~5 GPa. The densities of the simulated ice phases differ significantly, depending on the potential used, and are generally higher than experimental values. The importance of proper treatment of long-range electrostatic interactions is also discussed.",
author = "Jan Slov{\'a}k and Hideki Tanaka",
year = "2005",
month = "5",
day = "22",
doi = "10.1063/1.1906217",
language = "English",
volume = "122",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "20",

}

TY - JOUR

T1 - Computer simulation study of metastable ice VII and amorphous phases obtained by its melting

AU - Slovák, Jan

AU - Tanaka, Hideki

PY - 2005/5/22

Y1 - 2005/5/22

N2 - Molecular dynamics simulations of metastable ice VII and cubic ice Ic are carried out in order to examine (1) the ability of commonly used water interaction potentials to reproduce the properties of ices, and (2) the possibility of generating low-density amorphous (LDA) structures by heating ice VII, which is known to transform to LDA at ~135 K at normal pressure [S. Klotz, J. M. Besson, G. Hamel, R. J. Nelmes, J. S. Loveday, and W. G. Marshall, Nature (London) 398, 681 (1999)]. We test four simple empirical interaction potentials of water: TIP4P [W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983)], SPC/E [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. B 91, 6269 (1987)], TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000)], and ST2 [F. H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)]. We have found that TIP5P ice VII melts at 210 K, TIP4P at 90 K, and SPC/E at 70 K. Only TIP5P water after transition has a structure similar to that of LDA. TIP4P and SPC/E have almost identical structures, dissimilar to any known water or amorphous phases, but upon heating both slowly evolve towards LDA-like structure. ST2 ice VII is remarkably stable up to 430 K. TIP4P and SPC/E predict correctly the cubic ice collapse into a high-density amorphous ice (HDA) at ~1 GPa whereas TIP5P remains stable up to ~5 GPa. The densities of the simulated ice phases differ significantly, depending on the potential used, and are generally higher than experimental values. The importance of proper treatment of long-range electrostatic interactions is also discussed.

AB - Molecular dynamics simulations of metastable ice VII and cubic ice Ic are carried out in order to examine (1) the ability of commonly used water interaction potentials to reproduce the properties of ices, and (2) the possibility of generating low-density amorphous (LDA) structures by heating ice VII, which is known to transform to LDA at ~135 K at normal pressure [S. Klotz, J. M. Besson, G. Hamel, R. J. Nelmes, J. S. Loveday, and W. G. Marshall, Nature (London) 398, 681 (1999)]. We test four simple empirical interaction potentials of water: TIP4P [W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983)], SPC/E [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. B 91, 6269 (1987)], TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000)], and ST2 [F. H. Stillinger and A. Rahman, J. Chem. Phys. 60, 1545 (1974)]. We have found that TIP5P ice VII melts at 210 K, TIP4P at 90 K, and SPC/E at 70 K. Only TIP5P water after transition has a structure similar to that of LDA. TIP4P and SPC/E have almost identical structures, dissimilar to any known water or amorphous phases, but upon heating both slowly evolve towards LDA-like structure. ST2 ice VII is remarkably stable up to 430 K. TIP4P and SPC/E predict correctly the cubic ice collapse into a high-density amorphous ice (HDA) at ~1 GPa whereas TIP5P remains stable up to ~5 GPa. The densities of the simulated ice phases differ significantly, depending on the potential used, and are generally higher than experimental values. The importance of proper treatment of long-range electrostatic interactions is also discussed.

UR - http://www.scopus.com/inward/record.url?scp=20844451050&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=20844451050&partnerID=8YFLogxK

U2 - 10.1063/1.1906217

DO - 10.1063/1.1906217

M3 - Article

VL - 122

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 20

M1 - 204512

ER -