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.     

How much carbon dioxide can be stored in the structure H clathrate hydrates?: a molecular dynamics study

The stability of structure H (sH) carbon dioxide clathrate hydrates at three temperature-pressure conditions are determined by molecular dynamics simulations on a 3x3x3 sH unit cell replica. Simulations are performed at 100 K at ambient pressure, 273 K at 100 bars and also 300 K and 5.0 kbars. The small and medium cages of the sH unit cell are occupied by a single carbon dioxide guest and large cage guest occupancies of 1-5 are considered. Radial distribution functions are given for guests in the large cages and unit cell volumes and configurational energies are studied as a function of large cage CO(2) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies at three temperature/pressure conditions stated above. At the low temperature, large cage occupancy of 5 is the most stable while at the higher temperature, the occupancy of 3 is the most favored. Calculations are also performed to show that the CO(2) sH clathrate is more stable than the methane clathrate analog. Implications on CO(2) sequestration by clathrate formation are discussed.     

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.     

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.     

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.     

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.     

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.     

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.     

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 study of structure II hydrogen clathrates

Molecular-dynamics simulations are used to study the stability of structure II hydrogen clathrates with different H2 guest occupancies. Simulations are done at pressures of 2.5 kbars and 1.013 bars and for temperatures ranging from 100 to 250 K. For a structure II unit cell with 136 water molecules, H2 guest molecule occupancies of 0-64 are studied with uniform occupancies among each type of cage. The simulations show that at 100 K and 2.5 kbars, the most stable configurations have single occupancy in the small cages and quadruple occupancy in the large cages. The optimum occupancy for the large cages decreases as the temperature is raised. Double occupancy in the small cages increases the energy of the structures and causes tetragonal distortion in the unit cell. The spatial distribution of the hydrogen guest molecules in the cages is determined by studying the guest-water and guest-guest radial distribution functions at various temperatures.     

Neutron diffraction study of the type I clathrate Ba8Al(x)Si(46-x): site occupancies, cage volumes, and the interaction between the guest and the host framework

Samples with the type I clathrate structure and composition Ba(8)Al(x)Si(46-x), where x = 8, 10, 12, 14, and 15, were examined by neutron powder diffraction at 35 K. The clathrate type I structure contains Ba cations as guests in a framework derived from tetrahedrally coordinated Al/Si atoms. The framework is made up of five- and six-membered rings that form dodecahedral and tetrakaidecahedral cages. The change in distances between tetrahedral sites across the series is used to develop a model for the mixed Al/Si occupancy observed in the framework. The calculated volumes of the cages that contain the Ba atoms display a linear increase with increasing Al composition. In the smaller dodecahedral cages, the Ba atomic displacement parameter is symmetry constrained to be isotropic for all compositions. In the larger tetrakaidecahedral cages, the anisotropic atomic displacement of the Ba atom depends upon the composition: the displacement is perpendicular (x = 8) and parallel (x = 15) to the six-membered ring. This difference in direction of the displacement parameter is attributed to interaction with the Al in the framework and not to the size of the cage volume as x increases from 8 to 15. The influence of the site occupation of Al in the framework on displacement of the cation at the 6d site is demonstrated.     

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.     

NMR shielding constants for hydrogen guest molecules in structure II clathrates

Proton NMR shielding constants and chemical shifts for hydrogen guests in small and large cages of structure II clathrates are calculated using density-functional theory and the gauge-invariant atomic-orbital method. Shielding constants are calculated at the B3LYP level with the 6-311++G(d,p) basis set. The calculated chemical shifts are corrected with a linear regression to reproduce the experimental chemical shifts of a set of standard molecules. The calculated chemical shifts of single hydrogen molecules in the small and large structure II cages are 4.94 and 4.84 ppm, respectively, which show that within the error range of the method the H2 guest molecules in the small and large cages cannot be distinguished. Chemical shifts are also calculated for double occupancy of the hydrogen guests in small cages, and double, triple, and quadruple occupancy in large cages. Multiple occupancy changes the chemical shift of the hydrogen guests by approximately 0.2 ppm. The relative effects of other guest molecules and the cage on the chemical shift are studied for the cages with multiple occupancies.     

Theoretical Study of Clathrate Hydrates with Multiple Occupation

In this work, the electronic, structural, dynamic andthermodynamic properties of structure II, H and tetragonalAr clathrate hydrates have been calculated and the effectof multiple occupancy on their stability has been examinedusing first-principles and lattice dynamics calculations.The dynamic properties of these clathrates have beeninvestigated depending on the number of guest moleculesin a clathrate cage. It has been found that selectedhydrate structures are dynamically stable. The calculatedcell parameters are in agreement with experimental data.We also report the results of a systematic investigationof cage-like water structures using first-principles calculations. Ithas been observed that Ar clusters can be stabilized indifferent water cages and the stability is strongly dependenton the number of argon atoms inside the cages.     

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.     

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.     

Phase equilibrium measurements and the tuning behavior of new sII clathrate hydrates

We suggest two types of new amine-type sII formers: pyrrolidine and piperidine. These guest compounds fail to form clathrate hydrate structures with host water, but instead have to combine with light gaseous guest molecules (methane) for enclathration. First, two binary clathrate hydrates of (pyrrolidine + methane) and (piperidine + methane) were synthesized at various amine concentrations. ¹³C NMR and Raman analysis were done to identify the clathrate hydrate structure and guest distribution over sII-S and sII-L cages. XRD was also used to find the exact structure and corresponding cell parameters. At a dilute pyrrolidine concentration of less than 5.56 mol%, the tuning phenomenon is observed such that methane molecules surprisingly occupy sII-L cages. At the critical guest concentration of about 0.1 mol%, the cage occupancy ratio reaches the maximum of approximately 0.5. At very dilute guest concentration below 0.1 mol%, the methane molecules fail to occupy large cages on account of their rarefied distribution in the network. Direct-release experiments were performed to determine the actual guest compositions in the clathrate hydrate phases. Finally, we measured the clathrate hydrate phase equilibria of (pyrrolidine + methane) and (piperidine + methane).     

The thermodynamic stability of clathrate hydrate. Encaging non-spherical propane molecules

The thermodynamic stability of a clathrate hydrate encaging non-spherical molecules has been investigated by examining the free energy of cage occupancy. In the present study, a generalized van der Waals and Platteeuw theory is extended to treat the rotational motion of guest molecules in clathrate hydrate cages. The vibrational free energy of both guest and host molecules is divided into harmonic and anharmonic contributions. The anharmonic free energy associated with the non-spherical nature of the guest molecules is evaluated as a perturbation from the spherical guest. Predicted thermodynamic properties are compared with measured values. It is shown that this anharmonic contribution is important in the free energy of the hindered rotation of the guests.     

Topotactic transformations of sodalite cages: synthesis and NMR study of mixed salt-free and salt-bearing sodalites

A series of mixed sodalite samples, Na(8)[Al(6)Si(6)O(24)]Br(x).(H(3)O(2))(2-x), with the unit cell stoichiometries varying in the 0 < x <2 region, was made by hydrothermal synthesis and subsequently transformed into Na(6+x)[Al(6)Si(6)O(24)]Br(x).(4H(2)O)(2-x) and Na(6+x)[Al(6)Si(6)O(24)]Br(x).circle(2-x) sodalites. Here, circle refers to an empty sodalite cage. The three series, referred hereafter to as the Br/basic, Br/hydro, and Br/dry series, were characterized by powder diffraction X-ray and by (23)Na, (27)Al, and (81)Br magic angle spinning (MAS) NMR and high-resolution triple quantum (TQ) MAS NMR spectroscopy. We determined that incorporation of Br(-) anions is 130 times more preferred than incorporation of H(3)O(2)(-) anions during the formation of sodalite cages, which permitted precise control of the halide content in the solid. Monotonic trends in chemical shifts were observed as a function of cage occupancy, reflecting continuous changes in structural parameters. A linear correlation between (81)Br chemical shift and lattice constant with a slope of -86 ppm/A was observed for all three series. Likewise, (23)Na chemical shifts for Na(+) cations in salt-bearing sodalite cages correlate linearly with the lattice constant. Both results indicate a universal dependence of the (23)Na and (81)Br chemical shifts on the Na-Br distance. The (27)Al chemical shifts of Br/basic and Br/hydro sodalites obey an established relation between delta(cs) and the average T-O-T bond angle of 0.72 ppm/degrees. Br/dry sodalites show two aluminum resonances, characterized by significantly different chemical shifts and quadrupolar interaction parameters. In that series, local symmetry distortions are evident from strong quadrupolar perturbations in the NMR spectra. P(Q) values for (27)Al vary between 0.8 MHz in Br/basic sodalites and 4.4 MHz in the Br/dry series caused by deviations from the tetrahedral symmetry of the salt-free sodalite cages. For (23)Na, P(Q) values of 0.8, 0.8, 2.0, and 5.7 MHz were found for sodium in bromo, basic, hydro, and dry cages, respectively. In addition, both (23)Na and (81)Br spectra offer some evidence that the Br(-) anions in the Br/dry sodalite are displaced from the center of the expanded sodalite cage. For all three series, the spectral deconvolution of the (23)Na NMR line shapes permits an accurate determination of the mixed sodalite stoichiometry.     

High-pressure synthesis of a new silicon clathrate superconductor, Ba8Si46

A new silicon clathrate compound with a composition of Ba8Si46 was prepared under high-pressure and high-temperature conditions. The compound was isomorphous with Na8Si46 and became a superconductor with a transition temperature of 8.0 K. Barium atoms occupy all of the Si20 and Si24 cages of the clathrate structure. This is the first clathrate superconductor obtained as a bulk phase.