The effects of curvature and surface heterogeneity on the adsorption of water in finite length carbon nanopores: a computer simulation study

The effects of pore curvature and surface heterogeneity on the adsorption of water on a graphitic surface at 298 K were investigated using a Grand Canonical Monte Carlo (GCMC) simulation. Slit and cylindrical pores are used to study the curvature effects. To investigate the surface heterogeneity the functional group and the structural defect on the surface were specifically considered. The hydroxyl group (OH) is used as a model for the functional group and the water potential model proposed by Müller et al. is used to calculate the water interaction. For the homogeneous cylinder, the pore filling occurs at a pressure lower than the saturation pressure of the water model, while it is greater in the case of homogeneous slit pore. The size of hysteresis loop is more sensitive to the length of cylinder than that of the slit, and it increases with decreasing pore length. The isotherms of water in cylindrical pores are found to depend on the position and the concentration of the functional group. The pore filling pressure is lower with an increased number and/or with the position of the functional group. The structural defect shows significant effects on the adsorption isotherm in shifting to a lower pore filling pressure when it is located at a position away from the pore entrance. The adsorption of water on the heterogeneous surface was studied and it was found that the simulated isotherm can describe the behaviour of water on Graphitized Thermal Carbon Black (GTCB) satisfactorily. The water cluster grows mostly along the surface for the case of finite extent surface, while for the slit the pore grows in all directions but the preference is a direction perpendicular to the pore wall. Reasons for the direction of growth will be discussed.     

Explanation of the unusual peak of calorimetric heat in the adsorption of nitrogen, argon and methane on graphitized thermal carbon black

Heats of adsorption and adsorption isotherms of argon, nitrogen and methane on a perfect graphitic surface and a defective graphitic surface are studied with a Grand Canonical Monte Carlo Simulation (GCMC). For the perfect surface, the isosteric heat versus loading shows a typical pattern of adsorption of simple fluids on graphite. Depending on adsorbate, degree of graphitization and temperature, a spike in the heat curve versus loading is observed when the first layer is mostly covered with adsorbate molecules. The heat spike is observed for argon and nitrogen at 77 K while for argon at 87.3 K it is no longer present. These simulation results are consistent with the experimental data of J. Rouquerol, S. Partyka and F. Rouquerol, J. Chem. Soc., Faraday Trans. 1, 1977, 73, 306. In the case of methane we observe heat spikes at low temperatures, 84.5, 92.5 and 104 K. The heat spike shifts to higher loading with temperature and it then disappears at high temperatures. These observations are in qualitative agreement with the experimental data of A. Inaba, Y. Koga and J. A. Morrison, J. Chem. Soc., Faraday Trans. 2, 1986, 82, 1635. In all cases where heat spikes are observed, the GCMC simulation results indicate that the heat spike is associated with the squeezing of molecules into the already dense first layer, and the rearrangement of molecules to form a highly structured fluid of this layer. While this squeezing into the first layer is happening, molecules continue to adsorb onto the relatively sparse second layer.     

Monte Carlo Simulation and Experimental Studies of CO2, CH4 and Their Mixture Capture in Porous Carbons

Adsorption of carbon dioxide and methane in porous activated carbon and carbon nanotube was studied experimentally and by Grand Canonical Monte Carlo (GCMC) simulation. A gravimetric analyzer was used to obtain the experimental data, while in the simulation we used graphitic slit pores of various pore size to model activated carbon and a bundle of graphitic cylinders arranged hexagonally to model carbon nanotube. Carbon dioxide was modeled as a 3-center-Lennard-Jones (LJ) molecule with three fixed partial charges, while methane was modeled as a single LJ molecule. We have shown that the behavior of adsorption for both activated carbon and carbon nanotube is sensitive to pore width and the crossing of isotherms is observed because of the molecular packing, which favors commensurate packing for some pore sizes. Using the adsorption data of pure methane or carbon dioxide on activated carbon, we derived its pore size distribution (PSD), which was found to be in good agreement with the PSD obtained from the analysis of nitrogen adsorption data at 77 K. This derived PSD was used to describe isotherms at other temperatures as well as isotherms of mixture of carbon dioxide and methane in activated carbon and carbon nanotube at 273 and 300 K. Good agreement between the computed and experimental isotherm data was observed, thus justifying the use of a simple adsorption model.     

The role of accessibility in the characterization of porous solids and their adsorption properties

This paper addresses the role of accessibility for adsorption in porous solids on the adsorption properties including Henry constant, adsorption isotherms and isosteric heat of adsorption. The relevant parameters are the accessible volume, the accessible geometrical surface area and the accessible pore size and its associated volume. This concept will be demonstrated to be important and calls for the need to consider adsorption characteristics in the most coherent and consistent manner. It is particularly reinforced by the limitations inherent in the conventional ways in determining the void volume, surface area and pore size. We provide a number of examples to support this; the challenge that faces us is the development of consistent experimental procedures to determine these accessible quantities. We define the accessible pore size as the size of the largest sphere that rests on three closest solid atoms in such a manner that any probe particle residing in that sphere would have a non-positive solid-fluid potential energy. For each accessible pore size there is an associated accessible pore volume, giving rise to a new accessible pore size distribution (APSD). This is distinct from the classical pore size distribution commonly used in the literature, and in our definition of accessible pore size, a zero pore size is possible. It is also emphasized that the accessible quantities that we introduce here are dependent on the choice of molecular probe, which is entirely consistent with the concept of molecular sieving.     

Adsorption of water in finite length carbon slit pore: comparison between computer simulation and experiment

The effects of surface dimensions and topology on the adsorption of water on a graphite surface at 298 K were investigated using the grand canonical Monte Carlo (GCMC) simulation. Regarding the surface topology, we specifically considered the functional group and its position on the surface. The hydroxyl group (OH) is used as a model for the functional group. For describing the interaction of water, we used the potential model proposed by Muller et al., and the simulated isotherms of water in slit pores are found to depend on the position and concentration of the functional group. The onset of adsorption shifts to lower pressure when the concentration of functional group increases or when the functional group is positioned at the center of the graphene surface. The configuration of a group of functional groups also affects the adsorption isotherm. In all cases investigated, we have found that the hysteresis loop always exists, and the loop size depends on the concentration of the functional group and its position. Finally, we tested the molecular model of water adsorption on a functional graphite pore against the experimental data of a commercial activated carbon. The agreement is found to be satisfactory when the model porous solid is composed of pores having width in the range between 10 and 20 A and functional groups positioned at the center of the graphitic wall.     

A critical assessment of capillary condensation and evaporation equations: a computer simulation study

Nonaqueous reverse micelles of brij surfactants are prepared in benzene and ethylammonium nitrate (EAN). The effect of polar head group bulk on reverse micellar size was studied with brij-52, brij-56 and brij-58 whereas the effect of polarity of hydrocarbon chain was investigated taking brij-52 and brij-93 with varying W(s) (W(s)=[EAN]/[surfactant]). Dynamic light scattering (DLS) has been employed to reveal the size and shape of the reverse micelles. Micropolarities of these reverse micelles were investigated by visible spectroscopy using methylene blue (MB) and methyl orange (MO) as molecular optical probes. It has been revealed from the experimental results that with increase in polar head group size reverse micellar size increases. Moreover, it is also observed that with increasing polarity of the hydrocarbon chain the average size of the reverse micelles decreases. It can be concluded that polar head group size and polarity of hydrocarbon chain play important roles in determining reverse micellar size of the brij surfactants apart from the W(s) ratio, nature of the solvent medium, and concentration of the surfactants.     

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 to 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, shell 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 nm 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.