Simulated water adsorption in chemically heterogeneous carbon nanotubes

Grand canonical Monte Carlo simulations are used to study the adsorption of water in single-walled (10:10), (12:12), and (20:20) carbon nanotubes at 298 K. Water is represented by the extended simple point charge model and the carbon atoms as Lennard-Jones spheres. The nanotubes are decorated with different amounts of oxygenated sites, represented as carbonyl groups. In the absence of carbonyl groups the simulated isotherms are characterized by negligible amounts of water uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption-desorption hysteresis loops. In the presence of a few carbonyl groups the simulated adsorption isotherms are characterized by pore filling at lower pressures and by narrower adsorption-desorption hysteresis loops compared to the results obtained in the absence of carbonyl groups. Our results show that the distribution of the carbonyl groups has a strong effect on the adsorption isotherms. For carbonyl groups localized in a narrow section the adsorption of water may be gradual because a cluster of adsorbed water forms at low pressures and grows as the pressure increases. For carbonyl groups distributed along the nanotube the adsorption isotherm is of type V.     

Water in carbon nanotubes: adsorption isotherms and thermodynamic properties from molecular simulation

Grand canonical Monte Carlo simulations are performed to study the adsorption of water in single-walled (6:6), (8:8), (10:10), (12:12), and (20:20) carbon nanotubes in the 248-548 K temperature range. At room temperature the resulting adsorption isotherms in (10:10) and wider single-walled carbon nanotubes (SWCNs) are characterized by negligible water uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption/desorption hysteresis loops. The width of the hysteresis loops decreases as pore diameter narrows and it becomes negligible for water adsorption in (8:8) and (6:6) SWCNs. Results for the isosteric heat of adsorption, density profiles along the pore axis and across the pore radii, order parameter across the pore radii, and x-ray diffraction patterns are presented. Layered structures are observed when the internal diameter of the nanotubes is commensurate to the establishment of a hydrogen-bonded network. The structure of water in (8:8) and (10:10) SWCNs is ordered when the temperature is 298 and 248 K, respectively. By simulating adsorption isotherms at various temperatures, the hysteresis critical temperature, e.g., the lowest temperature at which no hysteresis can be detected, is determined for water adsorbed in (20:20), (12:12), and (10:10) SWCNs. The hysteresis critical temperature is lower than the vapor-liquid critical temperature for bulk Simple Point Charge-Extended (SPCE) water model.     

Water adsorption-desorption isotherms of two-dimensional hexagonal mesoporous silica around freezing point

Zr-doped mesoporous silica with a diameter of approximately 3.8 nm was synthesized via an evaporation-induced self-assembly process, and the adsorption-desorption isotherms of water vapor were measured in the temperature range of 263-298 K. The measured adsorption-desorption isotherms below 273 K indicated that water confined in the mesopores did not freeze at any relative pressure. All isotherms had a steep curve, resulting from capillary condensation/evaporation, and a pronounced hysteresis. The hysteresis loop, which is associated with a delayed adsorption process, increased with a decrease in temperature. Furthermore, the curvature radius where capillary evaporation/condensation occurs was evaluated by the combined Kelvin and Gibbs-Tolman-Koening-Buff (GTKB) equations for the modification of the interfacial tension due to the interfacial curvature. The thickness of the water adsorption layer for capillary condensation was slightly larger, whereas that for capillary evaporation was slightly smaller than 0.7 nm.     

Adsorption isotherms of linoleic and linolenic acids from solutions in carbon tetrachloride on surfaces of highly dispersed manganese and copper ferrites

A comparative study of adsorption-desorption isotherms of linoleic and linolenic acids from solutions in carbon tetrachloride on surfaces of manganese and copper ferrites is performed by means of equilibrium adsorption. Adsorption isotherms of fatty acids are described in terms of the theory of volume filling of micropores, and the values of the limiting adsorption, the characteristic adsorption energy, and the pore space volume are calculated. It is established that the limiting adsorption values of linoleic and linolenic acids from solutions in carbon tetrachloride on a copper ferrite surface are higher than on a manganese ferrite surface. It is shown that the adsorption-desorption isotherms have a hysteresis loop.     

Oleic acid adsorption-desorption isotherms on the surface of high-dispersity ferrites from a solution in carbon tetrachloride

The equilibrium adsorption method was used to comparatively study the adsorption-desorption isotherms of oleic acid on the surfaces of manganese and copper ferrites from a solution in carbon tetrachloride. The adsorption isotherms of the fatty acid were described in terms of the theory of volume filling of micropores. The theory was used to calculate the limiting adsorption values, characteristic energy, and porous space volumes. The isotherm of oleic acid adsorption on the surface of manganese ferrite from a solution in carbon tetrachloride was similar to the isotherms of fatty acid adsorption from solutions in heptane, whereas the isotherm of adsorption on the surface of copper ferrite was similar to the isotherms of fatty acid adsorption from hexane. The limiting adsorption from carbon tetrachloride was higher on the surface of manganese ferrite than on the surface of copper ferrite. The adsorption-desorption isotherms contained hysteresis loops.     

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.     

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.     

Gravimetric,Volumetric And Calorimetric Studies of The Surface Structure of Portland Cement

Dry cement powder and hardened cement paste were characterised by means of laser granulometry and volumetric measurement of the nitrogen adsorption isotherm at 77 K. Water sorption isotherms at ambient temperature were measured stepwise by means of a gravimetric apparatus. The isotherms show a very large hysteresis loop, reaching down to zero relative pressure, which reflects swelling of the cement gel. Thermoporometric measurements reveal that swelling of that gel consists in water take up within the continuously growing structure; no stable pore structure could be observed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Water adsorption in disordered mesoporous silica (Vycor) at 300 K and 650 K: a Grand Canonical Monte Carlo simulation study of hysteresis

This numerical simulation paper focuses on the adsorption/desorption of water in disordered mesoporous silica glasses (Vycor-like). The numerical adsorbent was previously obtained by off lattice method, and was shown to reproduce quite well the micro- and mesotextural properties of real Vycor, as well as morphological (pore size distribution) and topological (pore interconnections) disorder. The water-water interactions are described by the SPC model while water-silica interactions are calculated in the framework of the PN-TrAZ model. The water adsorption/desorption isotherms and the configurational energies are calculated by the Grand Canonical Monte Carlo simulation method. The low pressure results compare well with experiments, showing the good transferability of the intermolecular potential. It is shown that if the hysteresis loop observed in the adsorption/desorption isotherm is considered as a true phase transition (which is actually still an open question in the case of disordered porous materials), then it is possible to calculate the grand potential by applying the thermodynamic integration scheme. The grand potential is shown to be multivalued for low (subcritical) temperature, and continuous for high (supercritical) temperature. A coexistence point is found within the hysteresis loop, actually close to the vertical desorption line. Below the equilibrium chemical potential, the gaslike branch is stable whereas the liquidlike branch is metastable. The situation is reversed above the coexistence point.     

Calculation of the adsorption-desorption hysteresis curves in narrow-pore systems

The gas and liquid spinodal branches for an adsorbate located in narrow slit-shaped, cylindrical, and spherocylindrical pores were calculated. The adsorbate is modeled by Lennard-Jones spherical particles. The calculation was based on the lattice gas model taking into account the intermolecular interactions of nearest neighbors in the quasichemical approximation. The density-temperature diagrams for the gas and liquid spinodal branches in the pores are similar to the equilibrium vapor-liquid phase diagrams: they have a common critical point; the dense-phase branches are shifted to lower pore fillings, while the rarefied-phase branches are shifted toward higher pore fillings. The width of adsorption-desorption hysteresis loop in the adsorption isotherms for Lennard-Jones particles was analyzed as a function of the pore size and the interaction potential of the adsorbate with the pore walls. The effect of pore wall roughness and the accuracy of isotherm calculation on the width of the adsorption-desorption hysteresis loop in cylindrical pores is discussed Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 813–823, May, 2007.     

Morphology and surface properties of fumed silicas

Several series of fumed silicas and mixed fumed oxides produced and treated under different conditions were studied in gaseous and liquid media using nitrogen and water adsorption-desorption, mass spectrometry, FTIR, NMR, thermally stimulated depolarization current (TSDC), photon correlation spectroscopy (PCS), zeta potential, potentiometric titration, and Auger electron spectroscopy methods. Aggregation of primary particles and adsorption capacity (Vp) decrease and hysteresis loops of nitrogen adsorption-desorption isotherms becomes shorter with decreasing specific surface area (S(BET)). However, the shape of nitrogen adsorption-desorption isotherms can be assigned to the same type independent of S(BET) value. The main maximum of pore size distribution (gaps between primary nonporous particles in aggregates and agglomerates) shifts toward larger pore size and its intensity decreases with decreasing S(BET) value. The water adsorption increases with increasing S(BET) value; however, the opposite effect is observed for the content of surface hydroxyls (in mmol/m2). Associative desorption of water (2(SiOH)-->SiOSi+H2O) depends on both the morphology and synthesis conditions of fumed silica. The silica dissolution rate increases with increasing S(BET) and pH values. However, surface charge density and the modulus of zeta-potential increase with decreasing S(BET) value. The PCS, 1H NMR, and TSDC spectra demonstrate rearrangement of the fumed silica dispersion depending on the S(BET) value and the silica concentration (C(SiO2)) in the aqueous suspensions. A specific state of the dispersion is observed at the C(SiO2) values corresponding to the bulk density of the initial silica powder.     

Adsorption of simple gases in MCM-41 materials: the role of surface roughness

This paper reports the development and testing of atomistic models of silica MCM-41 pores. Model A is a regular cylindrical pore having a constant section. Model B has a surface disorder that reproduces the morphological features of a pore obtained from an on-lattice simulation that mimics the synthesis process of MCM-41 materials. Both models are generated using a similar procedure, which consists of carving the pore out of an atomistic silica block. The differences between the two models are analyzed in terms of small angle neutron scattering spectra as well as adsorption isotherms and isosteric heat curves for Ar at 87 K and Xe at 195 K. As expected for capillary condensation in regular nanopores, the Ar and Xe adsorption/desorption cycles for model A exhibit a large hysteresis loop having a symmetrical shape, i.e., with parallel adsorption and desorption branches. The features of the adsorption isotherms for model B strongly depart from those observed for model A. Both the Ar and Xe adsorption branches for model B correspond to a quasicontinuous pore filling that involves coexistence within the pore of liquid bridges and gas nanobubbles. As in the case of model A, the Ar adsorption isotherm for model B exhibits a significant hysteresis loop; however, the shape of the loop is asymmetrical with a desorption branch much steeper than the adsorption branch. In contrast, the adsorption/desorption cycle for Xe in model B is quasicontinuous and quasireversible. Comparison with adsorption and neutron scattering experiments suggests that model B is too rough at the molecular scale but reproduces reasonably the surface disorder of real MCM-41 at larger length scales. In contrast, model A is smooth at small length scales in agreement with experiments but seems to be too ordered at larger length scales.     

Analysis of van der Waals loops in the adsorption isotherms of porous systems

The properties of the van der Waals loops for isotherms of the adsorbate located in the slit-shaped and cylindrical pores of different width were analyzed. The adsorbate molecules are modeled by spherical Lennard-Jones type particles. The calculation was based on the latticegas model in the quasichemical approximation for intermolecular interactions. The accuracy of calculations of the adsorption hysteresis loops increases due to the use of distributed models, which reflect the spatial inhomogeneity of the distribution of molecules along the normal to the pore wall. The effect of the adsorbate-adsorbent interaction potential and pore width on the pattern of adsorption isotherms is considered. Taking into account the inhomogeneity of the spatial distribution of molecules changes the course of the spinodal sections of the adsorption isotherms.     

A grand canonical Monte Carlo study of capillary condensation in mesoporous media: effect of the pore morphology and topology

We study by means of Grand Canonical Monte Carlo simulations the condensation and evaporation of argon at 77 K in nanoporous silica media of different morphology or topology. For each porous material, our results are compared with data obtained for regular cylindrical pores. We show that both the filling and emptying mechanisms are significantly affected by the presence of a constriction. The simulation data for a constricted pore closed at one end reproduces the asymmetrical shape of the hysteresis loop that is observed for many real disordered porous materials. The adsorption process is a quasicontinuous mechanism that corresponds to the filling of the different parts of the porous material, cavity, and constriction. In contrast, the desorption branch for this pore closed at one end is brutal because the evaporation of Ar atoms confined in the largest cavity is triggered by the evaporation of the fluid confined in the constriction (which isolates the cavity from the gas reservoir). This evaporation process conforms to the classical picture of "pore blocking effect" proposed by Everett many years ago. We also simulate Ar adsorption in a disordered porous medium, which mimics a Vycor mesoporous silica glass. The adsorption isotherm for this disordered porous material having both topological and morphological defects presents the same features as that for the constricted pore (quasicontinuous adsorption and steep desorption process). However, the larger degree of disorder of the Vycor surface enhances these main characteristics. Finally, we show that the effect of the disorder, topological and/or morphological, leads to a significant lowering of the capillary condensation pressure compared to that for regular cylindrical nanopores. Also, our results suggest that confined fluids isolated from the bulk reservoir evaporate at a pressure driven by the smallest size of the pore.     

Interpretation of water isotherm hysteresis for an activated charcoal using stochastic pore networks

Water vapor adsorption equilibria on activated carbons typically exhibit hysteresis. The size and shape of the hysteresis loop which separates the adsorption and desorption branches is a strong function of the pore size and interconnectivity of the pores. Neither conventional pore filling models nor statistical thermodynamics approaches provide a means for predicting the extent of hysteresis from only adsorption measurements. This work uses the Kelvin Equation in conjunction with the structural concept of a stochastic pore network to describe measured water isotherms on BPL carbon. Using a pore segment distribution function determined from the adsorption branch, it is shown that totally random assemblies underestimate the extent of hysteresis. It is possible, however, to closely fit the measured BPL-water hysteresis loop using a patchy heterogeneity in which a proportion of the larger pores are preferentially located on the exterior, mid-range pores are concentrated in a sub-surface layer and some large pores form shielded voids behind much smaller pores.Nomenclature p vapor phase partial pressure of sorbate - p _sat saturation vapor pressure of sorbate - R gas constant - r pore radius - T absolute temperature - t adsorbed layer thickness - V _L molar volume of adsorbed phase - surface tension - contact angle     

Lower closure point of adsorption hysteresis in ordered mesoporous silicas

To examine the nature of the lower closure point of adsorption hysteresis in ordered mesoporous silicas, we measured the temperature dependence of the adsorption-desorption isotherm of nitrogen for three kinds of ordered silicas with cagelike pores and three kinds of ordered silicas with cylindrical pores. The lower closure point pressure of nitrogen in the cagelike pores with sufficiently small necks, that is, the cavitation pressure of a confined liquid, did not depend appreciably on the cage size in the temperature region far away from a hysteresis critical temperature (Tch) but its cage-size dependence was noticeable in the vicinity of Tch. The lower closure point in the cylindrical pores depended on the pore size, and its thermal behavior was totally different from that in the cagelike pores. Nevertheless, the hysteresis critical points of nitrogen in the ordered mesoporous silicas, which are defined as a threshold of temperatures (Tch) and pressure above which reversible capillary condensation takes place in a given size and shape of pores, fell on a common line in a temperature-pressure diagram regardless of the pore geometries. We consider this finding as evidence that capillary evaporation in the cylindrical pores follows a cavitation process in the vicinity of Tch in the same way as that in the cagelike pores and also that the low limit of the hysteresis loop that has been long recognized since 1965 is due to the occurrence of a vapor bubble in a stretched metastable liquid confined to the pores with decreasing pressure (cavitation).     

Adsorption of Carbon Tetrachloride by 3.4 nm Pore Diameter Siliceous MCM-41: Isotherms and Neutron Diffraction

Adsorption isotherms of carbon tetrachloride at temperatures between 273 and 323 K have been determined on the pure silica form of MCM-41 of pore diameter ca. 3.4 nm. All isotherms were of Type V, the isotherms at 273, 288 and 303 K showing hysteresis loops, whereas the isotherm at 323 K was completely reversible. Despite the questionable validity of the Kelvin equation when applied to narrow mesopores, changes in the relative pressure positions of capillary condensation and evaporation as a function of the temperature appear to be well described. Neutron diffraction measurements at 200 and 273 K show significant changes in the physical properties of the adsorbed CCl4 in the MCM-41 from those of bulk adsorbate. The results also suggest a highly heterogeneous surface and appear to show some flexibility in the pore walls upon pore filling. The conditions required for first order reversible capillary condensation are discussed.     

Cluster-Mediated Water Adsorption on Carbon Nanopores

The adsorption isotherms of water at 303 K and N2 at 77 K on various kinds of porous carbons were compared with each other. The saturated amounts of water adsorbed on carbons almost coincided with amounts of N2 adsorption in micropores. Although carbon aerogel samples have mesopores of the great pore volume, the saturated amount of adsorbed water was close to the micropore volume which is much small than the mesopore volume. These adsorption data on carbon aerogels indicated that the water molecules are not adsorbed in mesopores, but in micropores only. The adsorption isotherms of water on activated carbon having micropores of smaller than 0.7 nm in width had no clear adsorption hysteresis, while the water adsorption isotherms on micropores of greater than 0.7 nm had a remarkable adsorption hysteresis above P/P0 = 0.5. The disappearance of the clear hysteresis for smaller micropores suggested that the cluster of water molecules of about 0.7 nm in size gives rise to the water adsorption on the hydrophobic micropores; the formation and the structure of clusters of water molecules were associated with the adsorption mechanism. The cluster-mediated pore filling mechanism was proposed with a special relevance to the evidence on the formation of the ordered water molecular assembly in the carbon micropores by in situ X-ray diffraction.     

Adsorption hysteresis in ordered mesoporous silicas

We review some recent progress in experimental studies of the adsorption hysteresis of simple molecules in ordered mesoporous silicas. We show that the nature of the adsorption hysteresis due to capillary condensation can be examined with less ambiguity by measuring the hysteresis loop for the ordered mesoporous silicas with three types of pore geometries (cylindrical, interconnected cylindrical, and interconnected spherical) over a wide temperature range. The adsorption hysteresis arises from the metastability of a confined phase and the temperature at which the hysteresis disappears is lower than the critical temperature of vapor-liquid equilibrium in pores. The hysteresis occurs mainly on the desorption rather than adsorption branch, irrespective of the pore geometries.     

Adsorption properties of ordered mesoporous silicas synthesized in the presence of block copolymer Pluronic F127 under microwave irradiation

Ordered mesoporous silicas (OMSs) were prepared at different temperatures by using tetraethyl orthosilicate (TEOS) as a silica source, poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) triblock copolymer (Pluronic F127) as a structure directing agent and sodium chloride as an additive under acidic conditions and microwave irradiation. The small angle X-ray diffraction patterns of these samples indicate the presence of ordered mesopores, while adsorption studies show that they possess high volumes of pores, bimodal pore size distributions and large pore sizes. There is an interesting change in the hysteresis loop of nitrogen adsorption isotherms with increasing temperature of hydrothermal treatment; a delayed desorption characteristic for cage-like mesostructures is observed for the OMS samples treated at 100 and 120?°C, while the hydrothermal treatment at 140 and 160?°C leads to the samples having hysteresis loops characteristic for channel-like materials.