Adsorption/desorption hysteresis in inkbottle pores: a density functional theory and Monte Carlo simulation study

The mechanisms of adsorption and desorption in inkbottle-shaped pores are considered for lattice models using grand canonical mean field density functional theory and Monte Carlo simulation. We find that they depend significantly on the particular pore geometry, the nature of the fluid-solid interaction, and the temperature. We find two mechanisms for desorption. One mechanism involves the emptying of the main cavity even as the density of fluid in the necks remains high, a mechanism observed recently in studies of an off-lattice model of an inkbottle. The other is a simultaneous desorption from the entire pore space, behavior that is more closely related to the traditional picture of pore blocking in the inkbottle system.     

Capillary condensation of a binary mixture in slit-like pores

We investigate the capillary condensation of two model fluid mixtures in slit-like pores, which exhibit different demixing properties in the bulk phase. The interactions between adsorbate particles are modeled by using Lennard-Jones (12,6) potentials and the adsorbing potentials are of the Lennard-Jones (9,3) type. The calculations are performed for different pore widths and at different concentrations of the bulk gas, by means of density functional theory. We evaluate the capillary phase diagrams and discuss their dependence on the parameters of the model. Our calculations indicate that a binary mixture confined to a slit-like pore may exhibit rich phase behavior.     

Adsorption of binary mixtures on two-dimensional surfaces: theory and Monte Carlo simulations

The adsorption of binary mixtures on square lattices is studied by combining theoretical modeling and Monte Carlo (MC) simulations in grand canonical ensemble. The adsorption thermodynamics is analyzed through the total and partial isotherms. Two theoretical models have been used in the present study: (i) the first, which we called cluster approximation (CA), is based on exact calculations of configurations on finite cells. An efficient algorithm allows us to calculate the detailed structure of the configuration space for \(m = l \times l \) cells; and (ii) the second is a generalization of the classical quasi-chemical approximation (QCA) in which the adsorbate is a binary mixture of species \(a\) and \(b\) . Adsorbate–adsorbate lateral interactions are incorporated in the context of the two mentioned approximations. Results from CA and QCA are compared with MC simulations. Close agreement between simulated and theoretical data supports the validity of the theoretical models to describe the adsorption of mixed gases on two-dimensional surfaces.     

Kinetic Monte Carlo study of binary diffusion in silicalite

We report a Kinetic Monte Carlo (KMC) study of the diffusion of linear n-hexane (nC6) and 2,2-dimethylbutane (22DMB) mixture in zeolite silicalite. We first investigated the loading dependences of single component self- and corrected diffusivities of nC6 at 300 K. Anisotropic transition rates are implemented to account for the distribution of the molecules within the zeolite framework. Repulsive guest-guest interactions are modeled using the parameter introduced by Reed and Ehrlich (Surf. Sci. 102:588–601, 1981). The results are in good agreement with recent experimental Quasi Elastic Neutron Scattering data of Jobic et al. (J. Phys. Chem. B 110:2195–2201, 2006), although the influence of the adsorption isotherm inflection is not reproduced. The binary diffusion study of nC6/22DMB mixtures was performed by implementing the nC6 transition rates used for the single component study while 22DMB molecules propagate via intersection-intersection hops. This KMC model allows for different saturation capacities and accounts for interactions between molecules by introducing f _ij parameters. Results show the large impact of guest-guest interactions between nC6 and 22DMB on both self- and corrected diffusivities of the two components. Molecule-size effects are found to be predominant near 22DMB saturation capacity. Acceleration/deceleration effects already described in the literature are confirmed.     

Crystal nucleation in binary hard sphere mixtures: a Monte Carlo simulation study

We present calculations of the nucleation barrier during crystallization in binary hard sphere mixtures under moderate degrees of supercooling using Monte Carlo simulations in the isothermal-isobaric semigrand ensemble in conjunction with an umbrella sampling technique. We study both additive and negatively nonadditive binary hard sphere systems. The solid-fluid phase diagrams of such systems show a rich variety of behavior, ranging from simple spindle shapes to the appearance of azeotropes and eutectics to the appearance of substitutionally ordered solid phase compounds. We investigate the effect of these types of phase behavior upon the nucleation barrier and the structure of the critical nucleus. We find that the underlying phase diagram has a significant effect on the mechanism of crystal nucleation. Our calculations indicate that fractionation of the species upon crystallization increases the difficulty of crystallization of fluid mixtures and in the absence of fractionation (azeotropic conditions) the nucleation barrier is comparable to pure fluids. We also calculate the barrier to nucleation of a substitutionally ordered compound solid. In such systems, which also show solid-solid phase separation, we find that the phase that nucleates is the one whose equilibrium composition is closer to the composition of the fluid phase.     

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

Monte Carlo study of interfacial properties of associating fluids

Canonical Monte Carlo Simulations have been performed to calculate liquid-vapor properties of the associating square well and Lennard-Jones fluids with one and two sites. Simulations were carried out by using several values of reduced temperatures and association energies. The orthobaric densities, as well as the surface tension of associating square well fluids, were calculated and compared with those reported previously in literature; a good agreement was found among them. Results of surface tension of two-sites associating Lennard-Jones fluids are presented here for the first time.     

Monte Carlo study of four dimensional binary hard hypersphere mixtures

A multithreaded Monte Carlo code was used to study the properties of binary mixtures of hard hyperspheres in four dimensions. The ratios of the diameters of the hyperspheres examined were 0.4, 0.5, 0.6, and 0.8. Many total densities of the binary mixtures were investigated. The pair correlation functions and the equations of state were determined and compared with other simulation results and theoretical predictions. At lower diameter ratios the pair correlation functions of the mixture agree with the pair correlation function of a one component fluid at an appropriately scaled density. The theoretical results for the equation of state compare well to the Monte Carlo calculations for all but the highest densities studied.     

Surface excess free energy of simple fluids confined in cylindrical pores by isothermal-isobaric Monte Carlo: influence of pore size

Confined fluid properties are mainly determined by interfacial phenomena characterized by surface quantities. Based on a simple model of Lennard-Jones particles confined in a cylindrical pore, this study introduces a grand potential surface quantity to quantify the difference in the thermodynamic pressure between the bulk and the confined fluids. The usual surface tension gamma defined as this grand potential difference for the same chemical potential in both confined and bulk states is generally strongly dependent on both the chemical potential and temperature. It is proposed here to introduce another surface quantity zeta which measures the thermodynamic pressure difference between confined and bulk states for identical densities. It is shown that this quantity is much less dependent on confined fluid density or chemical potential. It is actually constant along the gas-like and liquid-like adsorption/desorption branches for an irreversible isotherm (hysteresis), with a different value for each branch. For reversible supercritical isotherms, zeta is shown to remain constant in the low and high density parts of the isotherm. This independence on chemical potential (or equivalently fluid density) is believed to be of great interest for practical applications when one desires to calculate thermodynamic quantities such as the usual surface tension gamma or the thermodynamic pressure of a confined fluid for any given chemical potential and temperature. Such calculations are required to determine fundamental properties such as metastability or coexistence. The effects of temperature, fluid/substrate interaction strength, and pore size are studied.     

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.     

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.     

Linear chain surfactants at a planar interface: a comparative Monte Carlo study of several lattice models

Linear chain surfactants in a densely packed arrangement (such as alkane chains in lipid monolayers in the “uniform tilt” structures) are described by a crude coarse-grained model where the endgroups grafted on the interface form a regular lattice and the chains are described by the bond fluctuation model with chains containing N = 4 effective monomers only. Square-well interactions between the monomers are studied for both the attractive and repulsive case for three choices of the interaction range. None of these models exhibits a structure with uniform tilt. For attractive interactions the last bond has a strong tendency to fold back thus leading to a very high density close to the interface. Only when an intrachain-potential favoring stiff chain configurations also is included one can obtain configurations with uniform tilt order. Although related models (with much longer chain lengths and lower grafting densities) are very useful for the study of polymer brushes, the present case of very short chains in a high-density state clearly is plagued by various lattice artefacts and it is concluded that for modelling linear chain surfactants one should use an off-lattice model even on a coarse-grained level.     

Adsorption of benzene and propene in zeolite MCM-22: a grand canonical Monte Carlo study

The GCMC (grand canonical Monte Carlo) simulation technique was used to predict the competition adsorption characteristics of benzene and propene in different pore systems of MCM-22. The nine-site model of benzene was used, which proved to be effective and efficient. The zeolite was divided into three adsorption sites following a simulated annealing method. It is found that benzene and propene have the same preferential adsorption site and a similar adsorption order in different sites. Moreover, the pure and mixture isotherms of the three sites are drawn. From the isotherms, we obtained a selectivity reversal of the mixture isotherms of benzene and propene in different sites. It is also noted that the competition adsorption in the three adsorption sites for the two adsorbates can fall into three successive steps and the adsorption order of propene in mixture in these three sites is S3→S1→S2. A new model is presented to predict the benzene and propene adsorption equilibrium in MCM-22. This approach yields better multicomponent equilibrium predictions than ideal adsorbed solution theory (IAST). Isotherms at different mole fraction of benzene in gas phase indicate an advantage to increase the feed radio of benzene and propene. Thus, this work is helpful for a better understanding of the adsorption mechanism of benzene and propene in MCM-22 and hence the relation of the catalytic properties of the zeolite to its structure.     

Adsorption and morphology transition of surfactants on hydrophobic surfaces: a lattice Monte Carlo study

In this work, we first show that there are only five independent interchange parameters in the surfactant-solvent-interface system in Larson's model, and then adsorption and morphology transition of surfactants on hydrophobic surfaces are studied by extensive lattice Monte Carlo simulations. In our simulations, we found that there exist six adsorbed morphologies: (1) premature admicelle, (2) hemisphere, (3) hemisphere-hemicylinder mixture, (4) wormlike hemicylinder, (5) perforated monolayer, and (6) monolayer. The surface morphologies and the amount of adsorption on hydrophobic surfaces are found to be affected obviously by two interchange parameters. One is the attractive interaction between tail groups and surface (chiTS), and the other is the solubility of head groups in bulk (chiHW). Phase diagrams in chiHW versus chiTS planes for surfactants with different hydrophobicities (chiTW) and for surfactants with different molecular structures are determined in this work, from which the transitions of surface morphologies and adsorption behaviors are discussed.     

Adsorption of fluids in slitlike pores containing a small amount of mobile ions

We apply density functional theory to investigate changes in the phase behavior of a fluid caused by the presence of mobile ions inside the pore. The approach has been based on the fundamental measure density functional theory and on the theory of nonuniform electrolytes developed recently by O. Pizio, A. Patrykiejew, S. Soko?owski [J. Chem. Phys. 121 (2005) 11,957]. We have evaluated capillary condensation phase diagrams for pores of different widths and for different concentrations of confined ions. The calculations have demonstrated that the presence of ions leads to lowering the critical temperature and to an increase of the value of the chemical potential at the capillary condensation point.     

Monte Carlo simulation of n-alkane adsorption isotherms in carbon slit pores

The single component adsorption of alkanes in carbon slit pores was studied using configurational-biased grand canonical Monte Carlo simulations. Wide ranges of temperature, pressure, alkane chain length, and slit height were studied to evaluate their effects on adsorption. Adsorption isotherms and density and orientation profiles were calculated. The behavior of long alkanes at high temperatures was found to be similar to short alkanes at lower temperatures. This suggests that the isotherms may be related through the Polanyi potential theory.     

Condensation phenomena in nanopores: a Monte Carlo study

The nonequilibrium dynamics of condensation phenomena in nanopores is studied via Monte Carlo simulations of a lattice-gas model. Hysteretic behavior of the particle density as a function of the density of a reservoir is obtained for various pore geometries in two and three dimensions. The shape of the hysteresis loops depend on the characteristics of the pore geometry. The evaporation of particles from a pore can be fitted to a stretched exponential decay of the particle density. Phase-separation dynamics inside the pore is effectively described by a random walk of the non-wetting phases. Domain evolution is significantly slowed down in the presence of a random wall-particle potential and gives rise to a temperature-dependent growth exponent. A geometric roughness of the pore wall only delays the onset of a pure domain growth.     

Adsorption of random copolymers by a selective layer: Monte Carlo studies

We use scaling arguments and computer simulations to investigate the adsorption of symmetric AB-random copolymers (RC) from a diluted solution onto a selective ABA layer. Depending on the ratio between the layer thickness and the size of excess blobs, d/xi, three regimes of RC adsorption are predicted. For large values of the layer thickness RC adsorption can be understood as adsorption on two selective interfaces where sequences of RC chains form bridges. When the layer thickness is of the order of xi, excess blobs are trapped in the layer and localize the copolymer chain strongly. If the layer thickness is very small a weak adsorption scenario is predicted where large loops are formed outside the layer. Our simulations using the bond fluctuation model are in good agreement with the scaling predictions. We show that chain properties display non-monotonous behavior with respect to the layer thickness with optimal values for d approximately xi. In particular, we discuss simulation results for density profiles, statistics of bridges, loops and tails formed by the adsorbed chains, as well as for the adsorption order parameter and free energy.     

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.