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.     

Order parameters for the multistep crystallization of clathrate hydrates

Recent reports indicate that the crystallization of clathrate hydrates occurs in multiple steps that involve amorphous intermediates and metastable clathrate crystals. The elucidation of the reaction coordinate for clathrate crystallization requires the use of order parameters able to identify the reactants, products, and intermediates in the crystallization pathway. Nevertheless, existing order parameters cannot distinguish between amorphous and crystalline clathrates or between different clathrate crystals. In this work, we present the first set of order parameters that discern between the sI and sII clathrate crystals, the amorphous clathrates, the blob of solvent-separated guests and the liquid solution. These order parameters can be used to monitor the advance of the crystallization and for the efficient implementation of methods to sample the rare clathrate nucleation events in molecular simulations. We illustrate the use of these order parameters in the analysis of the growth and the dissolution of clathrate crystals and the spontaneous nucleation and growth of clathrates under conditions of high supercooling.     

The nuclear magnetic resonance line shapes of Xe in the cages of clathrate hydrates

We report, for the first time, a prediction of the line shapes that would be observed in the (129)Xe nuclear magnetic resonance (NMR) spectrum of xenon in the cages of clathrate hydrates. We use the dimer tensor model to represent pairwise contributions to the intermolecular magnetic shielding tensor for Xe at a specific location in a clathrate cage. The individual tensor components from quantum mechanical calculations in clathrate hydrate structure I are represented by contributions from parallel and perpendicular tensor components of Xe-O and Xe-H dimers. Subsequently these dimer tensor components are used to reconstruct the full magnetic shielding tensor for Xe at an arbitrary location in a clathrate cage. The reconstructed tensors are employed in canonical Monte Carlo simulations to find the Xe shielding tensor component along a particular magnetic field direction. The shielding tensor component weighted according to the probability of finding a crystal fragment oriented along this direction in a polycrystalline sample leads to a predicted line shape. Using the same set of Xe-O and Xe-H shielding functions and the same Xe-O and Xe-H potential functions we calculate the Xe NMR spectra of Xe atom in 12 distinct cage types in clathrate hydrates structures I, II, H, and bromine hydrate. Agreement with experimental spectra in terms of the number of unique tensor components and their relative magnitudes is excellent. Agreement with absolute magnitudes of chemical shifts relative to free Xe atom is very good. We predict the Xe line shapes in two cages in which Xe has not yet been observed.     

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.     

Equilibrium conditions for clathrate hydrates formed from methane and aqueous propanol solutions

Equilibrium conditions for clathrate hydrates formed from methane and different concentrations of 1-propanol or 2-propanol aqueous solutions were experimentally determined at temperatures of 274.0–287.1 K and pressures up to 11.0 MPa. Each propanol has an inhibiting and/or promoting effect on hydrate formation depending on the propanol concentration. A structural transition from a structure I to a different hydrate structure occurred at concentrations between 3 and 5 mass% for 1-propanol and between 2 and 3 mass% for 2-propanol.     

Free energies of carbon dioxide sequestration and methane recovery in clathrate hydrates

Classical molecular dynamics simulations are used to compare the stability of methane, carbon dioxide, nitrogen, and mixed CO(2)N(2) structure I (sI) clathrates under deep ocean seafloor temperature and pressure conditions (275 K and 30 MPa) which were considered suitable for CO(2) sequestration. Substitution of methane guests in both the small and large sI cages by CO(2) and N(2) fluids are considered separately to determine the separate contributions to the overall free energy of substitution. The structure I clathrate with methane in small cages and carbon dioxide in large cages is determined to be the most stable. Substitutions of methane in the small cages with CO(2) and N(2) have positive free energies. Substitution of methane with CO(2) in the large cages has a large negative free energy and substitution of the methane in the large cages with N(2) has a small positive free energy. The calculations show that under conditions where storage is being considered, carbon dioxide spontaneously replaces methane from sI clathrates, causing the release of methane. This process must be considered if there are methane clathrates present where CO(2) sequestration is to be attempted. The calculations also indicate that N(2) does not directly compete with CO(2) during methane substitution or clathrate formation and therefore can be used as a carrier gas or may be present as an impurity. Simulations further reveal that the replacement of methane with CO(2) in structure II (sII) cages also has a negative free energy. In cases where sII CO(2) clathrates are formed, only single occupancy of the large cages will be observed.     

Nuclear magnetic resonance relaxation of methane protons in aqueous solutions

Proton nuclear magnetic resonance longitudinal relaxation rates are reported at 100MHz for CH₄ and CHD₃ in deuterium oxide solutions. The results demonstrate that methane reorients with a correlation time on the order of 0.1 psec, considerably faster than molecule correlation times in aqueous solutions.     

The cages, dynamics, and structuring of incipient methane clathrate hydrates

Interest in describing clathrate hydrate formation mechanisms spans multiple fields of science and technical applications. Here, we report findings from multiple molecular dynamics simulations of spontaneous methane clathrate hydrate nucleation and growth from fully demixed and disordered two-phase fluid systems of methane and water. Across a range of thermodynamic conditions and simulation geometries and sizes, a set of seven cage types comprises approximately 95% of all cages formed in the nucleated solids. This set includes the ubiquitous 5(12) cage, the 5(12)6(n) subset (where n ranges from 2-4), and the 4(1)5(10)6(n) subset (where n also ranges from 2-4). Transformations among these cages occur via water pair insertions/removals and rotations, and may elucidate the mechanisms of solid-solid structural rearrangements observed experimentally. Some consistency is observed in the relative abundance of cages among all nucleation trajectories. 5(12) cages are always among the two most abundant cage types in the nucleated solids and are usually the most abundant cage type. In all simulations, the 5(12)6(n) cages outnumber their 4(1)5(10)6(n) counterparts with the same number of water molecules. Within these consistent features, some stochasticity is observed in certain cage ratios and in the long-range ordering of the nucleated solids. Even when comparing simulations performed at the same conditions, some trajectories yield swaths of multiple adjacent sI unit cells and long-range order over 5 nm, while others yield only isolated sI unit cells and little long-range order. The nucleated solids containing long-range order have higher 5(12)6(2)/5(12) and 5(12)6(3)/4(1)5(10)6(2) cage ratios when compared to systems that nucleate with little long-range order. The formation of multiple adjacent unit cells of sI hydrate at high driving forces suggests an alternative or addition to the prevailing hydrate nucleation hypotheses which involve formation through amorphous intermediates.     

Molecular dynamics study of the stability of methane structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) methane clathrate hydrates in a 3 x 3 x 3 sH unit cell replica. Simulations are performed at experimental conditions of 300 K and 2 GPa for three methane intermolecular potentials. The five small cages of the sH unit cell are assigned methane guest occupancies of one and large cage guest occupancies of one to five are considered. Radial distribution functions, unit cell volumes, and configurational energies are studied as a function of large cage CH(4) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies. Large cage occupancy of five is the most stable configuration for a Lennard-Jones united-atom potential and the Tse-Klein-McDonald potential parametrized for condensed methane phases and two for the most stable configuation for the Murad and Gubbins potential.     

Stability and reactivity of methane clathrate hydrates: insights from density functional theory

The sI methane clathrate hydrate consists of methane gas molecules encapsulated as dodecahedron (5(12)CH(4)) and tetrakaidecahedron (5(12)6(2)CH(4)) water cages. The characterization of the stability of these cages is crucial to an understanding of the mechanism of their formation. In the present work, we perform calculations using density functional theory to calculate interaction energies, free energies, and reactivity indices of these cages. The contributions from polarization functions to interaction energies is more than diffuse functions from Pople basis sets, though both functions from the correlation-consistent basis sets contribute significantly to interaction energies. The interaction energies and free energies show that the formation of the 5(12)CH(4) cage (from the 5(12) cage) is more favored compared to the 5(12)6(2)CH(4) cage (from the 5(12)6(2) cage). The pressure-dependent study shows a spontaneous formation of the 5(12)CH(4) cage at 273 K (P ≥ 77 bar) and the 5(12)6(2)CH(4) cage (P = 100 bar). The reactivity of the 5(12)CH(4) cage is similar to that of the 5(12) cage, but the 5(12)6(2)CH(4) cage is more reactive than the 5(12)6(2) cage.     

2H and13C NMR study of guest molecule orientation in clathrate hydrates

At relatively high temperatures (200–270K), clathrate hydrate cages achieve their full crystallographic symmetry because of time averaging of different cage configurations which exist because of disorder in the water molecule orientations. The average orientation of guest molecules in the cages can be obtained from the NMR spectrum, in case of spin 1/2 nuclei from the nuclear shielding tensor, in case of spin 1 nuclei from the quadrupole coupling tensor. Guest molecules studied include carbon dioxide, carbonyl sulphide, methyl-d₃ fluoride, methyl-d₃ chloride, methyl-d₃ bromide, ethane-d₆, acetylene-d₂ in the structure I hydrates, and methyl-d₃ iodide in the structure II hydrate.For the slightly flattened large cage of structure I hydrate, the guest molecules rotate so that the plane which contains the long axis of the molecule is confined to be nearer to the equatorial plane of the cage than the axial regions.Since the structure II large cage has tetrahedral symmetry on time average, it exerts no orienting effect on guest molecules.NRCC no. 32722.     

Intriguing electric field effect on magnetic spin couplings in dielectron clathrate hydrates

Clathrate hydrates have appeared as promising icy materials as the radical, high-spin molecule, and even electron clathrate hydrates are found. In particular, dielectron clathrate hydrates are expected to develop as structural units for a novel class of icy magnetic materials because of not only possible spin coupling interaction, but also very sensitive response to electric field of the loosely bound electrons. However, electric field responses concerning the magnetic properties of such hydrates have not been reported so far. In this work, three representative dielectron clathrate hydrate model clusters (e₂@4⁶6⁸BB, e₂@5¹²6²BB, and e₂@4⁶6⁸AB) were considered for the exploration of their magnetic spin coupling properties, electron distributions, and energy responses to applied electric field. The results calculated at the density functional theory level show that the energies and electron spin coupling properties of these dielectron clathrate hydrate clusters are quite sensitive to applied electric field, presenting intriguing variations. Most importantly, applied electric field can regulate the strength of spin coupling between two trapped electrons, and even could realize the magnetic interconversion of such dielectron cluster structures between antiferromagnetic and paramagnetic or diamagnetic characteristics. Clearly, the intriguing variations should be attributed to the diffuse character, special mobility and polarizable properties of such trapped electrons, and especially the susceptible redistributions of two electrons (including the electron cloud shape and distance between two electron centers) to the electric field. This work opens up the possibility of designing novel icy magnetic materials with sensitive electric field responses of the magnetic properties.     

Thermodynamic modeling for clathrate hydrates of ozone

We report a theoretical study to predict the phase-equilibrium properties of ozone-containing clathrate hydrates based on the statistical thermodynamics model developed by van der Waals and Platteeuw. The Patel–Teja–Valderrama equation of state is employed for an accurate estimation of the properties of gas phase ozone. We determined the three parameters of the Kihara intermolecular potential for ozone as a = 6.815 · 10⁻² nm, σ = 2.9909 · 10⁻¹ nm, and ε · k_B⁻¹ = 184.00 K. An infinite set of ε–σ parameters for ozone were determined, reproducing the experimental phase equilibrium pressure–temperature data of the (O₃ + O₂ + CO₂) clathrate hydrate. A unique parameter pair was chosen based on the experimental ozone storage capacity data for the (O₃ + O₂ + CCl₄) hydrate that we reported previously. The prediction with the developed model showed good agreement with the experimental phase equilibrium data within ±2% of the average deviation of the pressure. The Kihara parameters of ozone showed slightly better suitability for the structure-I hydrate than CO₂, which was used as a help guest. Our model suggests the possibility of increasing the ozone storage capacity of clathrate hydrates (∼7% on a mass basis) from the previously reported experimental capacity (∼1%).     

Polymorphism in Br2 clathrate hydrates

The structure and composition of bromine clathrate hydrate has been controversial for more than 170 years due to the large variation of its observed stoichiometries. Several different crystal structures were proposed before 1997 when Udachin et al. (Udachin, K. A.; Enright, G. D.; Ratcliffe, C. I.; Ripmeester, J. A. J. Am. Chem. Soc. 1997, 119, 11481) concluded that Br2 forms only the tetragonal structure (TS-I). We show polymorphism in Br2 clathrate hydrates by identifying two distinct crystal structures through optical microscopy and resonant Raman spectroscopy on single crystals. After growing TS-I crystals from a liquid bromine-water solution, upon dropping the temperature slightly below -7 degrees C, new crystals of cubic morphology form. The new crystals, which have a limited thermal stability range, are assigned to the CS-II structure. The two structures are clearly distinguished by the resonant Raman spectra of the enclathrated Br2, which show long overtone progressions and allow the extraction of accurate vibrational parameters: omega(e) = 321.2 +/- 0.1 cm(-1) and omega(e)x(e) = 0.82 +/- 0.05 cm(-1) in TS-I and omega(e) = 317.5 +/- 0.1 cm(-1) and omega(e)x(e) = 0.70 +/- 0.1 cm(-1) in CS-II. On the basis of structural analysis, the discovery of the CS-II crystals implies stability of a large class of bromine hydrate structures and, therefore, polymorphism.     

Raman spectra of vibrational and librational modes in methane clathrate hydrates using density functional theory

The sI type methane clathrate hydrate lattice is formed during the process of nucleation where methane gas molecules are encapsulated in the form of dodecahedron (5(12)CH(4)) and tetrakaidecahedron (5(12)6(2)CH(4)) water cages. The characterization of change in the vibrational modes which occur on the encapsulation of CH(4) in these cages plays a key role in understanding the formation of these cages and subsequent growth to form the hydrate lattice. In this present work, we have chosen the density functional theory (DFT) using the dispersion corrected B97-D functional to characterize the Raman frequency vibrational modes of CH(4) and surrounding water molecules in these cages. The symmetric and asymmetric C-H stretch in the 5(12)CH(4) cage is found to shift to higher frequency due to dispersion interaction of the encapsulated CH(4) molecule with the water molecules of the cages. However, the symmetric and asymmetric O-H stretch of water molecules in 5(12)CH(4) and 5(12)6(2)CH(4) cages are shifted towards lower frequency due to hydrogen bonding, and interactions with the encapsulated CH(4) molecules. The CH(4) bending modes in the 5(12)CH(4) and 5(12)6(2)CH(4) cages are blueshifted, though the magnitude of the shifts is lower compared to modes in the high frequency region which suggests bending modes are less affected on encapsulation of CH(4). The low frequency librational modes which are collective motion of the water molecules and CH(4) in these cages show a broad range of frequencies which suggests that these modes largely contribute to the formation of the hydrate lattice.     

Molecular dynamics study of structure H clathrate hydrates of methane and large guest molecules

Methane storage in structure H (sH) clathrate hydrates is attractive due to the relatively higher stability of sH as compared to structure I methane hydrate. The additional stability is gained without losing a significant amount of gas storage density as happens in the case of structure II (sII) methane clathrate. Our previous work has showed that the selection of a specific large molecule guest substance (LMGS) as the sH hydrate former is critical in obtaining the optimum conditions for crystallization kinetics, hydrate stability, and methane content. In this work, molecular dynamics simulations are employed to provide further insight regarding the dependence of methane occupancy on the type of the LMGS and pressure. Moreover, the preference of methane molecules to occupy the small (5(12)) or medium (4(3)5(6)6(3)) cages and the minimum cage occupancy required to maintain sH clathrate mechanical stability are examined. We found that thermodynamically, methane occupancy depends on pressure but not on the nature of the LMGS. The experimentally observed differences in methane occupancy for different LMGS may be attributed to the differences in crystallization kinetics and/or the nonequilibrium conditions during the formation. It is also predicted that full methane occupancies in both small and medium clathrate cages are preferred at higher pressures but these cages are not fully occupied at lower pressures. It was found that both small and medium cages are equally favored for occupancy by methane guests and at the same methane content, the system suffers a free energy penalty if only one type of cage is occupied. The simulations confirm the instability of the hydrate when the small and medium cages are empty. Hydrate decomposition was observed when less than 40% of the small and medium cages are occupied.