Storage of hydrogen at 303 K in graphite slitlike pores from grand canonical Monte Carlo simulation

Grand canonical Monte Carlo (GCMC) simulations were used for the modeling of the hydrogen adsorption in idealized graphite slitlike pores. In all simulations, quantum effects were included through the Feynman and Hibbs second-order effective potential. The simulated surface excess isotherms of hydrogen were used for the determination of the total hydrogen storage, density of hydrogen in graphite slitlike pores, distribution of pore sizes and volumes, enthalpy of adsorption per mole, total surface area, total pore volume, and average pore size of pitch-based activated carbon fibers. Combining experimental results with simulations reveals that the density of hydrogen in graphite slitlike pores at 303 K does not exceed 0.014 g/cm(3), that is, 21% of the liquid-hydrogen density at the triple point. The optimal pore size for the storage of hydrogen at 303 K in the considered pore geometry depends on the pressure of storage. For lower storage pressures, p < 30MPa, the optimal pore width is equal to a 2.2 collision diameter of hydrogen (i.e., 0.65 nm), whereas, for p congruent with 50MPa, the pore width is equal to an approximately 7.2 collision diameter of hydrogen (i.e., 2.13 nm). For the wider pores, that is, the pore width exceeds a 7.2 collision diameter of hydrogen, the surface excess of hydrogen adsorption is constant. The importance of quantum effects is recognized in narrow graphite slitlike pores in the whole range of the hydrogen pressure as well as in wider ones at high pressures of bulk hydrogen. The enthalpies of adsorption per mole for the considered carbonaceous materials are practically constant with hydrogen loading and vary within the narrow range q(st) congruent with 7.28-7.85 kJ/mol. Our systematic study of hydrogen adsorption at 303 K in graphite slitlike pores gives deep insight into the timely problem of hydrogen storage as the most promising source of clean energy. The calculated maximum storage of hydrogen is equal to approximately 1.4 wt %, which is far from the United States Department of Energy (DOE) target (i.e., 6.5 wt %), thus concluding that the total storage amount of hydrogen obtained at 303 K in graphite slitlike pores of carbon fibers is not sufficient yet.     

State of hydrogen in idealized carbon slitlike nanopores at 77 K

The purpose of this letter is to clarify recent findings and answer to the question: "What is the state of hydrogen in carbon slitlike pores at 77 K?" For this purpose, we determined the volumetric density of hydrogen in idealized carbon pores of molecular dimension at 77 K and pressure up to 1 MPa. We used quantum corrected grand canonical Monte Carlo simulation. We recognized the highest volumetric density of confined hydrogen (around 71% of hydrogen liquid at boiling point) for effective pore width 5.6 angstroms (H* = 3.04) in the considered pressure range. Our computational results are in agreement with the calculations performed by Wang and Johnson and Rzepka et al. In contrast, we did not observe the high volumetric density of hydrogen in slitlike carbon pores exceeding the density of hydrogen liquid at the boiling point as was reported by Jagiello and Thommes. Moreover, we obtained qualitative agreement between the simulation results and some experimental findings reported by Nijkamp.     

Thermodynamics of hydrogen adsorption in slit-like carbon nanopores at 77 K. Classical versus path-integral Monte Carlo simulations

Hydrogen in slit-like carbon nanopores at 77 K represents a quantum fluid in strong confinement. We have used path-integral grand canonical Monte Carlo and classical grand canonical Monte Carlo simulations for the investigation of the "quantumness" of hydrogen at 77 K adsorbed in slit-like carbon nanopores up to 1 MPa. We find that classical simulations overpredict the hydrogen uptake in carbon nanopores due to neglect of the quantum delocalization. Such disagreement of both simulation methods depends on the slit-like carbon pore size. However, the differences between the final uptakes of hydrogen computed from both classical and quantum simulations are not large due to a similar effective size of quantum/classical hydrogen molecules in carbon nanospaces. For both types of molecular simulations, the volumetric density of stored energy in optimal carbon nanopores exceeds 6.4 MJ dm(-3) (i.e., 45 kg m(-3); Department of Energy target for 2010). In contrast to the hydrogen adsorption isotherms, we found a large reduction of isosteric enthalpy of adsorption computed from the quantum Feynman's path-integral simulations in comparison to the classical values at 77 K and pressures up to 1 MPa. Depression of the quantum isosteric enthalpy of adsorption depends on the slit-like carbon pore size. For the narrow pores (pore width H in [0.59-0.7] nm), the reduction of the quantum isosteric enthalpy of adsorption at zero coverage is around 50% in comparison to the classical one. We observed new phenomena called, by us, the quantum confinement-inducing polymer shrinking. In carbon nanospaces, the quantum cyclic polymers shrink, in comparison to its bulk-phase counterpart, due to a strong confinement effect. At considered storage conditions, this complex phenomenon depends on the size of the slit-like carbon nanopore and the density of hydrogen volumetric energy. For the smallest nanopores and a low density of hydrogen volumetric energy, the reduction of the polymer effective size is the highest, whereas an increase of the pore size and the density of hydrogen volumetric energy causes the polymer swelling up to a value slightly below the one computed from the bulk phase. Quantum confinement-inducing polymer shrinking is of great importance for realizing the potential of quantum molecular sieves.     

Calculation of the characteristics of ethanol and methanol adsorbed in a pore of activated carbon by the density functional method

Density functional theory was used to perform quantum-chemical calculations of changes in the energy and structure characteristics of methanol and ethanol molecules caused by their adsorption in model slitlike pores of activated carbon. The conclusion was made that changes in these characteristics (bond lengths, angles, charges on atoms, and harmonic vibration frequencies) is additional evidence of the validity of the Tolmachev thermodynamic model, in which adsorption is treated as a quasi-chemical reaction of the addition of adsorbate molecules to adsorption centers of an adsorbent. It was shown that the arrangement of alcohol molecules, when the C-O and C-C bonds were approximately parallel to pore walls and the hydrogen atom was directed toward a nearer pore wall, was most favorable energetically. Two adsorbed alcohol molecules are also arranged parallel to pore walls and form a hydrogen bond.     

Nanoscale tubular vessels for storage of methane at ambient temperatures

Novel carbon nanostructures can serve as effective storage media for methane, a source of "clean energy" for the future. We have used Grand Canonical Monte Carlo Simulation for the modeling of methane storage at 293 K and pressures up to 80 MPa in idealized bundles of (10,10) armchair-type single-walled carbon nanotubes and wormlike carbon pores. We have found that these carbon nanomaterials can be treated as the world's smallest high-capacity methane storage vessels. Our simulation results indicate that such novel carbon nanostructures can reach a high volumetric energy storage, exceeding the US FreedomCAR Partnership target of 2010 (5.4 MJ dm(-3)), at low to moderate pressures ranging from 1 to 7 MPa at 293 K. On the contrary, in the absence of these nanomaterials, methane needs to be compressed to approximately 13 MPa at 293 K to achieve the same target. The light carbon membranes composed of bundles of single-walled carbon nanotubes or wormlike pores efficiently physisorb methane at low to moderate pressures at 293 K, which we believe should be particularly important for automobiles and stationary devices. However, above 15-20 MPa at 293 K, all investigated samples of novel carbon nanomaterials are not as effective when compared with compression alone since the stored volumetric energy and power saturate at values below those of the bulk, compressed fluid.     

Tailoring of nanoscale porosity in carbide-derived carbons for hydrogen storage

The poor performance of hydrogen storage materials continues to hinder development of fuel cell-powered automobiles. Nanoscale carbons, in particular (activated carbon, exfoliated graphite, fullerenes, nanotubes, nanofibers, and nanohorns), have not fulfilled their initial promise. Here we show that carbon materials can be rationally designed for H2 storage. Carbide-derived carbons (CDC), a largely unknown class of porous carbons, are produced by high-temperature chlorination of carbides. Metals and metalloids are removed as chlorides, leaving behind a collapsed noncrystalline carbon with up to 80% open pore volume. The detailed nature of the porosity-average size and size distribution, shape, and total specific surface area (SSA)-can be tuned with high sensitivity by selection of precursor carbide (composition, lattice type) and chlorination temperature. The optimum temperature is bounded from below by thermodynamics and kinetics of chlorination reactions and from above by graphitization, which decreases SSA and introduces H2-sorbing surfaces with binding energies too low to be useful. Intuitively, pores of different size and shape should not contribute equally to hydrogen storage. By correlating pore properties with 77 K H2 isotherms from a wide variety of CDCs, we experimentally confirm that gravimetric hydrogen storage capacity normalized to total pore volume is optimized in materials with primarily micropores ( approximately 1 nm) rather than mesopores. Thus, in agreement with theoretical predictions, a narrow size distribution of small pores is desirable for storing hydrogen, while large pores merely degrade the volumetric storage capacity.     

Adsorption and selectivity of carbon dioxide with methane and nitrogen in slit-shaped carbonaceous micropores: Simulation and experiment

We have used the grand canonical Monte Carlo method to study the adsorption and selectivity of mixtures of carbon dioxide with methane and nitrogen at high (i.e., ambient) temperatures in model slit pores with graphitic surfaces. Experimental data, including new high pressure measurements for carbon dioxide and methane on a non-porous graphitic standard, were used to test the potential models. The mixture simulations predict that carbon dioxide is preferentially adsorbed in both systems. The results are discussed in terms of competing energetic and entropic effects and the underlying molecular mechanisms.     

On interfacial tension calculation from the test-area methodology in the grand canonical ensemble

We propose the extension of the test-area methodology, originally proposed to evaluate the surface tension of planar fluid-fluid interfaces along a computer simulation in the canonical ensemble, to deal with the solid-fluid interfacial tension of systems adsorbed on slitlike pores using the grand canonical ensemble. In order to check the adequacy of the proposed extension, we apply the method for determining the density profiles and interfacial tension of spherical molecules adsorbed in slitlike pore with different pore sizes and solid-fluid dispersive energy parameters along the same simulation. We also calculate the solid-fluid interfacial tension using the original test-area method in the canonical ensemble. Agreement between the results obtained from both methods indicate that both methods are fully equivalent. The advantage of the new methodology is that allows to calculate simultaneously the density profiles and the amount of molecules adsorbed onto a slitlike pore, as well as the solid-fluid interfacial tension. This ensures that the chemical potential at which all properties are evaluated during the simulation is exactly the same since simulations can be performed in the grand canonical ensemble, mimicking the conditions at which the adsorption experiments are most usually carried out in the laboratory.     

Computer simulation of the adsorption of methane, nitrogen, and mixtures thereof in pores of a lamellar carbon adsorbent

The Monte Carlo method in conjunction with the grand canonical ensemble was used to calculate the isotherms of adsorption of methane, nitrogen, and mixtures thereof in 1.34 × 3.02-nm rectangular-cross-section and 2.35 × 2.35-nm square pores in a lamellar carbon adsorbent. The phase diagrams of adsorbed methane were plotted, and the characteristics of the phases of the adsorbate were described. Modeling the adsorption of the binary mixture demonstrated that the square carbon pore is more selective with respect to methane.     

内置多孔介质卷式反应器富燃制氢实验研究

设计了一个新型内置多孔介质卷式制氢反应器,进行了天然气/空气预混气在当量比1.25~2.50,总流量60 L/min~120 L/min条件下的富燃实验,研究其自热重整制氢的燃烧特性。结果表明,富燃情况下可实现自维持燃烧,自稳定响应时间较快;反应器Swiss-roll结构的有效预热与中心多孔介质的蓄热,可提高反应温度至1 600 K以上,实现超绝热富燃制氢;固定流量条件下,甲烷的转化效率与能量利用率随着当量比的变化而改变。实验中,氢气生成率与能量利用率分别为30%~57%、50%~84%。流量为120 L/min,当量比为2.0时,尾气中H2浓度达22%,氢气生成率为57%,能量利用率为68%。     

Hydrogen storage in nanoporous carbon materials: myth and facts

We used Grand canonical Monte Carlo simulation to model the hydrogen storage in the primitive, gyroid, diamond, and quasi-periodic icosahedral nanoporous carbon materials and in carbon nanotubes. We found that none of the investigated nanoporous carbon materials satisfy the US Department of Energy goal of volumetric density and mass storage for automotive application (6 wt% and 45 kg H(2) m(-3)) at considered storage condition. Our calculations indicate that quasi-periodic icosahedral nanoporous carbon material can reach the 6 wt% at 3.8 MPa and 77 K, but the volumetric density does not exceed 24 kg H(2) m(-3). The bundle of single-walled carbon nanotubes can store only up to 4.5 wt%, but with high volumetric density of 42 kg H(2) m(-3). All investigated nanoporous carbon materials are not effective against compression above 20 MPa at 77 K because the adsorbed density approaches the density of the bulk fluid. It follows from this work that geometry of carbon surfaces can enhance the storage capacity only to a limited extent. Only a combination of the most effective structure with appropriate additives (metals) can provide an efficient storage medium for hydrogen in the quest for a source of "clean" energy.     

Molecular simulation of binary mixture adsorption in buckytubes and MCM-41

MCM-41 and buckytubes are novel porous materials with controllable pore sizes and narrow pore size distributions. Buckytubes are carbon tubes with internal diameters in the range 1–5 urn. The structure of each tube is thought to be similar to one or more graphite sheets rolled up in a helical manner. MCM-41 is one member of a new family of highly uniform mesoporous silicate materials produced by Mobil, whose pore size can be accurately controlled in the range 1.5–10 nm. We present grand canonical Monte Carlo (GCMC) simulations of single fluid and binary mixture adsorption in a model buckytube, and nonlocal density functional theory (DFT) calculations of trace pollutant separation in a range of buckytubes and MCM-41 pores. Three adsorbed fluids are considered; methane, nitrogen and propane. The GCMC studies show that the more strongly adsorbed pure fluid is adsorbed preferentially from an equimolar binary mixture. Ideal adsorbed solution theory (IAST) is shown to give good qualitative agreement with GCMC when predicting binary mixture separations. The DFT results demonstrate the very large increases in trace pollutant separation that can be achieved by tuning the pore size, structure, temperature and pressure of the MCM-41 and buckytube adsorbent systems to their optimal values.     

Effects of fluorination modification on pore size controlled electrospun activated carbon fibers for high capacity methane storage

Electrospun carbon fibers were prepared as a methane storage medium. Chemical activation was carried out using potassium carbonate to develop the pore structure, which can provide sites for the uptake of methane, and then fluorination surface modification was conducted to enhance the capacity of storage. Chemical activation provided a highly microporous structure, which is beneficial for methane storage, with a high specific surface area greater than 2500 m²/g. The pore size distribution showed that the prepared samples have pore sizes in the range of 0.7–1.6 nm. The effect of fluorination surface modification was also investigated. The functional groups, which were confirmed by XPS analysis, played an important role in guiding methane gas into the carbon silt pores via the attractive force felt by the electrons in the methane molecules due to the high electronegativity of fluorine. Eventually, the methane uptake increased up to 18.1 wt.% by the synergetic effects of the highly developed micropore structure and the guiding of methane to carbon pores by fluorine.     

巨正则系综Monte Carlo模拟方法确定活性炭的微孔尺寸

根据299K下甲烷在活性炭中的吸附实验数据,通过调节狭缝微孔的孔宽参数,利用巨正则系综MonteCarlo(GCEMC)方法得到不同孔宽下流体的微观结构以及吸附等温线.比较并拟合模拟结果和实验数据,确定了活性炭微孔的平均孔宽,为下一步求解微孔尺寸分布以及为预测吸附剂在不同温度下吸附不同吸附质分子时的吸附性能提供了基础与指导.模拟中,甲烷分子采用单点Lennard-Jones球型分子模型,活性炭用狭缝孔来近似表征,流体分子与单个狭缝墙的相互作用采用著名的Steele的10－4-3势能模型.模拟表明,此方法为考察介孔材料的微孔分布以及微孔平均孔宽提供了新的思路.     

Simulating of alcohol adsorption in slitlike micropores of active carbon by the molecular dynamics method

Benzene, methanol, ethanol, and 1,3-propanediol adsorption in model pores of active carbon is studied by the molecular dynamics method using the OPLSAA universal force field with allowance for the microheterogeneous structure of the carbon. It is shown that, in a slitlike pore with a width and a diameter of 0.7 and 3 nm, respectively, adsorbate molecules are, on average, located in parallel with pore walls and have a constant deviation from this position. When few molecules are contained in a limited pore space, they are located independently of one another in two different layers; however, as their concentration increases, initially, aggregates bonded (in the case of alcohols) through hydrogen bonds and, subsequently, dense layers are formed that strongly affect one another.     

Gas–solid interactions in the hydrogen/single-walled carbon nanotube system

Gas–solid interactions between hydrogen and single-walled carbon nanotubes (SWNTs) were investigated using highly purified SWNTs. The activation energy of hydrogen desorption, measured to be 0.21 eV, indicates that hydrogen is physisorbed in the pores and that inter-tube pores have an adsorption potential of about −0.21 eV, which induces hydrogen physisorption at ambient temperature. The total amount of adsorbed hydrogen, about 0.3 wt% at 9 MPa, shows that 38% of the inter-tube sites are occupied. These findings are interpreted in terms of the chemical potential of hydrogen and the adsorption potential of the inter-tube pores.     

Adsorption of Fluids in Pores Formed between Two Hard Cylinders

We study adsorption in pores with curved hard walls that are made of two uniaxial cylinders by using a density functional approach. Two cases are considered: adsorption of hard spheres and adsorption of a Lennard-Jones fluid. In the case of hard spheres, we perform a comparison with the results of grand canonical ensemble Monte Carlo data. This comparison indicates that the applied approach is capable of reproducing the fluid structure quite satisfactorily. For hard spheres, we also make a comparison of the total adsorption effect (expressed as the average density of a confined fluid) inside pores with curved walls with that evaluated for a slitlike pore. We have found that the differences between adsorption in pores with curved walls and in slits with the same wall-to-wall distance are quite low. The calculations for the Lennard-Jones fluid have been concerned with the investigation of the capillary evaporation and with the evaluation of phase diagrams for different pores, including slitlike pores. We have found that the curvature of the pore walls shifts the transition toward lower values of the chemical potential and increases slightly the value of the critical temperature in comparison with the values obtained for a slitlike pore. Copyright 2000 Academic Press.     

Density functional description of adsorption in slitlike pores modified with chain molecules: a simple model for pillaredlike materials

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in slitlike pores modified by chain molecules. Specifically, the chains are bonded by their ends to the opposite pore walls, so they can form pillaredlike structure. Two models are studied. In the first model, the nonterminating segments of chains can change their configuration inside the pore upon adsorption of spherical species. In the second model, the chains configuration remains fixed, so that the system is similar to a nonuniform quenched-annealed mixture. We study capillary condensation of fluid species inside such modified pores and compare the results obtained for two models.     

New nanosized catalytic membrane reactors for hydrogenation with stored hydrogen: Prerequisites and the experimental basis for their creation

The prerequisites and prospects for creating a new generation of nanosized membrane reactors are considered. For the first time, hydrogenation reactions take place in ceramic membrane pores with hydrogen adsorbed beforehand in mono- and multilayered oriented carbon nanotubes with graphene walls (OCNTGs) formed on the internal pore surface. It is shown for Trumem microfiltration membranes with D _avg ∼130 nm that oxidation reactions of CO on a Cu_0.03Ti_0.97O_{2 ± δ} catalyst and the oxidative conversion of methane into synthesis gas and light hydrocarbons on La + Ce/MgO are considerably enhanced when they occur in membranes. Regularities of hydrogen adsorption, storage, and desorption in nanosized membrane reactors are investigated through OCNTG formation in Trumem ultrafiltration membrane pores with D _avg = 50 and 90 nm and their saturation with hydrogen at a pressure of 10–13 MPa. It is shown that the amount of adsorbed hydrogen reaches 14.0% of OCNTG mass. Using thermogravimetric analysis in combination with mass-spectrometric analysis, hydrogen adsorption in OCNTG is first determined and its desorption is found to proceed at atmospheric pressure at a temperature of ∼175°C. It is shown that adsorbed hydrogen affects the transport properties of the membranes, reducing their efficiency with respect to liquids by 4–26 times. This is indirect confirmation of its high activity, due apparently the dissociative mechanism of adsorption.     

甲烷在活性炭上裂解制氢研究

在连续流动石英固定床反应器上研究了甲烷在活性炭上裂解制氢的反应，并对反应前后活性炭的比表面积以及孔径分布等的变化进行了测定。结果表明，甲烷在五种活性炭上的裂解行为基本相同，反应初期转化率最高，随着反应进行转化率逐渐降低直至一个平稳的状态；降低甲烷分压和增加甲烷与活性炭的接触时间可提高甲烷转化率；温度的升高有利于初始转化率的提高，但不利于活性炭的稳定性；反应后活性炭比表面积、孔容及微孔孔容都明显降低，平均孔径增大，孔径分布向中孔方向迁移，说明甲烷的裂解导致了活性炭孔特别是微孔内的炭沉积以及进一步的孔堵塞。