Simulation studies of fluid critical behaviour

We review and discuss recent advances in the simulation of bulk critical phenomena in model fluids. In particular we emphasise the extensions to finite-size scaling theory needed to cope with the lack of symmetry between coexisting fluid phases. The consequences of this asymmetry for simulation measurements of quantities such as the particle density and the heat capacity are pointed out and the relationship to experiment is discussed. A general simulation strategy based on the finite-size scaling theory is described and its utility illustrated via Monte-Carlo studies of the Lennard-Jones fluid and a two-dimensional spin fluid model. Recent applications to critical polymer blends and solutions are also briefly reviewed. Finally we consider the outlook for future simulation work in the field.Comment: 35 pages Revtex, 11 eps figures. Review article to appear in J. Phys.: Condens. Matte

Accurate Measurements of the Chemical Potential of Polymeric Systems by Monte-Carlo Simulation

We present a new Monte-Carlo method for estimating the chemical potential of model polymer systems. The method is based upon the gradual insertion of a penetrable `ghost' polymer into the system and is effective for large chain lengths and at high densities. Insertion of the ghost chain is facilitated by use of an expanded ensemble in which weighted transitions are permitted between states characterising the strength of the excluded volume and thermal interactions experienced by the ghost chain. We discuss the implementation and optimisation of the method within the framework of the bond fluctuation model, and demonstrate its precision by a calculation of the finite-size corrections to the chemical potential.Comment: 12 pages Latex, Report Number #IP-94.12

Liquid-vapour asymmetry in pure fluids: A Monte Carlo simulation study

Monte Carlo simulations within the grand canonical ensemble are used to obtain the joint distribution of density and energy fluctuations $p_L(\rho,u)$ for two model fluids: a decorated lattice gas and a polymer system. In the near critical region the form of $p_L(\rho,u)$ is analysed using a mixed field finite-size-scaling theory that takes account of liquid-vapour asymmetry. Field mixing transformations are performed that map $p_L(\rho,u)$ onto the joint distribution of critical scaling operators \ptMEstar\ appropriate to the Ising fixed point. Carrying out this procedure permits a very accurate determination of the critical point parameters. By forming various projections of \ptMEstar , the full universal finite-size spectrum of the critical density and energy distributions of fluids is also obtained. In the sub-critical coexistence region, an examination is made of the influence of field mixing on the asymmetry of the density distribution.Comment: 19 pages Latex, 15 Figures available on request. Report Number #IP-94.15

Errors in Monte Carlo simulations using shift register random number generators

We report large systematic errors in Monte Carlo simulations of the tricritical Blume-Capel model using single spin Metropolis updating. The error, manifest as a $20\%$ asymmetry in the magnetisation distribution, is traced to the interplay between strong triplet correlations in the shift register random number generator and the large tricritical clusters. The effect of these correlations is visible only when the system volume is a multiple of the random number generator lag parameter. No such effects are observed in related models.Comment: 7 pages Revtex, 4 ps figures (uuencoded). Paper also available from: http://moses.physik.uni-mainz.de/~wilding/home_wilding.htm

A non-equilibrium Monte Carlo approach to potential refinement in inverse problems

The inverse problem for a disordered system involves determining the interparticle interaction parameters consistent with a given set of experimental data. Recently, Rutledge has shown (Phys. Rev. E63, 021111 (2001)) that such problems can be generally expressed in terms of a grand canonical ensemble of polydisperse particles. Within this framework, one identifies a polydisperse attribute (`pseudo-species') $\sigma$ corresponding to some appropriate generalized coordinate of the system to hand. Associated with this attribute is a composition distribution $\bar\rho(\sigma)$ measuring the number of particles of each species. Its form is controlled by a conjugate chemical potential distribution $\mu(\sigma)$ which plays the role of the requisite interparticle interaction potential. Simulation approaches to the inverse problem involve determining the form of $\mu(\sigma)$ for which $\bar\rho(\sigma)$ matches the available experimental data. The difficulty in doing so is that $\mu(\sigma)$ is (in general) an unknown {\em functional} of $\bar\rho(\sigma)$ and must therefore be found by iteration. At high particle densities and for high degrees of polydispersity, strong cross coupling between $\mu(\sigma)$ and $\bar\rho(\sigma)$ renders this process computationally problematic and laborious. Here we describe an efficient and robust {\em non-equilibrium} simulation scheme for finding the equilibrium form of $\mu[\bar\rho(\sigma)]$. The utility of the method is demonstrated by calculating the chemical potential distribution conjugate to a specific log-normal distribution of particle sizes in a polydisperse fluid.Comment: 6 pages, 3 figure

Critical-point finite-size scaling in the microcanonical ensemble

We develop a scaling theory for the finite-size critical behavior of the microcanonical entropy (density of states) of a system with a critically-divergent heat capacity. The link between the microcanonical entropy and the canonical energy distribution is exploited to establish the former, and corroborate its predicted scaling form, in the case of the 3d Ising universality class. We show that the scaling behavior emerges clearly when one accounts for the effects of the negative background constant contribution to the canonical critical specific heat. We show that this same constant plays a significant role in determining the observed differences between the canonical and microcanonical specific heats of systems of finite size, in the critical region.Comment: 27 pages Revtex, 9 figure

A liquid state theory that remains successful in the critical region

A thermodynamically self-consistent Ornstein-Zernike approximation (SCOZA) is applied to a fluid of spherical particles with a pair potential given by a hard-core repulsion and a Yukawa attractive tail $w(r)=-\exp [-z(r-1)]/r$. This potential allows one to take advantage of the known analytical properties of the solution to the Ornstein-Zernike equation for the case in which the direct correlation function outside the repulsive core is given by a linear combination of two Yukawa tails and the radial distribution function $g(r)$ satisfies the exact core condition $g(r)=0$ for $r<1$. The predictions for the thermodynamics, the critical point, and the coexistence curve are compared here to other theories and to simulation results. In order to unambiguously assess the ability of the SCOZA to locate the critical point and the phase boundary of the system, a new set of simulations has also been performed. The method adopted combines Monte Carlo and finite-size scaling techniques and is especially adapted to deal with critical fluctuations and phase separation. It is found that the version of the SCOZA considered here provides very good overall thermodynamics and a remarkably accurate critical point and coexistence curve. For the interaction range considered here, given by $z=1.8$, the critical density and temperature predicted by the theory agree with the simulation results to about 0.6%.Comment: Prepared for the John Barker festschrift issue of Molecular Physics. 22 pages Latex, 6 ps figure

Quantifying density fluctuations in water at a hydrophobic surface: evidence for critical drying

Employing smart Monte Carlo sampling techniques within the grand canonical ensemble, we investigate the properties of water at a model hydrophobic substrate. By reducing the strength of substrate-water attraction we find that fluctuations in the local number density, quantified by a rigorous definition of the local compressibility $\chi(z)$, increase rapidly for distances $z$ within $1$ or $2$ molecular diameters from the substrate as the degree of hydrophobicity, measured by the macroscopic contact angle $\theta$, increases. Our simulations provide evidence for a continuous (critical) drying transition as the substrate-water interaction becomes very weak: $\cos(\theta)\to -1$. We speculate that the existence of such a transition might account for earlier simulation observations of strongly enhanced density fluctuations

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.

Liquid-vapour phase behaviour of a symmetrical binary fluid mixture

Using Monte-Carlo simulation and mean field calculations, we study the liquid-vapour phase diagram of a square well binary fluid mixture as a function of a parameter $\delta$ measuring the relative strength of interactions between particles of dissimilar and similar species. The results reveal a rich variety of liquid-vapour coexistence behaviour as $\delta$ is tuned. Specifically, we uncover critical end point behaviour, a triple point involving a vapour and two liquids of different density, and tricritical behaviour. For a certain range of $\delta$, the mean field calculations also predict a `hidden' (metastable) liquid-vapour binodal.Comment: Revtex, 8 figure