Directed percolation identified as equilibrium pre-transition towards non-equilibrium arrested gel states

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.

Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: phase behavior and structural properties

Competing short-range attractive (SA) and long range repulsive (LR) particle interactions can be used to describe three-dimensional charge-stabilized colloid or protein dispersions at low added salt concentrations, as well as membrane proteins with interaction contributions mediated by lipid molecules. Using Langevin dynamics (LD) simulations, we determine the generalized phase diagram, cluster shapes and size distributions of a generic quasi-two-dimensional (Q2D) dispersion of spherical SALR particles confined to in-plane motion inside a bulk fluid. The SA and LR interaction parts are modelled by a generalized Lennard-Jones potential and a screened Coulomb potential, respectively. The microstructures of the detected equilibrium and non-equilibrium Q2D phases are distinctly different from those observed in three-dimensional (3D) SALR systems, by exhibiting different levels of hexagonal ordering. We discuss a thermodynamic perturbation theory prediction for the metastable binodal line of a reference system of particles with SA interactions only, which in the explored Q2D-SALR phase diagram region separates cluster from non-clustered phases. The transition from the high-temperature (small SA) dispersed fluid (DF) phase to the lower-temperature equilibrium cluster (EC) fluid phase is characterised by a low-wavenumber peak height of the static structure factor (corresponding to a thermal correlation length of about twice the particle diameter) featuring a distinctly smaller value (≈1.4) than in 3D SALR systems. With decreasing temperature (increasing SA), the cluster morphology changes from disk-like shapes in the equilibrium cluster phase, to double-stranded anisotropic hexagonal cluster segments formed in a cluster-percolated (CP) gel-like phase. This transition can be quantified by a hexagonal order parameter distribution function. The mean cluster size and coordination number of particles in the CP phase are insensitive to changes in the attraction strength.

Percolation of nonequilibrium assemblies of colloidal particles in active chiral liquids

The growing interest in the non-equilibrium assembly of colloidal particles in active liquids is driven by the motivation to create novel structures endowed with tunable properties unattainable within the confines of equilibrium systems. Here, we present an experimental investigation of the structural features of colloidal assemblies in active liquids of chiral E. coli. The colloidal particles form dynamic clusters due to the effective interaction mediated by active media. The activity and chirality of the swimmers strongly influence the dynamics and local ordering of colloidal particles, resulting in clusters with persistent rotation, whose structure differs significantly from those in equilibrium systems with attractive interactions, such as colloid-polymer mixtures. Our colloid-bacteria mixture displays several hallmark features of a percolation transition at a critical density, where the clusters span the system size. A closer examination of the critical exponents associated with cluster size distribution, the average cluster size, and the correlation length in the vicinity of the critical density shows deviations from the prediction of the standard continuum percolation model. Therefore, our experiments reveal a richer phase behavior of colloidal assemblies in active liquids.

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.

Dilute gel networks vs. clumpy gels in colloidal systems with a competition between repulsive and attractive interactions

Using Brownian dynamics simulations we study gel-forming colloidal systems. The focus of this article lies on the differences of dense and dilute gel networks in terms of structure formation both on a local and a global level. We apply reduction algorithms and observe that dilute networks and dense gels differ in the way structural properties like the thickness of strands emerge. We also analyze the percolation behavior and find that two different regimes of percolation exist which might be responsible for structural differences. In dilute networks we confirm that solidity is mainly a consequence of pentagonal bipyramids forming in the network. In dense gels, tetrahedral structures also influence solidity.

Self-assembly of colloids with competing interactions confined in spheres

At low temperatures, colloidal particles with short-range attractive and long-range repulsive interactions can form various periodic microphases in bulk. In this paper, we investigate the self-assembly behaviour of colloids with competing interactions under spherical confinement by conducting molecular dynamics simulations. We find that the cluster, mixture, cylindrical, perforated lamellar and lamellar structures can be obtained, but the details of the ordered structures are different from those in bulk systems. Interestingly, the system tends to form more perforated structures when confined in smaller spheres. The mechanism behind this phenomenon is driven by the relationship between the energy of the ordered structures and the bending of the confinement wall, which is different from the mechanism in copolymer systems.

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, Pe_g, 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 Pe_g. 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.

Memory in aging colloidal gels with time-varying attraction

We report a combined rheology, x-ray photon correlation spectroscopy, and modeling study of gel formation and aging in suspensions of nanocolloidal spheres with volume fractions of 0.20 and 0.43 and with a short-range attraction whose strength is tuned by changing temperature. Following a quench from high temperature, where the colloids are essentially hard spheres, to a temperature below the gel point, the suspensions form gels that undergo aging characterized by a steadily increasing elastic shear modulus and slowing, increasingly constrained microscopic dynamics. The aging proceeds at a faster rate for stronger attraction strength. When the attraction strength is suddenly lowered during aging, the gel properties evolve non-monotonically in a manner resembling the Kovacs effect in glasses, in which the modulus decreases and the microscopic dynamics become less constrained for a period before more conventional aging resumes. Eventually, the properties of the gel following the decrease in attraction strength converge to those of a gel that has undergone aging at the lower attraction strength throughout. The time scale of this convergence increases as a power law with the age at which the attraction strength is decreased and decreases exponentially with the magnitude of the change in attraction. A model for gel aging in which particles attach and detach from the gel at rates that depend on their contact number reproduces these trends and reveals that the non-monotonic behavior results from the dispersion in the rates that the populations of particles with different contact number adjust to the new attraction strength.

Response to "Comment on 'Communication: Is directed percolation in colloid-polymer mixtures linked to dynamic arrest?' " [J. Chem. Phys. 148, 241101 (2018)]

None

Comment on "Communication: Is directed percolation in colloid-polymer mixtures linked to dynamic arrest?'

In their communication [J. Chem. Phys. 148, 241101 (2018)] Richard et al. misquote and criticize a work of Kohl et al.~[Nature Communications 7, 11817 (2016)] resting on simulations in a regime of the phase diagram that significantly differ from the one that Kohl et al.~consider. In this comment it is shown that both the effective density as well as the rescaled second virial coefficient indicate that the comparison presented by Richard et al.~is invalid. Therefore, the implications with respect to the work of Kohl et al.~are incorrect. There is no indication for a disagreement between the simulations of Richard et al.~and those of Kohl et al.~and I am confident that upon a consistent comparison and interpretation of the results both works can contribute to a more comprehensive picture of gel-forming systems.

Diminishing Interfacial Turbulence by Colloid-Polymer Electrolyte to Stabilize Zinc Ion Flux for Deep-Cycling Zn Metal Batteries

The fluidity of aqueous electrolytes and undesired H₂ evolution reaction (HER) can cause severe interfacial turbulence in aqueous Zn metal batteries (ZMBs) at deep cycling with high capacities and current densities, which would further perturb ion flux and aggravate Zn dendrite growth. In this study, a colloid-polymer electrolyte (CPE) with special colloidal phase and suppressed HER is designed to diminish interfacial turbulence and boost deep Zn electrochemistry. Density functional theory calculations confirm that the quantitative migratory barriers of Zn²⁺ along the transport pathway in CPE demonstrate much smaller fluctuations compared with normal aqueous electrolyte, indicating that CPE can effectively diminish interfacial disturbance. Benefitting from this, the Zn²⁺ ion flux can be homogenized and deposited evenly on the electrode, which is confirmed by finite element simulation and in situ Raman measurements. Consequently, CPE enables stable operation of Zn//Cu cells even with high capacity (up to 20 mAh cm^-2 ) and current density (up to 100 mA cm^-2 ) and Zn//Na₅ V₁₂ O₃₂ full-cell with N/P ratio as low as 1 (i.e., 100% Zn utilization). It is believed that this strategy opens a brand-new avenue based on CPE toward boosting deep-cycling electrochemistry in ZMBs and even other aqueous energy-storage applications.

Microscopic structural origin behind slowing down of colloidal phase separation approaching gelation

The gelation of colloidal particles interacting through a short-range attraction is widely recognized as a consequence of the dynamic arrest of phase separation into colloid-rich and solvent-rich phases. However, the microscopic origin behind the slowing down and dynamic arrest of phase separation remains elusive. In order to access microscopic structural changes through the entire process of gelation in a continuous fashion, we used core–shell fluorescent colloidal particles, laser scanning confocal microscopy, and a unique experimental protocol that allows us to initiate phase separation instantaneously and gently. Combining these enables us to track the trajectories of individual particles seamlessly during the whole phase-separation process from the early stage to the late arresting stage. We reveal that the enhancement of local packing and the resulting formation of locally stable rigid structures slow down the phase-separation process and arrest it to form a gel with an average coordination number of z = 6–7. This result supports a mechanical perspective on the dynamic arrest of sticky-sphere systems based on the microstructure, replacing conventional explanations based on the macroscopic vitrification of the colloid-rich phase. Our findings illuminate the microscopic mechanisms behind the dynamic arrest of colloidal phase separation, the emergence of mechanical rigidity, and the stability of colloidal gels.

Location of the gel-like boundary in patchy colloidal dispersions: Rigidity percolation, structure, and particle dynamics

During the past decade, there has been a hot debate about the physical mechanisms that determine when a colloidal dispersion approaches the gel transition. However, there is still no consensus on a possible unique route that leads to the conditions for the formation of a gel-like state. Based on gel states identified in experiments, Valadez-Pérez et al. [Phys. Rev. E 88, 060302(R) (2013)PLEEE81539-375510.1103/PhysRevE.88.060302] proposed rigidity percolation as the precursor of colloidal gelation in adhesive hard-sphere dispersions with coordination number 〈n_{b}〉 equal to 2.4. Although this criterion was originally established to describe mechanical transitions in network-forming molecular materials with highly directional interactions, it worked well to explain gel formation in colloidal suspensions with isotropic short-range attractive forces. Recently, this idea has also been used to account for the dynamical arrest experimentally observed in attractive spherocylinders. Then, by assuming that rigidity percolation also drives gelation in spherical colloids interacting with short-ranged and highly directional potentials, we locate the thermodynamic states where gelation seems to occur in dispersions made up of patchy colloids. To check whether the criterion 〈n_{b}〉=2.4 also holds in patchy colloidal systems, we apply the so-called bond-bending analysis to determine the fraction of floppy modes at some percolating clusters. This analysis confirms that the condition 〈n_{b}〉=2.4 is a good approximation to determine those percolating clusters that are either mechanically stable or rigid. Furthermore, our results point out that not all combinations of patches and coverages lead to a gel-like state. Additionally, we systematically study the structure and the cluster size distribution along those thermodynamic states identified as gels. We show that for high coverage values, the structure is very similar for systems that have the same coverage regardless the number or the position of the patches on the particle surface. Finally, by using dynamic Monte Carlo computer simulations, we calculate both the mean-square displacement and the intermediate scattering function at and in the neighborhood of the gel-like states.

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.

Real space analysis of colloidal gels: triumphs, challenges and future directions

Colloidal gels constitute an important class of materials found in many contexts and with a wide range of applications. Yet as matter far from equilibrium, gels exhibit a variety of time-dependent behaviours, which can be perplexing, such as an increase in strength prior to catastrophic failure. Remarkably, such complex phenomena are faithfully captured by an extremely simple model-'sticky spheres'. Here we review progress in our understanding of colloidal gels made through the use of real space analysis and particle resolved studies. We consider the challenges of obtaining a suitable experimental system where the refractive index and density of the colloidal particles is matched to that of the solvent. We review work to obtain a particle-level mechanism for rigidity in gels and the evolution of our understanding of time-dependent behaviour, from early-time aggregation to ageing, before considering the response of colloidal gels to deformation and then move on to more complex systems of anisotropic particles and mixtures. Finally we note some more exotic materials with similar properties.

Elasticity of colloidal gels: structural heterogeneity, floppy modes, and rigidity

Rheological measurements of model colloidal gels reveal that large variations in the shear moduli as colloidal volume-fraction changes are not reflected by simple structural parameters such as the coordination number, which remains almost a constant. We resolve this apparent contradiction by conducting a normal-mode analysis of experimentally measured bond networks of gels of colloidal particles with short-ranged attraction. We find that structural heterogeneity of the gels, which leads to floppy modes and a nonaffine-affine crossover as frequency increases, evolves as a function of the volume fraction and is key to understanding the frequency-dependent elasticity. Without any free parameters, we achieve good qualitative agreement with the measured mechanical response. Furthermore, we achieve universal collapse of the shear moduli through a phenomenological spring-dashpot model that accounts for the interplay between fluid viscosity, particle dissipation, and contributions from the affine and non-affine network deformation.

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.

Structural and dynamical properties of dilute gel networks in colloid-polymer mixtures

The competition of short-ranged depletion attraction and long-ranged repulsion between colloidal particles in colloid-polymer mixtures leads to the formation of heterogeneous gel-like structures. Our special focus will be on the states where the colloids arrange in thin strands that span the whole system and that we will refer to as dilute gel networks. These states occur at low packing fractions for attractions that are stronger than those at both the binodal line of the equilibrium gas-liquid phase separation and the directed percolation transition line. By using Brownian dynamics simulations, we explore the formation, structure, and aging dynamics of dilute gel networks. The essential connections in a dilute gel network are determined by constructing reduced networks. We compare the observed properties to those of clumpy gels or cluster fluids. Our results demonstrate that both the structure and the (often slow) dynamics of the stable or meta-stable heterogeneous states in colloid-polymer mixtures possess distinct features on various length and time scales and thus are richly diverse.

Growth transitions and critical behavior in the non-equilibrium aggregation of short, patchy nanorods

We have carried out Monte Carlo simulations to study the non-equilibrium aggregation of short patchy nanorods in two dimensions. Below a critical value of patch size ([Formula: see text]), the aggregates have finite sizes with small radii of gyration, [Formula: see text]. At [Formula: see text], the average radius of gyration shows a power law increase with time such that [Formula: see text], where [Formula: see text]. Above, [Formula: see text], the aggregates are fractal in nature and their fractal dimension depends on the value of patch size. These morphological differences are due to the fact that below the critical value of patch size ([Formula: see text]), the growth of the clusters is suppressed and the system reaches an 'absorbed state.' Above [Formula: see text], the system reaches an 'active state,' in which the cluster size keeps growing with a fixed rate at long times. Thus, the system encounters a non-equilibrium phase transition. Close to the transition, the growth rate scales as [Formula: see text], where [Formula: see text]. The long-time growth rate varies as [Formula: see text] where [Formula: see text]. These scaling exponents indicate that the transition belongs to the directed percolation universality class. The patchy nanorods also display a threshold patch size ([Formula: see text]), beyond which the long-time growth rate remains constant. We present geometric arguments for the existence of [Formula: see text]. The fractal dimension of the aggregates increases from 1.75, at [Formula: see text], to 1.81, at [Formula: see text]. It remains constant beyond [Formula: see text].

Critical properties of the susceptible-exposed-infected model with correlated temporal disorder

In this paper we study the critical properties of the nonequilibrium phase transition of the susceptible-exposed-infected (SEI) model under the effects of long-range correlated time-varying environmental noise on the Bethe lattice. We show that temporal noise is perturbatively relevant changing the universality class from the (mean-field) dynamical percolation to the exotic infinite-noise universality class of the contact process model. Our analytical results are based on a mapping to the one-dimensional fractional Brownian motion with an absorbing wall and is confirmed by Monte Carlo simulations. Unlike the contact process, our theory also predicts that it is quite difficult to observe the associated active temporal Griffiths phase in the long-time limit. Finally, we also show an equivalence between the infinite-noise and the compact directed percolation universality classes by relating the SEI model in the presence of temporal disorder to the Domany-Kinzel cellular automaton in the limit of compact clusters.

Aging and equilibration in bistable contagion dynamics

We analyze the late-time relaxation dynamics for a general contagion model. In this model, nodes are either active or failed. Active nodes can fail either "spontaneously" at any time or "externally" if their neighborhoods are sufficiently damaged. Failed nodes may always recover spontaneously. At late times, the breaking of time-translation invariance is a necessary condition for physical aging. We observe that time-translational invariance is lost for initial conditions that lie between the basins of attraction of the model's two stable stationary states. Based on corresponding mean-field predictions, we characterize the observed model behavior in terms of a phase diagram spanned by the fractions of spontaneously and externally failed nodes. For the square lattice, the phases in which the dynamics approaches one of the two stable stationary states are not linearly separable due to spatial correlation effects. Our results provide insights into aging and relaxation phenomena that are observable in a model of social contagion processes.

A unique route of colloidal phase separation yields stress-free gels

Phase separation often leads to gelation in soft and biomatter. For colloidal suspensions, we have a consensus that gels form by the dynamical arrest of phase separation. In this gelation, percolation of the phase-separated structure occurs before the dynamical arrest, leading to the generation of mechanical stress in the gel network. Here, we find a previously unrecognized type of gelation in dilute colloidal suspensions, in which percolation occurs after the local dynamical arrest, i.e., the formation of mechanically stable, rigid clusters. Thus, topological percolation generates little mechanical stress, and the resulting gel is almost stress-free when formed. We also show that the selection of these two types of gelation (stressed and stress-free) is determined solely by the volume fraction as long as the interaction is short-ranged. This universal classification of gelation of particulate systems may have a substantial impact on material and biological science.

Seaweed's Bioactive Candidate Compounds to Food Industry and Global Food Security

The world population is continuously growing, so it is important to keep producing food in a sustainable way, especially in a way that is nutritious and in a sufficient quantity to overcome global needs. Seaweed grows, and can be cultivated, in seawater and generally does not compete for arable land and freshwater. Thus, the coastal areas of the planet are the most suitable for seaweed production, which can be an alternative to traditional agriculture and can thus contribute to a reduced carbon footprint. There are evolving studies that characterize seaweed's nutritional value and policies that recognize them as food, and identify the potential benefits and negative factors that may be produced or accumulated by seaweed, which are, or can be, dangerous for human health. Seaweeds have a high nutritional value along with a low caloric input and with the presence of fibers, proteins, omega 3 and 6 unsaturated fatty acids, vitamins, and minerals. Moreover, several seaweed sub-products have interesting features to the food industry. Therefore, the focus of this review is in the performance of seaweed as a potential alternative and as a safe food source. Here described is the nutritional value and concerns relating to seaweed consumption, and also how seaweed-derived compounds are already commercially explored and available in the food industry and the usage restrictions to safeguard them as safe food additives for human consumption.

A computational study of the behavior of colloidal gel networks at low volume fraction

Colloidal gel networks appear in different scientific and industrial applications because of their unique properties. Molecular dynamics simulations could reveal the relation between molecular level and macroscopic properties of these systems. Nevertheless, the predictions of numerical simulations might depend on the specific form and parameters of interaction potentials. In this paper, a new effective interaction potential is used for characterizing the mechanical behavior of low volume fraction colloidal gels under large shear deformation. The findings are compared with those obtained from other available forms of interaction potentials in order to determine gel characteristics that are interaction potential independent. Furthermore, the macroscopic stress-strain behavior is discussed in terms of the behavior of different terms of the proposed interaction potential. The correlation between the stretch of interparticle bonds and their alignment in the direction of the maximum principal stress is also computed in order to provide microscopic explanations for the initial strain softening behavior. It is concluded that, in addition to topology, local mechanical interactions between colloidal particles are important in defining the mechanical response of soft gels.

Chemical vs. mechanical microstructure evolution in drying colloid and polymer coatings

Colloidal based films have been widely developed for a wide range of applications including chemical and electrical barrier coatings, photonic materials, biomaterials, and pharmaceutical oral drug delivery. Many previous studies investigate methods to generate uniformity or desired stratification of the final components with a desired microstructure. Few studies have been able to investigate this microstructure in-situ during drying. This experimental study directly tracks fluorescent colloids that are either stable in suspension or have attractive interactions during the drying process using high speed laser scanning confocal microscopy to obtain details of microstructural evolution during drying. The colloidal microstructure in stable suspensions evolves continuously during drying. Microstructures in these systems have a signature Voronoi polyhedra distribution that is defined by lognormal curve having a constant standard deviation that only depends on its chemical composition. Those formulations having strongly attractive constituents have microstructure that is heterogeneous and non-monotonic due to the mechanics associated with internal convection and capillary forces. Toward the end of drying, the influence of the mode of microstructure rearrangements remains evident.

Dynamic arrest of adhesive hard rod dispersions

The phenomenon of dynamic arrest, more commonly referred to as gel and glass formation, originates as particle motion slows significantly. Current understanding of gels and glasses stems primarily from dispersions of spherical particles, but much less is known about how particle shape affects dynamic arrest transitions. To better understand the effects of particle shape anisotropy on gel and glass formation, we systematically measure the rheology, particle dynamics, and static microstructure of thermoreversible colloidal dispersions of adhesive hard rods (AHR). First, the dynamic arrest transitions are mapped as a function of temperature T, aspect ratio L/D≈ 3 to 7, and volume fraction φ≈ 0.1 to 0.5. The critical gel temperature T_gel and glass volume fraction φ_g are determined from the particle dynamics and rheology. Second, an effective orientation-averaged, short-range attraction between rods is quantified from small-angle scattering measurements and characterized by a reduced temperature τ. Similar τ is found at low rod concentrations, indicating that rod gelation occurs at similar effective attraction strength independent of L/D. Monte Carlo simulations reveal a similar convergence in τ when rods cluster and percolate with an average bond coordination number 〈n_c〉≈ 2.4, supporting the link between physical gelation and rigidity percolation. Lastly, AHR results are mapped onto a dimensionless state diagram to compare with previous predictions of attraction-driven gels, repulsion-driven glasses, and liquid crystal phases.

Direct link between mechanical stability in gels and percolation of isostatic particles

Colloidal gels have unique mechanical and transport properties that stem from their bicontinuous nature, in which a colloidal network is intertwined with a viscous solvent, and have found numerous applications in foods, cosmetics, and construction materials and for medical applications, such as cartilage replacements. So far, our understanding of the process of colloidal gelation is limited to long-time dynamical effects, where gelation is viewed as a phase separation process interrupted by the glass transition. However, this purely out-of-equilibrium thermodynamic picture does not address the emergence of mechanical stability. With confocal microscopy experiments, we reveal that mechanical metastability is reached only after isotropic percolation of locally isostatic environments, establishing a direct link between the load-bearing ability of gels and the isostaticity condition. Our work suggests an operative description of gels based on mechanical equilibrium and isostaticity, providing the physical basis for the stability and rheology of these materials.

Temperature dependence of the transition packing fraction of thermal jamming in a harmonic soft sphere system

The glassy dynamics of soft harmonic spheres are often mapped onto the dynamics of hard spheres by considering an effective diameter for the soft particles and therefore an effective packing fraction. While in this approach the thermal fluctuations within valleys of the energy landscape are covered, the crossing of energy barriers from one valley into another usually is neglected. Here we argue-motivated by studies of the glass transition based on explorations of the energy landscape-that the crossing of energy barriers can be attributed by an effective decrease of the glass transition packing fraction with increasing temperature T according to T ^0.2. Furthermore, we reanalyzing data of soft sphere simulations. Since fitting scaling laws to simulation data always allows for some arbitrariness, we cannot prove based on the simulation data that our idea of a shift of the glass transition packing fraction due to barrier crossings is the only possible way to explain the discrepancies that have been observed previously. However, we show that a possible explanation of the simulation data is given by our approach to characterize the dynamics of soft spheres by both, the previously-considered temperature-dependent effective packing fraction due to the increase of the mean overlap between neighboring particles with stronger thermal fluctuations and the newly introduced increase of the glass transition packing with an increasing probability of barrier crossings.

The thermal jamming transition of soft harmonic disks in two dimensions

By exploring the properties of the energy landscape of a bidisperse system of soft harmonic disks in two dimensions we determine the thermal jamming transition. To be specific, we study whether the ground state of the system where the particles do not overlap can be reached within a reasonable time. Starting with random initial configurations, the energy landscape is probed by energy minimization steps as in case of athermal jamming and in addition steps where an energy barrier can be crossed with a small but non-zero probability. For random initial conditions we find that as a function of packing fraction the thermal jamming transition, i.e. the transition from a state where all overlaps can be removed to an effectively non-ergodic state where one cannot get rid of the overlaps, occurs at a packing fraction of [Formula: see text], which is smaller than the transition packing fraction of athermal jamming at [Formula: see text]. Furthermore, we show that the thermal jamming transition is in the universality class of directed percolation and therefore is fundamentally different from the athermal jamming transition.

Reversible Formation of Thermoresponsive Binary Particle Gels with Tunable Structural and Mechanical Properties

We investigate the collective behavior of suspended thermoresponsive microgels that expel solvent and subsequently decrease in size upon heating. Using a binary mixture of differently thermoresponsive microgels, we demonstrate how distinctly different gel structures form, depending on the heating profile used. Confocal laser scanning microscopy (CLSM) imaging shows that slow heating ramps yield a core-shell network through sequential gelation, while fast heating ramps yield a random binary network through homogelation. Here, secondary particles are shown to aggregate in a monolayer fashion upon the first gel, which can be qualitatively reproduced through Brownian dynamics simulations using a model based on a temperature-dependent interaction potential incorporating steric repulsion and van der Waals attraction. Through oscillatory rheology it is shown that secondary microgel deposition enhances the structural integrity of the previously formed single species gel, and the final structure exhibits higher elastic and loss moduli than its compositionally identical homogelled counterpart. Furthermore, we demonstrate that aging processes in the scaffold before secondary microgel deposition govern the final structural properties of the bigel, which allows a detailed control over these properties. Our results thus demonstrate how the temperature profile can be used to finely control the structural and mechanical properties of these highly tunable materials.

Coupling between criticality and gelation in "sticky" spheres: a structural analysis

We combine experiments and simulations to study the link between criticality and gelation in sticky spheres. We employ confocal microscopy to image colloid-polymer mixtures and Monte Carlo simulations of the square-well (SW) potential as a reference model. To this end, we map our experimental samples onto the SW model. We find an excellent structural agreement between experiments and simulations, both for locally favored structures at the single particle level and large-scale fluctuations at criticality. We follow in detail the rapid structural change in the critical fluid when approaching the gas-liquid binodal and highlight the role of critical density fluctuations for this structural crossover. Our results link the arrested spinodal decomposition to long-lived energetically favored structures, which grow even away from the binodal due to the critical scaling of the bulk correlation length and static susceptibility.

Vitrification and gelation in sticky spheres

Glasses and gels are the two dynamically arrested, disordered states of matter. Despite their importance, their similarities and differences remain elusive, especially at high density, where until now it has been impossible to distinguish them. We identify dynamical and structural signatures which distinguish the gel and glass transitions in a colloidal model system of hard and "sticky" spheres. It has been suggested that "spinodal" gelation is initiated by gas-liquid viscoelastic phase separation to a bicontinuous network and the resulting densification leads to vitrification of the colloid-rich phase, but whether this phase has sufficient density for arrest is unclear [M. A. Miller and D. Frenkel, Phys. Rev. Lett. 90, 135702 (2003) and P. J. Lu et al., Nature 435, 499-504 (2008)]. Moreover alternative mechanisms for arrest involving percolation have been proposed [A. P. R. Eberle et al., Phys. Rev. Lett. 106, 105704 (2011)]. Here we resolve these outstanding questions, beginning by determining the phase diagram. This, along with demonstrating that percolation plays no role in controlling the dynamics of our system, enables us to confirm spinodal decomposition as the mechanism for gelation. We are then able to show that gels can be formed even at much higher densities than previously supposed, at least to a volume fraction of ϕ = 0.59. Far from being networks, these gels apparently resemble glasses but are still clearly distinguished by the "discontinuous" nature of the transition and the resulting rapid solidification, which leads to the formation of inhomogeneous (with small voids) and far-from-equilibrium local structures. This is markedly different from the glass transition, whose continuous nature leads to the formation of homogeneous and locally equilibrated structures. We further reveal that the onset of the attractive glass transition in the form of a supercooled liquid is in fact interrupted by gelation. Our findings provide a general thermodynamic, dynamic, and structural basis upon which we can distinguish gelation from vitrification.

Hunting mermaids in real space: known knowns, known unknowns and unknown unknowns

We review efforts to realise so-called mermaid (or short-ranged attraction/long ranged repulsion) interactions in 3d real space. The repulsive and attractive contributions to these interactions in charged colloids and colloid-polymer mixtures, may be accurately realised, by comparing particle-resolved studies with colloids to computer simulation. However, when we review work where these interactions have been combined, despite early indications of behaviour consistent with predictions, closer analysis reveals that in the non-aqueous systems used for particle-resolved studies, the idea of summing the attractive and repulsive components leads to wild deviations with experiment. We suggest that the origin lies in the weak ion dissociation in these systems with low dielectric constant solvents. Ultimately this leads even to non-centro-symmetric interactions and a new level of complexity in these systems.

Ergodicity breaking transition in a glassy soft sphere system at small but non-zero temperatures

While the glass transition at non-zero temperature seems to be hard to access for experimental, theoretical, or simulation studies, jamming at zero temperature has been studied in great detail. Motivated by the exploration of the energy landscape that has been successfully used to investigate athermal jamming, we introduce a new method that includes the possibility of the thermally excited crossing of energy barriers. We then determine whether the ground state configurations of a soft sphere system are accessible or not and as a consequence whether the system is ergodic or effectively non-ergodic. Interestingly, we find an transition where the system becomes effectively non-ergodic if the density is increased. The transition density in the limit of small but non-zero temperatures is independent of temperature and below the transition density of athermal jamming. This confirms recent computer simulation studies where athermal jamming occurs deep inside the glass phase. In addition, we show that the ergodicity breaking transition is in the universality class of directed percolation. Therefore, our approach not only makes the transition from an ergodic to an effectively non-ergodic systems easily accessible and helps to reveal its universality class but also shows that it is fundamentally different from athermal jamming.

Lattice animals in diffusion limited binary colloidal system

In a soft matter system, controlling the structure of the amorphous materials has been a key challenge. In this work, we have modeled irreversible diffusion limited cluster aggregation of binary colloids, which serves as a model for chemical gels. Irreversible aggregation of binary colloidal particles leads to the formation of a percolating cluster of one species or both species which are also called bigels. Before the formation of the percolating cluster, the system forms a self-similar structure defined by a fractal dimension. For a one component system when the volume fraction is very small, the clusters are far apart from each other and the system has a fractal dimension of 1.8. Contrary to this, we will show that for the binary system, we observe the presence of lattice animals which has a fractal dimension of 2 irrespective of the volume fraction. When the clusters start inter-penetrating, we observe a fractal dimension of 2.5, which is the same as in the case of the one component system. We were also able to predict the formation of bigels using a simple inequality relation. We have also shown that the growth of clusters follows the kinetic equations introduced by Smoluchowski for diffusion limited cluster aggregation. We will also show that the chemical distance of a cluster in the flocculation regime will follow the same scaling law as predicted for the lattice animals. Further, we will also show that irreversible binary aggregation comes under the universality class of the percolation theory.

Shear-induced slab-like domains in a directed percolated colloidal gel

We explore the structural changes of a gel-forming colloid polymer mixture under shear by employing Brownian dynamics simulations of a colloidal system with short-ranged attractive depletion interaction in a linear flow profile. While the structure of unpercolated systems changes only slightly under shearing, we discover the formation of slab-like clusters in sheared directed percolated gel networks that are confined between two walls. These gel-slabs are stable over a long time and seem to be related to the syneresis phenomena that can be observed in directed percolated colloidal gels. Only at large shear strength the slabs are destroyed and a homogeneous state with many unbounded particles can be observed. We also quantitatively analyze our results by determining void volumes.

Phase behavior of colloid-polymer depletion mixtures with unary or binary depletants

Adding depletants to a colloidal suspension induces an attractive interparticle interaction that can be tuned to obtain desired structures or to probe phase behavior. When the depletant is not uniform in size, however, both the range and strength of the attraction become difficult to predict and hence control. We investigated the effects of depletant bidispersity on the non-equilibrium phase behavior of colloid-polymer mixtures. We added unary or binary mixtures of polystyrene as the depletant to suspensions of charged poly(methyl methacrylate) particles. The structure and dynamics of the particles were compared over three sets of samples with various mixtures of two different polystyrenes whose size varied by an order of magnitude. The structure and dynamics were nearly independent of depletant dispersity if the polymer concentration was represented as a sum of normalized concentrations of each species. Near the transition region between a fluid of clusters and an interconnected gel at intermediate volume fractions, partitioning of polymers in a binary mixture into colloid-rich and polymer-rich phase leads to a slightly different gelation pathway.

Why patchy diffusion-limited aggregation belongs to the directed-percolation universality class

We present a possible link between nonequilibrium phase transition observed in patchy diffusion-limited aggregation (DLA) [M. J. Kartha and A. Sayeed, Phys. Lett. A 380, 2791 (2016)10.1016/j.physleta.2016.06.036] and directed bond percolation (DP). A system of directed percolation with patchy particles (patchy DP) in which the bond connectivity is established depending on patch size p is analyzed. It is observed that patchy DP starting from a single seed shows a nonequilibrium phase transition. Below a critical value of the patch size p_{c}, the system reaches an absorbing state above which is a fluctuating active state as observed in the DP system. The value of this nonuniversal parameter p_{c} is observed to be slightly higher than the value observed in patchy DLA. Close to the critical value, the order parameter P(∞)∼(p-p_{c})^{β} where β=0.272±0.010, which is consistent with the directed-percolation universality class. Therefore the intrinsic nature of patchy DP is responsible for the phase transition in patchy DLA. This study reveals that the estimated critical value of patch size p_{c}=0.80625±0.00020 in patchy DP is different from the critical bond probability p_{c}=0.6447 in the DP system. This elucidates that the bond probability in DP is not equivalent to the patch probability of a particular site. Our work also gives an insight into the problem related with formation of an extended network of pentagon subunits in connection with the virus capsid.

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.