Phase diagram of the hard‐sphere/attractive‐Yukawa system

Modeling colloidal interactions that predict equilibrium and non-equilibrium states

Modulating the interaction potential between colloids suspended in a fluid can trigger equilibrium phase transitions as well as formation of non-equilibrium 'arrested states' such as gels and glasses. Faithful representation of such interactions are essential for using simulation to interrogate the microscopic details of non-equilibrium behavior, and for extrapolating observations to new regions of phase space that are difficult to explore in experiment. Although the extended law of corresponding states predicts equilibrium phases for systems with short-ranged interactions, it proves inadequate for equilibrium predictions of systems with longer-ranged interactions, and for predicting non-equilibrium phenomena in systems with either short-ranged or long-ranged interactions. These shortcomings highlight the need for new approaches to represent and disambiguate interaction potentials that replicate both equilibrium and non-equilibrium phase behavior. In this work, we use experiments and simulations to study a system with long-ranged thermoresponsive colloidal interactions and explore whether a resolution to this challenge can be found in regions of the phase diagram where temporal effects influence material state. We demonstrate that the conditions for non-equilibrium arrest by colloidal gelation are sensitive to both the shape of the interaction potential and the thermal quench rate. We exploit this sensitivity to propose a kinetics-based algorithm to extract distinct arrest conditions for candidate potentials that accurately selects between potentials that differ in shape but share the same predicted equilibrium structure. The method reveals that each potential has a quantitatively distinct arrest line, providing insight into how the shape of longer-ranged potentials influences the conditions for colloidal gelation.

Mean-field theory of inhomogeneous fluids

The Barker-Henderson perturbation theory is a bedrock of liquid-state physics, providing quantitative predictions for the bulk thermodynamic properties of realistic model systems. However, this successful method has not been exploited for the study of inhomogeneous systems. We develop and implement a first-principles "Barker-Henderson density functional," thus providing a robust and quantitatively accurate theory for classical fluids in external fields. Numerical results are presented for the hard-core Yukawa model in three dimensions. Our predictions for the density around a fixed test particle and between planar walls are in very good agreement with simulation data. The density profiles for the free liquid vapor interface show the expected oscillatory decay into the bulk liquid as the temperature is reduced toward the triple point, but with an amplitude much smaller than that predicted by the standard mean-field density functional.

Critical examination of the colloidal particle model of globular proteins

Recent studies of globular protein solutions have uniformly adopted a colloidal view of proteins as particles, a perspective that neglects the polymeric primary structure of these biological macromolecules, their intrinsic flexibility, and their ability to sample a large configurational space. While the colloidal perspective often serves as a useful idealization in many cases, the macromolecular identity of proteins must reveal itself under thermodynamic conditions in which the native state is no longer stable, such as denaturing solvents and high protein concentrations where macromolecules tend to have screened excluded volume, charge, and hydrodynamic interactions. Under extreme pH conditions, charge repulsion interactions within the protein chain can overcome the attractive hydrogen-bonding interactions, holding it in its native globular state. Conformational changes can therefore be expected to have great significance on the shear viscosity and other rheological properties of protein solutions. These changes are not envisioned in conventional colloidal protein models and we have initiated an investigation of the scattering and rheological properties of model proteins. We initiate this effort by considering bovine serum albumin because it is a globular protein whose solution properties have also been extensively investigated as a function of pH, temperature, ionic strength, and concentration. As we anticipated, near-ultraviolet circular dichroism measurements and intrinsic viscosity measurements clearly indicate that the bovine serum albumin tertiary structure changes as protein concentration and pH are varied. Our findings point to limited validity of the colloidal protein model and to the need for further consideration and quantification of the effects of conformational changes on protein solution viscosity, protein association, and the phase behavior. Small-angle Neutron Scattering measurements have allowed us to assess how these conformational changes influence protein size, shape, and interprotein interaction strength.

A study of the pair and triplet structures of the quantum hard-sphere Yukawa fluid

The pair and triplet structures of the quantum hard-sphere Yukawa fluid, evaluated for equilateral and isosceles correlations in both the r and the k spaces for a range of conditions and with a particular focus on a region where the onset of increasing number fluctuations takes place (for densities 0.4<or=rho(N) ( *)<or=0.5, along the isotherm lambda(B) ( *)=0.6), are computed via path-integral Monte Carlo simulations in the canonical ensemble and an appropriate Ornstein-Zernike framework. For a given type of correlation (instantaneous, continuous linear response, and centroids), the structural results in r space display how the correlation functions approach each other with decreasing densities as a result of the increasing fluctuations. An attempt at obtaining improved isothermal compressibilities by using a simplified grand-canonical correction to the canonical pair radial functions is also discussed in detail. The results for triplets in k space are based on triplet direct correlation function calculations and are restricted to the higher-density region of the interval studied. Complementary results report an assessment of the performances of the Kirkwood superposition and the Jackson-Feenberg convolution. Comparisons with results also obtained in this work for the bare quantum and the classical hard-sphere fluids are made, allowing one to draw conclusions on the interplay between the inclusion of Yukawa attractions and the quantum diffraction effects in hard-sphere fluids.

Phase diagram and surface tension of the hard-core attractive Yukawa model of variable range: Monte Carlo simulations

The liquid-vapor phase diagram and surface tension for hard-core Yukawa potential with 4<or=kappa<or=7 are calculated by applying canonical Monte Carlo simulation. The authors' new simulation results are more precise than those reported before, which allows them to affirm the accuracy of the previously reported self-consistent Ornstein-Zernike approximation calculations for kappa=5 and 7.

Phase diagrams of alkali halides using two interaction models: A molecular dynamics and free energy study

Phase diagrams for potassium and sodium chlorides are determined by molecular dynamics and free energy calculations. Two rigid-ion interaction models, namely, the Born-Mayer-Huggins (BMH) and Michielsen-Woerlee-Graaf (MWG) effective pair potentials, have been used. The critical and triple point properties are discussed and compared with available experimental and simulation data. The MWG model reproduces the experimental liquid-gas equilibria better than the BMH model, being the accordance very good in the lowest temperature region of the coexistent liquids, particularly for NaCl. However, both models underestimate the critical temperatures of KCl and NaCl. Relatively to the solid-gas equilibria, the models do not reproduce well the experimental data. As for the solid-liquid coexistences either the BMH or the MWG models appear unrealistic.

Phase diagram of the hard-core Yukawa fluid within the integral equation method

In this study, the integral equation method proposed recently by Sarkisov [J. Chem. Phys. 114, 9496 (2001).], which has proved accurate for continuous potentials, is extended successfully to the hard sphere potential plus an attractive Yukawa tail. By comparing the results of thermodynamic properties, including the liquid-vapor phase diagram, with available simulation data, it is found that this method remains reliable for this class of models of interaction often used in colloid science.

Phase and interfacial behavior of partially miscible symmetric Lennard-Jones binary mixtures

We have carried out extensive equilibrium molecular-dynamics simulations to study quantitatively the topology of the temperature versus density phase diagrams and related interfacial phenomena in a partially miscible symmetric Lennard-Jones binary mixture. The topological features are studied as a function of miscibility parameter, alpha = epsilonAB/epsilonAA. Here epsilonAA = epsilonBB and epsilonAB stand for the parameters related to the attractive part of the intermolecular interactions for similar and dissimilar particles, respectively. When the miscibility varies in the range 0 < alpha < 1, a continuous critical line of consolute points Tcons(rho)--critical demixing transition line--appears. This line intersects the liquid-vapor coexistence curve at different positions depending on the values of alpha, yielding mainly three different topologies for the phase diagrams. These results are in qualitative agreement to those found previously for square-well and hard-core Yukawa binary mixtures. The main contributions of the present paper are (i) a quantitative analysis of the phase behavior and (ii) a detailed study of the liquid-liquid interfacial and liquid-vapor surface tensions, as function of temperature and miscibility as well as its relationship to the topological features of the phase diagrams.

Effect of the range of attractive interactions on crystallization, metastable phase transition, and percolation in colloidal dispersions

The equilibrium as well as nonequilibrium phase behaviors of colloidal dispersions have been investigated using statistical-mechanical theories of fluids and solids in complement with the renormalization-group (RG) theory. It is shown that the osmotic second virial coefficient at the critical point of the fluid-fluid transition varies with the range of attractions and is sensitive to specific forms of the attractive potential in contrast to a common speculation that it remains practically constant. However, for colloids with short-ranged forces, the critical temperature of the fluid-fluid phase transition is well correlated with the range of attractions in good agreement with an earlier empirical correlation based on simulation results. A comparison of the relative positions of the fluid-fluid coexistence curve, freezing, melting, and percolation lines in the phase diagram indicates that the gelation in colloidal systems has significant effects on the equilibrium phase transitions and crystallization, especially when the attractions between colloidal particles are short ranged.

Phase equilibria and glass transition in colloidal systems with short-ranged attractive interactions: application to protein crystallization

We have studied a model of a complex fluid consisting of particles interacting through a hard-core and short-range attractive potential of both Yukawa and square-well form. Using a hybrid method, including a self-consistent and quite accurate approximation for the liquid integral equation in the case of the Yukawa fluid, perturbation theory to evaluate the crystal free energies, and mode-coupling theory of the glass transition, we determine both the equilibrium phase diagram of the system and the lines of equilibrium between the supercooled fluid and the glass phases. For these potentials, we study the phase diagrams for different values of the potential range, the ratio of the range of the interaction to the diameter of the repulsive core being the main control parameter. Our arguments are relevant to a variety of systems, from dense colloidal systems with depletion forces, through particle gels, nanoparticle aggregation, and globular protein crystallization.

Thermodynamically self-consistent theories of fluids interacting through short-range forces

The self-consistent Ornstein-Zernike approximation (SCOZA), the generalized mean spherical approximation (GMSA), the modified hypernetted chain (MHNC) approximation, and the hierarchical reference theory (HRT) are applied to the determination of thermodynamic and structural properties, and the phase diagram of the hard-core Yukawa fluid (HCYF). We investigate different Yukawa-tail screening lengths lambda, ranging from lambda=1.8 (a value appropriate to approximate the shape of the Lennard-Jones potential) to lambda=9 (suitable for a simple one-body modelization of complex fluids like colloidal suspensions and globular protein solutions). The comparison of the results obtained with computer simulation data shows that at relatively low lambda's all the theories are fairly accurate in the prediction of thermodynamic and structural properties; as far as the phase diagram is concerned, the SCOZA and HRT are able to predict the binodal line and the critical parameters in a quantitative manner. At lambda=4 some discrepancies begin to emerge in the performances of the different theoretical approaches: the MHNC remains, on the whole, reasonably accurate in predicting the energy and the contact value of the radial distribution function; the SCOZA predicts well the equation of state up to the highest lambda values investigated. The GMSA and the MHNC underestimate and overestimate, respectively, the liquid coexisting density, while the SCOZA and HRT yield liquid branches that fall between the two former theoretical predictions, although both appear to overestimate the critical temperature somewhat. At higher lambda's the GMSA and MHNC binodals further worsen, while the SCOZA appears to remain usefully predictive. In general, the predictions of all the theories tend to slightly worsen at low temperatures and high density. The determination of the freezing line, performed by means of a one-phase "freezing criterion" (due to other authors) is not particularly satisfactory within either the SCOZA or the MHNC; the GMSA prediction for the freezing line at lambda=7 and 9 is instead able to follow in a qualitative manner the pattern of the solid-vapor coexistence line as determined through computer simulation studies. The necessity of further assessments of the freezing predictions is also discussed. Finally, versions of the GMSA, SCOZA, and HRT that can be expected to be more accurate for interactions with extremely short-ranged attractions are identified.