-
Stationary particle currents in sedimenting active matter wetting a wall
Authors:
Matthieu Mangeat,
Shauri Chakraborty,
Adam Wysocki,
Heiko Rieger
Abstract:
Recently it was predicted, on the basis of a lattice gas model, that scalar active matter in a gravitational field would rise against gravity up a confining wall or inside a thin capillary - in spite of repulsive particle-wall interactions [Phys. Rev. Lett. 124, 048001 (2020)]. In this paper we confirm this prediction with sedimenting active Brownian particles (ABPs) in a box numerically and eluci…
▽ More
Recently it was predicted, on the basis of a lattice gas model, that scalar active matter in a gravitational field would rise against gravity up a confining wall or inside a thin capillary - in spite of repulsive particle-wall interactions [Phys. Rev. Lett. 124, 048001 (2020)]. In this paper we confirm this prediction with sedimenting active Brownian particles (ABPs) in a box numerically and elucidate the mechanism leading to the formation of a meniscus rising above the bulk of the sedimentation region. The height of the meniscus increases with the activity of the system, algebraically with the Péclet number. The formation of the meniscus is determined by a stationary circular particle current, a vortex, centered at the base of the meniscus, whose size and strength increase with the ABP activity. The origin of these vortices can be traced back to the confinement of the ABPs in a box: already the stationary state of ideal (non-interacting) ABPs without gravitation displays circular currents that arrange in a highly symmetric way in the eight octants of the box. Gravitation distorts this vortex configuration downward, leaving two major vortices at the two side walls, with a strong downward flow along the walls. Repulsive interactions between the ABPs change this situation only as soon as motility induced phase separation (MIPS) sets in and forms a dense, sedimented liquid region at the bottom, which pushes the center of the vortex upwards towards the liquid-gas interface. Self-propelled particles therefore represent an impressive realization of scalar active matter that forms stationary particle currents being able to perform visible work against gravity or any other external field, which we predict to be observable experimentally in active colloids under gravitation.
△ Less
Submitted 15 January, 2024; v1 submitted 18 September, 2023;
originally announced September 2023.
-
Jamming and flocking in the restricted active Potts model
Authors:
Mintu Karmakar,
Swarnajit Chatterjee,
Matthieu Mangeat,
Heiko Rieger,
Raja Paul
Abstract:
We study the active Potts model with either site occupancy restriction or on-site repulsion to explore jamming and kinetic arrest in a flocking model. The incorporation of such volume exclusion features leads to a surprisingly rich variety of self-organized spatial patterns. While bands and lanes of moving particles commonly occur without or under weak volume exclusion, strong volume exclusion alo…
▽ More
We study the active Potts model with either site occupancy restriction or on-site repulsion to explore jamming and kinetic arrest in a flocking model. The incorporation of such volume exclusion features leads to a surprisingly rich variety of self-organized spatial patterns. While bands and lanes of moving particles commonly occur without or under weak volume exclusion, strong volume exclusion along with low temperature, high activity, and large particle density facilitates jams due to motility-induced phase separation. Through several phase diagrams, we identify the phase boundaries separating the jammed and free-flowing phases and study the transition between these phases which provide us with both qualitative and quantitative predictions of how jamming might be delayed or dissolved. We further formulate and analyze a hydrodynamic theory for the restricted APM which predicts various features of the microscopic model.
△ Less
Submitted 18 July, 2023; v1 submitted 20 December, 2022;
originally announced December 2022.
-
Flocking of two unfriendly species: The two-species Vicsek model
Authors:
Swarnajit Chatterjee,
Matthieu Mangeat,
Chul-Ung Woo,
Heiko Rieger,
Jae Dong Noh
Abstract:
We consider the two-species Vicsek model (TSVM) consisting of two kinds of self-propelled particles, A and B, that tend to align with particles from the same species and to antialign with the other. The model shows a flocking transition that is reminiscent of the original Vicsek model: it has a liquid-gas phase transition and displays micro-phase-separation in the coexistence region where multiple…
▽ More
We consider the two-species Vicsek model (TSVM) consisting of two kinds of self-propelled particles, A and B, that tend to align with particles from the same species and to antialign with the other. The model shows a flocking transition that is reminiscent of the original Vicsek model: it has a liquid-gas phase transition and displays micro-phase-separation in the coexistence region where multiple dense liquid bands propagate in a gaseous background. The interesting features of the TSVM are the existence of two kinds of bands, one composed of mainly A particles and one mainly of B particles, the appearance of two dynamical states in the coexistence region: the PF (parallel flocking) state in which all bands of the two species propagate in the same direction, and the APF (antiparallel flocking) state in which the bands of species A and species B move in opposite directions. When PF and APF states exist in the low-density part of the coexistence region they perform stochastic transitions from one to the other. The system size dependence of the transition frequency and dwell times show a pronounced crossover that is determined by the ratio of the band width and the longitudinal system size. Our work paves the way for studying multispecies flocking models with heterogeneous alignment interactions.
△ Less
Submitted 15 February, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Polar flocks with discretized directions: the active clock model approaching the Vicsek model
Authors:
Swarnajit Chatterjee,
Matthieu Mangeat,
Heiko Rieger
Abstract:
We consider the off-lattice two-dimensional $q$-state active clock model (ACM) as a natural discretization of the Vicsek model (VM) describing flocking. The ACM consists of particles able to move in the plane in a discrete set of $q$ equidistant angular directions, as in the active Potts model (APM), with an alignment interaction inspired by the ferromagnetic equilibrium clock model. We find that…
▽ More
We consider the off-lattice two-dimensional $q$-state active clock model (ACM) as a natural discretization of the Vicsek model (VM) describing flocking. The ACM consists of particles able to move in the plane in a discrete set of $q$ equidistant angular directions, as in the active Potts model (APM), with an alignment interaction inspired by the ferromagnetic equilibrium clock model. We find that for a small number of directions, the flocking transition of the ACM has the same phenomenology as the APM, including macrophase separation and reorientation transition. For a larger number of directions, the flocking transition in the ACM becomes equivalent to the one of the VM and displays microphase separation and only transverse bands, i.e. no re-orientation transition. Concomitantly also the transition of the $q\to\infty$ limit of the ACM, the active XY model (AXYM), is in the same universality class as the VM. We also construct a coarse-grained hydrodynamic description for the ACM and AXYM akin to the VM.
△ Less
Submitted 9 May, 2022; v1 submitted 2 March, 2022;
originally announced March 2022.
-
How stickiness can speed up diffusion in confined systems
Authors:
Arthur Alexandre,
Matthieu Mangeat,
Thomas Guérin,
David S. Dean
Abstract:
The paradigmatic model for heterogeneous media used in diffusion studies is built from reflecting obstacles and surfaces. It is well known that the crowding effect produced by these reflecting surfaces slows the dispersion of Brownian tracers. Here, using a general adsorption desorption model with surface diffusion, we show analytically that making surfaces or obstacles attractive can accelerate d…
▽ More
The paradigmatic model for heterogeneous media used in diffusion studies is built from reflecting obstacles and surfaces. It is well known that the crowding effect produced by these reflecting surfaces slows the dispersion of Brownian tracers. Here, using a general adsorption desorption model with surface diffusion, we show analytically that making surfaces or obstacles attractive can accelerate dispersion. In particular, we show that this enhancement of diffusion can exist even when the surface diffusion constant is smaller than that in the bulk. Even more remarkably, this enhancement effect occurs when the effective diffusion constant, when restricted to surfaces only, is lower than the effective diffusivity with purely reflecting boundaries. We give analytical formulas for this intriguing effect in periodic arrays of spheres as well as undulating micro-channels. Our results are confirmed by numerical calculations and Monte Carlo simulations.
△ Less
Submitted 2 June, 2022; v1 submitted 10 December, 2021;
originally announced December 2021.
-
Effective diffusivity of Brownian particles in a two dimensional square lattice of hard disks
Authors:
M. Mangeat,
T. Guérin,
D. S. Dean
Abstract:
We revisit the classic problem of the effective diffusion constant of a Brownian particle in a square lattice of reflecting impenetrable hard disks. This diffusion constant is also related to the effective conductivity of non-conducting and infinitely conductive disks in the same geometry. We show how a recently derived Green's function for the periodic lattice can be exploited to derive a series…
▽ More
We revisit the classic problem of the effective diffusion constant of a Brownian particle in a square lattice of reflecting impenetrable hard disks. This diffusion constant is also related to the effective conductivity of non-conducting and infinitely conductive disks in the same geometry. We show how a recently derived Green's function for the periodic lattice can be exploited to derive a series expansion of the diffusion constant in terms of the disk's volume fraction $\varphi$. Secondly we propose a variant of the Fick-Jacobs approximation to study the large volume fraction limit. This combination of analytical results is shown to describe the behavior of the diffusion constant for all volume fractions.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
Steady state of overdamped particles in the non-conservative force field of a simple non-linear model of optical trap
Authors:
M. Mangeat,
T. Guérin,
D. S. Dean
Abstract:
Optically trapped particles are often subject to a non-conservative scattering force arising from radiation pressure. In this paper we present an exact solution for the steady state statistics of an overdamped Brownian particle subjected to a commonly used force field model for an optical trap. The model is the simplest of its kind that takes into account non-conservative forces. In particular, we…
▽ More
Optically trapped particles are often subject to a non-conservative scattering force arising from radiation pressure. In this paper we present an exact solution for the steady state statistics of an overdamped Brownian particle subjected to a commonly used force field model for an optical trap. The model is the simplest of its kind that takes into account non-conservative forces. In particular, we present exact results for certain marginals of the full three dimensional steady state probability distribution as well as results for the toroidal probability currents which are present in the steady state, as well as for the circulation of theses currents. Our analytical results are confirmed by numerical solution of the steady state Fokker-Planck equation.
△ Less
Submitted 8 October, 2021;
originally announced October 2021.
-
Narrow escape problem in two-shell spherical domains
Authors:
Matthieu Mangeat,
Heiko Rieger
Abstract:
Intracellular transport in living cells is often spatially inhomogeneous with an accelerated effective diffusion close to the cell membrane and a ballistic motion away from the centrosome due to active transport along actin filaments and microtubules, respectively. Recently it was reported that the mean first passage time (MFPT) for transport to a specific area on the cell membrane is minimal for…
▽ More
Intracellular transport in living cells is often spatially inhomogeneous with an accelerated effective diffusion close to the cell membrane and a ballistic motion away from the centrosome due to active transport along actin filaments and microtubules, respectively. Recently it was reported that the mean first passage time (MFPT) for transport to a specific area on the cell membrane is minimal for an optimal actin cortex width. In this paper we ask whether this optimization in a two-compartment domain can also be achieved by passive Brownian particles. We consider a Brownian motion with different diffusion constants in the two shells and a potential barrier between the two and investigate the narrow escape problem by calculating the MFPT for Brownian particles to reach a small window on the external boundary. In two and three dimensions, we derive asymptotic expressions for the MFPT in the thin cortex and small escape region limits confirmed by numerical calculations of the MFPT using the finite element method and stochastic simulations. From this analytical and numeric analysis we finally extract the dependence of the MFPT on the ratio of diffusion constants, the potential barrier height and the width of the outer shell. The first two are monotonous whereas the last one may have a minimum for a sufficiently attractive cortex, for which we propose an analytical expression of the potential barrier height matching very well the numerical predictions.
△ Less
Submitted 6 October, 2021; v1 submitted 27 April, 2021;
originally announced April 2021.
-
Flocking with a $q$-fold discrete symmetry: band-to-lane transition in the active Potts model
Authors:
Matthieu Mangeat,
Swarnajit Chatterjee,
Raja Paul,
Heiko Rieger
Abstract:
We study the $q$-state active Potts model (APM) on a two-dimensional lattice in which self-propelled particles have $q$ internal states corresponding to the $q$ directions of motion. A local alignment rule inspired by the ferromagnetic $q$-state Potts model and self-propulsion via biased diffusion according to the internal particle states leads to a collective motion at high densities and low nois…
▽ More
We study the $q$-state active Potts model (APM) on a two-dimensional lattice in which self-propelled particles have $q$ internal states corresponding to the $q$ directions of motion. A local alignment rule inspired by the ferromagnetic $q$-state Potts model and self-propulsion via biased diffusion according to the internal particle states leads to a collective motion at high densities and low noise. We formulate a coarse-grained hydrodynamic theory of the model, compute the phase diagrams of the APM for $q=4$ and $q=6$ and explore the flocking dynamics in the region, in which the high-density (polar liquid) phase coexists with the low-density (gas) phase and forms a fluctuating stripe of coherently moving particles. As a function of the particle self-propulsion velocity, a novel reorientation transition of the phase-separated profiles from transversal band motion to longitudinal lane formation is revealed, which is absent in the Vicsek model and APM for $q=2$ (active Ising model). The origin of this reorientation transition is obtained via a stability analysis: for large velocities, the transverse diffusion constant approaches zero and stabilizes lanes. Finally, we perform microscopic simulations that corroborate our analytical predictions about the flocking and reorientation transitions and validate the phase diagrams of the APM.
△ Less
Submitted 29 July, 2020;
originally announced July 2020.
-
Flocking and reorientation transition in the 4-state active Potts model
Authors:
Swarnajit Chatterjee,
Matthieu Mangeat,
Raja Paul,
Heiko Rieger
Abstract:
We study the active 4-state Potts model (APM) on the square lattice in which active particles have four internal states corresponding to the four directions of motion. A local alignment rule inspired by the ferromagnetic 4-state Potts model and self-propulsion via biased diffusion according to the internal particle states leads to flocking at high densities and low noise. We compute the phase diag…
▽ More
We study the active 4-state Potts model (APM) on the square lattice in which active particles have four internal states corresponding to the four directions of motion. A local alignment rule inspired by the ferromagnetic 4-state Potts model and self-propulsion via biased diffusion according to the internal particle states leads to flocking at high densities and low noise. We compute the phase diagram of the APM and explore the flocking dynamics in the region, in which the high-density (liquid) phase coexists with the low-density (gas) phase and forms a fluctuating band of coherently moving particles. As a function of the particle self-propulsion velocity, a novel reorientation transition of the phase-separated profiles from transversal to longitudinal band motion is revealed, which is absent in the Vicsek model and the active Ising model. We further construct a coarse-grained hydrodynamic description of the model which validates the results for the microscopic model.
△ Less
Submitted 2 March, 2020; v1 submitted 29 November, 2019;
originally announced November 2019.
-
The narrow escape problem in a circular domain with radial piecewise constant diffusivity
Authors:
Matthieu Mangeat,
Heiko Rieger
Abstract:
The stochastic motion of particles in living cells is often spatially inhomogeneous with a higher effective diffusivity in a region close to the cell boundary due to active transport along actin filaments. As a first step to understand the consequence of the existence of two compartments with different diffusion constant for stochastic search problems we consider here a Brownian particle in a circ…
▽ More
The stochastic motion of particles in living cells is often spatially inhomogeneous with a higher effective diffusivity in a region close to the cell boundary due to active transport along actin filaments. As a first step to understand the consequence of the existence of two compartments with different diffusion constant for stochastic search problems we consider here a Brownian particle in a circular domain with different diffusion constants in the inner and the outer shell. We focus on the narrow escape problem and compute the mean first passage time (MFPT) for Brownian particles starting at some pre-defined position to find a small region on the outer reflecting boundary. For the annulus geometry we find that the MFPT can be minimized for a specific value of the width of the outer shell. In contrast for the two-shell geometry we show that the MFPT depends monotonously on all model parameters, in particular on the outer shell width. Moreover we find that the distance between the starting point and the narrow escape region which maximizes the MFPT depends discontinuously on the ratio between inner and outer diffusivity.
△ Less
Submitted 26 August, 2019; v1 submitted 17 June, 2019;
originally announced June 2019.
-
The role of non-conservative scattering forces and damping on Brownian particles in optical traps
Authors:
Matthieu Mangeat,
Yacine Amarouchene,
Yann Louyer,
Thomas Guérin,
David S. Dean
Abstract:
We consider a model of a particle trapped in a harmonic optical trap but with the addition of a non-conservative radiation induced force. This model is known to correctly describe experimentally observed trapped particle statistics for a wide range of physical parameters such as temperature and pressure. We theoretically analyse the effect of non-conservative force on the underlying steady state d…
▽ More
We consider a model of a particle trapped in a harmonic optical trap but with the addition of a non-conservative radiation induced force. This model is known to correctly describe experimentally observed trapped particle statistics for a wide range of physical parameters such as temperature and pressure. We theoretically analyse the effect of non-conservative force on the underlying steady state distribution as well as the power spectrum for the particle position. We compute perturbatively the probability distribution of the resulting non-equilibrium steady states for all dynamical regimes, underdamped through to overdamped and give expressions for the associated currents in phase space (position and velocity). We also give the spectral density of the trapped particle's position in all dynamical regimes and for any value of the non-conservative force. Signatures of the presence non-conservative forces are shown to be particularly strong for in the underdamped regime at low frequencies.
△ Less
Submitted 21 December, 2018;
originally announced December 2018.
-
Nonequilibrium dynamics induced by scattering forces for optically trapped nanoparticles in strongly inertial regimes
Authors:
Yacine Amarouchene,
Matthieu Mangeat,
Benjamin Vidal Montes,
Lukas Ondic,
Thomas Guérin,
David. S. Dean,
Yann Louyer
Abstract:
As a first approximation, the forces acting on optically trapped particles are commonly assumed to be conservative. The influence of the nonconservative force has been shown to be negligible in overdamped liquid environments. However, its impact in the underdamped regime remains unexplored. Here, we experimentally study the combined effects of gradient and scattering forces on the nonlinear dynami…
▽ More
As a first approximation, the forces acting on optically trapped particles are commonly assumed to be conservative. The influence of the nonconservative force has been shown to be negligible in overdamped liquid environments. However, its impact in the underdamped regime remains unexplored. Here, we experimentally study the combined effects of gradient and scattering forces on the nonlinear dynamics of trapped particles at various air pressures, leading to significant low-frequency noise excess along the optical axis and the emergence of toroidal Brownian vortices. Our results supported, by numerical simulations and a theoretical model, provide fundamental insights into the functioning of optical tweezers and a means for investigating nonequilibrium steady states in the presence of nonconservative forces.
△ Less
Submitted 9 April, 2019; v1 submitted 14 December, 2018;
originally announced December 2018.
-
Dispersion in two-dimensional periodic channels with discontinuous profiles
Authors:
M. Mangeat,
T. Guérin,
D. S. Dean
Abstract:
The effective diffusivity of Brownian tracer particles confined in periodic micro-channels is smaller than the microscopic diffusivity due to entropic trapping. Here, we study diffusion in two-dimensional periodic channels whose cross-section presents singular points, such as abrupt changes of radius or the presence of thin walls, with openings, delimiting periodic compartments composing the chann…
▽ More
The effective diffusivity of Brownian tracer particles confined in periodic micro-channels is smaller than the microscopic diffusivity due to entropic trapping. Here, we study diffusion in two-dimensional periodic channels whose cross-section presents singular points, such as abrupt changes of radius or the presence of thin walls, with openings, delimiting periodic compartments composing the channel. Dispersion in such systems is analyzed using the Fick-Jacobs' approximation. This approximation assumes a much faster equilibration in the lateral than in the axial direction, along which the dispersion is measured. If the characteristic width $a$ of the channel is much smaller than the period $L$ of the channel, i.e. $\varepsilon = a/L$ is small, this assumption is clearly valid for Brownian particles. For discontinuous channels, the Fick-Jacobs' approximation is only valid at the lowest order in $\varepsilon$ and provides a rough, though on occasions rather accurate, estimate of the effective diffusivity. Here we provide formulas for the effective diffusivity in discontinuous channels that are asymptotically exact at the next-to-leading order in $\varepsilon$. Each discontinuity leads to a reduction of the effective diffusivity. We show that our theory is consistent with the picture of effective {\em trapping rates} associated with each discontinuity, for which our theory provides explicit and asymptotically exact formulas. Our analytical predictions are confirmed by numerical analysis. Our results provide a precise quantification of the kinetic entropic barriers associated with profile singularities.
△ Less
Submitted 14 July, 2018;
originally announced July 2018.
-
Dispersion in two dimensional channels - the Fick-Jacobs approximation revisited
Authors:
M. Mangeat,
T. Guérin,
D. S. Dean
Abstract:
We examine the dispersion of Brownian particles in a symmetric two dimensional channel, this classical problem has been widely studied in the literature using the so called Fick-Jacobs' approximation and its various improvements. Most studies rely on the reduction to an effective one dimensional diffusion equation, here we drive an explicit formula for the diffusion constant which avoids this redu…
▽ More
We examine the dispersion of Brownian particles in a symmetric two dimensional channel, this classical problem has been widely studied in the literature using the so called Fick-Jacobs' approximation and its various improvements. Most studies rely on the reduction to an effective one dimensional diffusion equation, here we drive an explicit formula for the diffusion constant which avoids this reduction. Using this formula the effective diffusion constant can be evaluated numerically without resorting to Brownian simulations. In addition a perturbation theory can be developed in $\varepsilon = h_0/L$ where $h_0$ is the characteristic channel height and $L$ the period. This perturbation theory confirms the results of Kalinay and Percus (Phys. Rev. E 74, 041203 (2006)), based on the reduction, to one dimensional diffusion are exact at least to ${\cal O}(\varepsilon^6)$. Furthermore, we show how the Kalinay and Percus pseudo-linear approximation can be straightforwardly recovered. The approach proposed here can also be exploited to yield exact results an appropriate limit $\varepsilon \to \infty$, we show that here the diffusion constant remains finite and show how the result can be obtained with a simple physical argument. Moreover we show that the correction to the effective diffusion constant is of order $1/\varepsilon$ and remarkably has a some universal characteristics. Numerically we compare the analytical results obtained with exact numerical calculations for a number of interesting channel geometries.
△ Less
Submitted 8 November, 2017; v1 submitted 7 October, 2017;
originally announced October 2017.
-
Geometry controlled dispersion in periodic corrugated channels
Authors:
M. Mangeat,
T. Guérin,
D. S. Dean
Abstract:
The effective diffusivity $D_e$ of tracer particles diffusing in periodically corrugated axisymmetric two and three dimensional channels is studied. The majority of previous studies of this class of problems are based on perturbative analyses about narrow channels, where the problem can be reduced to an effectively one dimensional one. Here we show how to analyze this class of problems using a muc…
▽ More
The effective diffusivity $D_e$ of tracer particles diffusing in periodically corrugated axisymmetric two and three dimensional channels is studied. The majority of previous studies of this class of problems are based on perturbative analyses about narrow channels, where the problem can be reduced to an effectively one dimensional one. Here we show how to analyze this class of problems using a much more general approach which even includes the limit of infinitely wide channels. Using the narrow and wide channel asymptotics, we provide a Pad{é} approximant scheme that is able to describe the dispersion properties of a wide class of channels. Furthermore, we systematically identify all the exact asymptotic scaling regimes of $D_e$ and the accompanying physical mechanisms that control dispersion, clarifying the distinction between smooth channels and compartmentalized ones, and identifying the regimes in which $D_e$ can be linked to first passage problems.
△ Less
Submitted 12 September, 2017;
originally announced September 2017.
-
Quantitative approximation schemes for glasses
Authors:
Matthieu Mangeat,
Francesco Zamponi
Abstract:
By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode-Coupling appro…
▽ More
By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode-Coupling approximation scheme. Our scheme becomes exact, by construction, in dimension $d\to\infty$ and it can be improved systematically by adding more terms in the expansion.
△ Less
Submitted 22 January, 2016; v1 submitted 13 October, 2015;
originally announced October 2015.