Direct Space Methods for Powder X-ray Diffraction for Guest−Host Materials: Applications to Cage Occupancies and Guest Distributions in Clathrate Hydrates

On the phase behaviors of CH4-CO2 binary clathrate hydrates: Two-phase and three-phase coexistences

We develop a statistical mechanical theory on clathrate hydrates in order to explore the phase behaviors of clathrate hydrates containing two kinds of guest species and apply it to CH4-CO2 binary hydrates. The two boundaries separating water and hydrate and hydrate and guest fluid mixtures are estimated, which are extended to the lower temperature and the higher pressure region far distant from the three-phase coexisting conditions. The chemical potentials of individual guest components can be calculated from free energies of cage occupations, which are available from intermolecular interactions between host water and guest molecules. This allows us to derive all thermodynamic properties pertinent to the phase behaviors in the whole space of thermodynamic variables of temperature, pressure, and guest compositions. It is found that the phase boundaries of CH4-CO2 binary hydrates with water and with fluid mixtures locate between simple CH4 and CO2 hydrates, but the composition ratios of CH4 guests in hydrates are disproportional to those in fluid mixtures. Such differences arise from the affinities of each guest species to the large and small cages of CS-I hydrates and significantly affect occupation of each cage type, which results in a deviation of the guest composition in hydrates from that in fluid on the two-phase equilibrium conditions. The present method provides a basis for the evaluation of the efficiency of the guest CH4 replacement to CO2 at the thermodynamic limit.

Computational investigations of stable multiple-cage-occupancy He clathrate-like hydrostructures

One of the several possibilities offered by the interesting clathrate hydrates is the opportunity to encapsulate several atoms or molecules, in such a way that more efficient storage materials could be explored or new molecules that otherwise do not exist could be created. These types of applications are receiving growing attention from technologists and chemists, given the future positive implications that they entail. In this context, we investigated the multiple cage occupancy of helium clathrate hydrates, to establish stable novel hydrate structures or ones similar to those predicted previously by experimental and theoretical studies. To this purpose, we analyzed the feasibility of including an increased number of He atoms inside the small (D) and large (H) cages of the sII structure through first-principles properly assessed density functional approaches. On the one hand, we have computed energetic and structural properties, in which we examined the guest-host and guest-guest interactions in both individual and two-adjacent clathrate-like sII cages by means of binding and evaporation energies. On the other hand, we have carried out a thermodynamical analysis on the stability of such He-containing hydrostructures in terms of changes in enthalpy, ΔH, Gibbs free energy, ΔG, and entropy, ΔS, during their formation process at various temperature and pressure values. In this way, we have been able to make a comparison with experiments, reaffirming the ability of computational DFT approaches to describe such weak guest-host interactions. In principle, the most stable structure involves the encapsulation of one and four He atoms inside the D and H sII cages, respectively; however, more He atoms could be entrapped under lower temperature and/or higher pressure thermodynamic conditions. We foresee such accurate computational quantum chemistry approaches contributing to the current emerging machine-learning model development.

Direct Measurement of the Four-Phase Equilibrium Coexistence Vapor-Aqueous Solution-Ice-Gas Hydrate in Water-Carbon Dioxide System

Precise data on the non-variant equilibrium of the four phases (vapor-aqueous solution-ice-gas hydrate) in P-T coordinates are highly desired for developing accurate thermodynamic models and can be used as reference points (similar to the triple point of water). Using the two-component hydrate-forming system CO₂-H₂O, we have proposed and validated a new express procedure for determining the temperature and pressure of the lower quadruple point Q₁. The essence of the method is the direct measurement of these parameters after the successive formation of the gas hydrate and ice phases in the initial two-phase gas-water solution system under intense agitation of the fluids. After relaxation, the system occurs in the same equilibrium state (T = 271.60 K, P = 1.044 MPa), regardless of the initial parameters and the order of crystallization of the CO₂ hydrate and ice phases. Considering the combined standard uncertainties (±0.023 K, ±0.021 MPa), the determined P and T values agree with the results of other authors obtained by a more sophisticated indirect method. Validating the developed approach for systems with other hydrate-forming gases is of great interest.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = -7(5).

Thermally Induced Phase Transition of Cubic Structure II Hydrate: Crystal Structures of Tetrahydropyran-CO2 Binary Hydrate

We report a thermally induced phase transition of cubic structure II hydrates of tetrahydropyran (THP) and CO₂ below about 140 K. The phase transition was characterized by powder X-ray diffraction measurements at variable temperatures. A dynamical ordering of the CO₂ guests in small pentagonal dodecahedral 5¹² host water cages, not previously observed in the simple CO₂ hydrate, occurs simultaneously with the symmetry lowering transition from a cubic structure II (space group Fd-3m with cell dimensions a = 17.3202(7) Å at 153 K) to a tetragonal (space group I4₁/amd with cell dimensions a = 17.484(4) Å and c = 12.145(1) Å at 138 K) unit cell. The effect of guest molecules on the phase transition at low temperatures is discussed, which demonstrates that the clathrate hydrate structures and thermodynamic properties can be modified by adjusting the size and chemical structure of larger and smaller guest molecules.

Modeling of Structure H Carbon Dioxide Clathrate Hydrates: Guest-Lattice Energies, Crystal Structure, and Pressure Dependencies

We performed first-principles computations to investigate the complex interplay of molecular interaction energies in determining the lattice structure and stability of CO₂@sH clathrate hydrates. Density functional theory computations using periodic boundary conditions were employed to characterize energetics and the key structural properties of the sH clathrate crystal under pressure, such as equilibrium lattice volume and bulk modulus. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. Structural relaxations of the fully CO₂-filled and empty sH unit cells yield crystal structure and lattice energies, while their compressibility parameters were derived by including the pressure dependencies. The present quantum chemistry computations suggest anisotropy in the compressibility of the sH clathrate hydrates, with the crystal being less compressible along the a-axis direction than along the c-axis one, in distinction from nearly isotropic sI and sII structures. The detailed results presented here give insight into the complex nature of the underlying guest-host interactions, checking earlier assumptions, providing critical tests, and improving estimates. Such entries may eventually lead to better predictions of thermodynamic properties and formation conditions, with a direct impact on emerging hydrate-based technologies.

Effect of Cage Occupancy on Stability and Decomposition of Methane Hydrate

Gas hydrates usually contain a certain number of empty cages that will both affect the hydrate stability and reduce the gas storage capacity. In this work, by MD simulations, we found that the hydrate stability is related to the cage occupancy, the empty cage types, and especially the distribution of empty cages. With the decrease of overall occupancy, the stability of hydrate becomes worse. Under the same overall occupancy, the more concentrated the empty cages are, the more unstable the hydrate is and hence the faster it decomposes. The methane molecules may migrate between distorted cages during the decomposition, resulting in a temporary increase in the stability of hydrate. Hydrates with different empty cage distributions show different decomposition mechanisms: when empty cages are concentrated, the melting rate is fast first due to the rapid decomposition of empty cages, but the remaining filled cages will reduce the melting rate; when empty cages are separated on the contrary, the early melting is slow because of the high local occupancy, and the following melting will be accelerated because of the high melting surface area. It indicates that the empty cage distribution plays a controlling role in hydrate decomposition kinetics at different stages.

Instrumental Methods for Cage Occupancy Estimation of Gas Hydrate

Studies revealed that gas hydrate cages, especially small cages, are incompletely filled with guest gas molecules, primarily associated with pressure and gas composition. The ratio of hydrate cages occupied by guest molecules, defined as cage occupancy, is a critical parameter to estimate the resource amount of a natural gas hydrate reservoir and evaluate the storage capacity of methane or hydrogen hydrate as an energy storage medium and carbon dioxide hydrate as a carbon sequestration matrix. As the result, methods have been developed to investigate the cage occupancy of gas hydrate. In this review, several instrument methods widely applied for gas hydrate analysis are introduced, including Raman, NMR, XRD, neutron diffraction, and the approaches to estimate cage occupancy are summarized.

The impact of the abnormal salinity enrichment in pore water on the thermodynamic stability of marine natural gas hydrates in the Arctic region

In this study, the thermodynamic and structural characteristics of natural gas hydrates (NGHs) retrieved from gas hydrate mounds (ARAON Mound 03 (AM03) and ARAON Mound 06 (AM06)) in the Chukchi Sea in the Arctic region were investigated. The gas compositions, crystalline structure, and cage occupancy of the NGHs at AM03 and AM06 were experimentally measured using gas chromatography (GC), ¹³C nuclear magnetic resonance (NMR), Raman spectroscopy, and powder X-ray diffraction (PXRD). In the NGHs from AM03 and AM06, a significantly large fraction of CH₄ (> 99%) and a very small amount of H₂S were enclathrated in small (5¹²) and large (5¹²6²) cages of sI hydrate. The NGHs from AM03 and AM06 were almost identical in composition, guest distributions, and existing environment to each other. The salinity of the residual pore water in the hydrate-bearing sediment (AM06) was measured to be 50.32‰, which was much higher than that of seawater (34.88‰). This abnormal salinity enrichment in the pore water of the low-permeability sediment might induce the dissociation of NGHs at a lower temperature than expected. The saturation changes in the NGHs that corresponded with an increase in the seawater temperature were also predicted on the basis of the salinity changes in the pore water. The experimental and predicted results of this study would be helpful for understanding the thermodynamic stability of NGHs and potential CH₄-releasing phenomena in the Arctic region.

On the size, shape and energetics of the hydration shell around alkanes

A large number of clathrate-like cages have been proposed as the very first hydration shell of alkanes. The cages include canonical structures commonly found in clathrate hydrates and many others, not previously reported, derived from the carbon fullerene cavities. These structures have a rich and variegated form, which can adapt to the shape and conformation of the solute. They avoid "wasting" hydrogen bonds, while minimizing the volume cage and maximizing the solute-solvent van der Waals interactions. DFT/M06-2X and MP2 ab initio calculations give comparable structural and energetic results although the latter predicts slightly larger cages for a given solute. It is shown that the van der Waals interactions are substantial and the large exoenergetic values found for isobutane and cyclopentane provide an explanation for the surprising high melting points of related hydrates at room pressure. The encaging enthalpy for various hydrocarbons is similar to the enthalpy of solution measured at a temperature just above the melting point of aqueous hydrocarbon solutions, thus indicating that water molecules should not deviate too much from the configuration with O-H bonds tangentially oriented with respect to the solute surface. The computed trend differs from the enthalpy of solution measured at room temperature, thus the very first hydration shell departs, up to a certain degree, from the clathrate-like structures.

Thermodynamic and structural features of chlorodifluoromethane (a sI-sII dual hydrate former) + external guest (N2 or CH4) hydrates and their significance for greenhouse gas separation

In this study, a new sI-sII dual hydrate former [chlorodifluoromethane (CHClF2); an important greenhouse gas with a global warming potential of 1810], which forms sI hydrate by itself and forms sII hydrate in the presence of external help guests such as CH4 and N2, was introduced and closely investigated for its potential significance in gas hydrate-based gas separation. The phase equilibria of CHClF2 hydrate, binary CHClF2 (5%) + N2 (95%) hydrate, and binary CHClF2 (5%) + CH4 (95%) hydrate were measured to examine the formation conditions and thermodynamic stability regions of CHClF2 + external guest hydrates. Nuclear magnetic resonance and in situ Raman spectroscopic results confirmed the formation of sII hydrates for CHClF2 + external guest (N2 or CH4) mixtures. Powder X-ray diffraction patterns clearly demonstrated a structural transition of sI to sII hydrates and a preferential incorporation of CHClF2 molecules in the hydrate phase when external guests (N2 or CH4) were involved in CHClF2 hydrate formation. The measured dissociation enthalpy values of CHClF2 hydrate, binary CHClF2 (5%) + N2 (95%) hydrate, and binary CHClF2 (5%) + CH4 (95%) hydrate using a high-pressure micro-differential scanning calorimeter also indicated preferential CHClF2 enclathration. The experimental results provide new insights into the thermodynamic and structural features of the CHClF2 (sI-sII dual hydrate former) + external guest hydrates for understanding and designing the hydrate-based CHClF2 separation process.

Evaluation of kinetic salt-enrichment behavior and separation performance of HFC-152a hydrate-based desalination using an experimental measurement and a thermodynamic correlation

Hydrate-based desalination (HBD), a type of freezing-based desalination, can concentrate salts of saline water and produce fresh water via hydrate crystal formation. In this study, the thermodynamic stability, crystallographic information, and kinetic growth behavior of HFC-152a hydrate were investigated to estimate the desalination efficiency of HBD. The phase equilibria revealed that the HFC-152a hydrate could be formed at a higher temperature in the presence of NaCl (0 wt%, 3.5 wt%, and 8.0 wt%) than the HFC-134a hydrate at 0.3 MPa. The hydration number of the HFC-152a hydrate (sI) was found to be 7.74 through the Rietveld refinement of the powder X-ray diffraction patterns, and it was also used to determine the dissociation enthalpy of the HFC-152a hydrate. The Hu-Lee-Sum correlation was employed to predict the equilibrium shift and hydrate depression temperature of both HFC-152a and HFC-134a hydrates in the presence of NaCl. Faster hydrate growth kinetics and higher hydrate conversion were observed for the HFC-152a hydrate in saline solutions despite the smaller initial driving force at 0.3 MPa and the subcooling temperature of 3 K. Additionally, to quantify the desalination efficiency of the HFC-152a HBD, the maximum achievable salinity and maximum water yield were examined using the HLS correlation. The salt-enrichment efficiency decreased with an increase in the initial salinity and increased with increasing the subcooling. The overall results indicate that HFC-152a is, potentially, a superior candidate for HBD. The novel approach examined in this study will be useful for assessing the desalination efficiency of the HBD process.

Exploring CO2 @sI Clathrate Hydrates as CO2 Storage Agents by Computational Density Functional Approaches

The formation of specific clathrate hydrates and their transformation at given thermodynamic conditions depends on the interactions between the guest molecule/s and the host water lattice. Understanding their structural stability is essential to control structure-property relations involved in different technological applications. Thus, the energetic aspects relative to CO₂ @sI clathrate hydrate are investigated through the computation of the underlying interactions, dominated by hydrogen bonds and van der Waals forces, from first-principles electronic structure approaches. The stability of the CO₂ @sI clathrate is evaluated by combining bottom-up and top-down approaches. Guest-free and CO₂ guest-filled aperiodic cages, up to the gradually CO₂ occupation of the entire sI periodic unit cells were considered. Saturation, cohesive and binding energies for the systems are determined by employing a variety of density functionals and their performance is assessed. The dispersion corrections on the non-covalent interactions are found to be important in the stabilization of the CO₂ @sI energies, with the encapsulation of the CO₂ into guest-free/empty cage/lattice being always an energetically favorable process for most of the functionals studied. The PW86PBE functional with XDM or D3(BJ) dispersion corrections predicts a lattice constant in accord to the experimental values available, and simultaneously provides a reliable description for the guest-host interactions in the periodic CO₂ @sI crystal, as well as the energetics of its progressive single cage occupancy process. It has been found that the preferential orientation of the single CO₂ in the large sI crystal cages has a stabilizing effect on the hydrate, concluding that the CO₂ @sI structure is favored either by considering the individual building block cages or the complete sI unit cell crystal. Such benchmark and methodology cross-check studies benefit new data-driven model research by providing high-quality training information, with new insights that indicate the underlying factors governing their structure-driven stability, and triggering further investigations for controlling the stabilization of these promising long-term CO₂ storage materials.

Mechanical creep instability of nanocrystalline methane hydrates

Mechanical creep behaviors of natural gas hydrates are of importance for understanding the mechanical instability of gas hydrate-bearing sediments on Earth. Limited by the experimental challenges, intrinsic creep mechanisms of nanocrystalline methane hydrates remain largely unknown yet at the molecular scale. Herein, using large-scale molecular dynamics simulations, mechanical creep behaviors of nanocrystalline methane hydrates are investigated. It is revealed that mechanical creep responses are greatly dictated by internal microstructures of crystalline grain size and external conditions of temperature and static stress. Interestingly, a long steady-state creep is observed in nanocrystalline methane hydrates, which can be described by a modified constitutive Bird-Dorn-Mukherjee model. Microstructural analysis shows that deformations of crystalline grains, grain boundary diffusion and grain boundary sliding collectively govern the mechanical creep behaviors of nanocrystalline methane hydrates. Furthermore, structural transformation also appears to be important in their mechanical creep behaviors. This study provides new insights into understanding the mechanical creep scenarios of gas hydrates.

Computational density-functional approaches on finite-size and guest-lattice effects in CO2@sII clathrate hydrate

We performed first-principles computations to investigate guest-host/host-host effects on the encapsulation of the CO₂ molecule in sII clathrate hydrates from finite-size clusters up to periodic 3D crystal lattice systems. Structural and energetic properties were first computed for the individual and first-neighbors clathrate-like sII cages, where highly accurate ab initio quantum chemical methods are available nowadays, allowing in this way the assessment of the density functional (DFT) theoretical approaches employed. The performance of exchange-correlation functionals together with recently developed dispersion-corrected schemes was evaluated in describing interactions in both short-range and long-range regions of the potential. On this basis, structural relaxations of the CO₂-filled and empty sII unit cells yield lattice and compressibility parameters comparable to experimental and previous theoretical values available for sII hydrates. According to these data, the CO₂ enclathration in the sII clathrate cages is a stabilizing process, either by considering both guest-host and host-host interactions in the complete unit cell or only the guest-water energies for the individual clathrate-like sII cages. CO₂@sII clathrates are predicted to be stable whatever the dispersion correction applied and in the case of single cage occupancy are found to be more stable than the CO₂@sI structures. Our results reveal that DFT approaches could provide a good reasonable description of the underlying interactions, enabling the investigation of formation and transformation processes as a function of temperature and pressure.

Characterization of the Clathrate Hydrate Formed with Fluoromethane and Pinacolone: The Thermodynamic Stability and Volumetric Behavior of the Structure H Binary Hydrate

To reveal the relation of guest dynamics within the structure H clathrate hydrate and its macroscopic physical properties, experimental and computational works have been conducted on the system of fluoromethane (HFC-41) and pinacolone coexisting with water. The phase boundaries of the hydrate formed from HFC-41 and pinacolone within the pressure range of (0.25-2.48) MPa and the temperature range of (277-293) K were measured. The equilibrium hydrate formation pressure incorporating HFC-41 was lowered by adding the pinacolone as a large guest molecule compound to form a sH phase compared to the HFC-41 single hydrate. Powder X-ray diffraction measurements confirmed the formation of the structure H hydrate with the HFC-41 and pinacolone binary hydrate. The lattice constants of the sH hydrate were also measured to see the effect of the help guest molecular size, which showed a different trend from that of the previous studies of sH pinacolone hydrates. Molecular dynamics simulations of the binary sH phase indicate weak hydrogen bonding of the pinacolone molecules with the water in the cages in the phase with HFC-41. The oblate HFC-41 molecules showed strong orientational preference to the equatorial planes of the D' cages, which may explain some of the trends in the behavior of this phase.

Dissociation kinetics of propane-methane and butane-methane hydrates below the melting point of ice

Understanding the dissociation mechanism of gas hydrates below the melting point of ice is crucial for expanding the practical applications of solid hydrates in gas storage. The kinetic processes for gas hydrates have not been clarified, except for those of pure CH4 hydrate and CO2 hydrates. In this study, using in situ X-ray diffraction analysis, the low-temperature onset of the dissociation of C3H8 and C4H10 hydrate fine particles encapsulating CH4 as a secondary guest was investigated during temperature ramping. At ∼200 K, the C3H8 + CH4 hydrate, n-C4H10 + CH4 hydrate, and iso-C4H10 + CH4 hydrate all dissociated in a single step, similar to pure C3H8 and pure iso-C4H10 hydrate. The dissociation of C3H8 hydrate was also found to accelerate the dissociation of CH4 hydrate. Based on the experimental results, it was confirmed that the C3H8 and C4H10 molecules released from the dissociating hydrates accelerated hydrate dissociation.

Managing hydrogen bonding in the clathrate hydrate of the 1-pentanol guest molecule

1-pentanol, long-chain alcohol molecule, can be encaged in the clathrate hydrate by managing the destabilizing influence of guest–host hydrogen bonding.

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.

Temperature effects on the C-H symmetric stretching vibrational frequencies of guest hydrocarbon molecules in 512, 51262 and 51264 cages of sI and sII clathrate hydrates

C–H symmetric stretching vibrational frequencies of CH₄, C₂H₄ and C₂H₆ molecules encapsulated in 5¹², 5¹²6² and 5¹²6⁴ cages of structures I (sI) and II (sII) clathrate hydrates measured by Raman spectroscopy in the temperature range of 93–183 K was analysed. The slopes of the symmetric stretch vibrational frequencies under changing temperatures (Δv/ΔT) for CH₄, C₂H₄ and C₂H₆ molecules encapsulated in sII 5¹²6⁴ cages were smaller than those for molecules in sI 5¹²6² cages, although sI 5¹²6² cages are smaller than sII 5¹²6⁴ cages. We compared the results of Δv/ΔT in this study with the geometrical properties of each host water cage, and these comparisons suggest that the geometry of host water cages affects Δv/ΔT.

Temperature effects on the C–H symmetric stretch of hydrocarbons in various cages of sI and sII clathrate hydrates were observed.

Graphene-Based Kinetic Promotion of Gas Hydrate Formation

Gas hydrate technology holds great potential in energy and environmental fields, and achieving efficient gas hydrate formation is critical for its industrial application. Graphene is a novel carbon-based nanostructured material with excellent thermal conductivity and a large specific surface area. Therefore, the use of graphene-based materials for the promotion of gas hydrate formation might be feasible and has aroused a lot of interests. Accordingly, to evaluate the current research on graphene-based promotion of gas hydrate formation, this work presents a review of existing studies involving graphene-based promoters of gas hydrate formation. Here, the studies applying various types of graphene-based promoters for gas hydrate formation are listed and detailed, the peculiar properties of graphene-based promoters are discussed, and the promotion mechanisms are analyzed. Through this review, comprehensive insight into graphene-based promotion of gas hydrate formation can be obtained, which can guide the design and applications of novel graphene-based promoters and might contribute to achieving efficient gas hydrate formation.

Effect of temperature and large guest molecules on the C-H symmetric stretching vibrational frequencies of methane in structure H and I clathrate hydrates

Large molecules such as 2-methylbutane (C₅H₁₂) or 2,2-dimethylbutane (C₆H₁₄) form structure H (sH) hydrates with methane (CH₄) as a help gas. In this study, the Raman spectra of the C–H symmetric stretch region of CH₄ enclathrated within various sH hydrates and structure I CH₄ hydrates were analyzed in the temperature range 137.7–205.4 K. Thermal expansions of these sH hydrate samples were also measured using powder X-ray diffraction. Symmetric stretch vibrational frequencies of CH₄ in host–water cages increased because of varying temperature, and the sizes of the host–water cages also increased; variation of CH₄ in small cages was less than in larger cages. Comparing the variations of the C–H symmetric stretching frequencies of CH₄ in gas hydrates with varying pressure and temperature, we suggest that the observed trend is caused by thermal vibrations of the CH₄ molecule in water cages.

Temperature effect on C–H symmetric stretching frequencies of CH₄ in water cages of sI and sH clathrate hydrates were clarified.

Effects of cage occupancy in hydrate by first-principles calculation and modification of the van der Waals-Platteeuw hypothesis

The changes of methane hydrate lattice with the decrease of cage occupancy were calculated by first-principles methods. The calculation results show that the decrease of the cages occupancy in sII and sH hydrate does not lead to large deformation in the lattice. Even if all the methane molecules are removed so that the hydrates have become new types of ice, the sII and sH lattices remain stable. The same conclusion is also true when the occupancy of the small cages in sI hydrate is reduced. However, the sI hydrate lattice will deform and almost collapse as the large cage occupancy decreases. These calculation results suggest that sI hydrate cannot exist with empty cages. Since the van der Waals-Platteeuw theory is based on the assumption that the stability of host lattice is independent of the occupancy of guest molecule, it would be applicable to sII and sH lattices, but not to sI hydrates. We propose a modification to the van der Waals-Platteeuw hypothesis so that the theory seems more reasonable.

Structure and Density Comparison of Noble Gas Hydrates Encapsulating Xenon, Krypton and Argon

Understanding the effect of guest species on the host framework is important for the development of structure-based properties of inclusion compounds. Herein, the crystal structures of the noble gas hydrates encapsulating Xe, Kr, and Ar were studied by powder X-ray diffraction measurements. The crystal structures and hydration numbers of these noble gas hydrates were solved by Rietveld refinements using optimized models with the direct-space technique. It was revealed that host cage size of these hydrates changed depending on the type of guest species even though their unit-cell parameters were the same. Based on the structure models obtained, the densities of Xe, Kr, and Ar gas hydrates were also determined to be 1.837, 1.445 and 1.097 g/cm³ at 93 K, respectively. Our findings, from a crystallographic point of view, may give insight into further understanding the thermodynamic stability and physical properties of gas hydrates encapsulating small guests.

A Systematic Protocol for Benchmarking Guest-Host Interactions by First-Principles Computations: Capturing CO2 in Clathrate Hydrates

Clathrate hydrates of CO₂ have been proposed as potential molecular materials in tackling important environmental problems related to greenhouse gases capture and storage. Despite the increasing interest in such hydrates and their technological applications, a molecular-level understanding of their formation and properties is still far from complete. Modeling interactions is a challenging and computationally demanding task, essential to reliably determine molecular properties. First-principles calculations for the CO₂ guest in all sI, sII, and sH clathrate cages were performed, and the nature of the guest-host interactions, dominated by both hydrogen-bond and van der Waals forces, was systematically investigated. Different families of density functionals, as well as pairwise CO₂ @H₂ O model potentials versus wavefunction-based quantum approaches were studied for CO₂ clathrate-like systems. Benchmark energies for new distance-dependent datasets, consisting of potential energy curves sampling representative configurations of the systems at the repulsive, near-equilibrium, and asymptotic/long-range regions of the full-dimensional surface, were generated, and a general protocol was proposed to assess the accuracy of such conventional and modern approaches at minimum and non-minimum orientations. Our results show that dispersion interactions are important in the guest-host stabilization energies of such clathrate cages, and the encapsulation of the CO₂ into guest-free clathrate cages is always energetically favorable. In addition, the orientation of CO₂ inside each cage was explored, and the ability of current promising approaches to accurately describe non-covalent CO₂ @H₂ O guest-host interactions in sI, sII, and sH clathrates was discussed, providing information for their applicability to future multiscale computer simulations.

Retracted Article: Effect of temperature and large guest molecules on the C–H symmetric stretching vibrational frequencies of methane in structure H and I clathrate hydrates

Large molecules such as 2-methylbutane (C₅H₁₂) or 2,2-dimethylbutane (C₆H₁₄) form structure H (sH) hydrates with methane (CH₄) as a help gas. In this study, the Raman spectra of the C–H symmetric stretch region of CH₄ enclathrated within various sH hydrates and structure I CH₄ hydrates were analyzed in the temperature range 83–183 K. Thermal expansions of these sH hydrate samples were also measured using powder X-ray diffraction. Symmetric stretch vibrational frequencies of CH₄ in host water cages increased because of varying temperature, and the sizes of the host water cages also increased; variation of CH₄ in small cages was less than in larger cages. Comparing the variations of the C–H symmetric stretching frequencies of CH₄ in gas hydrates with varying pressure and temperature, we suggest that the observed trend is caused by thermal vibrations of the CH₄ molecule in water cages.

Temperature effects on C–H symmetric stretching frequencies of CH₄ in water cages of sI and sH clathrate hydrates were clarified.

Elasticity and Stability of Clathrate Hydrate: Role of Guest Molecule Motions

Molecular dynamic simulations were performed to determine the elastic constants of carbon dioxide (CO₂) and methane (CH₄) hydrates at one hundred pressure–temperature data points, respectively. The conditions represent marine sediments and permafrost zones where gas hydrates occur. The shear modulus and Young’s modulus of the CO₂ hydrate increase anomalously with increasing temperature, whereas those of the CH₄ hydrate decrease regularly with increase in temperature. We ascribe this anomaly to the kinetic behavior of the linear CO₂ molecule, especially those in the small cages. The cavity space of the cage limits free rotational motion of the CO₂ molecule at low temperature. With increase in temperature, the CO₂ molecule can rotate easily, and enhance the stability and rigidity of the CO₂ hydrate. Our work provides a key database for the elastic properties of gas hydrates, and molecular insights into stability changes of CO₂ hydrate from high temperature of ~5 °C to low decomposition temperature of ~−150 °C.

Confinement of hydrogen and hydroxyl radicals in water cages: a density functional theory study

Radicals can be confined in water cages and exhibit similar structures and properties to their corresponding free forms.

Effects of the CO₂ Guest Molecule on the sI Clathrate Hydrate Structure

This paper analyzes the structural, energetic and mechanical properties of carbon dioxide hydrate clathrates calculated using finite cluster and periodic ab initio density-functional theory methodologies. Intermolecular interactions are described by the exchange-hole dipole moment method. The stability, gas saturation energetics, guest–host interactions, cage deformations, vibrational frequencies, and equation of state parameters for the low-pressure sI cubic phase of the CO₂@H₂O clathrate hydrate are presented. Our results reveal that: (i) the gas saturation process energetically favors complete filling; (ii) carbon dioxide molecules prefer to occupy the larger of the two cages in the sI structure; (iii) blue shifts occur in both the symmetric and antisymmetric stretching frequencies of CO₂ upon encapsulation; and (iv) free rotation of guest molecules is restricted to a plane parallel to the hexagonal faces of the large cages. In addition, we calculate the librational frequency of the hindered rotation of the guest molecule in the plane perpendicular to the hexagonal faces. Our calculated spectroscopic data can be used as signatures for the detection of clathrate hydrates in planetary environments.

CO₂ processing and hydration of fruit and vegetable tissues by clathrate hydrate formation

CO2 hydrate can be used to preserve fresh fruits and vegetables, and its application could contribute to the processing of carbonated frozen food. We investigated water transformation in the frozen tissue of fresh grape samples upon CO2 treatment at 2-3 MPa and 3°C for up to 46 h. Frozen fresh bean, radish, eggplant and cucumber samples were also investigated for comparison. X-ray diffraction indicated that after undergoing CO2 treatment for several hours, structure I CO2 hydrate formed within the grape tissue. Phase-contrast X-ray imaging using the diffraction-enhanced imaging technique revealed the presence of CO2 hydrate within the intercellular spaces of these tissues. The carbonated produce became effervescent because of the dissociation of CO2 hydrate through the intercellular space, especially above the melting point of ice. In addition, suppressed metabolic activity resulting from CO2 hydrate formation, which inhibits water and nutrient transport through intercellular space, can be expected.

Disorder of Hydrofluorocarbon Molecules Entrapped in the Water Cages of Structure I Clathrate Hydrate

Water versus fluorine: Clathrate hydrates encaging hydrofluorocarbons as guests show both isotropic and anisotropic distributions within host water cages, depending on the number of fluorine atoms in the guest molecule; this is caused by changes in intermolecular interactions to host water molecules in the hydrates.

Introducing new half-integer quadrupolar nuclei for solid state NMR of inclusion compounds

Broad developments in experimental NMR techniques have opened new and exciting opportunities for application of solid state nuclear magnetic resonance (SS NMR) in studies of gas hydrates and inclusion compounds in general. Perhaps the most important advance of the last 10 years was the extension into very high magnetic fields beyond 20 T. This progress is especially significant in studies concerned with low-γ, low natural abundance, and quadrupolar nuclei. This work reports our recent exploration of clathrate hydrates and other inclusion compounds (β-quinol, tert-Bu-Calix[4], and dodecasil-3C) with SS NMR of nuclei that were not so long ago completely out of reach for NMR, namely 131Xe, 83Kr, and 33S. Although 129Xe is a widely used NMR probe, applications of the low-γ isotope 131Xe were very scarce. Being a quadrupolar spin 3/2 nucleus, 131Xe provides an additional probe for sampling the electric field gradients in inclusion compounds. Another nucleus that has been seriously under-explored is 83Kr, with its very low γ being the main obstacle, and along with quadrupolar coupling we report the first detection of the chemical shift anisotropy in krypton. The relative values of the Sternheimer antishielding factors for 131Xe and 83Kr, obtained by comparison of the spectra of the two in identical cage environments, are also discussed. Though 33S NMR of solids is notoriously difficult due to its low γ, low natural abundance, and relatively large quadrupolar moment, working at the field of 21.1 T it was possible to acquire, in a reasonable time, natural abundance 33S SS NMR spectra of various H2S and SO2 gas hydrates and inclusion compounds. In most cases the spectra are dominated by the quadrupolar interactions, providing information on the symmetry of the cages encapsulating the guest molecules, and also show the effects of very rapid reorientation of the encaged H2S and SO2. The impact of the introduction of new NMR nuclei on hydrate research is discussed.

Crystal engineering and characterization of a structure-H ionic clathrate hydrate

Ionic hydrates are known to form numerous clathrate structures in which either the cations or anions sit in cages and the counterions are incorporated into the water framework. Due to the inclusion of the ionic species, such ionic clathrate hydrates not only show many peculiar features such as metal ion encagement and superoxide ion generation, but also exhibit notable physicochemical properties such as outstanding ionic conductivity and thermal stability. Thus, the ionic clathrate hydrates are considered for their potential applicability in various fields, including those that involve solid electrolytes, gas sensors, and energy storage. In this study, we report the design, synthesis, and characterization of the first ionic clathrate hydrate of the hexagonal structure-H (Str.H) crystal type. Diethyl-dimethyl-ammonium hydroxide hydrate was synthesized with CH4 and Xe as help gases, and the crystal structure was identified by powder X-ray diffraction analysis. Further confirmation of the formation of Str.H was obtained from Raman spectroscopy and 13C, 129Xe, and 2H solid-state NMR spectroscopy. From 13C NMR and ab initio calculations, it was shown that the quaternary ion occupies the large cage of Str.H with a conformation different from that in solution, due to constraints imposed by the dimensions of the cage. The H deficiency introduced by substitution of OH– for a water molecule appears, from 129Xe NMR, to be disordered over the framework, and, from 2H NMR, to substantially increase the rate of reorientational mobility of the D atoms in the framework, over that observed for a Str.I hydrate and for ice. The Str.H hydrates are commonly more stable than other structures, thus the present findings on the ionic Str.H clathrate hydrate may offer a new approach for improving the stability of ionic clathrate hydrates for their practical application.

A molecular dynamics study of guest-host hydrogen bonding in alcohol clathrate hydrates

Clathrate hydrates are typically stabilized by suitably sized hydrophobic guest molecules. However, it has been experimentally reported that isomers of amyl-alcohol C5H11OH can be enclosed into the 5(12)6(4) cages in structure II (sII) clathrate hydrates, even though the effective radii of the molecules are larger than the van der Waals radii of the cages. To reveal the mechanism of the anomalous enclathration of hydrophilic molecules, we performed ab initio and classical molecular dynamics simulations (MD) and analyzed the structure and dynamics of a guest-host hydrogen bond for sII 3-methyl-1-butanol and structure H (sH) 2-methyl-2-butanol clathrate hydrates. The simulations clearly showed the formation of guest-host hydrogen bonds and the incorporation of the O-H group of 3-methyl-1-butanol guest molecules into the framework of the sII 5(12)6(4) cages, with the remaining hydrophobic part of the amyl-alcohol molecule well accommodated into the cages. The calculated vibrational spectra of alcohol O-H bonds showed large frequency shifts due to the strong guest-host hydrogen bonding. The 2-methyl-2-butanol guests form strong hydrogen bonds with the cage water molecules in the sH clathrate, but are not incorporated into the water framework. By comparing the structures of the alcohols in the hydrate phases, the effect of the location of O-H groups in the butyl chain of the guest molecules on the crystalline structure of the clathrate hydrates is indicated.