Showing 1–50 of 101 results for author: Voituriez, R

Integer defects, flow localization, and bistability on curved active surfaces

Authors: Rushikesh Shinde, Raphaël Voituriez, Andrew Callan-Jones

Abstract: Biological surfaces, such as developing epithelial tissues, exhibit in-plane polar or nematic order and can be strongly curved. Recently, integer (+1) topological defects have been identified as morphogenetic hotspots in living systems. Yet, while +1 defects in active matter on flat surfaces are well-understood, the general principles governing curved active defects remain unknown. Here, we study… ▽ More Biological surfaces, such as developing epithelial tissues, exhibit in-plane polar or nematic order and can be strongly curved. Recently, integer (+1) topological defects have been identified as morphogenetic hotspots in living systems. Yet, while +1 defects in active matter on flat surfaces are well-understood, the general principles governing curved active defects remain unknown. Here, we study the dynamics of integer defects in an extensile or contractile polar fluid on two types of morphogenetically-relevant substrates : (1) a cylinder terminated by a spherical cap, and (2) a bump on an otherwise flat surface. Because the Frank elastic energy on a curved surface generically induces a coupling to $\textit{deviatoric}$ curvature, $\mathcal{D}$ (difference between squared principal curvatures), a +1 defect is induced on both surface types. We find that $\mathcal{D}$ leads to surprising effects including localization of orientation gradients and active flows, and particularly for contractility, to hysteresis and bistability between quiescent and flowing defect states. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Long-term memory induced correction to Arrhenius law

Authors: A. Barbier-Chebbah, O. Bénichou, R. Voituriez, T. Guérin

Abstract: The Kramers escape problem is a paradigmatic model for the kinetics of rare events, which are usually characterized by Arrhenius law. So far, analytical approaches have failed to capture the kinetics of rare events in the important case of non-Markovian processes with long-term memory, as occurs in the context of reactions involving proteins, long polymers, or strongly viscoelastic fluids. Here, b… ▽ More The Kramers escape problem is a paradigmatic model for the kinetics of rare events, which are usually characterized by Arrhenius law. So far, analytical approaches have failed to capture the kinetics of rare events in the important case of non-Markovian processes with long-term memory, as occurs in the context of reactions involving proteins, long polymers, or strongly viscoelastic fluids. Here, based on a minimal model of non-Markovian Gaussian process with long-term memory, we determine quantitatively the mean FPT to a rare configuration and provide its asymptotics in the limit of a large energy barrier $E$. Our analysis unveils a correction to Arrhenius law, induced by long-term memory, which we determine analytically. This correction, which we show can be quantitatively significant, takes the form of a second effective energy barrier $E'<E$ and captures the dependence of rare event kinetics on initial conditions, which is a hallmark of long-term memory. Altogether, our results quantify the impact of long-term memory on rare event kinetics, beyond Arrhenius law. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 7 pages (main text) + 9 pages (SI)

Exact propagators of one-dimensional self-interacting random walks

Authors: Julien Brémont, Olivier Bénichou, Raphaël Voituriez

Abstract: Self-interacting random walks (SIRWs) show long-range memory effects that result from the interaction of the random walker at time $t$ with the territory already visited at earlier times $t'<t$. This class of non-Markovian random walks has applications in contexts as diverse as foraging theory, the behaviour of living cells, and even machine learning. Despite this importance and numerous theoretic… ▽ More Self-interacting random walks (SIRWs) show long-range memory effects that result from the interaction of the random walker at time $t$ with the territory already visited at earlier times $t'<t$. This class of non-Markovian random walks has applications in contexts as diverse as foraging theory, the behaviour of living cells, and even machine learning. Despite this importance and numerous theoretical efforts, the propagator, which is the distribution of the walker's position and arguably the most fundamental quantity to characterize the process, has so far remained out of reach for all but a single class of SIRW. Here we fill this gap and provide an exact and explicit expression for the propagator of two important classes of SIRWs, namely, the once-reinforced random walk and the polynomially self-repelling walk. These results give access to key observables, such as the diffusion coefficient, which so far had not been determined. We also uncover an inherently non-Markovian mechanism that tends to drive the walker away from its starting point. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Forming 1D Periodic J-aggregates by Mechanical Bending of BNNTs: Evidence of Activated Molecular Diffusion

Authors: J. -B. Marceau, D. -M Ta, A. Aguilar, A. Loiseau, R. Martel, P. Bon, R. Voituriez, G. Recher, E. Gaufrès

Abstract: Driving molecular assembly into micrometer-scale patterns is key for defining advanced materials of interest in various fields, including life sciences, photovoltaics, and quantum photonics. However, the driving process competes with other forces, such as Brownian motion, ripening phenomena, capillary forces, and non-specific adsorption. Here we report on a guided diffusion mechanism of luminescen… ▽ More Driving molecular assembly into micrometer-scale patterns is key for defining advanced materials of interest in various fields, including life sciences, photovoltaics, and quantum photonics. However, the driving process competes with other forces, such as Brownian motion, ripening phenomena, capillary forces, and non-specific adsorption. Here we report on a guided diffusion mechanism of luminescent dye molecules encapsulated inside boron nitride nanotubes (BNNTs). Correlative measurements between BNNT bending and molecular position along the BNNT axis reveal an efficient and long-range migration of dyes from curved to straight regions of the nanotube. This curvature activated diffusion forms clusters of bright J-aggregates in periodic patterns of well-defined spacing and length. A phenomenological model of guided molecular transport in bended BNNTs is used to describe this directed 1D diffusion inside BNNT. It is shown to accurately predict the position and morphologies of a J-aggregate as a function of nanotube length. Coupling topological stimuli to 1D molecular diffusion at the nanoscale is here presented as an interesting tool capable of reconfiguring various emissive patterns of functional molecules at the mesoscopic scale. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: Supplementary added page 21

Evidence and quantification of memory effects in competitive first passage events

Authors: M. Dolgushev, T. V. Mendes, B. Gorin, K. Xie, N. Levernier, O. Bénichou, H. Kellay, R. Voituriez, T. Guérin

Abstract: Splitting probabilities quantify the likelihood of a given outcome out of competitive events for general random processes. This key observable of random walk theory, historically introduced as the Gambler's ruin problem for a player in a casino, has a broad range of applications beyond mathematical finance in evolution genetics, physics and chemistry, such as allele fixation, polymer translocation… ▽ More Splitting probabilities quantify the likelihood of a given outcome out of competitive events for general random processes. This key observable of random walk theory, historically introduced as the Gambler's ruin problem for a player in a casino, has a broad range of applications beyond mathematical finance in evolution genetics, physics and chemistry, such as allele fixation, polymer translocation, protein folding and more generally competitive reactions. The statistics of competitive events is well understood for memoryless (Markovian) processes. However, in complex systems such as polymer fluids, the motion of a particle should typically be described as a process with memory. Appart from scaling theories and perturbative approaches in one-dimension, the outcome of competitive events is much less characterized analytically for processes with memory. Here, we introduce an analytical approach that provides the splitting probabilities for general $d$-dimensional non-Markovian Gaussian processes. This analysis shows that splitting probabilities are critically controlled by the out of equilibrium statistics of reactive trajectories, observed after the first passage. This hallmark of non-Markovian dynamics and its quantitative impact on splitting probabilities are directly evidenced in a prototypical experimental reaction scheme in viscoelastic fluids. Altogether, these results reveal both experimentally and theoretically the importance of memory effects on competitive reactions. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Target search kinetics for random walkers with memory

Authors: Olivier Bénichou, Thomas Guérin, Nicolas Levernier, Raphaël Voituriez

Abstract: In this chapter, we consider the problem of a non-Markovian random walker (displaying memory effects) searching for a target. We review an approach that links the first passage statistics to the properties of trajectories followed by the random walker in the future of the first passage time. This approach holds in one and higher spatial dimensions, when the dynamics in the vicinity of the target i… ▽ More In this chapter, we consider the problem of a non-Markovian random walker (displaying memory effects) searching for a target. We review an approach that links the first passage statistics to the properties of trajectories followed by the random walker in the future of the first passage time. This approach holds in one and higher spatial dimensions, when the dynamics in the vicinity of the target is Gaussian, and it is applied to three paradigmatic target search problems: the search for a target in confinement, the search for a rarely reached configuration (rare event kinetics), or the search for a target in infinite space, for processes featuring stationary increments or transient aging. The theory gives access to the mean first passage time (when it exists) or to the behavior of the survival probability at long times, and agrees with the available exact results obtained perturbatively for examples of weakly non-Markovian processes. This general approach reveals that the characterization of the non-equilibrium state of the system at the instant of first passage is key to derive first-passage kinetics, and provides a new methodology, via the analysis of trajectories after the first-passage, to make it quantitative. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 22 pages, invited chapter for the book "The Target Problem" (Eds. D. S. Grebenkov, R. Metzler, G. Oshanin)

Aging Dynamics of $d-$dimensional Locally Activated Random Walks

Authors: Julien Brémont, Theresa Jakuszeit, Olivier Bénichou, Raphael Voituriez

Abstract: Locally activated random walks are defined as random processes, whose dynamical parameters are modified upon visits to given activation sites. Such dynamics naturally emerge in living systems as varied as immune and cancer cells interacting with spatial heterogeneities in tissues, or at larger scales animals encountering local resources. At the theoretical level, these random walks provide an expl… ▽ More Locally activated random walks are defined as random processes, whose dynamical parameters are modified upon visits to given activation sites. Such dynamics naturally emerge in living systems as varied as immune and cancer cells interacting with spatial heterogeneities in tissues, or at larger scales animals encountering local resources. At the theoretical level, these random walks provide an explicit construction of strongly non Markovian, and aging dynamics. We propose a general analytical framework to determine various statistical properties characterizing the position and dynamical parameters of the random walker on $d$-dimensional lattices. Our analysis applies in particular to both passive (diffusive) and active (run and tumble) dynamics, and quantifies the aging dynamics and potential trapping of the random walker; it finally identifies clear signatures of activated dynamics for potential use in experimental data. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Extreme Value Statistics of Jump Processes

Authors: Jérémie Klinger, Raphaël Voituriez, Olivier Bénichou

Abstract: We investigate extreme value statistics (EVS) of general discrete time and continuous space symmetric jump processes. We first show that for unbounded jump processes, the semi-infinite propagator $G_0(x,n)$, defined as the probability for a particle issued from $0$ to be at position $x$ after $n$ steps whilst staying positive, is the key ingredient needed to derive a variety of joint distributions… ▽ More We investigate extreme value statistics (EVS) of general discrete time and continuous space symmetric jump processes. We first show that for unbounded jump processes, the semi-infinite propagator $G_0(x,n)$, defined as the probability for a particle issued from $0$ to be at position $x$ after $n$ steps whilst staying positive, is the key ingredient needed to derive a variety of joint distributions of extremes and times at which they are reached. Along with exact expressions, we extract novel universal asymptotic behaviors of such quantities. For bounded, semi-infinite jump processes killed upon first crossing of zero, we introduce the \textit{strip probability} $μ_{0,\underline{x}}(n)$, defined as the probability that a particle issued from 0 remains positive and reaches its maximum $x$ on its $n^{\rm th}$ step exactly. We show that $μ_{0,\underline{x}}(n)$ is the essential building block to address EVS of semi-infinite jump processes, and obtain exact expressions and universal asymptotic behaviors of various joint distributions. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 5 pages + 8 pages SM

Imperfect Narrow Escape problem

Authors: T. Guérin, M. Dolgushev, O. Bénichou, R. Voituriez

Abstract: We consider the kinetics of the imperfect narrow escape problem, i.e. the time it takes for a particle diffusing in a confined medium of generic shape to reach and to be adsorbed by a small, imperfectly reactive patch embedded in the boundary of the domain, in two or three dimensions. Imperfect reactivity is modeled by an intrinsic surface reactivity $κ$ of the patch, giving rise to Robin boundary… ▽ More We consider the kinetics of the imperfect narrow escape problem, i.e. the time it takes for a particle diffusing in a confined medium of generic shape to reach and to be adsorbed by a small, imperfectly reactive patch embedded in the boundary of the domain, in two or three dimensions. Imperfect reactivity is modeled by an intrinsic surface reactivity $κ$ of the patch, giving rise to Robin boundary conditions. We present a formalism to calculate the exact asymptotics of the mean reaction time in the limit of large volume of the confining domain. We obtain exact explicit results in the two limits of large and small reactivities of the reactive patch, and a semi-analytical expression in the general case. Our approach reveals an anomalous scaling of the mean reaction time as the inverse square root of the reactivity in the large reactivity limit, valid for an initial position near the extremity of the reactive patch. We compare our exact results with those obtained within the ``constant flux approximation''; we show that this approximation turns out to give exactly the next-to-leading order term of the small reactivity limit, and provides a good approximation of the reaction time far from the reactive patch for all reactivities, but not in the vicinity of the boundary of the reactive patch due to the above mentioned anomalous scaling. These results thus provide a general framework to quantify the mean reaction times for the imperfect narrow escape problem. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Leftward, Rightward and Complete Exit Time Distributions of Jump Processes

Authors: Jérémie Klinger, Raphaël Voituriez, Olivier Bénichou

Abstract: First-passage properties of continuous stochastic processes confined in a 1--dimensional interval are well described. However, for jump processes (discrete random walks), the characterization of the corresponding observables remains elusive, despite their relevance in various contexts. Here we derive exact asymptotic expressions for the leftward, rightward and complete exit time distributions from… ▽ More First-passage properties of continuous stochastic processes confined in a 1--dimensional interval are well described. However, for jump processes (discrete random walks), the characterization of the corresponding observables remains elusive, despite their relevance in various contexts. Here we derive exact asymptotic expressions for the leftward, rightward and complete exit time distributions from the interval $[0,x]$ for symmetric jump processes starting from $x_0=0$, in the large $x$ and large time limit. We show that both the leftward probability $F_{\underline{0},x}(n)$ to exit through $0$ at step $n$ and rightward probability $F_{0,\underline{x}}(n)$ to exit through $x$ at step $n$ exhibit a universal behavior dictated by the large distance decay of the jump distribution parameterized by the Levy exponent $μ$. In particular, we exhaustively describe the $n\ll x^μ$ and $n\gg x^μ$ limits and obtain explicit results in both regimes. Our results finally provide exact asymptotics for exit time distributions of jump processes in regimes where continuous limits do not apply. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 5 pages, 4 figures, Supplementary Material(5 pages)

Splitting Probabilities of Jump Processes

Authors: Jérémie Klinger, Raphaël Voituriez, Olivier Bénichou

Abstract: We derive a universal, exact asymptotic form of the splitting probability for symmetric continuous jump processes, which quantifies the probability $ π_{0,\underline{x}}(x_0)$ that the process crosses $x$ before 0 starting from a given position $x_0\in[0,x]$ in the regime $x_0\ll x$. This analysis provides in particular a fully explicit determination of the transmission probability ($x_0=0$), in s… ▽ More We derive a universal, exact asymptotic form of the splitting probability for symmetric continuous jump processes, which quantifies the probability $ π_{0,\underline{x}}(x_0)$ that the process crosses $x$ before 0 starting from a given position $x_0\in[0,x]$ in the regime $x_0\ll x$. This analysis provides in particular a fully explicit determination of the transmission probability ($x_0=0$), in striking contrast with the trivial prediction $ π_{0,\underline{x}}(0)=0$ obtained by taking the continuous limit of the process, which reveals the importance of the microscopic properties of the dynamics. These results are illustrated with paradigmatic models of jump processes with applications to light scattering in heterogeneous media in realistic 3$d$ slab geometries. In this context, our explicit predictions of the transmission probability, which can be directly measured experimentally, provide a quantitative characterization of the effective random process describing light scattering in the medium. △ Less

Submitted 31 January, 2022; originally announced January 2022.

Comments: 5 pages + 7 pages for supplementary information

Intermittent Attractive Interactions Lead to Microphase Separation in Non-motile Active Matter

Authors: Henry Alston, Andrew O. Parry, Raphaël Voituriez, Thibault Bertrand

Abstract: Non-motile active matter exhibits a wide range of non-equilibrium collective phenomena yet examples are crucially lacking in the literature. We present a microscopic model inspired by the bacteria {\it Neisseria Meningitidis} in which diffusive agents feel intermittent attractive forces. Through a formal coarse-graining procedure, we show that this truly scalar model of active matter exhibits the… ▽ More Non-motile active matter exhibits a wide range of non-equilibrium collective phenomena yet examples are crucially lacking in the literature. We present a microscopic model inspired by the bacteria {\it Neisseria Meningitidis} in which diffusive agents feel intermittent attractive forces. Through a formal coarse-graining procedure, we show that this truly scalar model of active matter exhibits the time-reversal-symmetry breaking terms defining the {\it Active Model B+} class. In particular, we confirm the presence of microphase separation by solving the kinetic equations numerically. We show that the switching rate controlling the interactions provides a regulation mechanism tuning the typical cluster size, e.g. in populations of bacteria interacting via type IV pili. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 7 pages (4 figures) of main text plus 12 pages (2 figures) of supplementary information

Everlasting impact of initial perturbations on first-passage times of non-Markovian random walks

Authors: N. Levernier, T. V. Mendes, O. Bénichou, R. Voituriez, T. Guérin

Abstract: Persistence, defined as the probability that a fluctuating signal has not reached a threshold up to a given observation time, plays a crucial role in the theory of random processes. It quantifies the kinetics of processes as varied as phase ordering, reaction diffusion or interface relaxation dynamics. The fact that persistence can decay algebraically with time with non trivial exponents has trigg… ▽ More Persistence, defined as the probability that a fluctuating signal has not reached a threshold up to a given observation time, plays a crucial role in the theory of random processes. It quantifies the kinetics of processes as varied as phase ordering, reaction diffusion or interface relaxation dynamics. The fact that persistence can decay algebraically with time with non trivial exponents has triggered a number of experimental and theoretical studies. However, general analytical methods to calculate persistence exponents cannot be applied to the ubiquitous case of non-Markovian systems relaxing transiently after an imposed initial perturbation. Here, we introduce a theoretical framework that enables the non perturbative determination of persistence exponents of $d$-dimensional Gaussian non-Markovian processes with general non stationary dynamics relaxing to a steady state after an initial perturbation. Two prototypical classes of situations are analyzed: either the system is subjected to a temperature quench at initial time, or its past trajectory is assumed to have been observed and thus known. Altogether, our results reveal and quantify, on the basis of Gaussian processes, the deep impact of initial perturbations on first-passage statistics of non-Markovian processes. Our theory covers the case of spatial dimension higher than one, opening the way to characterize non-trivial reaction kinetics for complex systems with non-equilibrium initial conditions. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Crosslinking and depletion determine spatial instabilities in cytoskeletal active matter

Authors: Guillaume Sarfati, Ananyo Maitra, Raphael Voituriez, Jean-Christophe Galas, André Estevez-Torres

Abstract: Active gels made of cytoskeletal proteins are valuable materials with attractive non-equilibrium properties such as spatial self-organization and self-propulsion. At least four typical routes to spatial patterning have been reported to date in different types of cytoskeletal active gels: bending and buckling instabilities in extensile systems, and global and local contraction instabilities in cont… ▽ More Active gels made of cytoskeletal proteins are valuable materials with attractive non-equilibrium properties such as spatial self-organization and self-propulsion. At least four typical routes to spatial patterning have been reported to date in different types of cytoskeletal active gels: bending and buckling instabilities in extensile systems, and global and local contraction instabilities in contractile gels. Here we report the observation of these four instabilities in a single type of active gel and we show that they are controlled by two parameters: the concentrations of ATP and depletion agent. We demonstrate that as the ATP concentration decreases, the concentration of passive motors increases until the gel undergoes a gelation transition. At this point, buckling is selected against bending, while global contraction is favored over local ones. Our observations are coherent with a hydrodynamic model of a viscoelastic active gel where the filaments are crosslinked with a characteristic time that diverges as the ATP concentration decreases. Our work thus provides a unified view of spatial instabilities in cytoskeletal active matter. △ Less

Submitted 11 January, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Universal kinetics of imperfect reactions in confinement

Authors: Thomas Guérin, Maxim Dolgushev, Olivier Bénichou, Raphaël Voituriez

Abstract: Chemical reactions generically require that particles come into contact. In practice, reaction is often imperfect and can necessitate multiple random encounters between reactants. In confined geometries, despite notable recent advances, there is to date no general analytical treatment of such imperfect transport-limited reaction kinetics. Here, we determine the kinetics of imperfect reactions in c… ▽ More Chemical reactions generically require that particles come into contact. In practice, reaction is often imperfect and can necessitate multiple random encounters between reactants. In confined geometries, despite notable recent advances, there is to date no general analytical treatment of such imperfect transport-limited reaction kinetics. Here, we determine the kinetics of imperfect reactions in confining domains for any diffusive or anomalously diffusive Markovian transport process, and for different models of imperfect reactivity. We show that the full distribution of reaction times is obtained in the large confining volume limit from the knowledge of the mean reaction time only, which we determine explicitly. This distribution for imperfect reactions is found to be identical to that of perfect reactions upon an appropriate rescaling of parameters, which highlights the robustness of our results. Strikingly, this holds true even in the regime of low reactivity where the mean reaction time is independent of the transport process, and can lead to large fluctuations of the reaction time even in simple reaction schemes. We illustrate our results for normal diffusion in domains of generic shape, and for anomalous diffusion in complex environments, where our predictions are confirmed by numerical simulations. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: to appear in Communications Chemistry

Joint statistics of space and time exploration of $1d$ random walks

Authors: J. Klinger, A. Barbier-Chebbah, R. Voituriez, O. Bénichou

Abstract: The statistics of first-passage times of random walks to target sites has proved to play a key role in determining the kinetics of space exploration in various contexts. In parallel, the number of distinct sites visited by a random walker and related observables have been introduced to characterize the geometry of space exploration. Here, we address the question of the joint distribution of the fi… ▽ More The statistics of first-passage times of random walks to target sites has proved to play a key role in determining the kinetics of space exploration in various contexts. In parallel, the number of distinct sites visited by a random walker and related observables have been introduced to characterize the geometry of space exploration. Here, we address the question of the joint distribution of the first-passage time to a target and the number of distinct sites visited when the target is reached, which fully quantifies the coupling between kinetics and geometry of search trajectories. Focusing on 1-dimensional systems, we present a general method and derive explicit expressions of this joint distribution for several representative examples of Markovian search processes. In addition, we obtain a general scaling form, which holds also for non Markovian processes and captures the general dependence of the joint distribution on its space and time variables. We argue that the joint distribution has important applications to various problems, such as a conditional form of the Rosenstock trapping model, and the persistence properties of self-interacting random walks. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Self-interacting random walks : aging, exploration and first-passage times

Authors: Alex Barbier--Chebbah, Olivier Benichou, Raphael Voituriez

Abstract: Self-interacting random walks are endowed with long range memory effects that emerge from the interaction of the random walker at time $t$ with the territory that it has visited at earlier times $t'<t$. This class of non Markovian random walks has applications in a broad range of examples, ranging from insects to living cells, where a random walker modifies locally its environment -- leaving behin… ▽ More Self-interacting random walks are endowed with long range memory effects that emerge from the interaction of the random walker at time $t$ with the territory that it has visited at earlier times $t'<t$. This class of non Markovian random walks has applications in a broad range of examples, ranging from insects to living cells, where a random walker modifies locally its environment -- leaving behind footprints along its path, and in turn responds to its own footprints. Because of their inherent non Markovian nature, the exploration properties of self-interacting random walks have remained elusive. Here we show that long range memory effects can have deep consequences on the dynamics of generic self-interacting random walks ; they can induce aging and non trivial persistence and transience exponents, which we determine quantitatively, in both infinite and confined geometries. Based on this analysis, we quantify the search kinetics of self-interacting random walkers and show that the distribution of the first-passage time (FPT) to a target site in a confined domain takes universal scaling forms in the large domain size limit, which we characterize quantitatively. We argue that memory abilities induced by attractive self-interactions provide a decisive advantage for local space exploration, while repulsive self-interactions can significantly accelerate the global exploration of large domains. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 14 pages, 7 figures

Reply to Comment on "Inverse Square Lévy Walks are not Optimal Search Strategies for d \geq 2 "

Authors: Nicolas Levernier, Johannes Textor, Olivier Bénichou, Raphaël Voituriez

Abstract: We refute here the concernes raised in the Comment of our letter. This reply states clearly the validity range of our results and shows that the optimality of inverse-square Levy walks at the basis of the Levy flight foraging hypothesis is generically unfounded. We also give the precise set of conditions for which inverse-levy square Levy walks turn to be optimal, conditions which are unlikely to… ▽ More We refute here the concernes raised in the Comment of our letter. This reply states clearly the validity range of our results and shows that the optimality of inverse-square Levy walks at the basis of the Levy flight foraging hypothesis is generically unfounded. We also give the precise set of conditions for which inverse-levy square Levy walks turn to be optimal, conditions which are unlikely to be verified biologically. △ Less

Submitted 23 March, 2021; originally announced March 2021.

Journal ref: Phys. Rev. Lett. 126, 048902, 2021

Clustering and ordering in cell assemblies with generic asymmetric aligning interactions

Authors: Thibault Bertrand, Joseph d'Alessandro, Ananyo Maitra, Shreyansh Jain, Barbara Mercier, René-Marc Mège, Benoit Ladoux, Raphaël Voituriez

Abstract: Collective cell migration plays an essential role in various biological processes, such as development or cancer proliferation. While cell-cell interactions are clearly key determinants of collective cell migration -- in addition to individual cells self-propulsion -- the physical mechanisms that control the emergence of cell clustering and collective cell migration are still poorly understood. In… ▽ More Collective cell migration plays an essential role in various biological processes, such as development or cancer proliferation. While cell-cell interactions are clearly key determinants of collective cell migration -- in addition to individual cells self-propulsion -- the physical mechanisms that control the emergence of cell clustering and collective cell migration are still poorly understood. In particular, observations have shown that binary cell-cell collisions generally lead to anti-alignement of cell polarities and separation of pairs -- a process called contact inhibition of locomotion (CIL), which is expected to disfavor the formation of large scale cell clusters with coherent motion. Here, we adopt a joint experimental and theoretical approach to determine the large scale dynamics of cell assemblies from elementary pairwise cell-cell interaction rules. We quantify experimentally binary cell-cell interactions and show that they can be captured by a minimal equilibrium-like pairwise asymmetric aligning interaction potential that reproduces the CIL phenomenology. We identify its symmetry class, build the corresponding active hydrodynamic theory and show on general grounds that such asymmetric aligning interaction destroys large scale clustering and ordering, leading instead to a liquid-like microphase of cell clusters of finite size and short lived polarity, or to a fully dispersed isotropic phase. Finally, this shows that CIL-like asymmetric interactions in cellular systems -- or general active systems -- control cluster sizes and polarity, and can prevent large scale coarsening and long range polarity, except in the singular regime of dense confluent systems. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 11 pages, 5 figures

Chiral active hexatics: Giant number fluctuations, waves and destruction of order

Authors: Ananyo Maitra, Martin Lenz, Raphael Voituriez

Abstract: Active materials, composed of internally driven particles, have properties that are qualitatively distinct from matter at thermal equilibrium. However, the most spectacular departures from equilibrium phase behaviour are thought to be confined to systems with polar or nematic asymmetry. In this paper, we show that such departures are also displayed in more symmetric phases such as hexatics if, in… ▽ More Active materials, composed of internally driven particles, have properties that are qualitatively distinct from matter at thermal equilibrium. However, the most spectacular departures from equilibrium phase behaviour are thought to be confined to systems with polar or nematic asymmetry. In this paper, we show that such departures are also displayed in more symmetric phases such as hexatics if, in addition, the constituent particles have chiral asymmetry. We show that chiral active hexatics whose rotation rate does not depend on density have giant number fluctuations. If the rotation rate depends on density, the giant number fluctuations are suppressed due to a novel orientation-density sound mode with a linear dispersion which propagates even in the overdamped limit. However, we demonstrate that beyond a finite but large lengthscale, a chirality and activity-induced relevant nonlinearity invalidates the predictions of the linear theory and destroys the hexatic order. In addition, we show that activity modifies the interactions between defects in the active chiral hexatic phase, making them non-mutual. Finally, to demonstrate the generality of a chiral active hexatic phase we show that it results from the melting of chiral active crystals in finite systems. △ Less

Submitted 28 December, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Kinetics of rare events for non-Markovian stationary processes and application to polymer dynamics

Authors: Nicolas Levernier, Olivier Bénichou, Raphaël Voituriez, Thomas Guérin

Abstract: How much time does it take for a fluctuating system, such as a polymer chain, to reach a target configuration that is rarely visited -- typically because of a high energy cost ? This question generally amounts to the determination of the first-passage time statistics to a target zone in phase space with lower occupation probability. Here, we present an analytical method to determine the mean first… ▽ More How much time does it take for a fluctuating system, such as a polymer chain, to reach a target configuration that is rarely visited -- typically because of a high energy cost ? This question generally amounts to the determination of the first-passage time statistics to a target zone in phase space with lower occupation probability. Here, we present an analytical method to determine the mean first-passage time of a generic non-Markovian random walker to a rarely visited threshold, which goes beyond existing weak-noise theories. We apply our method to polymer systems, to determine (i) the first time for a flexible polymer to reach a large extension, and (ii) the first closure time of a stiff inextensible wormlike chain. Our results are in excellent agreement with numerical simulations and provide explicit asymptotic laws for the mean first-passage times to rarely visited configurations. △ Less

Submitted 4 February, 2020; originally announced February 2020.

Comments: Accepted in Physical Review Research (Rapid Communication) Main text + SI

Journal ref: Phys. Rev. Research 2, 012057 (2020)

Inverse square Lévy walks are not optimal search strategies for $d\ge 2$

Authors: Nicolas Levernier, Olivier Benichou, Johannes Textor, Raphael Voituriez

Abstract: The Lévy hypothesis states that inverse square Lévy walks are optimal search strategies because they maximise the encounter rate with sparse, randomly distributed, replenishable targets. It has served as a theoretical basis to interpret a wealth of experimental data at various scales, from molecular motors to animals looking for resources, putting forward the conclusion that many living organisms… ▽ More The Lévy hypothesis states that inverse square Lévy walks are optimal search strategies because they maximise the encounter rate with sparse, randomly distributed, replenishable targets. It has served as a theoretical basis to interpret a wealth of experimental data at various scales, from molecular motors to animals looking for resources, putting forward the conclusion that many living organisms perform Lévy walks to explore space because of their optimal efficiency. Here we provide analytically the dependence on target density of the encounter rate of Lévy walks for any space dimension $d$ ; in particular, this scaling is shown to be {\it independent} of the Lévy exponent $α$ for the biologically relevant case $d\ge 2$, which proves that the founding result of the Lévy hypothesis is incorrect. As a consequence, we show that optimizing the encounter rate with respect to $α$ is {\it irrelevant} : it does not change the scaling with density and can lead virtually to {\it any} optimal value of $α$ depending on system dependent modeling choices. The conclusion that observed inverse square Lévy patterns are the result of a common selection process based purely on the kinetics of the search behaviour is therefore unfounded. △ Less

Submitted 5 February, 2020; v1 submitted 1 February, 2020; originally announced February 2020.

Comments: Accepted in Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 124, 080601, 2020

Enhanced orientational ordering induced by active yet isotropic bath

Authors: Ananyo Maitra, Raphael Voituriez

Abstract: Can a bath of isotropic but active particles promote ordering of anisotropic but passive particles? In this paper, we uncover a fluctuation-driven mechanism by which this is possible. Somewhat counter-intuitively, we show that the passive particles tend to be more ordered upon increasing the noise-strength of the active isotropic bath. We first demonstrate this in a general dynamical model for a n… ▽ More Can a bath of isotropic but active particles promote ordering of anisotropic but passive particles? In this paper, we uncover a fluctuation-driven mechanism by which this is possible. Somewhat counter-intuitively, we show that the passive particles tend to be more ordered upon increasing the noise-strength of the active isotropic bath. We first demonstrate this in a general dynamical model for a non-conserved order parameter (model A) coupled to an active isotropic field and then concentrate on two examples, i. a collection of polar rods on a substrate in an active isotropic bath and ii. a passive apolar suspension in a momentum conserved, actively forced but isotropic fluid which is relevant for current research in active systems. Our theory, which is relevant for understanding ordering transitions in out-of-equilibrium systems can be tested in experiments, for instance, by introducing a low concentration of passive rod-like objects in active isotropic fluids and, since it is applicable to any non-conserved dynamical field, may have applications far beyond active matter. △ Less

Submitted 27 September, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Journal ref: Phys. Rev. Lett. 124, 048003 (2020)

Survival probability of stochastic processes beyond persistence exponents

Authors: N. Levernier, M. Dolgushev, O. Bénichou, R. Voituriez, T. Guérin

Abstract: For many stochastic processes, the probability $S(t)$ of not-having reached a target in unbounded space up to time $t$ follows a slow algebraic decay at long times, $S(t)\sim S_0/t^θ$. This is typically the case of symmetric compact (i.e. recurrent) random walks. While the persistence exponent $θ$ has been studied at length, the prefactor $S_0$, which is quantitatively essential, remains poorly ch… ▽ More For many stochastic processes, the probability $S(t)$ of not-having reached a target in unbounded space up to time $t$ follows a slow algebraic decay at long times, $S(t)\sim S_0/t^θ$. This is typically the case of symmetric compact (i.e. recurrent) random walks. While the persistence exponent $θ$ has been studied at length, the prefactor $S_0$, which is quantitatively essential, remains poorly characterized, especially for non-Markovian processes. Here we derive explicit expressions for $S_0$ for a compact random walk in unbounded space by establishing an analytic relation with the mean first-passage time of the same random walk in a large confining volume. Our analytical results for $S_0$ are in good agreement with numerical simulations, even for strongly correlated processes such as Fractional Brownian Motion, and thus provide a refined understanding of the statistics of longest first-passage events in unbounded space. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Journal ref: Nature Communications 2019

Tunable corrugated patterns in an active gel sheet

Authors: Anis Senoussi, Shunnichi Kashida, Ananyo Maitra, Raphael Voituriez, Jean-Christophe Galas, André Estevez-Torres

Abstract: Active matter locally converts chemical energy into mechanical work and, for this reason, it provides new mechanisms of pattern formation. In particular, active gels made of protein motors and filaments are far-from-equilibrium systems that exhibit spontaneous flow,[Kruse2004, Voituriez2005] leading to active turbulence in two and three dimensions[Sanchez2012, Kumar2018] and coherent flow in three… ▽ More Active matter locally converts chemical energy into mechanical work and, for this reason, it provides new mechanisms of pattern formation. In particular, active gels made of protein motors and filaments are far-from-equilibrium systems that exhibit spontaneous flow,[Kruse2004, Voituriez2005] leading to active turbulence in two and three dimensions[Sanchez2012, Kumar2018] and coherent flow in three dimensions[Wu2017] (3D). Although these dynamic flows reveal a characteristic length scale resulting from the interplay between active forcing and passive restoring forces, the observation of static and long-range spatial patterns in active gels has remained elusive. In this work, we demonstrate that a 2D free-standing nematic active gel, formed spontaneously by depletion forces from a 3D solution of kinesin motors and microtubule filaments, actively buckles out-of-plane into a centimeter-sized periodic corrugated sheet that is stable for several days at low activity. Importantly, the corrugations are formed in the absence of flow and their wavelength and stability are controlled by the motor concentration, in agreement with a hydrodynamic theory. At higher activities these patterns are transient with the gel becoming turbulent at longer times. Our results underline the importance of both passive and active forces in shaping active gels and indicate that a static material can be sculpted through an active mechanism. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Diffusion through nanopores in connected lipid bilayer networks

Authors: Manon Valet, Léa-Laetitia Pontani, Raphaël Voituriez, Elie Wandersman, Alexis M. Prevost

Abstract: A biomimetic model of cell-cell communication was developed to probe the passive molecular transport across ion channels inserted in synthetic lipid bilayers formed between contacting droplets arranged in a linear array. Diffusion of a fluorescent probe across the array was measured for different pore concentrations. The diffusion characteristic time scale is found to vary non-linearly with the po… ▽ More A biomimetic model of cell-cell communication was developed to probe the passive molecular transport across ion channels inserted in synthetic lipid bilayers formed between contacting droplets arranged in a linear array. Diffusion of a fluorescent probe across the array was measured for different pore concentrations. The diffusion characteristic time scale is found to vary non-linearly with the pore concentration. Our measurements are successfully modeled by a continuous time random walk description, whose waiting time is the first exit time from a droplet through a cluster of pores. The size of the cluster of pores is found to increase with their concentration. Our results provide a direct link between the mesoscopic permeation properties and the microscopic characteristics of the pores such as their number, size and spatial arrangement. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 123, 088101 (2019)

Tracer diffusion in crowded narrow channels. Topical review

Authors: O. Benichou, P. Illien, G. Oshanin, A. Sarracino, R. Voituriez

Abstract: We summarise different results on the diffusion of a tracer particle in lattice gases of hard-core particles with stochastic dynamics, which are confined to narrow channels -- single-files, comb-like structures and quasi-one-dimensional channels with the width equal to several particle diameters. We show that in such geometries a surprisingly rich, sometimes even counter-intuitive, behaviour emerg… ▽ More We summarise different results on the diffusion of a tracer particle in lattice gases of hard-core particles with stochastic dynamics, which are confined to narrow channels -- single-files, comb-like structures and quasi-one-dimensional channels with the width equal to several particle diameters. We show that in such geometries a surprisingly rich, sometimes even counter-intuitive, behaviour emerges, which is absent in unbounded systems. We also present a survey of different results obtained for a tracer particle diffusion in unbounded systems, which will permit a reader to have an exhaustively broad picture of the tracer diffusion in crowded environments. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 32 pages, 4 figures, Journal of Physics: Condensed Matter, to appear

Dynamics of Run-and-Tumble Particles in Dense Single-File Systems

Authors: Thibault Bertrand, Pierre Illien, Olivier Bénichou, Raphaël Voituriez

Abstract: We study a minimal model of self-propelled particle in a crowded single-file environment. We extend classical models of exclusion processes (previously analyzed for diffusive and driven tracer particles) to the case where the tracer particle is a run-and-tumble particle (RTP), while all bath particles perform symmetric random walks. In the limit of high density of bath particles, we derive exact e… ▽ More We study a minimal model of self-propelled particle in a crowded single-file environment. We extend classical models of exclusion processes (previously analyzed for diffusive and driven tracer particles) to the case where the tracer particle is a run-and-tumble particle (RTP), while all bath particles perform symmetric random walks. In the limit of high density of bath particles, we derive exact expressions for the full distribution $\mathcal{P}_n(X)$ of the RTP position $X$ and all its cumulants, valid for arbitrary values of the tumbling probability $α$ and time $n$. Our results highlight striking effects of crowding on the dynamics: even cumulants of the RTP position are increasing functions of $α$ at intermediate timescales, and display a subdiffusive anomalous scaling $\propto \sqrt{n}$ independent of $α$ in the limit of long times $n\to \infty$. These analytical results set the ground for a quantitative analysis of experimental trajectories of real biological or artificial microswimmers in extreme confinement. △ Less

Submitted 11 July, 2018; originally announced July 2018.

Comments: 6 pages, 4 figures

Signatures of motor susceptibility in the dynamics of a tracer particle in an active gel

Authors: Nitzan Razin, Raphael Voituriez, Nir S. Gov

Abstract: We study a model for the motion of a tracer particle inside an active gel, exposing the properties of the van Hove distribution of the particle displacements. Active events of a typical force magnitude give rise to non-Gaussian distributions, having exponential tails or side-peaks. The side-peaks appear when the local bulk elasticity of the gel is large enough and few active sources are dominant.… ▽ More We study a model for the motion of a tracer particle inside an active gel, exposing the properties of the van Hove distribution of the particle displacements. Active events of a typical force magnitude give rise to non-Gaussian distributions, having exponential tails or side-peaks. The side-peaks appear when the local bulk elasticity of the gel is large enough and few active sources are dominant. We explain the regimes of the different distributions, and study the structure of the peaks for active sources that are susceptible to the elastic stress that they cause inside the gel. We show how the van Hove distribution is altered by both the duty cycle of the active sources and their susceptibility, and suggest it as a sensitive probe to analyze microrheology data in active systems with restoring elastic forces. △ Less

Submitted 23 June, 2018; originally announced June 2018.

Comments: 4 pages, 4 figures and supplemental information (5 pages, 4 figures)

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

Optimized Diffusion of Run-and-Tumble Particles in Crowded Environments

Authors: Thibault Bertrand, Yongfeng Zhao, Olivier Bénichou, Julien Tailleur, Raphaël Voituriez

Abstract: We study the transport of self-propelled particles in dynamic complex environments. To obtain exact results, we introduce a model of run-and-tumble particles (RTPs) moving in discrete time on a $d$-dimensional cubic lattice in the presence of diffusing hard core obstacles. We derive an explicit expression for the diffusivity of the RTP, which is exact in the limit of low density of fixed obstacles… ▽ More We study the transport of self-propelled particles in dynamic complex environments. To obtain exact results, we introduce a model of run-and-tumble particles (RTPs) moving in discrete time on a $d$-dimensional cubic lattice in the presence of diffusing hard core obstacles. We derive an explicit expression for the diffusivity of the RTP, which is exact in the limit of low density of fixed obstacles. To do so, we introduce a generalization of Kac's theorem on the mean return times of Markov processes, which we expect to be relevant for a large class of lattice gas problems. Our results show the diffusivity of RTPs to be nonmonotonic in the tumbling probability for low enough obstacle mobility. These results prove the potential for optimization of the transport of RTPs in crowded and disordered environments with applications to motile artificial and biological systems. △ Less

Submitted 11 July, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Journal ref: Phys. Rev. Lett. 120, 198103 (2018)

Nonequilibrium fluctuations and enhanced diffusion of a driven particle in a dense environment

Authors: Pierre Illien, Olivier Bénichou, Gleb Oshanin, Alessandro Sarracino, Raphaël Voituriez

Abstract: We study the diffusion of a tracer particle driven out-of-equilibrium by an external force and traveling in a dense environment of arbitrary density. The system evolves on a discrete lattice and its stochastic dynamics is described by a master equation. Relying on a decoupling approximation that goes beyond the naive mean-field treatment of the problem, we calculate the fluctuations of the positio… ▽ More We study the diffusion of a tracer particle driven out-of-equilibrium by an external force and traveling in a dense environment of arbitrary density. The system evolves on a discrete lattice and its stochastic dynamics is described by a master equation. Relying on a decoupling approximation that goes beyond the naive mean-field treatment of the problem, we calculate the fluctuations of the position of the tracer around its mean value on a lattice of arbitrary dimension, and with different boundary conditions. We reveal intrinsically nonequilibrium effects, such as enhanced diffusivity of the tracer induced both by the crowding interactions and the external driving. We finally consider the high-density and low-density limits of the model and show that our approximation scheme becomes exact in these limits. △ Less

Submitted 8 April, 2018; v1 submitted 6 September, 2017; originally announced September 2017.

Journal ref: Phys. Rev. Lett. 120, 200606 (2018)

Forces in inhomogeneous open active-particle systems

Authors: Nitzan Razin, Raphael Voituriez, Jens Elgeti, Nir S. Gov

Abstract: We study the force that non-interacting point-like active particles apply to a symmetric inert object in the presence of a gradient of activity and particle sources and sinks. We consider two simple patterns of sources and sinks that are common in biological systems. We analytically solve a one dimensional model designed to emulate higher dimensional systems, and study a two dimensional model by n… ▽ More We study the force that non-interacting point-like active particles apply to a symmetric inert object in the presence of a gradient of activity and particle sources and sinks. We consider two simple patterns of sources and sinks that are common in biological systems. We analytically solve a one dimensional model designed to emulate higher dimensional systems, and study a two dimensional model by numerical simulations. We specify when the particle flux due to the creation and annihilation of particles can act to smooth the density profile that is induced by a gradient in the velocity of the active particles, and find the net resultant force due to both the gradient in activity and the particle flux. These results are compared qualitatively to observations of nuclear motion inside the oocyte, that is driven by a gradient in activity of actin-coated vesicles. △ Less

Submitted 4 March, 2018; v1 submitted 17 August, 2017; originally announced August 2017.

Comments: 14 pages, 8 figures

Journal ref: Phys. Rev. E 96, 052409 (2017)

Universal first-passage statistics of aging processes

Authors: N. Levernier, O. Bénichou, T. Guérin, R. Voituriez

Abstract: Many out of equilibrium phenomena, such as diffusion-limited reactions or target search processes, are controlled by first-passage events. So far the general determination of the mean first-passage time (FPT) to a target in confinement has left aside aging processes, involved in contexts as varied as glassy dynamics, tracer diffusion in biological membranes or transport of cold atoms in optical la… ▽ More Many out of equilibrium phenomena, such as diffusion-limited reactions or target search processes, are controlled by first-passage events. So far the general determination of the mean first-passage time (FPT) to a target in confinement has left aside aging processes, involved in contexts as varied as glassy dynamics, tracer diffusion in biological membranes or transport of cold atoms in optical lattices. Here we consider general non-Markovian scale-invariant processes in arbitrary dimension, displaying aging, and demonstrate that all the moments of the FPT obey universal scalings with the confining volume with non trivial exponents. Our analysis shows that a nonlinear scaling of the mean FPT with the volume is the hallmark of aging and provides a general tool to quantify its impact on first-passage kinetics in confinement. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Journal ref: Phys. Rev. E 98, 022125 (2018)

Generalized Archimedes' principle in active fluids

Authors: Nitzan Razin, Raphael Voituriez, Jens Elgeti, Nir S. Gov

Abstract: We show how a gradient in the motility properties of non-interacting point-like active particles can cause a pressure gradient that pushes a large inert object. We calculate the force on an object inside a system of active particles with position dependent motion parameters, in one and two dimensions, and show that a modified Archimedes' principle is satisfied. We characterize the system, both in… ▽ More We show how a gradient in the motility properties of non-interacting point-like active particles can cause a pressure gradient that pushes a large inert object. We calculate the force on an object inside a system of active particles with position dependent motion parameters, in one and two dimensions, and show that a modified Archimedes' principle is satisfied. We characterize the system, both in terms of the model parameters and in terms of experimentally measurable quantities: the spatial profiles of the density, velocity and pressure. This theoretical analysis is motivated by recent experiments, which showed that the nucleus of a mouse oocyte (immature egg cell) moves from the cortex to the center due to a gradient of activity of vesicles propelled by molecular motors; it more generally applies to artificial systems of controlled localized activity. △ Less

Submitted 29 September, 2017; v1 submitted 21 March, 2017; originally announced March 2017.

Comments: 16 pages, 9 figures

Journal ref: Phys. Rev. E 96, 032606 (2017)

Mean first-passage time of an anisotropic diffusive searcher

Authors: Nicolas Levernier, Olivier Bénichou, Raphaël Voituriez

Abstract: We consider an anisotropic needle-like Brownian particle with nematic symmetry confined in a $2D$ domain. For this system, the coupling of translational and rotational diffusion makes the process ${\bf x} (t)$ of the positions of the particle non Markovian. Using scaling arguments, a Gaussian approximation and numerical methods, we determine the mean first passage time $\langle\mathbf{T}\rangle$ o… ▽ More We consider an anisotropic needle-like Brownian particle with nematic symmetry confined in a $2D$ domain. For this system, the coupling of translational and rotational diffusion makes the process ${\bf x} (t)$ of the positions of the particle non Markovian. Using scaling arguments, a Gaussian approximation and numerical methods, we determine the mean first passage time $\langle\mathbf{T}\rangle$ of the particle to a target of radius $a$ and show in particular that $\langle\mathbf{T}\rangle \sim a^{-1/2}$ for $a\to 0$, in contrast with the classical logarithmic divergence obtained in the case of an isotropic $2D$ Brownian particle. △ Less

Submitted 1 February, 2017; v1 submitted 31 January, 2017; originally announced January 2017.

Journal ref: Journal of Physics A: Mathematical and Theoretical, Volume 50, Number 2 , 2016

Non-Markovian closure kinetics of flexible polymers with hydrodynamic interactions

Authors: Nicolas Levernier, Maxim Dolgushev, Olivier Bénichou, Alexander Blumen, Thomas Guérin, Raphaël Voituriez

Abstract: This paper presents a theoretical analysis of the closure kinetics of a polymer with hydrodynamic interactions. This analysis, which takes into account the non-Markovian dynamics of the end-to-end vector and relies on the preaveraging of the mobility tensor (Zimm dynamics), is shown to reproduce very accurately the results of numerical simulations of the complete non linear dynamics. It is found t… ▽ More This paper presents a theoretical analysis of the closure kinetics of a polymer with hydrodynamic interactions. This analysis, which takes into account the non-Markovian dynamics of the end-to-end vector and relies on the preaveraging of the mobility tensor (Zimm dynamics), is shown to reproduce very accurately the results of numerical simulations of the complete non linear dynamics. It is found that Markovian treatments, based on a Wilemski-Fixman approximation, significantly overestimate cyclization times (up to a factor 2), showing the importance of memory effects in the dynamics. In addition, this analysis provides scaling laws of the mean first cyclization time (MFCT) with the polymer size $N$ and capture radius $b$, which are identical in both Markovian and non-Markovian approaches. In particular, it is found that the scaling of the MFCT for large $N$ is given by $T\sim N^{3/2}\ln (N/b^2)$, which differs from the case of the Rouse dynamics where $T\sim N^{2}$. The extension to the case of the reaction kinetics of a monomer of a Zimm polymer with an external target in a confined volume is also presented. △ Less

Submitted 1 February, 2017; v1 submitted 31 January, 2017; originally announced January 2017.

Journal ref: The Journal of Chemical Physics > Volume 143, Issue 20, November 2015

Mean first-passage times of non-Markovian random walkers in confinement

Authors: T. Guérin, N. Levernier, O. Bénichou, R. Voituriez

Abstract: The first-passage time (FPT), defined as the time a random walker takes to reach a target point in a confining domain, is a key quantity in the theory of stochastic processes. Its importance comes from its crucial role to quantify the efficiency of processes as varied as diffusion-limited reactions, target search processes or spreading of diseases. Most methods to determine the FPT properties in c… ▽ More The first-passage time (FPT), defined as the time a random walker takes to reach a target point in a confining domain, is a key quantity in the theory of stochastic processes. Its importance comes from its crucial role to quantify the efficiency of processes as varied as diffusion-limited reactions, target search processes or spreading of diseases. Most methods to determine the FPT properties in confined domains have been limited to Markovian (memoryless) processes. However, as soon as the random walker interacts with its environment, memory effects can not be neglected. Examples of non Markovian dynamics include single-file diffusion in narrow channels or the motion of a tracer particle either attached to a polymeric chain or diffusing in simple or complex fluids such as nematics \cite{turiv2013effect}, dense soft colloids or viscoelastic solution. Here, we introduce an analytical approach to calculate, in the limit of a large confining volume, the mean FPT of a Gaussian non-Markovian random walker to a target point. The non-Markovian features of the dynamics are encompassed by determining the statistical properties of the trajectory of the random walker in the future of the first-passage event, which are shown to govern the FPT kinetics.This analysis is applicable to a broad range of stochastic processes, possibly correlated at long-times. Our theoretical predictions are confirmed by numerical simulations for several examples of non-Markovian processes including the emblematic case of the Fractional Brownian Motion in one or higher dimensions. These results show, on the basis of Gaussian processes, the importance of memory effects in first-passage statistics of non-Markovian random walkers in confinement. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: Submitted version. Supplementary Information can be found on the Nature website : http://www.nature.com/nature/journal/v534/n7607/full/nature18272.html

Journal ref: Nature {\bf 534}, 356 (2016)

Optimal search strategies of run-and-tumble walks

Authors: Jean-Francois Rupprecht, Olivier Bénichou, Raphael Voituriez

Abstract: The run-and-tumble walk, consisting in randomly reoriented ballistic excursions, models phenomena ranging from gas kinetics to bacteria motility. We evaluate the mean time required for this walk to find a fixed target within a 2D or 3D spherical confinement. We find that the mean search time admits a minimum as a function of the mean run duration for various types of boundary conditions and run du… ▽ More The run-and-tumble walk, consisting in randomly reoriented ballistic excursions, models phenomena ranging from gas kinetics to bacteria motility. We evaluate the mean time required for this walk to find a fixed target within a 2D or 3D spherical confinement. We find that the mean search time admits a minimum as a function of the mean run duration for various types of boundary conditions and run duration distributions (exponential, power-law, deterministic). Our result stands in sharp contrast to the pure ballistic motion, which is predicted to be the optimal search strategy in the case of Poisson distributed targets. △ Less

Submitted 15 July, 2016; v1 submitted 11 March, 2016; originally announced March 2016.

Comments: 10 pages, 7 figures

Journal ref: Physical Review E 94, 012117 (2016)

Active particles with soft and curved walls: Equation of state, ratchets, and instabilities

Authors: Nikolai Nikola, Alexandre P. Solon, Yariv Kafri, Mehran Kardar, Julien Tailleur, Raphaël Voituriez

Abstract: We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the non-uniform pressure along curved walls, an equation of state is recovered upon a proper spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, whic… ▽ More We study, from first principles, the pressure exerted by an active fluid of spherical particles on general boundaries in two dimensions. We show that, despite the non-uniform pressure along curved walls, an equation of state is recovered upon a proper spatial averaging. This holds even in the presence of pairwise interactions between particles or when asymmetric walls induce ratchet currents, which are accompanied by spontaneous shear stresses on the walls. For flexible obstacles, the pressure inhomogeneities lead to a modulational instability as well as to the spontaneous motion of short semi-flexible filaments. Finally, we relate the force exerted on objects immersed in active baths to the particle flux they generate around them. △ Less

Submitted 25 August, 2016; v1 submitted 17 December, 2015; originally announced December 2015.

Comments: 6 pages, 5 figures and 8 pages of supplementary information

Journal ref: Phys. Rev. Lett. 117, 098001 (2016)

Cover times of random searches

Authors: M. Chupeau, O. Bénichou, R. Voituriez

Abstract: How long does it take a random searcher to visit all sites of a given domain? This time, known as the cover time, is a key observable to quantify the efficiency of exhaustive searches, which require a complete exploration of an area and not only the discovery of a single target; examples range from immune system cells chasing pathogens to animals harvesting resources, robotized exploration by e.g.… ▽ More How long does it take a random searcher to visit all sites of a given domain? This time, known as the cover time, is a key observable to quantify the efficiency of exhaustive searches, which require a complete exploration of an area and not only the discovery of a single target; examples range from immune system cells chasing pathogens to animals harvesting resources, robotized exploration by e.g. automated cleaners or deminers, or algorithmics. Despite its broad relevance, the cover time has remained elusive and so far explicit results have been scarce and mostly limited to regular random walks. Here we determine the full distribution of the cover time for a broad range of random search processes, which includes the prominent examples of Lévy strategies, intermittent strategies, persistent random walks and random walks on complex networks, and reveal its universal features. We show that for all these examples the mean cover time can be minimized, and that the corresponding optimal strategies also minimize the mean search time for a single target, unambiguously pointing towards their robustness. △ Less

Submitted 19 November, 2015; originally announced November 2015.

Comments: submitted version

Journal ref: Nature Physics 11, 844-847 (2015)

Nonlinear response and emerging nonequilibrium micro-structures for biased diffusion in confined crowding environments

Authors: O. Bénichou, P. Illien, G. Oshanin, A. Sarracino, R. Voituriez

Abstract: We study analytically the dynamics and the micro-structural changes of a host medium caused by a driven tracer particle moving in a confined, quiescent molecular crowding environment. Imitating typical settings of active micro-rheology experiments, we consider here a minimal model comprising a geometrically confined lattice system -- a two-dimensional strip-like or a three-dimensional capillary-li… ▽ More We study analytically the dynamics and the micro-structural changes of a host medium caused by a driven tracer particle moving in a confined, quiescent molecular crowding environment. Imitating typical settings of active micro-rheology experiments, we consider here a minimal model comprising a geometrically confined lattice system -- a two-dimensional strip-like or a three-dimensional capillary-like -- populated by two types of hard-core particles with stochastic dynamics -- a tracer particle driven by a constant external force and bath particles moving completely at random. Resorting to a decoupling scheme, which permits us to go beyond the linear-response approximation (Stokes regime) for arbitrary densities of the lattice gas particles, we determine the force-velocity relation for the tracer particle and the stationary density profiles of the host medium particles around it. These results are validated a posteriori by extensive numerical simulations for a wide range of parameters. Our theoretical analysis reveals two striking features: a) We show that, under certain conditions, the terminal velocity of the driven tracer particle is a nonmonotonic function of the force, so that in some parameter range the differential mobility becomes negative, and b) the biased particle drives the whole system into a nonequilibrium steady-state with a stationary particle density profile past the tracer, which decays exponentially, in sharp contrast with the behavior observed for unbounded lattices, where an algebraic decay is known to take place. △ Less

Submitted 30 October, 2015; originally announced October 2015.

Comments: 14 pages, 8 figures

Journal ref: Phys. Rev. E 93, 032128 (2016)

Contact Kinetics in Fractal Macromolecules

Authors: Maxim Dolgushev, Thomas Guérin, Alexander Blumen, Olivier Bénichou, Raphaël Voituriez

Abstract: We consider the kinetics of first contact between two monomers of the same macromolecule. Relying on a fractal description of the macromolecule, we develop an analytical method to compute the Mean First Contact Time (MFCT) for various molecular sizes. In our theoretical description, the non-Markovian feature of monomer motion, arising from the interactions with the other monomers, is captured by a… ▽ More We consider the kinetics of first contact between two monomers of the same macromolecule. Relying on a fractal description of the macromolecule, we develop an analytical method to compute the Mean First Contact Time (MFCT) for various molecular sizes. In our theoretical description, the non-Markovian feature of monomer motion, arising from the interactions with the other monomers, is captured by accounting for the non-equilibrium conformations of the macromolecule at the very instant of first contact. This analysis reveals a simple scaling relation for the MFCT between two monomers, which involves only their equilibrium distance and the spectral dimension of the macromolecule, independently of its microscopic details. Our theoretical predictions are in excellent agreement with numerical stochastic simulations. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Journal ref: Phys. Rev. Lett. 115, 208301 (2015)

Diffusion and subdiffusion of interacting particles on comb-like structures

Authors: O. Bénichou, P. Illien, G. Oshanin, A. Sarracino, R. Voituriez

Abstract: We study the dynamics of a tracer particle (TP) on a comb lattice populated by randomly moving hard-core particles in the dense limit. We first consider the case where the TP is constrained to move on the backbone of the comb only, and, in the limit of high density of particles, we present exact analytical results for the cumulants of the TP position, showing a subdiffusive behavior… ▽ More We study the dynamics of a tracer particle (TP) on a comb lattice populated by randomly moving hard-core particles in the dense limit. We first consider the case where the TP is constrained to move on the backbone of the comb only, and, in the limit of high density of particles, we present exact analytical results for the cumulants of the TP position, showing a subdiffusive behavior $\sim t^{3/4}$. At longer times, a second regime is observed, where standard diffusion is recovered, with a surprising non analytical dependence of the diffusion coefficient on the particle density. When the TP is allowed to visit the teeth of the comb, based on a mean-field-like Continuous Time Random Walk description, we unveil a rich and complex scenario, with several successive subdiffusive regimes, resulting from the coupling between the inhomogeneous comb geometry and particle interactions. Remarkably, the presence of hard-core interactions speeds up the TP motion along the backbone of the structure in all regimes. △ Less

Submitted 2 July, 2015; originally announced July 2015.

Comments: 5 pages, 3 figures + supplemental material

Journal ref: Phys. Rev. Lett. 115, 220601 (2015)

Distribution of the position of a driven tracer in a hardcore lattice gas

Authors: Pierre Illien, Olivier Bénichou, Gleb Oshanin, Raphaël Voituriez

Abstract: We study the position of a biased tracer particle (TP) in a bath of hardcore particles moving on a lattice of arbitrary dimension and in contact with a reservoir. Starting from the master equation satisfied by the joint probability of the TP position and the bath configuration and resorting to a mean-field-type approximation, we presented a computation of the fluctuations of the TP position in a p… ▽ More We study the position of a biased tracer particle (TP) in a bath of hardcore particles moving on a lattice of arbitrary dimension and in contact with a reservoir. Starting from the master equation satisfied by the joint probability of the TP position and the bath configuration and resorting to a mean-field-type approximation, we presented a computation of the fluctuations of the TP position in a previous publication [Phys. Rev. E \textbf{87}, 032164 (2013)]. Counterintuitively, on a one-dimensional lattice, the diffusion coefficient of the TP was shown to be a non-monotonic function of the density of bath particles, and reaches a maximum for a nonzero value of the density. Here, we: (i) give the details of this computation and offer a physical insight into the understanding of the non-monotonicity of the diffusion coefficient; (ii) extend the mean-field-type approximation to decouple higher-order correlation functions, and obtain the evolution equation satisfied by the cumulant generating function of the position of the TP, valid in any space dimension. In the particular case of a one-dimensional lattice, we solve this equation and obtain the probability distribution of the TP position. We show that the position rescaled by its fluctuations is asymptotically distributed accordingly to a Gaussian distribution in the long-time limit. △ Less

Submitted 2 July, 2015; originally announced July 2015.

Comments: 53 pages, 12 figures

Microscopic theory for negative differential mobility in crowded environments

Authors: O. Bénichou, P. Illien, G. Oshanin, A. Sarracino, R. Voituriez

Abstract: We study the behavior of the stationary velocity of a driven particle in an environment of mobile hard-core obstacles. Based on a lattice gas model, we demonstrate analytically that the drift velocity can exhibit a nonmonotonic dependence on the applied force, and show quantitatively that such negative differential mobility (NDM), observed in various physical contexts, is controlled by both the de… ▽ More We study the behavior of the stationary velocity of a driven particle in an environment of mobile hard-core obstacles. Based on a lattice gas model, we demonstrate analytically that the drift velocity can exhibit a nonmonotonic dependence on the applied force, and show quantitatively that such negative differential mobility (NDM), observed in various physical contexts, is controlled by both the density and diffusion time scale of obstacles. Our study unifies recent numerical and analytical results obtained in specific regimes, and makes it possible to determine analytically the region of the full parameter space where NDM occurs. These results suggest that NDM could be a generic feature of biased (or active) transport in crowded environments. △ Less

Submitted 2 December, 2014; originally announced December 2014.

Comments: 5 pages, 2 figures + supplemental material

Journal ref: Phys. Rev. Lett. 113, 268002 (2014)

Exit time distribution in spherically symmetric two-dimensional domains

Authors: J. -F. Rupprecht, O. Bénichou, D. S. Grebenkov, R. Voituriez

Abstract: The distribution of exit times is computed for a Brownian particle in spherically symmetric two- dimensional domains (disks, angular sectors, annuli) and in rectangles that contain an exit on their boundary. The governing partial differential equation of Helmholtz type with mixed Dirichlet- Neumann boundary conditions is solved analytically. We propose both an exact solution relying on a matrix in… ▽ More The distribution of exit times is computed for a Brownian particle in spherically symmetric two- dimensional domains (disks, angular sectors, annuli) and in rectangles that contain an exit on their boundary. The governing partial differential equation of Helmholtz type with mixed Dirichlet- Neumann boundary conditions is solved analytically. We propose both an exact solution relying on a matrix inversion, and an approximate explicit solution. The approximate solution is shown to be exact for an exit of vanishing size and to be accurate even for large exits. For angular sectors, we also derive exact explicit formulas for the moments of the exit time. For annuli and rectangles, the approximate expression of the mean exit time is shown to be very accurate even for large exits. The analysis is also extended to biased diffusion. Since the Helmholtz equation with mixed boundary conditions is encountered in microfluidics, heat propagation, quantum billiards, and acoustics, the developed method can find numerous applications beyond exit processes. △ Less

Submitted 26 September, 2014; originally announced September 2014.

Comments: 34 pages, 12 figures in Journal of Statistical Physics, 2014

Mean cover time of one-dimensional persistent random walks

Authors: Marie Chupeau, Olivier Bénichou, Raphaël Voituriez

Abstract: The cover time is defined as the time needed for a random walker to visit every site of a confined domain. Here, we focus on persistent random walks, which provide a minimal model of random walks with short range memory. We derive the exact expression of the mean cover time of a one-dimensional lattice by such a persistent random walk, both for periodic and reflecting boundary conditions. The cover time is defined as the time needed for a random walker to visit every site of a confined domain. Here, we focus on persistent random walks, which provide a minimal model of random walks with short range memory. We derive the exact expression of the mean cover time of a one-dimensional lattice by such a persistent random walk, both for periodic and reflecting boundary conditions. △ Less

Submitted 6 June, 2014; originally announced June 2014.

Comments: To appear in Phys. Rev. E

Geometry-induced superdiffusion in driven crowded systems

Authors: Olivier Bénichou, Anna Bodrova, Dipanjan Chakraborty, Pierre Illien, Adam Law, Carlos Mejía-Monasterio, Gleb Oshanin, Raphaël Voituriez

Abstract: Recent Molecular Dynamics simulations of glass-forming liquids revealed superdiffusive fluctuations associated with the position of a tracer particle (TP) driven by an external force. Such anomalous response, whose mechanism remains elusive, has been observed up to now only in systems close to their glass transition, suggesting that this could be one of its hallmarks. Here, we show that the presen… ▽ More Recent Molecular Dynamics simulations of glass-forming liquids revealed superdiffusive fluctuations associated with the position of a tracer particle (TP) driven by an external force. Such anomalous response, whose mechanism remains elusive, has been observed up to now only in systems close to their glass transition, suggesting that this could be one of its hallmarks. Here, we show that the presence of superdiffusion is in actual fact much more general, provided that the system is crowded and geometrically confined. We present and solve analytically a minimal model consisting of a driven TP in a dense, crowded medium in which the motion of particles is mediated by the diffusion of packing defects, called vacancies. For such non glass-forming systems, our analysis predicts a long-lived superdiffusion which ultimately crosses over to giant diffusive behavior. We find that this trait is present in confined geometries, for example long capillaries and stripes, and emerges as a universal response of crowded environments to an external force. These findings are confirmed by numerical simulations of systems as varied as lattice-gases, dense liquids and granular fluids. △ Less

Submitted 30 October, 2013; originally announced October 2013.

Geometric friction directs cell migration

Authors: M. Le Berre, Yan-Jun Liu, J. Hu, P. Maiuri, O. Bénichou, R. Voituriez, Y. Chen, M. Piel

Abstract: In the absence of environmental cues, a migrating cell performs an isotropic random motion. Recently, the breaking of this isotropy has been observed when cells move in the presence of asymmetric adhesive patterns. However, up to now the mechanisms at work to direct cell migration in such environments remain unknown. Here, we show that a non-adhesive surface with asymmetric micro-geometry consisti… ▽ More In the absence of environmental cues, a migrating cell performs an isotropic random motion. Recently, the breaking of this isotropy has been observed when cells move in the presence of asymmetric adhesive patterns. However, up to now the mechanisms at work to direct cell migration in such environments remain unknown. Here, we show that a non-adhesive surface with asymmetric micro-geometry consisting of dense arrays of tilted micro-pillars can direct cell motion. Our analysis reveals that most features of cell trajectories, including the bias, can be reproduced by a simple model of active Brownian particle in a ratchet potential, which we suggest originates from a generic elastic interaction of the cell body with the environment. The observed guiding effect, independent of adhesion, is therefore robust and could be used to direct cell migration both in vitro and in vivo. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: To appear in PRL. Supplementary Information file available on request

Active polar fluid flow in finite droplets

Authors: Carl A. Whitfield, Davide Marenduzzo, Raphaël Voituriez, Rhoda J. Hawkins

Abstract: We present a continuum level analytical model of a droplet of active contractile fluid consisting of filaments and motors. We calculate the steady state flows that result from a splayed polarisation of the filaments. We account for the interaction with an arbitrary external medium by imposing a viscous friction at the fixed droplet boundary. We then show that the droplet has non-zero force dipole… ▽ More We present a continuum level analytical model of a droplet of active contractile fluid consisting of filaments and motors. We calculate the steady state flows that result from a splayed polarisation of the filaments. We account for the interaction with an arbitrary external medium by imposing a viscous friction at the fixed droplet boundary. We then show that the droplet has non-zero force dipole and quadrupole moments, the latter of which is essential for self-propelled motion of the droplet at low Reynolds' number. Therefore, this calculation describes a simple mechanism for the motility of a droplet of active contractile fluid embedded in a 3D environment, which is relevant to cell migration in confinement (for example, embedded within a gel or tissue). Our analytical results predict how the system depends on various parameters such as the effective friction coefficient, the phenomenological activity parameter and the splay of the imposed polarisation. △ Less

Submitted 19 February, 2014; v1 submitted 12 July, 2013; originally announced July 2013.

Comments: 9 pages, 5 figures + 4 appendices (15 pages total)

MSC Class: 92C17

Journal ref: Eur. Phys. J. E, 37 2 (2014) 8