Short-ranged Attractive Colloids Research Articles

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.

Emergent collective behavior of active Brownian particles with visual perception.

Systems comprised of self-steering active Brownian particles are studied via simulations for a minimal cognitive flocking model. The dynamics of the active Brownian particles is extended by an orientational response with limited maneuverability to an instantaneous visual input of the positions of neighbors within a vision cone and a cut-off radius. The system exhibits large-scale self-organized structures, which depend on selected parameter values, and, in particular, the presence of excluded-volume interactions. The emergent structures in two dimensions, such as worms, worm-aggregate coexistence, and hexagonally close-packed structures, are analysed and phase diagrams are constructed. The analysis of the particle's mean-square displacement shows ABP-like dynamics for dilute systems and the worm phase. In the limit of densely packed structures, the active diffusion coefficient is significantly smaller and depends on the number of particles in the cluster. Our analysis of the cluster-growth dynamics shows distinct differences to processes in systems of short-range attractive colloids in equilibrium. Specifically, the characteristic time for the growth and decay of clusters of a particular size is longer than that of isotropically attractive colloids, which we attribute to the non-reciprocal nature of the directed visual perception. Our simulations reveal a strong interplay between ABP-characteristic interactions, such as volume exclusion and rotational diffusion, and cognitive-based interactions and navigation.

The extended law of corresponding states when attractions meet repulsions.

Short-range attractive colloids show well-defined phase behaviour in the absence of repulsions, and highly intriguing equilibrium gelation in the presence of long-range repulsions. We present the state diagram of short-range attractive colloids with repulsions that range from fully screened to intermediately ranged, i.e. longer-ranged than the attractions, but shorter ranged than the colloid size. We demonstrate that although the macroscopic phase behaviour does not change perceptibly, there is a dramatic increase of inhomogeneities once the repulsions become longer-ranged than the attractions. The interaction potentials are characterized with small angle neutron scattering, and used to renormalize the state diagram with the minimum in the interaction potential, min[U(r)], and with the reduced second virial coefficient, B2*. We find that the extended law of corresponding states captures the onset of phase separation for shorter ranged repulsions, but fails for longer ranged repulsions. Instead, for a given model of U(r), the transition from visually homogeneous fluid to phase separation and/or gelation can be rescaled with min[U(r)] over the full range of repulsions. Finally, we suggest a generic state diagram to describe the effect of repulsions on short-range attractive systems.

Analyses of kinetic glass transition in short-range attractive colloids based on time-convolutionless mode-coupling theory.

For short-range attractive colloids, the phase diagram of the kinetic glass transition is studied by time-convolutionless mode-coupling theory (TMCT). Using numerical calculations, TMCT is shown to recover all the remarkable features predicted by the mode-coupling theory for attractive colloids: the glass-liquid-glass reentrant, the glass-glass transition, and the higher-order singularities. It is also demonstrated through the comparisons with the results of molecular dynamics for the binary attractive colloids that TMCT improves the critical values ofthe volume fraction. In addition, a schematic model of three control parameters is investigated analytically. It is thus confirmed that TMCT can describe the glass-glass transition and higher-order singularities even in such a schematic model.

Nucleation in Short-Range Attractive Colloids: Ordering and Symmetry of Clusters

Results from extensive Brownian dynamics simulations are presented for nucleation in a system of colloidal particles interacting via a short-range attractive potential. Our analysis shows that, even though the system is not in equilibrium, structures of small size clusters compare well with the theoretically predicted and experimentally observed ground state structures for short-range colloidal systems. In addition, the distribution of the symmetric structures in nucleation is comparable to the distribution seen in equilibrium. We also investigate how the shape and structure of fluctuating clusters in the prenucleation regime affect the formation of a stable nucleating cluster.

Arrested demixing opens route to bigels

Understanding and, ultimately, controlling the properties of amorphous materials is one of the key goals of material science. Among the different amorphous structures, a very important role is played by colloidal gels. It has been only recently understood that colloidal gels are the result of the interplay between phase separation and arrest. When short-ranged attractive colloids are quenched into the phase-separating region, density fluctuations are arrested and this results in ramified amorphous space-spanning structures that are capable of sustaining mechanical stress. We present a mechanism of aggregation through arrested demixing in binary colloidal mixtures, which leads to the formation of a yet unexplored class of materials--bigels. This material is obtained by tuning interspecies interactions. Using a computer model, we investigate the phase behavior and the structural properties of these bigels. We show the topological similarities and the geometrical differences between these binary, interpenetrating, arrested structures and their well-known monodisperse counterparts, colloidal gels. Our findings are supported by confocal microscopy experiments performed on mixtures of DNA-coated colloids. The mechanism of bigel formation is a generalization of arrested phase separation and is therefore universal.

Structure, compressibility factor, and dynamics of highly size-asymmetric binary hard-disk liquids

By using event-driven molecular dynamics simulation, we investigate effects of varying the area fraction of the smaller component on structure, compressibility factor, and dynamics of the highly size-asymmetric binary hard-disk liquids. We find that the static pair correlations of the large disks are only weakly perturbed by adding small disks. The higher-order static correlations of the large disks, by contrast, can be strongly affected. Accordingly, the static correlation length deduced from the bond-orientation correlation functions first decreases significantly and then tends to reach a plateau as the area fraction of the small disks increases. The compressibility factor of the system first decreases and then increases upon increasing the area fraction of the small disks and separating different contributions to it allows to rationalize this non-monotonic phenomenon. Furthermore, adding small disks can influence dynamics of the system in quantitative and qualitative ways. For the large disks, the structural relaxation time increases monotonically with increasing the area fraction of the small disks at low and moderate area fractions of the large disks. In particular, "reentrant" behavior appears at sufficiently high area fractions of the large disks, strongly resembling the reentrant glass transition in short-ranged attractive colloids and the inverted glass transition in binary hard spheres with large size disparity. By tuning the area fraction of the small disks, relaxation process for the small disks shows concave-to-convex crossover and logarithmic decay behavior, as found in other binary mixtures with large size disparity. Moreover, diffusion of both species is suppressed by adding small disks. Long-time diffusion for the small disks shows power-law-like behavior at sufficiently high area fractions of the small disks, which implies precursors of a glass transition for the large disks and a localization transition for the small disks. Therefore, our results demonstrate the generic dynamic features in highly size-asymmetric binary mixtures.

Pathways of Cluster Growth and Kinetic Slowing Down in a Model of Short-Range Attractive Colloids

We present the results of an extensive 3D Brownian dynamics simulation of the self-assembly of colloidal particles for a short-range attractive model that is quenched below its metastable critical point. In particular, results are obtained in the small-volume-fraction, low-temperature region in which we find so-called sticky beads that diffuse around the system, without reaching a final large cluster on the timescale of our simulation. For larger volume fractions in this low-temperature regime, a gel forms as the result of kinetically slowed down spinodal decomposition, as shown earlier for other short-range attractive models (Foffi, G.; De Michele, C.; Sciortino, F.; Tartaglia, P. Phys. Rev. Lett. 2005, 94, 078301. Zaccarelli, E. J. Phys.: Condens. Matter 2007, 19, 323101). We also show that for quenches below the critical point but above the intersection of the binodal with the glass line, two-step crystallization takes place. For sufficiently small volume fractions, the first step is the nucleation of dense fluid drops, followed by the second step of crystallization within these drops, as first proposed for a model of protein crystallization for quenches just above the metastable critical point (ten Wolde, P. R.; Frenkel, D. Science 1997, 277, 1975). For larger values of the volume fraction, the initial step is spinodal decomposition that leads to the formation of an interconnected network of low- and high-density fluids. The second step is crystallization that takes place within the dense fluid phase.

Kinetics and morphology of cluster growth in a model of short-range attractive colloids

We present results from detailed three-dimensional Brownian dynamics simulations of the self-assembly process in quenched short-range attractive colloids. Clusters obtained in the simulations range from dense faceted crystals to fractal aggregates which show ramified morphology on large length scales but close-packed crystalline morphology on short length scales. For low volume fractions of the colloids, the morphology and crystal structure of a nucleating cluster are studied at various times after the quench. As the volume fraction of the colloids is increased, growth of clusters is controlled by cluster diffusion and cluster-cluster interactions. For shallower quenches and low volume fractions, clusters are compact and the growth-law exponent agrees well with Binder-Stauffer predictions and with recent experimental results. As the volume fraction is increased, clusters do not completely coalesce when they meet each other and the kinetics crosses over to diffusion-limited cluster-cluster aggregation (DLCA) limit. For deeper quenches, clusters are fractals even at low volume fractions and the growth kinetics asymptotically reaches the irreversible DLCA case.

Elasticity of Arrested Short-Ranged Attractive Colloids: Homogeneous and Heterogeneous Glasses

We evaluate the elasticity of arrested short-ranged attractive colloids by combining an analytically solvable elastic model with a hierarchical arrest scheme. This new approach allows us to discriminate the microscopic (primary particle-level) from the mesoscopic (cluster-level) contribution to the macroscopic shear modulus. The results quantitatively predict experimental data in a wide range of volume fractions and indicate in which cases the relevant contribution is due to mesoscopic structures. On this basis we propose that different arrested states of short-ranged attractive colloids can be meaningfully distinguished as homogeneous or heterogeneous colloidal glasses in terms of the length scale which controls their elastic behavior.

Effective interaction between large colloidal particles immersed in a bidisperse suspension of short-ranged attractive colloids

The effective force between two large hard spheres mimicking lyophobic colloids (solute) immersed in an asymmetric two-component mixture of smaller particles (solvents), interacting via Baxter's sticky hard sphere (SHS) potential, was studied using integral equation theory and Monte Carlo simulation. The theoretical predictions were calculated from the analytic solution of the Percus-Yevick/Ornstein-Zernike integral equation for spatial correlations in a three-component mixture at vanishing solute concentration, while the simulation results were obtained by applying a special simulation technique developed for sampling the hard-sphere collision force. Due to layering of the solvent molecules, the effective force between the particles of the solute oscillates with periods equal to the molecular diameters of both solvent components. The attractive force between the solute particles in the SHS mixture comprising strongly attractive molecules of either component decays slower than that in the mixture with weaker interparticle attraction. Similar features are also observed when inspecting the separate contributions of individual components to the total solute-solute force. At sufficient strength of the interparticle stickiness, these oscillations disappear, the force becoming long ranged and attractive at all separations.

Crystallization and gelation in colloidal systems with short-ranged attractive interactions

We systematically study the relationship between equilibrium and nonequilibrium phase diagrams of a system of short-ranged attractive colloids. Using Monte Carlo and Brownian dynamics simulations we find a window of enhanced crystallization that is limited at high interaction strength by a slowing down of the dynamics and at low interaction strength by the high nucleation barrier. We find that the crystallization is enhanced by the metastable gas-liquid binodal by means of a two-stage crystallization process. First, the formation of a dense liquid is observed and second the crystal nucleates within the dense fluid. In addition, we find at low colloid packing fractions a fluid of clusters, and at higher colloid packing fractions a percolating network due to an arrested gas-liquid phase separation that we identify with gelation. We find that this arrest is due to crystallization at low interaction energy and it is caused by a slowing down of the dynamics at high interaction strength. Likewise, we observe that the clusters which are formed at low colloid packing fractions are crystalline at low interaction energy, but glassy at high interaction energy. The clusters coalesce upon encounter.

Short range attractive colloids: dynamics and energy landscape properties

We discuss the static and dynamical properties of a model of short range attractivecolloidal systems. Our model is a binary mixture of particles interacting with a modifiedLennard-Jones potential that we simulated by molecular dynamics. As previously found, weconfirm the existence of a maximum in diffusivity on varying the temperature at constanthigh density. The possibility to characterize such a maximum by studying the properties ofthe potential energy landscape (PEL) is also discussed. This has been achieved bycalculating the configurational entropy following two independent methods. It is shownthat this quantity presents a maximum in correspondence with that of diffusivity.

Kinetic arrest and glass-glass transition in short-ranged attractive colloids

A thermally reversible repulsive hard-sphere to sticky-sphere transition was studied in a model colloidal system over a wide volume fraction range. The static microstructure was obtained from high resolution small angle x-ray scattering, the colloid dynamics was probed by dynamic x-ray and light scattering, and supplementary mechanical properties were derived from bulk rheology. At low concentration, the system shows features of gas-liquid type phase separation. The bulk phase separation is presumably interrupted by a gelation transition at the intermediate volume fraction range. At high volume fractions, fluid-attractive glass and repulsive glass-attractive glass transitions are observed. It is shown that the volume fraction of the particles can be reliably deduced from the absolute scattered intensity. The static structure factor is modeled in terms of an attractive square-well potential, using the leading order series expansion of Percus-Yevick approximation. The ensemble-averaged intermediate scattering function shows different levels of frozen components in the attractive and repulsive glassy states. The observed static and dynamic behavior are consistent with the predictions of a mode-coupling theory and numerical simulations for a square-well attractive system.

Anomalous dynamic arrest in a mixture of large and small particles

We present molecular dynamics simulations of the slow dynamics of a mixture of large and small soft spheres with a large size disparity. The dynamics are investigated in a broad range of temperature and mixture composition. As a consequence of the large size disparity, large and small particles exhibit very different relaxation times. As previously reported for simple models of short-ranged attractive colloids and polymer blends, several anomalous dynamic features are observed: (i) sublinear behavior for mean-squared displacements, (ii) concave-to-convex crossover for density-density correlators, by varying the temperature or wave vector, and (iii) logarithmic decay for specific wave vectors of density-density correlators. These anomalous features are observed over time intervals extending up to four decades and strongly resemble predictions of the mode coupling theory (MCT) for state points close to higher-order MCT transitions, which originate from the competition between different mechanisms for dynamic arrest. For the large particles we suggest competition between soft-sphere repulsion and depletion effects induced by neighboring small particles. For the small particles we suggest competition between bulklike dynamics and confinement, respectively induced by neighboring small particles and by the slow matrix of large particles.

Glassy colloidal systems

This review focuses on recent developments in the theoretical, numerical and experimental study of slow dynamics in colloidal systems, with a particular emphasis on the glass transition phenomenon. Colloidal systems appear to be particularly suited for tackling the general problem of dynamic arrest, since they show a larger flexibility compared to atomic and molecular glasses because of their size and the possibility of manipulating the physical and chemical properties of the samples. Indeed, a wealth of new effects, not easily observable in molecular liquids, have been predicted and measured in colloidal systems. The slow dynamic behavior of three classes of colloidal suspension is reviewed – hard colloids, short-range attractive colloids and soft colloidal systems – selecting the model systems among the most prominent candidates for grasping the essential features of dynamic arrest. Emphasis is on the possibility of performing a detailed comparison between experimental data and theoretical predictions based on the mode coupling theory of the glass transition. Finally, the importance of understanding the system's kinetic arrest phase diagram, i.e. the regions in phase space where disordered arrested states can be expected, is stressed. When and how these states are kinetically stabilized with respect to the ordered lowest free energy phases is then examined in order to provide a framework for interpreting and developing new ideas in the study of new materials. Glassy colloidal systemsAll authorsF. Sciortino & P. Tartagliahttps://doi.org/10.1080/00018730500414570Published online:19 February 2007Table Download CSVDisplay Table

Scaling of Dynamics with the Range of Interaction in Short-Range Attractive Colloids

We numerically study the dependence of the dynamics on the range of interaction Delta for the short-range square well potential. We find that, for small Delta, dynamics scale exactly in the same way as thermodynamics, both for Newtonian and Brownian microscopic dynamics. For interaction ranges from a few percent down to the Baxter limit, the relative location of the attractive-glass line and the liquid-gas line does not depend on Delta. This proves that, in this class of potentials, disordered arrested states (gels) can be generated only as a result of a kinetically arrested phase separation.

Numerical study of the glass–glass transition in short-ranged attractive colloids

We report extensive numerical simulations in the glass region for a simple model ofshort-ranged attractive colloids, the square well model. We investigate the behaviour of thedensity autocorrelation function and of the static structure factor in the region oftemperatures and packing fractions where a glass–glass transition is expected according totheoretical predictions. We strengthen our observations by studying both waiting time andhistory dependence of the numerical results. We provide evidence supportingthe possibility that activated bond-breaking processes destabilize the attractiveglass, preventing the full observation of a sharp glass–glass kinetic transition.

Dynamical arrest in dense short-ranged attractive colloids

We study thermodynamic and dynamic properties of model colloidal systems interactingwith a hard core repulsion and a short-range attraction, and provide an overall picture oftheir phase diagrams which shows a very rich phenomenology. We focus on the slowdynamic properties of this model, investigating in detail the glass transition lines (bothrepulsive and attractive), the glass–glass transitions and the location of the higher ordersingularities. We discuss the relative location of the glass lines and of the metastableliquid–gas binodal, an issue relevant for the understanding of low density arrested states ofmatter.

Aging in short-ranged attractive colloids: a numerical study.

We study the aging dynamics in a model for dense simple liquids, in which particles interact through a hard-core repulsion complemented by a short-ranged attractive potential, of the kind found in colloidal suspensions. In this system, at large packing fractions, kinetically arrested disordered states can be created both on cooling (attractive glass) and on heating (repulsive glass). The possibility of having two distinct glasses, at the same packing fraction, with two different dynamics offers the unique possibility of comparing-within the same model-the differences in aging dynamics. We find that, while the aging dynamics of the repulsive glass is similar to the one observed in atomic and molecular systems, the aging dynamics of the attractive glass shows novel unexpected features.