Molecular dynamics simulation of the ice nucleation and growth process leading to water freezing

Upon cooling, water freezes to ice. This familiar phase transition occurs widely in nature, yet unlike the freezing of simple liquids, it has never been successfully simulated on a computer. The difficulty lies with the fact that hydrogen bonding between individual water molecules yields a disordered three-dimensional hydrogen-bond network whose rugged and complex global potential energy surface permits a large number of possible network configurations. As a result, it is very challenging to reproduce the freezing of ‘real’ water into a solid with a unique crystalline structure. For systems with a limited number of possible disordered hydrogen-bond network structures, such as confined water, it is relatively easy to locate a pathway from a liquid state to a crystalline structure. For pure and spatially unconfined water, however, molecular dynamics simulations of freezing are severely hampered by the large number of possible network configurations that exist. Here we present a molecular dynamics trajectory that captures the molecular processes involved in the freezing of pure water. We find that ice nucleation occurs once a sufficient number of relatively long-lived hydrogen bonds develop spontaneously at the same location to form a fairly compact initial nucleus. The initial nucleus then slowly changes shape and size until it reaches a stage that allows rapid expansion, resulting in crystallization of the entire system.

[1]  B. Ahrens,et al.  Bis(diphenylphosphino)methane and related ligands as hydrogen bond donors , 1998 .

[2]  Andersen,et al.  10(6)-particle molecular-dynamics study of homogeneous nucleation of crystals in a supercooled atomic liquid. , 1990, Physical review. B, Condensed matter.

[3]  Sastry,et al.  Singularity-free interpretation of the thermodynamics of supercooled water. , 1998, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  Y. Furukawa,et al.  Anisotropic growth kinetics of ice crystals from water studied by molecular dynamics simulation , 1996 .

[5]  W. Wernsdorfer,et al.  Quantum phase interference and parity effects in magnetic molecular clusters , 1999, Science.

[6]  Hideki Tanaka,et al.  Fluctuation, relaxations, and hydration in liquid water. Hydrogen-bond rearrangement dynamics , 1993 .

[7]  J Tejada Does macroscopic quantum coherence occur in ferritin? , 1996, Science.

[8]  Classical and Quantum Magnetization Reversal Studied in Nanometer-sized Particles and Clusters to Be Published In: Advances in Chemical Physics Invitation from Stuart A. Rice Second Version of 31 Dec. 2000 , 2000 .

[9]  Peter G. Wolynes,et al.  Large local energy fluctuations in water. II. Cooperative motions and fluctuations , 1988 .

[10]  A. Guy Orpen,et al.  Metal-bound chlorine often accepts hydrogen bonds , 1998 .

[11]  Smyth,et al.  Macroscopic quantum tunneling in magnetic proteins. , 1992, Physical review letters.

[12]  Peter D. W. Boyd,et al.  Photosynthetic water oxidation center: spin frustration in distorted cubane MnIVMnIII3 model complexes , 1992 .

[13]  R. Sessoli,et al.  Single-Molecule Magnets , 2000 .

[14]  Friedman,et al.  Macroscopic measurement of resonant magnetization tunneling in high-spin molecules. , 1996, Physical review letters.

[15]  S. Gider,et al.  Classical and quantum magnetic phenomena in natural and artificial ferritin proteins. , 1995, Science.

[16]  L. Thomas,et al.  Macroscopic quantum tunnelling of magnetization in a single crystal of nanomagnets , 1996, Nature.

[17]  M. Weissman 1/f noise and other slow, nonexponential kinetics in condensed matter. , 1988 .

[18]  G. Desiraju The C-H×××O Hydrogen Bond: Structural Implications and Supramolecular Design , 1996 .

[19]  J. Onuchic,et al.  Theory of protein folding: the energy landscape perspective. , 1997, Annual review of physical chemistry.

[20]  J. F. Stoddart,et al.  The magnitude of [C-H...O] hydrogen bonding in molecular and supramolecular assemblies. , 2001, Journal of the American Chemical Society.

[21]  T. Leisner,et al.  Homogeneous nucleation rates of supercooled water measured in single levitated microdroplets , 1999 .

[22]  P. Steinhardt,et al.  Bond-orientational order in liquids and glasses , 1983 .

[23]  R. Ramaswamy,et al.  Long time fluctuation of liquid water: 1/f spectrum of energy fluctuation in hydrogen bond network rearrangement dynamics , 1992 .

[24]  J. Klafter,et al.  The random walk's guide to anomalous diffusion: a fractional dynamics approach , 2000 .

[25]  A. C. Zettlemoyer,et al.  Homogeneous Nucleation Theory , 1974 .

[26]  F. Stillinger,et al.  Packing Structures and Transitions in Liquids and Solids , 1984, Science.

[27]  S. Saito,et al.  Translational and orientational dynamics of a water cluster (H2O)108 and liquid water: Analysis of neutron scattering and depolarized light scattering , 1995 .

[28]  L. Brammer,et al.  Synthesis, Molecular Structures, and Fluxional Behavior of dppm-Bridged Complexes of Platinum(II) with Linear Gold(I), Trigonal Silver(I), or Tetrahedral Mercury(II) Centers , 1996 .

[29]  D. Wales,et al.  From Topographies to Dynamics on Multidimensional Potential Energy Surfaces of Atomic Clusters , 1996, Science.

[30]  M Karplus,et al.  The fundamentals of protein folding: bringing together theory and experiment. , 1999, Current opinion in structural biology.

[31]  H. Stanley,et al.  The relationship between liquid, supercooled and glassy water , 1998, Nature.

[32]  C. Sangregorio,et al.  QUANTUM TUNNELING OF THE MAGNETIZATION IN AN IRON CLUSTER NANOMAGNET , 1997 .

[33]  L. Bartell Nucleation Rates in Freezing and Solid-State Transitions. Molecular Clusters as Model Systems , 1995 .

[34]  Kenichiro Koga,et al.  First-order transition in confined water between high-density liquid and low-density amorphous phases , 2000, Nature.

[35]  B. Monosmith,et al.  Collision-induced Raman scattering from water and aqueous solutions , 1989 .

[36]  I. Ohmine Liquid Water Dynamics: Collective Motions, Fluctuation, and Relaxation , 1995 .

[37]  M. Freytag,et al.  Hydrogen bonds C–H⋯Cl as astructure-determining factor in the gold(I) complexbis(3-bromopyridine)gold(I)dichloroaurate(I) , 2000 .

[38]  M. Sasai Instabilities of hydrogen bond network in liquid water , 1990 .

[39]  P. Kusalik,et al.  Electrofreezing of Liquid Water: A Microscopic Perspective , 1996 .

[40]  A. Caneschi,et al.  Magnetic bistability in a metal-ion cluster , 1993, Nature.

[41]  B. Barbara,et al.  Macroscopic quantum tunneling in antiferromagnets , 1990 .

[42]  Michael N. Leuenberger,et al.  Quantum computing in molecular magnets , 2000, Nature.

[43]  D. Frenkel,et al.  Enhancement of protein crystal nucleation by critical density fluctuations. , 1997, Science.

[44]  M. Matsumoto,et al.  A new approach to the dynamics of hydrogen bond network in liquid water , 1996 .

[45]  N. Dilley,et al.  Half-Integer-Spin Small Molecule Magnet Exhibiting Resonant Magnetization Tunneling , 1998 .