Effects of Confinement on Critical Adsorption: Absence of Critical Depletion for Fluids in Slit Pores

The adsorption of a near-critical fluid confined in a slit pore is investigated by means of density functional theory and by Monte Carlo simulation for a Lennard-Jones fluid. Our work was stimulated by recent experiments for SF_6 adsorbed in a mesoporous glass which showed the striking phenomenon of critical depletion, i.e. the adsorption excess "Gamma" first increases but then decreases very rapidly to negative values as the bulk critical temperature T_c is approached from above along near-critical isochores. By contrast, our density functional and simulation results, for a range of strongly attractive wall-fluid potentials, show Gamma monotonically increasing and eventually saturating as the temperature is lowered towards T_c along both the critical (rho=rho_c) and sub-critical isochores (rho<\rho_c). Such behaviour results from the increasingly slow decay of the density profile away from the walls, into the middle of the slit, as T->T_c. For rho < rho_c we find that in the fluid the effective bulk field, which is negative and which favours desorption, is insufficient to dominate the effects of the surface fields which favour adsorption. We compare this situation with earlier results for the lattice gas model with a constant (negative) bulk field where critical depletion was found. Qualitatively different behaviour of the density profiles and adsorption is found in simulations for intermediate and weakly attractive wall-fluid potentials but in no case do we observe the critical depletion found in experiments. We conclude that the latter cannot be accounted for by a single pore model.Comment: 21 pages Revtex. Submitted to Phys. Rev.

Effects of weak surface fields on the density profiles and adsorption of a confined fluid near bulk criticality

The density profile and Gibbs adsorption of a near-critical fluid confined between two identical planar walls is studied by means of Monte Carlo simulation and by density functional theory for a Lennard-Jones fluid. By reducing the strength of wall-fluid interactions relative to fluid-fluid interactions we observe a crossover from behaviour characteristic of the normal surface universality class, strong critical adsorption, to behaviour characteristic of a 'neutral' wall. The crossover is reminiscent of that which occurs near the ordinary surface transition in Ising films subject to vanishing surface fields. For the 'neutral' wall the density profile, away from the walls, is almost constant throughout the slit capillary and gives rise to an adsorption that is constant along the critical isochore. The same 'neutral' wall yields a line of capillary coexistence that is almost identical to the bulk coexistence line. In the crossover regime we observe features in the density profile similar to those found in the magnetisation profile of the critical Ising film subject to weak surface fields, namely two smooth maxima, located away from the walls, which merge into a single maximum at midpoint as the strength of the wall-fluid interaction is reduced or as the distance between walls is decreased. We discuss similarities and differences between the surface critical behaviour of fluids and of Ising magnets.Comment: 34 pages, 10 figures, submitted to the Journ. Chem. Phy

Accurate simulation estimates of phase behaviour in ternary mixtures with prescribed composition

This paper describes an isobaric semi-grand canonical ensemble Monte Carlo scheme for the accurate study of phase behaviour in ternary fluid mixtures under the experimentally relevant conditions of prescribed pressure, temperature and overall composition. It is shown how to tune the relative chemical potentials of the individual components to target some requisite overall composition and how, in regions of phase coexistence, to extract accurate estimates for the compositions and phase fractions of individual coexisting phases. The method is illustrated by tracking a path through the composition space of a model ternary Lennard-Jones mixture.Comment: 6 pages, 3 figure

Metastable liquid-liquid coexistence and density anomalies in a core-softened fluid

Linearly-sloped or `ramp' potentials belong to a class of core-softened models which possess a liquid-liquid critical point (LLCP) in addition to the usual liquid-gas critical point. Furthermore they exhibit thermodynamic anomalies in the density and compressibility, the nature of which may be akin to those occurring in water. Previous simulation studies of ramp potentials have focused on just one functional form, for which the LLCP is thermodynamically stable. In this work we construct a series of ramp potentials, which interpolate between this previously studied form and a ramp-based approximation to the Lennard-Jones (LJ) potential. By means of Monte Carlo simulation, we locate the LLCP, the first order high density liquid (HDL)-low density liquid (LDL) coexistence line, and the line of density maxima for a selection of potentials in the series. We observe that as the LJ limit is approached, the LLCP becomes metastable with respect to freezing into a hexagonal close packed crystalline solid. The qualitative nature of the phase behaviour in this regime shows a remarkable resemblance to that seen in simulation studies of accurate water models. Specifically, the density of the liquid phase exceeds that of the solid; the gradient of the metastable LDL-HDL line is negative in the pressure (p)-temperature (T) plane; while the line of density maxima in the p-T plane has a shape similar to that seen in water and extends well into the {\em stable} liquid region of the phase diagram. As such, our results lend weight to the `second critical point' hypothesis as an explanation for the anomalous behaviour of water.Comment: 7 pages, 8 figure

Free energies of crystalline solids: a lattice-switch Monte Carlo method

We present a method for the direct evaluation of the difference between the free energies of two crystalline structures, of different symmetry. The method rests on a Monte Carlo procedure which allows one to sample along a path, through atomic-displacement-space, leading from one structure to the other by way of an intervening transformation that switches one set of lattice vectors for another. The configurations of both structures can thus be sampled within a single Monte Carlo process, and the difference between their free energies evaluated directly from the ratio of the measured probabilities of each. The method is used to determine the difference between the free energies of the fcc and hcp crystalline phases of a system of hard spheres.Comment: 5 pages Revtex, 3 figure

Wetting of a symmetrical binary fluid mixture on a wall

We study the wetting behaviour of a symmetrical binary fluid below the demixing temperature at a non-selective attractive wall. Although it demixes in the bulk, a sufficiently thin liquid film remains mixed. On approaching liquid/vapour coexistence, however, the thickness of the liquid film increases and it may demix and then wet the substrate. We show that the wetting properties are determined by an interplay of the two length scales related to the density and the composition fluctuations. The problem is analysed within the framework of a generic two component Ginzburg-Landau functional (appropriate for systems with short-ranged interactions). This functional is minimized both numerically and analytically within a piecewise parabolic potential approximation. A number of novel surface transitions are found, including first order demixing and prewetting, continuous demixing, a tricritical point connecting the two regimes, or a critical end point beyond which the prewetting line separates a strongly and a weakly demixed film. Our results are supported by detailed Monte Carlo simulations of a symmetrical binary Lennard-Jones fluid at an attractive wall.Comment: submitted to Phys. Rev.

Self-trapping at the liquid vapor critical point

Experiments suggest that localization via self-trapping plays a central role in the behavior of equilibrated low mass particles in both liquids and in supercritical fluids. In the latter case, the behavior is dominated by the liquid-vapor critical point which is difficult to probe, both experimentally and theoretically. Here, for the first time, we present the results of path-integral computations of the characteristics of a self-trapped particle at the critical point of a Lennard-Jones fluid for a positive particle-atom scattering length. We investigate the influence of the range of the particle-atom interaction on trapping properties, and the pick-off decay rate for the case where the particle is ortho-positronium.Comment: 12 pages, 3 figures, revtex4 preprin