1
|
Chakraborty T, Pradhan P. Time-dependent properties of run-and-tumble particles. II. Current fluctuations. Phys Rev E 2024; 109:044135. [PMID: 38755901 DOI: 10.1103/physreve.109.044135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 03/19/2024] [Indexed: 05/18/2024]
Abstract
We investigate steady-state current fluctuations in two models of hardcore run-and-tumble particles (RTPs) on a periodic one-dimensional lattice of L sites, for arbitrary tumbling rate γ=τ_{p}^{-1} and density ρ; model I consists of standard hardcore RTPs, while model II is an analytically tractable variant of model I, called a long-ranged lattice gas (LLG). We show that, in the limit of L large, the fluctuation of cumulative current Q_{i}(T,L) across the ith bond in a time interval T≫1/D grows first subdiffusively and then diffusively (linearly) with T: 〈Q_{i}^{2}〉∼T^{α} with α=1/2 for 1/D≪T≪L^{2}/D and α=1 for T≫L^{2}/D, where D(ρ,γ) is the collective- or bulk-diffusion coefficient; at small times T≪1/D, exponent α depends on the details. Remarkably, regardless of the model details, the scaled bond-current fluctuations D〈Q_{i}^{2}(T,L)〉/2χL≡W(y) as a function of scaled variable y=DT/L^{2} collapse onto a universal scaling curve W(y), where χ(ρ,γ) is the collective particle mobility. In the limit of small density and tumbling rate, ρ,γ→0, with ψ=ρ/γ fixed, there exists a scaling law: The scaled mobility γ^{a}χ(ρ,γ)/χ^{(0)}≡H(ψ) as a function of ψ collapses onto a scaling curve H(ψ), where a=1 and 2 in models I and II, respectively, and χ^{(0)} is the mobility in the limiting case of a symmetric simple exclusion process; notably, the scaling function H(ψ) is model dependent. For model II (LLG), we calculate exactly, within a truncation scheme, both the scaling functions, W(y) and H(ψ). We also calculate spatial correlation functions for the current and compare our theory with simulation results of model I; for both models, the correlation functions decay exponentially, with correlation length ξ∼τ_{p}^{1/2} diverging with persistence time τ_{p}≫1. Overall, our theory is in excellent agreement with simulations and complements the prior findings [T. Chakraborty and P. Pradhan, Phys. Rev. E 109, 024124 (2024)1539-375510.1103/PhysRevE.109.024124].
Collapse
Affiliation(s)
- Tanmoy Chakraborty
- Department of Physics of Complex Systems, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700106, India
| | - Punyabrata Pradhan
- Department of Physics of Complex Systems, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700106, India
| |
Collapse
|
2
|
Khali SS, Peruani F, Chaudhuri D. When an active bath behaves as an equilibrium one. Phys Rev E 2024; 109:024120. [PMID: 38491633 DOI: 10.1103/physreve.109.024120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 01/22/2024] [Indexed: 03/18/2024]
Abstract
Active scalar baths consisting of active Brownian particles are characterized by a non-Gaussian velocity distribution, a kinetic temperature, and a diffusion coefficient that scale with the square of the active velocity v_{0}. While these results hold in overdamped active systems, inertial effects lead to normal velocity distributions, with kinetic temperature and diffusion coefficient increasing as ∼v_{0}^{α} with 1<α<2. Remarkably, the late-time diffusivity and mobility decrease with mass. Moreover, we show that the equilibrium Einstein relation is asymptotically recovered with inertia. In summary, the inertial mass restores an equilibriumlike behavior.
Collapse
Affiliation(s)
| | - Fernando Peruani
- LPTM, CY Cergy Paris Université, 2 Avenue A. Chauvin, 95302 Cergy-Pontoise Cedex, France
| | - Debasish Chaudhuri
- Institute of Physics, Sachivalaya Marg, Bhubaneswar 751005, India
- Max-Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, 01187 Dresden, Germany
| |
Collapse
|
3
|
Chakraborty T, Pradhan P. Time-dependent properties of run-and-tumble particles: Density relaxation. Phys Rev E 2024; 109:024124. [PMID: 38491605 DOI: 10.1103/physreve.109.024124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 01/03/2024] [Indexed: 03/18/2024]
Abstract
We characterize collective diffusion of hardcore run-and-tumble particles (RTPs) by explicitly calculating the bulk-diffusion coefficient D(ρ,γ) for arbitrary density ρ and tumbling rate γ, in systems on a d-dimensional periodic lattice. We study two minimal models of RTPs: Model I is the standard version of hardcore RTPs introduced in [Phys. Rev. E 89, 012706 (2014)10.1103/PhysRevE.89.012706], whereas model II is a long-ranged lattice gas (LLG) with hardcore exclusion, an analytically tractable variant of model I. We calculate the bulk-diffusion coefficient analytically for model II and numerically for model I through an efficient Monte Carlo algorithm; notably, both models have qualitatively similar features. In the strong-persistence limit γ→0 (i.e., dimensionless ratio r_{0}γ/v→0), with v and r_{0} being the self-propulsion speed and particle diameter, respectively, the fascinating interplay between persistence and interaction is quantified in terms of two length scales: (i) persistence length l_{p}=v/γ and (ii) a "mean free path," being a measure of the average empty stretch or gap size in the hopping direction. We find that the bulk-diffusion coefficient varies as a power law in a wide range of density: D∝ρ^{-α}, with exponent α gradually crossing over from α=2 at high densities to α=0 at low densities. As a result, the density relaxation is governed by a nonlinear diffusion equation with anomalous spatiotemporal scaling. In the thermodynamic limit, we show that the bulk-diffusion coefficient-for ρ,γ→0 with ρ/γ fixed-has a scaling form D(ρ,γ)=D^{(0)}F(ρav/γ), where a∼r_{0}^{d-1} is particle cross section and D^{(0)} is proportional to the diffusion coefficient of noninteracting particles; the scaling function F(ψ) is calculated analytically for model II (LLG) and numerically for model I. Our arguments are independent of dimensions and microscopic details.
Collapse
Affiliation(s)
- Tanmoy Chakraborty
- Department of Physics of Complex Systems, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700106, India
| | - Punyabrata Pradhan
- Department of Physics of Complex Systems, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700106, India
| |
Collapse
|
4
|
Zhang H, Zheng X. Invariable distribution of co-evolutionary complex adaptive systems with agent's behavior and local topological configuration. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2024; 21:3229-3261. [PMID: 38454726 DOI: 10.3934/mbe.2024143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
In this study, we developed a dynamical Multi-Local-Worlds (MLW) complex adaptive system with co-evolution of agent's behavior and local topological configuration to predict whether agents' behavior would converge to a certain invariable distribution and derive the conditions that should be satisfied by the invariable distribution of the optimal strategies in a dynamical system structure. To this end, a Markov process controlled by agent's behavior and local graphic topology configuration was constructed to describe the dynamic case's interaction property. After analysis, the invariable distribution of the system was obtained using the stochastic process method. Then, three kinds of agent's behavior (smart, normal, and irrational) coupled with corresponding behaviors, were introduced as an example to prove that their strategies converge to a certain invariable distribution. The results showed that an agent selected his/her behavior according to the evolution of random complex networks driven by preferential attachment and a volatility mechanism with its payment, which made the complex adaptive system evolve. We conclude that the corresponding invariable distribution was determined by agent's behavior, the system's topology configuration, the agent's behavior noise, and the system population. The invariable distribution with agent's behavior noise tending to zero differed from that with the population tending to infinity. The universal conclusion, corresponding to the properties of both dynamical MLW complex adaptive system and cooperative/non-cooperative game that are much closer to the common property of actual economic and management events that have not been analyzed before, is instrumental in substantiating managers' decision-making in the development of traffic systems, urban models, industrial clusters, technology innovation centers, and other applications.
Collapse
Affiliation(s)
- Hebing Zhang
- School of Intelligent Manufacture, Taizhou University, Jiaojiang 318000, Zhejiang, China
| | - Xiaojing Zheng
- School of Mathematical Sciences, Beihang University, Beijing 100191, China
| |
Collapse
|
5
|
Saintillan D. Dispersion of run-and-tumble microswimmers through disordered media. Phys Rev E 2023; 108:064608. [PMID: 38243487 DOI: 10.1103/physreve.108.064608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Accepted: 11/20/2023] [Indexed: 01/21/2024]
Abstract
Understanding the transport properties of microorganisms and self-propelled particles in porous media has important implications for human health as well as microbial ecology. In free space, most microswimmers perform diffusive random walks as a result of the interplay of self-propulsion and orientation decorrelation mechanisms such as run-and-tumble dynamics or rotational diffusion. In an unstructured porous medium, collisions with the microstructure result in a decrease in the effective spatial diffusivity of the particles from its free-space value. Here, we analyze this problem for a simple model system consisting of noninteracting point particles performing run-and-tumble dynamics through a two-dimensional disordered medium composed of a random distribution of circular obstacles, in the absence of Brownian diffusion or hydrodynamic interactions. The particles are assumed to collide with the obstacles as hard spheres and subsequently slide on the obstacle surface with no frictional resistance while maintaining their orientation, until they either escape or tumble. We show that the variations in the long-time diffusivity can be described by a universal dimensionless hindrance function f(ϕ,Pe) of the obstacle area fraction ϕ and Péclet number Pe, or ratio of the swimmer run length to the obstacle size. We analytically derive an asymptotic expression for the hindrance function valid for dilute media (Peϕ≪1), and its extension to denser media is obtained using stochastic simulations. As we explain, the model is also easily generalized to describe dispersion in three dimensions.
Collapse
Affiliation(s)
- David Saintillan
- Department of Mechanical and Aerospace Engineering, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093, USA
| |
Collapse
|
6
|
Smith NR. Nonequilibrium steady state of trapped active particles. Phys Rev E 2023; 108:L022602. [PMID: 37723780 DOI: 10.1103/physreve.108.l022602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Accepted: 07/26/2023] [Indexed: 09/20/2023]
Abstract
We consider an overdamped particle with a general physical mechanism that creates noisy active movement (e.g., a run-and-tumble particle or active Brownian particle, etc.), that is confined by an external potential. Focusing on the limit in which the correlation time τ of the active noise is small, we find the nonequilibrium steady-state distribution P_{st}(X) of the particle's position X. While typical fluctuations of X follow a Boltzmann distribution with an effective temperature that is not difficult to find, the tails of P_{st}(X) deviate from a Boltzmann behavior: In the limit τ→0, they scale as P_{st}(X)∼e^{-s(X)/τ}. We calculate the large-deviation function s(X) exactly for arbitrary trapping potential and active noise in dimension d=1, by relating it to the rate function that describes large deviations of the position of the same active particle in absence of an external potential at long times. We then extend our results to d>1 assuming rotational symmetry.
Collapse
Affiliation(s)
- Naftali R Smith
- Department of Environmental Physics, Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Sede Boqer Campus, 8499000, Israel
| |
Collapse
|
7
|
Rizkallah P, Sarracino A, Bénichou O, Illien P. Absolute Negative Mobility of an Active Tracer in a Crowded Environment. PHYSICAL REVIEW LETTERS 2023; 130:218201. [PMID: 37295085 DOI: 10.1103/physrevlett.130.218201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Revised: 02/17/2023] [Accepted: 04/11/2023] [Indexed: 06/12/2023]
Abstract
Absolute negative mobility (ANM) refers to the situation where the average velocity of a driven tracer is opposite to the direction of the driving force. This effect was evidenced in different models of nonequilibrium transport in complex environments, whose description remains effective. Here, we provide a microscopic theory for this phenomenon. We show that it emerges in the model of an active tracer particle submitted to an external force and which evolves on a discrete lattice populated with mobile passive crowders. Resorting to a decoupling approximation, we compute analytically the velocity of the tracer particle as a function of the different parameters of the system and confront our results to numerical simulations. We determine the range of parameters where ANM can be observed, characterize the response of the environment to the displacement of the tracer, and clarify the mechanism underlying ANM and its relationship with negative differential mobility (another hallmark of driven systems far from the linear response).
Collapse
Affiliation(s)
- Pierre Rizkallah
- Sorbonne Université, CNRS, Laboratoire de Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux (PHENIX), 4 Place Jussieu, 75005 Paris, France
| | - Alessandro Sarracino
- Dipartimento di Ingegneria, Università della Campania Luigi Vanvitelli, 81031 Aversa (CE), Italy
- Istituto dei Sistemi Complessi-CNR, P.le Aldo Moro 2, 00185, Rome, Italy
| | - Olivier Bénichou
- Sorbonne Université, CNRS, Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), 4 Place Jussieu, 75005 Paris, France
| | - Pierre Illien
- Sorbonne Université, CNRS, Laboratoire de Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux (PHENIX), 4 Place Jussieu, 75005 Paris, France
| |
Collapse
|
8
|
Plati A, Puglisi A, Sarracino A. Thermodynamic bounds for diffusion in nonequilibrium systems with multiple timescales. Phys Rev E 2023; 107:044132. [PMID: 37198828 DOI: 10.1103/physreve.107.044132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Accepted: 03/21/2023] [Indexed: 05/19/2023]
Abstract
We derive a thermodynamic uncertainty relation bounding the mean squared displacement of a Gaussian process with memory, driven out of equilibrium by unbalanced thermal baths and/or by external forces. Our bound is tighter with respect to previous results and also holds at finite time. We apply our findings to experimental and numerical data for a vibrofluidized granular medium, characterized by regimes of anomalous diffusion. In some cases our relation can distinguish between equilibrium and nonequilibrium behavior, a nontrivial inference task, particularly for Gaussian processes.
Collapse
Affiliation(s)
- A Plati
- Department of Physics, University of Rome Sapienza, Piazzale Aldo Moro 2, 00185, Rome, Italy
- Institute for Complex Systems-CNR, Piazzale Aldo Moro 2, 00185, Rome, Italy
- Université Paris-Saclay, CNRS, Laboratoire de Physique des Solides, 91405 Orsay, France
| | - A Puglisi
- Department of Physics, University of Rome Sapienza, Piazzale Aldo Moro 2, 00185, Rome, Italy
- Institute for Complex Systems-CNR, Piazzale Aldo Moro 2, 00185, Rome, Italy
- INFN, University of Rome Tor Vergata, Via della Ricerca Scientifica 1, 00133, Rome, Italy
| | - A Sarracino
- Institute for Complex Systems-CNR, Piazzale Aldo Moro 2, 00185, Rome, Italy
- Department of Engineering, University of Campania "Luigi Vanvitelli," 81031 Aversa (CE), Italy
| |
Collapse
|
9
|
Karan C, Chaudhuri D. Cooperation and competition in the collective drive by motor proteins: mean active force, fluctuations, and self-load. SOFT MATTER 2023; 19:1834-1843. [PMID: 36789956 DOI: 10.1039/d2sm01183b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
We consider the dynamics of a bio-filament under the collective drive of motor proteins. They are attached irreversibly to a substrate and undergo stochastic attachment-detachment with the filament to produce a directed force on it. We establish the dependence of the mean directed force and force correlations on the parameters describing the individual motor proteins using analytical theory and direct numerical simulations. The effective Langevin description for the filament motion gives mean-squared displacement, asymptotic diffusion constant, and mobility leading to an effective temperature. Finally, we show how competition between motor protein extensions generates a self-load, describable in terms of the effective temperature, affecting the filament motion.
Collapse
Affiliation(s)
- Chitrak Karan
- Institute of Physics, Sachivalaya Marg, Sainik School, Bhubaneswar, 751005, India.
- Homi Bhaba National Institute, Training School Complex, Anushakti Nagar, Mumbai, 400094, India.
| | - Debasish Chaudhuri
- Institute of Physics, Sachivalaya Marg, Sainik School, Bhubaneswar, 751005, India.
- Homi Bhaba National Institute, Training School Complex, Anushakti Nagar, Mumbai, 400094, India.
| |
Collapse
|