Adsorption and self-assembly of surfactant/supercritical CO2 systems in confined pores: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation has been carried out to study the adsorption and self-organization for a model surfactant/supercritical CO2 system confined in the slit-shape nanopores with amorphous silica-like surfaces. The solid surfaces were designed to be CO2-philic and CO2-phobic, respectively. For the CO2-philic surface, obviously surface adsorption is observed for the surfactant molecules. The various energy profiles were used to monitor the lengthy dynamics process of the adsorption and self-assembly for surfactant micelles or monomers in the confined spaces. The equilibrium properties, including the morphologies and micelle-size distributions of absorbed surfactants, were evaluated based on the equilibrium trajectory data. The interaction between the surfactant and the surface produces an obvious effect on the dynamics rate of surfactant adsorption and aggregation, as well as the final self-assembly equilibrium structures of the adsorbed surfactants. However, for the CO2-phobic surfaces, there are scarcely adsorption layers of surfactant molecules, meaning that the CO2-phobic surface repels the surfactant molecules. It seems to conclude that the CO2 solvent depletion near the interfaces determines the surface repellence to the surfactant molecules. The effect of the CO2-phobic surface confinement on the surfactant micelle structure in the supercritical CO2 has also been discussed. In summary, this study on the microscopic behaviors of surfactant/Sc-CO2 in confined pores will help to shed light on the surfactant self-assembly from the Sc-CO2 fluid phase onto solid surfaces and nanoporous media.     

Structural characterization of surfactant aggregates adsorbed in cylindrical silica nanopores

The self-assembly of nonionic surfactants in the cylindrical pores of SBA-15 silica with a pore diameter of 8 nm was studied by small-angle neutron scattering (SANS) at different solvent contrasts. The alkyl ethoxylate surfactants C(10)E(5) and C(12)E(5) exhibit strong aggregative adsorption in the pores as indicated by the sigmoidal shape of the adsorption isotherms. The SANS intensity profiles can be represented by a sum of two terms, one accounting for diffuse scattering from surfactant aggregates in the pores and the other for Bragg scattering from the pore lattice of the silica matrix. The Bragg reflections are analyzed with a form factor model in which the radial density profile of the surfactant in the pore is approximated by a two-step function. Diffuse scattering is represented by a Teubner-Strey-type scattering function which indicates a preferred distance between adsorbed surface aggregates in the pores. Our results suggest that adsorption starts with formation of discrete surface aggregates which increase in number and eventually merge to interconnected patches as the plateau value of the adsorption isotherm is approached. A grossly different behavior, viz. formation of micelles as in solution, is found for the maltoside surfactant C(10)G(2), in agreement with the observed weak adsorption of this surfactant in SBA-15.     

Structure of Lennard-Jones fluids confined in square nanoscale channels from density functional theory

The density distribution of Lennard-Jones fluids confined in square nanoscale channels with Lennard-Jones walls has been studied using the nonlocal density functional theory (DFT) based on the Tarazona model. The effect of channel lengths on the density profiles with various chemical potentials was discussed. It was found that there is an apparent layering phenomenon for the confined fluids due to the combining influences of the enhancing solid-fluid interaction and the excluded volume effect. The pronounced density peaks were observed at the corners of square channels due to the strong fluid-solid interactions. The grand canonical ensemble Monte Carlo simulation (GCEMC) was applied to test the nonlocal DFT results. The DFT calculations are in relatively good agreement with the GCEMC simulations. The adsorption isotherms in a series of square channels were evaluated based on the obtained density distributions. The adsorption mechanism within the square pores was investigated. A comparison between the adsorptions of the square pores with those of the corresponding slit-size pores has been given.     

Capillary condensation hysteresis in overlapping spherical pores: a monte carlo simulation study

The mechanisms of hysteretic phase transformations in fluids confined to porous bodies depend on the size and shape of pores, as well as their connectivity. We present a Monte Carlo simulation study of capillary condensation and evaporation cycles in the course of Lennard-Jones fluid adsorption in the system of overlapping spherical pores. This model system mimics pore shape and connectivity in some mesoporous materials obtained by templating cubic surfactant mesophases or colloidal crystals. We show different mechanisms of capillary hysteresis depending on the size of the window between the pores. For the system with a small window, the hysteresis cycle is similar to that in a single spherical pore: capillary condensation takes place upon achieving the limit of stability of adsorption film and evaporation is triggered by cavitation. When the window is large enough, the capillary condensation shifts to a pressure higher than that of the isolated pore, and the possibility for the equilibrium mechanism of desorption is revealed. These finding may have important implications for practical problems of assessment of the pore size distributions in mesoporous materials with cagelike pore networks.     

Adsorption isotherms of nonionic surfactants in SBA-15 measured by micro-column chromatography

In this work, micro-column liquid chromatography has been employed for the study of the adsorption of ethylene glucol monoether (C8E4) and pentaethylene glucol monoether (C10E5) nonionic surfactants into the nanopores of SBA-15 silica at two temperatures. The adsorption process for both the surfactants has been investigated in the range of concentration from very dilute solution to "just above" critical micelle concentration (CMC). The adsorption data for both the surfactants are characterized as typical LC-shaped isotherms with plateau near the CMC. A simple two-step adsorption model has been applied to represent an experimental data. An attempt to estimate the thickness of the surfactant layer in the pores, from the experimental data, has been made.     

模拟吸附在狭缝微孔中的丙烷的微观结构

用巨正则系综MonteCarlo（GCEMC）方法模拟了活性碳孔吸附丙烷时的微观结构．在GCEMC模拟中,非极性丙烷分子采用单点LJ球状分子模型,狭缝活性碳孔墙采用10－4-3势能模型．在温度T＝134．3K下,模拟并观察到了丙烷分子在狭缝活性碳孔中的吸附、脱附以及毛细凝聚现象,得到了吸附等温线和孔中流体的局部密度轮廓图．从分子水平出发,详细分析了吸附、毛细冷凝时孔中流体的微观结构,为认识、理解吸附的微观机理提供了工具与借鉴．     

Importance of molecular shape in the adsorption of nitrogen,carbon dioxide and methane on surfaces and in confined spaces

The importance of shape in the adsorption of nitrogen, carbon dioxide and methane (common molecular probes for solid characterization) on surfaces and in confined spaces is investigated for its effects on the adsorption capacity and isosteric heat. We study the possibility of using an equivalent pseudo-sphere model to describe the potential energy of interaction of these molecular probes. On a flat open surface, we find that the equivalent pseudo-sphere model describes adsorption of these species sufficiently well. However, in the confined space of pores, especially pores that accommodate three layers or less, the pseudo-sphere model describes the adsorption badly because of the geometrical constraint on the molecular packing. It is recommended that to study adsorption properly in small pores, potential models that correctly describe molecular shape should be used. In characterization, pseudo-sphere models are commonly used to generate the kernels (local isotherms) for the determination of pore size distribution which can lead to misleading results. We illustrate this with an example to show that the wrong pore size distribution results if pseudo-sphere kernels are used.     

Effect of temperature on the adsorption of water in porous carbons

We report experimental and simulation studies to investigate the effect of temperature on the adsorption isotherms for water in carbons. Adsorption isotherms are measured by a gravimetric technique in carbon-fiber monoliths at 378 and 423 K and studied by molecular simulation in ideal carbon pores in the temperature range 298-600 K. Experimental adsorption isotherms show a gradual water uptake, as the pressure increases, and narrow adsorption-desorption hysteresis loops. In contrast, simulated adsorption isotherms at room temperature are characterized by negligible uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption-desorption hysteresis loops. As the temperature increases, the relative pressure at which pore filling occurs increases and the size of the hysteresis loop decreases. Experimental adsorption-desorption hysteresis loops are narrower than those from simulation. Discrepancies between simulation and experimental results are attributed to heterogeneities in chemical composition, pore connectivity, and nonuniform pore-size distribution, which are not accounted for in the simulation model. The hysteresis phase diagram for confined water is obtained by recording the pressure-density conditions that bound the simulated hysteresis loop at each temperature. We find that the hysteresis critical temperature, i.e., the lowest temperature at which no hysteresis is detected, can be hundreds of degrees lower than the vapor-liquid critical temperature for bulk model water. The properties of confined water are discussed with the aid of simulation snapshots and by analyzing the structure of the confined fluid.     

Gibbs ensemble Monte Carlo simulation of supercritical CO2 adsorption on NaA and NaX zeolites

Adsorption of supercritical carbon dioxide on two kinds of zeolites with identical chemical composition but different pore structure (NaA and NaX) was studied using the Gibbs ensemble Monte Carlo simulation. The model frameworks for the two zeolites with SiAl ratio being unity have been chosen as the solid structures in the simulation. The adsorption behaviors of supercritical CO2 on the NaA and NaX zeolites, based on the adsorption isotherms and isosteric heats of adsorption, were discussed in detail and were compared with the available experimental results. A good agreement between the simulated and experimental results is obtained for both the adsorbed amount and the bulk phase density. The intermediate configurational snapshots and the radial distribution functions between zeolite and adsorbed CO2 molecules were collected in order to investigate the preferable adsorption locations and the confined structure behavior of CO2. The structure behaviors of the adsorbed CO2 molecules show various performances, as compared with the bulk phase, due to the confined effect in the zeolite pores.     

NMR studies of cooperative effects in adsorption

The conversion of gas adsorption isotherms into pore size distributions generally relies upon the assumption of thermodynamically independent pores. Hence, pore-pore cooperative adsorption effects, which might result in a significantly skewed pore size distribution, are neglected. In this work, cooperative adsorption effects in water adsorption on a real, amorphous, mesoporous silica material have been studied using magnetic resonance imaging (MRI) and pulsed-gradient stimulated-echo (PGSE) NMR techniques. Evidence for advanced adsorption can be seen directly using relaxation time weighted MRI. The number and spatial distributions of pixels containing pores of different sizes filled with condensate have been analyzed. The spatial distribution of filled pores has been found to be highly nonrandom. Pixels containing the largest pores present in the material have been observed to fill in conjunction with pixels containing much smaller pores. PGSE NMR has confirmed the spatially extensive nature of the adsorbed ganglia. Thus, long-range (≥40 μm) cooperative adsorption effects, between larger pores associated with smaller pores, occur within mesoporous materials. The NMR findings have also suggested particular types of pore filling mechanisms occur within the porous solid studied.     

巨正则系综Monte Carlo模拟方法确定活性炭的微孔尺寸

根据299K下甲烷在活性炭中的吸附实验数据,通过调节狭缝微孔的孔宽参数,利用巨正则系综MonteCarlo(GCEMC)方法得到不同孔宽下流体的微观结构以及吸附等温线.比较并拟合模拟结果和实验数据,确定了活性炭微孔的平均孔宽,为下一步求解微孔尺寸分布以及为预测吸附剂在不同温度下吸附不同吸附质分子时的吸附性能提供了基础与指导.模拟中,甲烷分子采用单点Lennard-Jones球型分子模型,活性炭用狭缝孔来近似表征,流体分子与单个狭缝墙的相互作用采用著名的Steele的10－4-3势能模型.模拟表明,此方法为考察介孔材料的微孔分布以及微孔平均孔宽提供了新的思路.     

Nonionic surfactant sorption onto the bacterial cell surface: a multi-interaction isotherm

The adsorption of linear polyoxyethylene (POE) alcohol surfactants of the form CxEy onto the surface of a Sphingomonas sp. has been examined. For this study, the alkyl chain length (x) was fixed at 12 and the POE chain length (y) was varied, with y = 4, 7, 9, 10, and 23 ethylene oxide units. Langmuirian isotherms were observed for C12E4 and C12E23, and more complex isotherms were observed for the three intermediate POE chain length surfactants, with C12E7 and C12E9 exhibiting strong S-shaped isotherms. All isotherms showed plateaus near the critical micelle concentration (CMC) with the plateau decreasing with increasing POE chain length. A simple multi-interaction isotherm is proposed that models the sorption isotherm as the sum of two interactions. The first interaction describes monolayer adsorption, whereas the second interaction describes lateral interactions between sorbed surfactant molecules and the formation of surface aggregates. Varying ratios of these two interactions as a function of POE chain length gives rise to the variety of observed isotherm shapes. Results of the isotherm analysis suggest that lateral interactions dominate for surfactants with low POE chain lengths, and the lateral interactions decrease as the POE chain length is increased.     

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.     

Does water condense in carbon pores?

Using grand canonical Monte Carlo (GCMC) simulations of molecular models, we investigate the nature of water adsorption and desorption in slit pores with graphitelike surfaces. Special emphasis is placed on the question of whether water exhibits capillary condensation (i.e., condensation when the external pressure is below the bulk vapor pressure). Three models of water have been considered. These are the SPC and SPC/E models and a model where the hydrogen bonding is described by tetrahedrally coordinated square-well association sites. The water-carbon interaction was described by the Steele 10-4-3 potential. In addition to determining adsorption/desorption isotherms, we also locate the states where vapor-liquid equilibrium occurs for both the bulk and confined states of the models. We find that for wider pores (widths >1 nm), condensation does not occur in the GCMC simulations until the pressure is higher than the bulk vapor pressure, P0. This is consistent with a physical picture where a lack of hydrogen bonding with the graphite surface destabilizes dense water phases relative to the bulk. For narrow pores where the slit width is comparable to the molecular diameter, strong dispersion interactions with both carbon surfaces can stabilize dense water phases relative to the bulk so that pore condensation can occur for P < P0 in some cases. For the narrowest pores studied--a pore width of 0.6 nm--pore condensation is again shifted to P > P0. The phase-equilibrium calculations indicate vapor-liquid coexistence in the slit pores for P < P0 for all but the narrowest pores. We discuss the implications of our results for interpreting water adsorption/desorption isotherms in porous carbons.     

Adsorption of Pure Nonionic Alkylethoxylated Surfactants down to Low Concentrations at a Silica/Water Interface as Determined Using a HPLC Technique

The adsorption of pure nonionic alkylethoxylated surfactants of the C₁₂E_nseries at silica/water interface has been determined using a very precise HPLC technique. The number of ethoxylated groups was varied from 2 to 9. The adsorption isotherms were constructed with special attention to the very low surface coverage domain. It is shown that at very low concentration, the adsorption amounts are higher as the number of ethoxylated groups increases but the reverse trend is found at higher surfactant concentration and above the critical micelle concentration. It is shown that this behavior is the consequence of the interplay of the primary and secondary adsorption mechanisms depending upon the length of the ethoxylated chain. The maximum adsorption quantities is not a linear function of the number of ethoxylated groups. This and other observations confirm the viewpoint that the behavior of nonionic surfactant aggregates adsorbed at a hydrophilic surface carries many similarities with the properties of this class of nonionic surfactant aggregates in bulk aqueous solutions.     

Porous Texture and Surface Character of Dehydroxylated and Rehydroxylated MCM-41 Mesoporous Silicas-Analysis of Adsorption Isotherms of Nitrogen Gas and Water Vapor

Four samples of MCM-41 mesoporous silicas whose average pore diameters are 2.4, 2.8, 3.2, and 3.6 nm were prepared using sodium orthosilicate and cationic surfactants of [CH(3)(CH(2))(n)N(CH(3))(3)]X (n=11, 13, 15, 17). These four samples were calcined at 1123 K in vacuo to obtain the dehydroxylated samples, which were further rehydroxylated at 298 K to obtain the rehydroxylated samples. The adsorption isotherms of nitrogen gas (77 K) for the 12 MCM-41 mesoporous silicas are of Type IVc, giving no adsorption hysteresis. On the other hand, the first adsorption isotherms of water vapor (298 K) for the dehydroxylated MCM-41 samples are quite different from those of nitrogen gas, giving the remarkable adsorption hysteresis. The second water isotherms for the rehydroxylated MCM-41 samples are of Type IV, showing slight hysteresis. Using the nitrogen isotherms, the relation between the pore size and carbon chain length of the surfactant has been determined, and the effect of dehydroxylation and rehydroxylation on the porous texture has been examined. Using the first and second water isotherms, the adsorption model of physisorbed waters adsorbed on the surface silanol groups has been proposed. From the pore size distribution curves of nitrogen and water, the presence of constrictions in the cylindrical pores has been predicted. Copyright 2000 Academic Press.     

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.     

交联聚苯乙烯型多孔吸附剂的中孔性质研究

采用77K温度下的氮气吸附方法,测定了经悬浮聚合制备的不同交联度的交联聚苯乙烯多孔吸附剂的吸附/脱附等温线.根据BET吸附模型计算了比表面,由吸附量计算了总的孔体积,由孔体积和比表面计算出平均孔径,并依据脱附等温线采用BJH方法计算孔径分布.结果表明,交联度对交联聚苯乙烯多孔吸附剂的孔结构均具有显著影响.随着交联聚苯乙烯多孔吸附剂的交联度升高,其孔径变小,比表面增大,而且低交联度吸附剂的中孔接近圆柱形,高交联吸附剂的中孔形状接近“墨水瓶”形.显然,交联度对孔性质的影响与孔结构在交联聚苯乙烯多孔吸附剂制备和后处理过程中的稳定性密切相关.交联度低时,初期形成的小孔不能保持稳定,在后续聚合及后处理过程中合并为大孔,结果造成低交联吸附剂大孔径、低比表面的现象.通过对孔径分布的研究,揭示了不同吸附剂在中孔范围内的孔特征,并对其形成机制进行了分析.