Phase Equilibria of Hydrogen Bonding Fluid in a Slit Pore with Broken Symmetry

Phase equilibria of hydrogen bonding (HB) fluid confined in a slit pore with broken symmetry were investigated by the density functional theory incorporated with modified fundamental measure theory, where the symmetry breaking originated from the distinct interactions between fluid molecules and two walls of the slit pore. In terms of adsorption-desorption isotherms and the corresponding grand potentials, phase diagrams of HB fluid under various conditions are presented. Furthermore, through phase coexistences of laying transition and capillary condensation, the effects of HB interaction, pore width, fluid-pore interaction and the broken symmetry on the phase equilibrium properties are addressed. It is shown that these factors can give rise to apparent influences on the phase equilibria of confined HB fluid because of the competition between intermolecular interaction and fluid-pore interaction. Interestingly, a significant influence of broken symmetry of the slit pore is found, and thus the symmetry breaking can provide a new way to regulate the phase behavior of various confined fluids.     

On the separation of nonadditive symmetric mixtures in nanoscopic slitlike pores: A simple model for racemic fluids

A grand canonical ensemble Monte Carlo simulation method is used to study the adsorption of nonadditive symmetric mixtures of Lennard-Jones spherical particles in nanoscopic slitlike pores. The walls of the pore are assumed to be formed by the parallel (100) planes of the model face centered cubic crystal of adjustable corrugation potential. It is demonstrated that depending on the nonadditivity effects in the mixture and the pore width the condensed phases formed inside the pore may have different structures. In particular, it is shown that the mixture may separate into layers containing only one component each and the stacking may depend on the pore width and properties of the mixture.     

The structure of fluids confined in crystalline slitlike nanoscopic pores: bilayers

Grand canonical and canonical ensemble Monte Carlo simulation methods are used to study the structure and phase behavior of Lennard-Jones fluids confined between the parallel (100) planes of the face centered cubic crystal. Ultra thin slit pores of the width allowing for the formation of only two adsorbate layers are considered. It is demonstrated that the structure of adsorbed phases is very sensitive to the wall-wall separation and to the strength of the fluid-wall potential. It is also shown that the structure of low temperature (solid) phases strongly depends on the fluid density. In particular, when the surface field is sufficiently strong, then the high density phases may exhibit a domain wall structure, quite the same as found in monolayer films adsorbed at a single substrate wall. On the other hand, the weakening of the surface potential leads to the regime in which only the hexagonally ordered bilayer structure is stable. The phase diagrams for a series of systems are estimated. It is shown that, depending on the pore width and the temperature, the condensation leads to the formation of the commensurate or incommensurate phases. The incommensurate phases may have the domain-wall or the hexagonal structure depending on the pore width and the strength of the fluid-wall potential.     

Phase behavior and local structure of a binary mixture in pores: mean-field lattice model calculations for analyzing neutron scattering data

We investigate the phase behavior of an asymmetric binary liquid A-W mixture confined between two planar homogenous substrates (slit pore). Molecules of species W interact preferentially with the solid walls via a long-range potential. Assuming nearest-neighbor attractions between the liquid molecules, we employ a lattice-gas model and a mean-field approximation for the grand potential. Minimization of this potential yields the density profiles of thermodynamically stable phases for fixed temperature, chemical potentials of both species, pore width and strengths of attraction. This model is used to analyze experimental small-angle neutron-scattering (SANS) data on the microscopic structure of the binary system isobutyric acid (iBA)+heavy water (D2O) inside a mesoscopic porous matrix (controlled-pore glass of about 10 nm mean pore width). Confinement-independent model parameters are adjusted so that the theoretical liquid-liquid coexistence curve in the bulk matches its experimental counterpart. By choosing appropriate values of the pore width and the attraction strength between substrates and water we analyze the effect of confinement on the phase diagram. In addition to a depression of the liquid-liquid critical point we observe surface induced phase transitions as well as water-film adsorption near the walls. The temperature dependence of the structure of water-rich and iBA-rich phases of constant composition are discussed in detail. The theoretical predictions are consistent with results of the SANS study and assist their interpretation.     

Phase behavior of a fluid confined in slitlike pores with walls modified by preadsorbed chain molecules

We have studied the microscopic structure, thermodynamics of adsorption, and phase behavior of Lennard-Jones fluid in slitlike pores with walls modified due to preadsorption of chain molecules. The chain species are grafted at the walls by terminating segments. Our theoretical considerations are based on a density functional approach in the semigrand canonical ensemble. The applied constraint refers to the constant number of grafted chain molecules in the pore without restriction of the number of chains at each of the walls. We have observed capillary condensation of Lennard-Jones fluid combined with the change of the distribution of chains from nonsymmetric to symmetric with respect to the pore walls. The phase diagrams of the model are analyzed in detail, dependent on the pore width, length of chains, and grafted density.     

On the structure of bilayer condensed phases confined between crystalline walls of triangular symmetry

Grand canonical ensemble Monte Carlo simulation method is used to study the structure of Lennard-Jones fluids confined between the parallel walls, formed by the (111) planes of the face centered cubic crystal. Thin slit pores with a width allowing for the formation of only two atomic layers are considered. It is shown that the structure of confined solidlike phases is very sensitive to the pore width, the misfit between the size of adsorbate atoms and the size of surface lattice, as well as the corrugation of the surface potential. In particular, when the misfit between the adsorbate atoms and the surface lattice is very small, even a weakly corrugated surface potential highly stabilizes solidlike phases of triangular symmetry. On the other hand, in the case of large misfit the packing effects become a dominating factor and lead to the appearance of solidlike phases of different structures and symmetries.     

Adsorption of Fluids in Pores Formed between Two Hard Cylinders

We study adsorption in pores with curved hard walls that are made of two uniaxial cylinders by using a density functional approach. Two cases are considered: adsorption of hard spheres and adsorption of a Lennard-Jones fluid. In the case of hard spheres, we perform a comparison with the results of grand canonical ensemble Monte Carlo data. This comparison indicates that the applied approach is capable of reproducing the fluid structure quite satisfactorily. For hard spheres, we also make a comparison of the total adsorption effect (expressed as the average density of a confined fluid) inside pores with curved walls with that evaluated for a slitlike pore. We have found that the differences between adsorption in pores with curved walls and in slits with the same wall-to-wall distance are quite low. The calculations for the Lennard-Jones fluid have been concerned with the investigation of the capillary evaporation and with the evaluation of phase diagrams for different pores, including slitlike pores. We have found that the curvature of the pore walls shifts the transition toward lower values of the chemical potential and increases slightly the value of the critical temperature in comparison with the values obtained for a slitlike pore. Copyright 2000 Academic Press.     

Molecular simulation of adsorption and separation of mixtures of short linear alkanes in pillared layered materials at ambient temperature

Grand canonical Monte Carlo and configurational-bias Monte Carlo techniques are carried out to simulate the adsorption of ternary and quaternary mixtures of short linear alkanes, involving methane, ethane, propane, and n-butane, in pillared layered materials at ambient temperature, T=300 K. In the simulation, a pillared layered pore is modeled by a uniform distribution of pillars between two layered walls built by making two separate talc lamellas parallel each other with a given size of interlayer distance. The interaction between fluid molecules and two layered walls is measured by storing potentials calculated in advance at a series of grid points. The interaction between fluid molecules and pillars is also calculated by a site-to-site method. The potential model proposed in this work is proved to be effective because of the simulation result being good agreement with the experimental data for the adsorption of nitrogen at 77 K. Then, the adsorption isotherms of mixtures of short linear alkanes in pillared layered pores with three different porosities psi=0.98, 0.93 and 0.85, and three pore widths H=1.02, 1.70 and 2.38 nm at 300 K are obtained by taking advantage of the model. The simulation results tell us that the longer chain component is preferentially adsorbed at low pressures, and its adsorption increases and then decreases as the pressure increases while the shorter chain component is still adsorbed at high pressures. Moreover, the sorption selectivity of pillared layered materials for the longest chain component in alkane mixtures increases as the mole fraction of methane in the gas phase increases. The selectivity of pillared layered materials for the longest chain component in alkane mixtures also increases as the pore width decreases and the porosity increases.     

Phase equilibria and interfacial tension of fluids confined in narrow pores

Correlation between phase behaviors of a Lennard-Jones fluid in and outside a pore is examined over wide thermodynamic conditions by grand canonical Monte Carlo simulations. A pressure tensor component of the confined fluid, a variable controllable in simulation but usually uncontrollable in experiment, is related with the pressure of a bulk homogeneous system in equilibrium with the confined system. Effects of the pore dimensionality, size, and attractive potential on the correlations between thermodynamic properties of the confined and bulk systems are clarified. A fluid-wall interfacial tension defined as an excess grand potential is evaluated as a function of the pore size. It is found that the tension decreases linearly with the inverse of the pore diameter or width.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

Solvation force, structure and thermodynamics of fluids confined in geometrically rough pores

The effect of periodic surface roughness on the behavior of confined soft sphere fluids is investigated using grand canonical Monte Carlo simulations. Rough pores are constructed by taking the prototypical slit-shaped pore and introducing unidirectional sinusoidal undulations on one wall. For the above geometry our study reveals that the solvation force response can be phase shifted in a controlled manner by varying the amplitude of roughness. At a fixed amplitude of roughness, a, the solvation force for pores with structured walls was relatively insensitive to the wavelength of the undulation, lambda for 2.3/=0.5. The predictions of the superposition approximation, where the solvation force response for the rough pores is deduced from the solvation force response of the slit-shaped pores, was in excellent agreement with simulation results for the structured pores and for lambda/sigma(ff)>/=7 in the case of smooth walled pores. Grand potential computations illustrate that interactions between the walls of the pore can alter the pore width corresponding to the thermodynamically stable state, with wall-wall interactions playing an important role at smaller pore widths and higher amplitudes of roughness.     

Adsorption of a binary mixture of adhesive fluids in planar pores: a Monte Carlo study

Grand canonical Monte Carlo simulation is used to study the adsorption of a binary sticky hard-sphere fluid mixture in planar pores. The wall-component 1 and wall-component 2 contact densities are determined to calculate the pressure as a function of the composition of the mixture and the separation between the walls. From these data dependence of the solvation force between the plates on pore width is estimated. The simulation results are compared with the predictions of the Percus-Yevick approximation for planar pores. The density profiles of particular components show interesting shapes stemming from the interplay between the steric effects and the competitive adhesion among all possible species pairs. It is shown that narrowing of the pore causes selective partitioning of individual components of the mixture between the bulk phase and the interior of the pore. The agreement between the two methods is better at wider pores and for the component comprised of weakly adhesive particles.     

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

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

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.     

Characteristic features of phase diagrams describing condensation of adsorbate in narrow pores

The phase diagrams describing condensation of adsorbate in micro- and mesoporous adsorbents having slit-shaped and cylindrical pores whose size varied from 1 to 20 monolayers were constructed. The study was performed using the lattice-gas model in the quasichemical approximation to take into account the intermolecular interactions. The phase diagrams for various values of the potential arising from different types of adsorbate--adsorbent interaction were analyzed for adsorption of helium, neon, methane, and carbon tetrachloride in graphite pores. Other adsorption systems are considered and the relationship between the pressure and temperature of adsorbate condensation is discussed. A nonmonotonic variation of the critical densities for pore widths from 3 to 10 molecular diameters was found. The pattern of this variation depends on the ratio of the energy of lateral interactions of the adsorbate molecules to the energy of interaction of the adsorbate molecules with pore walls. The critical temperature decreases monotonically with a decrease in the pore width. The stronger the adsorbate interaction with the pore walls, the greater the decrease in the critical temperature.     

Phase behavior of binary symmetric mixtures in pillared slit-like pores: a density functional approach

Density functional approach is applied to study the phase behavior of symmetric binary Lennard-Jones(12,6) mixtures in pillared slit-like pores. Our focus is in the evaluation of the first-order phase transitions in adsorbed phases and lines delimiting mixed and demixed adsorbed phases. The scenario of phase changes is sensitive to the pore width, to the energy of fluid-solid interaction, the amount, and the length of the pillars. Quantitative trends and qualitative changes of the phase diagrams topology are examined depending on the values of these parameters. The presence of pillars provides additional excluded volume effects, besides the confinement due to the pore walls. The effects of attraction between fluid species and pillars counteract this additional confinement. We have observed that both the increasing surface pillar density and the augmenting strength of fluid-solid interactions can qualitatively change the phase diagrams topology for the model with sufficiently strong trends for demixing. If the length of pillars is sufficiently large comparing to the pore width at low temperatures, we observe additional phase transitions of the first and second order due to the symmetry breaking of the distribution of chain segments and fluid species with respect to the slit-like pore center. Re-entrant symmetry changes and additional critical points then are observed.     

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.     

Computer simulation of Lennard-Jones fluid adsorption in wedge-shaped pores with smooth walls

The adsorption isotherms and local density distributions in wedge-shaped pores with smooth solid walls are obtained for a molecular fluid by the grand canonical Monte Carlo and multiple histogram reweighting methods. The effect of wedge opening angle on the phase behavior of the fluid is considered, and the peculiarities of the mechanism of wedge-shaped pore filling are established.     

Density functional theory of adsorption in pillared slit-like pores

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in pillared slit like pores. Specifically, the pillars are built of chains that are bonded by their ends to the opposite pore walls. The approach we propose combines theory of quenched-annealed systems and theory of nonuniform fluids involving chain molecules. We compare the results of theoretical predictions with grand canonical ensemble Monte Carlo simulations and compute theoretical capillary condensation phase diagrams for several model systems.