The chemical shifts of Xe in the cages of clathrate hydrate Structures I and II

We report, for the first time, a calculation of the isotropic NMR chemical shift of 129Xe in the cages of clathrate hydrates Structures I and II. We generate a shielding surface for Xe in the clathrate cages by quantum mechanical calculations. Subsequently this shielding surface is employed in canonical Monte Carlo simulations to find the average isotropic Xe shielding values in the various cages. For the two types of cages in clathrate hydrate Structure I, we find the intermolecular shielding values [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-214.0 ppm, and [sigma(Xe@5(12)6(2) cage)-sigma(Xe atom)]=-146.9 ppm, in reasonable agreement with the values -242 and -152 ppm, respectively, observed experimentally by Ripmeester and co-workers between 263 and 293 K. For the 5(12) and 5(12)6(4) cages of Structure II we find [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-206.7 ppm, and [sigma(Xe@5(12)6(4) cage)-sigma(Xe atom)]=-104.7 ppm, also in reasonable agreement with the values -225 and -80 ppm, respectively, measured in a Xe-propane type II mixed clathrate hydrate at 77 and 220-240 K by Ripmeester et al.     

Xe chemical shift tensor in silicalite and SSZ-24

We report, for the first time, a theoretical prediction of the (129)Xe nuclear magnetic resonance chemical shift tensor of xenon atom in a single crystal of silicalite at near-zero occupancy and the temperature dependence of the Xe NMR chemical shift tensor for the polycrystalline silicalite at maximum occupancy. The former is a measure of the sensitivity of the Xe tensor components to the local structure of the channels without Xe-Xe contributions. The latter is a measure of the sensitivity of the Xe-Xe tensor components to the Xe-Xe distributions, as determined by the Xe-Xe potential function in competition with the Xe-silicalite potential function. Both theoretical predictions can be compared against Xe NMR experiments: the first against the Xe spectra collected as a function of rotation of the single crystal about the three crystalline axes in a magnetic field, and the second against variable temperature Xe NMR studies (below room temperature) of polycrystalline silicalite at maximum Xe occupancy. With the same parameter set (Xe-O potential and shielding functions), we predict the line shapes of Xe in SSZ-24 zeolite under various conditions of occupancy and temperature.     

129Xe nuclear magnetic resonance in the clathrate hydrate of xenon

The sensitivity of the nuclear magnetic shielding of the Xe atom to the physical environment makes it possible to “see” well-separated. ¹²⁹Xe resonances from Xe atoms in the small and large clathrate deuteriohydrate cages as well as in the gas with which the hydrate is in equilibrium. The relative occupancy of the small cages is suostantially less than predicted oy existing models of guest-host interactions in Xe gas hydrate.     

Diastereomeric Xe chemical shifts in tethered cryptophane cages

Cryptophane cages serve as host molecules to a Xe atom. Functionalization of cryptophane-A has permitted the development of Xe as a biosensor. Synthetic routes used to prepare cryptophanes result in racemic mixtures of the chiral cages. In the preparation of a tethered cryptophane-A cage for biosensor applications, some achiral and chiral substituents such as left-handed amino acids have been used. When the substituent is achiral, the NMR signal of the Xe atom in the functionalized cage in solution is a single isotropic peak, since the Xe shielding tensor components in the R and L cages differ by no more than the signs of the off-diagonal elements. Chiral substituents can split the cage-encapsulated Xe NMR signal into one or more sets of doublets, depending on the number of asymmetric centers in the substituent. We carry out quantum mechanical calculations of Xe nuclear magnetic shielding for the Xe atom at the same strategic position within an L cryptophane-A cage, under the influence of chiral potentials that represent r or l substituents outside the cage. Calculations of the Xe shielding response in the Lr and Ll diastereomeric pairs permit the prediction of the relative order of the Xe chemical shifts in solutions containing the Rl and Ll diastereomers. Where the substituent itself possesses two chiral centers, comparison of the calculated isotropic shielding responses in the Llr, Lrl, Rll, and Lrr systems, respectively, permits the prediction of the Xe spectrum of diastereomeric systems in solutions containing Llr, Rlr, Lll, and Rll systems. Assignment of the peaks observed in the experimental Xe NMR spectra is therefore possible, without having to undertake the difficult synthetic route that produces a single optically pure enantiomer.     

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.     

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.     

ChemInform Abstract: NMR of 129Xe Trapped in Clathrates and Some Other Solids.

129Xe NMR chemical shifts, obtained for Xe trapped in the cages of clathrate hydrates and a clathrasil sample, together with shift data for solid Xe correlate linearly with the mean free radius of the cage.     

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 Structure II Clathrate Hydrates of Xenon and Large Hydrocarbon Guest Molecules

Molecular dynamics (MD) simulations of structure II clathrate hydrates are performed under canonical (NVT) and isobaric–isothermal (NPT) ensembles. The guest molecule as a small help gas is xenon and gases such as cyclopropane, isobutane and propane are used as large hydrocarbon guest molecule (LHGM). The dynamics of structure II clathrate hydrate is considered in two cases: empty small cages and small cages containing xenon. Therefore, the MD results for structure II clathrate hydrates of LHGM and LHGM + Xe are obtained to clarify the effects of guest molecules on host lattice structure. To understand the characteristic configurations of structure II clathrate hydrate the radial distribution functions (RDFs) are calculated for the studied hydrate system. The obtained results indicate the significance of interactions of the guest molecules on stabilizing the hydrate host lattice and these results is consistent with most previous experimental and theoretical investigations.     

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.     

Effect of Guest–Host Hydrogen Bonding on the Structures and Properties of Clathrate Hydrates

To provide improved understanding of guest–host interactions in clathrate hydrates, we present some correlations between guest chemical structures and observations on the corresponding hydrate properties. From these correlations it is clear that directional interactions such as hydrogen bonding between guest and host are likely, although these have been ignored to greater or lesser degrees because there has been no direct structural evidence for such interactions. For the first time, single‐crystal X‐ray crystallography has been used to detect guest–host hydrogen bonding in structure II (sII) and structure H (sH) clathrate hydrates. The clathrates studied are the tert‐butylamine (tBA) sII clathrate with H₂S/Xe help gases and the pinacolone + H₂S binary sH clathrate. X‐ray structural analysis shows that the tBA nitrogen atom lies at a distance of 2.64 Å from the closest clathrate hydrate water oxygen atom, whereas the pinacolone oxygen atom is determined to lie at a distance of 2.96 Å from the closest water oxygen atom. These distances are compatible with guest–water hydrogen bonding. Results of molecular dynamics simulations on these systems are consistent with the X‐ray crystallographic observations. The tBA guest shows long‐lived guest–host hydrogen bonding with the nitrogen atom tethered to a water HO group that rotates towards the cage center to face the guest nitrogen atom. Pinacolone forms thermally activated guest–host hydrogen bonds with the lattice water molecules; these have been studied for temperatures in the range of 100–250 K. Guest–host hydrogen bonding leads to the formation of Bjerrum L‐defects in the clathrate water lattice between two adjacent water molecules, and these are implicated in the stabilities of the hydrate lattices, the water dynamics, and the dielectric properties. The reported stable hydrogen‐bonded guest–host structures also tend to blur the longstanding distinction between true clathrates and semiclathrates.     

Grand canonical Monte Carlo simulations of the 129Xe NMR line shapes of xenon adsorbed in (+/-)-[Co(en)3]Cl3

The 129Xe NMR line shapes of xenon adsorbed in the nanochannels of the (+/-)-[Co(en)3]Cl3 ionic crystal have been calculated by grand canonical Monte Carlo (GCMC) simulations. The results of our GCMC simulations illustrate their utility in predicting 129Xe NMR chemical shifts in systems containing a transition metal. In particular, the nanochannels of (+/-)-[Co(en)3]Cl3 provide a simple, yet interesting, model system that serves as a building block toward understanding xenon chemical shifts in more complex porous materials containing transition metals. Using only the Xe-C and Xe-H potentials and shielding response functions derived from the Xe@CH4 van der Waals complex to model the interior of the channel, the GCMC simulations correctly predict the 129Xe NMR line shapes observed experimentally (Ueda, T.; Eguchi, T.; Nakamura, N.; Wasylishen, R. E. J. Phys. Chem. B 2003, 107, 180-185). At low xenon loading, the simulated 129Xe NMR line shape is axially symmetric with chemical-shift tensor components delta(parallel) = 379 ppm and delta(perpendicular) = 274 ppm. Although the simulated isotropic chemical shift, delta(iso) = 309 ppm, is overestimated, the anisotropy of the chemical-shift tensor is correctly predicted. The simulations provide an explanation for the observed trend in the 129Xe NMR line shapes as a function of the overhead xenon pressure: delta(perpendicular) increased from 274 to 292 ppm, while delta(parallel) changed by only 3 ppm over the entire xenon loading range. The overestimation of the isotropic chemical shifts is explained based upon the results of quantum mechanical 129Xe shielding calculations of xenon interacting with an isolated (+/-)-[Co(en)3]Cl3 molecule. The xenon chemical shift is shown to be reduced by about 12% going from the Xe@[Co(en)3]Cl3 van der Waals complex to the Xe@C2H6 fragment.     

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.     

Multiple H2 occupancy of cages of clathrate hydrate under mild conditions

Experiments were carried out by reacting H(2) gas with N(2) hydrate at a temperature of 243 K and a pressure of 15 MPa. The characterizations of the reaction products indicated that multiple H(2) molecules can be loaded into both large and small cages of structure II clathrate hydrates. The realization of multiple H(2) occupancy of hydrate cages under moderate conditions not only brings new insights into hydrogen clathrates but also refreshes the perspective of clathrate hydrates as hydrogen storage media.     

Cage occupancies in the high pressure structure H methane hydrate: a neutron diffraction study

A neutron diffraction study was performed on the CD(4) : D(2)O structure H clathrate hydrate to refine its CD(4) fractional cage occupancies. Samples of ice VII and hexagonal (sH) methane hydrate were produced in a Paris-Edinburgh press and in situ neutron diffraction data collected. The data were analyzed with the Rietveld method and yielded average cage occupancies of 3.1 CD(4) molecules in the large 20-hedron (5(12)6(8)) cages of the hydrate unit cell. Each of the pentagonal dodecahedron (5(12)) and 12-hedron (4(3)5(6)6(3)) cages in the sH unit cell are occupied with on average 0.89 and 0.90 CD(4) molecules, respectively. This experiment avoided the co-formation of Ice VI and sH hydrate, this mixture is more difficult to analyze due to the proclivity of ice VI to form highly textured crystals, and overlapping Bragg peaks of the two phases. These results provide essential information for the refinement of intermolecular potential parameters for the water-methane hydrophobic interaction in clathrate hydrates and related dense structures.     

Superexchange-like interaction of encaged molecular oxygen in nitrogen-doped water cages of clathrate hydrates

Clathrate hydrates are a highly prospective material in energy and environmental fields, but the inherent nature of inclusion phenomena occurring in the stacked water cages has not been completely resolved yet. Investigating the magnetism of guest molecules is a new experimental approach in clathrate hydrate research to open the possibility of icy magnetic applications as a novel material as well as to understand the unrevealed host-guest interactions in icy inclusion compounds. In this study, we observed an indirect spin coupling between encaged dioxygen molecules via a nonmagnetic water framework through the measurement of guest magnetization. This spin coupling is reminiscent of superexchange coupling between magnetic ions through intervening oxygens in antiferromagnetic oxides, such as MnO and CoO. Theoretical calculations revealed that OH(-) incorporated in the framework induced the mixing of perpendicular π* orbitals of two distant dioxygens and that ammonia doping into the hydrate cage leads to a longer lifetime of that orientation.     

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.     

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.