Inter-particle ratchet effect determines global current of heterogeneous particles diffusing in confinement

In a model of $N$ volume-excluding spheres in a $d$-dimensional tube, we consider how differences between particles in their drift velocities, diffusivities, and sizes influence the steady state distribution and axial particle current. We show that the model is exactly solvable when the geometrical constraints prevent any particle from overtaking every other -- a notion we term quasi-one-dimensionality. Then, due to a ratchet effect, the current is biased towards the velocities of the least diffusive particles. We consider special cases of this model in one dimension, and derive the exact joint gap distribution for driven tracers in a passive bath. We describe the relationship between phase space structure and irreversible drift that makes the quasi-one-dimensional supposition key to the model's solvability.Comment: 26 pages, 7 figure

A comparison of dynamical fluctuations of biased diffusion and run-and-tumble dynamics in one dimension

We compare the fluctuations in the velocity and in the fraction of time spent at a given position for minimal models of a passive and an active particle: an asymmetric random walker and a run-and-tumble particle in continuous time and on a 1D lattice. We compute rate functions and effective dynamics conditioned on large deviations for these observables. While generally different, for a unique and non-trivial choice of rates (up to a rescaling of time) the velocity rate functions for the two models become identical, whereas the effective processes generating the fluctuations remain distinct. This equivalence coincides with a remarkable parity of the spectra of the processes' generators. For the occupation-time problem, we show that both the passive and active particles undergo a prototypical dynamical phase transition when the average velocity is non-vanishing in the long-time limit.Comment: 27 pages, 10 figure

Exact spectral solution of two interacting run-and-tumble particles on a ring lattice

Exact solutions of interacting random walk models, such as 1D lattice gases, offer precise insight into the origin of nonequilibrium phenomena. Here, we study a model of run-and-tumble particles on a ring lattice interacting via hardcore exclusion. We present the exact solution for one and two particles using a generating function technique. For two particles, the eigenvectors and eigenvalues are explicitly expressed using two parameters reminiscent of Bethe roots, whose numerical values are determined by polynomial equations which we derive. The spectrum depends in a complicated way on the ratio of direction reversal rate to lattice jump rate, $\omega$. For both one and two particles, the spectrum consists of separate real bands for large $\omega$, which mix and become complex-valued for small $\omega$. At exceptional values of $\omega$, two or more eigenvalues coalesce such that the Markov matrix is non-diagonalizable. A consequence of this intricate parameter dependence is the appearance of dynamical transitions: non-analytic minima in the longest relaxation times as functions of $\omega$ (for a given lattice size). Exceptional points are theoretically and experimentally relevant in, e.g., open quantum systems and multichannel scattering. We propose that the phenomenon should be a ubiquitous feature of classical nonequilibrium models as well, and of relevance to physical observables in this context.Comment: 29 pages, 7 figures, revised submission to J. Stat. Mec

Parasites on parasites:Coupled fluctuations in stacked contact processes

We present a model for host-parasite dynamics which incorporates both vertical and horizontal transmission as well as spatial structure. Our model consists of stacked contact processes (CP), where the dynamics of the host is a simple CP on a lattice while the dynamics of the parasite is a secondary CP which sits on top of the host-occupied sites. In the simplest case, where infection does not incur any cost, we uncover a novel effect: a non-monotonic dependence of parasite prevalence on host turnover. Inspired by natural examples of hyperparasitism, we extend our model to multiple levels of parasites and identify a transition between the maintenance of a finite and infinite number of levels, which we conjecture is connected to a roughening transition in models of surface growth

From a microscopic solution to a continuum description of active particles with a recoil interaction in one dimension

We consider a model system of persistent random walkers that can jam, pass through each other or jump apart (recoil) on contact. In a continuum limit, where particle motion between stochastic changes in direction becomes deterministic, we find that the stationary inter-particle distribution functions are governed by an inhomogeneous fourth-order differential equation. Our main focus is on determining the boundary conditions that these distribution functions should satisfy. We find that these do not arise naturally from physical considerations, but need to be carefully matched to functional forms that arise from the analysis of an underlying discrete process. The inter-particle distribution functions, or their first derivatives, are generically found to be discontinuous at the boundaries.Comment: 16 pages; 5 figures; published in PR

Combinatorial mappings of exclusion processes

We review various combinatorial interpretations and mappings of stationary-state probabilities of the totally asymmetric, partially asymmetric and symmetric simple exclusion processes (TASEP, PASEP, SSEP respectively). In these steady states, the statistical weight of a configuration is determined from a matrix product, which can be written explicitly in terms of generalised ladder operators. This lends a natural association to the enumeration of random walks with certain properties. Specifically, there is a one-to-many mapping of steady-state configurations to a larger state space of discrete paths, which themselves map to an even larger state space of number permutations. It is often the case that the configuration weights in the extended space are of a relatively simple form (e.g., a Boltzmann-like distribution). Meanwhile, various physical properties of the nonequilibrium steady state - such as the entropy - can be interpreted in terms of how this larger state space has been partitioned. These mappings sometimes allow physical results to be derived very simply, and conversely the physical approach allows some new combinatorial problems to be solved. This work brings together results and observations scattered in the combinatorics and statistical physics literature, and also presents new results. The review is pitched at statistical physicists who, though not professional combinatorialists, are competent and enthusiastic amateurs.Comment: 56 pages, 21 figure