BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nigel B. Wilding

Nigel B. Wilding appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.soft

Catalog footprint

What is connected

27works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

Density depletion and enhanced fluctuations in water near hydrophobic solutes: identifying the underlying physics

We investigate the origin of the density depletion and enhanced density fluctuations that occur in water in the vicinity of an extended hydrophobic solute. We argue that both phenomena are remnants of the critical drying surface phase transition that occurs at liquid-vapor coexistence in the macroscopic planar limit, ie. as the solute radius $R_s\to\infty$. Focusing on the density profile $ρ(r)$ and a sensitive spatial measure of fluctuations, the local compressibility profile $χ(r)$, we develop a scaling theory which expresses the extent of the density depletion and enhancement in compressibility in terms of $R_s$, the strength of solute-water attraction $\varepsilon_s$, and the deviation from liquid-vapor coexistence $δμ$. Testing the predictions against results of classical density functional theory for a simple solvent and Grand Canonical Monte Carlo simulations of a popular water model, we find that the theory provides a firm physical basis for understanding how water behaves at a hydrophobe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Equilibrium phases and domain growth kinetics of calamitic liquid crystals

The anisotropic shape of calamitic LC particles results in distinct energy values when nematogens are placed side-by-side or end-to-end. The energy anisotropy governed by parameter K' has deep consequences on equilibrium & non-equilibrium properties. Using GB model, which shows Nm & low temperature Sm order, we undertake large-scale MC & MD simulations to probe effect of K' on the equilibrium phase diagram & the non-equilibrium domain growth following a quench in the temperature T (coarsening). There are 2 transitions in the model, I->Nm at Tc1 & Nm->Sm at Tc2<Tc1. K' decreases Tc1 significantly, but has relatively little effect on Tc2. Domain growth in Nm phase exhibits the well-known LAC law, L(t)~t^0.5 & evolution is via annihilation of string defects. The system exhibits dynamical scaling that is also robust with respect to K'. We find that Sm phase at quench temperatures T (T>Tc1->T<Tc2) that we consider has SmB order with a hexatic arrangement of the LC molecules in the layers. Coarsening in this phase exhibits a striking two-time-scale scenario: first the LC molecules align & develop orientational order, followed by emergence of characteristic layering along with the hexatic bond-orientational-order within layers. Consequently, the growth follows the LAC law L(t)~t^0.5 at early times & then shows a sharp crossover to a slower growth regime at later times. Our observations strongly suggest L(t)~t^0.25 in this regime. Interestingly, the correlation function shows dynamical scaling in both the regimes & the scaling function is universal. The dynamics is also robust with respect to changes in K', but the smecticity is more pronounced at larger values. Further, the early-time dynamics is governed by string defects, while the late-time evolution is dictated by interfacial defects. We believe this scenario is generic to Sm phase even with other kinds of local order within Sm layers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Measures of Fluctuations for a Liquid Near Critical Drying

We investigate density fluctuations in a liquid close to a solvophobic substrate at which a surface critical drying transition occurs. Using classical density functional theory, we determine three measures of the spatial extent and strength of the fluctuations, i.e., the local compressibility $χ_μ(z)$, the local thermal susceptibility $χ_T(z)$ and the reduced density $χ_*(z)$; $z$ is the distance from the substrate. Whilst the first measure is frequently used, the second and third were introduced very recently by Eckert et.al., Phys. Rev. Lett. 125, 268004 (2020). For state points in the critical drying regime, all three measures, each scaled by its bulk value, exhibit very similar forms and the ratio of $χ_T(z)$ to $χ_μ(z)$, for fixed $z$ in the vapour-liquid transition region, is constant. Using a scaling treatment of surface thermodynamics we show that such behaviour is to be expected on general grounds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase separation and multibody effects in three-dimensional active Brownian particles

Simulation studies of the phase diagram of repulsive active Brownian particles in three dimensions reveal that the region of motility-induced phase separation between a high and low density phase is enclosed by a region of gas-crystal phase separation. Near-critical loci and structural crossovers can additionally be identified in analogy with simple fluids. Motivated by the striking similarity to the behaviour of equilibrium fluids with short-ranged pair-wise attractions, we show that a direct mapping to pair potentials in the dilute limit implies interactions that are insufficiently attractive to engender phase separation. Instead, this is driven by the emergence of multi-body effects associated with particle caging that occurs at sufficiently high number density. We quantify these effects via information-theoretical measures of $n$-body effective interactions extracted from the configurational structure.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Correction of coarse-graining errors by a two-level method: application to the Asakura-Oosawa model

We present a method that exploits self-consistent simulation of coarse-grained and fine-grained models, in order to analyse properties of physical systems. The method uses the coarse-grained model to obtain a first estimate of the quantity of interest, before computing a correction by analysing properties of the fine system. We illustrate the method by applying it to the Asakura-Oosawa (AO) model of colloid-polymer mixtures. We show that the liquid-vapour critical point in that system is affected by three-body interactions which are neglected in the corresponding coarse-grained model. We analyse the size of this effect and the nature of the three-body interactions. We also analyse the accuracy of the method, as a function of the associated computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coarse-grained depletion potentials for anisotropic colloids: application to lock-and-key systems

When a colloid is mixed with a depletant such as a non-adsorbing polymer, one observes attractive effective interactions between the colloidal particles. If these particles are anisotropic, analysis of these effective interactions is challenging in general. We present a method for inference of approximate (coarse-grained) effective interaction potentials between such anisotropic particles. Using the example of indented (lock-and-key) colloids, we show how numerical solutions can be used to integrate out the (hard sphere) depletant, leading to a depletion potential that accurately characterises the effective interactions. The accuracy of the method is based on matching of contributions to the second virial coefficient of the colloids. The simplest version of our method yields a piecewise-constant effective potential; we also show how this scheme can be generalised to other functional forms, where appropriate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Critical Drying of Liquids

We report a detailed simulation and classical density functional theory study of the drying transition in a realistic model fluid at a smooth substrate. This transition (in which the contact angle $θ\to 180^\circ$) is shown to be critical for both short ranged and long-ranged substrate-fluid interaction potentials. In the latter case critical drying occurs at exactly zero attractive substrate strength. This observation permits the accurate elucidation of the character of the transition via a finite-size scaling analysis of the density probability function. We find that the critical exponent $ν_\parallel$ that controls the parallel correlation length, i.e. the extent of vapor bubbles at the wall, is over twice as large as predicted by mean field and renormalization group calculations. We suggest a reason for the discrepancy. Our findings shed new light on fluctuation phenomena in fluids near hydrophobic and solvophobic interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Improved grand canonical sampling of vapour-liquid transitions

Simulation within the grand canonical ensemble is the method of choice for accurate studies of first order vapour-liquid phase transitions in model fluids. Such simulations typically employ sampling that is biased with respect to the overall number density in order to overcome the free energy barrier associated with mixed phase states. However, at low temperature and for large system size, this approach suffers a drastic slowing down in sampling efficiency. The culprits are geometrically induced transitions (stemming from the periodic boundary conditions) which involve changes in droplet shape from sphere to cylinder and cylinder to slab. Since the overall number density doesn't discriminate sufficiently between these shapes, it fails as an order parameter for biasing through the transitions. Here we report two approaches to ameliorating these difficulties. The first introduces a droplet shape based order parameter that generates a transition path from vapour to slab states for which spherical and cylindrical droplet are suppressed. The second simply biases with respect to the number density in a tetragonal subvolume of the system. Compared to the standard approach, both methods offer improved sampling, allowing estimates of coexistence parameters and vapor-liquid surface tension for larger system sizes and lower temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 $χ(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 $θ$, increases. Our simulations provide evidence for a continuous (critical) drying transition as the substrate-water interaction becomes very weak: $\cos(θ)\to -1$. We speculate that the existence of such a transition might account for earlier simulation observations of strongly enhanced density fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unusual liquid phases for indented colloids with depletion interactions

We study indented spherical colloids, interacting via depletion forces. These systems exhibit liquid-vapor phase transitions whose properties are determined by a combination of strong "lock-and-key" bonds and weaker non-specific interactions. As the propensity for lock-and-key binding increases, the critical point moves to significantly lower density, and the coexisting phases change their structure. In particular, the liquid phase is "porous", exhibiting large percolating voids. The properties of this system depend strongly on the topological structure of an underlying bond network: we comment on the implications of this fact for the assembly of equilibrium states with controlled porous structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Demixing cascades in cluster crystals

In a cluster crystal, each lattice site is occupied by multiple soft-core particles. As the number density is increased at zero temperature, a `cascade' of isostructural phase transitions can occur between states whose site occupancy differs by unity. For low but finite temperature, each of these transitions terminates in a critical point. Using tailored Monte Carlo simulation techniques we have studied such demixing cascades in systems of soft particles interacting via potentials of the generalized exponential form $u(r)=ε\exp[-(r/σ)^n]$. We have estimated the critical parameters of the first few transitions in the cascade as a function of the softness parameter $n$. The critical temperature and pressure exhibit non-monotonic behaviour as $n$ is varied, although the critical chemical potential remains monotonic. The trends for the pressure and chemical potential are confirmed by cell model calculations at zero temperature. As $n\to 2^+$, all the transitions that we have observed are preempted by melting although we cannot rule out that clustering transitions survive at high density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantifying the effects of neglecting many-body interactions in coarse-grained models of complex fluids

We describe a general simulation scheme for assessing the thermodynamic consequences of neglecting many-body effects in coarse-grained models of complex fluids. The method exploits the fact that the asymptote of a simple-to-measure structural function provides direct estimates of virial coefficients. Comparing the virial coefficients of an atomistically detailed system with those of a coarse-grained version described by pair potentials, permits the role of many-body effects to be quantified. The approach is applied to two models: (i) a size-asymmetrical colloid-polymer mixture, and (ii) a solution of star polymers. In the latter case, coarse-graining to an effective fluid described by pair potentials is found to neglect important aspects of the true behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-assembly and crystallisation of indented colloids at a planar wall

We report experimental and simulation studies of the structure of a monolayer of indented ("lock and key") colloids, on a planar surface. On adding a non-absorbing polymer with prescribed radius and volume fraction, depletion interactions are induced between the colloids, with controlled range and strength. For spherical particles, this leads to crystallisation, but the indented colloids crystallise less easily than spheres, in both simulation and experiment. Nevertheless, simulations show that indented colloids do form plastic (rotator) crystals. We discuss the conditions under which this occurs, and the possibilities of lower-symmetry crystal states. We also comment on the kinetic accessibility of these states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-assembly of colloidal polymers via depletion-mediated lock and key binding

We study the depletion-induced self-assembly of indented colloids. Using state-of-the-art Monte Carlo simulation techniques that treat the depletant particles explicitly, we demonstrate that colloids assemble by a lock-and-key mechanism, leading to colloidal polymerization. The morphology of the chains that are formed depends sensitively on the size of the colloidal indentation, with smaller values additionally permitting chain branching. In contrast to the case of spheres with attractive patches, Wertheim's thermodynamic perturbation theory fails to provide a fully quantitative description of the polymerization transition. We trace this failure to a neglect of packing effects and we introduce a modified theory that accounts better for the shape of the colloids, yielding improved agreement with simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Monte Carlo method for chemical potential determination in single and multiple occupancy crystals

We describe a Monte Carlo scheme which, in a single simulation, yields a measurement of the chemical potential of a crystalline solid. Within the isobaric ensemble, this immediately provides an estimate of the system free energy, with statistical uncertainties that are determined precisely and transparently. An extension to multiple occupancy ("cluster") solids permits the direct determination of the cluster chemical potential and hence the equilibrium conditions. We apply the method to a model exhibiting cluster crystalline phases, where we find evidence for an infinite cascade of critical points terminating coexistence between crystals of differing site occupancies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transitions between imperfectly ordered crystalline structures: A phase switch Monte Carlo study

A model for two-dimensional colloids confined laterally by "structured boundaries" (i.e., ones that impose a periodicity along the slit) is studied by Monte Carlo simulations. When the distance D between the confining walls is reduced at constant particle number from an initial value D_0, for which a crystalline structure commensurate with the imposed periodicity fits, to smaller values, a succession of phase transitions to imperfectly ordered structures occur. These structures have a reduced number of rows parallel to the boundaries (from n to n-1 to n-2 etc.) and are accompanied by an almost periodic strain pattern, due to "soliton staircases" along the boundaries. Since standard simulation studies of such transitions are hampered by huge hysteresis effects, we apply the phase switch Monte Carlo method to estimate the free energy difference between the structures as a function of the misfit between D and D_0, thereby locating where the transitions occurs in equilibrium. For comparison, we also obtain this free energy difference from a thermodynamic integration method: the results agree, but the effort required to obtain the same accuracy as provided by phase switch Monte Carlo would be at least three orders of magnitude larger. We also show for a situation where several "candidate structures" exist for a phase, that phase switch Monte Carlo can clearly distinguish the metastable structures from the stable one. Finally, applying the method in the conjugate statistical ensemble (where the normal pressure conjugate to D is taken as an independent control variable) we show that the standard equivalence between the conjugate ensembles of statistical mechanics is violated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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=σ_s/σ_b. Small particles are treated grand canonically, their influence being parameterized in terms of their packing fraction in the reservoir, η_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 η_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 η_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 η_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 η_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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Polydispersity induced solid-solid transitions in model colloids

Specialized Monte Carlo simulation techniques and moment free energy method calculations, capable of treating fractionation exactly, are deployed to study the crystalline phase behaviour of an assembly of spherical particles described by a top-hat "parent" distribution of particle sizes. An increase in either the overall density or the degree of polydispersity is shown to generate a succession of phase transitions in which the system demixes into an ever greater number of face-centred cubic "daughter" phases. Each of these phases is strongly fractionated: it contains a much narrower distribution of particle sizes than is present in the system overall. Certain of the demixing transitions are found to be nearly continuous, accompanied by fluctuations in local particle size correlated over many lattice spacings. We explore possible factors controlling the stability of the phases and the character of the demixing transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystalline phases of polydisperse spheres

We use specialized Monte Carlo simulation methods and moment free energy calculations to provide conclusive evidence that dense polydisperse spheres at equilibrium demix into coexisting fcc phases, with more phases appearing as the spread of diameters increases. We manage to track up to four coexisting phases. Each of these is fractionated: it contains a narrower distribution of particle sizes than is present in the system overall. We also demonstrate that, surprisingly, demixing transitions can be nearly continuous, accompanied by fluctuations in local particle size correlated over many lattice spacings.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Grand canonical simulation of phase behaviour in highly size-asymmetrical binary fluids

We describe a Monte Carlo scheme for the grand canonical simulation study of fluid phase equilibria in highly size-asymmetrical binary mixtures. The method utilizes an expanded ensemble in which the insertion and deletion of large particles is accomplished {\em gradually} by traversing a series of states in which a large particle interacts only partially with the environment of small particles. Free energy barriers arising from interfacial coexistence states are surmounted with the aid of multicanonical preweighting, the associated weights being determined from the transition matrix. As an illustration, we present results for the liquid-vapour coexistence properties of a Lennard-Jones binary mixture having a 10:1 size ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase behaviour of polydisperse spheres: simulation strategies and an application to the freezing transition

The statistical mechanics of phase transitions in dense systems of polydisperse particles presents distinctive challenges to computer simulation and analytical theory alike. The core difficulty, namely dealing correctly with particle size fractionation between coexisting phases, is set out in the context of a critique of previous simulation work on such systems. Specialized Monte Carlo simulation techniques and moment free energy method calculations, capable of treating fractionation exactly, are then described and deployed to study the fluid-solid transition of an assembly of repulsive spherical particles described by a top-hat "parent" distribution of particle sizes. The cloud curve delineating the solid-fluid coexistence region is mapped as a function of the degree of polydispersity $δ$, and the properties of the incipient "shadow" phases are presented. The coexistence region is found to shift to higher densities as $δ$ increases, but does not exhibit the sharp narrowing predicted by many theories and some simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fluid phase coexistence and critical behaviour from simulations in the restricted Gibbs ensemble

The symmetrical restricted Gibbs ensemble (RGE) is a version of the Gibbs ensemble in which particles are exchanged between two boxes of fixed equal volumes. It has recently come to prominence because -- when combined with specialized algorithms -- it provides for the study of near-coexistence density fluctuations in highly size-asymmetric binary mixtures. Hitherto, however, a detailed framework for extracting accurate estimates of critical point and coexistence curve parameters from RGE density fluctuations has been lacking. Here we address this problem by exploiting an exact link between the RGE density fluctuations and those of the grand canonical ensemble. In the sub-critical region we propose and test a simple method for obtaining accurate estimates of coexistence densities. In the critical region we identify an observable that serves as a finite system size estimator for the critical point parameters, and present a finite-size scaling theory that allows extrapolation to the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

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 $μ_1$ and $μ_2$ of the two species differ. Attention is focused on the set of interparticle interaction strengths for which, when $μ_1=μ_2$, the phase diagram exhibits both a liquid-vapor critical point and a tricritical point. The corresponding phase behaviour for the case $μ_1\neμ_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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Critical end point behaviour in a binary fluid mixture

We consider the liquid-gas phase boundary in a binary fluid mixture near its critical end point. Using general scaling arguments we show that the diameter of the liquid-gas coexistence curve exhibits singular behaviour as the critical end point is approached. This prediction is tested by means of extensive Monte-Carlo simulations of a symmetrical Lennard-Jones binary mixture within the grand canonical ensemble. The simulation results show clear evidence for the proposed singularity, as well as confirming a previously predicted singularity in the coexistence chemical potential [Fisher and Upton, Phys. Rev. Lett. 65, 2402 (1990)]. The results suggest that the observed singularities, particularly that in the coexistence diameter, should also be detectable experimentally.

0 citations0 reviewsNo code linkedOpen access