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.     

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.     

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.     

Controlling nonclassical content of clathrate hydrates through the choice of molecular guests and temperature

Low-temperature, low-pressure studies of clathrate hydrates (CHs) have revealed that small ether and other proton-acceptor guests greatly enhance rates of clathrate hydrate nucleation and growth; rapid formation and transformations are enabled at temperatures as low as 110 K, and cool moist vapors containing small ether molecules convert to mixed-gas CHs on a subsecond time scale. More recently, FTIR spectroscopic studies of the tetrahydrofuran (THF)-HCN double clathrate hydrate revealed a sizable frequency shift accompanied by a four-fold intensification of the C-N stretch-mode absorption of the small cage HCN, behavior that is enhanced by cooling and which correlates precisely with similar significant changes of the ether C-O/C-C stretch modes. These temperature-dependent correlated changes in the infrared spectra have been attributed to equilibrated extensive hydrogen bonding of neighboring large- and small-cage guest molecules with water molecules of the intervening wall. An ether guest functions as a proton acceptor, particularly so when complemented by the action of a proton-donor (HCN)/electron-acceptor (SO(2)) small-cage guest. Because guest molecules of the classic clathrate hydrates do not participate in hydrogen bonds with the host water, this H-bonding of guests has been labeled "nonclassical". The present study has been enriched by comparing observed FTIR spectra with high-level molecular orbital computational results for guests and hydrogen-bonded guest-water dimers. Vibrational frequency shifts, from heterodimerization of ethers and water, correlate well with the corresponding observed classical to nonclassical shifts. The new spectroscopic data reveal that the nonclassical structures can contribute at observable levels to CH infrared spectra for a remarkable range of temperatures and choice of guest molecules. By the choice of guest molecules, it is now possible to select the abundance levels of nonclassical configurations, ranging from ～0 to 100%, for a given temperature. This ability is expected to hasten understanding of the role of guest-induced nonclassical structures in the acceleration or inhibition of the rates of CH formation and transformation.     

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

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

Water cavities of sH clathrate hydrate stabilized by molecular hydrogen

X-ray diffraction and Raman spectroscopic measurements confirm that molecular hydrogen can be contained within the small water cavities of a binary sH clathrate hydrate using large guest molecules that stabilize the large cavity. The potential increase in hydrogen storage could be more than 40% when compared with binary sII hydrates. This work demonstrates the stabilization of hydrogen in a hydrate structure previously unknown for encapsulating molecular hydrogen, indicating the potential for other inclusion compound materials with even greater hydrogen storage capabilities.     

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

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

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.     

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.     

Thermal expansivity of tetrahydrofuran clathrate hydrate with diatomic guest molecules

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

Time-resolved in situ neutron diffraction studies of gas hydrate: transformation of structure II (sII) to structure I (sI).

We report the in situ observation from diffraction data of the conversion of a gas hydrate with the structure II (sII) lattice to one with the structure I (sI) lattice. Initially, the in situ formation, dissociation, and reactivity of argon gas clathrate hydrate was investigated by time-of-flight neutron powder diffraction at temperatures ranging from 230 to 263 K and pressures up to 5000 psi (34.5 MPa). These samples were prepared from deuterated ice crystals and transformed to hydrate by pressurizing the system with argon gas. Complete transformation from D(2)O ice to sII Ar hydrate was observed as the sample temperature was slowly increased through the D(2)O ice melting point. The transformation of sII argon hydrate to sI hydrate was achieved by removing excess Ar gas and exposing the hydrate to liquid CO(2) by pressurizing the Ar hydrate with CO(2). Results suggest the sI hydrate formed from CO(2) exchange in argon sII hydrate is a mixed Ar/CO(2) hydrate. The proposed exchange mechanism is consistent with clathrate hydrate being an equilibrium system in which guest molecules are exchanging between encapsulated molecules in the solid hydrate and free molecules in the surrounding gas or liquid phase.     

Molecular hydrogen occupancy in binary THF-H2 clathrate hydrates by high resolution neutron diffraction

We have determined the time-space average filling of hydrogen molecules in a binary tetrahydrofuran (THF)-d(8) + D(2) sII clathrate hydrate using high resolution neutron diffraction. The filling of hydrogen in the lattice of a THF-d(8) clathrate hydrate occurred upon pressurization. The hydrogen molecules were localized in the small dodecahedral cavities at 20 K, with nuclear density from the hydrogen approximately spherically distributed and centered in the small cavity. With a formation pressure of 70 MPa, molecular hydrogen was found to only singly occupy the sII small cavity. This result helps explain discrepancies about the hydrogen occupancy in the THF binary hydrate system.     

Abnormal Thermal Expansion of Clathrate Hydrates Induced by Asymmetric Guest Molecules

We investigated for the first time the abnormal thermal expansion induced by an asymmetric guest structure using high‐resolution neutron powder diffraction. Three dihydrogen molecules (H₂, D₂, and HD) were tested to explore the guest dynamics and thermal behavior of hydrogen‐doped clathrate hydrates. We confirmed the restricted spatial distribution and doughnut‐like motion of the HD guest in the center of anisotropic sII‐S (sII‐S=small cages of structure II hydrates). However, we failed to observe a mass‐dependent relationship when comparing D₂ with HD. The use of asymmetric guest molecules can significantly contribute to tuning the cage dimension and thus can improve the stable inclusion of small gaseous molecules in confined cages.     

A new method for screening potential sII and sH hydrogen clathrate hydrate promoters with model potentials

A new predictive computational method for classifying clathrate hydrate promoter molecules is presented, based on the interaction energies between potential promoters and the water networks of sII and sH clathrates. The motivation for this work is identifying promoters for storing hydrogen compactly in clathrate hydrates. As a first step towards achieving this goal, we have developed a general method aimed at distinguishing between molecules that form sII clathrate hydrates and molecules that can-together with a weakly interacting help gas-form sH clathrate hydrates. The new computational method calculates differences in estimated formation energies of the sII and the sH clathrate hydrate. Model interaction potentials have been used, including the electrostatic interactions with newly calculated partial charges for all the considered potential promoter molecules. The methodology can discriminate between the clathrate structure types (sII or sH) formed by each potential promoter with good selectivity, i.e., better than achieved with a simple van der Waals diameter criterion.     

Water cavities of sH clathrate hydrate stabilized by molecular hydrogen: phase equilibrium measurements

In this experimental phase equilibrium study, we show for the first time that it is possible to stabilize structure sH of hydrogen clathrate hydrate with the help of some selected promoters. It was established that the formation pressures of these systems are significantly higher than that of structure sII of hydrogen clathrate hydrate when tetrahydrofuran (THF) is used as a promoter. Although no experimental evidence is available yet, it is estimated that the hydrogen storage capacity of structure sH can be as high as 1.4 wt % of H2, which is about 40% higher compared to the hydrogen storage capacity in structure sII.     

On the thermodynamic stability of hydrogen clathrate hydrates

The cage occupancy of hydrogen clathrate hydrate has been examined by grand canonical Monte Carlo (GCMC) simulations for wide ranges of temperature and pressure. The simulations are carried out with a fixed number of water molecules and a fixed chemical potential of the guest species so that hydrogen molecules can be created or annihilated in the clathrate. Two types of the GCMC simulations are performed; in one the volume of the clathrate is fixed and in the other it is allowed to adjust itself under a preset pressure so as to take account of compression by a hydrostatic pressure and expansion due to multiple cage occupancy. It is found that the smaller cage in structure II is practically incapable of accommodating more than a single guest molecule even at pressures as high as 500 MPa, which agrees with the recent experimental investigations. The larger cage is found to encapsulate at most 4 hydrogen molecules, but its occupancy is dependent significantly on the pressure of hydrogen.     

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).     

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.     

A tetragonal polymorph of caesium hydroxide monohydrate,CsOH·H2O,from X‐ray powder data

Tetragonal caesium hydro­xide monohydrate, CsOH·H₂O, a clathrate hydrate, is a polymorph of three known hexagonal or pseudo‐hexagonal modifications. It was obtained as a by‐product in a high‐pressure experiment. Whether it is a high‐pressure polymorph, however, remains to be verified. The Cs atoms are situated in cavities of the form of a bicapped pentagonal prism, within an infinite three‐dimensional hydrogen‐bonded oxy­gen framework that is locally identical to layers found in the hexagonal modifications. The Cs atom and one of the two H atoms are at sites with symmetry, the O atom has mm site symmetry and the second H atom has 2/m symmetry.     

On the thermodynamic stability of clathrate hydrates IV: double occupancy of cages

We have extended the van der Waals and Platteeuw theory to treat multiple occupancy of a single cage of clathrate hydrates, which has not been taken into account in the original theory but has been experimentally confirmed as a real entity. We propose a simple way to calculate the free energy of multiple cage occupancy and apply it to argon clathrate structure II in which a larger cage can be occupied by two argon atoms. The chemical potential of argon is calculated treating it as an imperfect gas, which is crucial to predict accurate pressure dependence of double occupancy expected at high pressure. It is found that double occupancy dominates over single occupancy when the guest pressure in equilibrium with the clathrate hydrate exceeds 270 MPa.