Source author record

Michael E. Cates

Michael E. Cates appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2025arXiv

Non-equilibrium pathways between cluster morphologies in active phase separation: necking, rupture and cavitation

We investigate the dynamical pathways of a geometric phase transition in a two-dimensional active lattice gas undergoing motility-induced phase separation. The transition is between metastable morphologies of the liquid cluster: a system-spanning "slab" and a compact "droplet". We generate trajectories of this transition in both directions using forward flux sampling. We find that the droplet-to-slab transition always follows a similar mechanism to its equilibrium counterpart, but the reverse (slab-to-droplet) transition depends on rare non-equilibrium fluctuations. At low Peclet numbers the equilibrium and non-equilibrium pathways compete, while at high Peclet numbers the equilibrium pathway is entirely suppressed, and the only allowed mechanism involves a large vapour bubble. We discuss the implications of these findings for active matter systems more generally.

preprint2022arXiv

Irreversibility and biased ensembles in active matter: Insights from stochastic thermodynamics

Active systems evade the rules of equilibrium thermodynamics by constantly dissipating energy at the level of their microscopic components. This energy flux stems from the conversion of a fuel, present in the environment, into sustained individual motion. It can lead to collective effects without any equilibrium equivalent, such as a phase separation for purely repulsive particles, or a collective motion (flocking) for aligning particles. Some of these effects can be rationalized by using equilibrium tools to recapitulate nonequilibrium transitions. An important challenge is then to delineate systematically to which extent the character of these active transitions is genuinely distinct from equilibrium analogs. We review recent works that use stochastic thermodynamics tools to identify, for active systems, a measure of irreversibility comprising a coarse-grained or informatic entropy production. We describe how this relates to the underlying energy dissipation or thermodynamic entropy production, and how it is influenced by collective behavior. Then, we review the possibility to construct thermodynamic ensembles out-of-equilibrium, where trajectories are biased towards atypical values of nonequilibrium observables. We show that this is a generic route to discovering unexpected phase transitions in active matter systems, which can also inform their design.

preprint2022arXiv

Rectification in a mixture of active and passive particles subject to a ratchet potential

We study by simulation a mixture of active (run-and-tumble) and passive (Brownian) particles with repulsive exclusion interactions in one dimension, subject to a ratchet (smoothed sawtooth) potential. Such a potential is known to rectify active particles at one-body level, creating a net current in the `easy direction'. This is the direction in which one encounters the lower maximum force en route to the top of a potential barrier. The exclusion constraint results in single-file motion, so the mean velocities of active and passive particles are identical; we study the effects of activity level, Brownian diffusivity, particle size, initial sequence of active and passive particles, and active/passive concentration ratio on this mean velocity (i.e., the current per particle). We show that in some parameter regimes the sign of the current is reversed. This happens when the passive particles are at high temperature and so would cross barriers relatively easily, and without rectification, except that they collide with `cold' active ones, which would otherwise be localized near the potential minima. In this case, the reversed current arises because hot passive particles push cold active ones preferentially in the direction with the lower spatial separation between the bottom and top of the barrier. A qualitatively similar mechanism operates in a mixture containing passive particles of two very different temperatures, although there is no quantitative mapping between that case and the systems studied here.

preprint2022arXiv

Stochastic Hydrodynamics of Complex Fluids: Discretisation and Entropy Production

Many complex fluids can be described by continuum hydrodynamic field equations, to which noise must be added in order to capture thermal fluctuations. In almost all cases, the resulting coarse-grained stochastic partial differential equations carry a short-scale cutoff -- which is also reflected in numerical discretisation schemes. We draw together our recent findings concerning the construction of such schemes and the interpretation of their continuum limits, focusing for simplicity on models with a purely diffusive scalar field, such as `Model B' which describes phase separation in binary fluid mixtures. We address the requirement that the steady state entropy production rate (EPR) must vanish for any stochastic hydrodynamic model in thermal equilibrium. Only if this is achieved can the given discretisation scheme be relied upon to correctly calculate the nonvanishing EPR for `active field theories' in which new terms are deliberately added to the fluctuating hydrodynamic equations that break detailed balance. To compute the correct probabilities of forward and time-reversed paths (whose ratio determines the EPR) we must make a careful treatment of so-called `spurious drift' and other closely related terms that depend on the discretisation scheme. We show that such subtleties can arise not only in the temporal discretisation (as is well documented for stochastic ODEs with multiplicative noise) but also from spatial discretisation, even when noise is additive, as most active field theories assume. We then review how such noise can become multiplicative, via off-diagonal couplings to additional fields that encode thermodynamically the underlying chemical processes responsible for activity. In this case the spurious drift terms need careful accounting, not just to evaluate correctly the EPR, but also to numerically implement the Langevin dynamics itself.

preprint2022arXiv

Tracer dynamics in one dimensional gases of active or passive particles

We consider one-dimensional systems comprising either active run-and-tumble particles (RTPs) or passive Brownian random walkers. These particles are either noninteracting or have hardcore exclusions. We study the dynamics of a single tracer particle embedded in such a system - this tracer may be either active or passive, with hardcore exclusion from environmental particles. In an active hardcore environment, both active and passive tracers show long-time subdiffusion: displacements scale as $t^{1/4}$ with a density-dependent prefactor that is independent of tracer type, and differs from the corresponding result for passive-in-passive subdiffusion. In an environment of noninteracting active particles, the passive-in-passive results are recovered at low densities for both active and passive tracers, but transient caging effects slow the tracer motion at higher densities, delaying the onset of any $t^{1/4}$ regime. For an active tracer in a passive environment, we find more complex outcomes, which depend on details of the dynamical discretization scheme. We interpret these results by studying the density distribution of environmental particles around the tracer. In particular, sticking of environment particles to the tracer cause it to move more slowly in noninteracting than in interacting active environments, while the anomalous behaviour of the active-in-passive cases stems from a `snowplough' effect whereby a large pile of diffusive environmental particles accumulates in front of a RTP tracer during a ballistic run.

preprint2021arXiv

Power fluctuations in sheared amorphous materials: A minimal model

The importance of mesoscale fluctuations in flowing amorphous materials is widely accepted, without a clear understanding of their role. We propose a mean-field elastoplastic model that admits both stress and strain-rate fluctuations, and investigate the character of its power distribution under steady shear flow. The model predicts the suppression of negative power fluctuations near the liquid-solid transition; the existence of a fluctuation relation in limiting regimes but its replacement in general by stretched-exponential power-distribution tails; and a crossover between two distinct mechanisms for negative power fluctuations in the liquid and the yielding solid phases. We connect these predictions with recent results from particle-based, numerical micro-rheological experiments.

preprint2020arXiv

Absorbing-state transitions in granular materials close to jamming

We consider a model for driven particulate matter in which absorbing states can be reached both by particle isolation and by particle caging. The model predicts a non-equilibrium phase diagram in which analogues of hydrodynamic and elastic reversibility emerge at low and high volume fractions respectively, partially separated by a diffusive, non-absorbing region. We thus find a single phase boundary that spans the onset of chaos in sheared suspensions to the onset of yielding in jammed packings. This boundary has the properties of a non-equilibrium second order phase transition, leading us to write a Manna-like mean-field description that captures the model predictions. Dependent on contact details, jamming marks either a direct transition between the two absorbing states, or occurs within the diffusive region.

preprint2020arXiv

Collective motion in large deviations of active particles

We analyse collective motion that occurs during rare (large deviation) events in systems of active particles, both numerically and analytically. We discuss the associated dynamical phase transition to collective motion, which occurs when the active work is biased towards larger values, and is associated with alignment of particles' orientations. A finite biasing field is needed to induce spontaneous symmetry breaking, even in large systems. Particle alignment is computed exactly for a system of two particles. For many-particle systems, we analyse the symmetry breaking by an optimal-control representation of the biased dynamics, and we propose a fluctuating hydrodynamic theory that captures the emergence of polar order in the biased state.

preprint2020arXiv

Constitutive model for time-dependent flows of shear-thickening suspensions

We develop a tensorial constitutive model for dense, shear-thickening particle suspensions subjected to time-dependent flow. Our model combines a recently proposed evolution equation for the suspension microstructure in rate-independent materials with ideas developed previously to explain the steady flow of shear-thickening ones, whereby friction proliferates among compressive contacts at large particle stresses. We apply our model to shear reversal, and find good qualitative agreement with particle-level, discrete-element simulations whose results we also present.

preprint2020arXiv

Inference, prediction and optimization of non-pharmaceutical interventions using compartment models: the PyRoss library

PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary structured models formulated stochastically (as master equations) or deterministically (as ODEs) and allows mid-run transitioning from one to the other. By supporting additional compartmental subdivision ad libitum, PyRoss can emulate time-since-infection models and allows medical stages such as hospitalization or quarantine to be modelled and forecast. The PyRoss library enables fitting to epidemiological data, as available, using Bayesian parameter inference, so that competing models can be weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact of idealized NPIs by convolving uncertainties arising from epidemiological data, model choice, parameters, and intrinsic stochasticity. Algorithms to optimize time-dependent NPI scenarios against user-defined cost functions are included. PyRoss's current age-structured compartment framework for well-mixed populations will in future reports be extended to include compartments structured by location, occupation, use of travel networks and other attributes relevant to assessing disease spread and the impact of NPIs. We argue that such compartment models, by allowing social data of arbitrary granularity to be combined with Bayesian parameter estimation for poorly-known disease variables, could enable more powerful and robust prediction than other approaches to detailed epidemic modelling. We invite others to use the PyRoss library for research to address today's COVID-19 crisis, and to plan for future pandemics.

preprint2020arXiv

Non-equilibrium phase separation with reactions: A canonical model and its behaviour

Materials undergoing both phase separation and chemical reactions (defined here as all processes that change particle type or number) form an important class of non-equilibrium systems. Examples range from suspensions of self-propelled bacteria with birth-death dynamics, to bio-molecular condensates, or 'membraneless organelles', within cells. In contrast to their passive counterparts, such systems have conserved and non-conserved dynamics that do not, in general, derive from a shared free energy. This mismatch breaks time-reversal symmetry and leads to new types of dynamical competition that are absent in or near equilibrium. We construct a canonical scalar field theory to describe such systems, with conserved and non-conserved dynamics obeying Model B and Model A respectively (in the Hohenberg-Halperin classification), chosen such that the two free energies involved are incompatible. The resulting minimal model is shown to capture the various phenomenologies reported previously for more complicated models with the same physical ingredients, including microphase separation, limit cycles and droplet splitting. We find a low-dimensional subspace of parameters for which time-reversal symmetry is accidentally recovered, and show that here the dynamics of the order parameter field (but not its conserved current) is exactly the same as an equilibrium system in which microphase separation is caused by long-range attractive interactions.

preprint2020arXiv

Self-propulsion of active droplets without liquid-crystalline order

The swimming of cells, far from any boundary, can arise in the absence of long-range liquid-crystalline order within the cytoplasm, but simple models of this effect are lacking. Here we present a two-dimensional model of droplet self-propulsion involving two scalar fields, representing the cytoplasm and a contractile cortex. An active stress results from coupling between these fields; self-propulsion results when rotational symmetry is spontaneously broken. The swimming speed is predicted, and shown numerically, to vary linearly with the activity parameter and with the droplet area fraction. The model exhibits a Crowley-like instability for an array of active droplets.

preprint2020arXiv

Steady state entropy production rate for scalar Langevin field theories

The entropy production rate (EPR) offers a quantitative measure of time reversal symmetry breaking in non-equilibrium systems. It can be defined either at particle level or at the level of coarse-grained fields such as density; the EPR for the latter quantifies the extent to which these coarse-grained fields behave irreversibly. In this work, we first develop a general method to compute the EPR of scalar Langevin field theories with additive noise. This large class of theories includes active versions of Model A (non-conserved density dynamics) and Model B (conserved) and also models where both types of dynamics are simultaneously present (such as Model AB). Treating the scalar field $ϕ$ (and its time derivative $\dotϕ$) as the sole observable(s), we arrive at an expression for the EPR that is non-negative for every field configuration and is quadratic in the time-antisymmetric component of the dynamics. Our general expression is a function of the quasipotential, which determines the full probability distribution for configurations, and is not generally calculable. To alleviate this difficulty, we present a small-noise expansion of the EPR, which only requires knowledge of the deterministic (mean-field) solution for the scalar field in steady state, which generally is calculable, at least numerically. We demonstrate this calculation for the case of Model AB. We then present a similar EPR calculation for Model AB with the conservative and non-conservative contributions to $\dotϕ= \dotϕ_{\rm A} + \dotϕ_{\rm B}$ viewed as separately observable quantities. The results are qualitatively different, confirming that the field-level EPR depends on the choice of coarse-grained information retained within the dynamical description.

preprint2020arXiv

Stealth entropy production in active field theories near Ising critical points

We address the steady-state entropy production rate (EPR) of active scalar $ϕ^4$ theories, which lack time-reversal symmetry, close to a phase-separation critical point. We consider both nonconserved (Model A) and conserved (Model B) dynamics at Gaussian level, and also address the former at leading order in $ε= 4-d$. In each case, activity is irrelevant in the RG sense: the active model lies in the same (dynamic Ising) universality class as its time-reversible counterpart. Hence one might expect that activity brings no new critical behavior. Here we show instead that the EPR per spacetime correlation volume either remains finite or diverges on approach to criticality. A nontrivial critical scaling for entropy production thus ranks among universal dynamic Ising-class properties.

preprint2020arXiv

Thermodynamic cycles with active matter

Active matter constantly dissipates energy to power the self-propulsion of its microscopic constituents. This opens the door to designing innovative cyclic engines without any equilibrium equivalent. We offer a consistent thermodynamic framework to characterize and optimize the performances of such cycles. Based on a minimal model, we put forward a protocol which extracts work by controlling only the properties of the confining walls at boundaries, and we rationalize the transitions between optimal cycles. We show that the corresponding power and efficiency are generally proportional, so that they reach their maximum values at same cycle time in contrast with thermal cycles, and we provide a generic relation constraining the fluctuations of the power.

preprint2016arXiv

Pattern Formation in Chemically Interacting Active Rotors with Self-Propulsion

We demonstrate that active rotations in chemically signalling particles, such as autochemotactic {\it E. coli} close to walls, create a route for pattern formation based on a nonlinear yet deterministic instability mechanism. For slow rotations, we find a transient persistence of the uniform state, followed by a sudden formation of clusters contingent on locking of the average propulsion direction by chemotaxis. These clusters coarsen, which results in phase separation into a dense and a dilute region. Faster rotations arrest phase separation leading to a global travelling wave of rotors with synchronized roto-translational motion. Our results elucidate the physics resulting from the competition of two generic paradigms in active matter, chemotaxis and active rotations, and show that the latter provides a tool to design programmable self-assembly of active matter, for example to control coarsening.

preprint2016arXiv

Testing the Foundations of Classical Entropy: Colloid Experiments

Defining the entropy of classical particles raises a number of paradoxes and ambiguities, some of which have been known for over a century. Several, such as Gibbs' paradox, involve the fact that classical particles are distinguishable, and in textbooks these are often `resolved' by appeal to the quantum-mechanical indistinguishability of atoms or molecules of the same type. However, questions then remain of how to correctly define the entropy of large poly-atomic particles such as colloids in suspension, of which no two are exactly alike. By performing experiments on such colloids, one can establish that certain definitions of the classical entropy fit the data, while others in the literature do not. Specifically, the experimental facts point firmly to an `informatic' interpretation that dates back to Gibbs: entropy is determined by the number of microstates that we as observers choose to treat as equivalent when we identify a macrostate. This approach, unlike some others, can account for the existence of colloidal crystals, and for the observed abundances of colloidal clusters of different shapes. We also address some lesser-known paradoxes whereby the physics of colloidal assemblies, which ought to be purely classical, seems to involve quantum mechanics directly. The experimental symptoms of such involvement are predicted to be `isotope effects' in which colloids with different inertial masses, but otherwise identical sizes and properties, show different aggregation statistics. These paradoxes are caused by focussing one's attention on some classical degrees while neglecting others; when all are treated equally, all isotope effects are found to vanish.

preprint2016arXiv

Unsteady flow and particle migration in dense, non-Brownian suspensions

We present experimental results on dense corn-starch suspensions as examples of non-Brownian, nearly-hard particles that undergo continuous and discontinuous shear thickening (CST and DST) at intermediate and high densities respectively. Our results offer strong support for recent theories involving a stress-dependent effective contact friction among particles. We show however that in the DST regime, where theory might lead one to expect steady-state shear bands oriented layerwise along the vorticity axis, the real flow is unsteady. To explain this, we argue that steady-state banding is generically ruled out by the requirement that, for hard non-Brownian particles, the solvent pressure and the normal-normal component of the particle stress must balance separately across the interface between bands. (Otherwise there is an unbalanced migration flux.) However, long-lived transient shear-bands remain possible.

preprint2015arXiv

Active Model H: Scalar Active Matter in a Momentum-Conserving Fluid

We present a continuum theory of self-propelled particles, without alignment interactions, in a momentum-conserving solvent. To address phase separation we introduce a scalar concentration field $ϕ$ with advective-diffusive dynamics. Activity creates a contribution $Σ_{ij}=-ζ((\partial_iϕ)(\partial_jϕ)-(\nablaϕ)^{2}δ_{ij}/d)$ to the deviatoric stress, where $ζ$ is odd under time reversal and $d$ is the number of spatial dimensions; this causes an effective interfacial tension contribution that is negative for contractile swimmers. We predict that domain growth then ceases at a length scale where diffusive coarsening is balanced by active stretching of interfaces, and confirm this numerically. Thus the interplay of activity and hydrodynamics is highly nontrivial, even without alignment interactions.

preprint2015arXiv

Activity-induced phase separation and self-assembly in mixtures of active and passive particles

We investigate the phase behavior and kinetics of a monodisperse mixture of active (\textit{i.e.}, self-propelled) and passive isometric Brownian particles through Brownian dynamics simulations and theory. As in a purely active system, motility of the active component triggers phase separation into a dense and a dilute phase; in the dense phase we further find active-passive segregation, with "rafts" of passive particles in a "sea" of active particles. We find that phase separation from an initially disordered mixture can occur with as little as 15 percent of the particles being active. Finally, we show that a system prepared in a suitable fully segregated initial state reproducibly self-assembles an active "corona" which triggers crystallization of the passive core by initiating a compression wave. Our findings are relevant to the experimental pursuit of directed self-assembly using active particles.

preprint2015arXiv

Clustering and pattern formation in chemorepulsive active colloids

We demonstrate that migration away from self-produced chemicals (chemorepulsion) generates a generic route to clustering and pattern formation among self-propelled colloids. The clustering instability can be caused either by anisotropic chemical production, or by a delayed orientational response to changes of the chemical environment. In each case, chemorepulsion creates clusters of a self-limiting area which grows linearly with self-propulsion speed. This agrees with recent observations of dynamic clusters in Janus colloids (albeit not yet known to be chemorepulsive). More generally, our results could inform design principles for the self-assembly of chemorepulsive synthetic swimmers and/or bacteria into nonequilibrium patterns.

preprint2014arXiv

Motility-Induced Phase Separation

Self-propelled particles include both self-phoretic synthetic colloids and various micro-organisms. By continually consuming energy, they bypass the laws of equilibrium thermodynamics. These laws enforce the Boltzmann distribution in thermal equilibrium: the steady state is then independent of kinetic parameters. In contrast, self-propelled particles tend to accumulate where they move more slowly. They may also slow down at high density, for either biochemical or steric reasons. This creates positive feedback which can lead to motility-induced phase separation (MIPS) between dense and dilute fluid phases. At leading order in gradients, a mapping relates variable-speed, self-propelled particles to passive particles with attractions. This deep link to equilibrium phase separation is confirmed by simulations, but generally breaks down at higher order in gradients: new effects, with no equilibrium counterpart, then emerge. We give a selective overview of the fast-developing field of MIPS, focusing on theory and simulation but including a brief speculative survey of its experimental implications.

preprint2014arXiv

Phase behaviour of active Brownian particles: The role of dimensionality

Recently, there has been much interest in activity-induced phase separations in concentrated suspensions of "active Brownian particles" (ABPs), self-propelled spherical particles whose direction of motion relaxes through thermal rotational diffusion. To date, almost all these studies have been restricted to 2 dimensions. In this work we study activity-induced phase separation in 3D and compare the results with previous and new 2D simulations. To this end, we performed state-of-the-art Brownian dynamics simulations of up to 40 million ABPs -- such very large system sizes are unavoidable to evade finite size effects in 3D. Our results confirm the picture established for 2D systems in which an activity-induced phase separation occurs, with strong analogies to equilibrium gas-liquid spinodal decomposition, in spite of the purely non-equilibrium nature of the driving force behind the phase separation. However, we also find important differences between the 2D and 3D cases. Firstly, the shape and position of the phase boundaries is markedly different for the two cases. Secondly, for the 3D coarsening kinetics we find that the domain size grows in time according to the classical diffusive $t^{1/3}$ law, in contrast to the nonstandard subdiffusive exponent observed in 2D.

preprint2014arXiv

Scalar $ϕ^4$ field theory for active-particle phase separation

Recent theories predict phase separation among orientationally disordered active particles whose propulsion speed decreases rapidly enough with density. Coarse-grained models of this process show time-reversal symmetry (detailed balance) to be restored for uniform states, but broken by gradient terms; hence detailed-balance violation is strongly coupled to interfacial phenomena. To explore the subtle generic physics resulting from such coupling we here introduce `Active Model B'. This is a scalar $ϕ^4$ field theory (or phase-field model) that minimally violates detailed balance via a leading-order square-gradient term. We find that this additional term has modest effects on coarsening dynamics, but alters the static phase diagram by creating a jump in (thermodynamic) pressure across flat interfaces. Both results are surprising, since interfacial phenomena are always strongly implicated in coarsening dynamics but are, in detailed-balance systems, irrelevant for phase equilibria.

preprint2013arXiv

A continuum theory of phase separation kinetics for active Brownian particles

Active Brownian particles (ABPs), when subject to purely repulsive interactions, are known to undergo activity-induced phase separation broadly resembling an equilibrium (attraction-induced) gas-liquid coexistence. Here we present an accurate continuum theory for the dynamics of phase-separating ABPs, derived by direct coarse-graining, capturing leading-order density gradient terms alongside an effective bulk free energy. Such gradient terms do not obey detailed balance; yet we find coarsening dynamics closely resembling that of equilibrium phase separation. Our continuum theory is numerically compared to large-scale direct simulations of ABPs and accurately accounts for domain growth kinetics, domain topologies and coexistence densities.

preprint2013arXiv

Rheology of Cubic Blue Phases

We study the behaviour of cubic blue phases under shear flow via lattice Boltzmann simulations. We focus on the two experimentally observed phases, Blue Phase I (BPI) and Blue Phase II (BPII). The disclination network of Blue Phase II continuously breaks and reforms under steady shear, leading to an oscillatory stress response in time. For larger shear rates, the structure breaks up into a Grandjean texture with a cholesteric helix lying along the flow gradient direction. Blue Phase I leads to a very different response. Here, oscillations are only possible for intermediate shear rates -- very slow flow causes a transition of the initially ordered structure into an amorphous network with an apparent yield stress. Larger shear rates lead to another amorphous state with different structure of the defect network. For even larger flow rates the same break-up into a Grandjean texture as for Blue Phase II is observed. At the highest imposed flow rates both cubic blue phases adopt a flow-aligned nematic state. Our results provide the first theoretical investigation of sheared blue phases in large systems, and are relevant to understanding the bulk rheology of these materials.

preprint2012arXiv

Field-Induced Breakup of Emulsion Droplets Stabilized by Colloidal Particles

We simulate the response of a particle-stabilized emulsion droplet in an external force field, such as gravity, acting equally on all $N$ particles. We show that the field strength required for breakup (at fixed initial area fraction) decreases markedly with droplet size, because the forces act cumulatively, not individually, to detach the interfacial particles. The breakup mode involves the collective destabilization of a solidified particle raft occupying the lower part of the droplet, leading to a critical force per particle that scales approximately as $N^{-1/2}$.

preprint2012arXiv

From compact to fractal crystalline clusters in concentrated systems of monodisperse hard spheres

We address the crystallization of monodisperse hard spheres in terms of the properties of finite- size crystalline clusters. By means of large scale event-driven Molecular Dynamics simulations, we study systems at different packing fractions ϕ ranging from weakly supersaturated state points to glassy ones, covering different nucleation regimes. We find that such regimes also result in different properties of the crystalline clusters: compact clusters are formed in the classical-nucleation-theory regime (ϕ \leq 0.54), while a crossover to fractal, ramified clusters is encountered upon increasing packing fraction (ϕ \geq 0.56), where nucleation is more spinodal-like. We draw an analogy between macroscopic crystallization of our clusters and percolation of attractive systems to provide ideas on how the packing fraction influences the final structure of the macroscopic crystals. In our previous work (Phys. Rev. Lett., 106, 215701, 2011), we have demonstrated how crystallization from a glass (at ϕ > 0.58) happens via a gradual (many-step) mechanism: in this paper we show how the mechanism of gradual growth seems to hold also in super-saturated systems just above freezing showing that static properties of clusters are not much affected by dynamics.

preprint2012arXiv

Size limit for particle-stabilized emulsion droplets under gravity

We demonstrate that emulsion droplets stabilized by interfacial particles become unstable beyond a size threshold set by gravity. This holds not only for colloids but supra-colloidal glass beads, using which we directly observe the ejection of particles near the droplet base. The number of particles acting together in these ejection events decreases with time until a stable acorn-like configuration is reached. Stability occurs when the weight of all remaining particles is less than the interfacial binding force of one particle. We also show the importance of the curvature of the droplet surface in promoting particle ejection.

preprint2012arXiv

Thermodynamic Interpretation of Soft Glassy Rheology Models

Mesoscopic models play an important role in our understanding of the deformation and flow of amorphous materials. One such description, based on the Shear Transformation Zone (STZ) theory, has recently been re-formulated within a non-equilibrium thermodynamics framework, and found to be consistent with it. We show here that a similar interpretation can be made for the Soft Glassy Rheology (SGR) model. Conceptually this means that the "noise temperature" x, proposed phenomenologically in the SGR model to control the dynamics of a set of slow mesoscopic degrees of freedom, can consistently be interpreted as their actual thermodynamic temperature. (Because such modes are slow to equilibrate, this generally does not coincide with the temperature of the fast degrees of freedom and/or heat bath.) If one chooses to make this interpretation, the thermodynamic framework significantly constrains extensions of the SGR approach to models in which x is a dynamical variable. We assess in this light some such extensions recently proposed in the context of shear banding.

preprint2010arXiv

Age-dependent transient shear banding in soft glasses

We study numerically the formation of long-lived transient shear bands during shear startup within two models of soft glasses (a simple fluidity model and an adapted `soft glassy rheology' model). The degree and duration of banding depends strongly on the applied shear rate, and on sample age before shearing. In both models the ultimate steady flow state is homogeneous at all shear rates, consistent with the underlying constitutive curve being monotonic. However, particularly in the SGR case, the transient bands can be extremely long lived. The banding instability is neither `purely viscous' nor `purely elastic' in origin, but is closely associated with stress overshoot in startup flow.

preprint2010arXiv

Hawkes process as a model of social interactions: a view on video dynamics

We study by computer simulation the "Hawkes process" that was proposed in a recent paper by Crane and Sornette (Proc. Nat. Acad. Sci. USA 105, 15649 (2008)) as a plausible model for the dynamics of YouTube video viewing numbers. We test the claims made there that robust identification is possible for classes of dynamic response following activity bursts. Our simulated timeseries for the Hawkes process indeed fall into the different categories predicted by Crane and Sornette. However the Hawkes process gives a much narrower spread of decay exponents than the YouTube data, suggesting limits to the universality of the Hawkes-based analysis.

preprint2009arXiv

Colloidal gels assembled via a temporary interfacial scaffold

The liquid-liquid phase separation of a binary solvent can be arrested by colloidal particles trapped at the interface [K. Stratford et al., Science 309, 5744 (2005)]. We show experimentally that the colloidal network so formed can remain stable after fully remixing the liquids, creating a new type of gel in which colloids in a single-phase solvent have locally planar coordination. We argue that this structure is likely maintained by primary-minimum DLVO bonding of our charged colloids, created under strong compression by capillary forces. We present simulation evidence that the combination of a short-ranged attraction with a repulsive barrier can strongly stabilize such locally planar gels.

preprint2006arXiv

Aging in attraction-driven colloidal glasses

Aging in an attraction-driven colloidal glass is studied by computer simulations. The system is equilibrated without attraction and instantaneously ``quenched'', at constant colloid volume fraction, to one of two states beyond the glass transition; one is close to the transition, and the other one deep in the glass. The evolution of structural properties shows that bonds form in the system, increasing the local density, creating density deficits (holes) elsewhere. This process slows down with the time elapsed since the quench. As a consequence of bond formation, there is a slowing down of the dynamics, as measured by the mean squared displacement and the density, bond, and environment correlation functions. The density correlations can be time-rescaled to collapse their long time (structural) decay. The time scale for structural relaxation shows for both quenches a super-linear dependence on waiting time; it grows faster than the bond lifetime, showing the collective origin of the transition. At long waiting times and high attraction strength, we observe {\rem completely} arrested dynamics for more than three decades in time, although individual bonds are not permanent on this time scale. The localization length decreases as the state moves deeper in the glass; the non-ergodicity parameter oscillates in phase with the structure factor. Our main results are obtained for systems with a barrier in the pair potential that inhibits phase separation. However, when this barrier is removed for the case of a deep quench, we find changes in the static structure but almost none in the dynamics. Hence our results for the aging behavior remain relevant to experiments in which the glass transition competes with phase separation.

preprint2002arXiv

Schematic models for dynamic yielding of sheared colloidal glasses

The nonlinear rheological properties of dense suspensions are discussed within simplified models, suggested by a recent first principles approach to the model of Brownian particles in a constant-velocity-gradient solvent flow. Shear thinning of colloidal fluids and dynamical yielding of colloidal glasses arise from a competition between a slowing down of structural relaxation, because of particle interactions, and enhanced decorrelation of fluctuations, caused by the shear advection of density fluctuations. A mode coupling approach is developed to explore the shear-induced suppression of particle caging and the resulting speed-up of the structural relaxation.