Intermittent dynamics and logarithmic domain growth during the spinodal decomposition of a glass-forming liquid

Density and inertia effects on two-dimensional active semiflexible filament suspensions

We examine the influence of density on the transition between chain and spiral structures in planar assemblies of active semiflexible filaments, utilizing detailed numerical simulations. We focus on how increased density, and higher Péclet numbers, affect the activity-induced transition spiral state in a semiflexible, self-avoiding active chain. Our findings show that increasing the density causes the spiral state to break up, reverting to a motile chain-like shape. This results in a density-dependent reentrant phase transition from spirals back to open chains. We attribute this phenomenon to an inertial effect observed at the single polymer level, where increased persistence length due to inertia has been shown in recent three-dimensional studies to cause polymers to open up. Our two-dimensional simulations further reveal that a reduction in the damping coefficient leads to partial unwinding of the spirals, forming longer arms. In suspension, interactions among these extended arms can trigger a complete unwinding of the spirals, driven by the combined effects of density and inertia.

The impact of viscosity asymmetry on phase separating binary mixtures with suspended colloids

The introduction of neutrally wetting colloidal particles into coarsening binary fluids is known to arrest the dynamics of the phase separation, as the colloids tend to be captured by the growing interfaces to reduce the free energy of the system. This phenomenon has often been studied in systems with symmetric fluid viscosities. In this study, we investigate the behavior of colloidal particles introduced into asymmetric binary fluids with a viscosity contrast. Our results show that due to the broken symmetry the colloidal particles more easily escape from the interface towards the more viscous fluid, which reduces the lifetime of the jammed phase. Moreover, the presence of colloidal particles near the interfaces promotes the formation of micro-droplets with typical sizes comparable to the colloids.

Necking and failure of a particulate gel strand: signatures of yielding on different length scales

"Sticky" spheres with a short-ranged attraction are a basic model of a wide range of materials from the atomic to the granular length scale. Among the complex phenomena exhibited by sticky spheres is the formation of far-from-equilibrium dynamically arrested networks which comprise "strands" of densely packed particles. The aging and failure of such gels under load is a remarkably challenging problem, given the simplicity of the model, as it involves multiple length- and time-scales, making a single approach ineffective. Here we tackle this challenge by addressing the failure of a single strand with a combination of methods. We study the mechanical response of a single strand of a model gel-former to deformation, both numerically and analytically. Under elongation, the strand breaks by a necking instability. We analyse this behaviour at three different length scales: a rheological continuum model of the whole strand; a microscopic analysis of the particle structure and dynamics; and the local stress tensor. Combining these different approaches gives a coherent picture of the necking and failure. The strand has an amorphous local structure and has large residual stresses from its initialisation. We find that neck formation is associated with increased plastic flow, a reduction in the stability of the local structure, and a reduction in the residual stresses; this indicates that the system loses its solid character and starts to behave more like a viscous fluid. These results will inform the development of more detailed models that incorporate the heterogeneous network structure of particulate gels.

Mechanical Slowing Down of Network-Forming Phase Separation of Polymer Solutions

Phase separation is a fundamental phenomenon leading to spatially heterogeneous material distribution, which is critical in nature, biology, material science, and industry. In ordinary phase separation, the minority phase always forms droplets. Contrary to this common belief, even the minority phase can form a network structure in viscoelastic phase separation (VPS). VPS can occur in any mixture with significant mobility differences between their components and is highly relevant to soft matter and biomatter. In contrast to classical phase separation, experiments have shown that VPS in polymer solutions lacks self-similar coarsening, resulting in the absence of a domain-coarsening scaling law. However, the underlying microscopic mechanism of this behavior remains unknown. To this end, we perform fluid particle dynamics simulations of bead-spring polymers, incorporating many-body hydrodynamic interactions between polymers through a solvent. We discover that polymers in the dense-network-forming phase are stretched and store elastic energy when the deformation speed exceeds the polymer dynamics. This self-generated viscoelastic stress mechanically interferes with phase separation and slows its dynamics, disrupting self-similar growth. We also highlight the essential role of many-body hydrodynamic interactions in VPS. The implications of our findings may hold importance in areas such as biological phase separation, porous material formation, and other fields where network structures play a pivotal role.

Two-step devitrification of ultrastable glasses

The discovery of ultrastable glasses raises novel challenges about glassy systems. Recent experiments studied the macroscopic devitrification of ultrastable glasses into liquids upon heating but lacked microscopic resolution. We use molecular dynamics simulations to analyze the kinetics of this transformation. In the most stable systems, devitrification occurs after a very large time, but the liquid emerges in two steps. At short times, we observe the rare nucleation and slow growth of isolated droplets containing a liquid maintained under pressure by the rigidity of the surrounding glass. At large times, pressure is released after the droplets coalesce into large domains, which accelerates devitrification. This two-step process produces pronounced deviations from the classical Avrami kinetics and explains the emergence of a giant lengthscale characterizing the devitrification of bulk ultrastable glasses. Our study elucidates the nonequilibrium kinetics of glasses following a large temperature jump, which differs from both equilibrium relaxation and aging dynamics, and will guide future experimental studies.

The interplay between chemo-phoretic interactions and crowding in active colloids

Many motile microorganisms communicate with each other and their environments via chemical signaling which leads to long-range interactions mediated by self-generated chemical gradients. However, consequences of the interplay between crowding and chemotactic interactions on their collective behavior remain poorly understood. In this work, we use Brownian dynamics simulations to investigate the effect of packing fraction on the formation of non-equilibrium structures in a monolayer of diffusiophoretic self-propelled colloids as a model for chemically active particles. Focusing on the case when a chemical field induces attractive positional and repulsive orientational interactions, we explore dynamical steady-states of active colloids of varying packing fractions and degrees of motility. In addition to collapsed, active gas, and dynamical clustering steady-states reported earlier for low packing fractions, a new phase-separated state emerges. The phase separation results from a competition between long-range diffusiophoretic interactions and motility and is observed at moderate activities and a wide range of packing fractions. Our analysis suggests that the fraction of particles in the largest cluster is a suitable order parameter for capturing the transition from an active gas and dynamical clustering states to a phase-separated state.

Structure and Mechanical Properties of a Porous Polymer Material via Molecular Dynamics Simulations

We characterize, using molecular dynamics simulations, the structure and mechanical response of a porous glassy system, obtained via arrested phase separation of a model polymer melt. In the absence of external driving, coarsening dynamics, with power-law time dependence, controls the slow structural evolution, in agreement with what was reported for other phase-separating systems. The mechanical response was investigated in athermal quasi-static conditions. In the elastic regime, low values for the Young's and shear modulus were found, as compared to dense glassy systems, which originate from the porous structure. For large deformations, stress-strain curves show a highly intermittent behavior, with avalanches of plastic events. The stress-drop distribution is characterized exploring a large set of parameters. This work goes beyond the previous numerical studies on atomic porous materials, as it first examines the role of chain connectivity in the elastic and plastic responses of materials of this type.

Gas-liquid phase separation at zero temperature: mechanical interpretation and implications for gelation

The relationship between glasses and gels has been intensely debated for decades; however, the transition between these two phases remains elusive. To investigate a gel formation process in the zero-temperature limit and its relation to the glass phase, we conducted numerical experiments on athermal quasistatic decompression. During decompression, the system experiences a cavitation event similar to phase separation and this is a gelation process at zero temperature. A normal mode analysis revealed that the phase separation is signaled by the vanishing of the lowest eigenenergy, similar to plastic events of glasses under shear. One primary difference from the shear-induced plasticity is that the vanishing mode experiences a qualitative change in its spatial energy distribution at the phase separation point. These findings enable us to define the glass-gel phase boundary based on mechanics.

Nucleation causes an actin network to fragment into multiple high-density domains

Actin networks rely on nucleation mechanisms to generate new filaments because spontaneous nucleation is kinetically disfavored. Branching nucleation of actin filaments by actin-related protein (Arp2/3), in particular, is critical for actin self-organization. In this study, we use the simulation platform for active matter MEDYAN to generate 2000 s long stochastic trajectories of actin networks, under varying Arp2/3 concentrations, in reaction volumes of biologically meaningful size (>20 μm3). We find that the dynamics of Arp2/3 increase the abundance of short filaments and increases network treadmilling rate. By analyzing the density fields of F-actin, we find that at low Arp2/3 concentrations, F-actin is organized into a single connected and contractile domain, while at elevated Arp2/3 levels (10 nM and above), such high-density actin domains fragment into smaller domains spanning a wide range of volumes. These fragmented domains are extremely dynamic, continuously merging and splitting, owing to the high treadmilling rate of the underlying actin network. Treating the domain dynamics as a drift-diffusion process, we find that the fragmented state is stochastically favored, and the network state slowly drifts toward the fragmented state with considerable diffusion (variability) in the number of domains. We suggest that tuning the Arp2/3 concentration enables cells to transition from a globally coherent cytoskeleton, whose response involves the entire cytoplasmic network, to a fragmented cytoskeleton, where domains can respond independently to locally varying signals.

Gelation Dynamics upon Pressure-Induced Liquid-Liquid Phase Separation in a Water-Lysozyme Solution

Employing X-ray photon correlation spectroscopy, we measure the kinetics and dynamics of a pressure-induced liquid-liquid phase separation (LLPS) in a water-lysozyme solution. Scattering invariants and kinetic information provide evidence that the system reaches the phase boundary upon pressure-induced LLPS with no sign of arrest. The coarsening slows down with increasing quench depths. The g₂ functions display a two-step decay with a gradually increasing nonergodicity parameter typical for gelation. We observe fast superdiffusive (γ ≥ 3/2) and slow subdiffusive (γ < 0.6) motion associated with fast viscoelastic fluctuations of the network and a slow viscous coarsening process, respectively. The dynamics age linearly with time τ ∝ t_w, and we observe the onset of viscoelastic relaxation for deeper quenches. Our results suggest that the protein solution gels upon reaching the phase boundary.

Phonon transport properties of particulate physical gels

Particulate physical gels are sparse, low-density amorphous materials in which clusters of glasses are connected to form a heterogeneous network structure. This structure is characterized by two length scales, ξ s and ξ G: ξ s measures the length of heterogeneities in the network structure and ξ G is the size of glassy clusters. Accordingly, the vibrational states (eigenmodes) of such a material also exhibit a multiscale nature with two characteristic frequencies, [Formula: see text] and ω G, which are associated with ξ s and ξ G, respectively: (i) phonon-like vibrations in the homogeneous medium at [Formula: see text], (ii) phonon-like vibrations in the heterogeneous medium at [Formula: see text], and (iii) disordered vibrations in the glassy clusters at ω > ω G. Here, we demonstrate that the multiscale characteristics seen in the static structures and vibrational states also extend to the phonon transport properties. Phonon transport exhibits two distinct crossovers at frequencies ω* and ω G (or at wavenumbers of [Formula: see text] and [Formula: see text]). In particular, both transverse and longitudinal phonons cross over between Rayleigh scattering at [Formula: see text] and diffusive damping at [Formula: see text]. Remarkably, the Ioffe–Regel limit is located at the very low frequency of ω*. Thus, phonon transport is localized above ω*, even where phonon-like vibrational states persist. This markedly strong scattering behavior is caused by the sparse, porous structure of the gel.

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.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

Mechanical properties and pore size distribution in athermal porous glasses

We study the mechanical properties and pore structure in a three-dimensional molecular dynamics model of porous glass under athermal quasistatic shear. The vitreous samples are prepared by rapid thermal quenching from a high-temperature molten state. The pore structures form via solid-gas phase separation. The quiescent samples exhibit a wide range of pore topography, from inter-connected pore networks to randomly distributed compact pores depending on the material density. We find that the shear modulus strongly depends on the density and porosity. Under mechanical loading, the pore structure rearranges which is reflected in the pore size distribution function. Our results show that with increase in strain the distribution widens as the adjacent pores coalesce and form larger pores. We also propose a universal scaling law for the pore size distribution function which offers excellent data collapse for highly porous materials in the undeformed case. From the data scaling, we identify a critical density that can be attributed to the transition point from a porous-type to bulk-type material. The validity of the scaling law under finite deformation is also analyzed.

Waiting-time dependent non-equilibrium phase diagram of simple glass- and gel-forming liquids

Under numerous circumstances, many soft and hard materials are present in a puzzling wealth of non-equilibrium amorphous states, whose properties are not stationary and depend on preparation. They are often summarized in unconventional "phase diagrams" that exhibit new "phases" and/or "transitions" in which time, however, is an essential variable. This work proposes a solution to the problem of theoretically defining and predicting these non-equilibrium phases and their time-evolving phase diagrams, given the underlying molecular interactions. We demonstrate that these non-equilibrium phases and the corresponding non-stationary (i.e., aging) phase diagrams can indeed be defined and predicted using the kinetic perspective of a novel non-equilibrium statistical mechanical theory of irreversible processes. This is illustrated with the theoretical description of the transient process of dynamic arrest into non-equilibrium amorphous solid phases of an instantaneously quenched simple model fluid involving repulsive hard-sphere plus attractive square well pair interactions.

Power-law coarsening in network-forming phase separation governed by mechanical relaxation

A space-spanning network structure is a basic morphology in phase separation of soft and biomatter, alongside a droplet one. Despite its fundamental and industrial importance, the physical principle underlying such network-forming phase separation remains elusive. Here, we study the network coarsening during gas-liquid-type phase separation of colloidal suspensions and pure fluids, by hydrodynamic and molecular dynamics simulations, respectively. For both, the detailed analyses of the pore sizes and strain field reveal the self-similar network coarsening and the unconventional power-law growth more than a decade according to ℓ ∝ t1/2, where ℓ is the characteristic pore size and t is the elapsed time. We find that phase-separation dynamics is controlled by mechanical relaxation of the network-forming dense phase, whose limiting process is permeation flow of the solvent for colloidal suspensions and heat transport for pure fluids. This universal coarsening law would contribute to the fundamental physical understanding of network-forming phase separation.

Interplay between Glass Formation and Liquid-Liquid Phase Separation Revealed by the Scattering Invariant

The interplay of the glass transition with liquid-liquid phase separation (LLPS) is a subject of intense debate. We use the scattering invariant Q to probe how approaching the glass transition affects the shape of LLPS boundaries in the temperature/volume fraction plane. Two protein systems featuring kinetic arrest with a lower and an upper critical solution temperature phase behavior, respectively, are studied varying the quench depth. Using Q we noninvasively identify system-dependent differences for the effect of glass formation on the LLPS boundary. The glassy dense phase appears to enter the coexistence region for the albumin-YCl₃ system, whereas it follows the equilibrium binodal for the γ-globulin-PEG system.

Relaxation in a phase-separating two-dimensional active matter system with alignment interaction

Via computer simulations, we study kinetics of pattern formation in a two-dimensional active matter system. Self-propulsion in our model is incorporated via the Vicsek-like activity, i.e., particles have the tendency of aligning their velocities with the average directions of motion of their neighbors. In addition to this dynamic or active interaction, there exists passive inter-particle interaction in the model for which we have chosen the standard Lennard-Jones form. Following quenches of homogeneous configurations to a point deep inside the region of coexistence between high and low density phases, as the systems exhibit formation and evolution of particle-rich clusters, we investigate properties related to the morphology, growth, and aging. A focus of our study is on the understanding of the effects of structure on growth and aging. To quantify the latter, we use the two-time order-parameter autocorrelation function. This correlation, as well as the growth, is observed to follow power-law time dependence, qualitatively similar to the scaling behavior reported for passive systems. The values of the exponents have been estimated and discussed by comparing with the previously obtained numbers for other dimensions as well as with the new results for the passive limit of the considered model. We have also presented results on the effects of temperature on the activity mediated phase separation.

Coarsening Kinetics of Complex Macromolecular Architectures in Bad Solvent

This study reports a general scenario for the out-of-equilibrium features of collapsing polymeric architectures. We use molecular dynamics simulations to characterize the coarsening kinetics, in bad solvent, for several macromolecular systems with an increasing degree of structural complexity. In particular, we focus on: flexible and semiflexible polymer chains, star polymers with 3 and 12 arms, and microgels with both ordered and disordered networks. Starting from a powerful analogy with critical phenomena, we construct a density field representation that removes fast fluctuations and provides a consistent characterization of the domain growth. Our results indicate that the coarsening kinetics presents a scaling behaviour that is independent of the solvent quality parameter, in analogy to the time-temperature superposition principle. Interestingly, the domain growth in time follows a power-law behaviour that is approximately independent of the architecture for all the flexible systems; while it is steeper for the semiflexible chains. Nevertheless, the fractal nature of the dense regions emerging during the collapse exhibits the same scaling behaviour for all the macromolecules. This suggests that the faster growing length scale in the semiflexible chains originates just from a faster mass diffusion along the chain contour, induced by the local stiffness. The decay of the dynamic correlations displays scaling behavior with the growing length scale of the system, which is a characteristic signature in coarsening phenomena.

Anisotropic viscoelastic phase separation in polydisperse hard rods leads to nonsticky gelation

Spinodal demixing into two phases having very different viscosities leads to viscoelastic networks-i.e., gels-usually as a result of attractive particle interactions. Here, however, we demonstrate demixing in a colloidal system of polydisperse, rod-like clay particles that is driven by particle repulsions instead. One of the phases is a nematic liquid crystal with a highly anisotropic viscosity, allowing flow along the director, but suppressing it in other directions. This phase coexists with a dilute isotropic phase. Real-space analysis and molecular-dynamics simulations both reveal a long-lived network structure that is locally anisotropic, yet macroscopically isotropic. We show that our system exhibits the characteristics of colloidal gelation, leading to nonsticky gels.

Arrested dynamics of the dipolar hard sphere model

We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (soft core) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: at low (η = 0.2) and intermediate (η = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations (η≥ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid.

Aging phenomena during phase separation in fluids: decay of autocorrelation for vapor-liquid transitions

We performed molecular dynamics simulations to study relaxation phenomena during vapor-liquid transitions in a single component Lennard-Jones system. Results from two different overall densities are presented: one in the neighborhood of the vapor branch of the coexistence curve and the other being close to the critical density. The nonequilibrium morphologies, growth mechanisms and growth laws in the two cases are vastly different. In the low density case growth occurs via diffusive coalescence of droplets in a disconnected morphology. On the other hand, the elongated structure in the higher density case grows via advective transport of particles inside the tube-like liquid domains. The objective in this work has been to identify how the decay of the order-parameter autocorrelation, an important quantity to understand aging dynamics, differs in the two cases. In the case of the disconnected morphology, we observe a very robust power-law decay, as a function of the ratio of the characteristic lengths at the observation time and at the age of the system, whereas the results for the percolating structure appear rather complex. To quantify the decay in the latter case, unlike the standard method followed in a previous study, here we have performed a finite-size scaling analysis. The outcome of this analysis shows the presence of a strong preasymptotic correction, while revealing that in this case also, albeit in the asymptotic limit, the decay follows a power-law. Even though the corresponding exponents in the two cases differ drastically, this study, combined with a few recent ones, suggests that power-law behavior of this correlation function is rather universal in coarsening dynamics.

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.

Numerical modelling of non-ionic microgels: an overview

Microgels are complex macromolecules. These colloid-sized polymer networks possess internal degrees of freedom and, depending on the polymer(s) they are made of, can acquire a responsiveness to variations of the environment (temperature, pH, salt concentration, etc.). Besides being valuable for many practical applications, microgels are also extremely important to tackle fundamental physics problems. As a result, these last years have seen a rapid development of protocols for the synthesis of microgels, and more and more research has been devoted to the investigation of their bulk properties. However, from a numerical standpoint the picture is more fragmented, as the inherently multi-scale nature of microgels, whose bulk behaviour crucially depends on the microscopic details, cannot be handled at a single level of coarse-graining. Here we present an overview of the methods and models that have been proposed to describe non-ionic microgels at different length-scales, from the atomistic to the single-particle level. We especially focus on monomer-resolved models, as these have the right level of details to capture the most important properties of microgels, responsiveness and softness. We suggest that these microscopic descriptions, if realistic enough, can be employed as starting points to develop the more coarse-grained representations required to investigate the behaviour of bulk suspensions.

The race to the bottom: approaching the ideal glass?

Key to resolving the scientific challenge of the glass transition is to understand the origin of the massive increase in viscosity of liquids cooled below their melting temperature (avoiding crystallisation). A number of competing and often mutually exclusive theoretical approaches have been advanced to describe this phenomenon. Some posit a bona fide thermodynamic phase to an 'ideal glass', an amorphous state with exceptionally low entropy. Other approaches are built around the concept of the glass transition as a primarily dynamic phenomenon. These fundamentally different interpretations give equally good descriptions of the data available, so it is hard to determine which-if any-is correct. Recently however this situation has begun to change. A consensus has emerged that one powerful means to resolve this longstanding question is to approach the putative thermodynamic transition sufficiently closely, and a number of techniques have emerged to meet this challenge. Here we review the results of some of these new techniques and discuss the implications for the existence-or otherwise-of the thermodynamic transition to an ideal glass.

Computational investigation of microgels: synthesis and effect of the microstructure on the deswelling behavior

We present computer simulations of a realistic model of microgels. Unlike the regular network frameworks usually assumed in the simulation literature, we model and simulate a realistic and efficient synthesis route, mimicking cross-linking of functionalized chains inside a cavity. This model is inspired, e.g., by microfluidic fabrication of microgels from macromolecular precursors and is different from standard polymerization routes. The assembly of the chains is mediated by a low fraction of interchain crosslinks. The microgels are polydisperse in size and shape but globally spherical objects. In order to deeply understand the microgel structure and eventually improve the synthesis protocol we characterize their conformational properties and deswelling kinetics, and compare them with the results found for microgels obtained via underlying regular (diamond-like) structures. For the same molecular weight, monomer concentration and effective degree of cross-linking, the specific microstructure of the microgel has no significant effect on the locus of the volume phase transition (VPT). However, it strongly affects the deswelling kinetics, as revealed by a consistent analysis of the domain growth during the microgel collapse. Though both the disordered and the regular networks exhibit a similar early growth of the domains, an acceleration is observed in the regular network at the late stage of the collapse. Similar trends are found for the dynamic correlations coupled to the domain growth. As a consequence, the fast late processes for the domain growth and the dynamic correlations in the regular network are compensated, and the dynamic correlations follow a power-law dependence on the growing length scale that is independent of the microgel microstructure.

Distributions of pore sizes and atomic densities in binary mixtures revealed by molecular dynamics simulations

We report on the results of a molecular dynamics simulation study of porous glassy media, formed in the process of isochoric rapid quenching from a high-temperature liquid state. The transition to a porous solid occurs due to the concurrent processes of phase separation and material solidification. The study is focused on topographies of the model porous structures and their dependence on temperature and average density. To quantify the pore-size distributions, we put forth a scaling relation that provides a satisfactory data collapse in systems with high porosity. We also find that the local density of the solid domains in the porous structures is broadly distributed, and, with increasing average density, a distinct peak in the local density distribution is displaced toward higher values.

Evolution of the pore size distribution in sheared binary glasses

Molecular dynamics simulations are carried out to investigate mechanical properties and porous structure of binary glasses subjected to steady shear. The model vitreous systems were prepared via thermal quench at constant volume to a temperature well below the glass transition. The quiescent samples are characterized by a relatively narrow pore size distribution whose mean size is larger at lower glass densities. We find that in the linear regime of deformation, the shear modulus is a strong function of porosity, and the individual pores become slightly stretched while their structural topology remains unaffected. By contrast, with further increasing strain, the shear stress saturates to a density-dependent plateau value, which is accompanied by pore coalescence and a gradual development of a broader pore size distribution with a discrete set of peaks at large length scales.

Arrested and temporarily arrested states in a protein-polymer mixture studied by USAXS and VSANS

We investigate the transition of the phase separation kinetics from a complete to an arrested liquid-liquid phase separation (LLPS) in mixtures of bovine γ-globulin with polyethylene glycol (PEG). The solutions feature LLPS with upper critical solution temperature phase behavior. At higher PEG concentrations or low temperatures, non-equilibrium, gel-like states are found. The kinetics is followed during off-critical quenches by ultra-small angle X-ray scattering (USAXS) and very-small angle neutron scattering (VSANS). For shallow quenches a kinetics consistent with classical spinodal decomposition is found, with the characteristic length (ξ) growing with time as ξ ∼ t^1/3. For deep quenches, ξ grows only very slowly with a growth exponent smaller than 0.05 during the observation time, indicating an arrested phase separation. For intermediate quench depths, a novel growth kinetics featuring a three-stage coarsening is observed, with an initial classical coarsening, a subsequent slowdown of the growth, and a later resumption of coarsening approaching again ξ ∼ t^1/3. Samples featuring the three-stage coarsening undergo a temporarily arrested state. We hypothesize that, while intermittent coarsening and collapse might contribute to the temporary nature of the arrested state, migration-coalescence of the minority liquid phase through the majority glassy phase may be the main mechanism underlying this kinetics, which is also consistent with earlier simulation results.

Structural arrest and dynamic localization in biocolloidal gels

Casein micelles interacting via an entropic intermediate-ranged depletion attraction exhibit a fluid-to-gel transition due to arrested spinodal decomposition. The bicontinuous networked structure of the gel freezes shortly after formation. We determine the timescales of structural arrest from the build-up of network rigidity after pre-shear rejuvenation, and find that the arrest time as well as the plateau elastic modulus of the gel diverge as a function of the volume fraction and interaction potential. Moreover, we show using scaling from naïve mode coupling theory that their mechanical properties are dictated by their microscopic dynamics rather than their heterogeneous large scale structure.

Ultra-long-range dynamic correlations in a microscopic model for aging gels

We use large-scale computer simulations to explore the nonequilibrium aging dynamics in a microscopic model for colloidal gels. We find that gelation resulting from a kinetically arrested phase separation is accompanied by "anomalous" particle dynamics revealed by superdiffusive particle motion and compressed exponential relaxation of time correlation functions. Spatiotemporal analysis of the dynamics reveals intermittent heterogeneities producing spatial correlations over extremely large length scales. Our study is a microscopically resolved model reproducing all features of the spontaneous aging dynamics observed experimentally in soft materials.

Cell division and death inhibit glassy behaviour of confluent tissues

We investigate the effects of cell division and apoptosis on collective dynamics in two-dimensional epithelial tissues. Our model includes three key ingredients observed across many epithelia, namely cell-cell adhesion, cell death and a cell division process that depends on the surrounding environment. We show a rich non-equilibrium phase diagram depending on the ratio of cell death to cell division and on the adhesion strength. For large apoptosis rates, cells die out and the tissue disintegrates. As the death rate decreases, however, we show, consecutively, the existence of a gas-like phase, a gel-like phase, and a dense confluent (tissue) phase. Most striking is the observation that the tissue is self-melting through its own internal activity, ruling out the existence of any glassy phase.

Effects of vertical confinement on gelation and sedimentation of colloids

We consider the sedimentation of a colloidal gel under confinement in the direction of gravity. The confinement allows us to compare directly experiments and computer simulations, for the same system size in the vertical direction. The confinement also leads to qualitatively different behaviour compared to bulk systems: in large systems gelation suppresses sedimentation, but for small systems sedimentation is enhanced relative to non-gelling suspensions, although the rate of sedimentation is reduced when the strength of the attraction between the colloids is strong. We map interaction parameters between a model experimental system (observed in real space) and computer simulations. Remarkably, we find that when simulating the system using Brownian dynamics in which hydrodynamic interactions between the particles are neglected, we find that sedimentation occurs on the same timescale as the experiments. An analysis of local structure in the simulations showed similar behaviour to gelation in the absence of gravity.

Kinetics of liquid-liquid phase separation in protein solutions exhibiting LCST phase behavior studied by time-resolved USAXS and VSANS

We study the kinetics of the liquid-liquid phase separation (LLPS) and its arrest in protein solutions exhibiting a lower critical solution temperature (LCST) phase behavior using the combination of ultra-small angle X-ray scattering (USAXS) and very-small angle neutron scattering (VSANS). We employ a previously established model system consisting of bovine serum albumin (BSA) solutions with YCl₃. We follow the phase transition from sub-second to 10⁴ s upon an off-critical temperature jump. After a temperature jump, the USAXS profiles exhibit a peak that grows in intensity and shifts to lower q values with time. Below 45 °C, the characteristic length scale (ξ) obtained from this scattering peak increases with time with a power of about 1/3 for different sample compositions. This is in good agreement with the theoretical prediction for the intermediate stage of spinodal decomposition where the growth is driven by interface tension. Above 45 °C, ξ follows initially the 1/3 power law growth, then undergoes a significant slowdown, and an arrested state is reached below the denaturation temperature of the protein. This growth kinetics may indicate that the final composition of the protein-rich phase is located close to the high density branch of the LLPS binodal when a kinetically arrested state is reached.

Exploring the dynamics of phase separation in colloid-polymer mixtures with long range attraction

We have studied the kinetics of phase separation and gel formation in a low-dispersity colloid - non-adsorbing polymer system with long range attraction using small-angle light scattering. This system exhibits two-phase and three-phase coexistence of gas, liquid and crystal phases when the strength of attraction is between 2 and 4kBT and gel phases when the strength of attraction is increased. For those samples that undergo macroscopic phase separation, whether to gas-crystal, gas-liquid or gas-liquid-crystal coexistence, we observe dynamic scaling of the structure factor and growth of a characteristic length scale that behaves as expected for phase separation in fluids. In samples that gel, the power law associated with the growth of the dominant length scale is not equal to 1/3, but appears to depend mainly on the strength of attraction, decreasing from 1/3 for samples near the coexistence region to 1/27 at 8kBT, over a wide range of colloid and polymer concentrations.

Designing disordered materials using DNA-coated colloids of bacteriophage fd and gold

DNA has emerged as an exciting binding agent for programmable colloidal self-assembly. Its popularity derives from its unique properties: it provides highly specific short-ranged interactions and at the same time it acts as a steric stabilizer against non-specific van der Waals and Coulomb interactions. Because complementary DNA strands are linked only via hydrogen bonds, DNA-mediated binding is thermally reversible: it provides an effective attraction that can be switched off by raising the temperature only by a few degrees. In this article we introduce a new binary system made of DNA-functionalized filamentous fd viruses of ∼880 nm length with an aspect ratio of ∼100, and 50 nm gold nanoparticles (gold NPs) coated with the complementary DNA strands. When quenching mixtures below the melt temperature Tm, at which the attraction is switched on, we observe aggregation. Conversely, above Tm the system melts into a homogenous particulate 'gas'. We present the aggregation behavior of three different gold NP to virus ratios and compare them to a gel made solely of gold NPs. In particular, we have investigated the aggregate structures as a function of cooling rate and determine how they evolve as function of time for given quench depths, employing fluorescence microscopy. Structural information was extracted in the form of an effective structure factor and chord length distributions. Rapid cooling rates lead to open aggregates, while slower controlled cooling rates closer to equilibrium DNA hybridization lead to more fine-stranded gels. Despite the different structures we find that for both cooling rates the quench into the two-phase region leads to initial spinodal decomposition, which becomes arrested. Surprisingly, although the fine-stranded gel is disordered, the overall structure and the corresponding length scale distributions in the system are remarkably reproducible. Such highly porous systems can be developed into new functional materials.

Depletion, melting and reentrant solidification in mixtures of soft and hard colloids

We present extensive experimental and theoretical investigations on the structure, phase behavior, dynamics and rheology of model soft-hard colloidal mixtures realized with large, multiarm star polymers as the soft component and smaller, compact stars as the hard one. The number and length of the arms in star polymers control their softness, whereas the size ratio, the overall density and the composition are additional parameters varied for the mixtures. A coarse-grained theoretical strategy is employed to predict the structure of the systems as well as their ergodicity properties on the basis of mode coupling theory, for comparison with rheological measurements on the samples. We discovered that dynamically arrested star-polymer solutions recover their ergodicity upon addition of colloidal additives. At the same time the system displays demixing instability, and the binodal of the latter meets the glass line in a way that leads, upon addition of a sufficient amount of colloidal particles, to an arrested phase separation and reentrant solidification. We present evidence for a subsequent solid-to-solid transition well within the region of arrested phase separation, attributed to a hard-sphere-mixture type of glass, due to osmotic shrinkage of the stars at high colloidal particle concentrations. We systematically investigated the interplay of star functionality and size ratio with glass melting and demixing, and rationalized our findings by the depletion of the big stars due to the smaller colloids. This new depletion potential in which, contrary to the classic colloid-polymer case, the hard component depletes the soft one, has unique and novel characteristics and allows the calculation of phase diagrams for such mixtures. This work covers a broad range of soft-hard colloidal mixture compositions in which the soft component exceeds the hard one in size and provides general guidelines for controlling the properties of such complex mixtures.

Monte Carlo simulation of kinetically slowed down phase separation

Supercooled colloidal or molecular systems at low densities are known to form liquid, crystalline or glassy drops, which may remain isolated for a long time before they aggregate. This paper analyses the properties of this large time window, and how it can be tackled by computer simulation. We use single-particle and virtual move Monte Carlo simulations of short-range attractive spheres which are undercooled to the temperature region, where the spinodal intersects the attractive glass line. We study two different systems and we report the following kinetic behavior. A low-density system is shown to exhibit universal linear growth regimes under single-particle Monte Carlo correlating the growth rate to the local structure. These regimes are suppressed under collective motion, where droplets aggregate into a single large disordered domain. It is shown that the aggregation can be avoided and linear regimes recovered, if long-range repulsion is added to the short-range attraction. The results provide an insight into the behavior of the virtual move algorithm generating cluster moves according to the local forcefields. We show that different choices of maximum Monte Carlo displacement affect the dynamical trajectories but lead to the same kinetically slowed down or arrested states.

Fragmentation and limits to dynamical scaling in viscous coarsening: an interrupted in situ x-ray tomographic study

X-ray microtomography was used to follow the coarsening of the structure of a ternary silicate glass experiencing phase separation in the liquid state. The volumes, surfaces, mean, and Gaussian curvatures of the domains of minority phase were measured after reconstruction of the 3D images and segmentation. A linear growth law of the characteristic length scale ℓ∼t was observed. A detailed morphological study was performed. While dynamical scaling holds for most of the geometrical observables under study, a progressive departure from scaling invariance of the distributions of local curvatures was evidenced. The latter results from a gradual fragmentation of the structure in the less viscous phase that also leads to a power-law size distribution of isolated domains.