Glasses in hard spheres with short-range attraction

Approaching the hard sphere limit in colloids suitable for confocal microscopy - the end of a decade lasting quest

PMMA-PHSA particles serve as the hard sphere model system since the 1980s. We investigate the fluid structure of fluorescent ones in three different solvents by laser scanning confocal microscopy: a decalin-tetrachloroethylene (TCE)-mixture and a decalin-cyclohexylbromide (CHB)-mixture with and without tetrabutylammoniumbromide (TBAB). The experimental 3D radial distribution functions are modeled by analytical theory and computer simulations taking polydispersity and the experimental position uncertainty into account. The quantitative comparison between experiment and simulation/theory establishes hard sphere like behavior for particles in decalin-TCE for a wide range of particle packing fractions. To the best of our knowledge, we present the first experimental dataset of a fluid structure that agrees convincingly with Percus-Yevick over a wide concentration range. Furthermore, charged sphere behavior is confirmed both for the decalin-CHB and the decalin-CHB-TBAB solvents, and it is demonstrated that a finite particle concentration reduces screening in the decalin-CHB-TBAB system compared to the bulk solvent.

Memory in aging colloidal gels with time-varying attraction

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

Colloidal Rod-Like Particles with Temperature-Driven Tunable Interactions

Temperature-sensitive rod-like colloidal particles were synthesized by grafting a temperature-responsive polymer, poly(2-(dimethylamino)ethyl methacrylate) (PDMA), on the surface of high aspect ratio silica rods by surface-initiated atom transfer radical polymerization. The stability of the grafted polymer on the surface of the particles in aqueous solutions was found to deteriorate with time, leading to a gradual decrease of the polymer content of the hybrid colloids, which was attributed to the mechanically activated hydrolysis of the labile bonds at the polymer-silica interface. The polymer degrafting was significantly suppressed by first growing a hydrophobic poly(methyl methacrylate) block onto the particle surface to act as a barrier layer for the penetration of water molecules at the polymer-particle interface, followed by chain-extension with the hydrophilic PDMA chains. Dynamic light scattering, microscopy, and rheological measurements revealed that the PDMA block conferred a temperature-responsive behavior to the rod-like particles, which formed aggregates at temperatures above the lower critical solution temperature (LCST) of the polymer. However, in contrast to their spherical counterparts, the polymer-grafted rod-like particles did not exhibit complete thermo-reversibility upon lowering the solution temperature below the LCST of PDMA, which was reflected by different values of the diffusion coefficient for the heating and cooling cycles, indicating an irreversible rod particle aggregation upon increasing the temperature.

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.

Depletion Induced Demixing and Crystallization in Binary Colloids Subjected to an External Potential Barrier

Depletion interaction plays an important role in determining the structural and dynamical properties of binary colloidal mixtures. We have investigated the effect of the attractive depletion interaction between an external potential barrier and larger species in the binary mixture on the phase behavior of a binary colloidal mixture using canonical-isokinetic ensemble molecular dynamics simulations. The demixing of the binary mixture due to this depletion interaction increases as the volume fraction increases, and a pure phase of larger particles forms in the region of the potential barrier. The local density of this pure phase is high enough that a face centered cubic crystalline domain is formed at this region. This crystalline phase diffuses perpendicular to the external potential barrier, indicating that moving crystals can be obtained in an equilibrium system. The temperature dependence of diffusivity of larger particles is non-Arrhenius and changes from sub-Arrhenius to super-Arrhenius as the volume fraction increases. This crossover from sub-Arrhenius to super-Arrhenius diffusion coincides with the crystalline formation near the potential barrier.

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.

Structural and short-time vibrational properties of colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition

We investigate the short-time vibrational properties and structure of two-dimensional, bidisperse, colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition, as a function of interparticle depletion attraction strength. The long-time spatiotemporal dynamics of the samples are measured to be non-monotonic, confirming that the suspensions evolve from repulsive glass to supercooled liquid to attractive glass with increasing depletion attraction. Here, we search for vibrational signatures of the re-entrant behavior in the short-time spatiotemporal dynamics, i.e., dynamics associated with particle motion inside its nearest-neighbor cage. Interestingly, we observe that the anharmonicity of these in-cage vibrations varies non-monotonically with increasing attraction strength, consistent with the non-monotonic long-time structural relaxation dynamics of the re-entrant glass. We also extract effective spring constants between neighboring particles; we find that spring stiffness involving small particles also varies non-monotonically with increasing attraction strength, while stiffness between large particles increases monotonically. Last, from study of depletion-dependent local structure and vibration participation fractions, we gain microscopic insight into the particle-size-dependent contributions to short-time vibrational modes in the glass and supercooled liquid states.

Glassy states in adsorbing surfactant-microgel soft nanocomposites

Mixtures of polymer-colloid hybrids such as star polymers and microgels with non-adsorbing polymeric additives have received a lot of attention. In these materials, the interplay between entropic forces and softness is responsible for a wealth of phenomena. By contrast, binary mixtures where one component can adsorb onto the other one have been far less studied. Yet real formulations in applications often contain low molecular weight additives that can adsorb onto soft colloids. Here we study the microstructure and rheology of soft nanocomposites made of surfactants and microgels using linear and nonlinear rheology, SAXS experiments, and cryo-TEM techniques. The results are used to build a dynamical state diagram encompassing various liquid, glassy, jammed, metastable, and reentrant liquid states, which results from a subtle interplay between enthalpic, entropic, and kinetic effects. We rationalize the rheological properties of the nanocomposites in each domain of the state diagram, thus providing exquisite solutions for designing new rheology modifiers at will.

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.

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.

Structural and rheological aging in model attraction-driven glasses by Rheo-SANS

Aging in a model colloidal suspension comprised of particles with a thermoreversible attraction is studied using Rheo-SANS techniques in the attractive-driven glass state. Multiple thermal pathways lead to a common rheological and microstructural aging trajectory, as was observed previously for a thermoreversible gel. SANS measurements of the colloidal glass microstructure as a function of temperature and time during various quench protocols are quantitatively characterized in terms of an effective interaction strength that becomes an order parameter defining the microstructural state of the glass. Using previously validated concepts of a fictive temperature, a semi-empirical, quantitative relationship similar to an Avrami relationship is established between the mechanical aging (elastic modulus) and microstructural aging (order parameter) that is independent of thermal history for the thermal profiles studied herein at long times. Furthermore, shear rejuvenation is studied, and while shear may only partially reduce the degree of structure in the glass, aging upon flow cessation is found to follow a common trajectory when viewed in terms of the microstructural order parameter.

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.

Scale invariance in a nonvibrating magnetic granular system

A nonvibrating magnetic granular system is studied by using a time series approach. The system consists of steel balls confined inside a circular wall that surrounds a glass plate. Kinetic energy is provided to the particles by the application of an external vertical time-dependent magnetic field of different amplitudes. We carried out a characterization of the system dynamics through the measurement of the correlations present in the time series of positions, in the x-direction, of each particle. In particular, by performing Fourier spectral analysis, we find that the time series are fractal and scale invariant, in such a way that the corresponding Fourier power spectra follow a power law \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P(f)\propto 1/f^\beta$$\end{document}P(f)∝1/fβ, with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0<\beta <2.5$$\end{document}0<β<2.5. More specifically, we find that the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β, and therefore the strength of the correlations, increase as the magnetic field also increases. In this way, the present system constitutes an experimental model to generate correlated random walks. Additionally, we show how the introduction of a constant magnetic field breaks down this scale invariance property in the positions of each particle. Finally, we confirm the above results by applying detrended fluctuation analysis.

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.

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 and reversible gelation in asymmetric binary mixtures: A study by mode coupling theory and molecular dynamics

The glass transition and the binodals of asymmetric binary mixtures are investigated from the effective fluid approach in the mode coupling theory and by molecular dynamics. Motivated by previous theoretical predictions, the hard-sphere mixture and the Asakura-Oosawa models are used to analyze experimental results from the literature, relative to polystyrene spheres mixed either with linear polymers or with dense microgel particles. In agreement with the experimental observations, the specificity of the depletant particles is shown to favor lower density gels. It further favors equilibrium gelation by reducing also the tendency of the system to phase separate. These results are confirmed by a phenomenological modification of the mode coupling theory in which the vertex functions are computed at an effective density lower than the actual one. A model effective potential in asymmetric mixtures of hard particles is used to further check this phenomenological modification against molecular dynamics simulation.

Dynamic Heterogeneities in Colloidal Supercooled Liquids: Experimental Tests of Inhomogeneous Mode Coupling Theory

The dynamics in supercooled liquids slow enormously upon approaching the glass transition, albeit without significant change of liquid structure. This empirical observation has stimulated development of many theoretical models which attempt to elucidate microscopic mechanisms in glasses and glass precursors. Here, quasi-two-dimensional colloidal supercooled liquids and glasses are employed to experimentally test predictions of widely used models: mode coupling theory (MCT) and its important extension, inhomogeneous MCT (IMCT). We measure two-point dynamic correlation functions in the glass forming liquids to determine structural relaxation times, τ_α, and mode coupling exponents, a, b, and γ; these parameters are then used to extract the mode coupling dynamic crossover packing area-fraction, ϕ _c. This information, along with our measurements of supercooled liquid spatiotemporal dynamics, permits characterization of dynamic heterogeneities in the samples and facilitates direct experimental tests of the scaling predictions of IMCT. The time scales at which dynamic heterogeneities are largest, and their spatial sizes, exhibit power law growth on approaching ϕ _c. Within experimental error, the exponents of the measured power laws are close to the predictions of IMCT.

Glassy dynamics in asymmetric binary mixtures of hard spheres

We perform a systematic and detailed study of the glass transition in highly asymmetric binary mixtures of colloidal hard spheres, combining differential dynamic microscopy experiments, event-driven molecular dynamics simulations, and theoretical calculations, exploring the whole state diagram and determining the self-dynamics and collective dynamics of both species. Two distinct glassy states involving different dynamical arrest transitions are consistently described, namely, a double glass with the simultaneous arrest of the self-dynamics and collective dynamics of both species, and a single glass of large particles in which the self-dynamics of the small species remains ergodic. In the single-glass scenario, spatial modulations in the collective dynamics of both species occur due to the structure of the large spheres, a feature not observed in the double-glass domain. The theoretical results, obtained within the self-consistent generalized Langevin equation formalism, are in agreement with both simulations and experimental data, thus providing a stringent validation of this theoretical framework in the description of dynamical arrest in highly asymmetric mixtures. Our findings are summarized in a state diagram that classifies the various amorphous states of highly asymmetric mixtures by their dynamical arrest mechanisms.

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.

Gravitational collapse of colloidal gels: non-equilibrium phase separation driven by osmotic pressure

Delayed gravitational collapse of colloidal gels is characterized by initially slow compaction that gives way to rapid bulk collapse, posing interesting questions about the underlying mechanistic origins. Here we study gel collapse utilizing large-scale dynamic simulation of a freely draining gel of physically bonded particles subjected to gravitational forcing. The hallmark regimes of collapse are recovered: slow compaction, transition to rapid collapse, and long-time densification. Microstructural changes are monitored by tracking particle positions, coordination number, and bond dynamics, along with volume fraction, osmotic pressure, and potential energy. Together these reveal the surprising result that collapse can occur with a fully intact network, where the tipping point arises when particle migration dissolves strands in a capillary-type instability. While it is possible for collapse to rupture a gel network into clusters that then sediment, and hydrodynamic interactions can make interesting contributions, neither is necessary. Rather, we find that the "delay" arises from gravity-enhanced coarsening, which triggers the re-emergence of phase separation. The mechanism of this transition is a leap toward lower potential energy of the gel, driven by bulk negative osmotic pressure that condenses the particle phase: the gel collapses in on itself under negative osmotic pressure allowing the gel, to tunnel through the equilibrium phase diagram to a higher volume fraction "state". Remarkably, collapse stops when condensation stops, when gravitational advection produces a positive osmotic pressure, re-arresting the gel.

Dynamics of fluids in quenched-random potential energy landscapes: a mode-coupling theory approach

Motivated by a number of recent experimental and computational studies of the dynamics of fluids plunged in quenched-disordered external fields, we report on a theoretical investigation of this topic within the framework of the mode-coupling theory, based on the simple model of the hard-sphere fluid in a Gaussian random field. The possible dynamical arrest scenarios driven by an increase of the disorder strength and/or of the fluid density are mapped, and the corresponding evolutions of time-dependent quantities typically used for the characterization of anomalous self-diffusion are illustrated with detailed computations. Overall, a fairly reasonable picture of the dynamics of the system at hand is outlined, which in particular involves a non-monotonicity of the self-diffusion coefficient with fluid density at fixed disorder strength, in agreement with experiments. The disorder correlation length is shown to have a strong influence on the latter feature.

Long-term aging behaviors in a model soft colloidal system

Colloidal and molecular systems share similar behaviors near to the glass transition volume fraction or temperature. Here, aging behaviors after volume fraction up-jump (induced by performing temperature down-jumps) conditions for a PS-PNIPAM/AA soft colloidal system were investigated using light scattering (diffusing wave spectroscopy, DWS). Both aging responses and equilibrium dynamics were investigated. For the aging responses, long-term experiments (100 000 s) were performed, and both equilibrium and non-equilibrium behaviors of the system were obtained. In the equilibrium state, as effective volume fraction increases (or temperature decreases), the colloidal dispersion displays a transition from the liquid to a glassy state. The equilibrium α-relaxation dynamics strongly depend on both the effective volume fraction and the initial mass concentration for the studied colloidal systems. Compared with prior results from our lab [X. Di, X. Peng and G. B. McKenna, J. Chem. Phys., 2014, 140, 054903], the effective volume fractions investigated spanned a wider range, to deeper into the glassy domain. The results show that the α-relaxation time τ_α of the samples aged into equilibrium deviate from the classical Vogel-Fulcher-Tammann (VFT)-type expectations and the super-Arrhenius signature disappears above the glass transition volume fraction. The non-equilibrium aging response shows that the time for the structural evolution into equilibrium and the α-relaxation time are decoupled. The DWS investigation of the aging behavior after different volume fraction jumps reveals a different non-equilibrium or aging behavior for the considered colloidal systems compared with either molecular glasses or the macroscopic rheology of a similar colloidal dispersions.

On tuning microgel character and softness of cross-linked polystyrene particles

Polystyrene (PS) microgel colloids have often been used successfully to model hard sphere behaviour even though the term "gel" invokes inevitably the notion of a more or less soft, deformable object. Here we systematically study the effect of reducing the cross-link density from 1 : 10 (1 cross-link per 10 monomers) to 1 : 100 on particle interactions and "softness". We report on the synthesis and purification of 1 : 10, 1 : 25, 1 : 50, 1 : 75 and 1 : 100 cross-linked PS particles and their characterization in terms of single particle properties, as well as the behaviour of concentrated dispersions. We are able to tune particle softness in the range between soft PNiPAM-microgels and hard PMMA particles while still allowing the mapping of the microgels onto hard sphere behavior with respect to phase diagram and static structure factors. This is mainly due to a rather homogeneous radial distribution of cross-links in contrast to PNiPAM microgels where the cross-link density decreases radially. We find that up to a cross-link density of 1 : 50 particle form factors are perfectly described by a homogeneous sphere model whereas 1 : 75 and 1 : 100 cross-linked spheres are slightly better described as fuzzy spheres. However the fuzziness is rather small compared to typical PNiPAM microgels so that a hard sphere mapping still holds even for these low cross-link densities. Finally, by varying the reaction conditions - changing from batch to semibatch emulsion polymerization and varying the feed rate or by adjusting the monomer to initiator ratio we can tune the fuzziness or significantly alter the inner structure to a more open, star-like architecture.

Connectivity, dynamics, and structure in a tetrahedral network liquid

We report a detailed computational study by Brownian dynamics simulations of the structure and dynamics of a liquid of patchy particles which forms an amorphous tetrahedral network upon decreasing the temperature. The highly directional particle interactions allow us to investigate the system connectivity by discriminating the total set of particles into different populations according to a penta-modal distribution of bonds per particle. With this methodology we show how the particle bonding process is not randomly independent but it manifests clear bond correlations at low temperatures. We further explore the dynamics of the system in real space and establish a clear relation between particle mobility and particle connectivity. In particular, we provide evidence of anomalous diffusion at low temperatures and reveal how the dynamics is affected by the short-time hopping motion of the weakly bounded particles. Finally we widely investigate the dynamics and structure of the system in Fourier space and identify two quantitatively similar length scales, one dynamic and the other static, which increase upon cooling the system and reach distances of the order of few particle diameters. We summarize our findings in a qualitative picture where the low temperature regime of the viscoelastic liquid is understood in terms of an evolving network of long time metastable cooperative domains of particles.

Nonvibrating granular model for a glass-forming liquid: Equilibration and aging

We studied experimentally a model of a glass-forming liquid on the basis of a nonvibrating magnetic granular system under an unsteady magnetic field. A sudden quenching was produced that drove the system from a liquid state to a different final state with lower temperature; the latter could be a liquid state or a solid state. We determined the mean-squared displacement in temporal windows to obtain the dynamic evolution of the system, and we determined the radial distribution function to obtain its structural characteristics. The results were analyzed using the intermediate scattering function and the effective potential between two particles. We observed that when quenching drives the system to a final state in the liquid phase far from the glass-transition temperature, equilibration occurs very quickly. When the final state has a temperature far below the glass-transition temperature, the system reaches its equilibrium state very quickly. In contrast, when the final state has an intermediate temperature but is below that corresponding to the glass transition, the system falls into a state that evolves slowly, presenting aging. The system evolves by an aging process toward more ordered states. However, after a waiting time, the dynamic behavior changes. It was observed that some particles get close enough to overpass the repulsive interactions and form small stable aggregates. In the effective potential curves, it was observed that the emergence of a second effective well due to the attraction quickly evolves and results in a deeper well than the initial effective well due to the repulsion. With the increase in time, more particles fall in the attractive well forming inhomogeneities, which produce a frustration in the aging process.

Structure of colloidal gels at intermediate concentrations: the role of competing interactions

Colloidal gels formed by colloid-polymer mixtures with an intermediate volume fraction (ϕ_c ≈ 0.4) are investigated by confocal microscopy. In addition, we have performed Monte Carlo simulations based on a simple effective pair potential that includes a short-range attractive contribution representing depletion interactions, and a longer-range repulsive contribution describing the electrostatic interactions due to the presence of residual charges. Despite neglecting non-equilibrium effects, experiments and simulations yield similar gel structures, characterised by, e.g., the pair, angular and bond distribution functions. We find that the structure hardly depends on the strength of the attraction if the electrostatic contribution is fixed, but changes significantly if the electrostatic screening is changed. This delicate balance between attractions and repulsions, which we quantify by the second virial coefficient, also determines the location of the gelation boundary.

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.

Tunable depletion potentials driven by shape variation of surfactant micelles

Depletion interaction potentials between micron-sized colloidal particles are induced by nanometer-scale surfactant micelles composed of hexaethylene glycol monododecyl ether (C_{12}E_{6}), and they are measured by video microscopy. The strength and range of the depletion interaction is revealed to arise from variations in shape anisotropy of the surfactant micelles. This shape anisotropy increases with increasing sample temperature. By fitting the colloidal interaction potentials to theoretical models, we extract micelle length and shape anisotropy as a function of temperature. This work introduces shape anisotropy tuning as a means to control interparticle interactions in colloidal suspensions, and it shows how the interparticle depletion potentials of micron-scale objects can be employed to probe the shape and size of surrounding macromolecules at the nanoscale.

Self-diffusion and structural properties of confined fluids in dynamic coexistence

Self-diffusion and radial distribution functions are studied in a strongly confined Lennard-Jones fluid. Surprisingly, in the solid-liquid phase transition region, where the system exhibits dynamic coexistence, the self-diffusion constants are shown to present up to three-fold variations from solid to liquid phases at fixed temperature, while the radial distribution function corresponding to both the liquid and the solid phases are essentially indistinguishable.

Non-monotonic effect of confinement on the glass transition

The relaxation dynamics of glass forming liquids and their structure are influenced in the vicinity of confining walls. This effect has mostly been observed to be a monotonic function of the slit width. Recently, a qualitatively new behaviour has been uncovered by Mittal and coworkers, who reported that the single particle dynamics in a hard-sphere fluid confined in a planar slit varies in a non-monotonic way as the slit width is decreased from five to roughly two particle diametres (Mittal et al 2008 Phys. Rev. Lett. 100 145901). In view of the great potential of this effect for applications in those fields of science and industry, where liquids occur under strong confinement (e.g. nano-technology), the number of researchers studying various aspects and consequences of this non-monotonic behaviour has been rapidly growing. This review aims at providing an overview of the research activity in this newly emerging field. We first briefly discuss how competing mechanisms such as packing effects and short-range attraction may lead to a non-monotonic glass transition scenario in the bulk. We then analyse confinement effects on the dynamics of fluids using a thermodynamic route which relates the single particle dynamics to the excess entropy. Moreover, relating the diffusive dynamics to the Widom's insertion probability, the oscillations of the local dynamics with density at moderate densities are fairly well described. At high densities belonging to the supercooled regime, however, this approach breaks down signaling the onset of strongly collective effects. Indeed, confinement introduces a new length scale which in the limit of high densities and small pore sizes competes with the short-range local order of the fluid. This gives rise to a non-monotonic dependence of the packing structure on confinement, with a corresponding effect on the dynamics of structural relaxation. This non-monotonic effect occurs also in the case of a cone-plate type channel, where the degree of confinement varies with distance from the apex. This is a very promising issue for future research with the possibility of uncovering the existence of alternating glassy and liquid-like domains.

Discontinuous nature of the repulsive-to-attractive colloidal glass transition

In purely repulsive colloidal systems a glass transition can be reached by increasing the particle volume fraction beyond a certain threshold. The resulting glassy state is governed by configurational cages which confine particles and restrict their motion. A colloidal glass may also be formed by inducing attractive interactions between the particles. When attraction is turned on in a repulsive colloidal glass a re-entrant solidification ensues. Initially, the repulsive glass melts as free volume in the system increases. As the attraction strength is increased further, this weakened configurational glass gives way to an attractive glass in which motion is hindered by the formation of physical bonds between neighboring particles. In this paper, we study the transition from repulsive-to-attractive glasses using three-dimensional imaging at the single-particle level. We show how the onset of cage weakening and bond formation is signalled by subtle changes in local structure. We then demonstrate the discontinuous nature of the solid-solid transition, which is marked by a critical onset at a threshold bonding energy. Finally, we highlight how the interplay between bonding and caging leads to complex and heterogeneous dynamics at the microscale.

Gelation in mixtures of polymers and bidisperse colloids

We investigated the effects of varying the volume fraction of large particles (r) on the linear rheology and microstructure of mixtures of polymers and bidisperse colloids, in which the ratio of the small and large particle diameters was α=0.31 or α=0.45. Suspensions formulated at a total volume fraction of ϕ_{T}=0.15 and a constant concentration of polymer in the free volume c/c^{*}≈0.7 contained solid-like gels for small r and fluids or fluids of clusters at large r. The solid-like rheology and microstructure of these suspensions changed little with r when r was small, and fluidized only when r>0.8. By contrast, dense suspensions with ϕ_{T}=0.40 and α=0.31 contained solid-like gels at all concentrations of large particles and exhibited only modest rheological and microstructural changes upon varying the volume fraction of large particles. These results suggest that the effect of particle-size dispersity on the properties of colloid-polymer mixtures are asymmetric in particle size and are most pronounced near a gelation boundary.

Avalanche-like fluidization of a non-Brownian particle gel

We report on the fluidization dynamics of an attractive gel composed of non-Brownian particles made of fused silica colloids. Extensive rheology coupled to ultrasonic velocimetry allows us to characterize the global stress response together with the local dynamics of the gel during shear startup experiments. In practice, after being rejuvenated by a preshear, the gel is left to age for a time tw before being subjected to a constant shear rate [small gamma, Greek, dot above]. We investigate in detail the effects of both tw and [small gamma, Greek, dot above] on the fluidization dynamics and build a detailed state diagram of the gel response to shear startup flows. The gel may display either transient shear banding towards complete fluidization or steady-state shear banding. In the former case, we unravel that the progressive fluidization occurs by successive steps that appear as peaks on the global stress relaxation signal. Flow imaging reveals that the shear band grows until complete fluidization of the material by sudden avalanche-like events which are distributed heterogeneously along the vorticity direction and correlated to large peaks in the slip velocity at the moving wall. These features are robust over a wide range of tw and [small gamma, Greek, dot above] values, although the very details of the fluidization scenario vary with [small gamma, Greek, dot above]. Finally, the critical shear rate [small gamma, Greek, dot above]* that separates steady-state shear-banding from steady-state homogeneous flow depends on the width of the shear cell and exhibits a nonlinear dependence with tw. Our work brings about valuable experimental data on transient flows of attractive dispersions, highlighting the subtle interplay between shear, wall slip and aging whose modeling constitutes a major challenge that has not been met yet.

Unusual dynamics of concentration fluctuations in solutions of weakly attractive globular proteins

The globular protein γB-crystallin exhibits a complex phase behavior, where liquid-liquid phase separation characterized by a critical volume fraction ϕc = 0.154 and a critical temperature Tc = 291.8 K coexists with dynamical arrest on all length scales at volume fractions around ϕ ≈ 0.3-0.35, and an arrest line that extends well into the unstable region below the spinodal. However, although the static properties such as the osmotic compressibility and the static correlation length are in quantitative agreement with predictions for binary liquid mixtures, this is not the case for the dynamics of concentration fluctuations described by the dynamic structure factor S(q,t). Using a combination of dynamic light scattering and neutron spin echo measurements, we demonstrate that the competition between critical slowing down and dynamical arrest results in a much more complex wave vector dependence of S(q,t) than previously anticipated.

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.

Bond lifetime and diffusion coefficient in colloids with short-range interactions

We use molecular dynamics simulations to study the influence of short-range structures in the interaction potential between hard-sphere-like colloidal particles. Starting from model potentials and effective potentials in binary mixtures computed from the Ornstein-Zernike equations, we investigate the influence of the range and strength of a possible tail beyond the usual core repulsion or the presence of repulsive barriers. The diffusion coefficient and mean "bond" lifetimes are used as indicators of the effect of this structure on the dynamics. The existence of correlations between the variations of these quantities with the physical parameters is discussed to assess the interpretation of dynamics slowing down in terms of long-lived bonds. We also discuss the question of a universal behaviour determined by the second virial coefficient B ((2)) and the interplay of attraction and repulsion. While the diffusion coefficient follows the B ((2)) law for purely attractive tails, this is no longer true in the presence of repulsive barriers. Furthermore, the bond lifetime shows a dependence on the physical parameters that differs from that of the diffusion coefficient. This raises the question of the precise role of bonds on the dynamics slowing down in colloidal gels.

Concentration-dependent effective attractions between PEGylated nanoparticles

Effective attractions between colloidal particles bearing a grafted PEG layer in water have been studied and, from a quantitative SANS analysis, are found to be concentration dependent.

Brownian motion of optically anisotropic spherical particles in polymeric suspensions

We studied the rotational and translational diffusion of optically anisotropic liquid crystal particles embedded in semidiluted polymer solutions of Poly-Ethylene-Oxide (PEO) at different concentrations and different molecular weights. The polymer radius of gyration was chosen to be similar to the size of the probe particles and the polymer concentrations used are just above the crossover concentration. Thus, the system consists of solid probe particles moving in a sea of overlapping particles of similar size. We found that the behavior of both particle dynamics, rotational and translational, is similar in the range of concentrations considered here. In both cases, two linear diffusive regimes are observed, separated by a subdiffusive time interval. The spatial scale at which this intermediate regime appears shows a dependence on both the polymer concentration and molecular weight, and has a value similar to the thickness of the polymer-depleted layer usually found in this kind of systems. Additionally, we observe that the colloidal dynamic scales with the overlapping degree of the polymer particles.

Multiple glass singularities and isodynamics in a core-softened model for glass-forming systems

We investigate the slow dynamics of a simple glass former whose interaction potential is the sum of a hard core and a square shoulder repulsion. According to mode coupling theory, the competition between the two repulsive length scales gives rise to a complex dynamic scenario: besides the fluid-glass line, the theory predicts a glass-glass line in the temperature-packing fraction plane with two end points. Interestingly, for critical values of the square-shoulder parameters, such end points can be accessed from the liquid phase. We verify, via extensive numerical simulations, the existence of both points through the observation of an unconventional subdiffusive/logarithmic dynamical behavior. Unexpectedly, we also discover that the simultaneous presence of two end points generates special loci in the state diagram along which the dynamics is identical at all length and time scales.

Crystallization kinetics of colloidal model suspensions: recent achievements and new perspectives

Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions, with rather well-characterized pair interactions and minimized external influences. In complementary approaches experiment, analytic theory and simulation have been employed to study colloidal solidification in great detail. These studies were based on advanced optical methods, careful system characterization and sophisticated numerical methods. Over the last decade, both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity have been addressed in a quantitative way. Key parameters of crystallization have been derived and compared to those of metal systems. These systematic investigations significantly contributed to an enhanced understanding of the crystallization processes in general. Further, new fundamental questions have arisen and (partially) been solved over the last decade: including, for example, a two-step nucleation mechanism in homogeneous nucleation, choice of the crystallization pathway, or the subtle interplay of boundary conditions in heterogeneous nucleation. On the other hand, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting microstructure be influenced. The present review attempts to cover the interesting developments that have occurred since the turn of the millennium and to identify important novel trends, with particular focus on experimental aspects.

Solidification of a colloidal hard sphere like model system approaching and crossing the glass transition

We investigated the process of vitrification and crystallization in a model system of colloidal hard spheres. The kinetics of the solidification process was measured using time resolved static light scattering, while the time evolution of the dynamic properties was determined using time resolved dynamic light scattering. By performing further analysis we confirm that solidification of hard sphere colloids is mediated by precursors. Analyzing the dynamic properties we can show that the long time dynamics and thus the shear rigidity of the metastable melt is highly correlated with the number density of solid clusters (precursors) nucleated. In crystallization these objects convert into highly ordered crystals whereas in the case of vitrification this conversion is blocked and the system is (temporarily) locked in the metastable precursor state. From the early stages of solidification one cannot clearly conclude whether the melt will crystallize or vitrify. Furthermore our data suggests that colloidal hard sphere glasses can crystallize via homogeneous nucleation.