Researcher profile

Michael E. Cates

Michael E. Cates contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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.

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.