Statistical mechanics provides a framework to understand how microscopic dynamics of active particles give rise to emergent macroscopic behavior in nonequilibrium systems. Unlike passive systems at equilibrium, active matter continuously consumes energy, leading to rich phenomena such as non-Boltzmann steady states, enhanced transport, and anomalous fluctuations.
In this project, we develop theoretical and computational approaches to characterize the statistical properties of active systems, with a focus on first-passage processes, transport, and fluctuations. Of particular interest are quantities such as mean first passage times, splitting probabilities, and effective transport coefficients, which govern how active particles explore and escape complex environments. These problems are studied using tools from stochastic processes, asymptotic analysis, and kinetic theory, with applications to confined geometries, external fields, and heterogeneous media.
References
2026
-
Splitting probabilities of confined active particles
Sarafa A. Iyaniwura, and Zhiwei Peng
New Journal of Physics (in press), Sep 2026
Active particles exhibit self-propulsion, leading to transport behavior that differs fundamentally from passive Brownian motion. In confined or structured domains, activity strongly influence escape probabilities and first-passage behavior. Understanding these effects is essential for describing transport in biological microenvironments, microfluidic devices, and heterogeneous media. In this work, leveraging the backward Fokker–Planck equation, we investigate the splitting probability of active particles in confined domains, focusing on both a one-dimensional interval and a two-dimensional corrugated channel. Analytical solutions are derived for the one-dimensional case in various asymptotic regimes. In corrugated channels with small aspect ratios, we develop a Fick–Jacobs reduction that yields effective transport equations along the axial direction, whereas for finite aspect ratios, the splitting dynamics are characterized numerically. We demonstrate how channel geometry, particle activity, and chirality modulate the likelihood of escape through different boundaries. Our results provide quantitative predictions for the transport of active matter in complex environments and highlight the interplay between confinement and activity.
2025
-
Mean first passage time of active Brownian particles in two dimensions
Sarafa A. Iyaniwura, and Zhiwei Peng
New Journal of Physics, Oct 2025
The mean first passage time (MFPT) is a key metric for understanding transport, search, and escape processes in stochastic systems. While well characterized for passive Brownian particles, its behavior in active systems-such as active Brownian particles (ABPs)-remains less understood due to their self-propelled, nonequilibrium dynamics. In this paper, we formulate and analyze an elliptic partial differential equation (PDE) to characterize the MFPT of ABPs in two-dimensional domains, including circular, annular, and elliptical regions. For annular regions, we analyze the MFPT of ABPs under various boundary conditions. Our results reveal rich behaviors in the MFPT of ABPs that differ fundamentally from those of passive particles. Notably, the MFPT exhibits non-monotonic dependence on the initial position and orientation of the particle, with maxima often occurring away from the domain center. We also find that increasing swimming speed can either increase or decrease the MFPT depending on the geometry and initial orientation. Asymptotic analysis of the PDE in the weak-activity regime provides insight into how activity modifies escape statistics of the particles in different geometries. Finally, our numerical solutions of the PDE are validated against Monte Carlo simulations.
2024
-
Macrotransport of active particles in periodic channels and fields: rectification and dispersion
Zhiwei Peng
The Journal of Chemical Physics, Oct 2024
Transport and dispersion of active particles in structured environments such as corrugated channels and porous media are important for the understanding of both natural and engineered active systems. Owing to their continuous selfpropulsion, active particles exhibit rectified transport under spatially asymmetric confinement. While good progress has been made in experiments and particle-based simulations, a theoretical understanding of the effective long-time transport dynamics in spatially periodic geometries remains less developed. In this paper, we apply generalized Taylor dispersion theory (GTDT) to analyze the long-time effective transport dynamics of active Brownian particles (ABPs) in periodic channels and fields. We show that the long-time transport behavior is governed by an effective advection-diffusion equation. The derived macrotransport equations allow us to characterize the average drift and effective dispersion coefficient. For the case of ABPs in periodic channels without external fields, we show that regardless of activity, the average drift is given by the net diffusive flux along the channel. For ABPs, their activity is the driving mechanism that sustains a density gradient, which ultimately leads to rectified motion along the channel. Our continuum theory is validated against direct Brownian dynamics simulations of the Langevin equations governing the motion of each ABP.
-
Asymptotic analysis and simulation of mean first passage time for active Brownian particles in 1-D
Sarafa Adewale Iyaniwura, and Zhiwei Peng
SIAM Journal on Applied Mathematics, May 2024
Active Brownian particles (ABPs) are a model for nonequilibrium systems in which the constituent particles are self-propelled in addition to their Brownian motion. Compared to the well-studied mean first passage time (MFPT) of passive Brownian particles, the MFPT of ABPs is much less developed. In this paper, we study the MFPT for ABPs in a 1-D domain with absorbing boundary conditions at both ends of the domain. To reveal the effect of swimming on the MFPT, we consider an asymptotic analysis in the weak-swimming or small Péclet (Pe) number limit. In particular, analytical expressions for the survival probability and the MFPT are developed up to O(Pe2). We explore the effects of the starting positions and starting orientations on the MFPT. Our analysis shows that if the starting orientations are biased towards one side of the domain, the MFPT as a function of the starting position becomes asymmetric about the center of the domain. The analytical results were confirmed by the numerical solutions of the full PDE model.
2022
-
Trapped-particle microrheology of active suspensions
Zhiwei Peng, and John F. Brady
The Journal of Chemical Physics, Sep 2022
In microrheology, the local rheological properties, such as the viscoelasticity of a complex fluid, are inferred from the free or forced motion of embedded colloidal probe particles. Theoretical machinery developed for forced-probe microrheology of colloidal suspensions focused on either constant-force (CF) or constant-velocity (CV) probes, while in experiments, neither the force nor the kinematics of the probe is fixed. More importantly, the constraint of CF or CV introduces a difficulty in the meaningful quantification of the fluctuations of the probe due to a thermodynamic uncertainty relation. It is known that, for a Brownian particle trapped in a harmonic potential well, the product of the standard deviations of the trap force and the particle position is dkBT in d dimensions, with kBT being the thermal energy. As a result, if the force (position) is not allowed to fluctuate, the position (force) fluctuation becomes infinite. To allow the measurement of fluctuations in theoretical studies, in this work, we consider a microrheology model in which the embedded probe is dragged along by a moving harmonic potential so that both its position and the trap force are allowed to fluctuate. Starting from the full Smoluchowski equation governing the dynamics of N hard active Brownian particles, we derive a pair Smoluchowski equation describing the dynamics of the probe as it interacts with one bath particle by neglecting hydrodynamic interactions among particles in the dilute limit. From this, we determine the mean and the variance (i.e., fluctuation) of the probe position in terms of the pair probability distribution. We then characterize the behavior of the system in the limits of both weak and strong trap. By taking appropriate limits, we show that our generalized model can be reduced to the well-studied CF or CV microrheology models.