Dissociation mechanism of gas hydrates (I,II, H) of alkane molecules: a comparative molecular dynamics simulation

Employing NPT molecular dynamics method with consistent valence force field, the dissociation processes of sI, sII and sH gas hydrates are simulated at different temperatures and at a constant pressure of 100 MPa. The dissociation mechanisms of gas hydrates are revealed by analysing the structural snapshots, radial distribution functions and diffusion coefficients at different temperatures. As temperature increases, the diffusion rates of water molecules and guest molecules increase; thus the clathrate skeleton formed by water molecules with hydrogen bonds distorts and breaks down; meanwhile the guest molecules encapsulated in the water cavities are released. The size of guest molecules affects the dissociation behaviour of gas hydrate. In addition, the dissociation behaviour also relies on the structural phase of gas hydrates.     

On the Stability of Clathrate Hydrates Encaging Polar Guest Molecules: Contrast in the Hydrogen Bonds of Methylamine and Methanol Hydrates

Abstract

The stability of clathrate hydrates encaging highly polar guests has been investigated in order to explain the experimental observation that some amines form clathrate hydrates but alcohols act as inhibitor to hydrate formation. We choose methylamine and methanol as guest species and examine the stable structure, at which the total potential energy has a minimum value. At the local minima of those two hydrates, the potential energies of water-water and guest-water, and their hydrogen bonded networks are compared. It is found that methanol does not retain the host lattice structure, while the host-network structure is kept in the presence of methylamine. It is shown that the difference in the magnitude of the partial charge on the hydrogen atom between the hydroxyl and amino groups plays a much more significant role on the stability of both clathrate hydrates than the difference in molecular geometry. This is supported from the result of a methylamine-like model that has the same partial charges on the atoms in the hydrophilic site as methanol.     

The Stability of Clathrate Hydrates: Temperature Dependence of Dissociation Pressure in Xe and Ar Hydrate

Abstract

The thermodynamic stability of clathrate hydrates I and II encaging xenon or argon has been investigated by examining the temperature dependence of the dissociation pressure. The evaluation of the stability is made based on the generalized van der Waals and Platteeuw theory developed by Tanaka and Kiyohara [J. Chem. Phys. 98, 4098 (1993)]. In the new treatment, the free energy of formation of hydrates in equilibrium with ice is calculated by taking the coupling of the host lattice vibrations with guests into consideration. The predicted dissociation pressures of Xe and Ar hydrates agree well with experiments in higher temperature range. A poor agreement between experiment and calculation for Ar clathrate hydrate at low temperature is improved by the use of a quantum mechanical partition function for a harmonic oscillator in evaluating the free energy difference between ice and empty hydrate.     

Monte Carlo study of sI hydrogen hydrates

In this study, we perform grand canonical Monte Carlo simulations to evaluate the hydrogen storage capacity of structure I (sI) hydrogen hydrates at pressures up to 500 MPa. Initially, we calculate the upper limit of H₂ content of sI hydrates by studying the hypothetical sI hydrate, where H₂ is the single guest component. It is found that the storage capacity of the hypothetical pure H₂ sI hydrate could reach 3.5 wt% at 500 MPa and 274 K. Depending on pressure, the large cavities of the pure H₂ hydrate can accommodate up to three H₂ molecules while the small ones are singly occupied at most, even at pressures as high as 500 MPa, without any double occupancy being observed. Subsequently, the binary H₂–ethylene oxide (EO) hydrate is examined. In this case, the large cavities are occupied by a single EO molecule while the small cavities can accommodate at most a single H₂ molecule. Such configuration results in a maximum H₂ content of only 0.37 wt%. The hydrogen storage capacity does not improve significantly even in case when EO is replaced by a component with smaller molecular weight.     

Vibrational spectra of deuterated methane and water molecules in structure I clathrate hydrate from ab initio MD simulation

The deuteration of the lattice molecules in clathrate hydrates is a widely used experimental technique to clearly separate the vibrational modes. However, the effect of the deuteration on the vibrational spectra and molecular motions is not fully understood. Since the guest–host coupling may change the vibrational spectra, a detailed analysis of the vibrational spectra of deuterated clathrate hydrate is significant in the understanding of the mechanism of the vibrational shift. In this study, the vibrational spectra of the deuterated methane hydrates were calculated by ab initio molecular dynamics simulation. The intramolecular vibrational frequency of the methane in D₂O lattice and deuterated methane in H₂O lattice was calculated and compared with the pure methane hydrate. The bending, rocking and overtone of the bending mode was also reported. The effect of coupling of the rattling motions of guest and host molecules on the vibrational spectra was revealed.     

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Microbial Stabilization and Kinetic Enhancement of Marine Methane Hydrates

Abstract

In clathrate hydrates, a water host lattice encages small guest molecules in cavities. Methane hydrates are the most widespread in-situ clathrate in the permafrost and continental-shelf ocean regions, constituting a significant energy resource, and prompting recent marine-hydrate gas-production trials. Despite exciting engineering advances and a few marine-mimicking laboratory studies of methane-hydrate kinetics and stabilization, from microbial perspectives, little is known about a potential microbial origin of marine hydrates, nor their possible formation kinetics or potential stabilization by microbial sources. Here, for the first time, we show that an exported, extra-cytoplasmic porin – produced by a marine methylotrophic bacterium culture – provides the basis for kinetic enhancement and stabilization of methane hydrates under conditions simulating the seabed environment. We then identify the key protein at play, and we therefore suggest microbe-based stabilization of marine hydrates is evidently a property likely to be found in many marine bacteria. Our research opens the possibility of managing marine-hydrate deposits using microbiological strategies for environmental and societal benefit.     

Molecular simulation of cage occupancy and selectivity of binary THF–H2 sII hydrate

The hydrogen capacity of the binary THF–H₂ sII hydrate is determined by the cage occupancy and by the selectivity of guest molecules. Grand canonical Monte Carlo (GCMC) simulation is used to study the cage occupancy and selectivity of guest molecules from the equilibrium configuration of the binary sII hydrate. The cage framework is regarded as a rigid body and the number of guest molecules is varied to preserve the grand canonical ensemble. The occupancy and selectivity were investigated at a temperature of 270 K for pressures ranging from 0.1 to 200 MPa. It was found that most large cages select THF as guest molecules while small cages include only hydrogen molecules. Multiple occupancy of hydrogen, up to four molecules in large cages and two molecules in small cages, was found as the pressure increases. GCMC results show that the hydrogen capacity is approximately 1.1 wt% at 200 MPa.     

Microhydration of caesium compounds: Cs,CsOH, CsI and Cs2I2 complexes with one to three H2O molecules of nuclear safety interest

The structure and properties of natural gas hydrates containing hydrocarbons, CO₂, and N₂ molecules were studied by using computational quantum chemistry methods via the density functional theory approach. All host cages involved in I, II, and H types structures where filled with hydrocarbons up to pentanes, CO₂ and N₂ molecules, depending on their size, and the structures of these host–guest systems optimized. Structural properties, vibrational spectra, and density of states were analyzed together with results from atoms-in-a-molecule and natural bond orbitals methods. The inclusion of dispersion terms in the used functional plays a vital role for obtaining reliable information, and thus, B97D functional was shown to be useful for these systems. Results showed remarkable interaction energies, not strongly affected by the type of host cage, with molecules tending to be placed at the center of the cavities when host cages and guest molecules cavities are of similar size, but with molecules approaching hexagonal faces for larger cages. Vibrational properties show remarkable features in certain regions, with shiftings rising from host-guest interactions, and useful patterns in the terahertz region rising from water surface vibrations strongly coupled with guest molecules. Likewise, calculations on crystal systems for the I and H types were carried out using a pseudopotential approach combined with Grimme’s method to take account of dispersion.

Mechanisms for Stabilising Water Clathrates

Molecular dynamics computer simulations are used to examine various possible mechanisms for stabilising water clathrates. It is found that the attractive interactions between guest and host molecules are not sufficient to explain the clathrate's stability; nor is there any significant dynamic coupling between the guest and host molecules in the clathrate. Instead, it appears that the repulsive guest-host interactions provide the mechanical stability of the clathrate, and that they are even sufficient to define its structure. Some of the thermodynamic implications of these observations are discussed.     

Electrodissociation of clathrate-like structures

In the present work, we develop molecular dynamics (MD) simulations in the NPT (isobaric–isothermic) ensemble to analyse the effect of an external electrostatic field over a cubic methane hydrate crystallite. The amplitude of the field is in the range 0.5–3.0 V/nm. For the simulations, we used the SPC/E rigid water model and a single-site model for methane at a temperature of 248 K and a pressure of 20 bar. When the external electrostatic field is applied, the water dipoles are oriented in such a way that the methane molecules can diffuse far away from the water cages, hence the clathrate dissociation takes place. This last phenomenon was observed for intensities above 1.5 V/nm. Taking the final configuration of each run as input, we develop a new set of MD simulations, and we observe that the stable clathrate is not recovered immediately when the external electrostatic field is turned off due to limitations in the simulation time.     

Effects of temperature and concentration on the structure of ethylene oxide–propylene oxide–ethylene oxide triblock copolymer (Pluronic P65) in aqueous solution: a molecular dynamics simulation study

The effects of temperature and solution concentration on the structure of triblock polymeric surfactant (ethylene oxide)₁₉(propylene oxide)₂₉(ethylene oxide)₁₉ (Pluronic P65) have been investigated by fully atomistic molecular dynamics simulations. The Flory–Huggins interaction parameter χ, hydrogen bonding and molecular mobility in the aqueous solution of P65 were investigated covering a composition range of 0.1–0.73 (water weight fraction) and a temperature range of 273–373 K. The Flory–Huggins parameters indicated that propylene oxide (PO) segments became hydrophobic with the increase in temperature, whereas ethylene oxide (EO) segments remained hydrophilic, which caused the increase in repulsion between EO and PO segments. The intermolecular hydrogen bonds in P65 solution including water–water hydrogen bonds and water–P65 hydrogen bonds increased with the increase in solution concentration and decreased with the increase in temperature. The critical micellar temperature of Pluronic P65 predicted by Flory–Huggins interaction parameter χ and hydrogen bonding was in good agreement with experimental data.     

Potential theory for gate adsorption on soft porous crystals

We demonstrate that an adsorption potential at the gate adsorption pressure of soft porous crystals (SPCs) based on the Polanyi's potential theory of adsorption shows a constancy to temperature. This was done using grand canonical Monte Carlo simulations and free energy analysis, which were carried out with a simplified stacked-layer SPC model. This finding implies that the characteristic curve obtained from an experimental gate adsorption isotherm on SPCs can be used to predict the temperature dependence of the gate-opening pressure, even though the potential theory of adsorption does not take into account the deformation of porous solids during the adsorption. We develop a modified potential theory for gate adsorption and show that the derived relation has a form that the Gibbs free energy change due to the host framework deformation per guest molecule, ? ΔG_host/N, and a correction term, C, are added to the expression of the original potential theory of adsorption. The term C is not an empirical correction factor but is the difference of intermolecular interaction potential and entropy between the bulk liquid phase at the saturated state and the adsorbed phase, originating from spatial constraint of adsorbed guest molecules in the host. By evaluating the modified expression for gate adsorption using the simulation results, we demonstrate that the constancy of the adsorption potential to temperature results from a compensation effect between three terms: guest–host interaction potential per guest molecule, ? ΔG_host/N and C, which have a temperature dependence.     

Ni‐Metalloid (B,Si, P,As, and Te) Alloys as Water Oxidation Electrocatalysts

Breakthroughs toward effective water‐splitting electrocatalysts for mass hydrogen production will necessitate material design strategies based on unexplored material chemistries. Herein, Ni‐metalloid (B, Si, P, As, Te) alloys are reported as an emergent class of highly promising electrocatalysts for the oxygen evolution reaction (OER) and insight is offered into the origin of activity enhancement on the premise of the surface electronic structure, the OER activation energy, influence of the guest metalloid elements on the lattice structure of the host metal (Ni), and surface‐oxidized metalloid oxoanions. The metalloids modify the lattice structure of Ni, causing changes in the nearest Ni–Ni interatomic distance (d_Ni–Ni). The activation energy E_a scales with d_Ni–Ni indicating an apparent dependence of the OER activity on lattice properties. During the OER, surface Ni atoms are oxidized to nickel oxyhydroxide, which is the active state of the catalyst, meanwhile, the surface metalloids are oxidized to the corresponding oxoanions that affect the interfacial electrode/electrolyte properties and hence the adsorption/desorption interaction energies of the reacting species.     

Structural properties of sH hydrate: a DFT study of anisotropy and equation of state

ABSTRACT

Structure-H (sH) hydrate is one of the canonical gas hydrates with significant potential applications and scarce characterised material properties despite the wide knowledge available on other gas hydrates. In this work we characterise some of the important physical properties of this hydrate at the atomistic level using Density Functional Theory. Two exchange-correlation functionals (revPBE and DRSLL) were used to simulate six sH hydrate systems encapsulating neohexane and different help gas molecules. The important role of dispersion forces is quantified. The density and isothermal bulk modulus of sH hydrate are higher when dispersion interactions are considered. The presence of those interactions imposes a direct relationship between the hydrate density and its bulk modulus, while their absence reveals the bulk modulus dependency on hydrogen bond density. Anisotropy is a distinguishing feature of this hydrate in distinction to nearly isotropic sI and sII hydrates. Structure-H hydrate experiences a compressional anisotropy in which the a-lattice and the c-lattice constants respond differently to applied pressure showing less compressibility along the c-axis. This compressional anisotropy was found dependant on the chemistry of help gas molecules. Taken together, these property characterisation results and analysis are a significant and novel contribution to the material physics of sH hydrates.     

Ab initio DFT study of structural and mechanical properties of methane and carbon dioxide hydrates

The structural and mechanical properties of methane and carbon dioxide hydrates were investigated using density functional theory simulations. Well-established equations of state of solids and exchange-correlation functionals were used for fitting the unit lattice total energy as a function of volume, and the full second-order elastic constants of these two gas hydrates were determined by energy–strain analyses. The polycrystalline elastic properties were also calculated from the unit lattice results. The final results for methane hydrate agree well with available experimental data and with other theoretical results. The two gas hydrates were found to be highly elastically isotropic, but they differed significantly in shear properties. The presented results for carbon dioxide hydrates are the first complete set reported so far. The results are a significant contribution to the ab initio material characterisation of gas hydrates required for ongoing fundamental studies and technological applications.     

Analysis of dissociation process for gas hydrates by molecular dynamics simulation

The dissociation processes of methane and carbon dioxide hydrates were investigated by molecular dynamics simulation. The simulations were performed with 368 water molecules and 64 gas molecules using NPT ensembles. The TraPPE (single-site) and 5-site models were adopted for methane molecules. The EPM2 (3-site) and SPC/E models were used for carbon dioxide and water molecules, respectively. The simulations were carried out at 270 K and 5.0 MPa for hydrate stabilisation. Then, temperature was increased up to 370 K. The temperature increasing rates were 0.1–20 TK/s. The gas hydrates dissociated during increasing temperature or at 370 K. The potential models of methane molecule did not much influence the dissociation process of methane hydrate. The mechanisms of dissociation process were analysed with the coordination numbers and mean square displacements. It was found that the water cages break down first, then the gas molecules escape from the water cages. The methane hydrate was more stable than the carbon dioxide hydrate at the calculated conditions.     

Stability analysis for binary sII hydrogen–promoter hydrates by molecular dynamics simulation

Molecular dynamics simulation was applied for the binary sII hydrogen–promoter hydrates to search the potential promoters to stabilise the hydrogen hydrates. The simulations were performed at 10.1 MPa. The simulation temperature was maintained at 260 K for 100 ps, and then it was increased at the rate of 0.1 TK/s. The cell volumes of the hydrates slowly increased with increasing temperature, and then the cell volumes rapidly increased. The temperature at which the cell volumes rapidly increased is identified as the simulated collapse temperature. The promoter which gives high simulated collapse temperature is judged to stabilise the hydrates. The simulated collapse temperature of the hydrate filled with cyclobutane is the highest among the promoters studied in this work.