Effect of pore constriction on adsorption behaviour in nanoporous carbon slit pore: a computer simulation

A Grand Canonical Monte Carlo simulation (GCMC) method is used to study the effects of pore constriction on the adsorption of argon at 87.3 K in carbon slit pores of infinite and finite lengths. It is shown that the pore constriction affects the pattern of adsorption isotherm. First, the isotherm of the composite pore is greater than that of the uniform pore having the same width as the larger cavity of the composite pore. Secondly, the hysteresis loop of the composite pore is smaller than and falls between those of uniform pores. Two types of hysteresis loops have been observed, irrespective of the absence or presence of constriction and their presence depend on pore width. One hysteresis loop is associated with the compression of adsorbed particles and this phenomenon occurs after pore has been filled with particles. The second hysteresis loop is the classical condensation-evaporation loop. The hysteresis loop of a composite pore depends on the sizes of the larger cavity and the constriction. Generally, it is found that the pore blocking effect is not manifested in composite slit pores, and this result does not support the traditional irkbottle pore hypothesis.

The effects of energy sites on adsorption of Lennard–Jones fluids and phase transition in carbon slit pore of finite length a computer simulation study

A Monte Carlo simulation method is Used 10 study the effects of adsorption strength and topology of sites on adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length. Argon is used as a model adsorbate, while the adsorbent is modeled as a finite carbon slit pore whose two walls composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. Impurities having well depth of interaction greater than that of carbon atom are assumed to be grafted onto the surface. Different topologies of the impurities; corner, centre, shelf and random topologies are studied. Adsorption isotherms of argon at 87.3 K are obtained for pore having widths of 1, 1.5 and 3 11111 using a Grand Canonical Monte Carlo simulation (GCMC). These results are compared with isotherms obtained for infinite pores. It is shown that the Surface heterogeneity affects significantly the overall adsorption isotherm, particularly the phase transition. Basically it shifts the onset of adsorption to lower pressure and the adsorption isotherms for these four impurity models are generally greater than that for finite pore. The positions of impurities on solid Surface also affect the shape of the adsorption isotherm and the phase transition. We have found that the impurities allocated at the centre of pore walls provide the greatest isotherm at low pressures. However when the pressure increases the impurities allocated along the edges of the graphene layers show the most significant effect on the adsorption isotherm. We have investigated the effect of surface heterogeneity on adsorption hysteresis loops of three models of impurity topology, it shows that the adsorption branches of these isotherms are different, while the desorption branches are quite close to each other. This suggests that the desorption branch is either the thermodynamic equilibrium branch or closer to it than the adsorption branch. (c) 2005 Elsevier Inc. All rights reserved.

Adsorption of Lennard-Jones fluids in carbon slit pores of a finite length. A computer simulation study

The adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length was studied with Canonical Ensemble (NVT) and Gibbs Ensemble Monte Carlo Simulations (GEMC). The Canonical Ensemble was a collection of cubic simulation boxes in which a finite pore resides, while the Gibbs Ensemble was that of the pore space of the finite pore. Argon was used as a model for Lennard-Jones fluids, while the adsorbent was modelled as a finite carbon slit pore whose two walls were composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. The Lennard-Jones (LJ) 12-6 potential model was used to compute the interaction energy between two fluid particles, and also between a fluid particle and a carbon atom. Argon adsorption isotherms were obtained at 87.3 K for pore widths of 1.0, 1.5 and 2.0 nm using both Canonical and Gibbs Ensembles. These results were compared with isotherms obtained with corresponding infinite pores using Grand Canonical Ensembles. The effects of the number of cycles necessary to reach equilibrium, the initial allocation of particles, the displacement step and the simulation box size were particularly investigated in the Monte Carlo simulation with Canonical Ensembles. Of these parameters, the displacement step had the most significant effect on the performance of the Monte Carlo simulation. The simulation box size was also important, especially at low pressures at which the size must be sufficiently large to have a statistically acceptable number of particles in the bulk phase. Finally, it was found that the Canonical Ensemble and the Gibbs Ensemble both yielded the same isotherm (within statistical error); however, the computation time for GEMC was shorter than that for canonical ensemble simulation. However, the latter method described the proper interface between the reservoir and the adsorbed phase (and hence the meniscus).

High-pressure adsorption capacity and structure of CO2 in carbon slit pores: Theory and simulation

We present new simulation results for the packing of single-center and three-center models of carbon dioxide at high pressure in carbon slit pores. The former shows a series of packing transitions that are well described by our density functional theory model developed earlier. In contrast, these transitions are absent for the three-center model. Analysis of the simulation results shows that alternations of flat-lying molecules and rotated molecules can occur as the pore width is increased. The presence or absence of quadrupoles has negligible effect on these high-density structures.

Comparative adsorption of spherical argon and flexible n-butane in carbon slit pores - a GCMC computer simulation study

In this paper we investigate the difference between the adsorption of spherical molecule argon (at 87.3 K) and the flexible normal butane (at an equivalent temperature of 150 K) in carbon slit pores. These temperatures are equivalent in the sense that they have the same relative distances between their respective triple points and critical points. Higher equivalent temperatures are also studied (122.67 K for argon and 303 K for n-butane) to investigate the effects of temperature on the 2D-transition in adsorbed density. The Grand Canonical Monte Carlo simulation is used to study the adsorption of these two model adsorbates. Beside the longer computation times involved in the computation of n-butane adsorption, n-butane exhibits many interesting behaviors such as: (i) the onset of adsorption occurs sooner (in terms of relative pressure), (ii) the hysteresis for 2D- and 3D-transitions is larger, (iii) liquid-solid transition is not possible, (iv) 2D-transition occurs for n-butane at 150 K while it does not happen for argon except for pores that accommodate two layers of molecules, (v) the maximum pore density is about four times less than that of argon and (vi) the sieving pore width is slightly larger than that for argon. Finally another feature obtained from the Grand Canonical Monte Carlo (GCMC) simulation is the configurational arrangement of molecules in pores. For spherical argon, the arrangement is rather well structured, while for n-butane the arrangement depends very much on the pore size. (C) 2004 Elsevier B.V. All rights reserved.

Effect of activated carbon xerogel pore size on the capacitance performance of ionic liquid electrolytes

The use of ionic liquid (IL) electrolytes promises to improve the energy density of electrochemical capacitors (ECs) by allowing for operation at higher voltages. Several studies have also shown that the pore size distribution of materials used to produce electrodes is an important factor in determining EC performance. In this research the capacitative, energy and power performance of ILs 1-ethyl-3- methylimidazolium tetrafluoroborate (EMImBF4), 1-ethyl-3-methylimidazolium dicyanamide (EMImN(CN)2), 1,2-dimethyl-3-propylimidazolium bis(trifluoromethylsulfonyl)imide (DMPImTFSI), and 1-butyl-3-methylpyrrolidinium tris(pentafluoroethyl)trifluorophosphate (BMPyT(F5Et)PF3) were studied and compared with the commercially utilised organic electrolyte 1M tetraethylammonium tetrafluoroborate solution in anhydrous propylene carbonate (Et4NBF4–PC 1 M). To assess the effect of pore size on IL performance, controlled porosity carbons were produced from phenolic resins activated in CO2. The carbon samples were characterised by nitrogen adsorption– desorption at 77 K and the relevant electrochemical behaviour was characterised by cyclic voltammetry, galvanostatic charge–discharge and electrochemical impedance spectroscopy. The best capacitance performance was obtained for the activated carbon xerogel with average pore diameter 3.5 nm, whereas the optimum rate performance was obtained for the activated carbon xerogel with average pore diameter 6 nm. When combined in an EC with IL electrolyte EMImBF4 a specific capacitance of 210 F g1 was obtained for activated carbon sample with average pore diameter 3.5 nm at an operating voltage of 3 V. The activated carbon sample with average pore diameter 6 nm allowed for maximum capacitance retention of approximately 70% at 64 mA cm2.

Adsorption of quadrupolar, diatomic nitrogen onto graphitized thermal carbon black and in slit-shaped carbon pores. Effects of surface mediation

In this paper, we study the effect of solid surface mediation on the intermolecular potential energy of nitrogen, and its impact on the adsorption of nitrogen on a graphitized carbon black surface and in carbon slit-shaped pores. This effect arises from the lower effective interaction potential energy between two particles close to the surface compared to the potential energy of the same two particles when they are far away from the surface. A simple equation is proposed to calculate the reduction factor and this is used in the Grand Canonical Monte Carlo (GCMC) simulation of nitrogen adsorption on graphitized thermal carbon black. With this modification, the GCMC simulation results agree extremely well with the experimental data over a wide range of pressure; the simulation results with the original potential energy (i.e. no surface mediation) give rise to a shoulder in the neighbourhood of monolayer coverage and a significant over-prediction of the second and higher layer coverages. The influence of this surface mediation on the dependence of the pore-filling pressure on the pore width is also studied. It is shown that such surface mediation has a significant effect on the pore-filling pressure. This implies that the use of the local isotherms obtained from the potential model without surface mediation could give rise to a serious error in the determination of the pore-size distribution.

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.

BET surface area of carbonaceous adsorbents: verification using geometric considerations and GCMC simulations on virtual porous carbon models

The applicability of BET model for calculation of surface area of activated carbons is checked by using molecular simulations. By calculation of geometric surface areas for the simple model carbon slit-like pore with the increasing width, and by comparison of the obtained values with those for the same systems from the VEGA ZZ package (adsorbate-accessible molecular surface), it is shown that the latter methods provide correct values. For the system where a monolayer inside a pore is created the ASA approach (GCMC, Ar, T = 87 K) underestimates the value of surface area for micropores (especially, where only one layer is observed and/or two layers of adsorbed Ar are formed). Therefore, we propose the modification of this method based on searching the relationship between the pore diameter and the number of layers in a pore. Finally BET; original andmodified ASA; and A, B and C-point surface areas are calculated for a series of virtual porous carbons using simulated Ar adsorption isotherms (GCMC and T = 87 K). The comparison of results shows that the BET method underestimates and not, as it was usually postulated, overestimates the surface areas of microporous carbons.

Dissolved inorganic carbon isotopes in pore waters of ODP Leg 204 sediments

Isotopic characterization of carbon in the dissolved inorganic carbon (DIC) pool is fundamental for a wide array of scientific studies directly related to gas hydrate research. In order to generate integrated and internally consistent data of d13C of DIC in pore waters from Hydrate Ridge, we used the modern continuous flow technology of a GasBench II automated sampler interfaced to a gas source stable isotope mass spectrometer for the rapid determination (~80 samples/day) of d13C DIC in small-volume water samples. The overall precision of this technique is conservatively estimated to be better than ±0.15 per mil (1 sigma), which is similar to the precision of methods in current use. Here we present the data generated from Ocean Drilling Program Leg 204 pore water samples.