Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study

Because of the increasing interest in studying the phenomenon exhibited by charge-stabilized colloidal suspensions in confining geometry, we present a density functional theory (DFT) for a hard-core multi-Yukawa fluid. The excess Helmholtz free-energy functional is constructed by using the modified fundamental measure theory and Rosenfeld's perturbative method, in which the bulk direct correlation function is obtained from the first-order mean spherical approximation. To validate the established theory, grand canonical ensemble Monte Carlo (GCMC) simulations are carried out to determine the density profiles and surface excesses of multi-Yukawa fluid in a slitlike pore. Comparisons of the theoretical results with the GCMC data suggest that the present DFT gives very accurate density profiles and surface excesses of multi-Yukawa fluid in the slitlike pore as well as the radial distribution functions of the bulk fluid. Both the DFT and the GCMC simulations predict the depletion of the multi-Yukawa fluid near a nonattractive wall, while the mean-field theory fails to describe this depletion in some cases. Because the simple form of the direct correlation function is used, the present DFT is computationally as efficient as the mean-field theory, but reproduces the simulation data much better than the mean-field theory.     

Accumulation of Macromolecules in Pores from Dilute Solutions in Poor Solvent

The phenomenon of spontaneous accumulation of macromolecules in pores from the dilute solutions of flexible polymers in poor solvent is studied using the density functional theory. It is shown that the partitioning of macromolecules between the bulk solution and slitlike pore, on whose walls the segments are predominantly adsorbed, depends on the ability of a semidilute polymer solution to wet the wall of a pore and on its size, as well as on the degree of undersaturation of bulk solution and temperature. Partitioning coefficient increases in a jumpwise manner when the surface first-order phase transition or capillary separation take place. It is revealed that the regularities of capillary separation of polymer solutions in a poor solvent are similar to those observed during the capillary condensation of unsaturated vapor. The applicability of the Kelvin equation to describe the conditions of capillary separation of flexible polymer solutions is analyzed.     

Structure and phase behavior of a confined nanodroplet composed of the flexible chain molecules

A polymer density functional theory has been employed for investigating the structure and phase behaviors of the chain polymer, which is modelled as the tangentially connected sphere chain with an attractive interaction, inside the nanosized pores. The excess free energy of the chain polymer has been approximated as the modified fundamental measure-theory for the hard spheres, the Wertheim's first-order perturbation for the chain connectivity, and the mean-field approximation for the van der Waals contribution. For the value of the chemical potential corresponding to a stable liquid phase in the bulk system and a metastable vapor phase, the flexible chain molecules undergo the liquid-vapor transition as the pore size is reduced; the vapor is the stable phase at small volume, whereas the liquid is the stable phase at large volume. The wide liquid-vapor coexistence curve, which explains the wide range of metastable liquid-vapor states, is observed at low temperature. The increase of temperature and decrease of pore size result in a narrowing of liquid-vapor coexistence curves. The increase of chain length leads to a shift of the liquid-vapor coexistence curve towards lower values of chemical potential. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. The equilibrium capillary phase transition occurs at a higher chemical potential than in the bulk phase.     

Phase behavior of ionic fluids in slitlike pores: a density functional approach for the restricted primitive model

We present a density functional theory of nonuniform ionic fluids. This theory is based on the application of the electrostatic contribution to the free energy functional arising from mean spherical approximation for a bulk restricted primitive model and from the energy route bulk equation of state. In order to employ this functional we define a reference fluid and additional averaged densities, according to the approach introduced by Gillespie, Nonner and Eisenberg [J. Phys.: Condens. Matter 14, 12129 (2002)]. In the case of bulk systems the proposed theory reduces to the mean spherical approximation equation of state, arising from the energy route and thus it predicts the first-order phase transition. We use this theory to investigate the effects of confinement on the liquid-vapor equilibria. Two cases are considered, namely an electrolyte confined to the pore with uncharged walls and with charged walls. The dependence of the capillary evaporation diagrams on the pore width and on the electrostatic potential is determined.     

受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

Structures and adsorption of binary hard-core Yukawa mixtures in a slitlike pore: grand canonical Monte Carlo simulation and density-functional study

The grand canonical ensemble Monte Carlo simulation and density-functional theory are applied to calculate the structures, local mole fractions, and adsorption isotherms of binary hard-core Yukawa mixtures in a slitlike pore as well as the radial distribution functions of bulk mixtures. The excess Helmholtz energy functional is a combination of the modified fundamental measure theory of Yu and Wu [J. Chem. Phys. 117, 10156 (2002)] for the hard-core contribution and a corrected mean-field theory for the attractive contribution. A comparison of the theoretical results with the results from the Monte Carlo simulations shows that the corrected theory improves the density profiles of binary hard-core Yukawa mixtures in the vicinity of contact over the original mean-field theory. Both the present corrected theory and the simulations suggest that depletion and desorption occur at low temperature, and the local segregation can be observed in most cases. For binary mixtures in the hard slitlike pore, the present corrected theory predicts more accurate surface excesses than the original one does, while in the case of the attractive pore, no improvement is found in the prediction of a surface excess of the smaller molecule.     

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.     

Density-functional theory and Monte Carlo simulation for the surface structure and correlation functions of freely jointed Lennard-Jones polymeric fluids

We present a nonlocal density-functional theory of polymeric fluids consisting of freely jointed Lennard-Jones chains with explicit consideration of the segment size, van der Waals attraction, and structural correlations due to chain connectivity. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the short-ranged repulsion and the first-order thermodynamic perturbation theory for chain connectivity. The contribution of the long-ranged attraction to the Helmholtz energy functional is taken into account using a quadratic density expansion with the direct correlation function obtained from the first-order mean-spherical approximation. The numerical performance of the density-functional theory is compared well with the simulation results from this work as well as those from the literature for the segment-level density profiles and correlation functions of Lennard-Jones chains in slit pores, near isolated nanoparticles, or in bulk.     

Casimir amplitudes and capillary condensation of near-critical fluids between parallel plates: renormalized local functional theory

We investigate the critical behavior of a near-critical fluid confined between two parallel plates in contact with a reservoir by calculating the order parameter profile and the Casimir amplitudes (for the force density and for the grand potential). Our results are applicable to one-component fluids and binary mixtures. We assume that the walls absorb one of the fluid components selectively for binary mixtures. We propose a renormalized local functional theory accounting for the fluctuation effects. Analysis is performed in the plane of the temperature T and the order parameter in the reservoir ψ(∞). Our theory is universal if the physical quantities are scaled appropriately. If the component favored by the walls is slightly poor in the reservoir, there appears a line of first-order phase transition of capillary condensation outside the bulk coexistence curve. The excess adsorption changes discontinuously between condensed and noncondensed states at the transition. With increasing T, the transition line ends at a capillary critical point T=T(c) (ca) slightly lower than the bulk critical temperature T(c) for the upper critical solution temperature. The Casimir amplitudes are larger than their critical point values by 10-100 times at off-critical compositions near the capillary condensation line.     

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.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Surface-induced phase transitions in ultrathin films of block copolymers

We study theoretically the lamellar-disorder-lamellar phase transitions of AB diblock and tetrablock copolymers confined in symmetric slitlike pores where the planar surface discriminatingly adsorbs A segments but repels B segments, mimicking the hydrophobic/hydrophilic effects that have been recently utilized for the fabrication of environmentally responsive "smart" materials. The effects of film thickness, polymer volume fraction, and backbone structure on the surface morphology have been investigated using a polymer density-functional theory. The surface-induced phase transition is manifested itself in a discontinuous switch of microdomains or a jump in the surface density dictated by the competition of surface adsorption and self-aggregation of the block copolymers. The surface-induced first-order phase transition is starkly different from the thickness-induced symmetric-asymmetric or horizontal-vertical transitions in thin films of copolymer melts reported earlier.     

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.     

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.     

Density functional description of adsorption in slitlike pores modified with chain molecules: a simple model for pillaredlike materials

We propose a density functional theory to describe adsorption of Lennard-Jones fluid in slitlike pores modified by chain molecules. Specifically, the chains are bonded by their ends to the opposite pore walls, so they can form pillaredlike structure. Two models are studied. In the first model, the nonterminating segments of chains can change their configuration inside the pore upon adsorption of spherical species. In the second model, the chains configuration remains fixed, so that the system is similar to a nonuniform quenched-annealed mixture. We study capillary condensation of fluid species inside such modified pores and compare the results obtained for two models.     

Investigation of excess adsorption, solvation force, and plate-fluid interfacial tension for Lennard-Jones fluid confined in slit pores

The excess Helmholtz free energy functional is formulated in terms of a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)] for a short ranged repulsion and a first-order mean-spherical approximation theory [Y. P. Tang, J. Chem. Phys. 118, 4140 (2003)] for a long ranged attraction. Within the framework of the density functional theory, the density profile, excess adsorption, solvation force, and plate-fluid interfacial tension of a Lennard-Jones fluid confined in slit pores are predicted, and the results agree well with the simulation data. The phase equilibria inside the slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal, and the plate-fluid interfacial tensions at equilibrium states are predicted consequently.     

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.     

Prediction of phase behavior of nanoconfined Lennard-Jones fluids with density functional theory based on the first-order mean spherical approximation

The recently proposed first-order mean spherical approximation (FMSA) [Y. Tang, J. Chem. Phys. 121, 10605 (2004)] for inhomogeneous fluids is extended to study the phase behavior of nanoconfined Lennard-Jones fluids, which is consistent with the phase equilibria calculation of the corresponding bulk fluid. With a combination of fundamental measure theory, FMSA provides Helmholtz free energy and direct correlation function to formulate density functional theory, which implementation is as easy as the mean-field theory. Following previous success in predicting density profiles inside slit pores, this work is focused specially on the vapor-liquid equilibrium of the Lennard-Jones fluids inside these pores. It is found that outside the critical region FMSA predicts well the equilibrium diagram of slit pores with the sizes of 5.0, 7.5, and 10 molecular diameters by comparing with available computer simulation data. As a quantitative method, FMSA can be treated as an extension from its bulk calculation, while the mean-field theory is only qualitative, as its bulk version.     

Phase behavior of polymer/nanoparticle blends near a substrate

We use the recent fluids density functional theory of Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005); J. Chem. Phys. 122, 094506 (2005)] to investigate the phase behavior of athermal polymer/nanoparticle blends near a substrate. The blends are modeled as a mixture of hard spheres and freely jointed hard chains, near a hard wall. There is a first order phase transition present in these blends in which the nanoparticles expel the polymer from the surface to form a monolayer at a certain nanoparticle concentration. The nanoparticle transition density depends on the length of the polymer, the nanoparticle diameter, and the overall bulk density of the system. The phase transition is due to both packing entropy effects related to size asymmetry between the components and to the polymer configurational entropy, justifying the so-called "entropic push" observed in experiments. In addition, a layered state is found at higher densities which resembles that in colloidal crystals, in which the polymer and nanoparticles form alternating discrete layers. We show that this laminar state has nearly the same free energy as the homogeneously mixed fluid in the bulk and is nucleated by the surface.     

球腔中氢键流体的吸附-解吸附转变:表面调控研究

集中讨论了球形微腔表面对腔中氢键流体相态结构的调控机制. 为了揭示微腔表面对腔中氢键流体相平衡的影响, 首先根据吸附-解吸附原理并利用经典流体的密度泛函理论计算了微腔中氢键流体的平衡密度分布, 进而通过吸附-解吸附等温线及巨势等温线绘制出体系的相图. 在此基础上, 重点考察了球腔尺寸、 表面作用强度和作用力程对氢键流体毛细凝聚及层化转变的影响. 结果表明, 这些因素可以有效地调控体系毛细凝聚和层化转变的临界约化温度、 临界密度和相区大小等特征, 从而阐明了表面调控的主要机制. 研究结果为设计相关吸附材料提供了理论参考.