## Asymmetric limit cycles within Lorenz chaos induce anomalous mobility for a memory-driven active particle

Rahil N. Valani 1 ∗ and Bruno S. Dandogbessi 2

1 Rudolf Peierls Centre for Theoretical Physics, Parks Road, University of Oxford, OX1 3PU, United Kingdom and 2 Ecole Nationale Sup´rieure de G´nie Math´matique et Mod´lisation (ENSGMM), e e e e Universit´ Nationale des Sciences, Technologies, Ing´nierie et Math´matiques, Abomey, Republique du B´nin e e e e (Dated: October 21, 2024)

On applying a small bias force, non-equilibrium systems may respond in paradoxical ways such as with giant negative mobility (GNM) - a large net drift opposite to the applied bias, or giant positive mobility (GPM) - an anomalously large drift in the same direction as the applied bias. Such behaviors have been extensively studied in idealized models of externally driven passive inertial particles. Here, we consider a minimal model of a memory-driven active particle inspired from experiments with walking and superwalking droplets, whose equation of motion maps to the celebrated Lorenz system. By adding a small bias force to this Lorenz model for the active particle, we uncover a dynamical mechanism for simultaneous emergence of GNM and GPM in the parameter space. Within the chaotic sea of the parameter space, a symmetric pair of coexisting asymmetric limit cycles separate and migrate under applied bias force, resulting in anomalous transport behaviors that are sensitive to the active particle's memory. Our work highlights a general dynamical mechanism for the emergence of anomalous transport behaviors for active particles described by low-dimensional nonlinear models.

Introduction. When a small bias force is applied to a system at equilibrium, it responds with a net drift in the direction of the applied bias. However, non-equilibrium systems can respond in counterintuitive ways to a small applied bias. A peculiar anomalous transport behavior that defies intuition is negative mobility, where a system with no net drift in the absence of a bias force, responds with a net drift in the direction opposite to an applied small bias force. Negative mobility has been observed in experiments of both classical as well as quantum systems, such as electrons in GaAs quantum wells [1], transport through semiconductor superlattices [2, 3] and particles moving in structured microfluidic channels [4]. Theoretically and numerically, negative mobility and related anomalous transport behaviors have been extensively investigated in various minimal models of a passive inertial particle, either deterministic or stochastic, with typical features of dissipation, spatially periodic potential and external time-periodic driving [5-17]. These studies have identified different mechanisms that give rise to negative mobility in these externally-driven non-autonomous systems, and the effect is attributed to at least one of the following ingredients: noise-induced fluctuations, chaotic attractors and periodic limit cycle attractors. Investigation of negative mobility in self-driven active systems has been relatively scarce but it is increasing with a recent surge of interest in active matter [18-23].

In this letter, we consider a simple Lorenz system model of a memory-driven active particle inspired from the experimental system of walking/superwalking droplets [24, 25]. Although the existence of absolute negative mobility in this model has been reported in an

∗ rahil.valani@physics.ox.ac.uk

FIG. 1. Lorenz system model for a memory-driven active particle. (a) Experimental image of a millimeter-sized superwalking droplet [25] and (b) a schematic (not to scale) of the corresponding memory-driven active particle model. A particle located at x d (blue circle) moves in one-dimension with horizontal velocity ˙ x d = X and experiences a memory force Y from its self-generated wave field (blue) built-up by superposition of individual cosine decaying waves generated by the particle at each instant of time (black and gray curves, with the higher intensity of the color indicating the waves created more recently). In addition, the particle experiences an effective drag force -X and an external constant applied bias force F . The wave field height at particle location is Z .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00168v2-figure-8d736434-d206-4312-8722-6a30cda25312.jpg)

earlier work by the first author [23], here we rationalize its emergence by identifying a general dynamical mechanism. We uncover a dynamical mechanism for simultaneous emergence of giant negative mobility (GNM) and giant positive mobility (GPM) rooted in asymmetric limit cycles of the Lorenz system [26].

Active particle model. A vertically vibrating liquid

FIG. 2. Correspondence between phase-space behaviors of the Lorenz system and the trajectories of the memory-driven active particle. Phase-plane projection in the ( X,Z ) plane (left) and space-time trajectory of the particle (right) for a fixed R = 2 5 and . τ = 0 5 (Type I - non-walking . state), τ = 2 (Type II - steady walking state), τ = 5 (Type III -chaotic walk) and τ = 8 9 (Type IV - back-and-forth os-. cillations with a net drift). Time evolution for two different initial conditions ( x d (0) , X (0) , Y (0) , Z (0)) = (0 , 1 0 0) (blue) , , and = (0 , -1 0 0) (red) is shown. , , In left panels, black filled circles represent stable equilibrium points whereas gray filled circles represent unstable equilibrium points of the Lorenz system.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00168v2-figure-0e9dcabb-e8c7-4907-8e39-67e659f1be93.jpg)

bath can support bouncing and walking droplets of the same liquid on its free surface [24, 25, 27]. Each impact in the periodic bouncing of the droplet generates a localized standing wave that decays in time. The droplet interacts with these self-generated waves on subsequent bounces to propel itself horizontally resulting in walking and superwalking droplets (see Fig. 1(a)). Such droplets have been shown to mimic hydrodynamic analogs of several quantum systems [28, 29]. Three key features to note about this system (i) the droplet and its underlying waves co-exist as a wave-particle entity; without the droplet, surface waves decay completely. (ii) There is memory in the system since at large amplitudes of bath vibrations, the droplet-generated waves decay very slowly in time. Hence the droplet motion is guided by the complex dynamical landscape sculpted by the history of waves along its trajectory. (iii) The droplet is self-driven or active [30] since the droplet locally extracts energy from the vibrating bath and converts it into persistent horizontal walking motion. A commonly used theoretical model that captures the key features of a droplet's walking dynamics is the stroboscopic model [31]. In this model, the vertical periodic bouncing motion of the droplet is averaged over and one gets a trajectory equation for the horizontal walking motion describing a particle driven by its wave-memory. As shown schematically in Fig. 1(b), we consider a simple setup where the particle can move only in one horizontal dimension, resulting in the following trajectory equation of motion (in dimensionless units) [31, 32]

$$\ddot { x } _ { d } + \dot { x } _ { d } = F _ { s e l f } + F _ { b i a s },$$

where x d is the position of the particle and an overdot denotes a time derivative. To non-dimensionalize the variables, the spatial length x is scaled by the inverse of the wavenumber k F = 2 π/λ F with λ F as the wavelength of particle-generated waves, and time t is scaled by m/D where m is the particle mass and D is an effective drag coefficient averaged over the periodic bouncing motion (see [31] and [32] for more details). The left-handside of Eq. (1) is composed of an inertial term ¨ x d and a dissipative force term proportional to particle velocity ˙ x d . The first term on the right-hand-side is the wave-memory force on the particle due to its self-generated waves. At each instant of time, the particle generates waves W x ( ) centered at the particle location and decaying exponentially in time. By integrating over the history of the particle's trajectory, one can calculate the overall wave field generated by a linear superposition of the individual waves. The wave-memory force is then proportional to the gradient of this overall wave field at the particle location and can be written as

$$\mathbb { I } _ { \mathbb { I } } \quad F _ { s e l f } = - R \int _ { - \infty } ^ { t } W ^ { \prime } \left ( ( x _ { d } ( t ) - x _ { d } ( s ) ) \right ) \, \mathrm e ^ { - \frac { ( t - s ) } { \tau } } \, \mathrm d s.$$

Here, the parameter R is a dimensionless wave-amplitude of particle-generated waves and the parameter τ is a dimensionless wave decay rate and hence determines the memory of the system (see [31] and [32] for more details on the parameters). The second term on the right hand side is the external bias force applied to the particle and we take it to be a constant force i.e. F bias = F .

Connection with Lorenz system. In experiments, droplet-generated waves have features of both spatial oscillations and spatial decay (see Fig. 1(a)). However, if we only retain spatial oscillations and consider an idealized cosine waveform i.e. W x ( ) = cos( x ), then this still

captures the key qualitative features of walking dynamics [33] and significantly simplifies the system by allowing us to map the integro-differential equation in Eq. (1) into a system of nonlinear ordinary differential equations [3336]. This is done by differentiating with time, the wavememory force

$$Y = F _ { s e l f } = R \int _ { - \infty } ^ { t } \sin ( x _ { d } ( t ) - x _ { d } ( s ) ) \, e ^ { - \frac { ( t - s ) } { \tau } } \, \mathrm d s,$$

and the wave field height at the particle location

$$Z = R \int _ { - \infty } ^ { t } \cos ( x _ { d } ( t ) - x _ { d } ( s ) ) \, e ^ { - \frac { ( t - s ) } { \tau } } \, d s,$$

using the Leibniz integral rule giving the following system of ordinary differential equations (see [32] for a derivation):

$$\dot { X } & = Y - X + F, \\ \dot { Y } & = - \frac { 1 } { \tau } Y + X Z, \quad \quad ( 2 ) \\ \dot { Z } & = R - \frac { 1 } { \tau } Z - X Y, \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

where X = ˙ x d is the particle velocity. With no applied bias force, F = 0, the system in Eq. (2) is equivalent to the celebrated Lorenz system [26] (see Supplemental Material). Hence, the motion of our memory-driven active particle can be directly mapped onto the dynamics of the Lorenz system. We numerically solve the dynamical system of Eq. (2) in MATLAB using the inbuilt ode45 solver.

With the above correspondence between the particle's trajectory equation and the Lorenz system, we can map the attractors and phase-space features of the Lorenz system onto the motion and trajectories of the active particle [32, 33]. We proceed by reviewing different attractors realized in the Lorenz model [32, 37] of Eq. (2) as a function of the memory parameter τ and map them onto the different types of trajectories for the particle. For small values of τ , there is only one equilibrium point of the system at ( X,Y,Z ) = (0 , 0 0) , which is stable. This corresponds to a (Type I) non-walking state for the particle since its velocity is ˙ x d = X = 0 (see Fig. 2(a)). With increasing τ , the Lorenz system undergoes a pitchfork bifurcation (at R = 1 /τ 2 ) where the non-walking state becomes unstable and a symmetric pair of stable equilibrium points emerge, ( X,Y,Z ) = ( ± √ R -1 /τ 2 , ± √ R -1 /τ 2 , 1 /τ ); this corresponds to steady walking motion (Type II) of the particle to right/left (see Fig. 2(b)). Further increase in τ leads to a complex set of homoclinic, heteroclinic and Hopf bifurcations in the Lorenz system [37, 38] and eventually chaos emerges at large τ with a strange attractor in the phase space. Phase-space dynamics on the Lorenz chaotic attractor corresponds to a (Type III) chaotic walk [33, 39] for the active particle where it unpredictably switches between left and right moving states

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00168v2-figure-6acae69b-3a2d-418e-b0c3-a2d6c22c7d2f.jpg)

FIG. 3. Splitting of asymmetric limit cycles in the parameter space induce anomalous transport under applied bias force (see also Supplemental Video). Time-average particle velocity (colorbar) in the ( R,τ ) parameter space of the Lorenz model with (a) no bias force, F = 0, and (b) a small applied bias force, F = 0 05. . The black dashed curve in (b) corresponds to a band of Type IV behavior in (a). These results are from simulations till t = 5000 where only the latter half of the time series was used to calculate the average velocity of the particle. Initial conditions were fixed to ( x d (0) , X (0) , Y (0) , Z (0)) = (0 , 1 0 0). , ,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00168v2-figure-e6924f7f-154b-480d-b654-5885013ee6be.jpg)

as the phase space trajectory jumps between left and right 'wing' of the Lorenz attractor (see Fig. 2(c)). In experiments with a walking droplet moving in free space, the non-walking state, steady walking state and hints of a Hopf bifurcation leading to a walking state with velocity oscillations have been observed, whereas access to the chaotic regime and beyond is not yet achievable [40, 41].

Dynamical mechanism for anomalous mobility. It is common for periodic windows to exist within a chaotic regime of parameter space. Since the Lorenz system obeys a symmetry i.e. the system is symmetric about

FIG. 4. Emergence of giant negative mobility (GNM) and giant positive mobility (GPM). (a) Average velocity (both time and initial-condition averaged) as a function of the bias force F for fixed R = 2 5 and, . τ = 6 9 (red) and . τ = 12 2 . (blue). (b) Average velocity as a function of the memory parameter τ for fixed R = 2 5 and, . F = 0 (gray) and F = 0 05 . (black). Shaded region indicates the spread in the average velocity for different initial conditions. The simulations were done till t = 2000 and only the latter half of the velocity time series was used for time averaging. To average over initial conditions, 50 different linearly-spaced initial conditions were chosen for X (0) ∈ -[ 5 5] with fixed , x d (0) = 0, Y (0) = 0 and Z (0) = 0.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00168v2-figure-1d9d1504-8035-41ab-bc1a-4ca4378f2f1d.jpg)

a rotation of 180 degrees around the Z axis, it restricts the types of periodic orbits that can exist. Two types of limit cycles can be found in the Lorenz system [37]: (i) a symmetric limit cycle about the Z axis which would correspond to a back-and-forth oscillation for the particle with no net drift, and (ii) a symmetric pair of asymmetric limit cycles, as shown in Fig. 2(d), that corresponds to the particle undergoing back-and-forth oscillations along with a net drift either towards the left or right (Type IV). A parameter space plot showing the timeaveraged particle velocity ⟨ ˙ x d ⟩ = ⟨ X ⟩ in the ( R,τ ) space is shown in Fig. 3(a) which helps identify these oscillating states that have a net drift. Region I corresponds to non-walking state ( X = 0), Region II is steady walking ( X = ± √ R -1 /τ 2 ), Region III is mostly filled with chaotic walks with no net drift ( ⟨ X ⟩ = 0). Of particular interest is Region IV bands corresponding to oscillating states with net particle drift. The mixture of positive (blue) and negative (red) drift states along this band indicates the sensitivity of these states to small changes in parameter values.

Now under the application of a small bias force of F = 0 05, . the corresponding parameter space plot in Fig. 3(a) transforms to Fig. 3(b). We clearly see that the band corresponding to Region IV splits into two in the parameter space, one corresponding to the asymmetric limit cycle with a net drift opposite to the applied bias (red) and the other corresponding to the asymmetric limit cycle with a net drift in the direction of the applied bias (blue). Since, on these Region IV bands, the magnitude of net particle drift velocity is around 5 times larger than neighboring parameter values, we identify the former as the giant negative mobility (GNM) band and the latter as the giant positive mobility (GPM) band. Further, on these separated bands, all initial conditions lead to the same limit cycle attractor, unlike the F = 0 band where the basin of attraction of left/right drifting states in the initial-condition space are intricately mixed. After complete separation, the two bands migrate away from each other (see Supplemental Video) with increasing small bias F such that the GNM band moves below the F = 0 band (black dashed curve) while the GPM band migrates above the F = 0 band [42]. We note that this band separation mechanism of coexisting asymmetric limit cycles is closely related to the separation mechanism of phase-locked periodic attractors reported by [14] and [15] for negative mobility of a passive particle in a periodic potential with external time-periodic driving. However, in contrast, the negative and positive mobility bands completely separate and migrate in the parameter space of our Lorenz model, whereas the two bands only partially separate in the system of [14] along with the negative mobility band shrinking and positive mobility band growing in the parameter space.

To demonstrate more clearly how this separation of bands manifests in GNM and GPM for our memorydriven active particle, let us consider the system at two different memory parameter τ values for a fixed R , one value below the region IV band with lower memory τ = 6 9 and the other above the band . τ = 12 2 with . higher memory. If we fixed the system parameter at these values, and investigate how the average velocity (both time and ensemble averaged) varies as a function of the bias force, we get the results shown in Fig. 4(a). For lower memory (red curve), the average velocity shows a strong dip of GNM near F = 0 05. This is because, the . GNM band after separation, migrates through the ( R,τ ) parameter space with increasing F and passes through the ( R,τ ) = (2 5 6 9) . , . point in parameter space near F = 0 05. . Similarly, for the higher memory value (blue curve), the average velocity shows a strong spike near F = 0 05 due to the passing of the GPM band. . Figure 4(b) further shows the sensitivity of these behaviors as a function of the system memory τ . With no bias force (gray curve), there is no net mobility of the particle (when ensemble averaged) at all values of memory τ . However, with a small applied bias, the splitting of the region IV band results in GNM regions for lower mem-

ory and GPM regions for higher memory. This highlights the sensitivity of anomalous transport behaviors to the active particle's memory under small applied bias.

Discussion and conclusions. Although in our simplified Lorenz model we considered an idealized cosine wave, W x ( ) = cos( x ), a wave form more closer to the experimental system of walking droplets that includes spatial decay shows a similar mechanism for GNM and GPM in a smaller region of the parameter space for weak spatial decay, whereas a strong spatial decay appears to suppress this behavior completely (see Supplemental Material). Further, since the present dynamical mechanism lies within the chaotic regime of the parameter space, it would be difficult to observe this mechanism in experiments with walking/superwalking droplets moving in free space as access to the chaotic regime for free walkers/superwalkers is currently not achievable [40]. However, it might be possible to experimentally demonstrate this effect in other active systems driven by lowdimensional chaos, such as autonomous mobile robots whose motion can be easily coupled to the Lorenz system [39, 43-46].

- [1] R. A. H¨ opfel, J. Shah, P. A. Wolff, and A. C. Gossard, Negative absolute mobility of minority electrons in GaAs quantum wells, Phys. Rev. Lett. 56 , 2736 (1986).
- [2] B. J. Keay, S. Zeuner, S. J. Allen, K. D. Maranowski, A. C. Gossard, U. Bhattacharya, and M. J. W. Rodwell, Dynamic localization, absolute negative conductance, and stimulated, multiphoton emission in sequential resonant tunneling semiconductor superlattices, Phys. Rev. Lett. 75 , 4102 (1995).
- [3] E. H. Cannon, F. V. Kusmartsev, K. N. Alekseev, and D. K. Campbell, Absolute negative conductivity and spontaneous current generation in semiconductor superlattices with hot electrons, Phys. Rev. Lett. 85 , 1302 (2000).
- [4] A. Ros, R. Eichhorn, J. Regtmeier, T. T. Duong, P. Reimann, and D. Anselmetti, Absolute negative particle mobility, Nature 436 , 928 (2005).
- [5] R. Eichhorn, P. Reimann, and P. H¨nggi, Brownian moa tion exhibiting absolute negative mobility, Phys. Rev. Lett. 88 , 190601 (2002).
- [6] R. Eichhorn, P. Reimann, and P. H¨ anggi, Paradoxical motion of a single brownian particle: Absolute negative mobility, Phys. Rev. E 66 , 066132 (2002).
- [7] B. S. Dandogbessi and A. Kenfack, Absolute negative mobility induced by potential phase modulation, Phys. Rev. E 92 , 062903 (2015).
- [8] L. Machura, M. Kostur, P. Talkner, J. glyph[suppress] Luczka, and P. H¨ anggi, Absolute negative mobility induced by thermal equilibrium fluctuations, Phys. Rev. Lett. 98 , 040601 (2007).
- [9] J. Spiechowicz and J. glyph[suppress] Luczka, Intriguing properties of transport at the microscales: Langevin equation approach, The European Physical Journal B 97 , 89 (2024).
- [10] M. Wi´ sniewski and J. Spiechowicz, Memory-induced ab-

In summary, we have shown a dynamical mechanism rooted in the celebrated Lorenz system model for the generation of giant negative mobility and giant positive mobility for a memory-driven active particle under a small applied bias. The degenerate pair of asymmetric limit cycles existing for the same parameter values completely separate and migrate in the parameter space upon small applied bias, resulting in the emergence of memory-dependent GNM and GPM for the active particle. We believe this dynamical mechanism to be general to many low-dimensional chaotic systems and it might be at play in generating anomalous transport behaviors for active particles described by low-dimensional nonlinear systems. Further, such mechanisms may be used to tune the internal parameters of the system in order to control a particle's transport behaviors, or in sorting of active particles with different intrinsic properties [47, 48].

Acknowledgments . This research was supported by the Leverhulme Trust [Grant No. LIP-2020-014]. Some of the numerical results were computed using the Hydra cluster in the Department of Physics at the University of Oxford.

solute negative mobility, Chaos 34 , 073101 (2024).

- [11] Y. Fang, Y. Luo, and C. Zeng, Dichotomous noiseinduced negative mass and mobility of inertial brownian particle, Chaos Solit. Fractals 155 , 111775 (2022).
- [12] W. Guo, L.-C. Du, Z.-Z. Liu, H. Yang, and D.-C. Mei, Uphill anomalous transport in a deterministic system with speed-dependent friction coefficient, Chin. Physics B 26 , 010502 (2017).
- [13] M. Wi´ sniewski and J. Spiechowicz, Anomalous transport in driven periodic systems: distribution of the absolute negative mobility effect in the parameter space, NJP 24 , 063028 (2022).
- [14] D. Speer, R. Eichhorn, and P. Reimann, Transient chaos induces anomalous transport properties of an underdamped brownian particle, Phys. Rev. E 76 , 051110 (2007).
- [15] A. Sglyph[suppress]apik, l J. Luczka, and J. Spiechowicz, Negative moglyph[suppress] bility of a brownian particle: Strong damping regime, Commun. Nonlinear Sci. Numer. Simul. 55 , 316 (2018).
- [16] M. Wi´ sniewski and J. Spiechowicz, Paradoxical nature of negative mobility in the weak dissipation regime, Chaos: An Interdisciplinary Journal of Nonlinear Science 33 , 063114 (2023).
- [17] J. Spiechowicz, P. H¨nggi, and J. glyph[suppress] uczka, Coexistence of a L absolute negative mobility and anomalous diffusion, New Journal of Physics 21 , 083029 (2019).
- [18] B.-q. Ai, W.-j. Zhu, Y.-f. He, and W.-r. Zhong, Giant negative mobility of inertial particles caused by the periodic potential in steady laminar flows, J. Chem. Phys. 149 , 164903 (2018).
- [19] F. Cecconi, A. Puglisi, A. Sarracino, and A. Vulpiani, Anomalous mobility of a driven active particle in a steady laminar flow, J. Condens. Matter Phys. 30 , 264002 (2018).

- [20] J.-C. Wu, F.-J. Lin, and B.-Q. Ai, Absolute negative mobility of active polymer chains in steady laminar flows, Soft Matter 18 , 1194 (2022).
- [21] P. K. Ghosh, P. H¨ anggi, F. Marchesoni, and F. Nori, Giant negative mobility of Janus particles in a corrugated channel, Phys. Rev. E 89 , 062115 (2014).
- [22] P. Rizkallah, A. Sarracino, O. B´nichou, and P. Illien, Abe solute negative mobility of an active tracer in a crowded environment, Phys. Rev. Lett. 130 , 218201 (2023).
- [23] R. N. Valani, Anomalous transport of a classical waveparticle entity in a tilted potential, Phys. Rev. E 105 , L012101 (2022).
- [24] Y. Couder, S. Proti`re, E. Fort, and A. Boudaoud, Dye namical phenomena: Walking and orbiting droplets, Nature 437 , 208 (2005).
- [25] R. N. Valani, A. C. Slim, and T. Simula, Superwalking droplets, Phys. Rev. Lett. 123 , 024503 (2019).
- [26] E. N. Lorenz, Deterministic nonperiodic flow, J. Atmos. Sci. 20 , 130 (1963).
- [27] Y. Couder, E. Fort, C.-H. Gautier, and A. Boudaoud, From bouncing to floating: noncoalescence of drops on a fluid bath, Phys. Rev. Lett. 94 , 177801 (2005).
- [28] J. W. M. Bush and A. U. Oza, Hydrodynamic quantum analogs, Rep. Prog. Phys. (2020).
- [29] J. W. M. Bush, V. Frumkin, and P. J. S´enz, Perspeca tives on pilot-wave hydrodynamics, Appl. Phys. Lett. 125 , 030503 (2024).
- [30] S. Ramaswamy, The mechanics and statistics of active matter, Annu. Rev. Condens. Matter Phys. 1 , 323 (2010).
- [31] A. U. Oza, R. R. Rosales, and J. W. M. Bush, A trajectory equation for walking droplets: hydrodynamic pilotwave theory, J. Fluid Mech. 737 , 552 (2013).
- [32] R. N. Valani, Infinite-memory classical wave-particle entities, attractor-driven active particles, and the diffusionless Lorenz equations, Chaos 34 , 013133 (2024).
- [33] R. N. Valani, A. Slim, D. Paganin, T. Simula, and T. Vo, Unsteady dynamics of a classical particle-wave entity, Phys. Rev. E 104 , 015106 (2021).
- [34] J. Mol´ aˇek, c Bouncing and walking droplets : towards a hydrodynamic pilot-wave theory , Ph.D. thesis, Massachusetts Institute of Technology (2013).
- [35] M. Durey, Bifurcations and chaos in a Lorenz-like pilotwave system, Chaos 30 , 103115 (2020).
- [36] R. N. Valani, Lorenz-like systems emerging from an integro-differential trajectory equation of a one-
- dimensional wave-particle entity, Chaos 32 , 023129 (2022).
- [37] C. Sparrow, The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors (Springer-Verlag, New York, 1982).
- [38] E. A. Jackson, Models based on third order differential systems, in Perspectives of Nonlinear Dynamics , Vol. 2 (Cambridge University Press, 1990) p. 125-230.
- [39] R. N. Valani and D. M. Paganin, Attractor-driven matter, Chaos 33 (2023), 023125.
- [40] V. Bacot, S. Perrard, M. Labousse, Y. Couder, and E. Fort, Multistable free states of an active particle from a coherent memory dynamics, Phys. Rev. Lett. 122 , 104303 (2019).
- [41] M. Hubert, S. Perrard, M. Labousse, N. Vandewalle, and Y. Couder, Tunable bimodal explorations of space from memory-driven deterministic dynamics, Phys. Rev. E 100 , 032201 (2019).
- [42] We note that with applied bias force and for small τ , negative steady walking states (Type II) (i.e. in the direction opposite to the applied bias) exist as well. However, we don't classify these as anomalous behaviors since these states exist even for F = 0, and furthermore, with increasing F , the magnitude of the negative steady velocity decreases for these states.
- [43] S. Paramanick, A. Pal, H. Soni, and N. Kumar, Programming tunable active dynamics in a self-propelled robot, The European Physical Journal E 47 , 34 (2024).
- [44] P. Sooraska and K. Klomkarn, 'No-CPU' chaotic robots: From classroom to commerce, IEEE Circuits Syst. Mag. 10 , 46 (2010).
- [45] Y. Nakamura and A. Sekiguchi, The chaotic mobile robot, IEEE trans. robot. autom. 17 , 898 (2001).
- [46] C. Li, Y. Song, F. Wang, Z. Wang, and Y. Li, A bounded strategy of the mobile robot coverage path planning based on Lorenz chaotic system, Int. J. Adv. Robot. Syst. 13 , 107 (2016).
- [47] A. Sglyph[suppress]apik, l J. glyph[suppress] Luczka, P. H¨nggi, and J. Spiechowicz, Tuna able mass separation via negative mobility, Phys. Rev. Lett. 122 , 070602 (2019).
- [48] A. Sglyph[suppress]apik, l J. Luczka, and J. Spiechowicz, Temperatureglyph[suppress] induced tunable particle separation, Phys. Rev. Appl. 12 , 054002 (2019).
Rahil N. Valani Bruno S. Dandogbessi
Asymmetric limit cycles within Lorenz chaos induce anomalous mobility for a memory-driven active particle
On applying a small bias force, non-equilibrium systems may respond in
paradoxical ways such as with giant negative mobility (GNM) -- a large net
drift opposite to the applied bias, or giant positive mobility (GPM) -- an
anomalously large drift in the same direction as the applied bias. Such
behaviors have been extensively studied in idealized models of externally
driven passive inertial particles. Here, we consider a minimal model of a
memory-driven active particle inspired from experiments with walking and
superwalking droplets, whose equation of motion maps to the celebrated Lorenz
system. By adding a small bias force to this Lorenz model for the active
particle, we uncover a dynamical mechanism for simultaneous emergence of GNM
and GPM in the parameter space. Within the chaotic sea of the parameter space,
a symmetric pair of coexisting asymmetric limit cycles separate and migrate
under applied bias force, resulting in anomalous transport behaviors that are
sensitive to the active particle's memory. Our work highlights a general
dynamical mechanism for the emergence of anomalous transport behaviors for
active particles described by low-dimensional nonlinear models.