Showing 1–50 of 111 results for author: Jack, R L

Effects of phase separation on extinction times in population models

Authors: Janik Schüttler, Robert L. Jack, Michael E. Cates

Abstract: We study the effect of phase separating diffusive dynamics on the mean time to extinction in several reaction-diffusion models with slow reactions. We consider a continuum theory similar to model AB, and a simple model where individual particles on two sites undergo on-site reactions and hopping between the sites. In the slow-reaction limit, we project the models' dynamics onto suitable one-dimens… ▽ More We study the effect of phase separating diffusive dynamics on the mean time to extinction in several reaction-diffusion models with slow reactions. We consider a continuum theory similar to model AB, and a simple model where individual particles on two sites undergo on-site reactions and hopping between the sites. In the slow-reaction limit, we project the models' dynamics onto suitable one-dimensional reaction coordinates, which allows derivation of quasi-equilibrium effective free energies. For weak noise, this enables characterisation of the mean time to extinction. This time can be enhanced or suppressed by the addition of phase separation, compared with homogeneous reference cases. We also discuss how Allee effects can be affected by phase separation. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 31 pages

Thermodynamically consistent flocking: From discontinuous to continuous transitions

Authors: Tal Agranov, Robert L. Jack, Michael E. Cates, Étienne Fodor

Abstract: We introduce a family of lattice-gas models of flocking, whose thermodynamically consistent dynamics admits a proper equilibrium limit at vanishing self-propulsion. These models are amenable to an exact coarse-graining which allows us to study their hydrodynamic behavior analytically. We show that the equilibrium limit here belongs to the universality class of Model C, and that it generically exhi… ▽ More We introduce a family of lattice-gas models of flocking, whose thermodynamically consistent dynamics admits a proper equilibrium limit at vanishing self-propulsion. These models are amenable to an exact coarse-graining which allows us to study their hydrodynamic behavior analytically. We show that the equilibrium limit here belongs to the universality class of Model C, and that it generically exhibits tricritical behavior. Self-propulsion has a non-perturbative effect on the phase diagram, yielding novel phase behaviors depending on the type of aligning interactions. For aligning interactions that increase monotonically with the density, the tricritical point diverges to infinite density reproducing the standard scenario of a discontinuous flocking transition accompanied by traveling bands. In contrast, for models where the aligning interaction is non-monotonic in density, the system can exhibit either (the nonequilibrium counterpart of) an azeotropic point, associated with a continuous flocking transition, or a state with counterpropagating bands. △ Less

Submitted 20 May, 2024; v1 submitted 18 January, 2024; originally announced January 2024.

Heterogeneous nucleation in the random field Ising model

Authors: Liheng Yao, Robert L. Jack

Abstract: We investigate the nucleation dynamics of the three-dimensional random field Ising model (RFIM) under an external field. We use umbrella sampling to compute the free-energy cost of a critical nucleus, and use forward flux sampling for the direct estimation of nucleation rates. For moderate to strong disorder, our results indicate that the size of the nucleating cluster is not a good reaction coord… ▽ More We investigate the nucleation dynamics of the three-dimensional random field Ising model (RFIM) under an external field. We use umbrella sampling to compute the free-energy cost of a critical nucleus, and use forward flux sampling for the direct estimation of nucleation rates. For moderate to strong disorder, our results indicate that the size of the nucleating cluster is not a good reaction coordinate, contrary to the pure Ising model. We rectify this problem by introducing a coordinate that also accounts for the location of the nucleus. Using the free energy barrier to predict the nucleation rate, we find reasonable agreement, although deviations become stronger as disorder increases. We attribute this effect to cluster shape fluctuations. We also discuss finite-size effects on the nucleation rate. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 11 pages

Partial and complete wetting of droplets of active Brownian particles

Authors: Francesco Turci, Robert L. Jack, Nigel B. Wilding

Abstract: We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a non… ▽ More We study wetting droplets formed of active Brownian particles in contact with a repulsive potential barrier, in a wedge geometry. Our numerical results demonstrate a transition between partially wet and completely wet states, as a function of the barrier height, analogous to the corresponding surface phase transition in passive fluids. We analyse partially wet configurations characterised by a nonzero contact angle $θしーた$ between the droplet surface and the barrier, including the average density profile and its fluctuations. These findings are compared with two equilibrium systems: a Lennard-Jones fluid and a simple contour model for a liquid-vapour interface. We locate the wetting transition where $\cos(θしーた)=1$, and the neutral state where $\cos(θしーた)=0$. We discuss the implications of these results for possible definitions of surface tensions in active fluids. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Unravelling Metastable Markovian Open Quantum Systems

Authors: Calum A. Brown, Katarzyna Macieszczak, Robert L. Jack

Abstract: We analyse the dynamics of metastable Markovian open quantum systems by unravelling their average dynamics into stochastic trajectories. We use quantum reset processes as examples to illustrate metastable phenomenology, including a simple three-state model whose metastability is of classical type, and a two-qubit model that features a metastable decoherence free subspace. In the three-state model,… ▽ More We analyse the dynamics of metastable Markovian open quantum systems by unravelling their average dynamics into stochastic trajectories. We use quantum reset processes as examples to illustrate metastable phenomenology, including a simple three-state model whose metastability is of classical type, and a two-qubit model that features a metastable decoherence free subspace. In the three-state model, the trajectories exhibit classical metastable phenomenology: fast relaxation into distinct phases and slow transitions between them. This extends the existing correspondence between classical and quantum metastability. It enables the computation of committors for the quantum phases, and the mechanisms of rare transitions between them. For the two-qubit model, the decoherence-free subspace appears in the unravelled trajectories as a slow manifold on which the quantum state has a continuous slow evolution. This provides a classical (non-metastable) analogue of this quantum effect. We discuss the general implications of these results, and we highlight the useful role of quantum reset processes for analysis of quantum trajectories in metastable systems. △ Less

Submitted 5 April, 2024; v1 submitted 27 August, 2023; originally announced August 2023.

Comments: 19 pages, 7 figures

Journal ref: Phys. Rev. A 109, 022244 (2024)

How do Cicadas Emerge Together? Thermophysical Aspects of Their Collective Decision-Making

Authors: Raymond E. Goldstein, Robert L. Jack, Adriana I. Pesci

Abstract: Certain periodical cicadas exhibit life cycles with durations of 13 or 17 years, and it is now generally accepted that such large prime numbers arise evolutionarily to avoid synchrony with predators. Less well explored is the question of how, in the face of intrinsic biological and environmental noise, insects within a brood emerge together in large successive swarms from underground during spring… ▽ More Certain periodical cicadas exhibit life cycles with durations of 13 or 17 years, and it is now generally accepted that such large prime numbers arise evolutionarily to avoid synchrony with predators. Less well explored is the question of how, in the face of intrinsic biological and environmental noise, insects within a brood emerge together in large successive swarms from underground during springtime warming. Here we consider the decision-making process of underground cicada nymphs experiencing random but spatially-correlated thermal microclimates like those in nature. Introducing short-range communication between insects leads to a model of consensus building that maps on to the statistical physics of an Ising model with a quenched, spatially correlated random magnetic field and annealed site dilution, which displays the kinds of collective swarms seen in nature. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 5 pages, 4 figures

Emerging mesoscale flows and chaotic advection in dense active matter

Authors: Yann-Edwin Keta, Juliane Klamser, Robert L. Jack, Ludovic Berthier

Abstract: We study two models of overdamped self-propelled disks in two dimensions, with and without aligning interactions. Active mesoscale flows leading to chaotic advection emerge in both models in the homogeneous dense fluid away from dynamical arrest, forming streams and vortices reminiscent of multiscale flow patterns in turbulence. We show that the characteristics of these flows do not depend on the… ▽ More We study two models of overdamped self-propelled disks in two dimensions, with and without aligning interactions. Active mesoscale flows leading to chaotic advection emerge in both models in the homogeneous dense fluid away from dynamical arrest, forming streams and vortices reminiscent of multiscale flow patterns in turbulence. We show that the characteristics of these flows do not depend on the specific details of the active fluids, and result from the competition between crowding effects and persistent propulsions. Our results suggest that dense active fluids present a type of `active turbulence' distinct from collective flows reported in other types of active systems. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 6 pages, 4 figs

Journal ref: Phys. Rev. Lett. 132, 218301 (2024)

Active-passive mixtures with bulk loading: a minimal active engine in one dimension

Authors: Jean-François Derivaux, Robert L. Jack, Michael E. Cates

Abstract: We study a one-dimensional mixture of active (run-and-tumble) particles and passive (Brownian) particles, with single-file constraint, in a sawtooth potential. The active particles experience a ratchet effect: this generates a current, which can push passive particles against an applied load. The resulting system operates as an active engine. Using numerical simulations, we analyse the efficiency… ▽ More We study a one-dimensional mixture of active (run-and-tumble) particles and passive (Brownian) particles, with single-file constraint, in a sawtooth potential. The active particles experience a ratchet effect: this generates a current, which can push passive particles against an applied load. The resulting system operates as an active engine. Using numerical simulations, we analyse the efficiency of this engine, and we discuss how it can be optimised. Efficient operation occurs when the active particles self-organise into teams, which can push the passive ones against large loads by leveraging collective behaviour. We discuss how the particle arrangement, conserved under the single-file constraint, affects the engine efficiency. We also show that relaxing this constraint still allows the engine to operate effectively. △ Less

Submitted 31 August, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 22 pages, 17 figures

Journal ref: J. Stat. Mech. (2023) 083212

Competition between lanes and transient jammed clusters in driven binary mixtures

Authors: Honghao Yu, Robert L. Jack

Abstract: We consider mixtures of oppositely driven particles, showing that their non-equilibrium steady states form lanes parallel to the drive, which coexist with transient jammed clusters where particles are temporarily immobilised. We analyse the interplay between these two types of non-equilibrium pattern formation, including their implications for macroscopic demixing perpendicular to the drive. Finit… ▽ More We consider mixtures of oppositely driven particles, showing that their non-equilibrium steady states form lanes parallel to the drive, which coexist with transient jammed clusters where particles are temporarily immobilised. We analyse the interplay between these two types of non-equilibrium pattern formation, including their implications for macroscopic demixing perpendicular to the drive. Finite-size scaling analysis indicates that there is no critical driving force associated with demixing, which appears as a crossover in finite systems. We attribute this effect to the disruption of long-ranged order by the transient jammed clusters. △ Less

Submitted 19 February, 2024; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: authors' accepted version

Journal ref: Phys. Rev. E 109, 024123 (2024)

Necking and failure of a colloidal gel arm: signatures of yielding on different length scales

Authors: Kristian Thijssen, T. B. Liverpool, C. Patrick Royall, Robert L. Jack

Abstract: Colloidal gels consist of percolating networks of interconnected arms. Their mechanical properties depend on the individual arms, and on their collective behaviour. We use numerical simulations to pull on a single arm, built from a model colloidal gel-former with short-ranged attractive interactions. Under elongation, the arm breaks by a necking instability. We analyse this behaviour at three diff… ▽ More Colloidal gels consist of percolating networks of interconnected arms. Their mechanical properties depend on the individual arms, and on their collective behaviour. We use numerical simulations to pull on a single arm, built from a model colloidal gel-former with short-ranged attractive interactions. Under elongation, the arm breaks by a necking instability. We analyse this behaviour at three different length scales: a rheological continuum model of the whole arm; a microscopic analysis of the particle structure and dynamics; and the local stress tensor. Combining these different measurements gives a coherent picture of the necking and failure: the neck is characterised by plastic flow that occurs for stresses close to the arm's yield stress. The arm has an amorphous local structure and has large residual stresses from its initialisation. We find that neck formation is associated with increased plastic flow, a reduction in the stability of the local structure, and a reduction in the residual stresses; this indicates that {the} system loses its solid character and starts to behave more like a viscous fluid. We discuss the implications of these results for the modelling of gel dynamics. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 15 pages

Intermittent relaxation and avalanches in extremely persistent active matter

Authors: Yann-Edwin Keta, Rituparno Mandal, Peter Sollich, Robert L. Jack, Ludovic Berthier

Abstract: We use numerical simulations to study the dynamics of dense assemblies of self-propelled particles in the limit of extremely large, but finite, persistence times. In this limit, the system evolves intermittently between mechanical equilibria where active forces balance interparticle interactions. We develop an efficient numerical strategy allowing us to resolve the statistical properties of elasti… ▽ More We use numerical simulations to study the dynamics of dense assemblies of self-propelled particles in the limit of extremely large, but finite, persistence times. In this limit, the system evolves intermittently between mechanical equilibria where active forces balance interparticle interactions. We develop an efficient numerical strategy allowing us to resolve the statistical properties of elastic and plastic relaxation events caused by activity-driven fluctuations. The system relaxes via a succession of scale-free elastic events and broadly distributed plastic events that both depend on the system size. Correlations between plastic events lead to emergent dynamic facilitation and heterogeneous relaxation dynamics. Our results show that dynamical behaviour in extremely persistent active systems is qualitatively similar to that of sheared amorphous solids, yet with some important differences. △ Less

Submitted 5 June, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Authors' accepted version for publication in Soft Matter

Journal ref: Soft Matter 19, 3871 (2023)

Tricritical behavior in dynamical phase transitions

Authors: Tal Agranov, Michael E. Cates, Robert L. Jack

Abstract: We identify a new scenario for dynamical phase transitions associated with time-integrated observables occurring in diffusive systems described by the macroscopic fluctuation theory. It is characterized by the pairwise meeting of first- and second-order bias-induced phase transition curves at two tricritical points. We formulate a simple, general criterion for its appearance and derive an exact La… ▽ More We identify a new scenario for dynamical phase transitions associated with time-integrated observables occurring in diffusive systems described by the macroscopic fluctuation theory. It is characterized by the pairwise meeting of first- and second-order bias-induced phase transition curves at two tricritical points. We formulate a simple, general criterion for its appearance and derive an exact Landau theory for the tricritical behavior. The scenario is demonstrated in three examples: the simple symmetric exclusion process biased by an activity-related structural observable; the Katz-Lebowitz-Spohn lattice gas model biased by its current; and in an active lattice gas biased by its entropy production. △ Less

Submitted 24 May, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: 21 pages, authors' accepted version

Dimensionality reduction of local structure in glassy binary mixtures

Authors: Daniele Coslovich, Robert L. Jack, Joris Paret

Abstract: We consider unsupervised learning methods for characterizing the disordered microscopic structure of supercooled liquids and glasses. Specifically, we perform dimensionality reduction of smooth structural descriptors that describe radial and bond-orientational correlations, and assess the ability of the method to grasp the essential structural features of glassy binary mixtures. In several cases,… ▽ More We consider unsupervised learning methods for characterizing the disordered microscopic structure of supercooled liquids and glasses. Specifically, we perform dimensionality reduction of smooth structural descriptors that describe radial and bond-orientational correlations, and assess the ability of the method to grasp the essential structural features of glassy binary mixtures. In several cases, a few collective variables account for the bulk of the structural fluctuations within the first coordination shell and also display a clear connection with the fluctuations of particle mobility. Fine-grained descriptors that characterize the radial dependence of bond-orientational order better capture the structural fluctuations relevant for particle mobility, but are also more difficult to parametrize and to interpret. We also find that principal component analysis of bond-orientational order parameters provides identical results to neural network autoencoders, while having the advantage of being easily interpretable. Overall, our results indicate that glassy binary mixtures have a broad spectrum of structural features. In the temperature range we investigate, some mixtures display well-defined locally favored structures, which are reflected in bimodal distributions of the structural variables identified by dimensionality reduction. △ Less

Submitted 23 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 22 pages, 16 figures. Data and workflow to reproduce the findings of this work are available at https://doi.org/10.5281/zenodo.7108316

Journal ref: J. Chem. Phys. 157, 204503 (2022)

Thermal vestiges of avalanches in the driven random field Ising model

Authors: Liheng Yao, Robert L. Jack

Abstract: We investigate the non-equilibrium behavior of the $3d$ random field Ising model at finite temperature, as an external field is increased through its coercive field. We show by numerical simulations that the phenomenology of avalanches -- which are sharply defined only at zero temperature -- also persists over a significant range of finite temperatures. We analyse the main differences between the… ▽ More We investigate the non-equilibrium behavior of the $3d$ random field Ising model at finite temperature, as an external field is increased through its coercive field. We show by numerical simulations that the phenomenology of avalanches -- which are sharply defined only at zero temperature -- also persists over a significant range of finite temperatures. We analyse the main differences between the thermal and zero-temperature systems, including an excess of small avalanches in the thermal case, whose behaviour is consistent with activated dynamical scaling. We also investigate the extent to which individual avalanches at finite temperature can be traced back to parent avalanches in the athermal system. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 16 pages

Entropy production and its large deviations in an active lattice gas

Authors: Tal Agranov, Michael E. Cates, Robert L. Jack

Abstract: Active systems are characterized by a continuous production of entropy at steady state. We study the statistics of entropy production within a lattice-based model of interacting active particles that is capable of motility-induced phase separation. Exploiting a recent formulation of the exact fluctuating hydrodynamics for this model, we provide analytical results for its entropy production statist… ▽ More Active systems are characterized by a continuous production of entropy at steady state. We study the statistics of entropy production within a lattice-based model of interacting active particles that is capable of motility-induced phase separation. Exploiting a recent formulation of the exact fluctuating hydrodynamics for this model, we provide analytical results for its entropy production statistics in both typical and atypical (biased) regimes. This complements previous studies of the large deviation statistics of entropy production in off-lattice active particle models that could only be addressed numerically. Our analysis uncovers an unexpectedly intricate phase diagram, with five different phases arising (under bias) within the parameter regime where the unbiased system is in its homogeneous state. Notably, we find the concurrence of first order and second order nonequilibrium phase transition curves at a bias-induced tricritical point, a feature not yet reported in previous studies of active systems. △ Less

Submitted 13 September, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Convergence of random-weight sequential Monte Carlo methods

Authors: Paul B. Rohrbach, Robert L. Jack

Abstract: We investigate the properties of a sequential Monte Carlo method where the particle weight that appears in the algorithm is estimated by a positive, unbiased estimator. We present broadly-applicable convergence results, including a central limit theorem, and we discuss their relevance for applications in statistical physics. Using these results, we show that the resampling step reduces the impact… ▽ More We investigate the properties of a sequential Monte Carlo method where the particle weight that appears in the algorithm is estimated by a positive, unbiased estimator. We present broadly-applicable convergence results, including a central limit theorem, and we discuss their relevance for applications in statistical physics. Using these results, we show that the resampling step reduces the impact of the randomness of the weights on the asymptotic variance of the estimator. In addition, we explore the limits of convergence of the sequential Monte Carlo method, with a focus on almost sure convergence. We construct an example algorithm where we can prove convergence in probability, but which does not converge almost surely, even in the non-random-weight case. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: 31 pages

MSC Class: 65C05; 62L10; 60F05; 60F15

Role of initial conditions in $1D$ diffusive systems: compressibility, hyperuniformity and long-term memory

Authors: Tirthankar Banerjee, Robert L. Jack, Michael E. Cates

Abstract: We analyse the long-lasting effects of initial conditions on fluctuations in one-dimensional diffusive systems. We consider both the fluctuations of current for non-interacting diffusive particles starting from a step-like initial density profile, and the mean-square displacement of tracers in homogeneous systems with single-file diffusion. For these two cases, we show analytically (via the propag… ▽ More We analyse the long-lasting effects of initial conditions on fluctuations in one-dimensional diffusive systems. We consider both the fluctuations of current for non-interacting diffusive particles starting from a step-like initial density profile, and the mean-square displacement of tracers in homogeneous systems with single-file diffusion. For these two cases, we show analytically (via the propagator and Macroscopic Fluctuation Theory, respectively) that the long-term memory of initial conditions is mediated by a single static quantity: a generalized compressibility that quantifies the density fluctuations of the initial state. We thereby identify a universality class of hyperuniform initial states whose dynamical variances coincide with the `quenched' cases studied previously; we also describe a continuous family of other classes among which equilibrated (or `annealed') initial conditions are but one family member. We verify our predictions through extensive Monte Carlo simulations. △ Less

Submitted 6 November, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: accepted in Physical Review E (Letter)

Journal ref: Phys. Rev. E 106, L062101 (2022)

Multilevel simulation of hard-sphere mixtures

Authors: Paul B. Rohrbach, Hideki Kobayashi, Robert Scheichl, Nigel B. Wilding, Robert L. Jack

Abstract: We present a multilevel Monte Carlo simulation method for analysing multi-scale physical systems via a hierarchy of coarse-grained representations, to obtain numerically-exact results, at the most detailed level. We apply the method to a mixture of size-asymmetric hard spheres, in the grand canonical ensemble. A three-level version of the method is compared with a previously-studied two-level vers… ▽ More We present a multilevel Monte Carlo simulation method for analysing multi-scale physical systems via a hierarchy of coarse-grained representations, to obtain numerically-exact results, at the most detailed level. We apply the method to a mixture of size-asymmetric hard spheres, in the grand canonical ensemble. A three-level version of the method is compared with a previously-studied two-level version. The extra level interpolates between the full mixture and a coarse-grained description where only the large particles are present -- this is achieved by restricting the small particles to regions close to the large ones. The three-level method improves the performance of the estimator, at fixed computational cost. We analyse the asymptotic variance of the estimator, and discuss the mechanisms for the improved performance. △ Less

Submitted 26 August, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: 24 pages, authors' accepted version

Mechanical analysis of a dynamical phase transition for particles in a channel

Authors: Jakub Dolezal, Robert L. Jack

Abstract: We analyse biased ensembles of trajectories for a two-dimensional system of particles, evolving by Langevin dynamics in a channel geometry. This bias controls the degree of particle clustering. On biasing to large clustering, we observe a dynamical phase transition where the particles break symmetry and accumulate at one of the walls. We analyse the mechanical properties of this symmetry-broken st… ▽ More We analyse biased ensembles of trajectories for a two-dimensional system of particles, evolving by Langevin dynamics in a channel geometry. This bias controls the degree of particle clustering. On biasing to large clustering, we observe a dynamical phase transition where the particles break symmetry and accumulate at one of the walls. We analyse the mechanical properties of this symmetry-broken state using the Irving-Kirkwood stress tensor. The biased ensemble is characterised by body forces which originate in random thermal noises, but have finite averages in the presence of the bias. We discuss the connection of these forces to Doob's transform and optimal control theory. △ Less

Submitted 2 September, 2022; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 10 pages, authors' accepted version

Journal ref: Phys. Rev. Research 4, 033134 (2022)

Perpendicular and Parallel Phase Separation in Two Species Driven Diffusive Lattice Gases

Authors: Honghao Yu, Kristian Thijssen, Robert L. Jack

Abstract: We study three different lattice models in which two species of diffusing particles are driven in opposite directions by an electric field. We focus on dynamical phase transitions that involve phase separation into domains that may be parallel or perpendicular to a driving field. In all cases, the perpendicular state appears for weak driving, consistent with previous work. For strong driving, we i… ▽ More We study three different lattice models in which two species of diffusing particles are driven in opposite directions by an electric field. We focus on dynamical phase transitions that involve phase separation into domains that may be parallel or perpendicular to a driving field. In all cases, the perpendicular state appears for weak driving, consistent with previous work. For strong driving, we introduce two models that support the parallel state. In one model, this state occurs because of the inclusion of dynamical rules that enhance lateral diffusion during collisions; in the other, it is a result of a nearest-neighbour attractive/repulsive interaction between particles of the same/opposite species. We discuss the connections between these results and the behaviour found in off-lattice systems, including laning and freezing by heating. △ Less

Submitted 4 August, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: authors' accepted version

Journal ref: Phys. Rev. E, 106, 024129 (2022)

Direct Imaging of Contacts and Forces in Colloidal Gels

Authors: Jun Dong, Francesco Turci, Robert L. Jack, Malcolm A. Faers, C. Patrick Royall

Abstract: Colloidal dispersions are prized as model systems to understand basic properties of materials, and are central to a wide range of industries from cosmetics to foods to agrichemicals. Among the key developments in using colloids to address challenges in condensed matter is to resolve the particle coordinates in 3D, allowing a level of analysis usually only possible in computer simulation. However i… ▽ More Colloidal dispersions are prized as model systems to understand basic properties of materials, and are central to a wide range of industries from cosmetics to foods to agrichemicals. Among the key developments in using colloids to address challenges in condensed matter is to resolve the particle coordinates in 3D, allowing a level of analysis usually only possible in computer simulation. However in amorphous materials, relating mechanical properties to microscopic structure remains problematic. This makes it rather hard to understand, for example, mechanical failure. Here we address this challenge by studying the contacts and the forces between particles, as well as their positions. To do so, we use a colloidal model system (an emulsion) in which the interparticle forces and local stress can be linked to the microscopic structure. We demonstrate the potential of our method to reveal insights into the failure mechanisms of soft amorphous solids by determining local stress in a colloidal gel. In particular, we identify `force chains' of load-bearing droplets, and local stress anisotropy, and investigate their connection with locally rigid packings of the droplets. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 12 pages

Macroscopic behaviour in a two-species exclusion process via the method of matched asymptotics

Authors: James Mason, Robert L Jack, Maria Bruna

Abstract: We consider a two-species simple exclusion process on a periodic lattice. We use the method of matched asymptotics to derive evolution equations for the two population densities in the dilute regime, namely a cross-diffusion system of partial differential equations for the two species densities. First, our result captures non-trivial interaction terms neglected in the mean-field approach, includin… ▽ More We consider a two-species simple exclusion process on a periodic lattice. We use the method of matched asymptotics to derive evolution equations for the two population densities in the dilute regime, namely a cross-diffusion system of partial differential equations for the two species densities. First, our result captures non-trivial interaction terms neglected in the mean-field approach, including a non-diagonal mobility matrix with explicit density dependence. Second, it generalises the rigorous hydrodynamic limit of Quastel [Commun. Pure Appl. Math. 45(6), 623--679 (1992)], valid for species with equal jump rates and given in terms of a non-explicit self-diffusion coefficient, to the case of unequal rates in the dilute regime. In the equal-rates case, by combining matched asymptotic approximations in the low- and high-density limits, we obtain a cubic polynomial approximation of the self-diffusion coefficient that is numerically accurate for all densities. This cubic approximation agrees extremely well with numerical simulations. It also coincides with the Taylor expansion up to the second-order in the density of the self-diffusion coefficient obtained using a rigorous recursive method. △ Less

Submitted 3 January, 2023; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 49 pages, 7 figures

MSC Class: 82C22; 41A60; 60J74

Disordered collective motion in dense assemblies of persistent particles

Authors: Yann-Edwin Keta, Robert L. Jack, Ludovic Berthier

Abstract: We explore the emergence of nonequilibrium collective motion in disordered non-thermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, char… ▽ More We explore the emergence of nonequilibrium collective motion in disordered non-thermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field △ Less

Submitted 23 June, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: Authors' accepted version for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 129, 048002 (2022)

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

Authors: Jean-François Derivaux, Robert L. Jack, Michael E. Cates

Abstract: 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… ▽ More 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. △ Less

Submitted 1 April, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: 19 pages, 11 figures

Journal ref: J. Stat. Mech. 2022, 043203 (2022)

Tracer dynamics in one dimensional gases of active or passive particles

Authors: Tirthankar Banerjee, Robert L. Jack, Michael E. Cates

Abstract: 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 e… ▽ More 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. △ Less

Submitted 9 February, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: 32 pages, 26 figures

Journal ref: J. Stat. Mech. (2022) 013209

Dynamical phase transition in the activity-biased fully-connected random field Ising model: connection with glass-forming systems

Authors: Jules Guioth, Robert L. Jack

Abstract: We analyse biased ensembles of trajectories for the random-field Ising model on a fully-connected lattice, which is described exactly by mean-field theory. By coupling the activity of the system to a dynamical biasing field, we find a range of dynamical phase transitions, including spontaneous symmetry breaking into ordered states. For weak bias, the phase behaviour is controlled by extrema of the… ▽ More We analyse biased ensembles of trajectories for the random-field Ising model on a fully-connected lattice, which is described exactly by mean-field theory. By coupling the activity of the system to a dynamical biasing field, we find a range of dynamical phase transitions, including spontaneous symmetry breaking into ordered states. For weak bias, the phase behaviour is controlled by extrema of the free energy, which may be local minima or saddle points. For large bias, the system tends to states of extremal activity, which may differ strongly from free energy minima. We discuss connections of these results to random first-order transition theory of glasses, which motivates an extension of the analysis to random-field Ising models where the dynamical activity is not symmetric under magnetisation reversal. △ Less

Submitted 26 October, 2021; v1 submitted 30 July, 2021; originally announced August 2021.

Comments: 26 pages, 10 figures

Journal ref: J. Stat. Mech. (2021) 103206

Cycle counts and affinities in stochastic models of non-equilibrium systems

Authors: Patrick Pietzonka, Jules Guioth, Robert L. Jack

Abstract: For non-equilibrium systems described by finite Markov processes, we consider the number of times that a system traverses a cyclic sequence of states (a cycle). The joint distribution of the number of forward and backward instances of any given cycle is described by universal formulae which depend on the cycle affinity, but are otherwise independent of system details. We discuss the similarities a… ▽ More For non-equilibrium systems described by finite Markov processes, we consider the number of times that a system traverses a cyclic sequence of states (a cycle). The joint distribution of the number of forward and backward instances of any given cycle is described by universal formulae which depend on the cycle affinity, but are otherwise independent of system details. We discuss the similarities and differences of this result to fluctuation theorems, and generalize the result to families of cycles, relevant under coarse-graining. Finally, we describe the application of large deviation theory to this cycle counting problem. △ Less

Submitted 9 January, 2022; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: authors' accepted version

Critical point for de-mixing of binary hard spheres

Authors: Hideki Kobayashi, Paul B. Rohrbach, Robert Scheichl, Nigel B. Wilding, Robert L. Jack

Abstract: We use a two-level simulation method to analyse the critical point associated with demixing of binary hard sphere mixtures. The method exploits an accurate coarse-grained model with two-body and three-body effective interactions. Using this model within the two-level methodology allows computation of properties of the full (fine-grained) mixture. The critical point is located by computing the prob… ▽ More We use a two-level simulation method to analyse the critical point associated with demixing of binary hard sphere mixtures. The method exploits an accurate coarse-grained model with two-body and three-body effective interactions. Using this model within the two-level methodology allows computation of properties of the full (fine-grained) mixture. The critical point is located by computing the probability distribution for the number of large particles in the grand canonical ensemble, and matching to the universal form for the $3d$ Ising universality class. The results have a strong and unexpected dependence on the size ratio between large and small particles, which is related to three-body effective interactions, and the geometry of the underlying hard sphere packings. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 11 pages

Journal ref: Phys. Rev. E 104, 044603 (2021)

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

Authors: Étienne Fodor, Robert L. Jack, Michael E. Cates

Abstract: 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… ▽ More 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. △ Less

Submitted 12 March, 2022; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 23 pages, 3 figures

Journal ref: Annu. Rev. Condens. Matter Phys 13, 215 (2022)

Large deviations at level 2.5 for Markovian open quantum systems: quantum jumps and quantum state diffusion

Authors: Federico Carollo, Juan P. Garrahan, Robert L. Jack

Abstract: We consider quantum stochastic processes and discuss a level 2.5 large deviation formalism providing an explicit and complete characterisation of fluctuations of time-averaged quantities, in the large-time limit. We analyse two classes of quantum stochastic dynamics, within this framework. The first class consists of the quantum jump trajectories related to photon detection; the second is quantum… ▽ More We consider quantum stochastic processes and discuss a level 2.5 large deviation formalism providing an explicit and complete characterisation of fluctuations of time-averaged quantities, in the large-time limit. We analyse two classes of quantum stochastic dynamics, within this framework. The first class consists of the quantum jump trajectories related to photon detection; the second is quantum state diffusion related to homodyne detection. For both processes, we present the level 2.5 functional starting from the corresponding quantum stochastic Schrödinger equation and we discuss connections of these functionals to optimal control theory. △ Less

Submitted 7 July, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: 45 pages, 4 figures

Journal ref: J. Stat. Phys. 184, 13 (2021)

Long-ranged correlations in large deviations of local clustering

Authors: Jakub Dolezal, Robert L. Jack

Abstract: In systems of diffusing particles, we investigate large deviations of a time-averaged measure of clustering around one particle. We focus on biased ensembles of trajectories, which realise large-deviation events. The bias acts on a single particle, but elicits a response that spans the whole system. We analyse this effect through the lens of Macroscopic Fluctuation Theory, focussing on the coupl… ▽ More In systems of diffusing particles, we investigate large deviations of a time-averaged measure of clustering around one particle. We focus on biased ensembles of trajectories, which realise large-deviation events. The bias acts on a single particle, but elicits a response that spans the whole system. We analyse this effect through the lens of Macroscopic Fluctuation Theory, focussing on the coupling of the bias to hydrodynamic modes. This explains that the dynamical free energy has non-trivial scaling relationships with the system size, in 1 and 2 spatial dimensions. We show that the long-ranged response to a bias on one particle also has consequences when biasing two particles. △ Less

Submitted 24 May, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 14 pages, 11 figures

Journal ref: Phys. Rev. E 103, 052132 (2021)

Efficient Bayesian inference of fully stochastic epidemiological models with applications to COVID-19

Authors: Yuting I. Li, Günther Turk, Paul B. Rohrbach, Patrick Pietzonka, Julian Kappler, Rajesh Singh, Jakub Dolezal, Timothy Ekeh, Lukas Kikuchi, Joseph D. Peterson, Austen Bolitho, Hideki Kobayashi, Michael E. Cates, R. Adhikari, Robert L. Jack

Abstract: Epidemiological forecasts are beset by uncertainties about the underlying epidemiological processes, and the surveillance process through which data are acquired. We present a Bayesian inference methodology that quantifies these uncertainties, for epidemics that are modelled by (possibly) non-stationary, continuous-time, Markov population processes. The efficiency of the method derives from a func… ▽ More Epidemiological forecasts are beset by uncertainties about the underlying epidemiological processes, and the surveillance process through which data are acquired. We present a Bayesian inference methodology that quantifies these uncertainties, for epidemics that are modelled by (possibly) non-stationary, continuous-time, Markov population processes. The efficiency of the method derives from a functional central limit theorem approximation of the likelihood, valid for large populations. We demonstrate the methodology by analysing the early stages of the COVID-19 pandemic in the UK, based on age-structured data for the number of deaths. This includes maximum a posteriori estimates, MCMC sampling of the posterior, computation of the model evidence, and the determination of parameter sensitivities via the Fisher information matrix. Our methodology is implemented in PyRoss, an open-source platform for analysis of epidemiological compartment models. △ Less

Submitted 13 August, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: Revised and accepted version, 24 pages

Journal ref: R. Soc. Open Sci. 8, 211065 (2021)

Collective motion in large deviations of active particles

Authors: Yann-Edwin Keta, Étienne Fodor, Frédéric van Wijland, Michael E. Cates, Robert L. Jack

Abstract: 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 spontaneo… ▽ More 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. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 28 pages

Journal ref: Phys. Rev. E 103, 022603 (2021)

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

Authors: R. Adhikari, Austen Bolitho, Fernando Caballero, Michael E. Cates, Jakub Dolezal, Timothy Ekeh, Jules Guioth, Robert L. Jack, Julian Kappler, Lukas Kikuchi, Hideki Kobayashi, Yuting I. Li, Joseph D. Peterson, Patrick Pietzonka, Benjamin Remez, Paul B. Rohrbach, Rajesh Singh, Günther Turk

Abstract: 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 s… ▽ More 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. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: Code and updates at https://github.com/rajeshrinet/pyross 75 pages, 14 figures, and 1 table

Giant leaps and long excursions: fluctuation mechanisms in systems with long-range memory

Authors: Robert L. Jack, Rosemary J. Harris

Abstract: We analyse large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_t of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behaviour, which differ from the Markovian case. In particular, the memo… ▽ More We analyse large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_t of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behaviour, which differ from the Markovian case. In particular, the memory can lead to large deviation principles with reduced speeds, and to non-analytic rate functions. We then explain how the mechanisms operating in these two models are generic for memory-dependent dynamics and show other examples including a non-Markovian symmetric exclusion process. △ Less

Submitted 16 June, 2020; v1 submitted 7 March, 2020; originally announced March 2020.

Comments: longer version (16 pages), with more detailed discussion

Journal ref: Phys. Rev. E 102, 012154 (2020)

Assessing the structural heterogeneity of supercooled liquids through community inference

Authors: Joris Paret, Robert L. Jack, Daniele Coslovich

Abstract: We present an information-theoretic approach inspired by distributional clustering to assess the structural heterogeneity of particulate systems. Our method identifies communities of particles that share a similar local structure by harvesting the information hidden in the spatial variation of two- or three-body static correlations. This corresponds to an unsupervised machine learning approach tha… ▽ More We present an information-theoretic approach inspired by distributional clustering to assess the structural heterogeneity of particulate systems. Our method identifies communities of particles that share a similar local structure by harvesting the information hidden in the spatial variation of two- or three-body static correlations. This corresponds to an unsupervised machine learning approach that infers communities solely from the particle positions and their species. We apply this method to three models of supercooled liquids and find that it detects subtle forms of local order, as demonstrated by a comparison with the statistics of Voronoi cells. Finally, we analyze the time-dependent correlation between structural communities and particle mobility and show that our method captures relevant information about glassy dynamics. △ Less

Submitted 17 April, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: 20 pages, 10 figures. Data relevant to this work are available at https://doi.org/10.5281/zenodo.3653942

Journal ref: J. Chem. Phys. 152, 144502 (2020)

Dynamical phase coexistence in the Fredrickson-Andersen model

Authors: Robert L. Jack, Takahiro Nemoto, Vivien Lecomte

Abstract: We analyse a first-order dynamical phase transition that takes place in the Fredrickson--Andersen (FA) model. We construct a two-dimensional spin system whose thermodynamic properties reproduce the dynamical large deviations of the FA model and we analyse this system numerically, comparing our results with finite-size scaling theory. This allows us to rationalise recent results for the FA model, i… ▽ More We analyse a first-order dynamical phase transition that takes place in the Fredrickson--Andersen (FA) model. We construct a two-dimensional spin system whose thermodynamic properties reproduce the dynamical large deviations of the FA model and we analyse this system numerically, comparing our results with finite-size scaling theory. This allows us to rationalise recent results for the FA model, including the exponential divergence of its susceptibility at phase coexistence. We also discuss a simple interfacial model that reproduces quantitatively the behaviour of the FA model at coexistence. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 32 pages

Journal ref: J. Stat. Mech. 053204 (2020)

Ergodicity and large deviations in physical systems with stochastic dynamics

Authors: Robert L. Jack

Abstract: In ergodic physical systems, time-averaged quantities converge (for large times) to their ensemble-averaged values. Large deviation theory describes rare events where these time averages differ significantly from the corresponding ensemble averages. It allows estimation of the probabilities of these events, and their mechanisms. This theory has been applied to a range of physical systems, where it… ▽ More In ergodic physical systems, time-averaged quantities converge (for large times) to their ensemble-averaged values. Large deviation theory describes rare events where these time averages differ significantly from the corresponding ensemble averages. It allows estimation of the probabilities of these events, and their mechanisms. This theory has been applied to a range of physical systems, where it has yielded new insights into entropy production, current fluctuations, metastability, transport processes, and glassy behaviour. We review some of these developments, identifying general principles. We discuss a selection of dynamical phase transitions, and we highlight some connections between large-deviation theory and optimal control theory. △ Less

Submitted 11 March, 2020; v1 submitted 22 October, 2019; originally announced October 2019.

Comments: accepted version, small additions and modifications. 23 pages, "colloquium" article

Correction of coarse-graining errors by a two-level method: application to the Asakura-Oosawa model

Authors: Hideki Kobayashi, Paul B. Rohrbach, Robert Scheichl, Nigel B. Wilding, Robert L. Jack

Abstract: We present a method that exploits self-consistent simulation of coarse-grained and fine-grained models, in order to analyse properties of physical systems. The method uses the coarse-grained model to obtain a first estimate of the quantity of interest, before computing a correction by analysing properties of the fine system. We illustrate the method by applying it to the Asakura-Oosawa (AO) model… ▽ More We present a method that exploits self-consistent simulation of coarse-grained and fine-grained models, in order to analyse properties of physical systems. The method uses the coarse-grained model to obtain a first estimate of the quantity of interest, before computing a correction by analysing properties of the fine system. We illustrate the method by applying it to the Asakura-Oosawa (AO) model of colloid-polymer mixtures. We show that the liquid-vapour critical point in that system is affected by three-body interactions which are neglected in the corresponding coarse-grained model. We analyse the size of this effect and the nature of the three-body interactions. We also analyse the accuracy of the method, as a function of the associated computational effort. △ Less

Submitted 18 July, 2019; originally announced July 2019.

Comments: 17 pages

Large deviations and optimal control forces for hard particles in one dimension

Authors: Jakub Dolezal, Robert L. Jack

Abstract: We analyse large deviations of the dynamical activity in one-dimensional systems of diffusing hard particles. Using an optimal-control representation of the large-deviation problem, we analyse effective interaction forces which can be added to the system, to aid sampling of biased ensembles of trajectories. We find several distinct regimes, as a function of the activity and the system size: we pre… ▽ More We analyse large deviations of the dynamical activity in one-dimensional systems of diffusing hard particles. Using an optimal-control representation of the large-deviation problem, we analyse effective interaction forces which can be added to the system, to aid sampling of biased ensembles of trajectories. We find several distinct regimes, as a function of the activity and the system size: we present approximate analytical calculations that characterise the effective interactions in several of these regimes. For high activity the system is hyperuniform and the interactions are long-ranged and repulsive. For low activity, there is a near-equilibrium regime described by macroscopic fluctuation theory, characterised by long-ranged attractive forces. There is also a far-from-equilibrium regime in which one of the interparticle gaps becomes macroscopic and the interactions depend strongly on the size of this gap. We discuss the extent to which transition path sampling of these ensembles is improved by adding suitable control forces. △ Less

Submitted 27 September, 2019; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: 37 pages

Journal ref: J. Stat. Mech. (2019) 123208

Large deviations in models of growing clusters with symmetry-breaking transitions

Authors: Robert L. Jack

Abstract: We analyse large deviations of the magnetisation in two models of growing clusters. The models have symmetry-breaking transitions, so the typical magnetisation of a growing cluster may be either positive or negative, with equal probability. For large clusters, the magnetisation obeys a large deviation principle. We show that the corresponding rate function is zero for values of the magnetisation t… ▽ More We analyse large deviations of the magnetisation in two models of growing clusters. The models have symmetry-breaking transitions, so the typical magnetisation of a growing cluster may be either positive or negative, with equal probability. For large clusters, the magnetisation obeys a large deviation principle. We show that the corresponding rate function is zero for values of the magnetisation that are intermediate between the two steady state values, which means that fluctuations with these values of the magnetisation are much less unlikely than previously thought. We show that their probabilities decay as power laws in the cluster size, instead of the exponential scaling that would be expected from the large deviation principle. We discuss how this observation is related to dynamical phase coexistence phenomena. We also comment on the typical size of magnetisation fluctuations. △ Less

Submitted 26 July, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 16 pages, authors accepted version

Journal ref: Phys. Rev. E 100, 012140 (2019)

Accelerated relaxation and suppressed dynamic heterogeneity in a kinetically constrained (East) model with swaps

Authors: Ricardo Gutierrez, Juan P. Garrahan, Robert L. Jack

Abstract: We introduce a kinetically constrained spin model with a local softness parameter, such that spin flips can violate the kinetic constraint with an (annealed) site-dependent rate. We show that adding MC swap moves to this model can dramatically accelerate structural relaxation. We discuss the connection of this observation with the fact that swap moves are also able to accelerate relaxation in stru… ▽ More We introduce a kinetically constrained spin model with a local softness parameter, such that spin flips can violate the kinetic constraint with an (annealed) site-dependent rate. We show that adding MC swap moves to this model can dramatically accelerate structural relaxation. We discuss the connection of this observation with the fact that swap moves are also able to accelerate relaxation in structural glasses. We analyse the rates of relaxation in the model. We also show that the extent of dynamical heterogeneity is strongly suppressed by the swap moves. △ Less

Submitted 27 September, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: 20 pages, accepted version

Journal ref: J. Stat. Mech. (2019) 094006

Microscopic analysis of thermo-orientation in systems of off-centre Lennard-Jones particles

Authors: Robert L. Jack, Peter Wirnsberger, Aleks Reinhardt

Abstract: When fluids of anisotropic molecules are placed in temperature gradients, the molecules may align themselves along the gradient: this is called thermo-orientation. We discuss the theory of this effect in a fluid of particles that interact by a spherically symmetric potential, where the particles' centres of mass do not coincide with their interaction centres. Starting from the equations of motion… ▽ More When fluids of anisotropic molecules are placed in temperature gradients, the molecules may align themselves along the gradient: this is called thermo-orientation. We discuss the theory of this effect in a fluid of particles that interact by a spherically symmetric potential, where the particles' centres of mass do not coincide with their interaction centres. Starting from the equations of motion of the molecules, we show how a simple assumption of local equipartition of energy can be used to predict the thermo-orientation effect, recovering the result of Wirnsberger et al. [Phys. Rev. Lett. 120, 226001 (2018)]. Within this approach, we show that for particles with a single interaction centre, the thermal centre of the molecule must coincide with the interaction centre. The theory also explains the coupling between orientation and kinetic energy that is associated with this non-Boltzmann distribution. We discuss deviations from this local equipartition assumption, showing that these can occur in linear response to a temperature gradient. We also present numerical simulations showing significant deviations from the local equipartition predictions, which increase as the centre of mass of the molecule is displaced further from its interaction centre. △ Less

Submitted 19 March, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 15 pages, authors' accepted version

Journal ref: J. Chem. Phys. 150, 134501 (2019)

Effects of random pinning on the potential energy landscape of a supercooled liquid

Authors: S. P. Niblett, V. K. de Souza, R. L. Jack, D. J Wales

Abstract: We use energy landscape methods to investigate the response of a supercooled liquid to random pinning. We classify the structural similarity of different energy minima using a measure of overlap. This analysis reveals a correspondence between distinct particle packings (which are characterised via the overlap) and funnels on the energy landscape (which are characterised via disconnectivity graphs)… ▽ More We use energy landscape methods to investigate the response of a supercooled liquid to random pinning. We classify the structural similarity of different energy minima using a measure of overlap. This analysis reveals a correspondence between distinct particle packings (which are characterised via the overlap) and funnels on the energy landscape (which are characterised via disconnectivity graphs). As the number of pinned particles is increased, we find a crossover from glassy behavior at low pinning to a structure-seeking landscape at high pinning, in which all thermally accessible minima are structurally similar. We discuss the consequences of these results for theories of randomly pinned liquids. We also investigate how the energy landscape depends on the fraction of pinned particles, including the degree of frustration and the evolution of distinct packings as the number of pinned particles is reduced. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Comments: Supporting research data may be found at https://doi.org/10.17863/CAM.30987

Journal ref: J. Chem. Phys. 149, 114503 (2018)

Unravelling the large deviation statistics of Markovian open quantum systems

Authors: Federico Carollo, Robert L. Jack, Juan P. Garrahan

Abstract: We analyse dynamical large deviations of quantum trajectories in Markovian open quantum systems in their full generality. We derive a {\em quantum level-2.5 large deviation principle} for these systems, which describes the joint fluctuations of time-averaged quantum jump rates and of the time-averaged quantum state for long times. Like its level-2.5 counterpart for classical continuous-time Markov… ▽ More We analyse dynamical large deviations of quantum trajectories in Markovian open quantum systems in their full generality. We derive a {\em quantum level-2.5 large deviation principle} for these systems, which describes the joint fluctuations of time-averaged quantum jump rates and of the time-averaged quantum state for long times. Like its level-2.5 counterpart for classical continuous-time Markov chains (which it contains as a special case) this description is both {\em explicit and complete}, as the statistics of arbitrary time-extensive dynamical observables can be obtained by contraction from the explicit level-2.5 rate functional we derive. Our approach uses an unravelled representation of the quantum dynamics which allows these statistics to be obtained by analysing a classical stochastic process in the space of pure states. For quantum reset processes we show that the unravelled dynamics is semi-Markov, and derive bounds on the asymptotic variance of the number of quantum jumps which generalise classical thermodynamic uncertainty relations. We finish by discussing how our level-2.5 approach can be used to study large deviations of non-linear functions of the state such as measures of entanglement. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 5+12 pages, 1 figure

Journal ref: Phys. Rev. Lett. 122, 130605 (2019)

Optimizing active work: Dynamical phase transitions, collective motion, and jamming

Authors: Takahiro Nemoto, Étienne Fodor, Michael E. Cates, Robert L. Jack, Julien Tailleur

Abstract: Active work measures how far the local self-forcing of active particles translates into real motion. Using Population Monte Carlo methods, we investigate large deviations in the active work for repulsive active Brownian disks. Minimizing the active work generically results in dynamical arrest; in contrast, despite the lack of aligning interactions, trajectories of high active work correspond to a… ▽ More Active work measures how far the local self-forcing of active particles translates into real motion. Using Population Monte Carlo methods, we investigate large deviations in the active work for repulsive active Brownian disks. Minimizing the active work generically results in dynamical arrest; in contrast, despite the lack of aligning interactions, trajectories of high active work correspond to a collectively moving, aligned state. We use heuristic and analytic arguments to explain the origin of dynamical phase transitions separating the arrested, typical, and aligned regimes. △ Less

Submitted 9 February, 2019; v1 submitted 8 May, 2018; originally announced May 2018.

Comments: 14 pages, 9 figures

Journal ref: Phys. Rev. E 99, 022605 (2019)

A Variational Structure for Interacting Particle Systems and their Hydrodynamic Scaling Limits

Authors: Marcus Kaiser, Robert L. Jack, Johannes Zimmer

Abstract: We consider hydrodynamic scaling limits for a class of reversible interacting particle systems, which includes the symmetric simple exclusion process and certain zero-range processes. We study a (non-quadratic) microscopic action functional for these systems. We analyse the behaviour of this functional in the hydrodynamic limit and we establish conditions under which it converges to the (quadratic… ▽ More We consider hydrodynamic scaling limits for a class of reversible interacting particle systems, which includes the symmetric simple exclusion process and certain zero-range processes. We study a (non-quadratic) microscopic action functional for these systems. We analyse the behaviour of this functional in the hydrodynamic limit and we establish conditions under which it converges to the (quadratic) action functional of Macroscopic Fluctuation Theory. We discuss the implications of these results for rigorous analysis of hydrodynamic limits. △ Less

Submitted 18 December, 2018; v1 submitted 3 May, 2018; originally announced May 2018.

Efficient characterisation of large deviations using population dynamics

Authors: Tobias Brewer, Stephen R. Clark, Russell Bradford, Robert L. Jack

Abstract: We consider population dynamics as implemented by the cloning algorithm for analysis of large deviations of time-averaged quantities. Using the simple symmetric exclusion process as a prototypical example, we investigate the convergence of the results with respect to the algorithmic parameters, focussing on the dynamical phase transition between homogeneous and inhomogeneous states, where converge… ▽ More We consider population dynamics as implemented by the cloning algorithm for analysis of large deviations of time-averaged quantities. Using the simple symmetric exclusion process as a prototypical example, we investigate the convergence of the results with respect to the algorithmic parameters, focussing on the dynamical phase transition between homogeneous and inhomogeneous states, where convergence is relatively difficult to achieve. We discuss how the performance of the algorithm can be optimised, and how it can be efficiently exploited on parallel computing platforms. △ Less

Submitted 9 May, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 23 pages, final version

Journal ref: J. Stat. Mech. (2018) 053204

Symmetries and Geometrical Properties of Dynamical Fluctuations in Molecular Dynamics

Authors: Robert L. Jack, Marcus Kaiser, Johannes Zimmer

Abstract: We describe some general results that constrain the dynamical fluctuations that can occur in non-equilibrium steady states, with a focus on molecular dynamics. That is, we consider Hamiltonian systems, coupled to external heat baths, and driven out of equilibrium by non-conservative forces. We focus on the probabilities of rare events (large deviations). First, we discuss a PT (parity-time) symmet… ▽ More We describe some general results that constrain the dynamical fluctuations that can occur in non-equilibrium steady states, with a focus on molecular dynamics. That is, we consider Hamiltonian systems, coupled to external heat baths, and driven out of equilibrium by non-conservative forces. We focus on the probabilities of rare events (large deviations). First, we discuss a PT (parity-time) symmetry that appears in ensembles of trajectories where a current is constrained to have a large (non-typical) value. We analyse the heat flow in such ensembles, and compare it with non-equilibrium steady states. Second, we consider pathwise large deviations that are defined by considering many copies of a system. We show how the probability currents in such systems can be decomposed into orthogonal contributions, that are related to convergence to equilibrium and to dissipation. We discuss the implications of these results for modelling non-equilibrium steady states. △ Less

Submitted 23 October, 2017; v1 submitted 14 September, 2017; originally announced September 2017.

Comments: final version, 21 pages

Canonical structure and orthogonality of forces and currents in irreversible Markov chains

Authors: Marcus Kaiser, Robert L. Jack, Johannes Zimmer

Abstract: We discuss a canonical structure that provides a unifying description of dynamical large deviations for irreversible finite state Markov chains (continuous time), Onsager theory, and Macroscopic Fluctuation Theory. For Markov chains, this theory involves a non-linear relation between probability currents and their conjugate forces. Within this framework, we show how the forces can be split into tw… ▽ More We discuss a canonical structure that provides a unifying description of dynamical large deviations for irreversible finite state Markov chains (continuous time), Onsager theory, and Macroscopic Fluctuation Theory. For Markov chains, this theory involves a non-linear relation between probability currents and their conjugate forces. Within this framework, we show how the forces can be split into two components, which are orthogonal to each other, in a generalised sense. This splitting allows a decomposition of the pathwise rate function into three terms, which have physical interpretations in terms of dissipation and convergence to equilibrium. Similar decompositions hold for rate functions at level 2 and level 2.5. These results clarify how bounds on entropy production and fluctuation theorems emerge from the underlying dynamical rules. We discuss how these results for Markov chains are related to similar structures within Macroscopic Fluctuation Theory, which describes hydrodynamic limits of such microscopic models.c Fluctuation Theory, which describes hydrodynamic limits of such microscopic models. △ Less

Submitted 16 February, 2018; v1 submitted 4 August, 2017; originally announced August 2017.

Comments: 36 page, 1 figure

MSC Class: 82C22; 82C35; 60J27; 60F10

Journal ref: Kaiser, M., Jack, R.L. & Zimmer, J. J Stat Phys (2018)