Hydrodynamics strongly affect the dynamics of colloidal gelation but not gel structure

The impact of colloid-solvent dynamic coupling on the coarsening rate of colloidal phase separation

Phase separation, a fundamental phenomenon in both natural and industrial settings, involves the coarsening of domains over time t to reduce interfacial energy. While well-understood for simple viscous liquid mixtures, the physical laws governing coarsening dynamics in complex fluids, such as colloidal suspensions, remain unclear. Here, we investigate colloidal phase separation through particle-based simulations with and without hydrodynamic interactions (HIs). The former incorporates many-body HIs through momentum conservation, while the latter simplifies their effects into a constant friction coefficient on a particle. In cluster-forming phase separation with HIs, the domain size ℓ grows as ℓ∝t1/3, aligning with the Brownian-coagulation mechanism. Without HIs, ℓ∝t1/5, attributed to an improper calculation of cluster thermal diffusion. For network-forming phase separation, ℓ∝t1/2 with HIs, while ℓ∝t1/3 without HIs. In both cases, network coarsening is governed by the mechanical stress relaxation of the colloid-rich phase, yet with distinct mechanisms: slow solvent permeation through densely packed colloids for the former and free draining for the latter. Our results provide a clear and concise physical picture of colloid-solvent dynamic coupling via momentum conservation, offering valuable insights into the self-organization dynamics of particles like colloids, emulsions, and globular proteins suspended in a fluid.

Tuning higher order structure in colloidal fluids

Colloidal particles self assemble into a wide range of structures under external AC electric fields due to induced dipolar interactions [Yethiraj and Van Blaaderen, Nature, 2003, 421, 513]. As a result of these dipolar interactions, at low volume fraction the system is modulated between a hard-sphere like state (in the case of zero applied field) and a "string fluid" upon application of the field. Using both particle-resolved experiments and computer simulations, we investigate the emergence of the string fluid with a variety of structural measures including two-body and higher-order correlations. We probe the higher-order structure using three-body spatial correlation functions and a many-body approach based on minimum energy clusters of a dipolar-Lennard-Jones system. The latter constitutes a series of geometrically distinct minimum energy clusters upon increasing the strength of the dipolar interaction, which are echoed in the higher-order structure of the colloidal fluids we study here. We find good agreement between experiment and simulation at the two-body level. Higher-order correlations exhibit reasonable agreement between experiment and simulation, again with more discrepancy at higher field strength for three-body correlation functions. At higher field strength, the cluster population in our experiments and simulations is dominated by the minimum energy clusters for all sizes 8 ≤ m ≤ 12.

Kinetic Pathways to Gelation and Effects of Flow-Induced Structuring in Depletion Gels

The kinetic pathways to gelation and the effects of flow-induced restructuring are studied here in depletion flocculated gels with short-ranged attractions, both experimentally and using computer simulations. In the experiments, we first carefully diffuse a screening organic salt to destabilize colloid-polymer mixtures and form a gel. We hence avoid flow history effects, typical of traditional mixing protocols. The initial gelation phases are then accessible and observed by time-resolved confocal microscopy. These insights show that quiescent gelation reduces heterogeneity and strand size with increasing attraction strength, with deeper quenches leading to earlier arrest. These findings are consistent with the simulations which include long-range hydrodynamic interactions. We then compare these results with gels formed by high-rate preshear followed by cessation of colloid-polymer-salt mixtures. The obtained microstructures do not seem in this case to depend on depletant concentration. Indeed, confocal images reveal that shear flow significantly impacts gel structure, from fluidization at high shear rates to dense heterogeneous aggregates formation at lower rates. We especially show how the heterogeneity is controlled by the strength of the flow relative to the attraction forces between the colloids. This study highlights the subtleties behind the preparation protocols of colloidal gels. In particular, it shows that differences in kinetic aggregation pathways can overshadow attraction effects, such as those caused by varying flow conditions during mixing at different attraction strengths. These insights provide a framework for understanding gelation kinetics and optimizing structural reproducibility in colloidal gel experiments.

Void Connectivity and Criticality in the Compression-Induced Gel-to-Glass Transition of Short-Range Attractive Colloids

Gels and attractive glasses-both dynamically arrested states formed through short-range attraction-are commonly found in a range of soft matter systems, including colloids, emulsions, proteins, and wet granular materials. Previous studies have revealed intriguing similarities and distinctions in their structural, dynamic, vibrational, and mechanical properties. However, the microstructural mechanisms underlying the gel-to-glass transition remain elusive. To address this, we investigate uniaxial compression-induced gel collapse using confocal microscopy, which provides experimental access to the relationship between mechanical stress and microstructure. Together with Brownian dynamics simulations, our study reveals two sequential transitions in void structure: from gels with percolating voids to isolated voids, and ultimately to voidless glasses. These transitions are closely linked to a shift from superlinear power-law scaling to explosive divergence toward the packing limit in both normal and deviatoric stresses as a function of volume fraction, particularly at lower temperatures. Understanding these mechanically self-organized transitions and their associated criticality deepens our insight into disordered solids, enabling better control over mechanical properties, interfacial characteristics, and transport behavior in porous materials.

Pairing-specific microstructure in depletion gels of bidisperse colloids

We report the ensemble-averaged and pairing-specific network microstructure formed by short-range depletion attractions in hard sphere-like colloidal systems. Gelation is induced by adding polystyrene molecules at a fixed concentration to colloids with different colloid bidispersity ratios (α = 1, 0.72, and 0.60) across a range of volume fractions (0.10 ≤ ϕ ≤ 0.40). 3D confocal microscopy imaging combined with a scale-invariant feature transform algorithm show that monodisperse colloids pack more efficiently, whereas increasing the size disparity leads to looser, more disordered, and sub-isostatic packings. Categorizing the structures formed by small and large particles reveal that certain cluster configurations may be favored due to the complex interplay between the differences in particle surface areas and attractive potentials. These pairwise bonds assemble to affect the density of tetrahedral and poly-tetrahedral clusters in bidisperse systems. With the exception of non-percolating samples at ϕ = 0.10, increasing the gel volume fraction leads to an increase in the number of nearest neighbors. However, the internal density within each cluster decreases, possibly due to kinetic arrest from the deeper potential wells of tetrahedral clusters at low volume fractions in which vertices are primarily made out of larger particles.

Unifying Atoms and Colloids near the Glass Transition through Bond-Order Topology

In this combined experimental and simulation study, we utilize bond-order topology to quantitatively match particle volume fraction in mechanically uniformly compressed colloidal suspensions with temperature in atomistic simulations. The obtained mapping temperature is above the dynamical glass transition temperature, indicating that the colloidal systems examined are structurally most like simulated undercooled liquids. Furthermore, the structural mapping procedure offers a unifying framework for quantifying relaxation in arrested colloidal systems.

In situ aggregation and early states of gelation of gold nanoparticle dispersions

The aggregation and onset of gelation of PEGylated gold nanoparticles dispersed in a glycerol-water mixture is studied by small-angle X-ray scattering and X-ray photon correlation spectroscopy. Tracking structural dynamics with sub-ms time resolution over a total experimental time of 8 hours corresponding to a time windows larger than 108 Brownian times and varying the temperature between 298 K and 266 K we can identify three regimes. First, while cooling to 275 K the particles show Brownian motion that slows down due to the increasing viscosity. Second, upon further cooling the static structure changes significantly, indicated by a broad structure factor peak. We attribute this to the formation of aggregates while the dynamics are still dominated by single-particle diffusion. Finally, the relaxation functions become more and more stretched accompanied by an increased slow down of the dynamics. At the same time the structure changes continuously indicating the onset of gelation. Our observations further suggest that the colloidal aggregation and gelation is characterized first by structural changes with a subsequent slowing down of the systems dynamics. The analysis also reveals that the details of the gelation process and the gel structure strongly depend on the thickness of the PEG-coating of the gold nanoparticles.

Hydrodynamic lubrication in colloidal gels

Colloidal gels are elasto-plastic materials composed of an out-of-equilibrium, self-assembled network of micron-sized (solid) particles suspended in a fluid. Recent work has shown that far-field hydrodynamic interactions do not change gel structure, only the rate at which the network forms and ages. However, during gel formation, the interplay between short-ranged attractions leading to gelation and equally short-ranged hydrodynamic lubrication interactions remains poorly understood. Here, we therefore study gelation using a range of hydrodynamic descriptions: from single-body (Brownian Dynamics), to pairwise (Rotne-Prager-Yamakawa), to (non-)lubrication-corrected many-body (Stokesian Dynamics). We confirm the current understanding informed by simulations accurate in the far-field. Yet, we find that accounting for lubrication can strongly impact structure at low colloid volume fraction. Counterintuitively, strongly dissipative lubrication interactions also accelerate the aging of a gel, irrespective of colloid volume fraction. Both elements can be explained by lubrication forces facilitating collective dynamics and therefore phase-separation. Our findings indicate that despite the computational cost, lubricated hydrodynamic modeling with many-body far-field interactions is needed to accurately capture the evolution of the gel structure.

Structuring colloidal gels via micro-bubble oscillations

Locally (re)structuring colloidal gels - micron-sized particles forming a connected network with arrested dynamics - can enable precise tuning of the micromechanical and -rheological properties of the system. A recent experimental study [B. Saint-Michel, G. Petekidis, and V. Garbin, Soft Matter, 2022, 18, 2092] showed that local ordering can be rapidly induced by acoustically modulating an embedded microbubble. Here, we perform Brownian dynamics simulations to understand the mechanical effect of an oscillating microbubble on the next-to-bubble structure of the embedding colloidal gel. Our simulations reveal hexagonal-close-packed structures over a range that is comparable to the amplitude of the oscillations. However, we were unable to reproduce the unexpectedly long-ranged modification of the gel structure - dozens of amplitudes - observed in experiment. This suggests including long-ranged effects, such as fluid flow, should be considered in future computational work.

Hydrodynamic stability criterion for colloidal gelation under gravity

Attractive colloids diffuse and aggregate to form gels, solidlike particle networks suspended in a fluid. Gravity is known to strongly impact the stability of gels once they are formed. However, its effect on the process of gel formation has seldom been studied. Here, we simulate the effect of gravity on gelation using both Brownian dynamics and a lattice-Boltzmann algorithm that accounts for hydrodynamic interactions. We work in a confined geometry to capture macroscopic, buoyancy-induced flows driven by the density mismatch between fluid and colloids. These flows give rise to a stability criterion for network formation, based on an effective accelerated sedimentation of nascent clusters at low volume fractions that disrupts gelation. Above a critical volume fraction, mechanical strength in the forming gel network dominates the dynamics: the interface between the colloid-rich and colloid-poor region moves downward at an ever-decreasing rate. Finally, we analyze the asymptotic state, the colloidal gel-like sediment, which we find not to be appreciably impacted by the vigorous flows that can occur during the settling of the colloids. Our findings represent the first steps toward understanding how flow during formation affects the life span of colloidal gels.

Impact of granular inclusions on the phase behavior of colloidal gels

Colloidal gels formed from small attractive particles are commonly used in formulations to keep larger components in suspension. Despite extensive work characterising unfilled gels, little is known about how the larger inclusions alter the phase behavior and microstructure of the colloidal system. Here we use numerical simulations to examine how larger 'granular' particles can alter the gel transition phase boundaries. We find two distinct regimes depending on both the filler size and native gel structure: a 'passive' regime where the filler fits into already-present voids, giving little change in the transition, and an 'active' regime where the filler no longer fits in these voids and instead perturbs the native structure. In this second regime the phase boundary is controlled by an effective colloidal volume fraction given by the available free volume.

Modeling colloidal interactions that predict equilibrium and non-equilibrium states

Modulating the interaction potential between colloids suspended in a fluid can trigger equilibrium phase transitions as well as formation of non-equilibrium 'arrested states' such as gels and glasses. Faithful representation of such interactions are essential for using simulation to interrogate the microscopic details of non-equilibrium behavior, and for extrapolating observations to new regions of phase space that are difficult to explore in experiment. Although the extended law of corresponding states predicts equilibrium phases for systems with short-ranged interactions, it proves inadequate for equilibrium predictions of systems with longer-ranged interactions, and for predicting non-equilibrium phenomena in systems with either short-ranged or long-ranged interactions. These shortcomings highlight the need for new approaches to represent and disambiguate interaction potentials that replicate both equilibrium and non-equilibrium phase behavior. In this work, we use experiments and simulations to study a system with long-ranged thermoresponsive colloidal interactions and explore whether a resolution to this challenge can be found in regions of the phase diagram where temporal effects influence material state. We demonstrate that the conditions for non-equilibrium arrest by colloidal gelation are sensitive to both the shape of the interaction potential and the thermal quench rate. We exploit this sensitivity to propose a kinetics-based algorithm to extract distinct arrest conditions for candidate potentials that accurately selects between potentials that differ in shape but share the same predicted equilibrium structure. The method reveals that each potential has a quantitatively distinct arrest line, providing insight into how the shape of longer-ranged potentials influences the conditions for colloidal gelation.

Direct Imaging of Contacts and Forces in Colloidal Gels

Colloidal dispersions are prized as model systems to understand basic properties of materials, and are central to a wide range of industries from cosmetics to foods to agrichemicals. Among the key developments in using colloids to address challenges in condensed matter is to resolve the particle coordinates in 3D, allowing a level of analysis usually only possible in computer simulation. However in amorphous materials, relating mechanical properties, and failure in particular to microscopic structure remains problematic. Here we address this challenge by studying the contacts and the forces between particles, as well as their positions. To do so, we use a colloidal model system (an emulsion) in which the interparticle forces and local stress can be linked to the microscopic structure. We demonstrate the potential of our method to reveal insights into the failure mechanisms of soft amorphous solids by determining local stress in a colloidal gel. In particular, we identify "force chains" of load--bearing droplets, and local stress anisotropy, and investigate their connection with locally rigid packings of the droplets.

ArGSLab: a tool for analyzing experimental or simulated particle networks

Microscopy and particle-based simulations are both powerful techniques to study aggregated particulate matter such as colloidal gels. The data provided by these techniques often contains information on a wide array of length scales, but structural analysis methods typically focus on the local particle arrangement, even though the data also contains information about the particle network on the mesoscopic length scale. In this paper, we present a MATLAB software package for quantifying mesoscopic network structures in colloidal samples. ArGSLab (Arrested and Gelated Structures Laboratory) extracts a network backbone from the input data, which is in turn transformed into a set of nodes and links for graph theory-based analysis. The routines can process both image stacks from microscopy as well as explicit coordinate data, and thus allows quantitative comparison between simulations and experiments. ArGSLab furthermore enables the accurate analysis of microscopy data where, e.g., an extended point spread function prohibits the resolution of individual particles. We demonstrate the resulting output for example datasets from both microscopy and simulation of colloidal gels, in order to showcase the capability of ArGSLab to quantitatively analyze data from various sources. The freely available software package can be used either with a provided graphical user interface or directly as a MATLAB script.

Diffusion, phase behavior, and gelation in a two-dimensional layer of colloids in osmotic equilibrium with a polymer reservoir

The addition of enough non-adsorbing polymers to an otherwise stable colloidal suspension gives rise to a variety of phase behaviors and kinetic arrest due to the depletion attraction induced between the colloids by the polymers. We report a study of these phenomena in a two-dimensional layer of colloids. The three-dimensional phenomenology of crystal-fluid coexistence is reproduced, but gelation takes a novel form, in which the strands in the gel structure are locally crystalline. We compare our findings with a previous simulation and theory and find substantial agreement.

Controlling Binder Adhesion to Impact Electrode Mesostructures and Transport

The complex three-phase composition of lithium-ion battery electrodes, containing an ion-conducting pore phase, a nanoporous electron-conducting carbon binder domain (CBD) phase, and an active material (AM) phase, provides several avenues of mesostructural engineering to enhance battery performance. We demonstrate a promising strategy for engineering electrode mesostructures by controlling the strength of adhesion between the AM and CBD phases. Using high-fidelity, physics-based colloidal and granular dynamics simulations, we predict that this strategy can provide significant control over electrochemical transport-relevant properties such as ionic conductivity, electronic conductivity, and available AM-electrolyte interface area. Importantly, the proposed strategy could be experimentally realized through surface functionalization of the AM and CBD phases and would be compatible with traditional electrode manufacturing methods.

Using Patchy Particles to Prevent Local Rearrangements in Models of Non-equilibrium Colloidal Gels

Simple models based on isotropic interparticle attractions often fail to capture experimentally observed structures of colloidal gels formed through spinodal decomposition and subsequent arrest: the resulting gels are typically denser and less branched than their experimental counterparts. Here, we simulate gels formed from soft particles with directional attractions ("patchy particles"), designed to inhibit lateral particle rearrangement after aggregation. We directly compare simulated structures with experimental colloidal gels made using soft attractive microgel particles, by employing a "skeletonization" method that reconstructs the three-dimensional backbone from experiment or simulation. We show that including directional attractions with sufficient valency leads to strongly branched structures compared to isotropic models. Furthermore, combining isotropic and directional attractions provides additional control over aggregation kinetics and gel structure. Our results show that the inhibition of lateral particle rearrangements strongly affects the gel topology and is an important effect to consider in computational models of colloidal gels.

Computer simulations of colloidal gels: how hindered particle rotation affects structure and rheology

The effects of particle roughness and short-ranged non-central forces on colloidal gels are studied using computer simulations in which particles experience a sinusoidal variation in energy as they rotate. The number of minima n and energy scale K are the key parameters; for large K and n, particle rotation is strongly hindered, but for small K and n particle rotation is nearly free. A series of systems are simulated and characterized using fractal dimensions, structure factors, coordination number distributions, bond-angle distributions and linear rheology. When particles rotate easily, clusters restructure to favor dense packings. This leads to longer gelation times and gels with strand-like morphology. The elastic moduli of such gels scale as G'∝ω^0.5 at high shear frequencies ω. In contrast, hindered particle rotation inhibits restructuring and leads to rapid gelation and diffuse morphology. Such gels are stiffer, with G'∝ω^0.35. The viscous moduli G'' in the low-barrier and high-barrier regimes scale according to exponents 0.53 and 0.5, respectively. The crossover frequency between elastic and viscous behaviors generally increases with the barrier to rotation. These findings agree qualitatively with some recent experiments on heterogeneously-surface particles and with studies of DLCA-type gels and gels of smooth spheres.

Magic number colloidal clusters as minimum free energy structures

Clusters in systems as diverse as metal atoms, virus proteins, noble gases, and nucleons have properties that depend sensitively on the number of constituent particles. Certain numbers are termed ‘magic’ because they grant the system with closed shells and exceptional stability. To this point, magic number clusters have been exclusively found with attractive interactions as present between atoms. Here we show that magic number clusters exist in a confined soft matter system with negligible interactions. Colloidal particles in an emulsion droplet spontaneously organize into a series of clusters with precisely defined shell structures. Crucially, free energy calculations demonstrate that colloidal clusters with magic numbers possess higher thermodynamic stability than those off magic numbers. A complex kinetic pathway is responsible for the efficiency of this system in finding its minimum free energy configuration. Targeting similar magic number states is a strategy towards unique configurations in finite self-organizing systems across the scales.

Magic number cluster with closed shells and increased stability often result from potential energy minimization between attractive atoms or particles. Here, Wang et al. show that such magic number clusters can also result from entropy maximization in colloidal systems with negligible interactions.