Structural and dynamical properties of methane clathrate hydrates

Equilibrium molecular dynamics (MD) simulations have been performed in both the NVT and NPT ensembles to study the structural and dynamical properties of fully occupied methane clathrate hydrates at 50, 125, and 200 K. Five atomistic potential models were used for water, ranging from fully flexible to rigid polarizable and nonpolarizable. A flexible and a rigid model were utilized for methane. The phonon densities of states were evaluated and the localized rattling modes for the methane molecules were found to couple to the acoustic phonons of the host lattice. The calculated methane density of states was found to be in reasonable agreement with available experimental data. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1569–1581, 2003

[1]  J. Bertie,et al.  Far infrared absorption and rotational vibrations of the guest molecules in structure I clathrate hydrates between 4.3 and 100 K , 1977 .

[2]  K. Kvenvolden.  Methane hydrate — A major reservoir of carbon in the shallow geosphere? , 1988 .

[3]  J. D. Bernal,et al.  A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions , 1933 .

[4]  I. R. Mcdonald,et al.  Molecular dynamics studies of ice Ic and the structure I clathrate hydrate of methane , 1983 .

[5]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[6]  C. Brooks.  Computer simulation of liquids , 1989 .

[7]  F. Trouw,et al.  Coupling of localized guest vibrations with the lattice modes in clathrate hydrates , 2001 .

[8]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[9]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[10]  C. Gutt,et al.  Quantum rotations in natural methane-clathrates from the Pacific sea-floor , 1999 .

[11]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[12]  Rahman,et al.  Molecular-dynamics study of atomic motions in water. , 1985, Physical review. B, Condensed matter.

[13]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[14]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

[15]  Peter A. Kollman,et al.  Simulation of the solvation free energies for methane, ethane, and propane and corresponding amino acid dipeptides : a critical test of the bond-PMF correction, a new set of hydrocarbon parameters, and the gas phase-water hydrophobicity scale , 1992 .

[16]  B. Berne,et al.  Combined fluctuating charge and polarizable dipole models: Application to a five-site water potential function , 2001 .

[17]  Pressure calculation in polar and charged systems using Ewald summation: Results for the extended simple point charge model of water , 1998, physics/9806038.

[18]  J. Tse,et al.  The low frequency vibrations in clathrate hydrates , 1997 .

[19]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[20]  John S. Tse,et al.  Thermal expansion of structure I ethylene oxide hydrate , 1987 .

[21]  G. MacDonald.  The Future of Methane as an Energy Resource , 1990 .

[22]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[23]  Michael L. Klein,et al.  Computer simulation studies of the structure I clathrate hydrates of methane, tetrafluoromethane, cyclopropane, and ethylene oxide , 1984 .

[24]  G. A. Jeffrey,et al.  Polyhedral Clathrate Hydrates. IX. Structure of Ethylene Oxide Hydrate , 1965 .

[25]  J. Bertie,et al.  Far‐infrared spectrum and x‐ray diffraction of polycrystalline structure I clathrate hydrate of xenon at 4.3 K , 1982 .

[26]  Franklin M. Orr,et al.  Gas Hydrate Formation in the Deep Sea: In Situ Experiments with Controlled Release of Methane, Natural Gas, and Carbon Dioxide , 1998 .

[27]  G. Ciccotti,et al.  Hoover NPT dynamics for systems varying in shape and size , 1993 .

[28]  J. Tse.  Dynamical Properties and Stability of Clathrate Hydrates , 1994 .

[29]  Niels Grønbech-Jensen.  Lekner Summation of Long Range Interactions in Periodic Systems , 1997 .

[30]  P. Rodger.  Lattice relaxation in type I gas hydrates , 1991 .

[31]  P. Cummings,et al.  Molecular dynamics study of the structure and thermophysical properties of model sI clathrate hydrates , 2002 .

[32]  Niall J. English,et al.  Atomistic simulations of liquid water using Lekner electrostatics , 2002 .

[33]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[34]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[35]  John Lekner,et al.  Summation of dipolar fields in simulated liquid-vapour interfaces , 1989 .

[36]  V. F. Sears,et al.  Rotational and Translational Motions of Trapped Methane. Incoherent Inelastic Neutron Scattering of Methane Hydrate , 1997 .

[37]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[38]  John Lekner,et al.  Summation of Coulomb fields in computer-simulated disordered systems , 1991 .