Dynamical and energetic properties of hydrogen and hydrogen-tetrahydrofuran clathrate hydrates

Classical equilibrium molecular dynamics (MD) simulations have been performed to investigate the dynamical and energetic properties in hydrogen and mixed hydrogen-tetrahydrofuran sII hydrates at 30 and 200 K and 0.05 kbar, and also at intermediate temperatures, using SPC/E and TIP4P-2005 water models. The potential model is found to have a large impact on overall density, with the TIP4P-2005 systems being on average 1% more dense than their SPC/E counterparts, due to the greater guest-host interaction energy. For the lightly-filled mixed H(2)-THF system, in which there is single H(2) occupation of the small cage (1s1l), we find that the largest contribution to the interaction energy of both types of guest is the van der Waals component with the surrounding water molecules in the constituent cavities. For the more densely-filled mixed H(2)-THF system, in which there is double H(2) occupation in the small cage (2s1l), we find that there is no dominant component (i.e., van der Waals or Coulombic) in the H(2) interaction energy with the rest of the system, but for the THF molecules, the dominant contribution is again the van der Waals interaction with the surrounding cage-water molecules; again, the Coulombic component increases in importance with increasing temperature. The lightly-filled pure H(2) hydrate (1s4l) system exhibits a similar pattern vis-à-vis the H(2) interaction energy as for the lightly-filled mixed H(2)-THF system, and for the more densely-filled pure H(2) system (2s4l), there is no dominant component of interaction energy, due to the multiple occupancy of the cavities. By consideration of Kubic harmonics, there is some evidence of preferential alignment of the THF molecules, particularly at 200 K; this was found to arise at higher temperatures due to transient hydrogen bonding of the oxygen atom in THF molecules with the surrounding cage-water molecules.     

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.     

Phase relations and binary clathrate hydrate formation in the system H2-THF-H2O

Experimentally determined equilibrium phase relations are reported for the system H2-THF-H2O as a function of aqueous tetrahydrofuran (THF) concentration from 260 to 290 K at pressures up to 45 MPa. Data are consistent with the formation of cubic structure-II (CS-II) binary H2-THF clathrate hydrates with a stoichiometric THF-to-water ratio of 1:17, which can incorporate modest volumes of molecular hydrogen at elevated pressures. Direct compositional analyses of the clathrate phase, at both low (0.20 mol %) and stoichiometric (5.56 mol %) initial THF aqueous concentrations, are consistent with observed phase behavior, suggesting full occupancy of large hexakaidecahedral (51264) clathrate cavities by THF, coupled with largely complete (80-90%) filling of small dodecahedral (512) cages by single H2 molecules at pressures of >30 MPa, giving a clathrate formula of (H2) < or =2.THF.17H2O. Results should help to resolve the current controversy over binary H2-THF hydrate hydrogen contents; data confirm recent reports that suggest a maximum of approximately 1 mass % H2, this contradicting values of up to 4 mass % previously claimed for comparable conditions.     

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.     

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.     

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.     

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.     

Experimental and theoretical analysis of the rotational Raman spectrum of hydrogen molecules in clathrate hydrates

The Raman spectra of H(2) and HD molecules in simple hydrogen and binary hydrogen-tetrahydrofuran clathrate hydrates have been measured at temperatures as low as 20 K. The rotational bands of trapped molecules in simple and binary hydrates have been analyzed, and the contributions originating from hydrogen molecules in the large cages have been separated from those in the small cages. A theoretical model, consisting in rigid cages enclosing interacting hydrogen molecules, has been exploited to calculate, on the basis of quantum mechanics, the Raman intensity of the rotational transitions for up to two interacting molecules in one cage. A comparison with experiment leads to a clear interpretation of sidebands appearing in the Raman rotational lines. The quantitative agreement between theory and experiment obtained in some cases clarifies the importance of the choice of the interaction potential, and of the proton disorder in the clathrate crystal.     

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.     

Pulsed EPR characterization of encapsulated atomic hydrogen in octasilsesquioxane cages

Hydrogen atoms encapsulated in molecular cages are potential candidates for quantum computing applications. They provide the simplest two-spin system where the 1s electron spin, S = 1/2, is hyperfine-coupled to the proton nuclear spin, I = 1/2, with a large isotropic hyperfine coupling (A = 1420.40575 MHz for a free atom). While hydrogen atoms can be trapped in many matrices at cryogenic temperatures, it has been found that they are exceptionally stable in octasilsesquioxane cages even at room temperature [Sasamori et al., Science, 1994, 256, 1691]. Here we present a detailed spin-lattice and spin-spin relaxation study of atomic hydrogen encapsulated in Si(8)O(12)(OSiMe(2)H)(8) using X-band pulsed EPR spectroscopy. The spin-lattice relaxation times T(1) range between 1.2 s at 20 K and 41.8 μs at room temperature. The temperature dependence of the relaxation rate shows that for T < 60 K the spin-lattice relaxation is best described by a Raman process with a Debye temperature of θ(D) = 135 K, whereas for T > 100 K a thermally activated process with activation energy E(a) = 753 K (523 cm(-1)) prevails. The phase memory time T(M) = 13.9 μs remains practically constant between 200 and 300 K and is determined by nuclear spin diffusion. At lower temperatures T(M) decreases by an order of magnitude and exhibits two minima at T = 140 K and T = 60 K. The temperature dependence of T(M) between 20 and 200 K is attributed to dynamic processes that average inequivalent hyperfine couplings, e.g. rotation of the methyl groups of the cage organic substituents. The hyperfine couplings of the encapsulated proton and the cage (29)Si nuclei are obtained through numerical simulations of field-swept FID-detected EPR spectra and HYSCORE experiments, respectively. The results are discussed in terms of existing phenomenological models based on the spherical harmonic oscillator and compared to those of endohedral fullerenes.     

Mechanisms for thermal conduction in hydrogen hydrate

Extensive equilibrium molecular dynamics simulations have been performed to investigate thermal conduction mechanisms via the Green-Kubo approach for (type II) hydrogen hydrate, at 0.05 kbar and between 30 and 250 K, for both lightly filled H(2) hydrates (1s4l) and for more densely filled H(2) systems (2s4l), in which four H(2) molecules are present in the large cavities, with respective single- and double-occupation of the small cages. The TIP4P water model was used in conjunction with a fully atomistic hydrogen potential along with long-range Ewald electrostatics. It was found that substantially less damping in guest-host energy transfer is present in hydrogen hydrate as is observed in common type I clathrates (e.g., methane hydrate), but more akin in to previous results for type II and H methane hydrate polymorphs. This gives rise to larger thermal conductivities relative to common type I hydrates, and also larger than type II and H methane hydrate polymorphs, and a more crystal-like temperature dependence of the thermal conductivity.     

Rotation and diffusion of H2 in hydrogen-ice clathrate by 1H NMR

A recently reported hydrogen-ice clathrate carries up to four H(2) in each large cage and one H(2) in each small cage. We report pulsed proton NMR line shape measurements on H(2)-D(2)O clathrate formed at 1500 bar and 250 K. The behavior of the two-pulse spin-echo amplitude with respect to the nutation angle of the refocusing pulse shows that intramolecular dipolar broadening, modulated by H(2) molecular reorientations, dominates the line width of the ortho-H(2). Dipolar interaction between H(2) guests and host D atoms explains the echo variation with the relative phases of the pulses. From 12 to 120 K, the line width varies as 1/T, demonstrating that the three sublevels of J = 1 are split by a constant energy, epsilon. The splitting arises from distortion in the otherwise high-symmetry cages from frozen-out D(2)O orientational disorder. Above 120 K, further line-narrowing signals the onset of H(2) diffusion from cage to cage. At the lowest temperature, 1.9 K, the spectrum has Pake powder doublet-like features; the doublet is not fully developed, indicating a broad distribution of order parameters and energies epsilon.     

Simulations of structure II H2 and D2 clathrates: potentials incorporating quantum corrections

Molecular dynamics simulations are used to study the stability of structure II H(2) and D(2) clathrates with different large and small guest occupancies at 160 and 250 K and 2.0 kbars. Simulations are performed with the recently proposed anisotropic site-site potentials of Wang for H2 and D2 [J. Quant. Spectrosc. Radiat. Transf. 76, 23 (2003)] which are parameterized to account for quantum corrections of order variant Planck's over 2pi(2) in the second virial coefficient. Occupancies of 0-2 in the small cages and 2-5 in the large cages are considered. Thermodynamic integration is used to determine the most stable guest occupancy at each temperature. Since lattice free energy and configurational energy differences are small for a number of different combinations of cage occupancies, one must expect that in bulk samples various combinations will indeed be observed. Special attention is given to the differences between H(2) and D(2) guests and implications on the hydrogen storage capacity of the clathrates are discussed.     

Intermolecular hydrogen transfer between guest species in small and large cages of methane + propane mixed gas hydrates

To investigate the molecular interaction between guest species inside of the small and large cages of methane + propane mixed gas hydrates, thermal stabilities of the methyl radical (possibly induced in small cages) and the normal propyl and isopropyl radicals (induced in large cages) were investigated by means of electron spin resonance measurements. The increase of the total amount of the normal propyl and isopropyl radicals reveals that the methyl radical in the small cage withdraws one hydrogen atom from the propane molecule enclathrated in the adjacent large cage of the structure-II hydrate. A guest species in a hydrate cage has the ability to interact closely with the other one in the adjacent cages. The clathrate hydrate may be utilized as a possible nanoscale reaction field.     

Structural Stability of the CO2@sI Hydrate: a Bottom-Up Quantum Chemistry Approach on the Guest-Cage and Inter-Cage Interactions

Through reliable first-principles computations, we have demonstrated the impact of CO₂ molecules enclathration on the stability of sI clathrate hydrates. Given the delicate balance between the interaction energy components (van der Waals, hydrogen bonds) present on such systems, we follow a systematic bottom-up approach starting from the individual 5¹² and 5¹²6² sI cages, up to all existing combinations of two-adjacent sI crystal cages to evaluate how such clathrate-like models perform on the evaluation of the guest-host and first-neighbors inter-cage effects, respectively. Interaction and binding energies of the CO₂ occupation of the sI cages were computed using DF-MP2 and different DFT/DFT−D electronic structure methodologies. The performance of selected DFT functionals, together with various semi-classical dispersion corrections schemes, were validated by comparison with reference ab initio DF-MP2 data, as well as experimental data from x-ray and neutron diffraction studies available. Our investigation confirms that the inclusion of the CO₂ in the cage/s is an energetically favorable process, with the CO₂ molecule preferring to occupy the large 5¹²6² sI cages compared to the 5¹² ones. Further, the present results conclude on the rigidity of the water cages arrangements, showing the importance of the inter-cage couplings in the cluster models under study. In particular, the guest-cage interaction is the key factor for the preferential orientation of the captured CO₂ molecules in the sI cages, while the inter-cage interactions seems to cause minor distortions with the CO₂ guest neighbors interactions do not extending beyond the large 5¹²6² sI cages. Such findings on these clathrate-like model systems are in accord with experimental observations, drawing a direct relevance to the structural stability of CO₂@sI clathrates.     

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.     

Electronic properties of mononuclear, dinuclear, and polynuclear cobaltacarboranes: electrochemical and spectroelectrochemical studies

Electronic interactions and metal-metal communication in a wide range of cobaltacarborane-hydrocarbon complexes containing one to six metal centers, and exhibiting a variety of modes of inter-cage connectivity and molecular architectures, have been investigated via cyclic voltammetry, controlled potential coulometry, and UV-visible spectroelectrochemistry. The properties of mixed-valent Co(III)/Co(IV) and Co(II)/Co(III) species that are generated on oxidation or reduction of dinuclear and polynuclear Co(III) complexes were examined and classified as Robin-Day Class I (localized), Class II (partially delocalized), or Class III (fully delocalized) systems. The extent of metal-metal communication between metallacarborane cage units is strongly influenced by the type of intercage connection (e.g., cage B-B or Cp-Cp); the vertexes involved (equatorial vs apical); the nature of the linking unit, if any; and the presence of substituents on the carborane cages. In multi-tripledecker complexes where three CpCo(C(2)B(3)H(4))CoCp units are linked through a central triethynyl benzene connector, the data suggest that Co-Co electronic communication is extensive (Class III) within individual sandwich units while intersandwich delocalization is weak or absent. An extended Hückel study of CpCoC(2)B(4)H(6) double-decker and CpCo(C(2)B(3)H(5))CoCp triple-decker sandwich model complexes shows significant differences in the orbital contributions involved in the HOMO and LUMO of the former vs the latter type. The calculations afford additional insight into the electronic structures and properties of these systems as elucidated by the experimental studies.     

Competitive cage occupancy of hydrogen and argon in structure-II hydrates

The cage occupancy of hydrogen in the single-crystals of simple hydrogen hydrates and hydrogen + argon mixed-gas hydrates was investigated by means of in situ Raman spectroscopy under the three-phase (hydrate + water + fluid) equilibrium condition. In the equilibrium pressure region higher than approximately 25 MPa, four hydrogen cluster and argon competitively occupied the large cages of structure-II hydrogen + argon mixed-gas hydrates. In addition, Raman spectroscopic analysis at liquid nitrogen temperature (77 K) supports that the clusters of two, three, or four hydrogen molecules occupy large cages.     

H2 molecules encapsulated in extended Ben cluster cages: toward light-metal nanofoams for hydrogen storage

Efficient storage of hydrogen is a bottleneck problem for hydrogen-based energy solutions. We demonstrate the feasibility of trapping a pair of hydrogen molecules in beryllium cluster cages. The systems are constructed by merging two smaller units with single molecules trapped, which are known to be stable in isolation. The resulting (H(2))(2)@Be(n) species can have hydrogen cores and beryllium shells of different shapes, and we report the calculated energy barriers for hydrogen exit from the cage. The relative stabilities are related to the molecular structure and charge distributions, and some initially counterintuitive features are explained. Aspects of the release of hydrogen from such structures, and of possible scaling up to larger extended systems of fused cages, are discussed in terms of hydrogen storage. The predicted capacity could potentially be sufficient for practical usage.