Ageing and Langmuir Behavior of the Cage Occupancy in the Nitrogen Gas Hydrate

Heat Transport in Clathrate Hydrates Controlled by Guest Frequency and Host-Guest Interaction

The underlying mechanism of common limited lattice thermal conductivity (κ) in energy-related host-guest crystalline compounds has been an ongoing topic in recent decades. Here, the guest-triggered intrinsic ultralow κ of the representative xenon clathrate hydrate was investigated using the time domain thermoreflectance technique and theoretical calculations. The localized guest modes were observed to hybridize with acoustic branches and severely limit the acoustic κ contribution. Besides, the strong mode coupling enables the reshaping of the overall lattice dynamics, especially for optical branches. More importantly, we identified that guest fillers prompt great phonon scattering in wide frequencies, which originates from both the guest-frequency-controlled enhancement of phase space and the host-guest-interaction-governed lattice anharmonicity. The extremely low guest frequency and strong host-guest interaction and coupling were thereby underlined to play vital but distinct roles in κ minimization. Our results unveil the dominant factors of guest reduction effects and facilitate the design of efficient thermoelectric or other thermal-related materials.

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

Unraveling the metastability of the SI and SII carbon monoxide hydrate with a combined DFT-neutron diffraction investigation

Clathrate hydrates are crystalline compounds consisting of water molecules forming cages (so-called "host") inside of which "guest" molecules are encapsulated depending on the thermodynamic conditions of formation (systems stable at low temperature and high pressure). These icelike systems are naturally abundant on Earth and are generally expected to exist on icy celestial bodies. Carbon monoxide hydrate might be considered an important component of the carbon cycle in the solar system since CO gas is one of the predominant forms of carbon. Intriguing fundamental properties have also been reported: the CO hydrate initially forms in the sI structure (kinetically favored) and transforms into the sII structure (thermodynamically stable). Understanding and predicting the gas hydrate structural stability then become essential. The aim of this work is, thereby, to study the structural and energetic properties of the CO hydrate using density functional theory (DFT) calculations together with neutron diffraction measurements. In addition to the comparison of DFT-derived structural properties with those from experimental neutron diffraction, the originality of this work lies in the DFT-derived energy calculations performed on a complete unit cell (sI and sII) and not only by considering guest molecules confined in an isolated water cage (as usually performed for extracting the binding energies). Interestingly, an excellent agreement (within less than 1% error) is found between the measured and DFT-derived unit cell parameters by considering the Perdew-Burke-Ernzerhof (denoted PBE) functional. Moreover, a strategy is proposed for evaluating the hydrate structural stability on the basis of potential energy analysis of the total nonbonding energies (i.e., binding energy and water substructure nonbonding energy). It is found that the sII structure is the thermodynamically stable hydrate phase. In addition, increasing the CO content in the large cages has a stabilizing effect on the sII structure, while it destabilizes the sI structure. Such findings are in agreement with the recent experimental results evidencing the structural metastability of the CO hydrate.