Three-body interactions in complex fluids: virial coefficients from simulation finite-size effects

A simulation technique is described for quantifying the contribution of three-body interactions to the thermodynamical properties of coarse-grained representations of complex fluids. The method is based on comparing the third virial coefficient $B_3$ for a complex fluid with that of an approximate coarse-grained model described by a pair potential. To obtain $B_3$ we introduce a new technique which expresses its value in terms of the measured volume-dependent asymptote of a certain structural function. The strategy is applicable to both Molecular Dynamics and Monte Carlo simulation. Its utility is illustrated via measurements of three-body effects in models of star polymer and highly size-asymmetrical colloid-polymer mixtures.Comment: 13 pages, 8 figure

Liquid-vapor interface of a polydisperse fluid

We report a Grand Canonical Monte Carlo simulation study of the liquid-vapor interface of a model fluid exhibiting polydispersity in terms of the particle size $\sigma$. The bulk density distribution, $\rho^0(\sigma)$, of the system is controlled by the imposed chemical potential distribution $\mu(\sigma)$. We choose the latter such that $\rho^0(\sigma)$ assumes a Schulz form with associated degree of polydispersity $\approx 14$. By introducing a smooth attractive wall, a planar liquid-vapor interface is formed for bulk state points within the region of liquid-vapor coexistence. Owing to fractionation, the pure liquid phase is enriched in large particles, with respect to the coexisting vapor. We investigate how the spatial non-uniformity of the density near the liquid-vapor interface affects the evolution of the local distribution of particle sizes between the limiting pure phase forms. We find (as previously predicted by density functional theory, Bellier-Castella {\em et al}, Phys. Rev. {\bf E65}, 021503 (2002)) a segregation of smaller particles to the interface. The magnitude of this effect is quantified for various $\sigma$ via measurements of the relative adsorption. Additionally, we consider the utility of various estimators for the interfacial width and highlight the difficulties of isolating the intrinsic contribution of polydispersity to this width.Comment: 9 pages, 10 Fig

Liquid-gas coexistence and critical point shifts in size-disperse fluids

Specialized Monte Carlo simulations and the moment free energy (MFE) method are employed to study liquid-gas phase equilibria in size-disperse fluids. The investigation is made subject to the constraint of fixed polydispersity, i.e. the form of the `parent' density distribution $\rho^0(\sigma)$ of the particle diameters $\sigma$, is prescribed. This is the experimentally realistic scenario for e.g. colloidal dispersions. The simulations are used to obtain the cloud and shadow curve properties of a Lennard-Jones fluid having diameters distributed according to a Schulz form with a large (40%) degree of polydispersity. Good qualitative accord is found with the results from a MFE method study of a corresponding van der Waals model that incorporates size-dispersity both in the hard core reference and the attractive parts of the free energy. The results show that polydispersity engenders considerable broadening of the coexistence region between the cloud curves. The principal effect of fractionation in this region is a common overall scaling of the particle sizes and typical inter-particle distances, and we discuss why this effect is rather specific to systems with Schulz diameter distributions. Next, by studying a family of such systems with distributions of various widths, we estimate the dependence of the critical point parameters on $\delta$. In contrast to a previous theoretical prediction, size-dispersity is found to raise the critical temperature above its monodisperse value. Unusually for a polydisperse system, the critical point is found to lie at or very close to the extremum of the coexistence region in all cases. We outline an argument showing that such behaviour will occur whenever size polydispersity affects only the range, rather than the strength of the inter-particle interactions.Comment: 14 pages, 12 figure

Phase behaviour and particle-size cutoff effects in polydisperse fluids

We report a joint simulation and theoretical study of the liquid-vapor phase behaviour of a fluid in which polydispersity in the particle size couples to the strength of the interparticle interactions. Attention is focussed on the case in which the particles 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-{\em independent} interaction strengths (N.B. Wilding, M. Fasolo and P. Sollich, J. Chem. Phys. {\bf 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 behaviour 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.Comment: 12 pages, 15 figure

Refusal of Medical Fees at Ass. Offices.

n/

Phase behaviour of a symmetrical binary fluid mixture

We have investigated the phase behaviour of a symmetrical binary fluid mixture for the situation where the chemical potentials $\mu_1$ and $\mu_2$ of the two species differ. Attention is focused on the set of interparticle interaction strengths for which, when $\mu_1=\mu_2$, the phase diagram exhibits both a liquid-vapor critical point and a tricritical point. The corresponding phase behaviour for the case $\mu_1\ne\mu_2$ is investigated via integral-equation theory calculations within the mean spherical approximation (MSA), and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes of phase behaviour can occur, these being distinguished by the relationship between the critical lines in the full phase diagram in the space of temperature, density, and concentration. We present the detailed form of the phase diagram for both subtypes and compare with the results from GCMC simulations, finding good overall agreement. The scenario via which one subtype evolves into the other, is also studied, revealing interesting features.Comment: 22 pages, 13 figure

Depletion potentials in highly size-asymmetric binary hard-sphere mixtures: Comparison of accurate simulation results with theory

We report a detailed study, using state-of-the-art simulation and theoretical methods, of the depletion potential between a pair of big hard spheres immersed in a reservoir of much smaller hard spheres, the size disparity being measured by the ratio of diameters q=\sigma_s/\sigma_b. Small particles are treated grand canonically, their influence being parameterized in terms of their packing fraction in the reservoir, \eta_s^r. Two specialized Monte Carlo simulation schemes --the geometrical cluster algorithm, and staged particle insertion-- are deployed to obtain accurate depletion potentials for a number of combinations of q\leq 0.1 and \eta_s^r. After applying corrections for simulation finite-size effects, the depletion potentials are compared with the prediction of new density functional theory (DFT) calculations based on the insertion trick using the Rosenfeld functional and several subsequent modifications. While agreement between the DFT and simulation is generally good, significant discrepancies are evident at the largest reservoir packing fraction accessible to our simulation methods, namely \eta_s^r=0.35. These discrepancies are, however, small compared to those between simulation and the much poorer predictions of the Derjaguin approximation at this \eta_s^r. The recently proposed morphometric approximation performs better than Derjaguin but is somewhat poorer than DFT for the size ratios and small sphere packing fractions that we consider. The effective potentials from simulation, DFT and the morphometric approximation were used to compute the second virial coefficient B_2 as a function of \eta_s^r. Comparison of the results enables an assessment of the extent to which DFT can be expected to correctly predict the propensity towards fluid fluid phase separation in additive binary hard sphere mixtures with q\leq 0.1.Comment: 16 pages, 9 figures, revised treatment of morphometric approximation and reordered some materia

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.

Monte Carlo cluster algorithm for fluid phase transitions in highly size-asymmetrical binary mixtures

Highly size-asymmetrical fluid mixtures arise in a variety of physical contexts, notably in suspensions of colloidal particles to which much smaller particles have been added in the form of polymers or nanoparticles. Conventional schemes for simulating models of such systems are hamstrung by the difficulty of relaxing the large species in the presence of the small one. Here we describe how the rejection-free geometrical cluster algorithm (GCA) of Liu and Luijten [Phys. Rev. Lett 92, 035504 (2004)] can be embedded within a restricted Gibbs ensemble to facilitate efficient and accurate studies of fluid phase behavior of highly size-asymmetrical mixtures. After providing a detailed description of the algorithm, we summarize the bespoke analysis techniques of Ashton et al. [J. Chem. Phys. 132, 074111 (2010)] that permit accurate estimates of coexisting densities and critical-point parameters. We apply our methods to study the liquid--vapor phase diagram of a particular mixture of Lennard-Jones particles having a 10:1 size ratio. As the reservoir volume fraction of small particles is increased in the range 0--5%, the critical temperature decreases by approximately 50%, while the critical density drops by some 30%. These trends imply that in our system, adding small particles decreases the net attraction between large particles, a situation that contrasts with hard-sphere mixtures where an attractive depletion force occurs.Comment: 11 pages, 10 figure

Phase behavior of a fluid with competing attractive and repulsive interactions

Fluids in which the interparticle potential has a hard core, is attractive at moderate separations, and repulsive at greater separations are known to exhibit novel phase behavior, including stable inhomogeneous phases. Here we report a joint simulation and theoretical study of such a fluid, focusing on the relationship between the liquid-vapor transition line and any new phases. The phase diagram is studied as a function of the amplitude of the attraction for a certain fixed amplitude of the long ranged repulsion. We find that the effect of the repulsion is to substitute the liquid-vapor critical point and a portion of the associated liquid-vapor transition line, by two first order transitions. One of these transitions separates the vapor from a fluid of spherical liquidlike clusters; the other separates the liquid from a fluid of spherical voids. At low temperature, the two transition lines intersect one another and a vapor-liquid transition line at a triple point. While most integral equation theories are unable to describe the new phase transitions, the Percus Yevick approximation does succeed in capturing the vapor-cluster transition, as well as aspects of the structure of the cluster fluid, in reasonable agreement with the simulation results.Comment: 15 pages, 20 figure