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

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Pore size characterization of mesoporous materials by a thermodynamic approach: A curvature-dependent solid-fluid potential

A thermodynamic analysis of nitrogen adsorption in cylindrical pores of MCM-41 and SBA-15 samples at 77 K is presented within the framework of the Broekhoff and de Boer (BdB) theory. We accounted for the effect of the solid surface curvature on the potential exerted by the pore walls. The developed model is in quantitative agreement with the non-local density functional theory (NLDFT) for pores larger than 2 tun. This modified BdB theory accounting for the Curvature Dependent Potential (CDP-BdB) was applied to determine the pore size distribution (PSD) of a number of MCM-41 and SBA-15 samples on the basis of matching the equilibrium theoretical isotherm against the adsorption branch of the experimental isotherm. In all cases investigated the PSDs determined with the new approach are very similar to those determined with the non-local density functional theory also using the same basis of matching of theoretical isotherm against the experimental adsorption branch. The developed continuum theory is very simple in its utilization, suggesting that CDP-BdB could be used as an alternative tool to obtain PSD for mesoporous solids from the analysis of adsorption branch of adsorption isotherms of any sub-critical fluids.

Improvement of the Derjaguin-Broekhoff-de Boer theory for capillary condensation/evaporation of nitrogen in mesoporous systems and its implications for pore size analysis of MCM-41 silicas and related materials

In this work, we propose an improvement of the classical Derjaguin-Broekhoff-de Boer (DBdB) theory for capillary condensation/evaporation in mesoporous systems. The primary idea of this improvement is to employ the Gibbs-Tolman-Koenig-Buff equation to predict the surface tension changes in mesopores. In addition, the statistical film thickness (so-called t-curve) evaluated accurately on the basis of the adsorption isotherms measured for the MCM-41 materials is used instead of the originally proposed t-curve (to take into account the excess of the chemical potential due to the surface forces). It is shown that the aforementioned modifications of the original DBdB theory have significant implications for the pore size analysis of mesoporous solids. To verify our improvement of the DBdB pore size analysis method (IDBdB), a series of the calcined MCM-41 samples, which are well-defined materials with hexagonally ordered cylindrical mesopores, were used for the evaluation of the pore size distributions. The correlation of the IDBdB method with the empirically calibrated Kruk-Jaroniec-Sayari (KJS) relationship is very good in the range of small mesopores. So, a major advantage of the IDBdB method is its applicability for small mesopores as well as for the mesopore range beyond that established by the KJS calibration, i.e., for mesopore radii greater than similar to4.5 nm. The comparison of the IDBdB results with experimental data reported by Kruk and Jaroniec for capillary condensation/evaporation as well as with the results from nonlocal density functional theory developed by Neimark et al. clearly justifies our approach. Note that the proposed improvement of the classical DBdB method preserves its original simplicity and simultaneously ensures a significant improvement of the pore size analysis, which is confirmed by the independent estimation of the mean pore size by the powder X-ray diffraction method.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Pore size distribution analysis of activated carbons: Application of density functional theory using nongraphitized carbon black as a reference system

The application of nonlocal density functional theory (NLDFT) to determine pore size distribution (PSD) of activated carbons using a nongraphitized carbon black, instead of graphitized thermal carbon black, as a reference system is explored. We show that in this case nitrogen and argon adsorption isotherms in activated carbons are precisely correlated by the theory, and such an excellent correlation would never be possible if the pore wall surface was assumed to be identical to that of graphitized carbon black. It suggests that pore wall surfaces of activated carbon are closer to that of amorphous solids because of defects of crystalline lattice, finite pore length, and the presence of active centers.. etc. Application of the NLDFT adapted to amorphous solids resulted in quantitative description of N-2 and Ar adsorption isotherms on nongraphitized carbon black BP280 at their respective boiling points. In the present paper we determined solid-fluid potentials from experimental adsorption isotherms on nongraphitized carbon black and subsequently used those potentials to model adsorption in slit pores and generate a corresponding set of local isotherms, which we used to determine the PSD functions of different activated carbons. (c) 2005 Elsevier Ltd. All rights reserved.

Pd/C catalysts based on synthetic carbons with bi- and tri-modal pore-size distribution:applications in flow chemistry

Two new types of phenolic resin-derived synthetic carbons with bi-modal and tri-modal pore-size distributions were used as supports for Pd catalysts. The catalysts were tested in chemoselective hydrogenation and hydrodehalogenation reactions in a compact multichannel flow reactor. Bi-modal and tri-modal micro-mesoporous structures of the synthetic carbons were characterised by N2 adsorption. HR-TEM, PXRD and XPS analyses were performed for characterising the synthesised catalysts. N2 adsorption revealed that tri-modal synthetic carbon possesses a well-developed hierarchical mesoporous structure (with 6.5 nm and 42 nm pores), contributing to a larger mesopore volume than the bi-modal carbon (1.57 cm3 g-1versus 1.23 cm3 g-1). It was found that the tri-modal carbon promotes a better size distribution of Pd nanoparticles than the bi-modal carbon due to presence of hierarchical mesopore limitting the growth of Pd nanoparticles. For all the model reactions investigated, the Pd catalyst based on tri-modal synthetic carbon (Pd/triC) show high activity as well as high stability and reproducibility. The trend in reactivities of different functional groups over the Pd/triC catalyst follows a general order alkyne ≫ nitro > bromo ≫ aldehyde.

Pore size distributions derived from adsorption isotherms, immersion calorimetry, and isosteric heats: A comparative study

We compare the pore size distribution of a well-characterized activated carbon derived from model-dependent, adsorption integral equation (AIE) methods with those from model-independent, immersion calorimetry and isosteric heat analyses. The AIE approach applied to nitrogen gave a mean pore width of 0.57 nm; the CO2 distribution exhibited wider dispersion. Spherical model application to CO2 and diffusion limitations for nitrogen and argon were proposed as primary reasons for inconsistency. Immersion enthalpy revealed a sharp decrease in available area equivalent to a cut-off due to molecular exclusion when the accessible surface was assessed against probe kinetic diameter. Mean pore width was identified as 0.58 ± 0.02 nm, endorsing the underlying assumptions for the nitrogen-based AIE approach. A comparison of the zero-coverage isosteric heat of adsorption for various non-polar adsorptives by the porous test sample was compared with the same adsorptives in contact with a non-porous reference adsorbent, leading to an energy ratio or adsorption enhancement factor. A linear relationship between the energy ratio and probe kinetic diameter indicated a primary pore size at 0.59 nm. The advantage of this enthalpy, model-independent methods over AIE were due to no assumptions regarding probe molecular shape, and no assumptions for pore shape and/or connectivity.

Large-pore mesoporous SBA-15 silica particles with submicrometer size as stationary phases for high-speed CEC separation

A new mesoporous sphere-like SBA-15 silica was synthesized and evaluated in terms of its suitability as stationary phases for CEC. The unique and attractive properties of the silica particle are its submicrometer particle size of 400 nm and highly ordered cylindrical mesopores with uniform pore size of 12 nm running along the same direction. The bare silica particles with submicrometer size have been successfully employed for the normal-phase electrochromatographic separation of polar compounds with high efficiency (e.g., 210 000 for thiourea), which is matched well with its submicrometer particle size. The Van Deemeter plot showed the hindrance to mass transfer because of the existence of pore structure. The lowest plate height of 2.0 mu m was obtained at the linear velocity of 1.1 mm/s. On the other hand, because of the relatively high linear velocity (e.g., 4.0 mm/s) can be generated, high-speed separation of neutral compounds, anilines, and basic pharmaceuticals in CEC with C-18-modified SBA-15 silica as stationary phases was achieved within 36, 60, and 34 s, respectively.

Volume, Surface, Connectivity and Size Distribution of Soil Pore Space in CT Images: Comparison of Samples at Different Depths from Nearby Natural and Tillage Areas

The study of soil structure, i.e., the pores, is of vital importance in different fields of science and technology. Total pore volume (porosity), pore surface, pore connectivity and pore size distribution are some (probably the most important) of the geometric measurements of pore space. The technology of X-ray computed tomography allows us to obtain 3D images of the inside of a soil sample enabling study of the pores without disturbing the samples. In this work we performed a set of geometrical measures, some of them from mathematical morphology, to assess and quantify any possible difference that tillage may have caused on the soil. We compared samples from tilled soil with samples from a soil with natural vegetation taken in a very close area. Our results show that the main differences between these two groups of samples are total surface area and pore connectivity per unit pore volume.

The correlation of pore morphology, interconnectivity and physical properties of 3D ceramic scaffolds with bone ingrowth

In the design of tissue engineering scaffolds, design parameters including pore size, shape and interconnectivity, mechanical properties and transport properties should be optimized to maximize successful inducement of bone ingrowth. In this paper we describe a 3D micro-CT and pore partitioning study to derive pore scale parameters including pore radius distribution, accessible radius, throat radius, and connectivity over the pore space of the tissue engineered constructs. These pore scale descriptors are correlated to bone ingrowth into the scaffolds. Quantitative and visual comparisons show a strong correlation between the local accessible pore radius and bone ingrowth; for well connected samples a cutoff accessible pore radius of approximately 100 microM is observed for ingrowth. The elastic properties of different types of scaffolds are simulated and can be described by standard cellular solids theory: (E/E(0))=(rho/rho(s))(n). Hydraulic conductance and diffusive properties are calculated; results are consistent with the concept of a threshold conductance for bone ingrowth. Simple simulations of local flow velocity and local shear stress show no correlation to in vivo bone ingrowth patterns. These results demonstrate a potential for 3D imaging and analysis to define relevant pore scale morphological and physical properties within scaffolds and to provide evidence for correlations between pore scale descriptors, physical properties and bone ingrowth.

The effects of pore architecture in silk fibroin scaffolds on the growth and differentiation of mesenchymal stem cells expressing BMP7

The pore architecture of scaffolds is known to play a critical role in tissue engineering as it provides the vital framework for seeded cells to organize into a functioning tissue. In this report we have investigated the effects of different concentrations of silk fibroin protein on three-dimensional (3D) scaffold pore microstructure. Four pore size ranges of silk fibroin scaffolds were made by the freeze drying technique, with the pore sizes ranging from 50 to 300 lm. The pore sizes of the scaffolds decreased as the concentration of fibroin protein increased. Human bone marrow mesenchymal stromal cells (BMSC) transfected with the BMP7 gene were cultured in these scaffolds. A cell viability colorimetric assay, alkaline phosphatase assay and reverse transcription-polymerase chain reaction were performed to analyze the effect of pore size on cell growth, the secretion of extracellular matrix (ECM) and osteogenic differentiation. Cell migration in 3D scaffolds was confirmed by confocal microscopy. Calvarial defects in SCID mice were used to determine the bone forming ability of the silk fibroin scaffolds incorporating BMSC expressing BMP7. The results showed that BMSC expressing BMP7 preferred a pore size between 100 and 300 lm in silk fibroin protein fabricated scaffolds, with better cell proliferation and ECM production. Furthermore, in vivo transplantation of the silk fibroin scaffolds combined with BMSC expressing BMP7 induced new bone formation. This study has shown that an optimized pore architecture of silk fibroin scaffolds can modulate the bioactivity of BMP7-transfected BMSC in bone formation.

The effects of pore architecture in silk fibroin scaffolds on the growth and differentiation of BMP7-expressing mesenchymal stem cells

Pore architecture of scaffolds is known to play a critical role in tissue engineering as it provides the vital framework for the seeded cells to organize into a functioning tissue. In this report, we investigated the effects of different concentration on silk fibroin protein 3D scaffold pore microstructure. Four pore size ranges of silk fibroin scaffolds were made by freeze-dry technique, with the pore sizes ranging from 50 to 300 µm. The pore size of the scaffold decreases as the concentration increases. Human mesenchymal stem cells were in vitro cultured in these scaffolds. After BMP7 gene transferred, DNA assay, ALP assay, hematoxylin–eosin staining, alizarin red staining and reverse transcription-polymerase chain reaction were performed to analyze the effect of the pore size on cell growth, differentiation and the secretion of extracellular matrix (ECM). Cell morphology in these 3D scaffolds was investigated by confocal microscopy. This study indicates mesenchymal stem cells prefer the group of scaffolds with pore size between 100 and 300 µm for better proliferation and ECM production

Pore formation during dehydration of polycrystalline gypsum observed and quantified in a time-series synchrotron radiation based X-ray micro-tomography experiment

We conducted an in-situ X-ray micro-computed tomography heating experiment at the Advanced Photon Source (USA) to dehydrate an unconfined 2.3 mm diameter cylinder of Volterra Gypsum. We used a purpose-built X-ray transparent furnace to heat the sample to 388 K for a total of 310 min to acquire a three-dimensional time-series tomography dataset comprising nine time steps. The voxel size of 2.2 μm3 proved sufficient to pinpoint reaction initiation and the organization of drainage architecture in space and time. We observed that dehydration commences across a narrow front, which propagates from the margins to the centre of the sample in more than four hours. The advance of this front can be fitted with a square-root function, implying that the initiation of the reaction in the sample can be described as a diffusion process. Novel parallelized computer codes allow quantifying the geometry of the porosity and the drainage architecture from the very large tomographic datasets (20483 voxels) in unprecedented detail. We determined position, volume, shape and orientation of each resolvable pore and tracked these properties over the duration of the experiment. We found that the pore-size distribution follows a power law. Pores tend to be anisotropic but rarely crack-shaped and have a preferred orientation, likely controlled by a pre-existing fabric in the sample. With on-going dehydration, pores coalesce into a single interconnected pore cluster that is connected to the surface of the sample cylinder and provides an effective drainage pathway. Our observations can be summarized in a model in which gypsum is stabilized by thermal expansion stresses and locally increased pore fluid pressures until the dehydration front approaches to within about 100 μm. Then, the internal stresses are released and dehydration happens efficiently, resulting in new pore space. Pressure release, the production of pores and the advance of the front are coupled in a feedback loop.