Journal of Chemical Physics, Vol.106, No.10, 4187-4195, 1997
Structure-H Clathrate Unit-Cell Coordinates and Simulation of the Structure-H Crystal Interface with Water
The Cartesian coordinates for the structure H clathrate unit cell have been determined to within 5% error. The unit cell, consisting of 34 water molecules and 6 hydrocarbon molecules, was replicated to create a system amenable to molecular dynamics simulation. Radial distribution functions were determined for this crystal. The method of molecular dynamics was employed to model the interface between the crystal and liquid water. It was found that the interfacial behavior spanned a region of 15 Angstrom on both sides of the interface. Both orientational and translational motion in the liquid phase became more vibrational as the interface was approached, while the temperature profile remained constant within about 15 degrees across the interface. In addition, there was a small tendency for the molecules in the liquid phase of the interfacial region to align themselves in an ordered, solidlike configuration. The solid region within about 10 Angstrom of the interface demonstrated persistent steady-state disorder compared with the solid bulk phase.