Gibbs ensemble Monte Carlo simulation of adsorption for model surfactant solution in confined slit pores

An isobaric-isothermal Gibbs ensemble Monte Carlo simulation has been carried out to study the adsorption of a model surfactant/solvent mixture in slit nanopores. The adsorption isotherms, the density distributions, and the configuration snapshots were simulated to illustrate the adsorption and self-assembly behaviors of the surfactant in the confined pores. The adsorption isotherms are stepwise: a two-step curve for the smaller (30 A) pore and a three-step one for the larger (50 A) pore. The adsorption isotherms and the interfacial aggregate structure of the surfactants in the pores with various sizes show a qualitatively consistent performance with the previous experimental observation. The micelle size distributions of the adsorbed surfactant aggregates have been analyzed in order to understand the adsorption mechanism, which suggests that the step rise in the surfactant adsorption is associated with the considerable formation of the micelle aggregates in the confined pores. The effect of the interaction between the pore surface and the surfactant on the adsorption behavior has also been investigated. The simulation results indicate that a change in the interaction can modify the shape of adsorption isotherms. A nonlinear mathematical model was used to represent the multistep adsorption isotherms. A good agreement between the model fitting and the simulation data was obtained for both the amount of adsorption and the jump point concentration.     

Monte Carlo modeling of chiral adsorption on nanostructured chiral surfaces and slit pores

Adsorptive separation of chiral molecules is a powerful technique that has long been used in the chemical and pharmaceutical industries. An important challenge in this field is to design and optimize new adsorbents to provide selective discrimination of enantiomers. In this article, we introduce an off-lattice model of chiral adsorption on nanostructured surfaces and slit pores with the aim of predicting their enantioslective properties. The concept presented here involves finding the optimal chiral pattern of active sites on the pore walls that maximizes the difference between the binding energies of the enantiomers. Our initial effort focuses on chiral molecules that do not have specific interactions with the pore surface. One candidate meeting this requirement is 1,2-dimethylcyclopropane (DMCP), a chiral hydrocarbon whose interaction with a model pore surface was described using the Lennard-Jones potential. To model the adsorption of DMCP, we used the Monte Carlo simulation method. It was demonstrated that the separation of the enantiomers of DMCP is hardly obtainable because of the smoothness of the potential energy surface for molecules physisorbed in the pore. However, the simulated results allowed the identification of key factors that influence the binding of the enantiomers of DMCP to the pore walls with a special distribution of active sites. This information will be useful in future considerations of the adsorption of more complex chiral molecules.     

Efficient simulation of binary adsorption isotherms using transition matrix Monte Carlo

Molecular simulations of binary adsorption in porous materials are a useful complement to experimental studies of mixture adsorption. Most molecular simulations of binary adsorption are performed using grand canonical Monte Carlo (GCMC) to independently examine a range of state points of interest. A disadvantage of this approach is that it only yields information at a discrete set of state points; therefore, if a complete isotherm is required for arbitrary conditions, some type of data fitting or interpolation must be used in combination with the GCMC data. We show that the transition matrix Monte Carlo (TMMC) method of Shen and Errington (Shen, V. K.; Errington, J. R. J. Chem.Phys. 2005, 122, 064508) is well-suited to simulation of binary adsorption in porous materials. At the completion of a TMMC simulation, the adsorption isotherm for all possible bulk phase compositions and pressures is available without data fitting or interpolation. It is also straightforward to use results from TMMC to compute derivatives of the isotherm such as the mixture thermodynamic correction factors, partial differential ln f(i)/partial differential ln c(j), again without data fitting or interpolation. This approach should be useful in contexts where information on the full adsorption isotherm is needed, such as the design of adsorption- or membrane-based separations.     

Monte Carlo simulation of polymer adsorption

We developed and employed the incremental gauge cell method to calculate the chemical potential (and thus free energies) of long, flexible homopolymer chains of Lennard-Jones beads with harmonic bonds. The free energy of these chains was calculated with respect to three external conditions: in the zero-density bulk limit, confined in a spherical pore with hard walls, and confined in a spherical pore with attractive pores, the latter case being an analog of adsorption. Using the incremental gauge cell method, we calculated the incremental chemical potential of free polymer chains before and after the globual-random coil transitions. We also found that chains confined in attractive pores exhibit behaviors typical of low temperature physisorption isotherms, such as layering followed by capillary condensation.     

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.     

Low temperature adsorption of CO2 in carbonaceous wedge pores: a Monte Carlo simulation study

Adsorption - A systematic study of carbon dioxide in wedge pores under subcritical conditions were conducted with Grand Canonical Monte Carlo simulation. The effects of various factors:...     

Grand canonical Monte Carlo simulation of the adsorption isotherms of water molecules on model soot particles

The grand canonical Monte Carlo method is used to simulate the adsorption isotherms of water molecules on different types of model soot particles. The soot particles are modeled by graphite-type layers arranged in an onionlike structure that contains randomly distributed hydrophilic sites, such as OH and COOH groups. The calculated water adsorption isotherm at 298 K exhibits different characteristic shapes depending both on the type and the location of the hydrophilic sites and also on the size of the pores inside the soot particle. The different shapes of the adsorption isotherms result from different ways of water aggregation in or/and around the soot particle. The present results show the very weak influence of the OH sites on the water adsorption process when compared to the COOH sites. The results of these simulations can help in interpreting the experimental isotherms of water adsorbed on aircraft soot.     

Monte Carlo simulation of cavitation in pores with nonwetting defects

We investigate the onset of cavitation in a metastable fluid confined to nanoscale pores with nonwetting defects present. Using grand canonical and gauge cell mesocanonical Monte Carlo simulations, we study the degree of metastability (relative vapor pressure), at which the critical bubble forms in a spherical pore with a circular nonwetting defect. It is shown that an increase of the defect size leads to a transition from homogeneous to heterogeneous nucleation of critical bubbles formed at the defect site. In this case, the desorption process may be initiated at larger relative vapor pressures than those predicted by the theories of homogeneous cavitation.     

Adsorption of argon from sub- to supercritical conditions on graphitized thermal carbon black and in graphitic slit pores: a grand canonical Monte Carlo simulation study

In this paper we consider the adsorption of argon on the surface of graphitized thermal carbon black and in slit pores at temperatures ranging from subcritical to supercritical conditions by the method of grand canonical Monte Carlo simulation. Attention is paid to the variation of the adsorbed density when the temperature crosses the critical point. The behavior of the adsorbed density versus pressure (bulk density) shows interesting behavior at temperatures in the vicinity of and those above the critical point and also at extremely high pressures. Isotherms at temperatures greater than the critical temperature exhibit a clear maximum, and near the critical temperature this maximum is a very sharp spike. Under the supercritical conditions and very high pressure the excess of adsorbed density decreases towards zero value for a graphite surface, while for slit pores negative excess density is possible at extremely high pressures. For imperfect pores (defined as pores that cannot accommodate an integral number of parallel layers under moderate conditions) the pressure at which the excess pore density becomes negative is less than that for perfect pores, and this is due to the packing effect in those imperfect pores. However, at extremely high pressure molecules can be packed in parallel layers once chemical potential is great enough to overcome the repulsions among adsorbed molecules.     

Molecular dynamics simulation of benzene in graphite and amorphous carbon slit pores

It is well known that confining a liquid into a pore strongly alters the liquid behavior. Investigations of the effect of confinement are of great importance for many scientific and technological applications. Here, we present a study of the behavior of benzene confined in carbon slit pores. Two types of pores are considered–graphite and amorphous carbon ones. We show that the effect of different pore structure is of crucial importance for the benzene behavior. © 2013 Wiley Periodicals, Inc.     

Grand canonical Monte Carlo simulation for determination of optimum parameters for adsorption of supercritical methane in pillared layered pores

A grand canonical Monte Carlo (GCMC) method is carried out to determine optimum adsorptive storage pressures of supercritical methane in pillared layered pores. In the simulation, the pillared layered pore is modeled by a uniform distribution of pillars between two solid walls. Methane is described as a spherical Lennard-Jones molecule, and Steele's 10-4-3 potential is used for representing the interaction between the fluid and a layered wall. The site-site interaction is also used for calculating the interaction energy between methane molecules and pillars. An effective potential model that reflects the characteristics of a real pillared layered material is proposed here. In the model, a binary interaction parameter, k(fw), is introduced into the combining rule for the cross-energy parameter for the interaction between the fluid and a layered wall. Based on the experimental results for the Zr-pillared material synthesized and characterized by Boksh, Kikkinides, and Yang, the binary interaction parameter, k(fw), is determined by fitting the simulation results to the experimental adsorption data of nitrogen at 77 K. Then, by taking it as a model of pillared layered material, a series of GCMC simulations have been carried out. The excess adsorption isotherms of methane in a pillared layered pore with three different pore widths and porosities are obtained at three supercritical temperatures T=207.3, 237.0, and 266.6 K. Based on the simulation results at different porosities, various pore widths and different supercritical temperatures, the pillared layered pore with porosity psi=0.94 and pore width hsigma(p)=1.02 nm is recommended as adsorption storage material of supercritical methane. Moreover, the optimum adsorption pressure is determined at a given temperature and a fixed width of the pillared layered pore. For example, at temperature T=207.3 K, the optimum adsorption pressures are 3.1, 3.7, and 4.5 M Pa at H=1.02, 1.70, and 2.38 nm, respectively. In summary, the GCMC method is a useful tool for optimizing adsorption storage of supercritical methane in pillared layered material.     

Molecular dynamic calculation of isotherms of alcohol adsorption in model pores of active carbon

The molecular dynamics method is used to calculate the isotherms of vapor adsorption in a model pore of active carbon taking into account the microheterogeneous structure of the adsorbent. In the calculations, an additional volume with a nonadsorbing surface (bulk phase) is added to the adsorbing volume of a slitlike model pore with a diameter and a width of ≈3 and 0.7 nm, respectively, and the trajectories of molecules are determined. After the dynamic equilibrium is established, the equilibrium numbers of molecules in the adsorption and bulk phases are estimated through averaging over time. Varying the total number of molecules, the necessary number of points is obtained for an isotherm, which, employing the theory of volume filling of micropores, is presented in ordinary coordinates and can be compared with the experimental data.     

Kinetics and isotherms of asphaltene adsorption in narrow pores

New data relating to the kinetics and adsorption isotherms of asphaltene in consolidated sandstone core samples are reported. The data were obtained from the measurements of electrokinetics of consolidated sandstone core samples in asphaltene/toluene solutions and petroleum oils. The numerical reduction in the (negative) zeta potential of the sandstone samples were attributed to the adsorption of positively charged molecules of asphaltenes. The hydrodynamics thickness δ of adsorption of asphaltene were followed by monitoring the pressure increase that occurred as the adsorbed layer restricted the rock pores and applying Poiseuille's equation. The flow rates indicated a plateau of asphaltene adsorption at a pore blocking thickness of about δ/r = 0.3, which was also the point at which the streaming current reached a plateau. After increasing to about 30% of the pore radius, the adsorbed layer thickness δ stopped growing either with time or with concentration of asphaltene in the flowing liquid. Alternative hypotheses involving asphaltene adsorption isotherms have been investigated. A theoretical treatment advanced describing particle adsorption in the same terms as molecular adsorption and the Langmuir isotherm, with the free energy of asphaltene adsorption on the rock surface (modeled on silica) calculated on the basis of van der Waals attraction. Acceptable agreement was obtained with the electrokinetic measurements.     

Model liquid crystalline dimers in a slit pore: Monte Carlo simulation study

ABSTRACT

In this work, we present results from (isobaric–isothermal) Monte Carlo Simulation studies of liquid crystalline dimer systems confined in a slit pore. Liquid crystalline dimer systems of various spacer numbers have been considered. Surface-induced conformational and alignment properties of these systems at different pressures under homeotropic anchoring condition have been investigated. We have used easily manageable coarse grained force fields to model both monomer–monomer and monomer–substrate interaction potentials. According to the simulated result, the anchoring of dimers to the surface and orientation of mesogenic units with respect to the surface normal seem to depend on the spacer number for messogen attractive confinement. Dimers with lower spacer number are able be adsorbed to the surface and most of their mesogens are oriented along the surface normal even at lower pressure. Those with larger spacer number are distributed throughout the volume at lower pressure. In the case of mesogen repulsive confinement, most of the dimers are adsorbed to the surface and most mesogens are randomly oriented at low pressure. As the pressure gets higher, the adsorption and orientability increase depending on the type of confinement and spacer number. As a result, clear submolecular partitioning and smectic A like structure have been identified.     

Polymolecular adsorption and capillary condensation in narrow slit pores

Capillary condensation and polymolecular adsorption in narrow slits has been calculated, where the fields of surface forces overlap one another. The calculations were carried out on the basis of macroscopic theory of dispersion forces and the isotherms of lone adsorption layers at the free surface. It has been shown that under the effect of mutual attraction through a gap, polymolecular adsorption films lose their stability long before their thickness has approached the half-width of a flat slit. This results in hysteresis of the capillary condensation in an ensemble of plane-parallel slits.

In the case of systems having strong adsorbate-adsorbate interaction, there has been detected the existence of the lower limit of sizes of slit pores, wherein the capillary meniscus can coexist with adsorption films. With a slit width smaller than the critical one, the meniscus is likely to form a finite contact angle with “dry” surfaces of a slit. Thus an explanation has been given of the lower limit of the capillary condensation in an ensemble of flat-surface, slit pores. In the case of strong adsorbate-adsorbent interaction, the coexistence of meniscus with adsorption films within the scope of the approach used is possible in slits of any width.

The value of corrections for the surface forces effect to be entered in the calculations of slit pores dimensions has been analyzed on the basis of the capillary condensation data obtained.

In wedge-shaped slits there also exists, besides lower limit the upper limit of capillary hysteresis.     

Branched versus linear alkane adsorption in carbonaceous slit pores

The presence of carbonaceous deposits on the internal surfaces of a spark ignition engine has been linked in the literature to impaired vehicle performance, as manifested by increased knocking, higher fuel consumption, higher emissions and other adverse effects. One of the proposed mechanisms, in which the deposits affect the processes in the engine, is the adsorption and desorption of fuel components in the pores within the deposit. In this article we investigate this mechanism in more detail by considering single component adsorption of normal and branched alkanes in selected model slit pores representing the structure of the deposits. We further extend these studies to the binary mixture of normal heptane and isooctane, corresponding to a primary reference fuel blend. In particular, we show that in larger pores adsorption selectivity towards isooctane is about 1.2 on average throughout the pressure range. However, in the smaller 10 Å pore selectivity towards isooctane can be in excess of three as a result of packing effects. These results are then placed in the context of engine performance issues.     

Monte Carlo simulations of hydrogen adsorption in alkali-doped single-walled carbon nanotubes

Monte Carlo simulations and Widom's test particle insertion method have been used to calculate the solubility coefficients (S) and the adsorption equilibrium constants (K) in single-walled (10,10) armchair carbon nanotubes including single nanotubes, and nanotube bundles with various configurations with and without alkali dopants. The hydrogen adsorption isotherms at room temperature were predicted by following the Langmuir adsorption model using the calculated constants S and K. The simulation results were in good agreement with experimental data as well as the grand canonical Monte Carlo simulation results reported in the literature. The simulations of nanotube bundle configurations suggest that the gravimetric hydrogen adsorption increases with internanotube gap size. It may be attributed to favorable hydrogen-nanotube interactions outside the nanotubes. The effect of alkali doping on hydrogen adsorption was studied by incorporating K+ or Li+ ions into nanotube arrays using a Monte Carlo simulation. The results on hydrogen adsorption isotherms indicate hydrogen adsorption of 3.95 wt% for K-doping, and 4.21 wt% for Li-doping, in reasonable agreement with the experimental results obtained at 100 atm and room temperature.     

Stockmayer流体在活性炭孔中的吸附的分子模拟

应用巨正则系综monteCarlo方法模拟Stockmayer流体[以一氯二氟甲烷(R22)为代表]在活性炭孔中的吸附。模拟中R22分子采用等效Stockmayer势能模型,狭缝碳孔墙采用10-4-3模型。通过模拟得到了最佳孔径,并在最佳孔径下,针对不同的主体压力及活性基团密度,得到了吸附等温线、孔中流体的局部密度分布图和较为直观的孔内流体分子的瞬时构象,分析了吸附等温线的特征及孔内流体的吸附结构,认为在0.0,1.0sites/nm^2的活性基团密度下的碳孔内分别发生物理及化学吸附,并确定了最佳操作压力,为工业设计合适的催化剂提供依据。     

Effect of adsorbate vibrations on adsorption isotherms in slit-like pores

A molecular model of adsorbate melting in slit-like pores [1] is used to calculate adsorption isotherms with regard to the contribution from vibrational motions of the adsorbate. The equations are based on discrete distribution functions (the lattice gas model). Molecular distributions are calculated in a quasi-chemical approximation reflecting the effects of direct correlations of interacting particles using the Lennard-Jones potential. The vibrational motion of molecules is taken into account using a modified quasi-dimer Mie model. It is shown that considering vibrations shifts the adsorption isotherms in the chemical potential-density coordinates to higher vapor pressures in the region of high filling. This is due to the need for additional compression of the vapor to transfer it into the adsorbed state with increased kinetic energy.