Using Newtonian and Brownian dynamics simulations, the structural and transport properties of hard and soft spheres have been studied. The soft spheres were modeled using inverse power potentials (V approximately r(-n), with 1n the potential softness). Although, at constant density, the pressure, diffusion coefficient, and viscosity depend on the particle softness up to extremely high values of n, we show that scaling the density with the freezing point for every system effectively collapses these parameters for n > or = 18 (including hard spheres) for large densities. At the freezing points, the long range structure of all systems is identical, when length is measured in units of the interparticle distance, but differences appear at short distances (due to the different shapes of the interaction potential). This translates into differences at short times in the velocity and stress autocorrelation functions, although they concur to give the same value of the corresponding transport coefficient (for the same density to freezing ratio); the microscopic dynamics also affects the short time behavior of the correlation functions and absolute values of the transport coefficients, but the same scaling with the freezing density works for Newtonian or Brownian dynamics. For hard spheres, the short time behavior of the stress autocorrelation function has been studied in detail, confirming quantitatively the theoretical forms derived for it.

The equilibrium phase diagram of a colloidal system composed of 1:1 mixture of positive and negative particles with equal charge is studied by means of Monte Carlo simulations. The system is the colloidal analog of the restricted primitive model (RPM) for ionic fluids. A liquid-gas transition is found in the low-temperature-low-density region, similar to the liquid-gas transition in the RPM. The fluid-crystal transition is also studied, and the liquid phase is shown to be stable in a narrow range of temperatures. In the liquid, the pair distribution function shows alternating layers of particles with opposite sign of charge surrounding every particle. In the vapor phase, clusters of particles are observed, again in agreement with the RPM. However, a decreasing distribution of clusters is obtained, instead of the discrimination between charged and neutral clusters found in the RPM.

We report on theoretical results concerning the relation between the liquid-liquid transition and the density anomaly for a family of ramp potentials (hard-core plus linear short range repulsion and linear long range attraction). Using first order perturbation, we have studied the influence of the range of the attractive interactions, taking the repulsive part of the interaction as the reference system. Two different mechanisms of liquid-liquid coexistence have been predicted: attraction and compression. The attractive case is attributed to long ranged potentials, while the second one is obtained when the interaction is shortened. The density anomaly appears linked to regions where the temperature derivative of the density derivative of the energy is bigger (in absolute value) than a limit. This condition is fulfilled when the range of the attractive part of the potential is short enough.

Articles you may be interested inThe liquid-gas transition of an electroneutral mixture of oppositely charged colloids, studied by Monte Carlo simulations, is found in the low-temperature-low-density region. The critical temperature shows a nonmonotonous behavior as a function of the interaction range, −1 , with a maximum at Ϸ 10, implying an island of coexistence in the -plane. The system is arranged in such a way that each particle is surrounded by shells of particles with alternating charge. In contrast with the electrolyte primitive model, both neutral and charged clusters are obtained in the vapor phase.

We computed the complete phase diagram of the symmetrical colloidal electrolyte by means of Monte Carlo simulations. Thermodynamic integration, together with the Einstein-crystal method, and Gibbs-Duhem integration were used to calculate the equilibrium phase behavior. The system was modeled via the linear screening theory, where the electrostatic interactions are screened by the presence of salt in the medium, characterized by the inverse Debye length, kappa (in this work kappasigma=6). Our results show that at high temperature, the hard-sphere picture is recovered, i.e., the liquid crystallizes into a fcc crystal that does not exhibit charge ordering. In the low temperature region, the liquid freezes into a CsCl structure because charge correlations enhance the pairing between oppositely charged colloids, making the liquid-gas transition metastable with respect to crystallization. Upon increasing density, the CsCl solid transforms into a CuAu-like crystal and this one, in turn, transforms into a tetragonal ordered crystal near close packing. Finally, we have studied the ordered-disordered transitions finding three triple points where the phases in coexistence are liquid-CsCl-disordered fcc, CsCl-CuAu-disordered fcc, and CuAu-tetragonal-disordered fcc.