Cage occupancies, lattice constants, and guest chemical potentials for structure II hydrogen clathrate hydrate from Gibbs ensemble Monte Carlo simulations

Dissociation line and driving force for nucleation of the nitrogen hydrate from computer simulation. II. Effect of multiple occupancy

In this work, we determine the dissociation line of the nitrogen (N2) hydrate by computer simulation using the TIP4P/Ice model for water and the TraPPE force field for N2. This work is the natural extension of Paper I, in which the dissociation temperature of the N2 hydrate has been obtained at 500, 1000, and 1500 bar [Algaba et al., J. Chem. Phys. 159, 224707 (2023)] using the solubility method and assuming single occupancy. We extend our previous study and determine the dissociation temperature of the N2 hydrate at different pressures, from 500 to 4500 bar, taking into account the single and double occupancy of the N2 molecules in the hydrate structure. We calculate the solubility of N2 in the aqueous solution as a function of temperature when it is in contact with a N2-rich liquid phase and when in contact with the hydrate phase with single and double occupancy via planar interfaces. Both curves intersect at a certain temperature that determines the dissociation temperature at a given pressure. We observe a negligible effect of occupancy on the dissociation temperature. Our findings are in very good agreement with the experimental data taken from the literature. We have also obtained the driving force for the nucleation of the hydrate as a function of temperature and occupancy at several pressures. As in the case of the dissociation line, the effect of occupancy on the driving force for nucleation is negligible. To the best of our knowledge, this is the first time that the effect of the occupancy on the driving force for nucleation of a hydrate that exhibits sII crystallographic structure is studied from computer simulation.

Mesomorphology of clathrate hydrates from molecular ordering

Clathrate hydrates are crystals formed by guest molecules that stabilize cages of hydrogen-bonded water molecules. Whereas thermodynamic equilibrium is well described via the van der Waals and Platteeuw approach, the increasing concerns with global warming and energy transition require extending the knowledge to non-equilibrium conditions in multiphase, sheared systems, in a multiscale framework. Potential macro-applications concern the storage of carbon dioxide in the form of clathrates, and the reduction of hydrate inhibition additives currently required in hydrocarbon production. We evidence porous mesomorphologies as key to bridging the molecular scales to macro-applications of low solubility guests. We discuss the coupling of molecular ordering with the mesoscales, including (i) the emergence of porous patterns as a combined factor from the walk over the free energy landscape and 3D competitive nucleation and growth and (ii) the role of molecular attachment rates in crystallization-diffusion models that allow predicting the timescale of pore sealing. This is a perspective study that discusses the use of discrete models (molecular dynamics) to build continuum models (phase field models, crystallization laws, and transport phenomena) to predict multiscale manifestations at a feasible computational cost. Several advances in correlated fields (ice, polymers, alloys, and nanoparticles) are discussed in the scenario of clathrate hydrates, as well as the challenges and necessary developments to push the field forward.

Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.

Modification of the van der Waals and Platteeuw model for gas hydrates considering multiple cage occupancy

This study reviews available van der Waals- and Platteeuw-based hydrate models considering multiple occupancy of cavities. Small guest molecules, such as hydrogen and nitrogen, are known to occupy lattice cavities multiple times. This phenomenon has a significant impact on hydrate stability and thermodynamic properties of the hydrate phase. The objective of this work is to provide a comprehensive overview and required correlations for the implementation of a computationally sufficient cluster model that considers up to five guest molecules per cavity. Two methodologies for cluster size estimation are evaluated by existing nitrogen hydrate models showing accurate results for phase equilibria calculations. Furthermore, a preliminary hydrogen hydrate model is introduced and compared with the results of other theoretical studies, indicating that double occupancy of small sII cavities is improbable and four-molecule clusters are predominant in large sII cavities for pressures above 300 MPa. This work lays the foundation for further exploration and optimization of hydrate-based technologies for small guest molecules, e.g., storage and transportation, emphasizing their role in the future landscape of sustainable energy solutions.

Analysing the stability of He-filled hydrates: how many He atoms fit in the sII crystal?

Clathrate hydrates have the ability to encapsulate atoms and molecules within their cavities, and thus they could be potentially large storage capacity materials. The present work studies the multiple cage occupancy effects in the recently discovered He@sII crystal. On the basis of previous theoretical and experimental findings, the stability of He(1/1)@sII, He(1/4)@sII and He(2/4)@sII crystals was analysed in terms of structural, mechanical and energetic properties. For this purpose, first-principles DFT/DFT-D computations were performed by using both semi-local and non-local functionals, not only to elucidate which configuration is the most energetically favoured, but also to scrutinize the relevance of the long-range dispersion interactions in these kinds of compounds. We have encountered that dispersion interactions play a fundamental role in describing the underlying interactions, and different tendencies were observed depending on the choice of the functional. We found that PW86PBE-XDM shows the best performance, while the non-local functionals tested here were not able to correctly account for them. The present results reveal that the most stable configuration is the one presenting singly occupied D cages and tetrahedrally occupied H cages (He(1/4)@sII) in line with the experimental observation.

Phase equilibria molecular simulations of hydrogen hydrates via the direct phase coexistence approach

We report the three-phase (hydrate-liquid water-vapor) equilibrium conditions of the hydrogen-water binary system calculated with molecular dynamics simulations via the direct phase coexistence approach. A significant improvement of ∼10.5 K is obtained in the current study, over earlier simulation attempts, by using a combination of modifications related to the hydrogen model that include (i) hydrogen Lennard-Jones parameters that are a function of temperature and (ii) the water-guest energy interaction parameters optimized further by using the Lorentz-Berthelot combining rules, based on an improved description of the solubility of hydrogen in water.

Hydrogen Inter-Cage Hopping and Cage Occupancies inside Hydrogen Hydrate: Molecular-Dynamics Analysis

The inter-cage hopping in a type II clathrate hydrate with different numbers of H2 and D2 molecules, from 1 to 4 molecules per large cage, was studied using a classical molecular dynamics simulation at temperatures of 80 to 240 K. We present the results for the diffusion of these guest molecules (H2 or D2) at all of the different occupations and temperatures, and we also calculated the activation energy as the energy barrier for the diffusion using the Arrhenius equation. The average occupancy number over the simulation time showed that the structures with double and triple large-cage H2 occupancy appeared to be the most stable, while the small cages remained with only one guest molecule. A Markov model was also calculated based on the number of transitions between the different cage types.

Gas hydrates in sustainable chemistry

This review includes the current state of the art understanding and advances in technical developments about various fields of gas hydrates, which are combined with expert perspectives and analyses.

Molecular Insights into Cage Occupancy of Hydrogen Hydrate: A Computational Study

Density functional theory calculations and molecular dynamics simulations were performed to investigate the hydrogen storage capacity in the sII hydrate. Calculation results show that the optimum hydrogen storage capacity is ~5.6 wt%, with the double occupancy in the small cage and quintuple occupancy in the large cage. Molecular dynamics simulations indicate that these multiple occupied hydrogen hydrates can occur at mild conditions, and their stability will be further enhanced by increasing the pressure or decreasing the temperature. Our work highlights that the hydrate is a promising material for storing hydrogen.