Patchy colloidal gels under the influence of gravity.

In this contribution, gravitational effects in gel-forming patchy colloidal systems are studied. We focus on how the gel structure is modified by gravity. Through Monte Carlo computer simulations of gel-like states recently identified by the rigidity percolation criterion [J. A. S. Gallegos et al., Phys. Rev. E 104, 064606 (2021)], the influence of the gravitational field, characterized by the gravitational Péclet number, Pe, on patchy colloids is studied in terms of the patchy coverage, χ. Our findings point out that there exists a threshold Péclet number, Peg, that depends on χ above which the gravitational field enhances the particle bonding and, in consequence, promotes the aggregation or clustering of particles; the smaller the χ value, the higher the Peg. Interestingly, when χ ∼ 1 (near the isotropic limit), our results are consistent with an experimentally determined threshold Pe value where gravity affects the gel formation in short-range attractive colloids. In addition, our results show that the cluster size distribution and the density profile undergo variations that lead to changes in the percolating cluster, i.e., gravity is able to modify the structure of the gel-like states. These changes have an important impact on the structural rigidity of the patchy colloidal dispersion; the percolating cluster goes from a uniform spatially network to a heterogeneous percolated structure, where an interesting structural scenario emerges, namely, depending on the Pe value, the new heterogeneous gel-like states can coexist with both diluted and dense phases or they simply reach a crystalline-like state. In the isotropic case, the increase in the Pe number can shift the critical temperature to higher temperatures; however, when Pe > 0.01, the binodal disappears and the particles fully sediment at the bottom of the sample cell. Furthermore, gravity moves the rigidity percolation threshold to lower densities. Finally, we also note that within the values of the Péclet number here explored, the cluster morphology is barely altered.

Potential-invariant network structures in Asakura-Oosawa mixtures with very short attraction range.

We systematically investigated the structure and aggregate morphology of gel networks formed by colloid-polymer mixtures with a moderate colloid volume fraction and different values of the polymer-colloid size ratio, always in the limit of short-range attraction. Using the coordinates obtained from confocal microscopy experiments, we determined the radial, angular, and nearest-neighbor distribution functions together with the cluster radius of gyration as a function of size ratio and polymer concentration. The analysis of the structural correlations reveals that the network structure becomes increasingly less sensitive to the potential strength with the decreasing polymer-colloid size ratio. For the larger size ratios, compact clusters are formed at the onset of network formation and become progressively more branched and elongated with increasing polymer concentration/attraction strength. For the smallest size ratios, we observe that the aggregate structures forming the gel network are characterized by similar morphological parameters for different values of the size ratio and the polymer concentration, indicating a limited evolution of the gel structure with variations of the parameters that determine the interaction potential between colloids.

Reversible Aggregation and Colloidal Cluster Morphology: The Importance of the Extended Law of Corresponding States.

Cluster morphology of spherical particles interacting with a short-range attraction has been extensively studied due to its relevance to many applications, such as the large-scale structure in amorphous materials, phase separation, protein aggregation, and organelle formation in cells. Although it was widely accepted that the range of the attraction solely controls the fractal dimension of clusters, recent experimental results challenged this concept by also showing the importance of the strength of attraction. Using Monte Carlo simulations, we conclusively demonstrate that it is possible to reduce the dependence of the cluster morphology to a single variable, namely, the reduced second virial coefficient, B_{2}^{*}, linking the local properties of colloidal systems to the extended law of corresponding states. Furthermore, the cluster size distribution exhibits two well-defined regimes: one identified for small clusters, whose fractal dimension, d_{f}, does not depend on the details of the attraction, i.e., small clusters have the same d_{f}, and another related to large clusters, whose morphology depends exclusively on B_{2}^{*}, i.e., d_{f} of large aggregates follows a master curve, which is only a function of B_{2}^{*}. This physical scenario is confirmed with the reanalysis of experimental results on colloidal-polymer mixtures.

Structure of colloidal gels at intermediate concentrations: the role of competing interactions.

Colloidal gels formed by colloid-polymer mixtures with an intermediate volume fraction (ϕc ≈ 0.4) are investigated by confocal microscopy. In addition, we have performed Monte Carlo simulations based on a simple effective pair potential that includes a short-range attractive contribution representing depletion interactions, and a longer-range repulsive contribution describing the electrostatic interactions due to the presence of residual charges. Despite neglecting non-equilibrium effects, experiments and simulations yield similar gel structures, characterised by, e.g., the pair, angular and bond distribution functions. We find that the structure hardly depends on the strength of the attraction if the electrostatic contribution is fixed, but changes significantly if the electrostatic screening is changed. This delicate balance between attractions and repulsions, which we quantify by the second virial coefficient, also determines the location of the gelation boundary.

Extended law of corresponding states for protein solutions.

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.

Generalized phase behavior of cluster formation in colloidal dispersions with competing interactions.

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.

Phase behavior of the modified-Yukawa fluid and its sticky limit.

Simple model systems with short-range attractive potentials have turned out to play a crucial role in determining theoretically the phase behavior of proteins or colloids. However, as pointed out by D. Gazzillo [J. Chem. Phys. 134, 124504 (2011)], one of these widely used model potentials, namely, the attractive hard-core Yukawa potential, shows an unphysical behavior when one approaches its sticky limit, since the second virial coefficient is diverging. However, it is exactly this second virial coefficient that is typically used to depict the experimental phase diagram for a large variety of complex fluids and that, in addition, plays an important role in the Noro-Frenkel scaling law [J. Chem. Phys. 113, 2941 (2000)], which is thus not applicable to the Yukawa fluid. To overcome this deficiency of the attractive Yukawa potential, D. Gazzillo has proposed the so-called modified hard-core attractive Yukawa fluid, which allows one to correctly obtain the second and third virial coefficients of adhesive hard-spheres starting from a system with an attractive logarithmic Yukawa-like interaction. In this work we present liquid-vapor coexistence curves for this system and investigate its behavior close to the sticky limit. Results have been obtained with the self-consistent Ornstein-Zernike approximation (SCOZA) for values of the reduced inverse screening length parameter up to 18. The accuracy of SCOZA has been assessed by comparison with Monte Carlo simulations.

Phase behavior of colloids and proteins in aqueous suspensions: theory and computer simulations.

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.

Clusters in colloidal dispersions with a short-range depletion attraction: Thermodynamic identification and morphology.

HYPOTHESIS: Particle aggregation is ubiquitous for many colloidal systems, and drives the phase separation or the formation of materials with a highly heterogeneous large-scale structure, such as gels, porous media and attractive glasses. While the macroscopic properties of such materials strongly depend on the shape and size of these particle aggregates, the morphology and underlining aggregation physical mechanisms are far from being fully understood. Recently, it has been proposed that for reversible colloidal aggregation, the cluster morphology in the case of colloids interacting with short-range attractive forces is determined by a single variable, namely, the reduced second virial coefficient, B2∗. EXPERIMENTS: We examined this proposal by performing confocal microscopy experiments and computer simulations on a large collection of short-ranged attractive colloidal systems with different values of the attraction strength and range. FINDINGS: We show that in all cases a connection between the colloidal cluster morphology and B2∗ can be established both in experiments and simulations. This physical scenario holds at all investigated thermodynamic conditions, namely, in the fluid state, in the metastable region and in non-equilibrium conditions. Our findings support the connection between reversible colloidal aggregation and the so-called extended law of corresponding states.

Thermodynamics, static properties and transport behaviour of fluids with competing interactions.

Competing interaction fluids have become ideal model systems to study a large number of phenomena, for example, the formation of intermediate range order structures, condensed phases not seen in fluids driven by purely attractive or repulsive forces, the onset of particle aggregation under in- and out-of-equilibrium conditions, which results in the birth of reversible and irreversible aggregates or clusters whose topology and morphology depend additionally on the thermodynamic constrictions, and a particle dynamics that has a strong influence on the transport behaviour and rheological properties of the fluid. In this contribution, we study a system of particles interacting through a potential composed by a continuous succession of a short-ranged square-well (SW), an intermediate-ranged square-shoulder and a long-ranged SW. This potential model is chosen to systematically analyse the contribution of every component of the interaction potential on the phase behaviour, the microstructure, the morphology of the resulting aggregates and the transport phenomena of fluids described by competing interactions. Our results indicate that the inclusion of a barrier and a second well leads to new and interesting effects, which in addition result in variations of the physical properties associated to the competition among interactions.

Efficient Determination of Slip-Link Parameters from Broadly Polydisperse Linear Melts.

We investigate the ability of a coarse-grained slip-link model and a simple double reptation model to describe the linear rheology of polydisperse linear polymer melts. Our slip-link model is a well-defined mathematical object that can describe the equilibrium dynamics and non-linear rheology of flexible polymer melts with arbitrary polydispersity and architecture with a minimum of inputs: the molecular weight of a Kuhn step, the entanglement activity, and Kuhn step friction. However, this detailed model is computationally expensive, so we also examine predictions of the cheaper double reptation model, which is restricted to only linear rheology near the terminal zone. We report the storage and loss moduli for polydisperse polymer melts and compare with experimental measurements from small amplitude oscillatory shear. We examine three chemistries: polybutadiene, polypropylene and polyethylene. We also use a simple double reptation model to estimate parameters for the slip-link model and examine under which circumstances this simplified model works. The computational implementation of the slip-link model is accelerated with the help of graphics processing units, which allow us to simulate in parallel large ensembles made of up to 50,000 chains. We show that our simulation can predict the dynamic moduli for highly entangled polymer melts over nine decades of frequency. Although the double reptation model performs well only near the terminal zone, it does provide a convenient and inexpensive way to estimate the entanglement parameter for the slip-link model from polydisperse data.

Dynamical arrest in adhesive hard-sphere dispersions driven by rigidity percolation.

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.