Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

Graphene Oxide: Study of Pore Size Distribution and Surface Chemistry Using Immersion Calorimetry

In this work, the textural parameters of graphene oxide (GO) and graphite (Gr) samples were determined. The non-local density functional theory (NLDFT) and quenched solid density functional theory (QSDFT) kernels were used to evaluate the pore size distribution (PSD) by modeling the pores as slit, cylinder and slit-cylinder. The PSD results were compared with the immersion enthalpies obtained using molecules with different kinetic diameter (between 0.272 nm and 1.50 nm). Determination of immersion enthalpy showed to track PSD for GO and graphite (Gr), which was used as a comparison solid. Additionally, the functional groups of Gr and GO were determined by the Boehm method. Donor number (DN) Gutmann was used as criteria to establish the relationship between the immersion enthalpy and the parameter of the probe molecules. It was found that according to the Gutmann DN the immersion enthalpy presented different values that were a function of the chemical groups of the materials. Finally, the experimental and modeling results were critically discussed.

The Structural Properties and Diffusion of a Three-Dimensional Isotropic Core-Softened Model Fluid in Disordered Porous Media. Molecular Dynamics Simulation

The microscopic structure and dynamic properties of an isotropic three-dimensional core-softened model fluid in disordered matrices of Lennard-Jones particles have been studied. Molecular dynamics computer simulations in Grand Canonical ensemble were used as the methodological tools. It was shown that the microscopic structure of the fluid is characterized by anomalies similar to those found in a bulk model, but that it is affected by the fluid-matrix interactions. The dynamic properties also exhibit anomalous dependence on fluid density, but the magnitude of these anomalies is suppressed in comparison to the bulk fluid model. The anomalous behaviour of the diffusion coefficient is attributed to structural changes in the first coordination shell of a given fluid particle. It seems that the anomalies can only be suppressed at matrix densities which are higher than those studied in the present work.

Assessment of MOF's Quality: Quantifying Defect Content in Crystalline Porous Materials

A quantitative method for assessment of defects in metal-organic framework (MOF) is presented based on isotherms calculated using Grand Canonical Monte Carlo (GCMC) simulations. Defects in MOF structures generated during the synthesis and sample preparation can lead to large variations in experimentally measured adsorption isotherms but are difficult to quantify. We use as a case study CO2 adsorption on Cu3(BTC)2 MOF (BTC = benzene-1,3,5-tricarboxylic acid) to show that different samples reported in the literature have various proportions of principal pores blocked or side pores blocked, resulting in isotherms with different capacity and affinity toward CO2. The approach presented is easily generalized to other materials, showing that simulation results combined with experimentally measured gas adsorption isotherms can be used to quantitatively identify key defective features of the material.

Capillary condensation in pores with rough walls: A density functional approach

The effect of surface roughness of slit-like pore walls on the capillary condensation of a spherical particles and short chains is studied. The gas molecules interact with the substrate by a Lennard-Jones (9,3) potential. The rough layer at each pore wall has a variable thickness and density and consists of a disordered quenched matrix of spherical particles. The system is described in the framework of a density functional approach and using computer simulations. The contribution due to attractive van der Waals interactions between adsorbate molecules is described by using first-order mean spherical approximation and mean-field approximation.

Features of nitrogen adsorption on nonporous carbon and silica surfaces in the framework of classical density functional theory

Equilibrium adsorption data of nitrogen on a series of nongraphitized carbon blacks and nonporous silica at 77 K were analyzed by means of classical density functional theory to determine the solid-fluid potential. The behavior of this potential profile at large distance is particularly considered. The analysis of nitrogen adsorption isotherms seems to indicate that the adsorption in the first molecular layer is localized and controlled mainly by short-range forces due to the surface roughness, crystalline defects, and functional groups. At distances larger than approximately 1.3-1.5 molecular diameters, the adsorption is nonlocalized and appears as a thickening of the adsorbed film with increasing bulk pressure in a relatively weak adsorption potential field. It has been found that the asymptotic decay of the potential obeys the power law with the exponent being -3 for carbon blacks and -4 for silica surface, which signifies that in the latter case the adsorption potential is mainly exerted by surface oxygen atoms. In all cases, the absolute value of the solid-fluid potential is much smaller than that predicted by the Lennard-Jones pair potential with commonly used solid-fluid molecular parameters. The effect of surface heterogeneity on the heat of adsorption is also discussed.

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.