We present an analysis of the stationary distributions of run-and-tumble particles trapped in external potentials in terms of a thermophoretic potential that emerges when trapped active motion is mapped to trapped passive Brownian motion in a fictitious inhomogeneous thermal bath. We elaborate on the meaning of the non-Boltzmann-Gibbs stationary distributions that emerge as a consequence of the persistent motion of active particles. These stationary distributions are interpreted as a class of distributions in nonequilibrium statistical mechanics known as superstatistics. Our analysis provides an original insight on the link between the intrinsic nonequilibrium nature of active motion and the well-known concept of local equilibrium used in nonequilibrium statistical mechanics and contributes to the understanding of the validity of the concept of effective temperature. Particular cases of interest, regarding specific trapping potentials used in other theoretical or experimental studies, are discussed. We point out as an unprecedented effect, the emergence of new modes of the stationary distribution as a consequence of the coupling of persistent motion in a trapping potential that varies highly enough with position.

The diffusion of chiral active Brownian particles in three-dimensional space is studied analytically, by consideration of the corresponding Fokker-Planck equation for the probability density of finding a particle at position $\mathit{x}$ and moving along the direction $\widehat{\mathit{v}}$ at time $t$, and numerically, by the use of Langevin dynamics simulations. The analysis is focused on the marginal probability density of finding a particle at a given location and at a given time (independently of its direction of motion), which is found from an infinite hierarchy of differential-recurrence relations for the coefficients that appear in the multipole expansion of the probability distribution, which contains the whole kinematic information. This approach allows the explicit calculation of the time dependence of the mean-squared displacement and the time dependence of the kurtosis of the marginal probability distribution, quantities from which the effective diffusion coefficient and the “shape” of the positions distribution are examined. Oscillations between two characteristic values were found in the time evolution of the kurtosis, namely, between the value that corresponds to a Gaussian and the one that corresponds to a distribution of spherical shell shape. In the case of an ensemble of particles, each one rotating around a uniformly distributed random axis, evidence is found of the so-called effect “anomalous, yet Brownian, diffusion,” for which particles follow a non-Gaussian distribution for the positions yet the mean-squared displacement is a linear function of time.

We investigate a one-dimensional model of active motion, which takes into account the effects of persistent self-propulsion through a memory function in a dissipative-like term of the generalized Langevin equation for particle swimming velocity. The proposed model is a generalization of the active Ornstein-Uhlenbeck model introduced by G. Szamel [ Phys. Rev. E 90, 012111 (2014)]. We focus on two different kinds of memory which arise in many natural systems: an exponential decay and a power law, supplemented with additive colored noise. We provide analytical expressions for the velocity autocorrelation function and the mean-squared displacement, which are in excellent agreement with numerical simulations. For both models, damped oscillatory solutions emerge due to the competition between the memory of the system and the persistence of velocity fluctuations. In particular, for a power-law model with fractional Brownian noise, we show that long-time active subdiffusion occurs with increasing long-term memory.

We study the free diffusion in two dimensions of active Brownian swimmers subject to passive fluctuations on the translational motion and to active fluctuations on the rotational one. The Smoluchowski equation is derived from a Langevin-like model of active swimmers and analytically solved in the long-time regime for arbitrary values of the Péclet number; this allows us to analyze the out-of-equilibrium evolution of the positions distribution of active particles at all time regimes. Explicit expressions for the mean-square displacement and for the kurtosis of the probability distribution function are presented and the effects of persistence discussed. We show through Brownian dynamics simulations that our prescription for the mean-square displacement gives the exact time dependence at all times. The departure of the probability distribution from a Gaussian, measured by the kurtosis, is also analyzed both analytically and computationally. We find that for the inverse of Péclet numbers $\lesssim 0.1$, the distance from Gaussian increases as $\sim t^{-2}$ at short times, while it diminishes as $\sim t^{-1}$ in the asymptotic limit.

Starting from a Langevin description of active particles that move with constant speed in infinite two-dimensional space and its corresponding Fokker-Planck equation, we develop a systematic method that allows us to obtain the coarse-grained probability density of finding a particle at a given location and at a given time in arbitrary short-time regimes. By going beyond the diffusive limit, we derive a generalization of the telegrapher equation. Such generalization preserves the hyperbolic structure of the equation and incorporates memory effects in the diffusive term. While no difference is observed for the mean-square displacement computed from the two-dimensional telegrapher equation and from our generalization, the kurtosis results in a sensible parameter that discriminates between both approximations. We carry out a comparative analysis in Fourier space that sheds light on why the standard telegrapher equation is not an appropriate model to describe the propagation of particles with constant speed in dispersive media.

A theoretical analysis of active motion on curved surfaces is presented in terms of a generalization of the telegrapher equation. Such a generalized equation is explicitly derived as the polar approximation of the hierarchy of equations obtained from the corresponding Fokker-Planck equation of active particles diffusing on curved surfaces. The general solution to the generalized telegrapher equation is given for a pulse with vanishing current as initial data. Expressions for the probability density and the mean squared geodesic displacement are given in the limit of weak curvature. As an explicit example of the formulated theory, the case of active motion on the sphere is presented, where oscillations observed in the mean squared geodesic displacement are explained.