Unusual species of methane hydrate detected in nanoporous media using solid state 13C NMR.

Methane is considered to be a cubic structure I (CS-I) clathrate hydrate former, although in a number of instances, small amounts of structure II (CS-II) clathrate hydrate have been transiently observed as well. In this work, solid-state magic angle spinning 13C NMR spectra of methane hydrate formed at low temperatures inside silica-based nanoporous materials with pores in the range of 3.8-20.0 nm (CPG-20, Vycor, and MCM-41) show methane in several different environments. In addition to methane encapsulated in the dodecahedral 512 (D) and tetrakaidecahedral 51262 (T) cages typical of the CS-I clathrate hydrate phase, methane guests in pentakaidecahedral 51263 (P) and hexakaidecahedral 51264 (H) cages are also identified, and these appear to be stabilized for extended periods of time. The ratio of methane guests among the D and T cages determined from the line intensities is significantly different from that of bulk CS-I samples and indicates that both CS-I and CS-II are present as the dominant species. This is the first observation of methane in P cages, and the possible structures in which they could be present are discussed. Broad and relatively strong methane peaks, which are also observed in the spectra, can be related to methane dissolved in an amorphous component of water adjacent to the pore walls. Nanoconfinement and interaction with the pore walls clearly have a strong influence on the hydrate formed and may reflect species present in the early stages of hydrate growth.

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 CO2 below about 140 K. The phase transition was characterized by powder X-ray diffraction measurements at variable temperatures. A dynamical ordering of the CO2 guests in small pentagonal dodecahedral 512 host water cages, not previously observed in the simple CO2 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 I41/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.

Molecular dynamics simulations of interfacial structure, dynamics, and interfacial tension of tetrabutylammonium bromide aqueous solution in the presence of methane and carbon dioxide.

The interfacial behavior of tetrabutylammonium bromide (TBAB) aqueous solutions in the absence of gas and the presence of methane and carbon dioxide gases is studied by molecular dynamics simulations. The aqueous TBAB phase, at concentrations similar to the solid semiclathrate hydrate (1:38 mol ratio), has a smaller interfacial tension and an increase in the gas molecules adsorbed at the interface compared to that in pure water. Both these factors may contribute to facilitating the uptake of the gases into the solid phase during the process of semiclathrate hydrate formation. At similar gas pressures, CO2 is adsorbed preferentially compared to CH4, giving it a higher surface density, due to the stronger intermolecular interactions of CO2 molecules of the solution at the interface. The increase in relative adsorption of CH4 at the solution surface compared to that in pure water surface is due to the hydrophobic interactions between the n-alkyl chains of the TBA+ cation and methane gas.

Comment on "Cage occupancy of methane clathrate hydrates in the ternary H2O-NH3-CH4 system" by C. Petuya, M. Choukroun, T. H. Vu, A. Desmedt, A. G. Davies, and C. Sotin, Chem. Commun., 2020, 56, 12391.

Alternative interpretations of the experimental results given in the Communication of Petuya et al.1 are presented. There is evidence that under certain conditions, ammonia can be incorporated into clathrate hydrate cages.

Methane Clathrate Formation is Catalyzed and Kinetically Inhibited by the Same Molecule: Two Facets of Methanol.

Here, we perform molecular dynamics simulations to provide atomic-level insights into the dual roles of methanol in enhancing and delaying the rate of methane clathrate hydrate nucleation. Consistent with experiments, we find that methanol slows clathrate hydrate nucleation above 250 K but promotes clathrate formation at temperatures below 250 K. We show that this behavior can be rationalized by the unusual temperature dependence of the methane-methanol interaction in an aqueous solution, which emerges due to the hydrophobic effect. In addition to its antifreeze properties at temperatures above 250 K, methanol competes with water to interact with methane prior to the formation of clathrate nuclei. Below 250 K, methanol encourages water to occupy the space between methane molecules favoring clathrate formation and it may additionally promote water mobility.

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.

Thermophysical Property Measurements of Tetrabutylphosphonium Oxalate (TBPOx) Ionic Semiclathrate Hydrate as a Media for the Thermal Energy Storage System.

With increasing global power demand, thermal energy storage technology could play a role ensuring a sustainable energy supply in power generation from renewable energy sources and power demand concentration. Hydrates have high potential as phase change materials (PCMs) for the use as a thermal energy storage medium. To develop thermal energy storage technology using a hydrate-based material, further investigation of thermophysical properties and the selection of a suitable hydrate are required. Tetrabutylphosphonium oxalate (TBPOx) ionic semiclathrate hydrate contains oxalic acid in salt form, as a guest compound, which is classified as carboxylic acid group with low environmental impact. In the present study, the phase equilibrium temperature and the dissociation heat of TBPOx hydrate were measured. The highest equilibrium temperature of the solid hydrate formed was 9.4°C at the mass fraction 0.35 of TBPOx in aqueous solution. The largest dissociation heat was 186.0 ± 0.5 kJ·kg-1 at the mass fraction of 0.35. Comparing with other PCMs with close phase equilibrium temperatures, TBPOx hydrate is superior in safety and sustainability. These results indicate that TBPOx hydrate would be suitable as the thermal storage medium for the general air conditioning systems.

Molecular dynamics simulations of interfacial properties of the CO2-water and CO2-CH4-water systems.

Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO2 and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO2-water interface has a lower interfacial tension than the CH4-water interface. A greater number of CO2 molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO2 and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO2. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO2-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO2 capture and sequestration processes are discussed.

Molecular dynamics simulations of nano-confined methanol and methanol-water mixtures between infinite graphite plates: Structure and dynamics.

Molecular dynamics simulations are used to investigate microscopic structures and dynamics of methanol and methanol-water binary mixture films confined between hydrophobic infinite parallel graphite plate slits with widths, H, in the range of 7-20 Å at 300 K. The initial geometric densities of the liquids were chosen to be the same as bulk methanol at the same temperature. For the two narrowest slit widths, two smaller initial densities were also considered. For the nano-confined system with H = 7 Å and high pressure, a solid-like hexagonal arrangement of methanol molecules arranged perpendicular to the plates is observed which reflects the closest packing of the molecules and partially mirrors the structure of the underlying graphite structure. At lower pressures and for larger slit widths, in the contact layer, the methanol molecules prefer having the C-O bond oriented parallel to the walls. Layered structures of methanol parallel to the wall were observed, with contact layers and additional numbers of central layers depending on the particular slit width. For methanol-water mixtures, simulations of solutions with different composition were performed between infinite graphite slits with H = 10 and 20 Å at 300 K. For the nanoslit with H = 10 Å, in the solution mixtures, three layers of molecules form, but for all mole fractions of methanol, methanol molecules are excluded from the central fluid layer. In the nanopore with H = 20 Å, more than three fluid layers are formed and methanol concentrations are enhanced near the confining plates walls compared to the average solution stoichiometry. The self-diffusion coefficients of methanol and water molecules in the solution show strong dependence on the solution concentration. The solution mole fractions with minimal diffusivity are the same in confined and non-confined bulk methanol-water mixtures.

Interfacial properties of hydrocarbon/water systems predicted by molecular dynamic simulations.

The presence of small hydrocarbons is known to reduce the interfacial tension of the gas-water interface, and this phenomenon can affect the formation of the clathrate hydrates of these gases. In this work, the interfacial behavior of the pure methane-, ethane-, and propane-water, and the ternary 90:7:3 mol. % gas mixture of (methane + ethane + propane)-water were studied with molecular dynamics simulations. The interfacial tension, γ, and z-density profiles for the gases and water from simulations of the gas-water systems were determined at the temperatures of 275.15 and 298.15 K, and pressures up to 10 MPa for methane and up to near the experimental saturation pressures of ethane and propane. The goal is to accurately calculate the interfacial tension for the hydrocarbon/water systems and to analyze the molecular behaviors at the interfaces which lead to the observed trends. At the same hydrostatic gas phase pressure, propane, ethane, and methane reduce the gas-water interfacial tension in that order. The local density of the gas molecules at the interface is enhanced relative to the bulk gas, and it was determined that about 13%-20%, 33%-40%, and 54%-59% of the gas molecules in the simulation congregated at the interfaces for the CH4-, C2H6-, and C3H8-water systems, respectively, at the different simulated hydrostatic pressure ranges. For all gases in the pressure range studied, a complete monolayer of gas had not formed at the water interface. Furthermore, a dynamic equilibrium with fast exchange between molecules at the interface and in the gas phase was observed. For the gas mixture, deviations were observed between total calculated interfacial tension, γmix, and the "ideal mixture" value, ∑xiγi,pure, calculated from the interfacial tensions of the pure gases, where xi is the mole fraction of each substance in the simulation. Some possible implications of the results on the mechanism of clathrate hydrate formation are discussed.

The anomalous halogen bonding interactions between chlorine and bromine with water in clathrate hydrates.

Clathrate hydrate phases of Cl2 and Br2 guest molecules have been known for about 200 years. The crystal structure of these phases was recently re-determined with high accuracy by single crystal X-ray diffraction. In these structures, the water oxygen-halogen atom distances are determined to be shorter than the sum of the van der Waals radii, which indicates the action of some type of non-covalent interaction between the dihalogens and water molecules. Given that in the hydrate phases both lone pairs of each water oxygen atom are engaged in hydrogen bonding with other water molecules of the lattice, the nature of the oxygen-halogen interactions may not be the standard halogen bonds characterized recently in the solid state materials and enzyme-substrate compounds. The nature of the halogen-water interactions for the Cl2 and Br2 molecules in two isolated clathrate hydrate cages has recently been studied with ab initio calculations and Natural Bond Order analysis (Ochoa-Resendiz et al. J. Chem. Phys. 2016, 145, 161104). Here we present the results of ab initio calculations and natural localized molecular orbital analysis for Cl2 and Br2 guests in all cage types observed in the cubic structure I and tetragonal structure I clathrate hydrates to characterize the orbital interactions between the dihalogen guests and water. Calculations with isolated cages and cages with one shell of coordinating molecules are considered. The computational analysis is used to understand the nature of the halogen bonding in these materials and to interpret the guest positions in the hydrate cages obtained from the X-ray crystal structures.

Understanding decomposition and encapsulation energies of structure I and II clathrate hydrates.

When compressed with water or ice under high pressure and low temperature conditions, some gases form solid gas hydrate inclusion compounds which have higher melting points than ice under those pressures. In this work, we study the balance of the guest-water and water-water interaction energies that lead to the formation of the clathrate hydrate phases. In particular, molecular dynamics simulations with accurate water potentials are used to study the energetics of the formation of structure I (sI) and II (sII) clathrate hydrates of methane, ethane, and propane. The dissociation enthalpy of the clathrate hydrate phases, the encapsulation enthalpy of methane, ethane, and propane guests in the corresponding phases, and the average bonding enthalpy of water molecules are calculated and compared with accurate calorimetric measurements and previous classical and quantum mechanical calculations, when available. The encapsulation energies of methane, ethane, and propane guests stabilize the small and large sI and sII hydrate cages, with the larger molecules giving larger encapsulation energies. The average water-water interactions are weakened in the sI and sII phases compared to ice. The relative magnitudes of the van der Waals potential energy in ice and the hydrate phases are similar, but in the ice phase, the electrostatic interactions are stronger. The stabilizing guest-water "hydrophobic" interactions compensate for the weaker water-water interactions and stabilize the hydrate phases. A number of common assumptions regarding the guest-cage water interactions are used in the van der Waals-Platteeuw statistical mechanical theory to predict the clathrate hydrate phase stability under different pressure-temperature conditions. The present calculations show that some of these assumptions may not accurately reflect the physical nature of the interactions between guest molecules and the lattice waters.

Phase Transition of a Structure†II Cubic Clathrate Hydrate to a Tetragonal Form.

The crystal structure and phase transition of cubic structure II (sII) binary clathrate hydrates of methane (CH4 ) and propanol are reported from powder X-ray diffraction measurements. The deformation of host water cages at the cubic-tetragonal phase transition of 2-propanol+CH4 hydrate, but not 1-propanol+CH4 hydrate, was observed below about 110 K. It is shown that the deformation of the host water cages of 2-propanol+CH4 hydrate can be explained by the restriction of the motion of 2-propanol within the 5(12) 6(4) host water cages. This result provides a low-temperature structure due to a temperature-induced symmetry-lowering transition of clathrate hydrate. This is the first example of a cubic structure of the common clathrate hydrate families at a fixed composition.

Selective occupancy of methane by cage symmetry in TBAB ionic clathrate hydrate.

Methane trapped in the two distinct dodecahedral cages of the ionic clathrate hydrate of TBAB was studied by single crystal XRD and MD simulation. The relative CH4 occupancies over the cage types were opposite to those of CO2, which illustrates the interplay between the cage symmetry and guest shape and dynamics, and thus the gas selectivity.

Molecular simulations and density functional theory calculations of bromine in clathrate hydrate phases.

Bromine forms a tetragonal clathrate hydrate structure (TS-I) very rarely observed in clathrate hydrates of other guest substances. The detailed structure, energetics, and dynamics of Br2 and Cl2 in TS-I and cubic structure I (CS-I) clathrate hydrates are studied in this work using molecular dynamics and quantum chemical calculations. X-ray diffraction studies show that the halogen-water-oxygen distances in the cages of these structures are shorter than the sum of the van der Waals radii of halogen and oxygen atoms. This suggests that the stabilizing effects of halogen bonding or other non-covalent interactions (NCIs) may contribute to the formation of the unique tetragonal bromine hydrate structure. We performed molecular dynamics simulations of Br2 and Cl2 clathrate hydrates using our previously developed five-site charge models for the dihalogen molecules [Dureckova et al. Can. J. Chem. 93, 864 (2015)] which reproduce the computed electrostatic potentials of the dihalogens and account for the electropositive σ-hole of the halogen bond donor (the dihalogen). Analysis of the radial distribution functions, enthalpies of encapsulation, velocity and orientation autocorrelation functions, and polar angle distributions are carried out for Br2 and Cl2 guests in various cages to contrast the properties of these guests in the TS-I and CS-I phases. Quantum chemical partial geometry optimizations of Br2 and Cl2 guests in the hydrate cages using the M06-2X functional give short halogen-water distances compatible with values observed in X-ray diffraction experiments. NCI plots of guest-cage structures are generated to qualitatively show the relative strength of the non-bonding interactions between dihalogens and water molecules. The differences between behaviors of Br2 and Cl2 guests in the hydrate cages may explain why bromine forms the unique TS-I phase.

Formation of methane nano-bubbles during hydrate decomposition and their effect on hydrate growth.

Molecular dynamic simulations are performed to study the conditions for methane nano-bubble formation during methane hydrate dissociation in the presence of water and a methane gas reservoir. Hydrate dissociation leads to the quick release of methane into the liquid phase which can cause methane supersaturation. If the diffusion of methane molecules out of the liquid phase is not fast enough, the methane molecules agglomerate and form bubbles. Under the conditions of our simulations, the methane-rich quasi-spherical bubbles grow to become cylindrical with a radius of ∼11 Å. The nano-bubbles remain stable for about 35 ns until they are gradually and homogeneously dispersed in the liquid phase and finally enter the gas phase reservoirs initially set up in the simulation box. We determined that the minimum mole fraction for the dissolved methane in water to form nano-bubbles is 0.044, corresponding to about 30% of hydrate phase composition (0.148). The importance of nano-bubble formation to the mechanism of methane hydrate formation, growth, and dissociation is discussed.

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.