Influence of the glass transition on the liquid-gas spinodal decomposition

Modular RNA motifs for orthogonal phase separated compartments

Recent discoveries in biology have highlighted the importance of protein and RNA-based condensates as an alternative to classical membrane-bound organelles. Here, we demonstrate the design of pure RNA condensates from nanostructured, star-shaped RNA motifs. We generate condensates using two different RNA nanostar architectures: multi-stranded nanostars whose binding interactions are programmed via linear overhangs, and single-stranded nanostars whose interactions are programmed via kissing loops. Through systematic sequence design, we demonstrate that both architectures can produce orthogonal (distinct and immiscible) condensates, which can be individually tracked via fluorogenic aptamers. We also show that aptamers make it possible to recruit peptides and proteins to the condensates with high specificity. Successful co-transcriptional formation of condensates from single-stranded nanostars suggests that they may be genetically encoded and produced in living cells. We provide a library of orthogonal RNA condensates that can be modularly customized and offer a route toward creating systems of functional artificial organelles for the task of compartmentalizing molecules and biochemical reactions.

Athermal quasistatic cavitation in amorphous solids: Effect of random pinning

Amorphous solids are known to fail catastrophically via fracture, and cavitation at nano-metric scales is known to play a significant role in such a failure process. Micro-alloying via inclusions is often used as a means to increase the fracture toughness of amorphous solids. Modeling such inclusions as randomly pinned particles that only move affinely and do not participate in plastic relaxations, we study how the pinning influences the process of cavitation-driven fracture in an amorphous solid. Using extensive numerical simulations and probing in the athermal quasistatic limit, we show that just by pinning a very small fraction of particles, the tensile strength is increased, and also the cavitation is delayed. Furthermore, the cavitation that is expected to be spatially heterogeneous becomes spatially homogeneous by forming a large number of small cavities instead of a dominant cavity. The observed behavior is rationalized in terms of screening of plastic activity via the pinning centers, characterized by a screening length extracted from the plastic-eigenmodes.

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.

Active Brownian particles in random and porous environments

The transport of active particles may occur in complex environments, in which it emerges from the interplay between the mobility of the active components and the quenched disorder of the environment. Here, we explore the structural and dynamical properties of active Brownian particles (ABPs) in random environments composed of fixed obstacles in three dimensions. We consider different arrangements of the obstacles. In particular, we consider two particular situations corresponding to experimentally realizable settings. First, we model pinning particles in (non-overlapping) random positions and, second, in a percolating gel structure and provide an extensive characterization of the structure and dynamics of ABPs in these complex environments. We find that the confinement increases the heterogeneity of the dynamics, with new populations of absorbed and localized particles appearing close to the obstacles. This heterogeneity has a profound impact on the motility induced phase separation exhibited by the particles at high activity, ranging from nucleation and growth in random disorder to a complex phase separation in porous environments.

Microscopic observation of two-level systems in a metallic glass model

The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements.

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.

Universal mechanical instabilities in the energy landscape of amorphous solids: Evidence from athermal quasistatic expansion

Using numerical simulations, we study the failure of an amorphous solid under athermal quasistatic expansion starting from a homogeneous high-density state. During the expansion process, plastic instabilities occur, manifested via sudden jumps in pressure and energy, with the largest event happening via cavitation leading to the material's yielding. We demonstrate that all these plastic events are characterized by saddle-node bifurcation, during which the smallest nonzero eigenvalue of the Hessian matrix vanishes via a square-root singularity. We find that after yielding and prior to complete fracture, the statistics of pressure or energy jumps corresponding to the plastic events show subextensive system-size scaling, similar to the case of simple shear but with different exponents. Thus, overall, our paper reveals universal features in the fundamental characteristics during mechanical failure in amorphous solids under any quasistatic deformation protocol.

Structure induced laminar vortices control anomalous dispersion in porous media

Natural porous systems, such as soil, membranes, and biological tissues comprise disordered structures characterized by dead-end pores connected to a network of percolating channels. The release and dispersion of particles, solutes, and microorganisms from such features is key for a broad range of environmental and medical applications including soil remediation, filtration and drug delivery. Yet, owing to the stagnant and opaque nature of these disordered systems, the role of microscopic structure and flow on the dispersion of particles and solutes remains poorly understood. Here, we use a microfluidic model system that features a pore structure characterized by distributed dead-ends to determine how particles are transported, retained and dispersed. We observe strong tailing of arrival time distributions at the outlet of the medium characterized by power-law decay with an exponent of 2/3. Using numerical simulations and an analytical model, we link this behavior to particles initially located within dead-end pores, and explain the tailing exponent with a hopping across and rolling along the streamlines of vortices within dead-end pores. We quantify such anomalous dispersal by a stochastic model that predicts the full evolution of arrival times. Our results demonstrate how microscopic flow structures can impact macroscopic particle transport.

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.

Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions

We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.

Electric Fields Enhance Ice Formation from Water Vapor by Decreasing the Nucleation Energy Barrier

Video images of ice formation from moist air under temperature and electric potential gradients reveal that ambient electricity enhances ice production rates while changing the habit of ice particles formed under low supersaturation. The crystals formed under an electric field are needles and dendrites instead of the isometric ice particles obtained within a Faraday cage. Both a non-classical mechanism and classical nucleation theory independently explain the observed mutual feedback between ice formation and its electrification. The elongated shapes result from electrostatic repulsion at the crystal surfaces, opposing the attractive intermolecular forces and thus lowering the ice-air interfacial tension. The video images allow for the estimation of ice particle dimensions, weight, and speed within the electric field. Feeding this data on standard equations from electrostatics shows that the ice surface charge density attains 0.62–1.25 × 10−6 C·m−2, corresponding to 73–147 kV·m−1 potential gradients, reaching the range measured within thunderstorms. The present findings contribute to a better understanding of natural and industrial processes involving water phase change by acknowledging the presence and effects of the pervasive electric fields in the ambient environment.

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.

Nonequilibrium Processes in Polymer Membrane Formation: Theory and Experiment

Porous polymer and copolymer membranes are useful for ultrafiltration of functional macromolecules, colloids, and water purification. In particular, block copolymer membranes offer a bottom-up approach to form isoporous membranes. To optimize permeability, selectivity, longevity, and cost, and to rationally design fabrication processes, direct insights into the spatiotemporal structure evolution are necessary. Because of a multitude of nonequilibrium processes in polymer membrane formation, theoretical predictions via continuum models and particle simulations remain a challenge. We compiled experimental observations and theoretical approaches for homo- and block copolymer membranes prepared by nonsolvent-induced phase separation and highlight the interplay of multiple nonequilibrium processes─evaporation, solvent-nonsolvent exchange, diffusion, hydrodynamic flow, viscoelasticity, macro- and microphase separation, and dynamic arrest─that dictates the complex structure of the membrane on different scales.

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.

Effect of softness on glass melting and re-entrant solidification in mixtures of soft and hard colloids

We present a systematic investigation of the structure and dynamic properties of model soft-hard colloidal mixtures. Results of a coarse-grained theoretical model are contrasted with rheological data, where the soft and hard colloids are mimicked by large star polymers with high functionality as the soft component and smaller stars with ultrahigh functionality as the hard one. Previous work by us revealed the recovery of the ergodicity of glassy soft star solutions and subsequent arrested phase separation and re-entrant solid transition upon progressive addition of small hard depletants. Here, we use different components to show that a small variation in softness has a significant impact on the state diagram of such mixtures. In particular, we establish that rendering the soft component more penetrable and modifying the size ratio bring about a remarkable shift in both the phase separation region and the glass-melting line so that the region of restored ergodicity can be notably enhanced and extended to much higher star polymer concentrations than for pure systems. We further rationalize our findings by analyzing the features of the depletion interaction induced by the smaller component that result from the interplay between the size ratio and the softness of the large component. These results demonstrate the great sensitivity of the phase behavior of entropic mixtures to small changes in the molecular architecture of the soft stars and point to the importance of accounting for details of the internal microstructure of soft colloidal particles for tailoring the flow properties of soft composites.

Dynamic Crossover in Fluids: From Hard Spheres to Molecules

We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.

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.

Glassy Behavior of Sticky Spheres: What Lies beyond Experimental Timescales?

We use the swap Monte Carlo algorithm to analyze the glassy behavior of sticky spheres in equilibrium conditions at densities where conventional simulations and experiments fail to reach equilibrium, beyond predicted phase transitions and dynamic singularities. We demonstrate the existence of a unique ergodic region comprising all the distinct phases previously reported, except for a phase-separated region at strong adhesion. All structural and dynamic observables evolve gradually within this ergodic region, the physics evolving smoothly from well-known hard sphere glassy behavior at small adhesions and large densities, to a more complex glassy regime characterized by unusually broad distributions of relaxation timescales and length scales at large adhesions.

Violation of Stokes-Einstein and Stokes-Einstein-Debye relations in polymers at the gas-supercooled liquid coexistence

Molecular dynamics simulations are performed on a system of model linear polymers to look at the violations of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations near the mode coupling theory transition temperatureT_cat three (one higher and two lower) densities. At low temperatures, both lower density systems show stable gas-supercooled-liquid coexistence whereas the higher density system is homogeneous. We show that monomer density relaxation exhibits SE violation for all three densities, whereas molecular density relaxation shows a weak violation of the SE relation nearT_cin both lower density systems. This study identifies disparity in monomer mobility and observation of jumplike motion in the typical monomer trajectories resulting in the SE violations. In addition to the SE violation, a weak SED violation is observed in the gas-supercooled-liquid coexisting domains of the lower densities. Both lower density systems also show a decoupling of translational and rotational dynamics in this polymer system.

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.

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.

Multi-component colloidal gels: interplay between structure and mechanical properties

We present a detailed numerical study of multi-component colloidal gels interacting sterically and obtained by arrested phase separation. Under deformation, we found that the interplay between the different intertwined networks is key. Increasing the number of components leads to softer solids that can accommodate progressively larger strains before yielding. The simulations highlight how this is the direct consequence of the purely repulsive interactions between the different components, which end up enhancing the linear response of the material. Our work provides new insight into mechanisms at play for controlling the material properties and opens a road to new design principles for soft composite solids.

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.

Formation of Ultrastable Glasses via Precipitation: A Modeling Study

The precipitation of a glass forming solute from solution is modeled using a lattice model previously introduced to study dissolution kinetics of amorphous materials. The model includes the enhancement of kinetics at the surface of a glass in contact with a plasticizing solvent. We demonstrate that precipitation can produce a glass substantially more stable than that produced by very long time annealing of the bulk glass former. The energy of these ultrastable amorphous precipitates is found to be dominated by residual solvent rather than high energy glass configurations.

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.

Formation of porous crystals via viscoelastic phase separation

Viscoelastic phase separation of colloidal suspensions can be interrupted to form gels either by glass transition or by crystallization. With a new confocal microscopy protocol, we follow the entire kinetics of phase separation, from homogeneous phase to different arrested states. For the first time in experiments, our results unveil a novel crystallization pathway to sponge-like porous crystal structures. In the early stages, we show that nucleation requires a structural reorganization of the liquid phase, called stress-driven ageing. Once nucleation starts, we observe that crystallization follows three different routes: direct crystallization of the liquid phase, the Bergeron process, and Ostwald ripening. Nucleation starts inside the reorganized network, but crystals grow past it by direct condensation of the gas phase on their surface, driving liquid evaporation, and producing a network structure different from the original phase separation pattern. We argue that similar crystal-gel states can be formed in monatomic and molecular systems if the liquid phase is slow enough to induce viscoelastic phase separation, but fast enough to prevent immediate vitrification. This provides a novel pathway to form nanoporous crystals of metals and semiconductors without dealloying, which may be important for catalytic, optical, sensing, and filtration applications.

Optically transparent dense colloidal gels

Traditionally it has been difficult to study the porous structure of dense colloidal gels and (macro) molecular transport through them simply because of the difference in refractive index between the colloid material and the continuous fluid phase surrounding it, rendering the samples opaque even at low colloidal volume fractions. Here, we demonstrate a novel colloidal gel that can be refractive index-matched in aqueous solutions owing to the low refractive index of fluorinated latex (FL)-particles (n = 1.37). Synthesizing them from heptafluorobutyl methacrylate using emulsion polymerization, we demonstrate that they can be functionalized with short DNA sequences via a dense brush-layer of polystyrene-b-poly(ethylene oxide) block-copolymers (PS-PEO). The block-copolymer, holding an azide group at the free PEO end, was grafted to the latex particle utilizing a swelling-deswelling method. Subsequently, DNA was covalently attached to the azide-end of the block copolymer via a strain-promoted alkyne-azide click reaction. For comparison, we present a structural study of single gels made of FL-particles only and composite gels made of a percolating FL-colloid gel coated with polystyrene (PS) colloids. Further we demonstrate that the diffusivity of tracer colloids dispersed deep inside a refractive index matched FL-colloidal gel can be measured as function of the local confinement using Dynamic Differential Microscopy (DDM).

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.

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.

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.

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.

Microdynamics and arrest of coarsening during spinodal decomposition in thermoreversible colloidal gels

Coarsening and kinetic arrest of colloidal systems undergoing spinodal decomposition (SD) is a conserved motif for forming hierarchical, bicontinuous structures. Although the thermodynamic origins of SD in colloids are widely known, the microstructural processes responsible for its coarsening and associated dynamics en route to arrest remain elusive. To better elucidate the underlying large-scale microdynamical processes, we study a colloidal system with moderate-range attractions which displays characteristic features of arrested SD, and study its dynamics during coarsening through a combination of differential dynamic microscopy and real-space tracking. Using these recently developed imaging techniques, we reveal directly that the coarsening arises from collective dynamics of dense domains, which undergo slow, intermittent, and ballistic motion. These collective motions indicate interfacial effects to be the driving force of coarsening. The nature of the gelation enables control of the arrested length scale of coarsening by the depths of quenching into the spinodal regime, which we demonstrate to provide an effective means to control the elasticity of colloidal gels.

Probing Colloidal Gels at Multiple Length Scales: The Role of Hydrodynamics

Colloidal gels are out-of-equilibrium structures, made up of a rarefied network of colloidal particles. Comparing experiments to numerical simulations, with hydrodynamic interactions switched off, we demonstrate the crucial role of the solvent for gelation. Hydrodynamic interactions suppress the formation of larger local equilibrium structures of closed geometry, and instead lead to the formation of highly anisotropic threads, which promote an open gel network. We confirm these results with simulations which include hydrodynamics. Based on three-point correlations, we propose a scale-resolved quantitative measure for the anisotropy of the gel structure. We find a strong discrepancy for interparticle distances just under twice the particle diameter between systems with and without hydrodynamics, quantifying the role of hydrodynamics from a structural point of view.

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.

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

We use large-scale molecular dynamics simulations of a simple glass-forming system to investigate how its liquid-gas phase separation kinetics depends on temperature. A shallow quench leads to a fully demixed liquid-gas system whereas a deep quench makes the dense phase undergo a glass transition and become an amorphous solid. This glass has a gel-like bicontinuous structure that evolves very slowly with time and becomes fully arrested in the limit where thermal fluctuations become negligible. We show that the phase separation kinetics changes qualitatively with temperature, the microscopic dynamics evolving from a surface tension-driven diffusive motion at high temperature to a strongly intermittent, heterogeneous, and thermally activated dynamics at low temperature, with a logarithmically slow growth of the typical domain size. These results elucidate the microscopic mechanisms underlying a specific class of viscoelastic phase separation.

Morphological evolution of domains in spinodal decomposition

Domain growth in spinodal decomposition is usually described by a single time-evolving length scale. We show that the evolution of morphology of domains is nonmonotonic. The domains elongate rapidly at first and then, with the help of hydrodynamics, return to a more circular shape. The initial elongation phase does not alter with hydrodynamics. A small deviation from critical composition changes the morphology dramatically.

Self-assembly and cooperative dynamics of a model colloidal gel network

We study the assembly into a gel network of colloidal particles, via effective interactions that yield local rigidity and make dilute network structures mechanically stable. The self-assembly process can be described by a Flory-Huggins theory, until a network of chains forms, whose mesh size is on the order of, or smaller than, the persistence length of the chains. The localization of the particles in the network, akin to some extent to caging in dense glasses, is determined by the network topology, and the network restructuring, which takes place via bond breaking and recombination, is characterized by highly cooperative dynamics. We use NVE and NVT molecular dynamics as well as Langevin dynamics and find a qualitatively similar time dependence of time correlations and of the dynamical susceptibility of the restructuring gel. This confirms that the cooperative dynamics emerge from the mesoscale organization of the network.

Aggregation dynamics, structure, and mechanical properties of bigels

Recently we have introduced bigels, inter-penetrating gels made of two different colloidal species. Even if particles with simple short-range isotropic potential are employed, the selective interactions enable the tunability of the self-assembly, leading to the formation of complex structures. In the present paper, we explore the non-equilibrium dynamics and the phenomenology underlying the kinetic arrest under quench and the formation of bigels. We demonstrate that the peculiar bigel kinetics can be described through an arrested spinodal decomposition driven by demixing of the colloidal species. The role played by the presence of a second colloidal species on the phase diagram, as expanded to account for the increased number of parameters, is clarified both via extensive numerical simulations and experiments. We provide details on the realisation of bigels, by means of DNA-coated colloids (DNACCs), and the consequent imaging techniques. Moreover we evidence, by comparison with the usual one-component gel formation, the emergence of controllable timescales in the aggregation of the bigels, whose final stages are also experimentally studied to provide morphological details. Finally, we use numerical models to simulate the bigel response to mechanical strain, highlighting how such a new material can bear significantly higher stress compared to the usual one-component gel. We conclude by discussing possible technological uses and by providing insights on the viable research steps to undertake for more complex and yet tuneable multi-component colloidal systems.