Colloidal liquids interacting with short range attraction and long range repulsion, such as proposed for some protein solutions, have been found to exhibit novel states consisting of equilibrium particle clusters. Monte Carlo simulations are performed for two physically meaningful inter-particle potentials across a broad range of interaction parameters, temperatures and volume fractions to locate the conditions where clustered states are found. A corresponding states phase behavior is identified when normalized by the critical point of an appropriately selected reference attractive fluid. Clustered fluid states and cluster percolated states are found exclusively within the two phase region of the state diagram for a reference attractive fluid, confirming the underlying intrinsic relation between clustered states and bulk phase separation. Clustered and cluster percolated states consistently exhibit an intermediate range order peak in their structure factors with a magnitude above 2.7, leading to a semi-empirical rule for identifying clustered fluids in scattering experiments.

The fluid phase behavior of colloidal suspensions with short-range attractive interactions is studied by means of Monte Carlo computer simulations and two theoretical approximations, namely, the discrete perturbation theory and the so-called self-consistent Ornstein-Zernike approximation. The suspensions are modeled as hard-core attractive Yukawa (HCAY) and Asakura-Oosawa (AO) fluids. A detailed comparison of the liquid-vapor phase diagrams obtained through different routes is presented. We confirm Noro-Frenkel's extended law of scaling according to which the properties of a short-ranged fluid at a given temperature and density are independent of the detailed form of the interaction, but just depend on the value of the second virial coefficient. By mapping the HCAY and AO fluids onto an equivalent square-well fluid of appropriate range at the critical point we show that the critical temperature as a function of the effective range is independent of the interaction potential, i.e., all curves fall in a master curve. Our findings are corroborated with recent experimental data for lysozyme proteins.

The so-called extended law of corresponding states, as proposed by Noro and Frenkel [J. Chem. Phys. 113, 2941 (2000)], involves a mapping of the phase behaviors of systems with short-range attractive interactions. While it has already extensively been applied to various model potentials, here we test its applicability to protein solutions with their complex interactions. We successfully map their experimentally determined metastable gas-liquid binodals, as available in the literature, to the binodals of short-range square-well fluids, as determined by previous as well as new Monte Carlo simulations. This is achieved by representing the binodals as a function of the temperature scaled with the critical temperature (or as a function of the reduced second virial coefficient) and the concentration scaled by the cube of an effective particle diameter, where the scalings take into account the attractive and repulsive contributions to the interaction potential, respectively. The scaled binodals of the protein solutions coincide with simulation data of the adhesive hard-sphere fluid. Furthermore, once the repulsive contributions are taken into account by the effective particle diameter, the temperature dependence of the reduced second virial coefficients follows a master curve that corresponds to a linear temperature dependence of the depth of the square-well potential. We moreover demonstrate that, based on this approach and cloud-point measurements only, second virial coefficients can be estimated, which we show to agree with values determined by light scattering or by Derjaguin-Landau-Verwey-Overbeek (DLVO)-based calculations.

One major goal in condensed matter is identifying the physical mechanisms that lead to arrested states of matter, especially gels and glasses. The complex nature and microscopic details of each particular system are relevant. However, from both scientific and technological viewpoints, a general, consistent and unified definition is of paramount importance. Through Monte Carlo computer simulations of states identified in experiments, we demonstrate that dynamical arrest in adhesive hard-sphere dispersions is the result of rigidity percolation with coordination number 〈n(b)〉 equal to 2.4. This corresponds to an established mechanism leading to mechanical transitions in network-forming materials [Phys. Rev. Lett. 54, 2107 (1985)]. Our findings connect the concept of critical gel formation in colloidal suspensions with short-range attractive interactions to the universal concept of rigidity percolation. Furthermore, the bond, angular, and local distributions along the gelation line are explicitly studied in order to determine the topology of the structure at the critical gel state.