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.     

A diffusion model for the fluids confined in micropores

The self-diffusion coefficients were calculated by molecular dynamics simulations and the effects of pore width, temperature, and fluid density on diffusion behavior of simple fluid argon and polar fluid water confined in micropores were analyzed and studied. A mathematical model describing diffusion behavior of fluids confined in micropores was proposed from the theories of molecular dynamics and molecular kinematics, and validated on the basis of the simulation results at various conditions. The model indicates that the diffusion coefficient is proportional to the square root of the pore width and to the temperature divided by the density squared. It is applicable to either liquid or gas states and only two parameters are required.     

Phase behavior in thin films of confined colloid-polymer mixtures

Using self-consistent-field and density-functional theories, we first investigate colloidal self-assembly of colloid-polymer films confined between two soft surfaces grafted by polymers. With increasing colloidal concentrations, the film undergoes a series of transitions from disordered liquid --> sparse square --> hexagonal (or mixed square-hexagonal) --> dense square --> cylindrical structures in a plane, which results from the competition between the entropic elasticity of polymer brushes and the steric packing effect of colloidal particles. A phase diagram displays the stable regions of different in-layer ordering structures as the colloidal concentration is varied and layering transitions as the polymer-grafted density is decreased. Our results show a new control mechanism to stabilize the ordering of structures within the films.     

Phase behavior of undecane-tetradecane mixtures confined in SBA-15

Phase behavior of undecane-tetradecane (n-C11H24-C14H30, C11-C14) mixtures in bulk and confined in SBA-15 have been studied using differential scanning calorimetry. The bulk C11-C14 system shows multiple phase regions due to rotator phase. Confined in the pores of SBA-15 (pore diameters 3.8-7.8 nm), the mixtures only show a melting boundary of a straight line and a curve, respectively. In SBA-15 (17.2 nm), phase behavior of themixtures has some similarity to that of the bulk. Under confinement, the phase diagrams of the mixtures vary with the pore size, temperature, and compositions.     

Phase behavior of a simple lattice model with a two-scale repulsive interaction

The properties of a simple one-dimensional lattice model with two repulsive ranges are studied in terms of its analytic solution. Its phase behavior is characterized by the presence of a disorder-order-disorder transition (or a fluid-solid-fluid transition in lattice gas language). A similar situation was discussed by Hemmer and Stell [Phys. Rev. Lett. 24, 1284 (1970)] when considering the purely repulsive version of their ramp potential. The melting of the solid phase, when pressure is increased along an isotherm, is a feature common to both models and one of the characteristic features of water.     

Phase behavior of the restricted primitive model of ionic fluids with association in slitlike pores. Density-functional approach

We present results of investigations of condensation of restricted primitive model of electrolyte solutions with association between oppositely charged ions confined to slitlike pores. The associative interaction leads to the formation of ionic pairs. It is accounted for by incorporating the first-order thermodynamic perturbation theory into the free energy functional. In order to elucidate the role of association, the phase diagrams are compared with those obtained by us recently [O. Pizio et al., J. Chem. Phys. 121, 11957 (2004)] for the restricted primitive model. The inclusion of the association into the theory leads to lowering the critical temperature for the fluid confined to pores with uncharged and with charged walls. We have observed that the average fraction of bonded ions is high along the coexistence envelope.     

Association model of fluids. Phase equilibria in sulphur dioxide mixtures

The vapor-liquid and liquid-liquid equilibria of binary mixtures formed by sulphur dioxide with organic components are reproduced well using a new associated-solution model whose association and solvation constants are defined in terms of the modified segment fractions of chemical species. The model shows a good performance in predicting ternary vapor-liquid and liquid-liquid equilibria of sulphur dioxide mixtures from only binary parameters.     

Phase behavior of magnetic nanoparticles dispersions in bulk and confined geometries

Dispersions of magnetic nanoparticles (ferrofluids) have been for a long time taken as examples of dipolar fluids. Theoretical papers conclude now that an important parameter is the ratio of the anisotropic attractions (dipolar ones) to the isotropic attractions (Van der Waals ones). It is confirmed by the recent experimental results concerning the behavior of small magnetic particles in bulk or confined 2D geometries.     

Viscosity of confined inhomogeneous nonequilibrium fluids

We use the nonlocal linear hydrodynamic constitutive model, proposed by Evans and Morriss [Statistical Mechanics of Nonequilibrium Liquids (Academic, London, 1990)], for computing an effective spatially dependent shear viscosity of inhomogeneous nonequilibrium fluids. The model is applied to a simple atomic fluid undergoing planar Poiseuille flow in a confined channel of several atomic diameters width. We compare the spatially dependent viscosity with a local generalization of Newton's law of viscosity and the Navier-Stokes viscosity, both of which are known to suffer extreme inaccuracies for highly inhomogeneous systems. The nonlocal constitutive model calculates effective position dependent viscosities that are free from the notorious singularities experienced by applying the commonly used local constitutive model. It is simple, general, and has widespread applicability in nanofluidics where experimental measurement of position dependent transport coefficients is currently inaccessible. In principle the method can be used to predict approximate flow profiles of any arbitrary inhomogeneous system. We demonstrate this by predicting the flow profile for a simple fluid undergoing planar Couette flow in a confined channel of several atomic diameters width.     

Local linear viscoelasticity of confined fluids

In this paper the authors propose a novel method to study the local linear viscoelasticity of fluids confined between two walls. The method is based on the linear constitutive equation and provides details about the real and imaginary parts of the local complex viscosity. They apply the method to a simple atomic fluid undergoing zero mean oscillatory flow using nonequilibrium molecular dynamics simulations. The method shows that the viscoelastic properties of the fluid exhibit dramatic spatial changes near the wall-fluid boundary due to the high density in this region. It is also shown that the real part of the viscosity converges to the frequency dependent local shear viscosity sufficiently far away from the wall. This also provides valuable information about the transport properties in the fluid, in general. The viscosity is compared with predictions from the local average density model. The two methods disagree in that the local average density model predicts larger viscosity variations near the wall-fluid boundary than what is observed through the method presented here.     

Phase behavior of a single polyethylene chain confined between two adsorption walls

The phase behavior of a single polyethylene chain confined between two adsorption walls is investigated by using molecular dynamics simulations. In the free space, it is confirmed in our calculation that the isolated polymer chain exhibits a disordered coil state at high temperatures, and collapses into a condensed state at low temperatures, that is, the coil‐to‐globule transition, and the finite chain length effects are considered since the critical region depends on chain lengths. When the chain is confined between two attractive walls, however, the equilibrium properties not only depend on the chain length but also depend on the adsorption energy and the confinement. Mainly, we focus on the influence of polymer chain length, confinement, and adsorption interaction on the equilibrium thermodynamic properties of the polyethylene chains. Chain lengths of N = 40, 80, and 120 beads, distances between the two walls of D = 10, 20, 30, 50, and 90 Å, and adsorption energies of w = 1.5, 2.5, 3.5, 6.5, and 8.5 kcal/mol are considered here. By considering the confinement–adsorption interactions, some new folding structures are found, that is, the hairpin structure for short chain of N = 40 beads, and the enhanced hairpin or crystal like structures for long chains of N = 80 and 120 beads. The results obtained in our simulations may provide some insights into the phase behaviors of confined polymers, which can not be obtained by previous studies without considering confinement–adsorption interactions. © 2008 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 46: 370–387, 2008     

Phase behavior and particle size cutoff effects in polydisperse fluids

We report a joint simulation and theoretical study of the liquid-vapor phase behavior of a fluid in which polydispersity in the particle size couples to the strength of the interparticle interactions. Attention is focused on the case in which the particle diameters are distributed according to a fixed Schulz form with degree of polydispersity delta = 14%. The coexistence properties of this model are studied using grand canonical ensemble Monte Carlo simulations and moment free energy calculations. We obtain the cloud and shadow curves as well as the daughter phase density distributions and fractional volumes along selected isothermal dilution lines. In contrast to the case of size-independent interaction [N. B. Wilding et al., J. Chem. Phys. 121, 6887 (2004)], the cloud and shadow curves are found to be well separated, with the critical point lying significantly below the cloud curve maximum. For densities below the critical value, we observe that the phase behavior is highly sensitive to the choice of upper cutoff on the particle size distribution. We elucidate the origins of this effect in terms of extremely pronounced fractionation effects and discuss the likely appearance of new phases in the limit of very large values of the cutoff.     

Phase behavior of short-range square-well model

Various Monte Carlo techniques are used to determine the complete phase diagrams of the square-well model for the attractive ranges lambda = 1.15 and lambda = 1.25. The results for the latter case are in agreement with earlier Monte Carlo simulations for the fluid-fluid coexistence curve and yield new results for the liquidus-solidus lines. Our results for lambda = 1.15 are new. We find that the fluid-fluid critical point is metastable for both cases, with the case lambda = 1.25 being just below the threshold value for metastability. We compare our results with prior studies and with experimental results for the gamma(II)-crystallin.     

Phase behavior of model lipid bilayers

We investigated the phase behavior of double-tail lipids, as a function of temperature, headgroup interaction and tail length. At low values of the head-head repulsion parameter a(hh), the bilayer undergoes with increasing temperature the transitions from the subgel phase L(c) via the flat gel phase L(beta) to the fluid phase L(alpha). For higher values of a(hh), the transition from the L(c) to the L(alpha) phase occurs via the tilted gel phase L(beta)(') and the rippled phase P(beta)('). The occurrence of the L(beta)(') phase depends on tail length. We find that the rippled structure (P(beta)(')) occurs if the headgroups are sufficiently surrounded by water and that the ripple is a coexistence between the L(c) or L(beta)(') phase and the L(alpha) phase. The anomalous swelling, observed at the P(beta)(') --> L(alpha) transition, is not directly related to the rippled phase, but a consequence of conformational changes of the tails.     

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive 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. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.     

Freezing of fluids confined between mica surfaces

Using grand ensemble simulations, we show that octamethyl-cyclo-tetra-siloxane (OMCTS) confined between two mica surfaces can form a variety of frozen phases which undergo solid-solid transitions as a function of the separation between the surfaces. For atomically smooth mica surfaces, the following sequence of transitions 1[triangle up] --> 1[triangle up]b --> 2B --> 2 square --> 2[triangle up] are observed in the one- and two-layered regimes, where n[triangle up], n[square], and nB denote triangular, square, and buckled phases, respectively, with the prefix n denoting the number of confined layers. The presence of potassium on mica is seen to have a strong influence on the degree of order induced in the fluid. The sequence of solid-solid transitions that occurs with the smooth mica surface is no longer observed. When equilibrated with a state point near the liquid-solid transition, a counterintuitive freezing scenario is observed in the presence of potassium. Potassium disrupts in-plane ordering in the fluid in contact with the mica surface, and freezing is observed only in the inner confined layers. The largest mica separations at which frozen phases were observed ranged from separations that could accommodate six to seven fluid layers. The extent of freezing and the square-to-triangular lattice transition was found to be sensitive to the presence of potassium as well as the thermodynamic conditions of the bulk fluid. The implications of our results on interpretation of surface force experiments as well as the generic phase behavior of confined soft spheres is discussed.     

Dielectric permittivity profiles of confined polar fluids

The dielectric response of a simple model of a polar fluid near neutral interfaces is examined by a combination of linear response theory and extensive molecular dynamics simulations. Fluctuation expressions for a local permittivity tensor epsilon(r) are derived for planar and spherical geometries, based on the assumption of a purely local relationship between polarization and electric field. While the longitudinal component of epsilon exhibits strong oscillations on the molecular scale near interfaces, the transverse component becomes ill defined and unphysical, indicating nonlocality in the dielectric response. Both components go over to the correct bulk permittivity beyond a few molecular diameters. Upon approaching interfaces from the bulk, the permittivity tends to increase, rather than decrease as commonly assumed, and this behavior is confirmed for a simple model of water near a hydrophobic surface. An unexpected finding of the present analysis is the formation of "electrostatic double layers" signaled by a dramatic overscreening of an externally applied field inside the polar fluid close to an interface. The local electric field is of opposite sign to the external field and of significantly larger amplitude within the first layer of polar molecules.     

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.     

Phase behavior of self-associating fluids with weaker dispersion interactions between bonded particles

In this study, we explore the global phase behavior of a simple model for self-associating fluids where association reduces the strength of the dispersion interactions between bonded particles. Recent research shows that this type of behavior likely explains the thermodynamic properties of strongly polar fluids and certain micellar solutions. Based on Wertheim's theory of associating liquids [M. S. Wertheim, J. Stat. Phys. 42, 459 (1986); 42, 477 (1986)], our model takes into account the effect that dissimilar particle interactions have on the equilibrium constant for self-association in the system. We find that weaker interactions between bonded molecules tend to favor the dissociation of chains at any temperature and density. This effect stabilizes a monomeric liquid phase at high densities, enriching the global phase behavior of the system. In particular, for systems in which the energy of mixing between bonded and unbonded species is positive, we find a triple point involving a vapor, a dense phase of chain aggregates, and a monomeric liquid. Phase coexistence between the vapor and the monomeric fluid is always more stable at temperatures above the triple point, but a highly associated fluid may exist as a metastable phase under these conditions. The presence of this metastable phase may explain the characteristic nucleation behavior of the liquid phase in strongly dipolar fluids.