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

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

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.     

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.     

Spectroscopic signatures of halogens in clathrate hydrate cages. 1. Bromine

We report the first UV-vis spectroscopic study of bromine molecules confined in clathrate hydrate cages. Bromine in its natural hydrate occupies 51262 and 51263 lattice cavities. Bromine also can be encapsulated into the larger 51264 cages of a type II hydrate formed mainly from tetrahydrofuran or dichloromethane and water. The visible spectra of the enclathrated halogen molecule retain the spectral envelope of the gas-phase spectra while shifting to the blue. In contrast, spectra of bromine in liquid water or amorphous ice are broadened and significantly more blue-shifted. The absorption bands shift by about 360 cm-1 for bromine in large 51264 cages of type II clathrate, by about 900 cm-1 for bromine in a combination of 51262 and 51263 cages of pure bromine hydrate, and by more than 1700 cm-1 for bromine in liquid water or amorphous ice. The dramatic shift and broadening in water and ice is due to the strong interaction of the water lone-pair orbitals with the halogen sigma* orbital. In the clathrate hydrates, the oxygen lone-pair orbitals are all involved in the hydrogen-bonded water lattice and are thus unavailable to interact with the halogen guest molecule. The blue shifts observed in the clathrate hydrate cages are related to the spatial constraints on the halogen excited states by the cage walls.     

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

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

Molecular Dynamics 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.     

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.     

Participation of xenon guest in hydrogen bond network of β-hydroquinone crystal

In the β-hydroquinone (β-HQ)-Xe crystal, the Xe guest is placed between two hexagonal rings of coupled [···O-H···O-](6) H-bonds. This clathrate is treated as the model for monitoring the H-bonding system with the Xe participation. Three kinds of isotope effects due to the H/D substitution in the [···O-H···O-](6) bonds are considered: (i) structural changes in the clathrate (X-ray diffraction), (ii) variations of (129)Xe NMR signal of the guest (CP MAS), and (iii) variations of selected vibrations of the host (IR). This study predicts subtle inclination of every other hydroxyl group of the [···O-H···O-](6) rings into the Xe atom and formation of six Xe···H-O pairs in every cage, the frequency shift of the γOH mode due to these contacts, -ΔγOH(Xe···H) > 74 cm(-1), as well as the enthalpy formation, -ΔH(Xe···H) > 6-8 kJ mol(-1). Our IR results reveal a tendency of the Xe atom to form the H-bond-like network inside its cage and much weaker Xe···D-O interactions in the H/D substituted crystal. The (129)Xe NMR results do not reflect this kind of interactions due to averaging of the (129)Xe shielding phenomena, probably. We also predict elongation of the O···O distances due to the β-HQ-Xe crystal heating and the Xe escape.     

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.     

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.     

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.     

A water-soluble Xe@cryptophane-111 complex exhibits very high thermodynamic stability and a peculiar (129)Xe NMR chemical shift

The known xenon-binding (±)-cryptophane-111 (1) has been functionalized with six [(η(5)-C(5)Me(5))Ru(II)](+) ([Cp*Ru](+)) moieties to give, in 89% yield, the first water-soluble cryptophane-111 derivative, namely [(Cp*Ru)(6)1]Cl(6) ([2]Cl(6)). [2]Cl(6) exhibits a very high affinity for xenon in water, with a binding constant of 2.9(2) × 10(4) M(-1) as measured by hyperpolarized (129)Xe NMR spectroscopy. The (129)Xe NMR chemical shift of the aqueous Xe@[2](6+) species (308 ppm) resonates over 275 ppm downfield of the parent Xe@1 species in (CDCl(2))(2) and greatly broadens the practical (129)Xe NMR chemical shift range made available by xenon-binding molecular hosts. Single crystal structures of [2][CF(3)SO(3)](6)·xsolvent and 0.75H(2)O@1·2CHCl(3) reveal the ability of the cryptophane-111 core to adapt its conformation to guests.     

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.     

Communication: Single crystal x-ray diffraction observation of hydrogen bonding between 1-propanol and water in a structure II clathrate hydrate

Single crystal x-ray crystallography is used to detect guest-host hydrogen bonding in structure II (sII) binary clathrate hydrate of 1-propanol and methane. X-ray structural analysis shows that the 1-propanol oxygen atom is at a distance of 2.749 and 2.788 ? from the closest clathrate hydrate water oxygen atoms from a hexagonal face of the large sII cage. The 1-propanol hydroxyl hydrogen atom is disordered and at distances of 1.956 and 2.035 ? from the closest cage water oxygen atoms. These distances are compatible with guest-water hydrogen bonding. The C-C-C-O torsional angle in 1-propanol in the cage is 91.47° which corresponds to a staggered conformation for the guest. Molecular dynamics studies of this system demonstrated guest-water hydrogen bonding in this hydrate. The molecular dynamics simulations predict most probable distances for the 1-propanol-water oxygen atoms to be 2.725 ?, and the average C-C-C-O torsional angle to be ~59° consistent with a gauche conformation. The individual cage distortions resulting from guest-host hydrogen bonding from the simulations are rather large, but due to the random nature of the hydrogen bonding of the guest with the 24 water molecules making up the hexagonal faces of the large sII cages, these distortions are not observed in the x-ray structure.     

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.     

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.     

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.     

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.     

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

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