A molecular dynamics study of ethanol-water hydrogen bonding in binary structure I clathrate hydrate with CO2

Guest-host hydrogen bonding in clathrate hydrates occurs when in addition to the hydrophilic moiety which causes the molecule to form hydrates under high pressure-low temperature conditions, the guests contain a hydrophilic, hydrogen bonding functional group. In the presence of carbon dioxide, ethanol clathrate hydrate has been synthesized with 10% of large structure I (sI) cages occupied by ethanol. In this work, we use molecular dynamics simulations to study hydrogen bonding structure and dynamics in this binary sI clathrate hydrate in the temperature range of 100-250 K. We observe that ethanol forms long-lived (>500 ps) proton-donating and accepting hydrogen bonds with cage water molecules from both hexagonal and pentagonal faces of the large cages while maintaining the general cage integrity of the sI clathrate hydrate. The presence of the nondipolar CO(2) molecules stabilizes the hydrate phase, despite the strong and prevalent alcohol-water hydrogen bonding. The distortions of the large cages from the ideal form, the radial distribution functions of the guest-host interactions, and the ethanol guest dynamics are characterized in this study. In previous work through dielectric and NMR relaxation time studies, single crystal x-ray diffraction, and molecular dynamics simulations we have observed guest-water hydrogen bonding in structure II and structure H clathrate hydrates. The present work extends the observation of hydrogen bonding to structure I hydrates.     

Molecular dynamics simulation of NMR powder lineshapes of linear guests in structure I clathrate hydrates

We perform molecular dynamics simulations (up to 6 ns) for the structure I clathrate hydrates of linear molecules CS, CS(2), OCS, and C(2)H(2) in large cages at different temperatures in the stability range to determine the angular distribution and dynamics of the guests in the large cages. The long axes of linear guest molecules in the oblate large structure I clathrate hydrate cages are primarily confined near the equatorial plane of the cage rather than axial regions. This non-uniform spatial distribution leads to well-known anisotropic lineshapes in the solid-state NMR spectra of the guest species. We use the dynamic distribution of guest orientations in the cages during the MD simulations at different temperatures to predict the (13)C NMR powder lineshapes of the guests in the large cages. The length of the guests and intermolecular interactions of the guests in the water cages determine the angular distribution and the mobility of the guests in the sI large cages at different temperatures. At low temperatures the range of motion of the guests in the cages are limited and this is reflected in the skew of the predicted (13)C lineshapes. As the guest molecules reach the fast motion limit at higher temperatures, the lineshapes for CS, OCS, and C(2)H(2) are predicted to have the "standard" powder lineshapes of guest molecules.     

Determination of NMR Lineshape Anisotropy of Guest Molecules within Inclusion Complexes from Molecular Dynamics Simulations

Nonspherical cages in inclusion compounds can result in non‐uniform motion of guest species in these cages and anisotropic lineshapes in NMR spectra of the guest. Herein, we develop a methodology to calculate lineshape anisotropy of guest species in cages based on molecular dynamics simulations of the inclusion compound. The methodology is valid for guest atoms with spin 1/2 nuclei and does not depend on the temperature and type of inclusion compound or guest species studied. As an example, the nonspherical shape of the structure I (sI) clathrate hydrate large cages leads to preferential alignment of linear CO₂ molecules in directions parallel to the two hexagonal faces of the cages. The angular distribution of the CO₂ guests in terms of a polar angle θ and azimuth angle ? and small amplitude vibrational motions in the large cage are characterized by molecular dynamics simulations at different temperatures in the stability range of the CO₂ sI clathrate. The experimental ¹³C NMR lineshapes of CO₂ guests in the large cages show a reversal of the skew between the low temperature (77 K) and the high temperature (238 K) limits of the stability of the clathrate. We determine the angular distributions of the guests in the cages by classical MD simulations of the sI clathrate and calculate the ¹³C NMR lineshapes over a range of temperatures. Good agreement between experimental lineshapes and calculated lineshapes is obtained. No assumptions regarding the nature of the guest motions in the cages are required.     

Effect of small cage guests on hydrogen bonding of tetrahydrofuran in binary structure II clathrate hydrates

Molecular dynamics simulations of the pure structure II tetrahydrofuran clathrate hydrate and binary structure II tetrahydrofuran clathrate hydrate with CO(2), CH(4), H(2)S, and Xe small cage guests are performed to study the effect of the shape, size, and intermolecular forces of the small cages guests on the structure and dynamics of the hydrate. The simulations show that the number and nature of the guest in the small cage affects the probability of hydrogen bonding of the tetrahydrofuran guest with the large cage water molecules. The effect on hydrogen bonding of tetrahydrofuran occurs despite the fact that the guests in the small cage do not themselves form hydrogen bonds with water. These results indicate that nearest neighbour guest-guest interactions (mediated through the water lattice framework) can affect the clathrate structure and stability. The implications of these subtle small guest effects on clathrate hydrate stability are discussed.     

Thermal expansivity of tetrahydrofuran clathrate hydrate with diatomic guest molecules

The guest dynamics and thermal behavior occurring in the cages of clathrate hydrates appear to be too complex to be clearly understood through various structural and spectroscopic approaches, even for the well-known structures of sI, sII, and sH. Neutron diffraction studies have recently been carried out to clarify the special role of guests in expanding the host water lattices and have contributed to revealing the influence factors on thermal expansivity. Through this letter we attempt to address three noteworthy features occurring in guest inclusion: (1) the effect of guest dimension on host water lattice expansion; (2) the effect of thermal history on host water lattice expansion; and (3) the effect of coherent/incoherent scattering cross sections on guest thermal patterns. The diatomic guests of H 2, D 2, N 2, and O 2 have been selected for study, and their size and mass dependence on the degree of lattice expansion have been examined, and four sII clathrate hydrates with tetrahydrofuran (THF) have been synthesized in order to determine their neutron powder diffraction patterns. After thermal cycling, the THF + H 2 clathrate hydrate is observed to exhibit an irreversible plastic deformation-like pattern, implying that the expanded lattices fail to recover the original state by contraction. The host-water cage dimension after degassing the guest molecules remains as it was expanded, and thus host-guest as well as guest-guest interactions will be altered if guest uptake reoccurs.     

Spectroscopic identification of the mixed hydrogen and carbon dioxide clathrate hydrate

In this contribution, we first found the novel clathrate hydrate containing two gaseous guests of hydrogen and carbon dioxide by spectroscopic analysis. X-ray powder diffraction and NMR spectroscopy were used to identify structure and guest distribution of the mixed H2 + CO2 hydrate. X-ray diffraction result confirmed that the unit cell parameter was 11.8602 +/- 0.0010 A, and the formed hydrate was identified as structure I hydrate. 1H magic angle spinning (MAS) NMR and 13C cross-polarization (CP) NMR spectroscopy were used to examine the distribution of hydrogen and carbon dioxide molecules in the cages of structure I, respectively. These NMR spectra showed that carbon dioxide molecules occupied both small 512 cages and large 51262 cages, and hydrogen molecules only were occluded in small 512 cages of structure I. The new finding of the mixed hydrogen hydrate is expected to contribute toward the development of hydrogen production technology and, particularly, inclusion chemistry.     

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.     

Dynamical properties of hydrogen sulphide motion in its clathrate hydrate from ab initio and classical isobaric-isothermal molecular dynamics

Equilibrium ab initio (AI) and classical constant pressure-constant temperature molecular dynamics (MD) simulations have been performed to investigate the dynamical properties in (type I) hydrogen sulphide hydrate at 150 and 300 K and 1 bar, and also lower temperatures, with particular scrutiny of guest motion. The rattling motions of the guests in the large and small cavities were around 45 and 75-80 cm(-1), respectively, from AIMD, with the corresponding classical MD modes being 10-12 cm(-1) less at 150 K and around 5 cm(-1) lower at 300 K. The rattling motion in the small cavity overlapped somewhat with the translational motion of the host lattice (with modes at circa 85 and 110 cm(-1)), due in part to a smaller cage radius and more frequent occurrences of guest-host hydrogen bonding leading to greater coupling in the motion. The experimentally determined H-S stretch and H-S-H bending frequencies, in the vicinity of 2550-2620 and 1175 cm(-1) [H.R. Richardson et al., J. Chem. Phys.1985, 83, 4387] were reproduced successfully in the AIMD simulations. Consideration of Kubic harmonics for the guest molecules from AIMD revealed that a preferred orientation of the dipole-vector (or C(2)-axis) exists at 150 K vis--vis the [100] cube axis in both the small and large cavities, but is markedly more significant for the small cavity, while there is no preferred orientation at 300 K. In comparison, classical MD did not reveal any preferred orientation at either temperature, or at 75 K (closer to the AIMD simulation at 150 K vis-à-vis that approach's estimated melting point). Probing rotational dynamics of the guests reinforced this temperature effect, revealing more rapid rotational time scales at 300 K with faster decay times of dipole-vector (C(2)) and H-H-vector (C(2)(y)) being similar for each cage, at around 0.25 and 0.2 ps, respectively, versus approximately 0.45 and 0.5 ps (large) and 0.8 ps (small) at 150 K. It was found that the origin of the observed preferred orientations, especially in the small cages, at 150 K via AIMD was attributable to optimization of the dipolar interaction between the guest and outward-pointing water dipoles in the cavity, with guests "flitting" rotationally between various such configurations, forming occasionally hydrogen bonds with the host molecules.     

13C chemical shifts of propane molecules encaged in structure II clathrate hydrate

Experimental NMR measurements for (13)C chemical shifts of propane molecules encaged in 16-hedral cages of structure II clathrate hydrate were conducted to investigate the effects of guest-host interaction of pure propane clathrate on the (13)C chemical shifts of propane guests. Experimental (13)C NMR measurements revealed that the clathrate hydration of propane reverses the (13)C chemical shifts of methyl and methylene carbons in propane guests to gaseous propane at room temperature and atmospheric pressure or isolated propane, suggesting a change in magnetic environment around the propane guest by the clathrate hydration. Inversion of the (13)C chemical shifts of propane clathrate suggests that the deshielding effect of the water cage on the methyl carbons of the propane molecule encaged in the 16-hedral cage is greater than that on its methylene carbon.     

Stability of rare gas structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) clathrate hydrates with the rare gases Ne, Ar, Kr, and Xe. Simulations on a 3 x 3 x 3 sH unit cell replica are performed at ambient pressure at 40 and 100 K temperatures. The small and medium (s+m) cages of the sH unit cell are assigned rare gas guest occupancies of 1 and for large (l) cages guest occupancies of 1-6 are considered. Radial distribution functions for guest pairs with occupancies in the l-l, l-(s+m), and (s+m)-(s+m) cages are presented. The unit cell volumes and configurational energies are studied as a function of large cage occupancy for the rare gases. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at 100 K and 1 bar and 20 kbar pressures. These studies show that the most stable argon clathrate has five guests in the large cages. For krypton and xenon the most stable configurations have three and two guests in the large cages, respectively.     

13C NMR studies of hydrocarbon guests in synthetic structure H gas hydrates: experiment and computation

(13)C NMR chemical shifts were measured for pure (neat) liquids and synthetic binary hydrate samples (with methane help gas) for 2-methylbutane, 2,2-dimethylbutane, 2,3-dimethylbutane, 2-methylpentane, 3-methylpentane, methylcyclopentane, and methylcyclohexane and ternary structure H (sH) clathrate hydrates of n-pentane and n-hexane with methane and 2,2-dimethylbutane, all of which form sH hydrates. The (13)C chemical shifts of the guest atoms in the hydrate are different from those in the free form, with some carbon atoms shifting specifically upfield. Such changes can be attributed to conformational changes upon fitting the large guest molecules in hydrate cages and/or interactions between the guests and the water molecules of the hydrate cages. In addition, powder X-ray diffraction measurements revealed that for the hexagonal unit cell, the lattice parameter along the a-axis changes with guest hydrate former molecule size and shape (in the range of 0.1 ?) but a much smaller change in the c-axis (in the range of 0.01 ?) is observed. The (13)C NMR chemical shifts for the pure hydrocarbons and all conformers were calculated using the gauge invariant atomic orbital method at the MP2/6-311+G(2d,p) level of theory to quantify the variation of the chemical shifts with the dihedral angles of the guest molecules. Calculated and measured chemical shifts are compared to determine the relative contribution of changes in the conformation and guest-water interactions to the change in chemical shift of the guest upon clathrate hydrate formation. Understanding factors that affect experimental chemical shifts for the enclathrated hydrocarbons will help in assigning spectra for complex hydrates recovered from natural sites.     

Molecular dynamics simulations of the mechanisms of thermal conduction in methane hydrates

The thermal conductivity of methane hydrate is an important physical parameter affecting the processes of methane hydrate exploration,mining,gas hydrate storage and transportation as well as other applications.Equilibrium molecular dynamics simulations and the Green-Kubo method have been employed for systems from fully occupied to vacant occupied sI methane hydrate in order to estimate their thermal conductivity.The estimations were carried out at temperatures from 203.15 to 263.15 K and at pressures from 3 to 100 MPa.Potential models selected for water were TIP4P,TIP4P-Ew,TIP4P/2005,TIP4P-FQ and TIP4P/Ice.The effects of varying the ratio of the host and guest molecules and the external thermobaric conditions on the thermal conductivity of methane hydrate were studied.The results indicated that the thermal conductivity of methane hydrate is essentially determined by the cage framework which constitutes the hydrate lattice and the cage framework has only slightly higher thermal conductivity in the presence of the guest molecules.Inclusion of more guest molecules in the cage improves the thermal conductivity of methane hydrate.It is also revealed that the thermal conductivity of the sI hydrate shows a similar variation with temperature.Pressure also has an effect on the thermal conductivity,particularly at higher pressures.As the pressure increases,slightly higher thermal conductivities result.Changes in density have little impact on the thermal conductivity of methane hydrate.     

甲烷水合物导热机理的分子动力学模拟

甲烷水合物导热系数是甲烷水合物勘探、开采、储运以及其他应用过程中一个十分重要的物理参数.我们采用平衡分子动力学（EMD）方法Green-Kubo理论计算温度203.15～263.15K、压力范围3～100MPa、晶穴占有率为0～1的sI甲烷水合物的导热系数,采用的水分子模型包括TIP4P、TIP4P-Ew、TIP4P-FQ、TIP4P/2005、TIP4P/Ice.研究了主客体分子、外界温压条件等对甲烷水合物导热性能的影响.研究结果显示甲烷水合物的低导热性能由主体分子构建的sI笼型结构决定,而客体分子进入笼型结构后,使得笼型结构导热性能增强,同时进入笼型结构的客体分子越多,甲烷水合物导热性能越强.研究结果还显示在高温区域（T〉TDebye/3）内不同温度作用下,所有sI水合物具有相似的导热规律.压力对导热系数有一定影响,尤其是在较高压力条件下,压力越高,导热系数越大.而在不同温度和不同压力作用过程中,密度的改变对导热系数的增大或减小几乎没有影响.     

Structures of hydrocarbon hydrates during formation with and without inhibitors

The formation of hydrates from a methane-ethane-propane mixture is more complex than with single gases. Using nuclear magnetic resonance (NMR) and high-pressure powder X-ray diffraction (PXRD), we have investigated the structural properties of natural gas hydrates crystallized in the presence of kinetic hydrate inhibitors (KHIs), two commercial inhibitors and two biological ice inhibitors, or antifreeze proteins (AFPs). NMR analyses indicated that hydrate cage occupancy was at near saturation for controls and most inhibitor types. Some exceptions were found in systems containing a new commercial KHI (HIW85281) and a recombinant plant AFP, suggesting that these two inhibitors could impact the kinetics of cavity formation. NMR analysis confirmed that the hydrate composition varies during crystal growth by kinetic effects. Strikingly, the coexistence of both structures I (sI) and II (sII) were observed in NMR spectra and PXRD profiles. It is suggested that sI phases may form more readily from liquid water. Real time PXRD monitoring showed that sI hydrates were less stable than sII crystals, and there was a conversion to the stable phase over time. Both commercial KHIs and AFPs had an impact on hydrate metastability, but transient sI PXRD intensity profiles indicated significantly different modes of interaction with the various inhibitors and the natural gas hydrate system.     

Molecular dynamics simulations of binary structure H hydrogen and methyl-tert-butylether clathrate hydrates

Binary structure H (sH) hydrogen and methyl-tert-butylether (MTBE) clathrate hydrates are studied with molecular dynamics simulations. Simulations on a 3 x 3 x 3 sH unit cell with up to 4.7 mass % hydrogen gas are run at pressures of 100 bars and 2 kbars at 100 and 273 K. For the small and medium cages of the sH unit cell, H2 guest molecule occupancies of 0, 1 (single occupancy), and 2 (double occupancy) are considered with the MTBE molecule occupying all of the large cages. An increase of the small and medium cage occupancies from 1 to 2 leads to a jump in the unit cell volume and configurational energy. Calculations are also set up with 13, 23, and 89 of the MTBE molecules in the large cages replaced by sets of three to six H2 molecules, and the effects on the configurational energy and volume of the simulation cell are determined. As MTBE molecules are replaced with sets of H2 guests in the large cages, the configurational energy of the unit cell increases. At the lower temperature, the energy and volume of the clathrate are not sensitive to the number of hydrogen guests in the large cages; however, at higher temperatures the repulsions among the H2 guest molecules in the large cages cause an increase in the system energy and volume.     

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.     

Molecular dynamics simulation for studying the stability of structure H clathrate-hydrates of argon and large guest molecules

The aim of present paper is to study the stability of (argon + large guest molecules) structure H clathrate-hydrates by using molecular dynamics simulations and with employing the COMPASS force field to consider the molecular interactions. The simulations are performed by embedding the structure H clathrate-hydrates in a simulation cell under isobaric-isothermal (NPT) conditions. The obtained equilibrium lattice parameters are compared with the experimental data, where a good consistency is observed. The results show that the size and dipole moment of the guest molecules enclosed in the hydrate cages play the main role in the interactions between the guest molecules and the water molecules, which constitute the surrounding walls of the hydrate cage and these interactions would stabilize the hydrate structure. The characteristics of the clathrate-hydrate structure are analyzed by evaluating the radial distribution function, where the agreement between the results obtained in this work and other similar theoretical and experimental investigations validates the simulation procedure and related interpretations.     

Microscopic observation and in-situ Raman scattering studies on high-pressure phase transformations of Kr hydrate

Direct observations through a microscope and in-situ Raman scattering measurements of synthesized single-crystalline Kr hydrate have been performed at pressures up to 5.2 GPa and 296 K. We have observed that the initial cubic structure II (sII) of Kr hydrate successively transforms to a cubic structure I (sI), a hexagonal structure, and an orthorhombic structure (sO) called "filled ice" at 0.45, 0.75, and 1.8 GPa, respectively. The sO phase exists at least up to 5.2 GPa. In addition to these transformations, we have also found the new phase behavior at 1.0 GPa, which is most likely caused by the change of cage occupancy of host water cages by guest Kr atoms without structural change. Raman scattering measurements for observed phases have shown that the lattice vibrational peak at around 130 cm(-1) disappears in the pressure region of sI, which enables us to distinguish the sI phase from sII and sH phases.     

Molecular-dynamics simulations of binary structure II hydrogen and tetrahydrofurane clathrates

The binary structure II hydrogen and tetrahydrofurane (THF) clathrates are studied with molecular-dynamics simulations. Simulations are done at pressures of 120 and 1.013 bars for temperatures ranging from 100 to 273 K. For the small cages of the structure II unit cell, H2 guest molecule occupancies of 0, 16 (single occupancy), and 32 (double occupancy) are considered. THF occupancies of 0-8 in the large cages are studied. For cases in which THF does not occupy all large cages in a unit cell, the remaining large cages can be occupied with sets of four H2 guest molecules. The unit-cell volumes and configurational energies are compared in the different occupancy cases. Increasing the small cage occupancy leads to an increase in the unit-cell volume and thermal-expansion coefficient. Among simulations with the same small cage occupancy, those with the large cages containing 4H2 guests have the largest volumes. The THF guest molecules have a stabilizing effect on the clathrate and the configurational energy of the unit cell decreases linearly as the THF content increases. For binary THF + H2 clathrates, the substitution of the THF molecules in the large cages with sets of 4H2 molecules increases the configurational energy. For the binary clathrates, various combinations of THF and H2 occupancies have similar configurational energies.     

Structure transition and tuning pattern in the double (tetramethylammonium hydroxide + gaseous guests) clathrate hydrates

In this study, we present an extraordinary structural transition accompanying the occurrence of more than two coexisting clathrate hydrate phases in the double (CH4 + tetramethylammonium hydroxide (Me(4)NOH)) and (H2 + Me(4)NOH) ionic clathrate hydrates using solid-state NMR spectroscopy (high-powered decoupling and CP/MAS) and powder X-ray diffraction. It was confirmed that structure-I (sI) and structure-II (sII) hydrates coexist as the water concentration increases. In the Me(4)NOH-depleted region, the unique tuning phenomenon was first observed at a chemical shift of -8.4 ppm where relatively small gaseous CH4 molecules partly occupy the sII large cages (sII-L), pulling out large cationic Me(4)N+ that is considered to be strongly bound with the surrounding host lattices. Moreover, we note that, while pure Me(4)NOH.16H(2)O clathrate hydrates melted at 249 K under atmospheric pressure conditions, the double (CH4 + Me(4)NOH) clathrate hydrate maintained a solid state up to approximately 283 K under 120 bar of CH4 with a conductivity of 0.065 S cm(-1), suggesting its potential use as a solid electrolyte. The present results indicate that ionic contributions must be taken into account for ionic clathrate hydrate systems because of their distinctive guest dynamic behavior and structural patterns. In particular, microscopic analyses of ionic clathrate hydrates for identifying physicochemical characteristics are expected to provide new insights into inclusion chemistry.