Grand canonical Monte Carlo simulation study of hydrogen storage by Li-decorated pha-graphene

Grand canonical Monte Carlo simulation(GCMCs)is utilized for studying hydrogen storage gravimetric density by pha-graphene at different metal densities,temperatures and pressures.It is demonstrated that the optimum adsorbent location for Li atoms is the center of the seven-membered ring of pha-graphene.The binding energy of Li-decorated phagraphene is larger than the cohesive energy of Li atoms,implying that Li can be distributed on the surface of pha-graphene without forming metal clusters.We fitted the force field parameters of Li and C atoms at different positions and performed GCMCs to study the absorption capacity of H_(2).The capacity of hydrogen storage was studied by the differing density of Li decoration.The maximum hydrogen storage capacity of 4Li-decorated pha-graphene was 15.88 wt%at 77 K and100 bar.The enthalpy values of adsorption at the three densities are in the ideal range of 15 kJ·mol^(-1)-25 kJ·mol^(-1).The GCMC results at different pressures and temperatures show that with the increase in Li decorative density,the hydrogen storage gravimetric ratio of pha-graphene decreases but can reach the 2025 US Department of Energy's standard(5.5 wt%).Therefore,pha-graphene is considered to be a potential hydrogen storage material.

Adsorption of Benzene and Propylene in Zeolite ZSM-5:Grand Canonical Monte Carlo Simulations

The adsorption behavior of benzene and propylene in zeolite ZSM-5 was studied by Grand Canonical Monte Carlo（GCMC） simulations. It could be found that benzene and propylene molecules showed different adsorption behavior in the zeolite cavities. The loadings of propylene were significantly larger than those of benzene at 100 kPa. From the figures of potential energy distribution, the potential energy of benzene/zeolite was more negative than that of propylene/zeolite, so benzene could be adsorbed more stably than propylene. When the temperature was in- creased from 298 to 443 K at 100 kPa, the loading ofpropylene was reduced from 99 to 82 molecules, whereas that of benzene changed little. When benzene and propylene were adsorbed in zeolite simultaneously, the competitive adsorption of them occurred; therefore, the potential energy distribution could be changed significantly. Besides, the adsorption isotherms of benzene and propylene in ZSM-5 at 298 and 443 K were simulated. The results exhibit that the different factors influenced the molecular adsorption at various temperatures and pressures, leading to the diffe- rent rules for the adsorption of benzene and propylene molecules in the zeolite. At a low pressure, the unfavorable energy of propylene/zeolite and the ＂commensurate freezing＂ phenomenon of benzene would make the loadings of itself higher than those of propylene. When pressure was higher than 5 kPa, the adsorption of benzene in ZSM-5 would nearly reach saturation.

Diffusion and Adsorption of Benzene and Propylene in MFI,MWW and BEA Zeolites:Molecular Dynamics and Grand Canonical Monte Carlo Simulations

The diffusion and adsorption behaviors of benzene and propylene in zeolites MFI, MWW and BEA have been studied by molecular dynamics（MD） and grand canonical Monte Carlo（GCMC） simulations. The diffusion coefficients of benzene and propylene in MFI, MWW and BEA zeolites were calculated by simulating the mean-square displacements（MSD） at 298 and 600 K. Benzene and propylene showed the different adsorption rules in the channels of the three zeolites. For propylene, the molecular loadings decreased in the order： BEA（linear channel）〉BEA （tortuous channel）〉MFI（linear channel）〉MWW（12-membered rings, 12MR channel）〉MFI（tortuous channel）〉MWW （10-membered rings, 10MR channel）; for benzene, the molecular loadings decreased in the order： BEA（linear chan-nel）〉BEA（tortuous channel）〉MWW（12MR channel）〉MFI（linear channel）〉MFI（tortuous channel）〉MWW（10MR channel）. Besides, the adsorption isotherms of benzene and propylene in the three zeolites at 298 and 443 K were simulated. The results show that the different factors influenced the molecular adsorption at various temperatures and pressures, leading to the different rules for the adsorption of benzene and propylene molecules in the zeolites. At a low pressure, the unfavorable energy would make the loadings of propylene lower than those of benzene. When pressure was higher than 0.25 kPa, the adsorption of benzene in MFI would nearly reach saturation.

Sorption and permeation of gaseous molecules in amorphous and crystalline PPX C membranes: molecular dynamics and grand canonical Monte Carlo simulation studies

Amorphous and crystalline poly （chloro-p-xytylene） （PPX C） membranes are constructed by using a novel com- putational technique, that is, a combined method of NVT＋NPT-molecular dynamics （MD） and gradually reducing the size （GRS） methods. The related free volumes are defined as homology clusters. Then the sorption and the permeation of gases in PPX C polymers are studied using grand canonical Monte Carlo （GCMC） and NVT-MD methods. The results show that the crystalline PPX C membranes provide smaller free volumes for absorbing or transferring gases relative to the amorphous PPX C area. The gas sorption in PPX C membranes mainly belongs to the physical one, and H bonds can appear obviously in the amorphous area. By cluster analyzing on the mean square displacement of gases, we find that gases walk along the x axis in the crystalline area and walk randomly in the amorphous area. The calculated permeability coefficients are close to the experimental data.

The Grand Canonical Monte Carlo Simulations of Benzene and Propylene in ITQ-1 Zeolite

Grand Canonical Monte Carlo (GCMC) simulations have been performed to study the localization and adsorption behavior of benzene and propylene, in purely siliceous MWW zeolite (ITQ-1). By analyzing the locations of benzene and propylene in ITQ-1, it can be deduced that the alkylation of benzene and propylene will mainly happen in 12-MR supercages at the external surface or close to the external surface. The adsorption isotherms of benzene and propylene at 315K and 0 similar to 3.5kPa are predicted, and the results for benzene generally coincide with the trend from the experiments of a series of aromatic compounds.

Physical Experiments and Mechanism Study on the Occurrence State of Hydrogen in Clay Minerals

In recent years,natural hydrogen has been discovered in various geological environments both domestically and internationally,which has sparked a global interest among geologists and led to a surge in the exploration of hydrogen gas(Klein et al.,2019;Prinzhofer et al.,2019;Moretti and Webber,2021;Scott,2021;Bezruchko,2022).However,there is a lack of research on the occurrence state of natural hydrogen gas,which hinders a deeper understanding of its behavior in underground storage and migration.

Molecular Simulations of Adsorption and Diffusion Behaviors of Benzene Molecules in NaY Zeolite

In the article the Grand Canonical Monte Carlo （GCMC）, molecular dynamics（MD）, and kinetic Monte Carlo （KMC） simulations with particular focus on ascertaining the loading dependence of benzene diffusion in the zeolite were performed. First, a realistic representation of the structure of the sorbate-sorbent system was obtained based on GCMC simulation. The simulation clearly shows the characteristics of the adsorption sites of the benzene-NaY system, from which two kinds of preferably adsorbing sites for benzene molecules, called SⅡ and W sites, are identified. The structure thus obtained was then used as a basis for KMC and MD simulations. A compara-tive study by introducing and comparing two different mechanisms underlying jump diffusion in the zeolite of in-terest shows that the.MS diffusivity values predicted by the KMC and MD methods are fairly close to each other,leading to the conclusion that for benzene diffusion in NaY, the SⅡ→W→SⅡ jumps of benzene molecules are dominated,while the W→Wjumps do not exist in the process. These findings provide further support to our previous conclusion about the absence of the W→W jumps in the process of benzene diffusion in NaY. Finally, to relations, for predicting the self-and MS difthsivities were derived and found to be in fair agreement with the KMC and MD simulations.

Multiscale Simulation Starting at the Molecular Level for Adsorption Process Development

This article presents a multiscale simulation approach starting at the molecular level for the adsorption process development. A grand canonical Monte Carlo method is used for the prediction of adsorption isotherms of methanol on an activated carbon at the molecular level. The adsorption isotherms obtained in the linear region （or adsorption constant） are exploited as a model parameter required for the adsorption process simulation. The adsorption process model described by a set of partial differential equations （PDEs） is solved by using the conservation element and solution element method, which produces a fast and an accurate numerical solution to PDEs. The simulation results obtained from the adsorption constant estimated at the molecular level are in good agreement with the experimental results of the pulse response. The systematical multiscale simulation approach addressed in this study may be useful to accelerate the adsorption process development by reducing the number of experiments.

Theoretical study of molecular hydrogen and spiltover hydrogen storage on two-dimensional covalent-organic frameworks

Molecular hydrogen and spiltover hydrogen storages on five two-dimensional （2D） covalent-organic frameworks （COFs） （PPy-COF, TP-COF, BTP-COF, COF-18 A, and HHTP-DPB COF） are investigated using the grand canonical Monte Carlo （GCMC） simulations and the density functional theory （DFT）, respectively. The GCMC simulated results show that HHTP-DPB COF has the best performance for hydrogen storage, followed by BTP-COF, TP-COF, COF-18 A, and PPy-COE However, their adsorption amounts at room temperature are all too low to meet the uptake target set by US Department of Energy （US-DOE） and enable practical applications. The effects of pore size, surface area, and isosteric heat of hydrogen on adsorption amount are considered, which indicate that these three factors are all the important factors for determining the H2 adsorption amount. The chemisorptions of spiltover hydrogen atoms on these five COFs represented by the cluster models are investigated using the DFT method. The saturation cluster models are constructed by considering all possible adsorption sites for these cluster models. The average binding energy of a hydrogen atom and the saturation hydrogen storage density are calculated. The large average binding energy indicates that the spillover process may pro- ceed smoothly and reversibly. The saturation hydrogen storage density is much larger than the physisorption uptake of H2 molecules at 298 K and 100 bar （1 bar = 105 Pa）, and is close to or exceeds the 2010 US-DOE target of 6 wt% for hydrogen storage. This suggests that the hydrogen storage capacities of these COFs by spillover may be significantly enhanced. Thus 2D COFs studied in this paper are suitable hydrogen storage media by spillover.

Molecular Simulation of Carbon Dioxide Adsorbed in a Slit Carbon Pore

Both the grand canonical Monte Carlo and molecular dynamics simulation methods are used to investigate the adsorption and diffusion of carbon dioxide confined in a 1.86 nm slit carbon pore at 4 temperatures from subcritical (120 K) to supercritical (313 K) conditions. Layering transition, capillary condensation and adsorption hysteresis are found at 120 K. The microstructure of carbon dioxide fluid in the slit carbon pore is analyzed. The diffusion coefficients of carbon dioxide parallel to the slit wall are significantly larger than those normal to the slit wall.

Molecular simulation on carbon dioxide capture performance for carbons doped with various elements

Among the different types of CO_(2)capture technologies for post-combustion,sorption CO_(2)capture technology with carbon-based sorbents have been extensively explored with the purpose of enhancing their sorption perfor-mance by doping hetero elements due to the rapid reaction kinetics and low costs.Herein,sorption capacity and selectivity for CO_(2)and N 2 on carbon-based sorbents doped with elements such as nitrogen,sulfur,phosphorus,and boron,are evaluated and compared using the grand canonical Monte Carlo(GCMC)method,the universal force field(UFF),and transferable potentials for phase equilibria(TraPPE).The sorption capacities of N-doped porous carbons(PCs)at 50℃were 76.1%,70.7%,50.6%,and 35.7%higher than those of pure PCs,S-doped PCs,P-doped PCs,and B-doped PCs,respectively.Its sorption selectivity at 50℃was approximately 14.0,nearly twice that of pure PCs or other hetero-element-doped PCs.The N-doped PCs showed the largest sorption heat at 50℃among all the PCs,approximately 20.6 kJ·mol^(−1),which was 9.7%−25.5%higher than that of the pure PCs under post-combustion conditions.Additionally,with the product purity of 41.7 vol.%−75.9 vol.%for vacuum pressure swing sorption,and 53.4 vol.%−83.6 vol.%for temperature swing sorption,the latter is more suitable for post-combustion conditions than pressure-swing sorption.

Effects of coal molecular structure and pore morphology on methane adsorption and accumulation mechanism

The adsorption,diffusion,and aggregation of methane from coal are often studied based on slit or carbon nanotube models and isothermal adsorption and thermodynamics theories.However,the pore morphology of the slit model involves a single slit,and the carbon nanotube model does not consider the molecular structure of coal.The difference of the adsorption capacity of coal to methane was determined without considering the external environmental conditions by the molecular structure and pore morphology of coal.The study of methane adsorption by coal under single condition cannot reveal its mechanism.In view of this,elemental analysis,FTIR spectrum,XPS electron energy spectrum,13C NMR,and isothermal adsorption tests were conducted on the semi-anthracite of Changping mine and the anthracite of Sihe Mine in Shanxi Province,China.The grand canonical Monte Carlo(GCMC)and molecular dynamics simulation method was used to establish the coal molecular structure model.By comparing the results with the experimental test results,the accuracy and practicability of the molecular structure model are confirmed.Based on the adsorption potential energy theory and aggregation model,the adsorption force of methane on aromatic ring structure,pyrrole nitrogen structure,aliphatic structure,and oxygen-containing functional group was calculated.The relationship between pore morphology,methane aggregation morphology,and coal molecular structure was revealed.The results show that the adsorption force of coal molecular structure on methane is as follows:aromatic ring structure(1.96 kcal/mol)>pyridine nitrogen(1.41 kcal/mol)>pyrrorole nitrogen(1.05 kcal/mol)>aliphatic structure(0.29 kcal/mol)>oxygen-containing functional group(0.20 kcal/mol).In the long and narrow regular pores of semi-anthracite and anthracite,methane aggregates in clusters at turns and aperture changes,and the adsorption and aggregation positions are mainly determined by the aromatic ring structure,the positions of pyrrole nitrogen and pyridine nitrogen.The degree of aggregation is controlled by the interaction energy and pore morphology.The results pertaining to coal molecular structure and pore morphology on methane adsorption and aggregation location and degree are conducive to the evaluation of the adsorption mechanism of methane in coal.

Hydrogen storage in Li-doped fullerene-intercalated hexagonal boron nitrogen layers

New materials for hydrogen storage of Li-doped fullerene （C20, C28, C36, C50, C60, C70）-intercalated hexagonal boron nitrogen （h-BN） frameworks were designed by using density functional theory （DFT） calculations. First-principles molecular dynamics （MD） simulations showed that the struc- tures of the Cn-BN （n = 20, 28, 36, 50, 60, and 70） frameworks were stable at room temperature. The interlayer distance of the h-BN layers was expanded to 9.96-13.59 A° by the intercalated fullerenes. The hydrogen storage capacities of these three-dimensional （3D） frameworks were studied using grand canonical Monte Carlo （GCMC） simulations. The GCMC results revealed that at 77 K and 100 bar （10 MPa）, the C50-BN framework exhibited the highest gravimetric hydrogen uptake of 6.86 wt% and volumetric hydrogen uptake of 58.01 g/L. Thus, the hydrogen uptake of the Li-doped Cn-intercalated h-BN frameworks was nearly double that of the non-doped framework at room temperature. Furthermore, the isosteric heats of adsorption were in the range of 10-21 kJ/mol, values that are suitable for adsorbing/desorbing the hydrogen molecules at room temperature. At 193 K （-80 ℃） and 100 bar for the Li-doped C50-BN framework, the gravimetric and volumetric uptakes of H2 reached 3.72 wt% and 30.08 g/L, respectively.