Higher-order glass-transition singularities in colloidal systems with attractive interactions

Microscopic theory of the elastic shear modulus and length-scale-dependent dynamic re-entrancy phenomena in very dense sticky particle fluids

We apply the hybrid projectionless dynamic theory (hybrid PDT) formulation of the elastically collective nonlinear Langevin equation (ECNLE) activated dynamics approach to study dense fluids of sticky spheres interacting with short range attractions. Of special interest is the problem of non-monotonic evolution with short range attraction strength of the elastic modulus ("re-entrancy") at very high packing fractions far beyond the ideal mode coupling theory (MCT) nonergodicity boundary. The dynamic force constraints explicitly treat the bare attractive forces that drive transient physical bond formation, while a projection approximation is employed for the singular hard-sphere potential. The resultant interference between repulsive and attractive forces contribution to the dynamic vertex results in the prediction of localization length and elastic modulus re-entrancy, qualitatively consistent with experiments. The non-monotonic evolution of the structural (alpha) relaxation time predicted by the ECNLE theory with the hybrid PDT approach is explored in depth as a function of packing fraction, attraction strength, and attraction range. Isochronal dynamic arrest boundaries based on activated relaxation display the classic non-monotonic glass melting form. Comparisons of these results with the corresponding predictions of ideal MCT, and also the ECNLE and NLE activated theories based on projection, reveal large qualitative differences. The consequences of stochastic trajectory fluctuations on intra-cage single particle dynamics with variable strength of attractions are also studied. Large dynamical heterogeneity effects in attractive glasses are properly captured. These include a rapidly increasing amplitude of the non-Gaussian parameter with packing fraction and a non-monotonic evolution with attraction strength, in qualitative accord with recent simulations. Extension of the microscopic theoretical approach to treat double yielding in attractive glass nonlinear rheology is possible.

The distinguishable-particle lattice model of glasses in three dimensions

The nature of glassy states in realistic finite dimensions is still under fierce debate. Lattice models can offer valuable insights and facilitate deeper theoretical understanding. Recently, a disordered-interacting lattice model with distinguishable particles in two dimensions (2D) has been shown to produce a wide range of dynamical properties of structural glasses, including the slow and heterogeneous characteristics of the glassy dynamics, various fragility behaviors of glasses, and so on. These findings support the usefulness of this model for modeling structural glasses. An important question is whether such properties still hold in the more realistic three dimensions. In this study, we aim to extend the distinguishable-particle lattice model (DPLM) to three dimensions (3D) and explore the corresponding glassy dynamics. Through extensive kinetic Monte Carlo simulations, we found that the 3D DPLM exhibits many typical glassy behaviors, such as plateaus in the mean square displacement of particles and the self-intermediate scattering function, dynamic heterogeneity, variability of glass fragilities, and so on, validating the effectiveness of the DPLM in a broader realistic setting. The observed glassy behaviors of the 3D DPLM appear similar to those of its 2D counterpart, in accordance with recent findings in molecular models of glasses. We further investigate the role of void-induced motions in dynamical relaxations and discuss their relation to dynamic facilitation. As lattice models tend to keep the minimal but important modeling elements, they are typically much more amenable to analysis. Therefore, we envisage that the DPLM will benefit future theoretical developments, such as the configuration tree theory, towards a more comprehensive understanding of structural glasses.

Glassy phases of the Gaussian core model

We present results from molecular dynamics simulations exploring the supercooled dynamics of the Gaussian Core Model in the low- and intermediate-density regimes. In particular, we analyse the transition from the low-density hard-sphere-like glassy dynamics to the high-density one. The dynamics at low densities is well described by the caging mechanism, giving rise to intermittent dynamics. At high densities, the particles undergo a more continuous motion in which the concept of cage loses its meaning. We elaborate on the idea that these different supercooled dynamics are in fact the precursors of two different glass states.

Comparative study of the dynamics of colloidal glass and gel

We investigate and compare the difference in the dynamics of two arrested states: colloidal glass and colloidal gel. Real-space experiments reveal two distinct nonergodicity origins for their slow dynamics, namely, cage effects for the glass and attractive bondings for the gel. Such distinct origins lead to a faster decay of the correlation function and a smaller nonergodicity parameter of the glass than those of the gel. We also find that the gel exhibits stronger dynamical heterogeneity compared with the glass due to the greater correlated motions in the gel. Moreover, a logarithmic decay in the correlation function is observed as the two nonergodicity origins merge, consistent with the mode coupling theory.

Complex optical transport, dynamics, and rheology of intermediately attractive emulsions

Introducing short-range attractions in Brownian systems of monodisperse colloidal spheres can substantially impact their structures and consequently their optical transport and rheological properties. Here, for size-fractionated colloidal emulsions, we show that imposing an intermediate strength of attraction, well above but not much larger than thermal energy ([Formula: see text] [Formula: see text], through micellar depletion leads to a striking notch in the measured inverse mean free path of optical transport, [Formula: see text], as a function of droplet volume fraction, [Formula: see text]. This notch, which appears between the hard-sphere glass transition, [Formula: see text], and maximal random jamming, [Formula: see text], implies the existence of a greater population of compact dense clusters of droplets, as compared to tenuous networks of droplets in strongly attractive emulsion gels. We extend a prior decorated core-shell network model for strongly attractive colloidal systems to include dense non-percolating clusters that do not contribute to shear rigidity. By constraining this extended model using the measured [Formula: see text], we improve and expand the microrheological interpretation of diffusing wave spectroscopy (DWS) experiments made on attractive colloidal systems. Our measurements and modeling demonstrate richness and complexity in optical transport and shear rheological properties of dense, disordered colloidal systems having short-range intermediate attractions between moderately attractive glasses and strongly attractive gels.

The Hydrophobic Effect Studied by Using Interacting Colloidal Suspensions

Interactions between nanoparticles (NPs) determine their self-organization and dynamic processes. In these systems, a quantitative description of the interparticle forces is complicated by the presence of the hydrophobic effect (HE), treatable only qualitatively, and due to the competition between the hydrophobic and hydrophilic forces. Recently, instead, a sort of crossover of HE from hydrophilic to hydrophobic has been experimentally observed on a local scale, by increasing the temperature, in pure confined water and studying the occurrence of this crossover in different water-methanol solutions. Starting from these results, we then considered the idea of studying this process in different nanoparticle solutions. By using photon correlation spectroscopy (PCS) experiments on dendrimer with OH terminal groups (dissolved in water and methanol, respectively), we show the existence of this hydrophobic-hydrophilic crossover with a well defined temperature and nanoparticle volume fraction dependence. In this frame, we have used the mode coupling theory extended model to evaluate the measured time-dependent density correlation functions (ISFs). In this context we will, therefore, show how the measured spectra are strongly dependent on the specificity of the interactions between the particles in solution. The observed transition demonstrates that just the HE, depending sensitively on the system thermodynamics, determines the hydrophobic and hydrophilic interaction properties of the studied nanostructures surface.

Brownian dynamics simulations of hard rods in external fields and with contact interactions

We propose a simulation method for Brownian dynamics of hard rods in one dimension for arbitrary continuous external force fields. It is an event-driven procedure based on the fragmentation and mergers of clusters formed by particles in contact. It allows one to treat particle interactions in addition to the hard-sphere exclusion as long as the corresponding interaction forces are continuous functions of the particle coordinates. We furthermore develop a treatment of sticky hard spheres as described by Baxter's contact interaction potential.

The Interplay between the Theories of Mode Coupling and of Percolation Transition in Attractive Colloidal Systems

In the recent years a considerable effort has been devoted to foster the understanding of the basic mechanisms underlying the dynamical arrest that is involved in glass forming in supercooled liquids and in the sol-gel transition. The elucidation of the nature of such processes represents one of the most challenging unsolved problems in the field of material science. In this context, two important theories have contributed significantly to the interpretation of these phenomena: the Mode-Coupling theory (MCT) and the Percolation theory (PT). These theories are rooted on the two pillars of statistical physics, universality and scale laws, and their original formulations have been subsequently modified to account for the fundamental concepts of Energy Landscape (EL) and of the universality of the fragile to strong dynamical crossover (FSC). In this review, we discuss experimental and theoretical results, including Molecular Dynamics (MD) simulations, reported in the literature for colloidal and polymer systems displaying both glass and sol-gel transitions. Special focus is dedicated to the analysis of the interferences between these transitions and on the possible interplay between MCT and PT. By reviewing recent theoretical developments, we show that such interplay between sol-gel and glass transitions may be interpreted in terms of the extended F13 MCT model that describes these processes based on the presence of a glass-glass transition line terminating in an A3 cusp-like singularity (near which the logarithmic decay of the density correlator is observed). This transition line originates from the presence of two different amorphous structures, one generated by the inter-particle attraction and the other by the pure repulsion characteristic of hard spheres. We show here, combining literature results with some new results, that such a situation can be generated, and therefore experimentally studied, by considering colloidal-like particles interacting via a hard core plus an attractive square well potential. In the final part of this review, scaling laws associated both to MCT and PT are applied to describe, by means of these two theories, the specific viscoelastic properties of some systems.

Motile Topological Defects Hinder Dynamical Arrest in Dense Liquids of Active Ellipsoids

Recent numerical studies have identified the persistence time of active motion as a critical parameter governing glassy dynamics in dense active matter. Here we studied dynamics in liquids of granular active ellipsoids with tunable persistence and velocity. We show that increasing the persistence time at moderate supercooling is equivalent to increasing the strength of attraction in equilibrium liquids and results in reentrant dynamics not just in the translational degrees of freedom, as anticipated, but also in the orientational ones. However, at high densities, motile topological defects, unique to active liquids of elongated particles, hindered dynamical arrest. Most remarkably, for the highest activity, we observed intermittent dynamics due to the jamming-unjamming of these defects for the first time.

Toward a flow-dependent phase-stability criterion: Osmotic pressure in sticky flowing suspensions

Equilibrium phase instability of colloids is robustly predicted by the Vliegenthart-Lekkerkerker (VL) critical value of the second virial efficient, but no such general criterion has been established for suspensions undergoing flow. A transition from positive to negative osmotic pressure is one mechanical hallmark of a change in phase stability in suspensions and provides a natural extension of the equilibrium osmotic pressure encoded in the second virial coefficient. Here, we propose to study the non-Newtonian rheology of an attractive colloidal suspension using the active microrheology framework as a model for focusing on the pair trajectories that underlie flow stability. We formulate and solve a Smoluchowski relation to understand the interplay between attractions, hydrodynamics, Brownian motion, and flow on particle microstructure in a semi-dilute suspension and utilize the results to study the viscosity and particle-phase osmotic pressure. We find that an interplay between attractions and hydrodynamics leads to dramatic changes in the nonequilibrium microstructure, which produces a two-stage flow-thinning of viscosity and leads to pronounced flow-induced negative osmotic pressure. We summarize these findings with an osmotic pressure heat map that predicts where hydrodynamic enhancement of attractive bonds encourages flow-induced aggregation or phase separation. We identify a critical isobar-a flow-induced critical pressure consistent with phase instability and a nonequilibrium extension of the VL criterion.

Glassy dynamics of sticky hard spheres beyond the mode-coupling regime

Sticky hard spheres, i.e., hard particles decorated with a short-ranged attractive interaction potential, constitute a relatively simple model with highly non-trivial glassy dynamics. The mode-coupling theory of the glass transition (MCT) offers a qualitative account of the complex reentrant dynamics of sticky hard spheres, but the predicted glass transition point is notoriously underestimated. Here we apply an improved first-principles-based theory, referred to as generalized mode-coupling theory (GMCT), to sticky hard spheres. This theoretical framework seeks to go beyond MCT by hierarchically expanding the dynamics in higher-order density correlation functions. We predict the phase diagrams from the first few levels of the GMCT hierarchy and the dynamics-related critical exponents, all of which are much closer to the empirical observations than MCT. Notably, the prominent reentrant glassy dynamics, the glass-glass transition, and the higher-order bifurcation singularity classes (A3 and A4) of sticky hard spheres are found to be preserved within GMCT at arbitrary order. Moreover, we demonstrate that when the hierarchical order of GMCT increases, the effect of the short-ranged attractive interactions becomes more evident in the dynamics. This implies that GMCT is more sensitive to subtle microstructural differences than MCT, and that the framework provides a promising first-principles approach to systematically go beyond the MCT regime.

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.

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.

Self-motion and heterogeneous droplet dynamics in moderately attractive dense emulsions

We show that diffusing wave spectroscopy (DWS) is sensitive to the presence of a moderate short-range attraction between droplets in uniform fractionated colloidal emulsions near and below the jamming point associated with monodisperse hard spheres. This moderate interdroplet attraction, induced by micellar depletion, has an energy of about ∼2.4k_BT, only somewhat larger than thermal energy. Although changes in the mean free path of optical transport caused by this moderate depletion attraction are small, DWS clearly reveals an additional secondary decay-to-plateau in the intensity autocorrelation function at long times that is not present when droplet interactions are nearly hard. We hypothesize that this secondary decay-to-plateau does not reflect the average self-motion of individual droplets experiencing Brownian excitations, but instead results from heterogeneous dynamics involving a sub-population of droplets that still experience bound motion yet with significantly larger displacements than the average. By effectively removing the contribution of this secondary decay-to-plateau, which is linked to greater local heterogeneity in droplet structure caused by the moderate attraction, we obtain self-motion mean square displacements (MSDs) of droplets that reflect only the initial primary decay-to-plateau. Moreover, we show that droplet self-motion primary plateau MSDs can be interpreted using the generalized Stokes-Einstein relation of passive microrheology, yielding quantitative agreement with plateau elastic shear moduli measured mechanically.

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.

Generalized mode-coupling theory of the glass transition. II. Analytical scaling laws

Generalized mode-coupling theory (GMCT) constitutes a systematically correctable, first-principles theory to study the dynamics of supercooled liquids and the glass transition. It is a hierarchical framework that, through the incorporation of increasingly many particle density correlations, can remedy some of the inherent limitations of the ideal mode-coupling theory (MCT). However, despite MCT's limitations, the ideal theory also enjoys several remarkable successes, notably including the analytical scaling laws for the α- and β-relaxation dynamics. Here, we mathematically derive similar scaling laws for arbitrary-order multi-point density correlation functions obtained from GMCT under arbitrary mean-field closure levels. More specifically, we analytically derive the asymptotic and preasymptotic solutions for the long-time limits of multi-point density correlators, the critical dynamics with two power-law decays, the factorization scaling laws in the β-relaxation regime, and the time-density superposition principle in the α-relaxation regime. The two characteristic power-law-divergent relaxation times for the two-step decay and the non-trivial relation between their exponents are also obtained. The validity ranges of the leading-order scaling laws are also provided by considering the leading preasymptotic corrections. Furthermore, we test these solutions for the Percus-Yevick hard-sphere system. We demonstrate that GMCT preserves all the celebrated scaling laws of MCT while quantitatively improving the exponents, rendering the theory a promising candidate for an ultimately quantitative first-principles theory of glassy dynamics.

Diffusing wave microrheology of strongly attractive dense emulsions

We advance the microrheological interpretation of optical diffusing wave spectroscopy (DWS) measurements of strongly attractive emulsions at dense droplet volume fractions, ϕ. Beyond accounting for collective scattering, we show that measuring the mean free path of optical transport over a wide range of ϕ is necessary to quantify the effective size of the DWS probes, which we infer to be local dense clusters of droplets through a decorated core-shell network model. This approach yields microrheological elastic shear moduli that are in quantitative agreement with mechanical rheometry.

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules (2019 J. Chem. Phys. 150 024903). In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.

Complex modulus and compliance for airway smooth muscle cells

A cell can be described as a complex viscoelastic material with structural relaxations that is modulated by thermal and chemically nonequilibrium processes. Tissue morphology and function rely upon cells' physical responses to mechanical force. We measured the frequency-dependent mechanical relaxation response of adherent human airway smooth muscle cells under adenosine triphosphate (ATP) depletion and normal ATP conditions. The frequency dependence of the complex compliance J^{*} and modulus G^{*} was measured over the frequencies 10^{-1}<f<10^{3} Hz at selected temperatures between 4<T<54^{∘}C. Our results show characteristic relaxation features which can be interpreted by the mode-coupling theory (MCT) of viscoelastic liquids. We analyze the shape of the spectra in terms of a so-called A_{4} scenario with logarithmic scaling laws. Characteristic timescales τ_{β} and τ_{α} appear with corresponding energy barriers E_{β}≈(10-20)k_{B}T and E_{α}≈(20-30)k_{B}T. We demonstrate that cells are close to a glass transition. We find that the cell becomes softer around physiological temperatures, where its surface structure is more liquid-like with a plateau modulus around 0.1-0.8 kPa compared with the more solid-like interior cytoskeletal structures with a plateau modulus 1-15 kPa. Corresponding values for the viscosity are 10^{2}-10^{3} Pa s for the surface structures closer to the membrane and 10^{4}-10^{6} Pa s for the core cytoskeletal structures.

Temperature expansions in the square-shoulder fluid. II. Thermodynamics

In Paper I [O. Coquand and M. Sperl, J. Chem. Phys. 152, 124112 (2020)], we derived analytical expressions for the structure factor of the square-shoulder potential in a perturbative way around the high- and low-temperature regimes. Here, various physical properties of these solutions are derived. In particular, we investigate the large wave number sector and relate it to the contact values of the pair-correlation function. Then, the thermoelastic properties of the square-shoulder fluids are discussed.

Active topological glass

The glass transition in soft matter systems is generally triggered by an increase in packing fraction or a decrease in temperature. It has been conjectured that the internal topology of the constituent particles, such as polymers, can cause glassiness too. However, the conjecture relies on immobilizing a fraction of the particles and is therefore difficult to fulfill experimentally. Here we show that in dense solutions of circular polymers containing (active) segments of increased mobility, the interplay of the activity and the topology of the polymers generates an unprecedented glassy state of matter. The active isotropic driving enhances mutual ring threading to the extent that the rings can relax only in a cooperative way, which dramatically increases relaxation times. Moreover, the observed phenomena feature similarities with the conformation and dynamics of the DNA fibre in living nuclei of higher eukaryotes.

Slowing down supercooled liquids by manipulating their local structure

Glasses remain an elusive and poorly understood state of matter. It is not clear how we can control the macroscopic dynamics of glassy systems by tuning the properties of their microscopic building blocks. In this paper, we propose a simple directional colloidal model that reinforces the optimal icosahedral local structure of binary hard-sphere glasses. We show that this specific symmetry results in a dramatic slowing down of the dynamics. Our results open the door to controlling the dynamics of dense glassy systems by selectively promoting specific local structural environments.

Glass-transition asymptotics in two theories of glassy dynamics: Self-consistent generalized Langevin equation and mode-coupling theory

We contrast the generic features of structural relaxation close to the idealized glass transition that are predicted by the self-consistent generalized Langevin equation theory (SCGLE) against those that are predicted by the mode-coupling theory of the glass transition (MCT). We present an asymptotic solution close to conditions of kinetic arrest that is valid for both theories, despite the different starting points that are adopted in deriving them. This in particular provides the same level of understanding of the asymptotic dynamics in the SCGLE as was previously done only for MCT. We discuss similarities and different predictions of the two theories for kinetic arrest in standard glass-forming models, as exemplified through the hard-sphere system. Qualitative differences are found for models where a decoupling of relaxation modes is predicted, such as the generalized Gaussian core model, or binary hard-sphere mixtures of particles with very disparate sizes. These differences, which arise in the distinct treatment of the memory kernels associated to self- and collective motion of particles, lead to distinct scenarios that are predicted by each theory for partially arrested states and in the vicinity of higher-order glass-transition singularities.

Dynamically asymmetric binary glass formers studied by dielectric and NMR spectroscopy

We investigate the component dynamics in asymmetric binary glass formers. Focusing on the dielectric spectra of the high-T_g components m-tricresyl phosphate and quinaldine mixed with toluene as low-T_g component, the broadend spectra cannot be described by Kohlrausch or Cole-Davidson (CD) functions. Instead, we apply a generalized CD function which allows to control the width of the susceptibility independently of its high-frequency flank. The spectra show a common broadening and failure of the frequency-temperature superposition with increasing toluene concentration. This is confirmed by stimulated echo experiments showing an increased stretching of the probed orientational correlation function. In analogy to the definition of T_g, we consider "isodynamic points". For each component, a different but linear concentration dependence of 1/T_iso is revealed, indicating different time scales. Qualitativly, we do not find significant differences for the present mixtures with T_g-contrasts of 63-89K compared to those with larger T_g-contrast ( [Formula: see text] K): Whereas the high-T_g component shows relaxation features similar to those of neat glass formers, yet, with "atypical" weak relaxation broadening, the faster low-T_g component displays pronounced dynamic heterogeneities. This is supported by scrutinizing NMR relaxation data of several mixtures investigated previously as a function of concentration. A universal evolution of the dynamics of the high-T_g as well as the low-T_g component is suggested for mixtures with high [Formula: see text]T_g .

α-relaxation, shear viscosity, and elastic moduli of hard-particle fluids from a mode-coupling theory with a retarded vertex

The recently proposed modification of the mode-coupling theory (MCT) in which the static structure used in the vertex is computed at a lower density than the actual one is tested on several dynamics-related properties. The predictions from this modified version of MCT calibrated on the one-component hard-sphere fluid are found in very good agreement with simulation data for one-component and binary hard-sphere fluids. They are also relevant for the stress moduli for models with attractive tails beyond the hard core. The clear improvement observed on several properties should give a new impetus to the use of MCT as a quantitative tool.

Slow dynamics coupled with cluster formation in ultrasoft-potential glasses

We numerically investigate the slow dynamics of a binary mixture of ultrasoft particles interacting with the generalized Hertzian potential. If the softness parameter, α, is small, the particles at high densities start penetrating each other, form clusters, and eventually undergo the glass transition. We find multiple cluster-glass phases characterized by a different number of particles per cluster, whose boundary lines are sharply separated by the cluster size. Anomalous logarithmic slow relaxation of the density correlation functions is observed in the vicinity of these glass-glass phase boundaries, which hints the existence of the higher-order dynamical singularities predicted by the mode-coupling theory. Deeply in the cluster glass phases, it is found that the dynamics of a single particle is decoupled from that of the collective fluctuations.

Microscopic Theory of Two-Step Yielding in Attractive Colloids

Attractive colloids display two distinct amorphous solid phases: the attractive glass due to particle bonding and the repulsive glass due to the hard-core repulsion. By means of a microscopic mean field approach, we analyze their response to a quasistatic shear strain. We find that the presence of two distinct interaction length scales may result in a sharp two-step yielding process, which can be associated with a hysteretic stress response or with a reversible but nonmonotonic stress-strain curve. We derive a generic phase diagram characterized by two distinct yielding lines, an inverse yielding and a critical point separating the hysteretic and reversible regimes. Our results should be applicable to a large class of glassy materials characterized by two distinct interaction length scales.

Tracer Transport Probes Relaxation and Structure of Attractive and Repulsive Glassy Liquids

Dynamic coupling of small penetrants to slow, cooperative relaxations within crowded cells, supercooled liquids, and polymer matrices has broad consequences for applications ranging from drug delivery to nanocomposite processing. Interactions between the constituents of these and other disordered media alter the cooperative relaxations, but their effect on penetrant dynamics remains incompletely understood. We use molecular dynamics simulations to show that the motions of hard-sphere tracer particles probe differences in local structure and cooperative relaxation processes in attractive and repulsive glassy liquid matrices with equal bulk packing fractions and long-time diffusivities. Coupling of the tracer dynamics to collective relaxations in each matrix affects the shape of tracer trajectories, which are fractal within the repulsive matrix and more compact in the attractive. These results reveal that the structure of relaxations controls penetrant transport and dispersion in cooperatively relaxing systems and provide insight into dynamical heterogeneity within glassy liquids.

New insights into re-entrant melting of microgel particles by polymer-induced aggregation experiments

While microgels are in general described as soft particles, polystyrene (PS) microgels can be synthesized in a way that cross-link density has only a minor influence on their physical properties. Even though the particles swell in a good solvent, the imparted slight softness still allows a mapping on hard sphere behaviour for a large range of cross-link densities [Schneider et al., Soft Matter, 2017, 13, 445]. Nevertheless, the hard sphere analogy breaks down as soon as polymer chains are added to these systems. Quantitative differences between PS microgels and true hard sphere systems appear and the differences between stronger and weaker cross-linked PS microgels can be observed. To gain deeper insights into the origin of these deviations from true hard sphere behaviour, this work is addressed to a systematic study of the colloid-polymer interactions in PS microgel-polymer mixtures. We investigated the aggregation behaviour (namely aggregation concentration and cluster structure) as a function of colloid size, cross-link density and colloid-polymer size ratio in very dilute colloidal suspensions. Our results show that the interplay of cross-link density and polymer size is a key parameter for the strength of the colloid-polymer interactions. Furthermore, the centre-to-centre distance of the colloidal particles in the formed clusters decreases if the cross-link density is decreased, allowing for a higher packing density. This may also explain the unusually high fluid packing fractions observed in the re-entry region of the phase diagram of PS microgel-free PS polymer mixtures.

Structure formation in soft nanocolloids: liquid-drop model

Using a model where soft nanocolloids such as spherical polymer brushes and star polymers are viewed as compressible liquid drops, we theoretically explore contact interactions between such particles. By numerically minimizing the phenomenological free energy consisting of bulk and surface terms, we find that at small deformations the drop-drop interaction is pairwise additive and described by a power law. We also propose a theory to describe the small-deformation regime, and the agreement is very good at all drop compressibilities. The large-deformation regime, which is dominated by many-body interactions, is marked by a rich phase diagram which includes the face- and body-centered-cubic, σ, A15, and simple hexagonal lattice as well as isostructural and re-entrant transitions. Most of these features are directly related to the non-convex deformation free energy emerging from many-body effects in the partial-faceting regime. The phase diagram, which depends on just two model parameters, contains many of the condensed phases observed in experiments. We also provide statistical-mechanical arguments that relate the two model parameters to the molecular architecture of the polymeric nanocolloids, chain rigidity, and solvent quality. The model represents a generic framework for the overarching features of the phase behavior of polymeric nanocolloids at high compressions.

Interplay between Caging and Bonding in Binary Concentrated Colloidal Suspensions

When a liquid becomes dynamically arrested, a gel, a repulsive glass, or an attractive glass state will form. Bonding and caging mechanisms decide their static structures and dynamic properties. To better understand their interplay, the competition between bonding and caging in a binary mixture of polystyrene core/poly( N-isopropylacrylamide) shell (CS) microgels and sulfonated polystyrene (PSS) particles is studied. CS microgels have short-range attraction above the volume phase transition temperature, whereas PSS species experiences relatively long-range electrostatic repulsion. Adding more PSS into the binary mixture will, of course, increase the total effective volume fraction but lead to different properties in gel or glass states. For instance, in gels, it increases the localization length and weakens the gel, whereas in glass, it decreases the localization length and strengthens the glass. This thus implies that the static and dynamic properties of gels are mainly controlled by bonding and those of both repulsive and attractive glasses are governed by caging. On the other hand, increasing the temperature will decrease the effective volume fraction because of the volume phase transition of the CS microgels. A discontinuous repulsive glass-to-liquid-to-gel transition can be observed when the PSS concentration is low, but a continuous repulsive glass-to-gel transition can also be observed with the increase of the PSS concentration. This may hint that glass transition and physical gelation share a similar mechanism, whereas the former has a longer relaxation time.

Static structure of active Brownian hard disks

We explore the changes in static structure of a two-dimensional system of active Brownian particles (ABP) with hard-disk interactions, using event-driven Brownian dynamics simulations. In particular, the effect of the self-propulsion velocity and the rotational diffusivity on the orientationally-averaged fluid structure factor is discussed. Typically activity increases structural ordering and generates a structure factor peak at zero wave vector which is a precursor of motility-induced phase separation. Our results provide reference data to test future statistical theories for the fluid structure of active Brownian systems. This manuscript was submitted for the special issue of the Journal of Physics: Condensed Matter associated with the Liquid Matter Conference 2017.

Describing magnetorheology under a colloidal glass approach

The equilibrium structure and dynamics of magnetorheological (MR) fluids are studied in this work by simulations, where particles are modeled as dipoles with a quasihard spherical core. Upon increasing the interaction strength, controlled experimentally by the magnetic field, elongated clusters grow and, for intense fields, thick columns form, aligned with the field. The dynamics of the system is monitored by the mean-squared displacement and density correlation functions, which show an increasing slowing down with the attraction strength. The correlation function shows a two-step decay, with a separation between microscopic and long time dynamics, a typical hallmark of undercooled fluids. We have therefore analyzed the dynamics of this MR fluid using the typical concepts for undercooled fluids. Thus, the second decay of the density correlation function is fitted with a stretched exponential, and the wave-vector dependence of the fitting parameters studied. Both the amplitude and the time scale oscillate in phase with the structure factor. Our results support the idea that the magnetorheological effect is in fact the manifestation of a colloidal system approaching an attractive glass transition (or gel transition).

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

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

Glass transition in hard-core fluids and beyond, using an effective static structure in the mode coupling theory

The dynamical arrest in classical fluids is studied using a simple modification of the mode coupling theory (MCT) aimed at correcting its overestimation of the tendency to glass formation while preserving its overall structure. As in previous attempts, the modification is based on the idea of tempering the static pair correlations used as input. It is implemented in this work by computing the static structure at a different state point than the one used to solve the MCT equation for the intermediate scattering function, using the pure hard-sphere glass for calibration. The location of the glass transition predicted from this modification is found to agree with simulations data for a variety of systems --pure fluids and mixtures with either purely repulsive interaction potentials or ones with attractive contributions. Besides improving the predictions in the long-time limit, and so reducing the non-ergodicity domain, the same modification works as well for the time-dependent correlators.

Cluster Glass Transition of Ultrasoft-Potential Fluids at High Density

Using molecular dynamics simulation, we investigate the slow dynamics of a supercooled binary mixture of soft particles interacting with a generalized Hertzian potential. At low density, it displays typical slow dynamics near its glass transition temperature. At higher densities, particles bond together, forming clusters, and the clusters undergo the glass transition. The number of particles in a cluster increases one by one as the density increases. We demonstrate that there exist multiple cluster-glass phases characterized by a different number of particles per cluster, each of which is separated by distinct minima. Surprisingly, a so-called higher order singularity of the mode-coupling theory signaled by a logarithmic relaxation is observed in the vicinity of the boundaries between monomer and cluster glass phases. The system also exhibits rich and anomalous dynamics in the cluster glass phases, such as the decoupling of the self- and collective dynamics.

Gelation transitions of colloidal systems with bridging attractions

Gelation transitions in a colloidal system, where there is a strong reversible attraction between small, soft microgels and large, hard spheres, are systematically investigated. Different from widely studied depletion attraction systems that are also two-component systems, the strong attraction between small solvent and large solute particles introduces bridging attractions between large solute particles. We conclusively demonstrate that the formation of physical gels at the intermediate volume fraction of our bridging attraction system follows more closely with the percolation line that is in stark contrast to what is observed in depletion attraction systems, where the gelation transition is related with the frustrated spinodal separation, not a purely kinetic phenomenon. Our results introduce a different way to control gelation transitions in spherical colloidal systems, and imply that people need to be prudent when generalizing the physical picture of the gelation transitions obtained from systems with different origins of effective attraction as the solvent molecule may play important roles.

Vibrational properties of quasi-two-dimensional colloidal glasses with varying interparticle attraction

We measure the vibrational modes and particle dynamics of quasi-two-dimensional colloidal glasses as a function of interparticle interaction strength. The interparticle attractions are controlled via a temperature-tunable depletion interaction. Specifically, the interparticle attraction energy is increased gradually from a very small value (nearly hard-sphere) to moderate strength (∼4k_{B}T), and the variation of colloidal particle dynamics and vibrations are concurrently probed. The particle dynamics slow monotonically with increasing attraction strength, and the particle motions saturate for strengths greater than ∼2k_{B}T, i.e., as the system evolves from a nearly repulsive glass to an attractive glass. The shape of the phonon density of states is revealed to change with increasing attraction strength, and the number of low-frequency modes exhibits a crossover for glasses with weak compared to strong interparticle attraction at a threshold of ∼2k_{B}T. This variation in the properties of the low-frequency vibrational modes suggests a new means for distinguishing between repulsive and attractive glass states.

Microstructure, local dynamics, and flow behavior of colloidal suspensions with weak attractive interactions

We present a comprehensive micro- and macrorheological study of the effect of weak depletion attraction (Ψdep ≈ 1-10 kBT) on dense colloidal suspensions stabilized by short-range repulsive interactions. We used aqueous polymer dispersions as model system and demonstrated the unique capabilities of multiple particle tracking (MPT) to disclose structural changes in such technically important systems exhibiting many characteristic features of hard sphere systems. Below the hard sphere freezing point ϕc, viscosity increases monotonically with increasing Ψdep due to the transition from a fluid to a fluid/crystalline and finally to a gel state. Above ϕc, increasing attraction strength first results in a viscosity reduction corresponding to the formation of large, permeable crystals and then in a viscosity increase when a network of dense, small crystals forms. The fraction of the fluid and crystal phase, particle concentration in each phase as well as the modulus of the micro-crystals are obtained, the latter decreases with Ψdep. Above the colloidal glass transition strong heterogeneities and different local particle mobility in the repulsive and attractive arrested states are found. Particles are trapped in the cage of neighboring particles rather than in an attractive potential well. The intermediate ergodic state exhibits uniform tracer diffusivity.

Glassy dynamics of Brownian particles with velocity-dependent friction

We consider a two-dimensional model system of Brownian particles in which slow particles are accelerated while fast particles are damped. The motion of the individual particles is described by a Langevin equation with Rayleigh-Helmholtz velocity-dependent friction. In the case of noninteracting particles, the time evolution equations lead to a non-Gaussian velocity distribution. The velocity-dependent friction allows negative values of the friction or energy intakes by slow particles, which we consider active motion, and also causes breaking of the fluctuation dissipation relation. Defining the effective temperature proportional to the second moment of velocity, it is shown that for a constant effective temperature the higher the noise strength, the lower the number of active particles in the system. Using the Mori-Zwanzig formalism and the mode-coupling approximation, the equations of motion for the density autocorrelation function are derived. The equations are solved using the equilibrium structure factors. The integration-through-transients approach is used to derive a relation between the structure factor in the stationary state considering the interacting forces, and the conventional equilibrium static structure factor.

Scaling and universality in glass transition

Kinetic facilitated models and the Mode Coupling Theory (MCT) model B are within those systems known to exhibit a discontinuous dynamical transition with a two step relaxation. We consider a general scaling approach, within mean field theory, for such systems by considering the behavior of the density correlator 〈q(t)〉 and the dynamical susceptibility 〈q(2)(t)〉 - 〈q(t)〉(2). Focusing on the Fredrickson and Andersen (FA) facilitated spin model on the Bethe lattice, we extend a cluster approach that was previously developed for continuous glass transitions by Arenzon et al. (Phys. Rev. E 90, 020301(R) (2014)) to describe the decay to the plateau, and consider a damage spreading mechanism to describe the departure from the plateau. We predict scaling laws, which relate dynamical exponents to the static exponents of mean field bootstrap percolation. The dynamical behavior and the scaling laws for both density correlator and dynamical susceptibility coincide with those predicted by MCT. These results explain the origin of scaling laws and the universal behavior associated with the glass transition in mean field, which is characterized by the divergence of the static length of the bootstrap percolation model with an upper critical dimension dc = 8.