Effective diffusion of confined active Brownian swimmers

Turning catalytically active pores into active pumps

We develop a semi-analytical model of self-diffusioosmotic transport in active pores, which includes advective transport and the inverse chemical reaction that consumes solute. In previous work [Antunes et al., Phys. Rev. Lett. 129, 188003 (2022)], we have demonstrated the existence of a spontaneous symmetry breaking in fore-aft symmetric pores that enables them to function as a micropump. We now show that this pumping transition is controlled by three timescales. Two timescales characterize advective and diffusive transport. The third timescale corresponds to how long a solute molecule resides in the pore before being consumed. Introducing asymmetry to the pore (either via the shape or the catalytic coating) reveals a second type of advection-enabled transition. In asymmetric pores, the flow rate exhibits discontinuous jumps and hysteresis loops upon tuning the parameters that control the asymmetry. This work demonstrates the interconnected roles of shape and catalytic patterning in the dynamics of active pores and shows how to design a pump for optimum performance.

Ratchet effect in an asymmetric two-dimensional system of Janus particles

We consider a disk-like Janus particle self-driven by a force of constant magnitude f, but an arbitrary direction depending on the stochastic rotation of the disk. The particle diffuses in a two-dimensional channel of varying width 2h(x). We applied the procedure mapping the 2+1-dimensional Fokker-Planck equation onto the longitudinal coordinate x; the result is the Fick-Jacobs equation extended by the spatially dependent effective diffusion constant D(x) and an additional effective potential -γ(x), derived recursively within the mapping procedure. Unlike the entropic potential ∼lnh(x), γ(x) becomes an increasing or decreasing function also in periodic channels, depending on the asymmetry of h(x) and thus it visualizes the net force driving the ratchet current. We demonstrate the appearance of the ratchet effect on a trial asymmetric channel; our theory is verified by a numerical solution of the corresponding Fokker-Planck equation. Isotropic driving force f results in the monotonic decrease of the ratchet current with a growing ratio α=D_{R}/D_{T} of the rotation and the translation diffusion constants; asymptotically going ∼1/α^{2}. If we allow anisotropy of the force, we can observe the current reversal depending on α.

Boundary design regulates the diffusion of active matter in heterogeneous environments

Physical boundaries play a key role in governing the overall transport properties of nearby self-propelled particles. In this work, we develop dispersion theories and conduct Brownian dynamics simulations to predict the coupling between surface accumulation and effective diffusivity of active particles in boundary-rich media. We focus on three models that are well-understood for passive systems: particle transport in (i) an array of fixed volume-excluding obstacles; (ii) a pore with spatially heterogeneous width; and (iii) a tortuous path with kinks and corners. While the impact of these entropic barriers on passive particle transport is well established, we find that these classical models of porous media flows break down due to the unique interplay between activity and the microstructure of the internal geometry. We study the activity-induced slowdown of effective diffusivity by formulating a Smoluchowski description of long-time self diffusivity which contains contributions from the density and fluctuation fields of the active particles. Particle-based and finite element simulations corroborate this perspective and reveal important nonequilibrium considerations of active transport.

Closed Formula for Transport across Constrictions

In the last decade, the Fick-Jacobs approximation has been exploited to capture transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that, for fore-aft symmetric channels, the flux of noninteracting systems is fully captured by its linear response regime. For this case, we derive a very simple formula that captures the correct trends and can be exploited as a simple tool to design experiments or simulations. Lastly, we show that higher-order corrections in the flux may appear for nonsymmetric channels.

Fick-Jacobs description and first passage dynamics for diffusion in a channel under stochastic resetting

The transport of particles through channels is of paramount importance in physics, chemistry, and surface science due to its broad real world applications. Much insight can be gained by observing the transition paths of a particle through a channel and collecting statistics on the lifetimes in the channel or the escape probabilities from the channel. In this paper, we consider the diffusive transport through a narrow conical channel of a Brownian particle subject to intermittent dynamics, namely, stochastic resetting. As such, resetting brings the particle back to a desired location from where it resumes its diffusive phase. To this end, we extend the Fick-Jacobs theory of channel-facilitated diffusive transport to resetting-induced transport. Exact expressions for the conditional mean first passage times, escape probabilities, and the total average lifetime in the channel are obtained, and their behavior as a function of the resetting rate is highlighted. It is shown that resetting can expedite the transport through the channel-rigorous constraints for such conditions are then illustrated. Furthermore, we observe that a carefully chosen resetting rate can render the average lifetime of the particle inside the channel minimal. Interestingly, the optimal rate undergoes continuous and discontinuous transitions as some relevant system parameters are varied. The validity of our one-dimensional analysis and the corresponding theoretical predictions is supported by three-dimensional Brownian dynamics simulations. We thus believe that resetting can be useful to facilitate particle transport across biological membranes-a phenomenon that can spearhead further theoretical and experimental studies.

First-passage times in conical varying-width channels biased by a transverse gravitational force: Comparison of analytical and numerical results

We study the crossing time statistic of diffusing point particles between the two ends of expanding and narrowing two-dimensional conical channels under a transverse external gravitational field. The theoretical expression for the mean first-passage time for such a system is derived under the assumption that the axial diffusion in a two-dimensional channel of smoothly varying geometry can be approximately described as a one-dimensional diffusion in an entropic potential with position-dependent effective diffusivity in terms of the modified Fick-Jacobs equation. We analyze the channel crossing dynamics in terms of the mean first-passage time, combining our analytical results with extensive two-dimensional Brownian dynamics simulations, allowing us to find the range of applicability of the one-dimensional approximation. We find that the effective particle diffusivity decreases with increasing amplitude of the external potential. Remarkably, the mean first-passage time for crossing the channel is shown to assume a minimum at finite values of the potential amplitude.

Transverse dichotomic ratchet in a two-dimensional corrugated channel

A particle diffusing in a two-dimensional channel of varying width h(x) is considered. It is driven by a force of constant magnitude f, but random orientation across the channel. We suggest the projection technique to study the ratchet effect appearing in this system. Reducing the transverse coordinate, as well as the orientation of the force in the full-dimensional Fokker-Planck equation, we arrive at the generalized Fick-Jacobs equation, describing dynamics of the system in the longitudinal coordinate x only. The additional effective potential -γ(x), calculated within the mapping procedure, exhibits an increasing or decreasing part in the channel shaped by an asymmetric periodic h(x), which determines the appearing ratchet current. As shown on a specific example, random driving in the transverse direction is much more effective than that in the longitudinal direction, at least for quickly flipping orientation of the force. Also, the transverse and the longitudinal driving push the rectified current in opposite directions along the same channel.

Dichotomic ratchet in a two-dimensional corrugated channel

We consider a particle diffusing in a two-dimensional (2D) channel of varying width h(x). It is driven by a force of constant magnitude f but random orientation there or back along the channel. We derive the effective generalized Fick-Jacobs equation for this system, which describes the dynamics of such a particle in the longitudinal coordinate x. Aside from the effective diffusion coefficient D(x), our mapping also generates an additional effective potential -γ(x) added to the entropic potential -log[h(x)]. It acquires an increasing or decreasing component in asymmetric periodic channels, and thus it explains appearance of the ratchet current. We study this effect on a trial example and compare the results of our true 2D theory with a commonly used effective one-dimensional description; the data are verified by the numerical solution of the full 2D problem.

Two-dimensional diffusion biased by a transverse gravitational force in an asymmetric channel: Reduction to an effective one-dimensional description

We focus on the derivation of a general position-dependent effective diffusion coefficient to describe two-dimensional (2D) diffusion in a narrow and smoothly asymmetric channel of varying width under a transverse gravitational external field, a generalization of the symmetric channel case using the projection method introduced earlier by Kalinay and Percus [P. Kalinay and J. K. Percus, J. Chem. Phys. 122, 204701 (2005)10.1063/1.1899150]. To this end, we project the 2D Smoluchowski equation into an effective one-dimensional generalized Fick-Jacobs equation in the presence of constant force in the transverse direction. The expression for the diffusion coefficient given in Eq. (34) is our main result. This expression is a more general effective diffusion coefficient for narrow 2D channels in the presence of constant transverse force, which contains the well-known previous results for a symmetric channel obtained by Kalinay, as well as the limiting cases when the transverse gravitational external field goes to zero and infinity. Finally, we show that diffusivity can be described by the interpolation formula proposed by Kalinay, D_{0}/[1+(1/4)w^{'2}(x)]^{-η}, where spatial confinement, asymmetry, and the presence of a constant transverse force can be encoded in η, which is a function of channel width (w), channel centerline, and transverse force. The interpolation formula also reduces to well-known previous results, namely, those obtained by Reguera and Rubi [D. Reguera and J. M. Rubi, Phys. Rev. E 64, 061106 (2001)10.1103/PhysRevE.64.061106] and by Kalinay [P. Kalinay, Phys. Rev. E 84, 011118 (2011)10.1103/PhysRevE.84.011118].

Reduced dynamics of a one-dimensional Janus particle

A Janus particle diffusing on a line is considered. Aside from its own driving force f acting forward or backward according to its stochastic orientation, it moves in a position-dependent potential U(x). We propose here the mapping scheme generating the effective generalized Fick-Jacobs equation, describing motion of the particle in the spatial coordinate x only; the orientation is understood as the transverse coordinate. The self-propulsion, driving the system out of equilibrium, is reflected as an additional effective potential -γ(x) in the reduced picture. It enables us to understand peculiarities of this system in a handy way. The additionally corrected potential redistributes the confined particles in quasiequilibrium causing their piling at the walls. In periodic asymmetric channels, it acquires a growing contribution, responsible for driving the ratchet effect.

Local pressure for confined systems

We derive a general closed expression for the local pressure exerted onto the corrugated walls of a channel confining a fluid medium. When the fluid medium is at equilibrium, the local pressure is a functional of the shape of the walls. It is shown that, due to the intrinsic nonlocal character of the interactions among the particles forming the fluid, the applicability of approximate schemes such as the concept of a surface of tension or morphometric thermodynamics is limited to wall curvatures that are small compared to the range of particle-particle interactions.

Motility of Escherichia coli in a quasi-two-dimensional porous medium

Bacterial migration through confined spaces is critical for several phenomena, such as biofilm formation, bacterial transport in soils, and bacterial therapy against cancer. In the present work, E. coli (strain K12-MG1655 WT) motility was characterized by recording and analyzing individual bacterium trajectories in a simulated quasi-two-dimensional porous medium. The porous medium was simulated by enclosing, between slide and cover slip, a bacterial-culture sample mixed with uniform 2.98-μm-diameter spherical latex particles. The porosity of the medium was controlled by changing the latex particle concentration. By statistically analyzing several trajectory parameters (instantaneous velocity, turn angle, mean squared displacement, etc.), and contrasting with the results of a random-walk model developed ad hoc, we were able to quantify the effects that different obstacle concentrations have upon bacterial motility.

Unbiased diffusion in two-dimensional channels with corrugated walls

This paper deals with diffusion of point particles in linearly corrugated two-dimensional channels. Such geometry allows one to obtain an approximate analytical expression that gives the particle effective diffusivity as a function of the geometric parameters of the channel. To establish its accuracy and the range of applicability, the expression is tested against Brownian dynamics simulation results. The test shows that the expression works very well for long channel periods, but fails when the period is not long enough compared to the minimum width of the channel. To fix this deficiency, we propose a simple empirical correction to the analytical expression. The resulting corrected expression for the effective diffusivity is in excellent agreement with the simulation results for all values of the channel period.

Two-dimensional motion of Brownian swimmers in linear flows

The motion of viruses and bacteria and even synthetic microswimmers can be affected by thermal fluctuations and by external flows. In this work, we study the effect of linear external flows and thermal fluctuations on the diffusion of those swimmers modeled as spherical active (self-propelled) particles moving in two dimensions. General formulae for their mean-square displacement under a general linear flow are presented. We also provide, at short and long times, explicit expressions for the mean-square displacement of a swimmer immersed in three canonical flows, namely, solid-body rotation, shear and extensional flows. These expressions can now be used to estimate the effect of external flows on the displacement of Brownian microswimmers. Finally, our theoretical results are validated by using Brownian dynamics simulations.