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).

Equilibrium adsorption in cylindrical mesopores: A modified Broekhoff and de boer theory versus density functional theory

This paper presents a thermodynamic analysis of capillary condensation phenomena in cylindrical pores. Here, we modified the Broekhoff and de Boer (BdB) model for cylindrical pores accounting for the effect of the pore radius on the potential exerted by the pore walls. The new approach incorporates the recently published standard nitrogen and argon adsorption isotherm on nonporous silica LiChrospher Si-1000. The developed model is tested against the nonlocal density functional theory (NLDFT), and the criterion for this comparison is the condensation/evaporation pressure versus the pore diameter. The quantitative agreement between the NLDFT and the refined version of the BdB theory is ascertained for pores larger than 2 nm. The modified BdB theory was applied to the experimental adsorption branch of adsorption isotherms of a number of MCM-41 samples to determine their pore size distributions (PSDs). It was found that the PSDs determined with the new BdB approach coincide with those determined with the NLDFT (also using the experimental adsorption branch). As opposed to the NLDFT, the modified BdB theory is very simple in its utilization and therefore can be used as a convenient tool to obtain PSDs of all mesoporous solids from the analysis of the adsorption branch of adsorption isotherms of any subcritical fluids.

Modeling of adsorption on nongraphitized carbon surface: GCMC simulation studies and comparison with experimental data

We model nongraphitized carbon black surfaces and investigate adsorption of argon on these surfaces by using the grand canonical Monte Carlo simulation. In this model, the nongraphitized surface is modeled as a stack of graphene layers with some carbon atoms of the top graphene layer being randomly removed. The percentage of the surface carbon atoms being removed and the effective size of the defect ( created by the removal) are the key parameters to characterize the nongraphitized surface. The patterns of adsorption isotherm and isosteric heat are particularly studied, as a function of these surface parameters as well as pressure and temperature. It is shown that the adsorption isotherm shows a steplike behavior on a perfect graphite surface and becomes smoother on nongraphitized surfaces. Regarding the isosteric heat versus loading, we observe for the case of graphitized thermal carbon black the increase of heat in the submonolayer coverage and then a sharp decline in the heat when the second layer is starting to form, beyond which it increases slightly. On the other hand, the isosteric heat versus loading for a highly nongraphitized surface shows a general decline with respect to loading, which is due to the energetic heterogeneity of the surface. It is only when the fluid-fluid interaction is greater than the surface energetic factor that we see a minimum-maximum in the isosteric heat versus loading. These simulation results of isosteric heat agree well with the experimental results of graphitization of Spheron 6 (Polley, M. H.; Schaeffer, W. D.; Smith, W. R. J. Phys. Chem. 1953, 57, 469; Beebe, R. A.; Young, D. M. J. Phys. Chem. 1954, 58, 93). Adsorption isotherms and isosteric heat in pores whose walls have defects are also studied from the simulation, and the pattern of isotherm and isosteric heat could be used to identify the fingerprint of the surface.

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.

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.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

High-pressure adsorption of methane and carbon dioxide on coal

Adsorption isotherms of methane and carbon dioxide on two kinds of Australian coals have been measured at three temperatures up to pressures of 20 MPa. The adsorption behavior is described by three isotherm equations: extended three-parameter, Langmuir, and Toth. Among these, the Toth equation is found to be the most suitable, yielding the most realistic values of pore volume of the coals and the adsorbed phase density. Also, the surface area of coals obtained from CO2 adsorption at 273 K is found to be the meaningful parameter which captures the CO2 adsorption capacity. A maximum in the excess amount adsorbed of each gas appears at a lower pressure with a decrease in temperature. For carbon dioxide, after the appearance of the maximum, an inflection point in the excess amount adsorbed is observed close to the critical density at each temperature, indicating that the decrease in the gas-phase density change with pressure influences the behavior of the excess amount adsorbed. In the context of CO2 sequestration, it is found that CO2 injection pressures of lower than 10 MPa may be desirable for the CH4 recovery process and CO2-holding capacity.

The adsorption of phenolic and organotin compounds by clays and cation exchanged clays

Quaternary ammonium exchanged laponites (Quat-laponites) show selectivity in the adsorption of phenols and chlorinated phenols. Strong adsorbate-adsorbent interactions are indicated by adsorption isotherms. Adsorption of phenols and chlorinated phenols by Quat-smectites is greater than that by the Bi Quat-Smectites prepared in this study. It is thought that the quaternary ammonium exchanged smectite components of the Bi Quat-smectites interact with each other (adsorbent-adsorbent interactions) reducing the number of sites available for adsorbate-adsorbent interactions. Solidification/stabilisation studies of 2-chlorophenol show that a blend of ground granulated blast furnace slag and ordinary Portland cement attenuates 2-chlorophenol more effectively than ordinary Portland cement alone. Tetramethyl ammonium- (TMA-) and tetramethyl phosphonium- (TMP-) montmorillonites were exposed to solutions of phenol or chlorinated phenols. TMP- montmorillonite was the better adsorbent and preferentially adsorbed 4-chlorophenol over phenol. Hydration of the interlayer cations occurs to a greater extent in the TMA-montmorillonite than the TMP-montmorillonite restricting interlayer adsorption. Contrary to that observed for phenols and chlorinated phenols, the Quat-smectites were ineffective as adsorbents for triphenyltin hydroxide and bis(tributyltin) oxide at room temperature. Under microwave conditions, only bis(tributyltin) oxide was adsorbed by the quaternary ammonium exchanged smectites. Bis(tributyltin) oxide was adsorbed from ethanol on the surface of the smectite clays at room temperature and under microwave conditions. The adsorbate-adsorbent interactions were weak. Adsorption is accompanied by conversion of bis(tributyltin) oxide to a different tin(IV) species and the release of sodium cations from the montmorillonite interlayer region. Attempts to introduce conditions suitable for charge transfer interactions between synthesised quaternary ammonium compounds and 2,4,6-trichlorophenol are documented. Transition metal complex exchanged clays adsorb 2,4,6-trichlorophenol and phenol. Strong adsorbate-adsorbent interactions (Type I isotherms) occur when the adsorbate is 2,4,6-trichlorophenol and when the adsorbent is [Fe(bipy)3]2+ exchanged montmorillonite or [Co(bipy)3]3+ exchanged montmorillonite. The 2,2'-bipyridyl ligands of the adsorbents are electron rich and the 2,4,6-trichlorophenol is electron deficient. This may have enhanced adsorbate-adsorbent interactions.

Adsorption of dimethyl ether (DME) on zeolite molecular sieves

In recent years there has been growing interest in the use of dimethyl ether (DME) as an alternative fuel. In this study, the adsorption of DME on molecular sieves 4Å (Mol4A) and 5Å (Mol5A) has been experimentally investigated using the volumetric adsorption method. Data on the adsorption isotherms, heats of adsorption, and adsorption kinetic have been obtained and used to draw conclusions and compare the performance of the two adsorbents. Within the conditions considered, the adsorption capacity of Mol5A was found to be around eight times higher than the capacity of Mol4A. Low temperature adsorption and thermal pre-treatment of the adsorbents in vacuum were observed to be favourable for increased adsorption capacity. The adsorption isotherms for both adsorbent were fitted to the Freundlich model and the corresponding model parameters are proposed. The adsorption kinetic analysis suggest that the DME adsorption on Mol5A is controlled by intracrystalline diffusion resistance, while on Mol4A it is mainly controlled by surface layering resistance with the diffusion only taking place at the start of adsorption and for a very limited short time. The heats of adsorption were calculated by a calorimetric method based on direct temperature measurements inside the adsorption cell. Isosteric heats, calculated by the thermodynamic approach (Clasius-Clapeyron equation), have consistently shown lower values. The maximum heat of adsorption was found to be 25.9kJmol-1 and 20.1kJmol-1 on Mol4A and Mol5A, respectively; thus indicating a physisorption type of interactions. © 2014 Elsevier B.V.

The investigation of photocatalysts and iron based materials in the oxidation and adsorption of toxic organic and chromium materials

The presences of heavy metals, organic contaminants and natural toxins in natural water bodies pose a serious threat to the environment and the health of living organisms. Therefore, there is a critical need to identify sustainable and environmentally friendly water treatment processes. In this dissertation, I focus on the fundamental studies of advanced oxidation processes and magnetic nano-materials as promising new technologies for water treatments. Advanced oxidation processes employ reactive oxygen species (ROS) which can lead to the mineralization of a number of pollutants and toxins. The rates of formation, steady-state concentrations, and kinetic parameters of hydroxyl radical and singlet oxygen produced by various TiO2 photocatalysts under UV or visible irradiations were measured using selective chemical probes. Hydroxyl radical is the dominant ROS, and its generation is dependent on experimental conditions. The optimal condition for generation of hydroxyl radical by of TiO2 coated glass microspheres is studied by response surface methodology, and the optimal conditions are applied for the degradation of dimethyl phthalate. Singlet oxygen (1O2) also plays an important role for advanced processes, so the degradation of microcystin-LR by rose bengal, an 1O2 sensitizer was studied. The measured bimolecular reaction rate constant between MC-LR and 1O2 is ∼ 106 M-1 s-1 based on competition kinetics with furfuryl alcohol. The typical adsorbent needs separation after the treatment, while magnetic iron oxides can be easily removed by a magnetic field. Maghemite and humic acid coated magnetite (HA-Fe3O4) were synthesized, characterized and applied for chromium(VI) removal. The adsorption of chromium(VI) by maghemite and HA-Fe3O4 follow a pseudo-second-order kinetic process. The adsorption of chromium(VI) by maghemite is accurately modeled using adsorption isotherms, and solution pH and presence of humic acid influence adsorption. Humic acid coated magnetite can adsorb and reduce chromium(VI) to non-toxic chromium (III), and the reaction is not highly dependent on solution pH. The functional groups associated with humic acid act as ligands lead to the Cr(III) complex via a coupled reduction-complexation mechanism. Extended X-ray absorption fine structure spectroscopy demonstrates the Cr(III) in the Cr-loaded HA-Fe 3O4 materials has six neighboring oxygen atoms in an octahedral geometry with average bond lengths of 1.98 Å.

Porosity in ion-exchanged and acid activated clays evaluated using n-nonane pre-adsorption

The applicability of the n-nonane pre-adsorption method for characterising the porosity in clays is presented. Na-SD, a Na+-exchanged purified bentonite, and materials obtained by Al3+-exchange and acid treatments of Na-SD and SAz-1 were used. Nitrogen adsorption isotherms, at -196 ºC, were determined before and after n-nonane pre-adsorption on each of the samples. In all materials, n-nonane remained adsorbed in ultramicropores after outgassing at 25 ºC. Outgassing at higher temperatures (50, 75 and 200 ºC) removed nonane and ultramicropores became available for nitrogen adsorption. All treatments on Na-SD led to increase in micropore volume. Larger ultramicropore and supermicropore volumes were obtained for Na-SD acid activated with HCl at 95 ºC than for treatments at 25 ºC with HCl or following Al3+-exchange (Al-SD), and increased with increasing acid concentration to 3 M. Activation with 4 M HCl led to the largest pore volume with contribution from mesopores. However, the specific external surface area was the same as that obtained for Na-SD, Al-SD and for most of the other acid activated samples. Treatments at 95 ºC with 1 M and 6 M HCl promoted increase in specific external surface area. The micropore volumes and specific external surface area for SAz-1 treated with 1 M HCl at 95 ºC were larger than those of Al-SAz-1, but lower than those obtained for corresponding materials derived from Na-SD. The n-nonane pre-adsorption method enabled micropore volumes and specific external surface areas to be obtained for all samples.

Removal of paracetamol on biomass-derived activated carbon: Modeling the fixed bed breakthrough curves using batch adsorption experiments

The remediation of paracetamol (PA), an emerging contaminant frequently found in wastewater treatment plants, has been studied in the low concentration range (0.3–10 mg L−1) using as adsorbent a biomass-derived activated carbon. PA uptake of up to 100 mg g−1 over the activated carbon has been obtained, with the adsorption isotherms being fairly explained by the Langmuir model. The application of Reichemberg and the Vermeulen equations to the batch kinetics experiments allowed estimating homogeneous and heterogeneous diffusion coefficients, reflecting the dependence of diffusion with the surface coverage of PA. A series of rapid small-scale column tests were carried out to determine the breakthrough curves under different operational conditions (temperature, PA concentration, flow rate, bed length). The suitability of the proposed adsorbent for the remediation of PA in fixed-bed adsorption was proven by the high PA adsorption capacity along with the fast adsorption and the reduced height of the mass transfer zone of the columns. We have demonstrated that, thanks to the use of the heterogeneous diffusion coefficient, the proposed mathematical approach for the numerical solution to the mass balance of the column provides a reliable description of the breakthrough profiles and the design parameters, being much more accurate than models based in the classical linear driving force.

Adsorption of pesticides onto activated carbons from wood composites

A series of activated carbon was produced from particleboard and medium-density fibreboard monoliths, which are waste originated from the industry, and then characterized and evaluated for potential application for phenoxyacetic acids removals, such 2,4-dichlorophenoxyacetic acid (2,4-D), 2-methyl-4-chlorophenoxy acetic acid (MCPA) and 3-(3,4-dichlorophenyl)-1,1-dimethylurea (diuron), from the liquid phase. All AC retain the shape of the precursor, and displays a microporous structure well-developed, reaching 0.58 cm 3 g -1. The adsorption isotherms for three pesticides were obtained in the optimal conditions and the AC with high superficial area and micropore volume exhibited better performance, allowing to state that, this AC could be a great substitute of those habitually used for this purpose. The pesticides adsorption data were linearized using the Langmuir and Freundlich equation, being the first a very good fit to the experimental data.