Rates and Mechanisms of Conversion of Ice Nanocrystals to Ether Clathrate Hydrates: Guest-Molecule Catalytic Effects at ∼120 K

NH3 as simple clathrate-hydrate catalyst: Experiment and theory

The catalytic action of NH₃ within the all-vapor approach for instant clathrate hydrate (CH) formation is studied using both FTIR spectroscopy and ab initio molecular dynamics simulations. A unique property of NH₃, namely, the rapid abundant penetration and occupation of the water network, creates defects, particularly Bjerrum D-defects, in the hydrate frame that are generally stabilized by guest NH₃ molecules in the cages. Furthermore, insertion of NH₃ seriously disturbs the hydrate network where the guest NH₃ molecules also make fluxional H-bonds with the host water molecules. These defects strongly facilitate a sub-second formation of the simple NH₃ s-II gas hydrate at 160 K. FTIR spectra of aerosols of the NH₃ s-II CH have been measured, and the displacement of both small and large cage NH₃ guests by CO₂ and tetrahydrofuran is examined.

Kinetics of Trifluoromethane Clathrate Hydrate Formation from CHF3 Gas and Ice Particles

We report the formation kinetics of trifluoromethane clathrate hydrate (CH) from less than 75 μm diameter ice particles and CHF₃ gas. As previously observed for difluoromethane and propane hydrate formation, the initial stages of the reaction exhibit a strong negative correlation of the reaction rate with temperature, consistent with a negative activation energy of formation. The values obtained for trifluoromethane, ca. -6 kJ/mol (H₂O), are similar to those for difluoromethane, even though the two molecules have different intermolecular interactions and sizes. The activation energy is lesser per mole of H₂O, but greater per mole of guest molecule, than for propane hydrate, which has a different crystal structure. We propose a possible explanation for the negative activation barrier based on the stabilization of metastable structures at low temperature. A pronounced dependence of the formation kinetics on the gas flow rate into the cell is observed. At 253 K and a flow rate of 15 mmol/h, the stage II enclathration of trifluoromethane proceeds so quickly that the overpressure, the difference between the gas cell pressure and the hydrate vapor pressure, is only 0.06 MPa.

Molecular Modes and Dynamics of HCl and DCl Guests of Gas Clathrate Hydrates

Recent years have yielded advances in the placement of unusual molecules as guests within clathrate hydrates (CHs) without severe distortion of the classic lattice structures. Reports describing systems for which observable but limited distortion does occur are available for methanol, ammonia, acetone, and small ether molecules. In these particular examples, the large-cage molecules often participate as non-classical guests H-bonded to the cage walls. Here, we expand the list of such components to include HCl/DCl and HBr as small-cage guests. Based on FTIR spectra of nanocrystalline CHs from two distinct preparative methods combined with critical insights derived from on-the-fly molecular dynamics and ab initio computational data, a coherent argument emerges that these strong acids serve as a source of molecular small-cage guests, ions, and orientational defects. Depending on the HCl/DCl content the ions, defects and molecular guests determine the CH structures, some of which form in sub-seconds via an all-vapor preparative method.

Catalytic activity of methanol in all-vapor subsecond clathrate-hydrate formation

Methanol's property as a catalyst in the formation of gas clathrate hydrates has been recognized for several years and was recently employed in a broad ranging study [K. Shin, K. A. Udachin, I. L. Moudrakovski, D. M. Leek, S. Alavi, C. I. Ratcliffe, and J. A. Ripmeester, Proc. Natl. Acad. Sci. U.S.A. 110, 8437 (2013)]. A new measure of that activity is offered here from comparative rates of formation of methanol (MeOH) clathrate hydrates within our all-vapor aerosol methodology for which tetrahydrofuran (THF) and other small ethers have set a standard for catalytic action. We have previously described numerous examples of the complete conversion of warm all-vapor mixtures to aerosols of gas clathrate hydrates on a sub-second time scale, generally with the catalyst confined primarily to the large cage of either structure-I (s-I) or structure-II (s-II) hydrates. THF has proven to be the most versatile catalyst for the complete subsecond conversion of water to s-II hydrate nanocrystals that follows pulsing of appropriate warm vapor mixtures into a cold chamber held in the 140-220 K range. Here, the comparative ability of MeOH to catalyze the formation of s-I hydrates in the presence of a small-cage help-gas, CO2 or acetylene, is examined. The surprising result is that, in the presence of either help gas, CH-formation rates appear largely unchanged by a complete replacement of THF by MeOH in the vapor mixtures for a chamber temperature of 170 K. However, as that temperature is increased, the dependence of effective catalysis by MeOH on the partial pressure of help gases also increases. Nevertheless, added MeOH is shown to markedly accelerate the s-II THF-CO2 CH formation rate at 220 K.

CO2 and C2H2 in cold nanodroplets of oxygenated organic molecules and water

Recent demonstrations of subsecond and microsecond timescales for formation of clathrate hydrate nanocrystals hint at future methods of control of environmental and industrial gases such as CO2 and methane. Combined results from cold-chamber and supersonic-nozzle [A. S. Bhabhe, "Experimental study of condensation and freezing in a supersonic nozzle," Ph.D. thesis (Ohio State University, 2012), Chap. 7] experiments indicate extremely rapid encagement of components of all-vapor pre-mixtures. The extreme rates are derived from (a) the all-vapor premixing of the gas-hydrate components and (b) catalytic activity of certain oxygenated organic large-cage guests. Premixing presents no obvious barrier to large-scale conditions of formation. Further, from sequential efforts of the groups of Trout and Buch, a credible defect-based model of the catalysis mechanism exists for guidance. Since the catalyst-generated defects are both mobile and abundant, it is often unnecessary for a high percentage of the cages to be occupied by a molecular catalyst. Droplets represent the liquid phase that bridges the premixed vapor and clathrate hydrate phases but few data exist for the droplets themselves. Here we describe a focused computational and FTIR spectroscopic effort to characterize the aerosol droplets of the all-vapor cold-chamber methodology. Computational data for CO2 and C2H2, hetero-dimerized with each of the organic catalysts and water, closely match spectroscopic redshift patterns in both magnitude and direction. Though vibrational frequency shifts are an order of magnitude greater for the acetylene stretch mode, both CO2 and C2H2 experience redshift values that increase from that for an 80% water-methanol solvent through the solvent series to approximately doubled values for tetrahydrofuran and trimethylene oxide (TMO) droplets. The TMO solvent properties extend to a 50 mol.% solution of CO2, more than an order of magnitude greater than for the water-methanol solvent mixture. The impressive agreement between heterodimer and experimental shift values throughout the two series encourages speculation concerning local droplet structures while the stable shift patterns appear to be useful indicators of the gas solubilities.

Antifreezes act as catalysts for methane hydrate formation from ice

Contrary to the thermodynamic inhibiting effect of methanol on methane hydrate formation from aqueous phases, hydrate forms quickly at high yield by exposing frozen water-methanol mixtures with methanol concentrations ranging from 0.6-10 wt% to methane gas at pressures from 125 bars at 253 K. Formation rates are some two orders of magnitude greater than those obtained for samples without methanol and conversion of ice is essentially complete. Ammonia has a similar catalytic effect when used in concentrations of 0.3-2.7 wt%. The structure I methane hydrate formed in this manner was characterized by powder X-ray diffraction and Raman spectroscopy. Steps in the possible mechanism of action of methanol were studied with molecular dynamics simulations of the Ih (0001) basal plane exposed to methanol and methane gas. Simulations show that methanol from a surface aqueous layer slowly migrates into the ice lattice. Methane gas is preferentially adsorbed into the aqueous methanol surface layer. Possible consequences of the catalytic methane hydrate formation on hydrate plug formation in gas pipelines, on large scale energy-efficient gas hydrate formation, and in planetary science are discussed.

Tracking all-vapor instant gas-hydrate formation and guest molecule populations: a possible probe for molecules trapped in water nanodroplets

Quantitative Fourier-transform infrared spectra for low-temperature (160-200 K) aerosols of clathrate-hydrate nanoparticles that contain large-cage catalysts and small-cage nonpolar guests have been extended to a broad range of vapor compositions and sampling conditions. The data better reveal the stages by which room-temperature vapor mixtures, when cooled below ∼220 K, instantly generate aerosols with particles composed exclusively of the corresponding clathrate hydrates. In particular the quantitative data help relate the nature of the hydrates that form to the composition of the aqueous nanodroplets of the first stages of the rapid transition from the all-vapor mixture. The overall transition from an all-vapor mixture to "gas"-hydrate nanocrystals is a multistage one that has been characterized as homogeneous nucleation and growth of solution nanodroplets (∼240 K) followed by nucleation and growth of the gas-hydrate particles (∼220 K); all occurring within a subsecond that follows pulsing of the warm vapor into a sampling cold chamber. This may serve well as a general description of the instantaneous generation of the gas-hydrate aerosols, but closer consideration of the nature of the sampling method, in context with recent computation-based insights to (a) gas-hydrate nucleation stages∕rates and (b) the lifetimes of trapped small nonpolar molecules in cold aqueous nanodroplets, suggests a more complex multistage transition. The simulated lifetimes and extensive new quantitative infrared data significantly broaden the knowledge base in which the instantaneous transition from vapor to crystalline hydrate particles is viewed. The apparent need for a high occupancy of large-cage catalytic guest molecules currently limits the practical value of the all-vapor method. Only through greater clarity in the molecular-level description of the transition will the ultimate limits be defined.

Spectroscopic methods in gas hydrate research

Gas hydrates are crystalline structures comprising a guest molecule surrounded by a water cage, and are particularly relevant due to their natural occurrence in the deep sea and in permafrost areas. Low molecular weight molecules such as methane and carbon dioxide can be sequestered into that cage at suitable temperatures and pressures, facilitating the transition to the solid phase. While the composition and structure of gas hydrates appear to be well understood, their formation and dissociation mechanisms, along with the dynamics and kinetics associated with those processes, remain ambiguous. In order to take advantage of gas hydrates as an energy resource (e.g., methane hydrate), as a sequestration matrix in (for example) CO(2) storage, or for chemical energy conservation/storage, a more detailed molecular level understanding of their formation and dissociation processes, as well as the chemical, physical, and biological parameters that affect these processes, is required. Spectroscopic techniques appear to be most suitable for analyzing the structures of gas hydrates (sometimes in situ), thus providing access to such information across the electromagnetic spectrum. A variety of spectroscopic methods are currently used in gas hydrate research to determine the composition, structure, cage occupancy, guest molecule position, and binding/formation/dissociation mechanisms of the hydrate. To date, the most commonly applied techniques are Raman spectroscopy and solid-state nuclear magnetic resonance (NMR) spectroscopy. Diffraction methods such as neutron and X-ray diffraction are used to determine gas hydrate structures, and to study lattice expansions. Furthermore, UV-vis spectroscopic techniques and scanning electron microscopy (SEM) have assisted in structural studies of gas hydrates. Most recently, waveguide-coupled mid-infrared spectroscopy in the 3-20 μm spectral range has demonstrated its value for in situ studies on the formation and dissociation of gas hydrates. This comprehensive review summarizes the importance of spectroscopic analytical techniques to our understanding of the structure and dynamics of gas hydrate systems, and highlights selected examples that illustrate the utility of these individual methods.

Nucleation pathways of clathrate hydrates: effect of guest size and solubility

Understanding the microscopic mechanism of nucleation of clathrate hydrates is important for their use in hydrogen storage, CO(2) sequestration, storage and transport of natural gas, and the prevention of the formation of hydrate plugs in oil and gas pipelines. These applications involve hydrate guests of varied sizes and solubility in water that form different hydrate crystal structures. Nevertheless, molecular studies of the mechanism of nucleation of hydrates have focused on the single class of small hydrophobic guests that stabilize the sI crystal. In this work, we use molecular dynamics simulations with a very efficient coarse-grained model to elucidate the mechanisms of nucleation of clathrate hydrates of four model guests that span a 2 orders of magnitude range in solubility in water and that encompass sizes which stabilize each one a different hydrate structure (sI and sII, with and without occupancy of the dodecahedral cages). We find that the overall mechanism of clathrate nucleation is similar for all guests and involves a first step of formation of blobs, dense clusters of solvent-separated guest molecules that are the birthplace of the clathrate cages. Blobs of hydrophobic guests are rarer and longer-lived than those for soluble guests. For each guest, we find multiple competing channels to form the critical nuclei, filled dodecahedral (5(12)) cages, empty 5(12) cages, and a variety of filled large (5(12)6(n) with n = 2, 3, and 4) clathrate cages. Formation of empty dodecahedra is an important nucleation channel for all but the smallest guest. The empty 5(12) cages are stabilized by the presence of guests from the blob in their first solvation shell. Under conditions of high supercooling, the structure of the critical and subcritical nuclei is mainly determined by the size of the guest and does not reflect the cage composition or ordering of the stable or metastable clathrate crystals.

Clathrate hydrates with hydrogen-bonding guests

Clathrate hydrates (CHs) are inclusion compounds in which "tetrahedrally" bonded H(2)O forms a crystalline host lattice composed of a periodic array of cages. The structure is stabilized by guest particles which occupy the cages and interact with cage walls via van der Waals interactions. A host of atoms or small molecules can act as guests; here the focus is on guests that are capable of strong to intermediate H-bonding to water (small ethers, H(2)S, etc.) but nevertheless "choose" this hydrate crystal form in which H-bonding is absent from the equilibrium crystal structure. These CHs can form by exposure of ice to guest molecules at temperatures as low as 100-150 K, at the (low) guest saturation pressure. This is in contrast to the "normal" CHs whose formation typically requires temperatures well above 200 K and at least moderate pressures. The experimental part of this study addresses formation kinetics of CHs with H-bonding guests, as well as transformation kinetics between different CH forms, studied by CH infrared spectroscopy. The accompanying computational study suggests that the unique properties of this family of CHs are due to exceptional richness of the host lattice in point defects, caused by defect stabilization by H-bonding of water to the guests.

Evidence for the surface origin of point defects in ice: Control of interior proton activity by adsorbates

Spectroscopic studies are presented of H-D isotopic exchange in the interior of ice nanocrystals. The exchange process is dominated by ionic and orientational defects long viewed as governing the electrical properties of ice. A new finding that interior exchange rates can be controlled by acidic and basic adsorbates is evidence that the defects originate at the ice surface. In particular, it is argued that interior isotopic exchange is a reflection of proton concentrations equilibrated at the ice surface.

Spectroscopic Signatures of Halogens in Clathrate Hydrate Cages. 1. Bromine

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