A powerful computational crystallography method to study ice polymorphism
Abstract
Classical molecular dynamics (MD) simulations are employed as a tool to investigate structural properties of ice crystals under several temperature and pressure conditions. All ice crystal phases are analyzed by means of a computational protocol based on a clustering approach following standard MD simulations. The MD simulations are performed by using a recently published classical interaction potential for oxygen and hydrogen in bulk water, derived from neutron scattering data, able to successfully describe complex phenomena such as proton hopping and bond formation/breaking. The present study demonstrates the ability of the interaction potential model to well describe most ice structures found in the phase diagram of water and to estimate the relative stability of 16 known phases through a cluster analysis of simulated powder diagrams of polymorphs obtained from MD simulations. The proposed computational protocol is suited for automated crystal structure identification.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- May 2011
- DOI:
- 10.1063/1.3593200
- arXiv:
- arXiv:1105.1114
- Bibcode:
- 2011JChPh.134t4506C
- Keywords:
-
- bonds (chemical);
- computerised tomography;
- crystal structure;
- ice;
- molecular dynamics method;
- neutron diffraction;
- phase diagrams;
- polymorphism;
- 61.66.Fn;
- 81.30.Dz;
- Inorganic compounds;
- Phase diagrams of other materials;
- Condensed Matter - Materials Science
- E-Print:
- RevTex 4.1, 7 figures - to be published in the Journal of Chemical Physics