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.     

Molecular dynamics simulation on the decomposition of type SII hydrogen hydrate and the performance of tetrahydrofuran as a stabiliser

Molecular dynamics simulation is used to study the decomposition and stability of SII hydrogen and hydrogen/tetrahydrofuran (THF) hydrates at 150 K, 220 K and 100 bar. The modelling of the microscopic decomposition process of hydrogen hydrate indicates that the decomposition of hydrogen hydrate is led by the diffusive behaviour of H₂ molecules. The hydrogen/THF hydrate presents higher stability, by comparing the distributions of the tetrahedral angle of H₂O molecules, radial distribution functions of H₂O molecules and mean square displacements or diffusion coefficients of H₂O and H₂ molecules in hydrogen hydrate with those in hydrogen/THF hydrate. It is also found that the resistance of the diffusion behaviour of H₂O and H₂ molecules can be enhanced by encaging THF molecules in the (5¹²6⁴) cavities. Additionally, the motion of THF molecules is restricted due to its high interaction energy barrier. Accordingly, THF, as a stabiliser, is helpful in increasing the stability of hydrogen hydrate.     

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.     

The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation

In this work, we present temperature dependence of lattice parameter and normalised lattice parameter in the atmospheric pressure and 120 bar and also pressure dependence of unit cell volume and normalised unit cell volume at 150 and 250?K for variety guests with different size, polarity and guest–host hydrogen bonding capability such as trimethylene oxide (TMO), ethylene oxide (EO), formaldehyde (FA), cyclobutane (CB), cyclopropane (CP) and ethane (Et) in the large cages with CH₄ in small cages of sI clathrate hydrates by molecular dynamics simulations. The obtained values of lattice parameters for the guest species are compatible with the experimental values. These clathrate hydrates are simulated with TIP4P/ice four-site water potential. Herein, isobaric thermal expansivity and isothermal compressibility are calculated at a temperature range of 50–250?K and a wide pressure range. These structural properties have been compared for guests which they are isoelectronic and have similar masses but with different size and polarity. We use molecular dynamics simulations to relate microscopic guest properties, like guest–host hydrogen bonding to macroscopic sI clathrate hydrate properties. The temperature dependence of thermodynamic properties such as constant-volume and constant-pressure heat capacity is presented in the atmospheric pressure for these guest species.     

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.     

Open carbon frameworks - a search for optimal geometry for hydrogen storage

Properties of a new class of hypothetical high-surface-area porous carbons (open carbon frameworks) have been discussed. The limits of hydrogen adsorption in these carbon porous structures have been analyzed in terms of competition between increasing surface accessible for adsorption and the lowering energy of adsorption. From an analysis of an analytical model and simulations of adsorption the physical limits of hydrogen adsorption have been defined: (i) higher storage capacities in slit-shaped pores can be obtained by fragmentation/truncation of graphene sheets into nano-metric elements which creates surface areas in excess of 2600 m²/g, the surface area for infinite graphene sheets; (ii) the positive influence of increasing surface area is compensated by the decreasing energy of adsorption in the carbon scaffolds of nano-metric sizes; (iii) for open carbon frameworks (OCF) built from coronene and benzene molecules with surface areas 6500 m² g^-1, we find an impressive excess adsorption of 75–110 g?H₂/kg C at 77 K, and high storage capacity of 110–150 g?H₂/kg C at 77 K and 100 bar; (iv) the new OCF, if synthesized and optimized, could lead to required hydrogen storage capacity for mobile applications.     

Molecular dynamics simulation of replacement of methane hydrate with carbon dioxide

Molecular dynamics simulation was performed to analyse the phenomena of replacement of methane hydrate with carbon dioxide (CO₂) at 270 K and 5.0 MPa for 5300 ps. The methane hydrate phase was constructed with 16 unit cells of hydrate. Every cage in the hydrate was occupied by one methane molecule. The methane hydrate phase was sandwiched between two CO₂ phases. During the simulation the hydrate partially melted and liquid water phase appeared, and CO₂ dissolved in the liquid water phase. The replacements were observed three times at the hydrate–liquid water interface during the simulation. In the first case, the replacement occurred at a S-cage without changing the structure. In the second case, an M-cage of methane hydrate partially collapsed, and methane and CO₂ molecules exchanged. After the exchange, the cage occupied by CO₂ remained in the M-cage structure. In the third case, a S-cage of methane hydrate partially collapsed, and methane and CO₂ molecules exchanged. After the exchange, the cage occupied by CO₂ changed to an M-cage-like structure.     

Catalytic Influence of Various Cerium Precursors on the Hydrogen Sorption Properties of NaAlH4

Sodium alanate (NaAlH₄) is one of the metal complex hydrides most often investigated for use as a hydrogen‐storage material. Doped with transition or rare earth metal compounds, NaAlH₄ can absorb and release hydrogen in low and medium temperature ranges with excellent reversibility and cycling stability. The properties of NaAlH₄ doped with CeCl₃ differ from materials with other dopants, with faster sorption kinetics and a more stable capacity. In this paper, various precursors of Ce are applied to investigate their catalytic effects on the sorption performance of this material. The re‐hydrogenation is found to be completed in approximately 10 min. Although all the Ce precursors investigated in this work result in reversible hydrogen storage materials, desorption kinetics are enhanced upon formation of cerium aluminide (CeAl₄) in the composites. While the use of CeAl₄ instead of CeCl₃ can increase the hydrogen capacity by bypassing the formation of the ineffective NaCl, the highest capacity of 4.9 wt%—close to the theoretical value—is obtained from NaAlH₄ doped directly with metallic cerium. Furthermore, dehydriding under back pressures is also investigated to evaluate the H₂ desorption rates under practical conditions. At 3 bar H₂ pressure, the second desorption step of NaAlH₄ is fully suppressed at 150 °C and only 2.5 wt% H was released, whereas at 180 °C the capacity is not much affected, which is interesting for combination in a system with a high‐temperature PEM fuel cell.     

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 investigations on planar (MgO)n clusters (n = 1–5) and their hydrogen adsorption behaviour

Ab initio calculations (B3LYP and PBE-D3) of the structures, stabilities, vibrational, electronic and hydrogen adsorption behaviour of (MgO)_n clusters are performed using 6-311+ + G(d,p) basis set. The planar (MgO)_n clusters are found to be global minima for n ≤ 3 and local minima for n = 4 and 5. In addition, we have also analysed global minimum structures of (MgO)₄ and (MgO)₅. The binding energies suggest that their stabilities increase successively. Vibrational frequencies and IR intensities further support the enhanced stability with an increase in the size of (MgO)_n clusters. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) surfaces are used to explain and discuss the electronic properties. Finally, we have demonstrated hydrogen storage capacity of (MgO)_n clusters, considering hydrogen adsorption on planar as well as global minimum (MgO)₄ and (MgO)₅ clusters. We have noticed that four and five H₂ molecules can be easily adsorbed by (MgO)₄ and (MgO)₅ clusters having adsorption energy of 0.13–0.14 eV with mass ratio of 4.76%. Thus, the present study is expected to motivate further the applications of small clusters for efficient hydrogen energy storage.     

Influences of lithium doping and fullerene impregnation on hydrogen storage in metal organic frameworks

Using the grand canonical ensemble Monte Carlo method, two similar metal organic frameworks (isoreticular MOFs [IRMOF]-12 and -14) and their modified structures by doping lithium (Li) atoms above the organic units and/or impregnating with fullerenes in their cavities have been employed to investigate the capacities of H₂ storage. Our simulations show that the H₂ uptakes of Li-C₆₀@Li-IRMOF-12 and Li-C₆₀@Li-IRMOF-14 achieve the U.S. Department of Energy targets before 2017 both in gravimetric density and in volumetric density at 243 K and 100 bar. Combining the results of IRMOF-10-based structures, we further study the relationships between the H₂ uptakes and the physical properties of the materials to identify the influence factors on the H₂ storage at room temperature.     

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.

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.     

Kinetic modeling of light limitation and sulfur deprivation effects in the induction of hydrogen production with Chlamydomonas reinhardtii. Part II: Definition of model‐based protocols and experimental validation

Photosynthetic hydrogen production under light by the green microalga Chlamydomonas reinhardtii was investigated in a torus‐shaped PBR in sulfur‐deprived conditions. Culture conditions, represented by the dry biomass concentration of the inoculum, sulfate concentration, and incident photon flux density (PFD), were optimized based on a previously published model (Fouchard et al., 2009. Biotechnol Bioeng 102:232–245). This allowed a strictly autotrophic production, whereas the sulfur‐deprived protocol is usually applied in photoheterotrophic conditions. Experimental results combined with additional information from kinetic simulations emphasize effects of sulfur deprivation and light attenuation in the PBR in inducing anoxia and hydrogen production. A broad range of PFD was tested (up to 500 µmol photons m⁻² s⁻¹). Maximum hydrogen productivities were 1.0 ± 0.2 mL H₂/h/L (or 25 ± 5 mL H₂/m² h) and 3.1 mL ± 0.4 H₂/h L (or 77.5 ± 10 mL H₂/m² h), at 110 and 500 µmol photons m⁻² s⁻¹, respectively. These values approached a maximum specific productivity of approximately 1.9 mL ± 0.4 H₂/h/g of biomass dry weight, clearly indicative of a limitation in cell capacity to produce hydrogen. The efficiency of the process and further optimizations are discussed. Biotechnol. Bioeng. 2011;108: 2288–2299. © 2011 Wiley Periodicals, Inc.     

Prediction of thermodynamically reversible hydrogen storage reactions utilizing Ca–M(M = Li, Na, K)–B–H systems: a first-principles study

Calcium borohydride is a potential candidate for onboard hydrogen storage because it has a high gravimetric capacity (11.5 wt.%) and a high volumetric hydrogen content (～130 kg m^?3). Unfortunately, calcium borohydride suffers from the drawback of having very strongly bound hydrogen. In this study, Ca(BH₄)₂ was predicted to form a destabilized system when it was mixed with LiBH₄, NaBH₄, or KBH₄. The release of hydrogen from Ca(BH₄)₂ was predicted to proceed via two competing reaction pathways (leading to CaB₆ and CaH₂ or CaB₁₂H₁₂ and CaH₂) that were found to have almost equal free energies. Using a set of recently developed theoretical methods derived from first principles, we predicted five new hydrogen storage reactions that are among the most attractive of those presently known. These combine high gravimetric densities (>6.0 wt.% H₂) with have low enthalpies [approximately 35 kJ/(mol^?1 H₂)] and are thermodynamically reversible at low pressure within the target window for onboard storage that is actively being considered for hydrogen storage applications. Thus, the first-principles theoretical design of new materials for energy storage in future research appears to be possible.

Tuning hydrogen storage of carbon nanotubes by mechanical bending: theoretical study

The effects of mechanical bending on tuning the hydrogen storage of titanium functionalised (4,0) carbon nanotube have been assessed using density functional theory calculations with reference to the ultimate targets of the US Department of Energy (DOE). The assessment has been carried out in terms of physisorption, gravimetric capacity, projected densities of states, statistical thermodynamic stability and reaction kinetics. The Ti atom binds at the hollow site of the hexagonal ring. The average adsorption energies (?0.54 eV) per hydrogen molecule meet the DOE target for physisorption (?0.20 to ?0.60 eV). The curvature attributed to the bending angle has no effect on the average adsorption energies per H₂ molecule. With no metal clustering, the system gravimetric capacities are expected to be as large as 9.0 wt%. The reactions of the deformed (bent) carbon nanotube have higher probabilities of occurring than those of the un-deformed carbon nanotube. The Gibbs free energies, enthalpies and entropies meet the ultimate targets of the DOE for all temperatures and pressures. The closest reactions to zero free energy occur at (378.15 K/2.961 atm.) and reverse at (340 and 360 K/1 atm.). The translational component is found to exact a dominant effect on the total entropy change with temperature. Favourable kinetics of the reactions at the temperatures targeted by DOE are reported regardless of the applied pressure. The more preferable thermodynamic properties assigned to the bending nanotube imply that hydrogen storage can be improved compared to the nonbending nanotube.     

Molecular dynamics simulation of diffusion of hydrogen in binary hydrogen–tetrahydrofuran hydrate

The binary structure II hydrogen–tetrahydrofuran (THF) hydrate was studied with molecular dynamics simulation. The simulations were carried out at 300, 310 K and 10.1 MPa, and with various contents of hydrogen and THF. The migrations of hydrogen molecules from cage to cage were observed. The migration process of hydrogen was also analysed, and the diffusion coefficients of hydrogen in the hydrate were calculated. The calculated diffusion coefficients qualitatively agreed with the experimental data. Double and quintet occupancies of hydrogen molecules were observed in the small and large cages, respectively, without changing the hydrate structure.     

Treatment of Diesel-Contaminated Soil by Fenton and Persulfate Oxidation with Zero-Valent Iron

The objective of this research is to investigate Fenton and persulfate oxidation with zero-valent iron [Fe(0)] as a batch type ex-situ remediation technology for the treatment of diesel-contaminated soil. Results from batch experiments indicate that Fe(0) is a better catalyst for H₂O₂ and persulfate than Fe²⁺ for the enhancement of Fenton and persulfate oxidation in a batch system. Maximum removal was obtained after 12 h when 1 and 2 g of Fe(0) were added to hydrogen peroxide (250 mg/L) and persulfate (250 mg/L), respectively, in a soil-water system. As the amounts of Fe(0) and persulfate were increased three times at the optimal ratio, the removal of total petroleum hydrocarbon (TPH) was enhanced accordingly. More than 90% of the TPH was removed in 3 h, and the treated soil met the Korean regulation level (500 mg/kg) for TPH. Increased amounts of Fe(0) and hydrogen peroxide (up to 10 g and 1250 mg/L, respectively) also significantly enhanced degradation under the optimal conditions. The results of our study suggest that Fe(0)-assisted Fenton and persulfate oxidation in a batch reactor may be an alternative option to treat diesel-contaminated soil.     

Sufficient intake of high amylose cornstarch maintains high colonic hydrogen production for 24 h in rats

Colonic hydrogen (H₂) can suppress oxidative stress and damage in the body. We examined the minimum requirement of high amylose cornstarch (HAS) to maintain high colonic H₂ production for 24 h. Ileorectostomized and sham-operated rats were fed a control diet supplemented with or without 20% HAS for 7 days. Colonic starch utilization was determined. Next, rats were fed the control diet with or without 10% or 20% HAS for 14 or 28 days, respectively. Breath and flatus H₂ excretion for 24 h was measured. 1.04 g of resistant fraction in HAS was utilized for 24 h by colonic bacteria. High H₂ excretion was not maintained for 24 h in rats fed the 10% HAS diet, from which only 0.89 g of resistant starch was estimated to be delivered. High colonic H₂ production for 24 h would be maintained by delivering more HAS to the large intestine than is utilized.