Biomedical subjects

Nigel B Wilding

Publications and source records attributed to Nigel B Wilding.

13 recordsLinked to original sources

PubMed2006

Phase behavior of a symmetrical binary fluid mixture.

We have investigated the phase behavior 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 behavior for the case mu(1) not equalmu(2) is investigated via integral-equation theory calculations within the mean spherical approximation and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes of phase behavior can occur, these being distinguished by the relationship between the triple 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.

Journal Article↗

PubMed2006

Wetting transitions in polydisperse fluids.

The properties of the coexisting bulk gas and liquid phases of a polydisperse fluid depend not only on the prevailing temperature but also on the overall parent density. As a result, a polydisperse fluid near a wall will exhibit density-driven wetting transitions inside the coexistence region. We propose a likely topology for the wetting phase diagram, which we test using Monte Carlo simulations of a model polydisperse fluid at an attractive wall, tracing the wetting line inside the cloud curve and identifying the relationship to prewetting.

Journal Article↗

PubMed2006

Simulation of phase transitions in highly asymmetric fluid mixtures.

We present a novel method for the accurate numerical determination of the phase behavior of fluid mixtures having large particle-size asymmetries. By incorporating the recently developed geometric cluster algorithm within a restricted Gibbs ensemble, we are able to probe directly the density and concentration fluctuations that drive phase transitions, but that are inaccessible to conventional simulation algorithms. We develop a finite-size scaling theory that relates these density fluctuations to those of the grand-canonical ensemble, thereby enabling accurate location of critical points and coexistence curves of multicomponent fluids. Several illustrative examples are presented.

Journal Article↗

PubMed2006

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.

Journal Article↗

PubMed2006

Simulation estimates of cloud points of polydisperse fluids.

We describe two distinct approaches to obtaining the cloud-point densities and coexistence properties of polydisperse fluid mixtures by Monte Carlo simulation within the grand-canonical ensemble. The first method determines the chemical potential distribution mu(sigma) (with the polydisperse attribute) under the constraint that the ensemble average of the particle density distribution rho(sigma) match a prescribed parent form. Within the region of phase coexistence (delineated by the cloud curve) this leads to a distribution of the fluctuating overall particle density n, p(n), that necessarily has unequal peak weights in order to satisfy a generalized lever rule. A theoretical analysis shows that as a consequence, finite-size corrections to estimates of coexistence properties are power laws in the system size. The second method assigns mu(sigma) such that an equal-peak-weight criterion is satisfied for p(n) for all points within the coexistence region. However, since equal volumes of the coexisting phases cannot satisfy the lever rule for the prescribed parent, their relative contributions must be weighted appropriately when determining mu(sigma). We show how to ascertain the requisite weight factor operationally. A theoretical analysis of the second method suggests that it leads to finite-size corrections to estimates of coexistence properties which are exponentially small in the system size. The scaling predictions for both methods are tested via Monte Carlo simulations of a polydisperse lattice-gas model near its cloud curve, the results showing excellent quantitative agreement with the theory.

Journal Article↗

PubMed2006

Influence of polydispersity on the critical parameters of an effective-potential model for asymmetric hard-sphere mixtures.

We report a Monte Carlo simulation study of the properties of highly asymmetric binary hard-sphere mixtures. This system is treated within an effective fluid approximation in which the large particles interact through a depletion potential [R. Roth, Phys. Rev. E 62 5360 (2000)] designed to capture the effects of a virtual sea of small particles. We generalize this depletion potential to include the effects of explicit size dispersity in the large particles and consider the case in which the particle diameters are distributed according to a Schulz form having a degree of polydispersity 14%. The resulting alteration (with respect to the monodisperse limit) of the metastable fluid-fluid critical point parameters is determined for two values of the ratio of the diameters of the small and large particles: q(triple bond)sigma(s)/(-)sigma(b)=0.1 and q=0.05. We find that the inclusion of polydispersity moves the critical point to lower reservoir volume fractions of the small particles and high volume fractions of the large ones. The estimated critical point parameters are found to be in good agreement with those predicted by a generalized corresponding states argument which provides a link to the known critical adhesion parameter of the adhesive hard-sphere model. Finite-size scaling estimates of the cluster percolation line in the one phase fluid region indicate that inclusion of polydispersity moves the critical point deeper into the percolating regime. This suggests that phase separation is more likely to be preempted by dynamical arrest in polydisperse systems.

Journal Article↗

PubMed2006

Freezing line of the Lennard-Jones fluid: a phase switch Monte Carlo study.

We report a phase switch Monte Carlo (PSMC) method study of the freezing line of the Lennard-Jones (LJ) fluid. Our work generalizes to soft potentials the original application of the method to hard-sphere freezing and builds on a previous PSMC study of the LJ system by Errington [J. Chem. Phys. 120, 3130 (2004)]. The latter work is extended by tracing a large section of the Lennard-Jones freezing curve, the results of which we compare with a previous Gibbs-Duhem integration study. Additionally, we provide new background material regarding the statistical-mechanical basis of the PSMC method and extensive implementation details.

Journal Article↗

PubMed2005

Finite-size scaling and particle-size cutoff effects in phase-separating polydisperse fluids.

We study the liquid-vapor phase behavior of a polydisperse fluid using grand canonical simulations and moment free energy calculations. The strongly nonlinear variation of the fractional volume of liquid across the coexistence region prevents naive extrapolation from detecting the cloud point. We describe a finite-size scaling method which, nevertheless, permits accurate determination of cloud points from simulations of a single system size. By varying a particle-size cutoff, we find that the cloud point density is highly sensitive to the presence of rare large particles; this could affect the reproducibility of experimentally measured phase behavior in colloids and polymers.

Journal Article↗

PubMed2005

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, rho0(sigma), of the system is controlled by the imposed chemical potential distribution mu(sigma), the form of which is specified such that rho0(sigma) assumes a Schulz form with associated degree of polydispersity approximately = 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 variation 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, Phys. Rev. E 65, 021503 (2002)] a segregation of smaller particles to the interface. The magnitude of this effect as a function of sigma is quantified via measurements of the relative adsorption. Additionally, we consider the utility of various estimators for the interfacial width and highlight the difficulty of isolating the intrinsic contribution of polydispersity to this width.

Journal Article↗

PubMed2004

Polydisperse hard spheres at a hard wall.

The structural properties of polydisperse hard spheres in the presence of a hard wall are investigated via Monte Carlo simulation and density functional theory (DFT). Attention is focused on the local density distribution rho(sigma,z), measuring the number density of particles of diameter sigma at a distance z from the wall. Estimates of rho(sigma,z) are obtained for bulk volume fractions eta(b)=0.2 and eta(b)=0.4, and for two choices of the bulk parent distribution: a top-hat form, which we study for degrees of polydispersity delta=11.5% and delta=40.4%, and a truncated Schulz form having delta=40.7%. Excellent overall agreement is found between the DFT and simulation results, particularly at eta(b)=0.2. A detailed analysis of rho(sigma,z) confirms the presence of oscillatory size segregation effects, as observed in a previous DFT study [I. Pagonabarraga, M. E. Cates, and G. J. Ackland, Phys. Rev. Lett. 84, 911 (2000)]. For large delta, the character of these oscillation is observed to depend strongly on the shape of the parent distribution. In the vicinity of the wall, attractive sigma-dependent depletion interactions are found to greatly enhance the density of the largest particles. The local degree of polydispersity delta(z) is suppressed in this region, while further from the wall it exhibits oscillations.

Journal Article↗

PubMed2004

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 (delta approximately 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 interparticle 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 behavior will occur whenever polydispersity affects only the range, rather than the strength of the interparticle interactions.

Journal Article↗

PubMed2003

Continuous demixing at liquid-vapor coexistence in a symmetrical binary fluid mixture.

We report a Monte Carlo finite-size scaling study of the demixing transition of a symmetrical Lennard-Jones binary fluid mixture. For equal concentration of species, and for a choice of the unlike-to-like interaction ratio delta=0.7, this transition is found to be continuous at liquid-vapor coexistence. The associated critical end point exhibits an Ising-like universality. These findings confirm those of earlier smaller scale simulation studies of the same model, but contradict the findings of recent integral equation and hierarchical reference theory investigations.

Journal Article↗

PubMed2002

Phase behavior and thermodynamic anomalies of core-softened fluids.

We report extensive simulation studies of phase behavior in single component systems of particles interacting via a core-softened interparticle potential. Two recently proposed examples of such potentials are considered; one in which the hard core exhibits a shoulder [Sadr-Lahijany et al., Phys. Rev. Lett. 81, 4895 (1998)], and the other in which the softening takes the form of a linear ramp [Jagla, Phys. Rev. E 63, 061501 (2001)]. Using a combination of state-of-the-art Monte Carlo methods, we obtain the gas, liquid, and solid phase behavior of the shoulder model in two dimensions. We then focus on the thermodynamic anomalies of the liquid phase, namely, maxima in the density and compressibility as a function of temperature. Analysis of the finite-size behavior of these maxima suggests that, rather than stemming from a metastable liquid-liquid critical point, as previously supposed, they are actually induced by the quasicontinuous nature of the two dimensional freezing transition. For the ramp model in three dimensions, we confirm the existence of a stable liquid-liquid ("second") critical point occurring at higher pressure and lower temperature than the liquid-gas critical point. Both these critical points and portions of their associated coexistence curves are located to high precision. In contrast to the shoulder model, the observed thermodynamic anomalies of this model are found to be authentic, i.e., they are not engendered by an incipient new phase. We trace the locus of density and compressibility maxima, the former of which appears to terminate close to the second critical point.

Journal Article↗