Abstract
Information storage is a key element of autonomous, out-of-equilibrium dynamics, especially for biological and synthetic active matter. In synthetic active matter however, the implementation of internal memory in self-propelled systems is often absent, limiting our understanding of memory-driven dynamics. Recently, a system comprised of a droplet generating its guiding wavefield appeared as a prime candidate for such investigations. Indeed, the wavefield, propelling the droplet, encodes information about the droplet trajectory and the amount of information can be controlled by a single scalar experimental parameter. In this work, we show numerically and experimentally that the accumulation of information in the wavefield induces the loss of time correlations, where the dynamics can then be described by a memory-less process. We rationalize the resulting statistical behavior by defining an effective temperature for the particle dynamics where the wavefield acts as a thermostat of large dimensions, and by evidencing a minimization principle of the generated wavefield.
Similar content being viewed by others
Introduction
Many simple active biological systems possess memory mechanisms which are commonly thought to play a key role in their statistical dynamical behaviors. However, assessing the influence of memory in such active systems is an ill-defined task1. Indeed, what is commonly called biological memory involves mechanisms acting at different time scales from allosteric switching (~10−5−10−3 s) to biochemical circuits (~10−2−1 s). In contrast, in synthetic active matter, already crafting a system able to self-propel, such as light-activated colloids2, colloidal rollers3 or self-propelled disks4, is by itself an experimental tour de force. As a consequence, these systems are intentionally designed to be minimalistic so there is a chance to rationalize them. Implementing a memory repository in these physical systems would raise complex experimental issues, althought, some recent robotic5 or electronically back-controlled6 strategies could solve this issue. Correlation times may potentially play the role for a memory of a physical system7, but in practice, their tunabilty and controllabilty upon a variation of experimentally controllable parameters are often limited. For all these theoretical and experimental reasons, physicists, including the authors themselves, are often ill-at-ease to rationalize the influence of memory in physical or biological active systems. As a consequence, this question has often been evaded in the field of active matter8,9,10, even if this multiscale and ill-posed concept has fascinated and puzzled for long11.
In the last two decades, a candidate for such an investigation appeared with walking droplets, or walkers for short12,13,14. Experimentally, a droplet bounces periodically on a vertically and sinusoidally oscillating oil surface. The result is the emergence of a complex standing wavefield15,16 which propels the droplet and stores information about its past positions17,18,19,20,21,22,23,24,25,26. Indeed, each bounce of the droplet imprints a standing long-lasting axisymmetric wave centered at the point of impact, which lifetime is controlled by the vertical acceleration magnitude of the bath. In this system, because the droplet slides down the gradient of the local liquid surface, the wavefield acts as a memory that drives the droplet motion. Crucially, the amount of information encoded in the wavefield is completely controllable in a continuous fashion through the magnitude of the acceleration of the interface.
This unique feedback between the droplet and the wavefield dynamics is at the core of a stream of research mainly motivated by analogies with quantum systems27,28,29,30,31,32,33,34,35,36,37,38,39. In addition, recent numerical, theoretical and experimental studies40,41,42,43,44,45,46,47 have shown that the memory of the walker leads to run-and-tumble-like chaotic dynamics48,49,50, similar to Marangoni-driven drops51, or particles in in-silico superfluids52. All these studies40,41,42,43,44,45,46,47 identify the key role of the wavefield memory in the emergent statistical behavior. Nevertheless, the dynamics and statistical properties of the wavefield are poorly understood. Indeed, most of the dynamical descriptions of the walker wavefield focus on the regime of so-called “low memory” where memory time is small33,53, or in the infinite memory limit where correlations become irrelevant42. In between, in the high memory regime, the complex and highly correlated interaction between the droplet and the wavefield makes theoretical analysis a daunting task, while numerical investigations are requiring large resources to obtain statistically relevant measures.
In this article, using simulations and experiments, we show that a walker trapped in a weak harmonic potential with large memory reaches an active statistical limit where the wavefield becomes a memory-controlled thermostat. We show that an excess of memory leads to an effective memory-less particle dynamics, paving the way for further understanding of highly correlated memory-driven dynamics.
Results
Walker as a memory-driven agent
A walker is the association of a sub-millimetric oil droplet, periodically bouncing on a vertically-vibrated oil surface, and the guiding standing waves generating by the drop bounces (see Fig. 1a)17,18,19,20,21,22,23,24,25,26 (see Methods). In this article, for practical reasons, the space available to the walker is bounded by confining the particle with an applied external potential. Experimentally, the external potential is applied to the ferrofluid core of the droplet through a combination of magnetic fields originating from electric coils and a permanent magnet (see Fig. 1a and Methods). The harmonic potential per unit mass is noted U = ω2∣r∣2/2, with ω/2π ~ 0.01−0.4 Hz its frequency. It is worth noticing that, since the liquid surface has no ferromagnetic susceptibility, it is not altered by the presence of the magnetic field. In addition, previous studies have focused on the different ways to confine a walker, such as circular or elliptic confinements31,37, or quasi one-dimensional channels34,35. Nevertheless, these confinement strongly alter the shape of waves and add a new layer of complexity to the walker dynamics. For these reasons, in this article, the confinement is only originating from a weak external force to be as close as possible to a free walker dynamics.
The dynamics is simulated by solving the recurrent coupled equations relating the speed v of the walker and wavefield ζ through [Eqs. 1, 2]29
and
β is the effective damping coefficient applied to the walker, τF is the bouncing period, U is the external confining potential per unit mass, c is the wave coupling coefficient, J0 is the bessel function of first kind and order zero, λF is the wavelength and δ is a spatial damping factor. Me is the memory parameter which controls the amount of information in the wavefield. Experimentally, Me is controlled by the temporal dampening of the wave τ, controlled by the proximity of the applied vertical acceleration magnitude γm of the surface to the Faraday threshold γF, i.e., \({{{{{{{\rm{Me}}}}}}}}=\tau /{\tau }_{F}={(1-{\gamma }_{m}/{\gamma }_{F})}^{-1}\). One should note that several models24,25,29,54, varying in their strategies to solve the walker’s equation, have been proposed, and share the same core ingredients. They are in good qualitative agreement with the model we use. In all models, the control parameters is the amount of information, and not the duration of a single delay, in contrast with recent electronic systems with temporal feedbacks6. Finally, to be precise, the wave coupling slightly depends on the speed at each impact which is denoted by \({{{{{{{\mathcal{O}}}}}}}}\left({{{{{{{\bf{v}}}}}}}}\left({t}_{N}\right){{{{{{{\boldsymbol{\nabla }}}}}}}}\zeta {({{{{{{{{\bf{r}}}}}}}}}_{N},{t}_{N})}^{2}\right)\) (see SI and Methods).
Figure 1b and c illustrate the dynamics from the waves and particle point of view. The wavefield stores information from a chain of standing-waves sources following the walker trajectory. This chain characteristic length scales as ~MeVτF, with V being the mean particle speed. We highlight that it is not possible to reduce this physical system (droplet and waves) to a point-like particle whose motion depends only on its current position and velocity. The droplet is able to read its memory (the two last terms in Eq. (1)) and to edit its memory (through the Bessel function in Eq. (2)), similarly to a Turing machine55. Erasing the memory is also possible with a specific protocol55. The concept of memory is therefore justified since the walker can therefore write, read and also edit the information it inscribes onto the liquid interface.
Dynamics of the particle: overview of the regimes with the memory
In Eq. (2), the memory parameter Me acts as a control parameter which allows to numerically span several and different dynamical regimes (Fig. 1d, e and f, see also Fig. SI 2 for complementary numerical results at large Me). At low memory (typically Me < 20), only circular trajectories are observed56 and the walker speed increases monotonously with Me. At intermediate values of the memory parameter (typically 20 < Me < 135, Fig 1d), a quantized set of close-looped trajectories are observed and result from a wave-energy minimization33,53. As the memory is further increased in this intermediate regime (Fig 1e), an intermittent and chaotic dynamics is triggered where the trajectory navigates between the many possible eigenmodes of the dynamics, as investigated experimentally in32 and theoretically in26,42. While a chaotic behavior has been measured and characterized in this regime, the dynamics still shows strong auto-correlation, as proved by the appearance of circle, lemniscates or loops in the trajectory. This coherence fades away as the memory is increased even further (Fig. 1f). In the high memory regime (Me > 135), the corresponding dynamics does not present any signature of patterns reported in the literature33 and instead the dynamics shows a fully developed chaos without apparent underlying structure which results from a Shil’nikov bifurcation40. In the limit Me → ∞, the temporal decay of the waves is removed and each past impact contributes with the same intensity to the wavefield, which is only investigated neatly in simulations and theory57,58. In the intermediate, high and infinite memory regimes, the mean walker speed is relatively constant while the persistence length of the dynamics decreases with the memory parameter40. In what follows, we focus only on the mainly-unexplored high memory regime.
We report in Fig. 1g a typical trajectory obtained experimentally in the high memory regime. Both experimentally and numerically, the trajectories are disordered with the presence of many loopy paths without apparent underlying structure. Over a long period of time, the overall trajectory shares the same symmetry of the confining potential (Fig. 1h) and the radial probability density function is Gaussian (Fig. 1i). Interestingly, the trajectories for lower memory parameters do not lead to Gaussian probability density functions (PDF) (Fig. SI 1). In this memory regime, the self-organization of the overall trajectory (Fig. 1d) and the intermittent dynamics (Fig. 1e) leads to structured statistics, which are not Gaussian.
High-memory dynamics of the wavefield
The dynamical rules of the walker evolution are mediated by the information stored in a wavefield, so that the statistical properties of the trajectory are directly related to those of the wavefield through the force Fw ∝ − ∇ζ(rN, tN) (see Eq. (1)). We expand the wavefield onto a cylindrical Bessel frame of reference by using Graf’s addition theorem59 (see SI). The complex quantity
is the weight of the eigenmode invariant by a rotation of 2π/n, with \(n\in {\mathbb{Z}}\). While the simulations assume δ = 2.5λF (see Methods), the limit δ → ∞ has been used to compute the modes an (Eq. (3)). This choice does not remove the fundamental features of the wavefield and walkers dynamics. Even though there is an infinity of eigenmodes an, only a handful number of modes are necessary, which are determined by the confinement applied to the walker. Indeed, for n ≥ 1, Jn(2πrp/λF) ≈ 0 around rp = 0, as well as all their (n − 1)-th derivatives. Consequently the interval of rp such that Jn(2πrp/λF) ≪ 1 increases as the index n increases too. As a results only a limited number of modes contributes effectively to the wave field.
Typical time series for the amplitudes ∣an∣ are erratic as a consequence of the chaotic dynamics of the droplet40 (Fig. 2a and Fig. SI 3a, real part of the associated eigenmodes are shown in Fig. 2b). Given the apparent lack of simple structure in their time series, we investigate their associated statistical PDF. The PDF P(∣an∣2) shows an exponential decrease (Fig. 2c), indicating that P(∣an∣) follows a Gaussian distribution as would be a sum of uncorrelated memory-less events. This evolution is especially valid for large ∣an∣. In addition, the distribution for ∣a0∣ differ from the others, presumably because this mode plays a special role since it shares the same spatial symmetry as the harmonic potential containing the walker. Note that even though the modes n > 0 do not share the potential symmetry, they do not vanish. This is due to the temporal decay of the waves, which lead to an axial symmetry breaking of the source positions. Furthermore, the different wave modes of amplitude an are weakly correlated to each other (Fig. 2d, Fig. SI 3d). This correlation decreases with the memory parameter Me.
The standard deviation of P(∣an∣2) is found to be identical for 0 < n < 7, which hints toward an equipartition of energy within the eigenmodes. Measuring the root mean squared (rms) width \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}=\sqrt{\langle | {a}_{n}{| }^{2}\rangle }\) of the PDF (Fig 2e, Fig. SI 3e), we observe that \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) decreases slowly with n (barely a 10% decrease for the first 25 modes) indicating that many modes fluctuate strongly. Furthermore, the modes PDF does not change significantly with the memory parameter Me. Indeed, as Me changes from 500 to 10,000, the value of \({a}_{0}^{{{{{{{{\rm{rms}}}}}}}}}\) changes approximately from 4ζ0 to 6ζ0.
We further analyze the wavefield dynamics by computing the field intensity E which we defined as (see SI) [Eq. 4]
The PDF P(E) is well fitted by a Gamma distribution, especially for large Me (Fig. 2f, Fig. SI 3f, see SI for the fitting parameters). Such a choice is guided by the fact that a sum of independent variables whose distribution follows the same exponential distribution follows a Gamma distribution. Given the moderate correlations between the effectively-contributing modes an, the exponential distribution of the \(| {a}_{n}^{2}|\) and the weak variation of \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) with n, this choice is expected to approximate correctly the wave-intensity distribution. Also we observe that larger Me lead to an increase and widening of P(E), which is a consequence of the increasing \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) (Fig. 2e). Finally, we note that the numerical PDFs are in quantitative agreement with the experimentally-highest-reachable Me (Fig. SI 4). Three independent experiments at similar Me and different frequencies also show Gaussian P(∣an∣), a slow decrease of \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) with n, and a Gamma-like PDF of the wavefield intensity.
To isolate the influence of the correlations between the successive walker positions, we compare the wave-driven dynamics of the walker with a random wavefield generated by the superposition of 5 × Me random sources. The radial PDF of the random sources is chosen to be equal to the radial PDF of the walker position for Me = 2500. This random field is then statistically equivalent to the wavefield generated from randomized positions of the particle. Similarly to the case of walkers, the PDFs P(∣an∣) are Gaussian (Fig. 2c, Fig. SI 3c), \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) decreases over one hundred units of n (Fig. 2e, Fig. SI 3e), and P(E) can be fitted by a Gamma function (Fig. 2f, Fig. SI 3f).
Yet, important differences exist. Contrarily to the walkers case, the correlation matrix shows no correlation between modes (Fig. 2d, Fig. SI 3d), a feature which results from our randomly-generated field. As Me increases, the different \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) for the random field increase by a common factor. Indeed, for a random distribution of source positions P(r) with standard deviation σ, \({a}_{n}^{{{{{{{{\rm{rms}}}}}}}}}\) reads [Eq. 5]
where In is the modified Bessel function of first kind of order n (see SI). This result plotted in Fig. 2e implies that the number of wave modes in the case of a random distribution is only determined by the standard deviation of the distribution of sources. On the contrary, the walker case shows an increase of the number of modes storing energy while keeping a roughly identical distribution of wave sources P(∣r∣) as discussed in further details in the following sections. The increase of modes effectively contributing to the wavefield can be understood as effective degrees of freedom (DOF) which we define by [Eq. 6]
The DOF does not depend on Me for a random wavefield, while the DOF increases logarithmically with the memory in the case of a walker, as if the wavefield were acting as a reservoir of increasing dimension (Fig. 2g, Fig. SI 3g). As a conclusion, the memory parameter gives a control on the properties of the reservoir surrounding the walker.
A last difference between the randomized and walker wavefield lies in the evolution of the mean wave intensity E with the memory Me. The average wave intensity increases differently in the case of the walker and a random wavefield (Fig. 2h, Fig. SI 3h). The mean intensity of the random field is proportional to Me (R2 = 0.999) as expected from our theoretical prediction Eq. (5), (see SI) while the correlated chain-like distribution created by walkers presents a sublinear scaling Mep with p = 0.382 ± 0.024 (R2 = 0.999). For a random field, the intensity per unit memory is E/Me ~ 1 which means that each source contributes in average equally. In contrast, the intensity per unit memory for the walker wavefield yields E/Me ~ Me−0.62. This decaying evolution indicates that the system tends to decrease significantly its wave intensity by selecting trajectories which promote destructive interference. A significant decrease of the wavefield intensity has already been observed for lower values of the Me with quantified trajectories33,53. Here we show that via destructive interference, this wave minimization mechanism extends to more complex and chaotic trajectories at high Me. Our analysis proves that the averaged modes depend not only on the probability density function but also of higher-order correlation functions, a difference with the result of Durey et al.42 which we expect to hold when the memory time exceeds the averaging period of time required to insure ergodicity.
Markovian walker dynamics from an overload of memory
We now consider the influence of the dimensional increase of the dynamics resulting from the additional wave DOF triggered by the walker, with increasing Me. The time series of the force Fw,x(t) experienced by the walker solely originating from the waves along the x direction is erratic, as expected from previous studies40 (Fig. 3a, Fig. SI 5a). The y direction is statistically identical given the axisymmetry of the confining potential. For the smallest memory parameter illustrated (Me = 500), we observe a correlated signal (Fig. 3b, Fig. SI 5b), especially for low frequency of the harmonic potential. This correlation exists because of the intermittent dynamics of the walker55.
As the memory parameter Me is increased, the correlation is lost and the auto-correlation \({C}_{x}(t)={\langle {F}_{w,x}({t}_{0}){F}_{w,x}(t+{t}_{0})\rangle }_{{t}_{0}}\) converges toward a sharp peak at t = 0, which can be approximated by a Dirac function, namely Cx(t) ≃ 2Dδ(t) where D defines an effective diffusion coefficient. This feature is even more pronounced when renormalising the correlation time by the memory time (Fig 3c, Fig. SI 5c). Indeed, the time-correlation within Fw,x becomes practically non-existent when compared to the amount of information stored into the wavefield for increasing Me. At memory parameter Me = 500 the correlation time is of the order of 10−2τ, while for memory parameter Me = 25,000 the correlation time is of the order of 10−4τ.
Along the tangential direction however, the auto-correlation for the force Fw,T does not converge to zero. Indeed, \({C}_{T}(t)={\langle {F}_{w,T}({t}_{0}){F}_{w,T}(t+{t}_{0})\rangle }_{{t}_{0}}\) presents a sharp peak at the origin and a plateau at longer time which ensures the average propulsion of the particle at mean speed v0. Small oscillations at high frequency are observed at short time scales, in particular at small Me. They correspond to the oscillations of the walker speed, which can trigger chaotic behaviors as described in40. As a consequence, the force exerted by the wavefield on the particle can be divided into two contributions, a constant tangential component ensuring self-propulsion, and a random component.
The PDF of the force Fw,x, P(Fw,x), is Gaussian for all values of Me > 500 (Fig. 3d, Fig. SI 5d). It is worth noticing that the standard deviation of this probability density function shows a very weak increase with Me. The associated power spectral density S(Fw,x) (Fig. 3e, Fig. SI 5e) shows that for the highest memory investigated, Me = 25,000, S(Fw,x) is flat over three orders of magnitude in frequency, equivalent to a white noise. Observed deviations from the flat power spectrum have two origins. First, for smaller values of Me and low frequency (Fig. SI 5e), a small bump is observed around \({t}^{-1} \sim 4. 1{0}^{-3}{\tau }_{F}^{-1}\), which corresponds to the characteristic orbital period of a walker at low memory parameter33,53,56. This deviation vanishes for large memory parameters. The second source of deviation is observed at high frequency for all Me. It can be attributed to the non-vanishing correlations over a few (≈10) bounces as discussed in40.
As a consequence of those observations, from the particle point-of-view, the wave reservoir eventually preserves the self-propulsion, and the fluctuating component acts as a white noise force. The wave force can be described as a combination of a simple deterministic propelling force Fp(v) aligned with the velocity which contains the correlations at short time scales, and a memory-less white noise η(t). It is surprising to lose all correlations in a memory-driven dynamics such as it becomes approximated by a Markovian process.
Effective temperature induced by a memory
As a result of the previous observations, the probability P(r, v) to find the particle at a position r and a velocity v can be approximated with good accuracy by a Fokker–Planck equation60. Neither in the experimental analysis, nor in the numerical simulations, were correlations between position and velocity identified. Hence we infer the ansatz P(r, v) = P(r)P(v). For practical purposes, we define the time-average kinetic energy per unit mass \(K=\frac{1}{2}\langle {v}^{2}\rangle\) which is a quantity independent of ω. K is also found to vary less than 0.6% with Me in the range [200:25,000]. The position PDF for the position P(r) is well described by a Boltzmann-Gibbs probability density function \({{{{{{{\mathcal{P}}}}}}}}({{{{{{{\bf{r}}}}}}}})=\alpha {\omega }^{2}| {{{{{{{\bf{r}}}}}}}}| \exp \left(-\beta {\omega }^{2}| {{{{{{{\bf{r}}}}}}}}{| }^{2}/2\right)\) with α a normalization factor and β−1 the equivalent temperature of our system (per unit mass) (Fig. 4a). This observation holds for both numerical simulations and experiments, and for all frequencies investigated (ω/2π ∈ [0.01;0.25]) and Me > 200.
As suggested in61, we expect the velocity PDF P(v) to be given by \(\alpha ^{\prime} \exp \left(-{{\Phi }}({{{{{{{\bf{v}}}}}}}})/D\right)\), with Fp = −∇vΦ, D an effective diffusion coefficient and Φ(v) a velocity potential to be determined. The velocity potential \({{\Phi }}({{{{{{{\bf{v}}}}}}}})={\phi }_{0}{\left(| {{{{{{{\bf{v}}}}}}}}| -{v}_{0}\right)}^{2}/2\) shows a good agreement between the theoretical prediction and the numerical speed PDF (Fig. 4b). The experimental data also collapses adequately onto the master curve. It is interesting to compare the velocity potential used here, namely \({{\Phi }}({{{{{{{\bf{v}}}}}}}})={\phi }_{0}{\left(| {{{{{{{\bf{v}}}}}}}}| -{v}_{0}\right)}^{2}/2\), with the one used in previous investigations performed at lower memories36,56 which is stiffer with the presence of v4 terms. This suggest that the constrain on the self-propulsion speed v0 is softer at high memory than at short memory.
However, the properties of P(r) and P(v) related to the walker dynamics do not change strongly with memory. This can be measured by computing the standard deviation of the PDFs shown in Fig. 4a, which also corresponds to the effective system temperature of the walker. We measure this quantity over two decades of values of the memory parameter by fitting the PDF P(r) with a Gaussian (Fig. 4c). We obtain β−1 ≃ 3.1 ± 0.5K−1. A very weak evolution with the memory Me could be argued, in which case a power law fitting β−1 ∝ Meν gives a very small exponent ν = 0.096 ± 0.021 (95% confidence interval) when the fit is applied to all frequencies (0.010–0.25 Hz) at once to increase the precision. Mostly all the numerical and experimental data fall onto the same master curve, at the exception of the numerical points at ω/2π = 0.25 Hz and experimental data at ω/2π = 0.2 Hz where an inflection is observed, leading to a deviation of β with respect to the other frequencies. This behavior might be related to the external potential which is strong and for which the chaotic dynamics are at the edge of a strongly-developed chaos. Finally, while β−1 increases very weakly with Me, β−1 strongly correlates with the wave intensity E linking back the features of the wavefield to the dynamics of the walker (Fig. SI 6).
Discussion
The analysis conducted in this article reveals a complex interplay between the agent and its wavefield. If Eqs. (1) and (2) together describe a deeply non-Markovian dynamics, the build-up of the memory wavefield with increasing Me breaks the correlations in the dynamics of the walker and leads to a stochastic dynamics described in average by a white-noise-driven Markovian dynamics. In addition, the observations related to the wavefield statistical behavior allow us to conclude that the erratic nature of the droplet dynamics does not arise from a chaotic mode mixing. Indeed, as evidenced in Fig. 2d and also by Fig. SI 3d, the eigenmodes of the wavefield are still moderately correlated. Therefore, the main statistical properties cannot result from a chaotic mixing of the wave modes but rather from an increase, up to a factor seven, of wave DOF (see Fig SI 3g) which stems from the increasing memory parameters and amount of information.
We finish by discussing our results in the context of cortical waves following a thought-provoking and inspirational review article by Muller et al.62. Although there are differences between the two systems, and our investigation is not motivated by neuroscience, the parallel between the two systems is very intriguing. In the visual cortex, it has been shown that stimulus-evoked responses can be described by a stationary bump and a propagative wave over a domain of several millimeters. These two responses encode separately the position and the initial time of the stimulus. As a consequence, two or more stimuli generate two or more responses centered at different loci which superpose onto each other. Although, the superposition mechanism is more complex than the simple linear superposition of density waves, it has been proposed that the wave state resulting for the superposition of several cortical waves may serve as computational principles. The case of many spatio-temporally-separated stimuli, which, in our analogy, would correspond to many secondary sources, would be of particular importance. The correlation between these stimuli may be of crucial importance in information processing, as the statistical features of a field resulting either from spatio-temporally-correlated or from uncorrelated sources strongly differ.
In this article, we numerically implemented a deterministic dynamical system which stores information and showed an experimental proof of principle. In this system, the past trajectory of the particle is encoded into a standing wavefield built by the drop previous bounces, which in return propels the particle. The unique property of this system is the control of the amount of information stored. Here, this information storage shapes the properties of the wavefield which acts as a indirect controllable thermal reservoir for the particle. In the long memory regime, we have shown that the wave reservoir conserves a short-term correlation, which is sufficient to maintain a propulsion. In the long-term limit, the wavefield possesses the properties of a white noise. The intensity stored in the waves does not diverge with memory and is self-regulated by means of destructive interference. It is very striking to observe that a system with multiple readable memories can be used to embed the properties of a thermal reservoir.
Methods
Experiments
In a 14 cm-circular tank filled with 5 mm of silicon oil with viscosity η = 20 cp and surface tension σ = 20.9 × 10−2 N/m, a droplet of diameter D = 700 ± 50 μm of the same silicon oil and with a core of ferrofluid is deposited. The coalescence of the droplet in the tank is prevented by vertically and sinusoidally oscillating the tank at 80 Hz and with an acceleration magnitude around γm = 3.8g, g being the gravitational acceleration (see Fig. 1a). This leads to a droplet bouncing period which is twice the tank oscillation period, i.e., the Faraday period τF = 2/f = 2.5 × 10−2 s and a wavelength which is the Faraday wavelength λF = 4.75 mm. A small amount (≈5% in volume) of ferrofluid (iron-cobalt nanoparticles in glycerol, magnetic susceptibility χ = 2.6) is encapsulated inside the drop, so that it becomes paramagnetic. Using a combination of a homogeneous magnetic field generated by two coils in the Helmholtz configuration, and a radial gradient parallel to the bath surface generated by a permanent magnet, we confine the drop in a harmonic well, whose minimum is located at the center of the tank (see Fig. 1a). Details of the calculation for the magnetic confinement can be found in ref. 33. The frequency of the potential well can be tuned by the distance between the permanent magnet and the oil surface. Note that the magnetic fields applied do not interfere with wave generation as the silicon oil has no magnetic properties. The magnetic force acting on the drop was calibrated by recording the motion of the drop on circular trajectories for various values of the potential well frequency. The procedure is described in details in33. The particle trajectory is tracked using image processing, and each drop is recorded for 1 h, corresponding to 1.5 × 105τF and a travelled distance of 7.2 × 103λF. We analyze the trajectory statistics for a memory parameter Me = 250 ± 50, which is directly related to the magnitude of the oil surface acceleration by the formula \({{{{{{{\rm{Me}}}}}}}}={(1-{\gamma }_{m}/{\gamma }_{F})}^{-1}=\tau /{\tau }_{F}\), γF being the threshold of the Faraday instability. The value Me = 250 ± 50 corresponds to the largest memory value that we were able to reach with our experimental setup. We also focus on magnetic potential frequency ω/2π < 0.30 Hz, allowing for large trajectories. Special care has been given to the homogeneity of the vertical vibration, by looking at the regularity and the homogeneity of the Faraday waves above the instability threshold γF.
Numerical simulations
The results presented in this article have been obtained via a discrete step algorithm29,63. In essence, the algorithm modelling the walking droplet dynamics consists of two phases, alternating periodically. The first phase is the “bouncing phase” where the droplet is considered as a perfectly inelastic ball, bouncing on a vertically oscillating rigid surface. If the surface oscillates with a dimensionless acceleration Γ = γm/g = A(2πf)2/g, where A is the amplitude and f the frequency of oscillation, the dynamics of the ball is uniquely determined64. The dimensionless acceleration Γ = 4.12 is chosen so that the bath oscillates twice as fast as the drop. As a consequence, the relative speed at impact and the duration of contact with the interface can be computed. The former information is related to the wave intensity and therefore the amplitude of each wave, and the latter to the duration of interaction with the fluid interface. The second phase is the period during which the ball sits on the interface. During that period, the droplet get a kick of momentum from the standing wavefield in the direction normal to the liquid interface, then the wavefield is updated by adding a new wave source at the drop impact position and finally, the drop loses kinetic energy via friction with the interface. For the kick of momentum, an increase of horizontal speed is applied and reads [Eq. 7]
where V = (v, vz) is the 3D velocity of the droplet, and N = (n, nz) is the 3D vector normal to the surface. Lowercase symbols corresponds to the horizontal component of the 3D vectors. The normal vector N reads [Eq. 8]
Following the kick of momentum, a new standing wave is created on the interface at the impact position while the droplet gets a kick of momentum in the direction normal to the wavefield. The shape of the standing wave and its time evolution obeys Eq. (2), with τF = 2.5 × 102 s, λF = 4.75 mm and δ = 2.5λF. Finally, the droplet speed decreases exponentially with a characteristic time scale τv = 4.5 × 10−2 s during the contact time. The stability and reproducibility of the numerics have been tested and validated in29,40. The initial conditions of the walking dynamics are identical for each simulation. No waves exist on the interface prior to the particle motion. The particle starts its motion with (x0, y0) = (0, λ/2) and (vx,0, vy,0) = (6.66,3.33) mm/s. These values were chosen to break the symmetry of the potential, while being close to the equilibrium speed value. For statistical investigations in Figs. 3 and 4, the algorithm was integrated over 5 × 106 bounces, i.e., period of the waves τF. In order to remove transient dynamics, the first 10% of the dynamics were not considered. In Fig. 2, the algorithm was used over 2.5 × 106 bounces and 125 waves modes were considered. As previously, the first 10% of the dynamics were not considered. PDFs presented in this article are normalized such as the sum of all obtained probability equals to one.
Data availability
All the important data sets that support the findings of our study are available at the link: https://mycore.core-cloud.net/index.php/s/bZsJUfvld9MxYbT.
Code availability
The codes are available at the following link: https://mycore.core-cloud.net/index.php/s/bZsJUfvld9MxYbT.
References
Libchaber, A. & Tlusty, T. Walking droplets, swimming microbes: on memory in physics and life. Comptes Rendus. Mécanique 348, 545–554 (2020).
Palacci, J., Sacanna, S., Steinberg, A. P., Pine, D. J. & Chaikin, P. M. Living crystals of light-activated colloidal surfers. Science 339, 936–940 (2013).
Bricard, A., Caussin, J. B., Desreumaux, N., Dauchot, O. & Bartolo, D. Emergence of macroscopic directed motion in populations of motile colloids. Nature 503, 95–98 (2013).
Deseigne, J., Dauchot, O. & Chaté, H. Collective motion of vibrated polar disks. Phys. Rev. Lett. 105, 098001 (2010).
Rubenstein, M., Cornejo, A. & Nagpal, R. Programmable self-assembly in a thousand-robot swarm. Science 345, 795–799 (2014).
Khadka, U., Holubec, V., Yang, H. & Cichos, F. Active particles bound by information flows. Nat. Commun. 9, 1–9 (2018).
Keim, N. C., Paulsen, J. D., Zeravcic, Z., Sastry, S. & Nagel, S. R. Memory formation in matter. Rev. Modern Phys. 91, 035002 (2019).
Marchetti, M. C. et al. Hydrodynamics of soft active matter. Rev. Modern Phys. 85, 1143–1189 (2013).
Bechinger, C. et al. Active Particles in Complex and Crowded Environments. Rev. Modern Phys. 88, 045006 (2016).
Gompper, G. et al. The 2020 motile active matter roadmap. J. Phys.: Condensed Matter 32, 193001 (2020).
Schrödinger, E.What is Life? The Physical Aspect of the Living Cell (Cambridge University Press, 1944).
Couder, Y., Protière, S., Fort, E. & Boudaoud, A. Walking and orbiting droplets. Nature 437, 208–208 (2005).
Bush, J. W. M. Pilot-Wave Hydrodynamics. Ann. Rev. Fluid Mech. 47, 269–292 (2015).
Bush, J. W. M. & Oza, A. U. Hydrodynamic quantum analogs. Report Prog. Phys. 84, 017001 (2021).
Faraday, M. On a Peculiar Class of Acoustical Figures; and on Certain Forms Assumed by Groups of Particles upon Vibrating Elastic Surfaces. Philosoph. Transact. Royal Soc. London 121, 299–340 (1831).
Miles, J. & Henderson, D. Parametrically forced surface waves. Annu. Rev. Fluid. Mech. 22, 143–165 (1990).
Walker, J. Drops of liquid can be made to float on the liquid. What enables them to do so? Sci. Am. 238, 151–158 (1978).
Protiere, S., Boudaoud, A. & Couder, Y. Particle wave association on a fluid interface. J. Fluid Mech. 554, 85–108 (2006).
Couder, Y., Fort, E., Gautier, C. & Boudaoud, A. From bouncing to floating: non-coalescence of drops on a fluid bath. Phys. Rev. Lett. 94, 177801 (2005).
Vandewalle, N., Terwagne, D., Mulleners, K., Gilet, T. & Dorbolo, S. Dynamics of a bouncing droplet onto a vertically vibrated surface. Phys. Rev. Lett. 100, 167802 (2008).
Eddi, A. et al. Information stored in Faraday waves: the origin of a path memory. J. Fluid Mech. 674, 433–463 (2011).
Moláček, J. & Bush, J. W. M. Drops walking on a vibrating bath: towards a hydrodynamic pilot-wave theory. J. Fluid Mech. 727, 612–647 (2013).
Moláček, J. & Bush, J. W. M. Drops bouncing on a vibrating bath. J. Fluid Mech. 727, 582–611 (2013).
Oza, A. U., Rosales, R. R. & Bush, J. W. M. A trajectory equation for walking droplets: hydrodynamic pilot-wave theory. J. Fluid Mech.737, 552–570 (2013).
Milewski, P. A., Galeano-Rios, C. A., Nachbin, A. & Bush, J. W. M. Faraday pilot-wave dynamics: Modelling and computation. J. Fluid Mech. 778, 361–388 (2015).
Durey, M. & Milewski, P. A. Faraday wave-droplet dynamics: discrete-time analysis. J. Fluid Mech. 821, 296–329 (2017).
Couder, Y. & Fort, E. Single-particle diffraction and interference at macroscopic scale. Phys. Rev. Lett. 97, 1–4 (2006).
Eddi, A., Fort, E., Moisy, F. & Couder, Y. Unpredictable tunneling of a classical wave-particle association. Phys. Rev. Lett. 102, 240401 (2009).
Fort, E., Eddi, A., Boudaoud, A., Moukhtar, J. & Couder, Y. Path-memory induced quantization of classical orbits. Proc. Natl Acad. Sci. 107, 17515–17520 (2010).
Oza, A. U., Harris, D. M., Rosales, R. R. & Bush, J. W. M. Pilot-wave dynamics in a rotating frame: on the emergence of orbital quantization. J. Fluid Mech. 744, 404–429 (2014).
Harris, D. M., Moukhtar, J., Fort, E., Couder, Y. & Bush, J. W. M. Wavelike statistics from pilot-wave dynamics in a circular corral. Phys. Rev. E 88, 011001 (2013).
Perrard, S., Labousse, M., Fort, E. & Couder, Y. Chaos Driven by Interfering Memory. Phys. Rev. Lett. 113, 104101 (2014).
Perrard, S., Labousse, M., Miskin, M., Fort, E. & Couder, Y. Self-organization into quantized eigenstates of a classical wave-driven particle. Nat. Commun. 5, 3219 (2014).
Filoux, B., Hubert, M. & Vandewalle, N. Strings of droplets propelled by coherent waves. Phys. Rev. E 92, 041004 (2015).
Filoux, B., Hubert, M., Schlagheck, P. & Vandewalle, N. Walking droplets in linear channels. Phys. Rev. Fluids 2, 013601 (2017).
Hubert, M., Labousse, M. & Perrard, S. Self-propulsion and crossing statistics under random initial conditions. Phys. Rev. E 95, 062607 (2017).
Sáenz, P. J., Cristea-Platon, T. & Bush, J. W. M. Statistical projection effects in a hydrodynamic pilot-wave system. Nat. Phys. 14, 315–319 (2018).
Sáenz, P. J., Cristea-Platon, T. & Bush, J. W. M. A hydrodynamic analog of Friedel oscillations. Sci. Adv. 6, eaay9234 (2020).
Papatryfonos, K. et al. Hydrodynamic superradiance in wave-mediated cooperative tunneling. Commun Phys. 5, 142 (2022).
Hubert, M., Perrard, S., Labousse, M., Vandewalle, N. & Couder, Y. Tunable bimodal explorations of space from memory-driven deterministic dynamics. Phys. Rev. E 100, 032201 (2019).
Bacot, V., Perrard, S., Labousse, M., Couder, Y. & Fort, E. Multistable free states of an active particle from a coherent memory dynamics. Phys. Rev. Lett. 122, 104303 (2019).
Durey, M., Milewski, P. A. & Bush, J. W. M. Dynamics, emergent statistics, and the mean-pilot-wave potential of walking droplets. Chaos: Interdisciplinary J. Nonlinear Sci. 28, 096108 (2018).
Durey, M. Bifurcations and chaos in a lorenz-like pilot-wave system. Chaos:Interdisciplinary J. Nonlinear Sci. 30, 103115 (2020).
Durey, M. & Bush, J. W. M. Pilot-wave dynamics: the free particle. Chaos: Interdisciplinary J. Nonlinear Sci. 31, 033136 (2021).
Durey, M., Turton, S. & Bush, J. W. M. Speed oscillations in classical pilot-wave dynamics. Proc. Royal Soc. A 476, 2239 (2020).
Devauchelle, O., Lajeunesse, E., James, F., Josserand, C. & Lagrée, P. Walkers in a wave field with memory. Comptes Rendus. Mécanique 438, 591–611 (2020).
Valani, R. N., Slim, A. C., Paganin, D. M., Simula, T. P. & Vo, T. Unsteady dynamics of a classical particle-wave entity. Phys. Rev. E 104, 015106 (2021).
Berg, H. & Brown, D. Chemotaxis in escherichia coli analysed by three-dimensional tracking. Nature 239, 500–504 (1972).
Tailleur, J. & Cates, M. E. Statistical mechanics of interacting run-and-tumble bacteria. Phys. Rev. Lett. 100, 218103 (2008).
Patteson, A. E., Gopinath, A., Goulian, M. & Arratia, P. E. Running and tumbling with E. coli in polymeric solutions. Sci. Rep. 5, 15761 (2015).
Hokmabad, B. V. et al. Emergence of bimodal motility in active droplets. Phys. Rev. X 11, 011043 (2021).
Kolmakov, G. V. & Aranson, I. S. Superfluid swimmers. Phys. Rev. Res. 3, 013188 (2021).
Labousse, M., Perrard, S., Couder, Y. & Fort, E. Build-up of macroscopic eigenstates in a memory-based constrained system. N. J. Phys. 16, 113027 (2014).
Galeano-Rios, C. A., Milewski, P. A. & Vanden-Broeck, J.-M. Quasi-normal free-surface impacts, capillary rebounds and application to faraday walkers. J. Fluid Mech. 873, 856 (2019).
Perrard, S., Fort, E. & Couder, Y. Wave-based turing machine: time reversal and information erasing. Phys. Rev. Lett. 117, 094502 (2016).
Labousse, M. & Perrard, S. Non-Hamiltonian features of a classical pilot-wave dynamics. Phys. Rev. E 90, 022913 (2014).
Gilet, T. Dynamics and statistics of wave-particle interactions in a confined geometry. Phys. Rev. E 90, 052917 (2014).
Gilet, T. Quantumlike statistics of deterministic wave-particle interactions in a circular cavity. Phys. Rev. E 93 (2016).
Olver, F., Lozier, D., R.F., B. & Clark, C. NIST Handbook of Mathematical Functions (Cambridge University Press, 2010).
Romanczuk, P., Bär, M., Ebeling, W., Lindner, B. & Schimansky-Geier, L. Active Brownian particles: from individual to collective stochastic dynamics: From individual to collective stochastic dynamics. Eur. Phys. J.: Special Topics 202, 1–162 (2012).
Erdmann, U., Ebeling, W., Schimansky-Geier, L. & Schweitzer, F. Brownian particles far from equilibrium. Eur. Phys. J. B 15, 105–113 (2000).
Muller, L., Chavane, F., Reynolds, J. & Sejnowski, T. J. Cortical travelling waves: mechanisms and computational principles. Nat. Rev. Neurosci. 19, 255–268 (2018).
Labousse, M. Etude d’une dynamique à mémoire de chemin: une expérimentation théorique. PhD Thesis, Université Pierre et Marie Curie-Paris VI (2014).
Gilet, T., Vandewalle, N. & Dorbolo, S. Completely inelastic ball. Phys. Rev. E 79, 055201 (2009).
Acknowledgements
We would like to thank Yves Couder for past fruitful insights. The authors thank warmly Vincent Bacot and Emmanuel Fort for insightful discussions. M.H. acknowledges financial support by the Actions de Recherches Concertées (ARC) of the Belgium Wallonia-Brussels Federation under Contract No. 12-17/02. M.L. and S.P. acknowledge the financial support of the French Agence Nationale de la Recherche, through the project ANR Freeflow, LABEX WIFI (Laboratory of Excellence ANR-10-LABX-24), within the French Program Investments for the Future under reference ANR-10IDEX-0001-02 PSL. Computational resources have been provided by the Consortium des Équipements de Calcul Intensif (CÉCI), funded by the Fonds de la Recherche Scientifique de Belgique (F.R.S.-FNRS) under Grant No. 2.5020.11.
Author information
Authors and Affiliations
Contributions
M.H. performed the simulations, S.P. performed the experiments, N.V. and M.L. supervised the study. All authors participate in the redaction of the paper.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Carlos Antonio Galeano Ríos and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hubert, M., Perrard, S., Vandewalle, N. et al. Overload wave-memory induces amnesia of a self-propelled particle. Nat Commun 13, 4357 (2022). https://doi.org/10.1038/s41467-022-31736-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-31736-z