Dino Costa

The structural properties of a two-Yukawa fluid: Simulation and analytical results.

Standard Monte Carlo simulations are carried out to assess the accuracy of theoretical predictions for the structural properties of a model fluid interacting through a hard-core two-Yukawa potential composed of a short-range attractive well next to a hard repulsive core, followed by a smooth, long-range repulsive tail. Theoretical calculations are performed in the framework provided by the Ornstein-Zernike equation, solved either analytically with the mean spherical approximation (MSA) or iteratively with the hypernetted-chain (HNC) closure. Our analysis shows that both theories are generally accurate in a thermodynamic region corresponding to a dense vapor phase around the critical point. For a suitable choice of potential parameters, namely, when the attractive well is deep and/or large enough, the static structure factor displays a secondary low-Q peak. In this case HNC predictions closely follow the simulation results, whereas MSA results progressively worsen the more pronounced this low-Q peak is. We discuss the appearance of such a peak, also experimentally observed in colloidal suspensions and protein solutions, in terms of the formation of equilibrium clusters in the homogeneous fluid.

Theory and simulation of short-range models of globular protein solutions

We report theoretical and simulation studies of phase coexistence in model globular protein solutions, based on short-range, central, pair potential representations of the interaction among macro-particles. After reviewing our previous investigations of hard-core Yukawa and generalized Lennard-Jones potentials, we report more recent results obtained within a DLVO-like description of lysozyme solutions in water and added salt. We show that a one-parameter fit of this model, based on static light scattering and self-interaction chromatography data in the dilute protein regime, yields demixing and crystallization curves in good agreement with experimental protein-rich?protein-poor and solubility envelopes. The dependence of cloud and solubility point temperatures of the model on the ionic strength is also investigated. Our findings highlight the minimal assumptions on the properties of the microscopic interaction sufficient for a satisfactory reproduction of the phase diagram topology of globular protein solutions.

Phase coexistence in a DLVO model of globular protein solutions

Globular protein solutions of lysozyme in water and added salt are modelled according to the Derjaguin–Landau–Verwey–Overbeek (DLVO) theory, in order to determine their fluid–fluid and fluid–solid coexistence lines. Calculations are based on both computer simulations and theoretical approaches. Our results indicate that, when the potential parameters are obtained by fitting physical properties directly deducible from either static or dynamic light scattering data, the fluid–fluid phase coexistence predictions agree quite well with the experiments. Our description of the solid phase allows only a qualitative reproduction of the experimental solubility boundaries. The overall accuracy of our predictions is discussed in view of the well known limitations of the DLVO representation of protein solutions.

Structure and equation of state of interaction site models for disc-shaped lamellar colloids

We apply RISM (reference interaction site model) and PRISM (polymer-RISM) theories to calculate the site–site pair structure and the osmotic equation of state of suspensions of circular or hexagonal platelets (lamellar colloids) over a range of ratios of the particle diameter over thickness Dσ. Despite the neglect of edge effects, the simpler PRISM theory yields results in good agreement with the more elaborate RISM calculations, provided the correct form factor, characterizing the intramolecular structure of the platelets, is used. The RISM equation of state is sensitive to the number n of sites used to model the platelets, but saturates when the hard spheres, associated with the interaction sites, nearly touch; the limiting equation of state agrees reasonably well with available simulation data for all densities up to the isotropic–nematic transition. When properly scaled with the second virial coefficient, the equations of state of platelets with different aspect ratios Dσ nearly collapse on a single master curve.

Temperature study of cluster formation in two-Yukawa fluids

An accurate thermodynamically self-consistent integral equation theory of the liquid state is used to investigate model fluids with competing attractive interaction at short distances and long-range repulsion, focusing on thermodynamic states where the formation of clusters is expected to occur. We find a remarkable accuracy of theoretical predictions, through a detailed assessment against results of Monte Carlo simulations. The behavior of theoretical radial distribution functions and structure factors faithfully follows the onset and growth of cluster aggregates in the homogeneous dense-vapor phase. The thermodynamic properties of the system sensitively depend on the ratio between the repulsive barrier and the attraction strength. We elucidate the role of accurate theoretical tools to investigate the properties of fluids with complex phase behaviors.

Communication: Thermodynamic signatures of cluster formation in fluids with competing interactions

Convergent theoretical evidence, based on self-consistent integral equations for the pair structure and on Monte Carlo simulations, is presented for the existence of small simultaneous jump discontinuities of several thermodynamic and structural properties of systems of colloidal particles with competing short-range attractive and long-range repulsive interactions, under physical conditions close to the onset of particle clustering. The discontinuities thus provide a signature of the transition from a homogeneous fluid phase to a locally inhomogeneous cluster phase.

A comprehensive study of the phase diagram of symmetrical hard-core Yukawa mixtures

The phase diagrams of hard-core Yukawa mixtures (HCYM), constituted of equal sized hard spheres interacting through an attractive Yukawa tail, are determined by means of Gibbs Ensemble Monte Carlo (GEMC) simulations, Semi-grand Canonical Monte Carlo (SGCMC) simulations, and through the modified hypernetted-chain (MHNC) theory. Freezing lines are obtained according to an approach recently proposed by Giaquinta and co-workers [Physica A 187, 145 (1992); Phys Rev. A 45, 6966 (1992)] in which an analysis of multiparticle contributions to the excess entropy, Δs, is performed, with the determination of the Δs=0 locus. Liquid–vapor coexistence, determined through GEMC simulations, turns out to be favored when the strength ratio ν of unlike to like particle interaction, is close to 1. For lower ν’s, liquid–vapor coexistence is favored at low densities, and liquid–liquid coexistence, determined through SGCMC simulations, at high densities. The liquid–vapor binodal shifts downward in temperature and flattens when ν...

Simulation and theory of a model for tetrahedral colloidal particles.

We study the thermodynamic and structural properties of a five-site tetrahedral molecular model by means of different Monte Carlo simulation techniques, and the reference interaction site model (RISM) theory of molecular fluids. Simulations and theory signal the onset, at sufficiently low temperatures, of two different tetrahedral molecular arrangements, with a more open topology progressively giving place to a fully bonded one, as the temperature decreases. The RISM theory reproduces the splitting of the static structure factor at low temperatures, a feature intimately related to the onset of the tetrahedral ordering. Less accurate predictions are obtained for the liquid-vapor coexistence and the short-range correlations.

Free energy determination of phase coexistence in model C60: A comprehensive Monte Carlo study

The free energy of the solid and fluid phases of the Girifalco C60 model are determined through extensive Monte Carlo simulations. In this model the molecules interact through a spherical pair potential, characterized by a narrow and attractive well, adjacent to a harshly repulsive core. We have used the Widom test particle method and a mapping from an Einstein crystal in order to estimate the absolute free energy in the fluid and solid phases, respectively; we have then determined the free energy along several isotherms, and the whole phase diagram, by means of standard thermodynamic integrations. The dependence of the simulation’s results on the size of the sample is also monitored in a number of cases. We highlight how the interplay between the liquid–vapor and the liquid–solid coexistence conditions determines the existence of a narrow liquid pocket in the phase diagram, whose stability is assessed and confirmed in agreement with previous studies. In particular, the critical temperature follows closel...

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid–fluid equilibrium, with a metastable fluid–fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76×106 MD steps. By progressively reducing the system temperature below the solid–fluid line, we first observe the metastable fluid–fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the pre...