S Perrard

M Labousse

Transition to chaos in wave memory dynamics in a harmonic well : deterministic and noise-driven behaviour

A walker is the association of a sub-millimetric bouncing drop moving along with a co-evolving Faraday wave. When confined in a harmonic potential, its stable trajectories are periodic and quantised both in extension and mean angular momentum. In this article we present the rest of the story, specifically the chaotic paths. They are chaotic and show intermittent behaviours between unstable quantised set of attractors. First we present the two possible situations we find experimentally. Then we emphasise theoretically two mechanisms that lead to unstable situations. It corresponds either to noise-driven chaos or low-dimensional deterministic chaos. Finally we characterise experimentally each of these distinct situations. This article aims at presenting a comprehensive investigation of the unstable paths in order to complete the picture of walkers in a two dimensional harmonic potential.

A walker is the association of a bouncing drop and its associated guiding Faraday wave field. This dual system is one of the rare experimental situation encoding a pilot wave dynamics. It has shown multiple behaviours driven by temporal non local effects and reminiscent of one particle quantum-like behaviours. When confined by a harmonic potential, its stable solutions form a set of attractors that are quantised both in extension and mean angular momentum. This striking results has been observed experimentally and numerically by different groups. Here we investigate the case of unstable solutions presenting intermit behaviours between the footprint of the aforementioned set of attractors. We present here a detailed and unified picture of these chaotic trajectories. We show that at least two distinct mechanisms arise and investigate their theoretical origin. The first mechanism is the signature of a low-dimensional chaos as previously reported by several groups. The second mechanism is driven by a multiplicative noise encoded within the multiple wave degrees of freedom. Finally we characterise experimentally these two mechanisms and present a unified and synthetic picture of all intermittent transitions that have been documented so far.

I. INTRODUCTION

The transition from deterministic dynamics toward highly disordered, unpredictable dynamics has been thoroughly studied in physics. The dynamics of a walker [1], the association of a bouncing drop with the wave field it generates, also shows a transition from regular, periodic trajectories toward disordered states under confinement. The specificity of walker dynamics is that it exhibits wave-like behaviours and quantised sets of attractors both in the deterministic regime [2][3][4] or in the highly disordered regime [5,6]. The study of the transient from one regime to another has been less studied, but it has now been reported experimentally [7] and numerically [8]. Other studies have also reported the arising of more complex trajectories [9] that are located in this intermediate regime. The origin of the instability mechanisms driving the system from stable to disorder are still unclear. However the use of standard dynamical system tools (bifurcation diagram, first return map, period doubling research) have drawn a transition to low-dimensional chaos which implies a loss of determinism by increasing lack of predictability. Tambasco et al. [8] in particular have compared three confinement configurations (Coriolis force, central force and Coulomb potential) to show that several scenarii of transition to chaos can be observed, depending on the confinement type.

In the present article, we go beyond our first report of low dimensional chaos [7] both experimentally and theoretically. Experimentally, we base our analysis on data obtained with a central force field. We provide evidence of a transition to chaos from periodic orbits to unpredictable dynamics of two types. One type follows the previous observation of Perrard et al. [7] and the previous detailed analysis of Tambasco et al. [8]. The transition is in particular now thoroughly characterised and consistent between authors for small extension orbits. In contrast, for larger orbit size, the experiments show a total loss of determinism while the transition between reminiscent periodic trajectories are still deterministic. It is in favour of a noise amplification mechanism throughout information storage in the wave field. This article is organised as follows. In Sec. II, we present the experimental setup and summarise the main results. In Sec. III A, we derive theoretically a set of equation in the continuous limit, that enables us to compute the global Lyapunov exponent of the walker dynamics. The result is not specific to the central force field: we show that any volume of phase space globally contracts at all memory but it converges toward zero in the infinite memory limit. This constrain roots a route towards noise-driven chaos. Then we provide in Sec. III B an analysis of periodic orbits and show that a constraint on the mode amplitudes applies. It evidences the number of degree of freedom that are relevant to describe the dynamics, depending on the spatial extension of walker trajectory. Finally we characterise experimentally these two scenarios in Sec. IV. For the sake of simplicity, we start with low-dimensional chaotic type in Sec. IV A and emphasise the noise amplification mechanism in Sec. IV B.

II. EXPERIMENTAL SET UP AND EVIDENCE OF CHAOS

A. Methods

The experimental setup consists in a bath filled with silicon oil of viscosity ν = 20 cp shaken vertically at an acceleration γ = γ 0 cos(2πf t) at a frequency f = 80 Hz. The acceleration amplitude γ 0 of the bath is tuned in the vicinity below the Faraday acceleration threshold γ F 4.5 g above which standing Faraday waves are spontaneously generated [10][11][12]. The control parameter for the damping time of the drop generated waves is M = γ 0 /(γ Fγ 0 ) which ranges typically from 10 to 200. A sub-millimetric drop of mass m D bounces on the bath on a doubling period regimes at the Faraday period T F , as described in [13]. In this situation the drop is self-propelled by a dynamically [3]. c) Zoom in view on the mean spatial extension R of trajectories as a function of the confinement Λ for 50 < M < 100. Regions of stable periodic orbits (dark grey) and regions of intermittent behaviour (light grey). We focus here on the light grey regions. Figure 1b adapted from [3] with courtesy of Perrard et al.

co-moving localised standing waves of wavelength λ F = 4.5 mm [1]. The dynamics of the drop for increasing memory is explored in a confined situations, in which the drop is submitted to an external force. This force is generated by the use of magnetic fields. Two coils placed in a Helmholtz configuration around the liquid bath generate a homogeneous vertical magnetic field B 0 in the plane of the unperturbed liquid surface. Magnetic field gradients are generated with a cylindrical magnet made of an Neodyme Bore alloy placed at a tunable distance d of the liquid surface. The total magnetic field is therefore B = B 0 e z + B m where B m is the magnetic field generated by the magnet. The drop is made magnetic sensitive by encapsulating inside a droplet of ferrofluid of mass ≈0.05 m D . The potential energy of interaction writes E m ∝ -m.B where m ∝ B is the induced magnetic dipole. Our ferrofluid was previously characterised in [START_REF] Browaeys | Les ferrofluides : ondes de surface, résistance de vague et simulation de la convection dans le manteau terrestre[END_REF] and the entire magnetic set up calibrated using a magnet oscillating horizontally. The full experimental details are given in [3]. In the central region of interest (r < 3λ F ), we are left with a harmonic potential for the drop where the magnetic force per unit mass F writes :

F = -ω 2 r (1)
Here ω = k/m W is the magnetic frequency and r is the radial vector whose length is the distance to the center of the system. A correction to the drop inertia originating from the wave field has been identified experimentally in the central force set up as well as on a rotating bath [2,[START_REF] Oza | [END_REF] and theoretically rationalised [16]. The deviation from m W = m D was then measured throughout the radius R of the stable circular motion of the confined drop. Eventually, we calibrate the system directly on circular orbits. In the following we will use the non dimensional confinement parameter Λ = V 0 /(ωλ F ) that sets the orbit size where V 0 is the mean speed of the walker. For circular motions and no significant additional force generated by the wave field we observe R/λ F = Λ. The dynamics is further explored as a function of two main parameters. The memory parameter M (or memory time τ = M T F ) sets the temporal coherence of the system. It is tuned throughout the vertical acceleration of the path. The confinement parameter Λ sets indirectly the mean extension of the trajectories. Λ is varied by changing the distance d from the magnet to the bath. We then explore different type trajectories of increasing complexity by increasing Λ.

B. Stable & unstable periodic orbits in a harmonic potential

The stable trajectories generated by a walker in a harmonic well have been reported experimentally using the set up previously described [3,7], numerically using discrete time models [17,18] or continuous time model [8,19] and investigated theoretically [17][18][19][20]. The observation of stable trajectories of different symmetries is a common feature, robust to changes of the model details. A quantisation of mean radius R and mean angular momentum L for stable trajectories is also observed by all authors, although differences may arise between the models in particular for the exact quantisation selection rule [19] or the accuracy of the angular momentum quantisation [18]. We introduce here the reader to walker's dynamics in a harmonic well throughout the experimental data first reported in [3]. The mean radius of the stable orbits observed experimentally in the range 50 < M < 100 are represented in figure 1a as a function of the confinement parameter Λ. The colour codes the symmetry of the various trajectories: black for circular symmetry, red for two fold symmetry (lemniscate) and blue for 3-fold symmetry (trifolium). All the stable solutions that we found experimentally are classified in a diagram ( R, Lz ) (figure 1b) and each cluster of states is labelled accordingly by a set of two integers (n, m). However, these stable orbits are only one side of the full story. As the memory parameter M is increased, most parameter values Λ lead to unstable orbits with intermittent shift between trajectories of different symmetry. This feature has been reported experimentally [7] and thoroughly studied numerically by Tambasco et al. [8]. We propose here to describe the transition from stable to unstable orbit in a more systematic manner than it was previously reported. We will highlight in particular the two types of intermittency that have been experimentally observed.

C. Intermittent regimes :definition & first observation

The first destabilisation of a non periodic orbit emerges at intermediate memory (M ≈ 50) for highly confined trajectories (Λ ≈ 0.5). This regime corresponds to a dimensionless memory length S Me = V 0 τ /λ F of order unity. Figure 1c shows a close up of mean radius trajectory as a function of the confinement parameter Λ for the modes n = 1 and n = 2 for M = 50. The stable zones corresponding to simple periodic orbits have been highlighted in dark grey. We observe that between two stable regions, no stable orbits are observed. We will here detail the analysis for the observation of the first two regions of non periodic orbits. For this purpose, we define the critical values Λ ± (n,m) of the confinement parameter Λ that delimits the stability range of each mode (n, m). Each mode characterised by an integer couple (n, m) is therefore observed in the range Λ - (n,m) < Λ < Λ + (n,m) . The first zone of interest is located in Λ values between circular orbits n = 1 and lemniscate mode (n, m) = (2, 0), hence Λ + 1,1 < Λ < Λ - 2,0 . The second zone is located between lemniscate mode (n, m) = (2, 0) and circular orbits n = 2, hence Λ + 2,0 < Λ < Λ - 2,2 . The width of unstable zones increases with M , such that the limit of stability Λ ± n,m is a function of memory.

Four experimental trajectories obtained for increasing values of Λ are represented in figure 2 for M = 50. As the confinement parameter Λ increases the circular orbits become unstable. This phenomenon is better observed from a temporal record of the angular momentum L z , as illustrated in figure 2. It shows a first reversal of rotation direction after a long duration of periodical motion (typically 50 orbital periods T o ). For higher values of Λ, the average time spent in each direction of rotation decreases. The trajectory is then composed of gently diverging oscillations of angular momentum (henceforth laminar phase) and abrupt transitions in between (henceforth chaotic phase). In the laminar phase, the angular momentum Lz averaged on the period of the orbit is still equal to the unperturbed mode. Time intervals of mean angular momentum Lz = L (1) z has been coloured in yellow, while Lz = L (-1) z zones were coloured in green. For higher values of Λ, a third type of motion emerge corresponding to a large oscillation of L. The direct trajectory observation reveals that the transition is mediated by a third building block of motion corresponding to the lemniscate mode (red zone of fig 2). Yellow, green and red zones of variable length are then sequentially observed. The total time spend in this three types of motion covers most of the recording. The interval spent in one of these modes is called laminar phase while the short, erratic motion between these phases is called chaotic burst. It is already remarkable that a succession of several eigenmodes of different symmetries appears in an apparent erratic dynamics.

The arising of chaos from a stable limit cycle that destabilises spontaneously can be observed in numerous nonlinear dynamical systems as soon as the number of degree of freedom is equal or greater than 3 for continuous non linear systems [START_REF] Strogatz | Nonlinear dynamics and chaos : with applications to physics, biology and chemistry[END_REF]. In this route to chaos several scenario can be identified each of them being associated to different properties of the associated intermittency. They can be sorted out in particular by looking at the distribution of time spend in laminar phases. These routes are described in details in [START_REF] Bergé | Order within Chaos: Towards a Deterministic Approach to Turbulence[END_REF]. For the dynamics of a walker, the temporal records of angular momentum show similarity with intermittency of type II and III, in which duration of the laminar x/λ F 16 08 13 3-34 L=0.48616, V=9.5, M=74, w=4.9, R= 0.45 phases are unbound. For walkers dynamics it has been observed and characterised numerically by Tombasco et al. [8].

t(s) t(s) y/λ F y/λ F -1 0 1 L -1 0 1 L -1 0 1 L -1 0 1 L 50 100 150 t(s) t(s)
The appearance of non commensurate frequency may be associated to type III intermittency which is in practice rarely observed experimentally [START_REF] Manneville | Instabilités, Chaos et Turbulence[END_REF].

III. THEORETICAL ORIGIN OF WAVE SELF-ORGANISATION AND DIFFERENT ROUTE TOWARDS WAVE CHAOS.

In this Section, we emphasize the theoretical mechanisms leading to the two distinct chaotic transitions. A numerical investigation has been performed by Tombasco et al. [8] and indicates also these different routes towards chaos. We present here a complementary point of view by pointing out general mechanisms supporting these different classes of instabilities. We show that these mechanisms arise from fundamental reasons intrinsically due to the wave nature of the system. The first mechanism described in Sec. III A is due the way a wave stores information within an infinite number of degrees of freedom. The second mechanism described in Sec. III B is a direct consequence of the mathematical constraints arising from periodic orbits.

A. Phase space contraction and its relation to the memory parameter: a route towards a noise-driven chaotic scenario

Instabilities and chaotic transitions are usually characterised by analysing the rate of divergence in each direction of phase space or the Lyapunov spectrum {λ i } i∈N . In general this is a difficult task for integro-differential systems as is the walker dynamics where the definition of a phase space can be by itself a tedious task. First, we reformulate the equations of motion into a normal form and we elucidate the mathematical structure of the phase space. Then we show that it is possible to calculate the divergence rate of a small volume of phase space; We relate it to the memory parameter of the system. Finally we calculate the global Lyapunov λ. We emphasize important consequences on system stability originating from the discreteness of the system.

Normal form of the path memory dynamics

We first recall the equations of motion and then reformulate them into a local normal form Ż = F (Z) where Z is the state vector of the dynamics and F a local evolution function to specify. The dynamics is driven by a pilot-wave originating from the modulation of the interface. Its footprint is a standing surface h field co-evolving with the drop horizontal motion [2]

h = h 0 N k=-∞ J 0 (k F r N -r k ) e -(t N -t k )/τ , (2)
where h 0 the wave amplitude, k F = 2π/λ F is the Faraday wave vector, r N is the drop position at the N th impact at time t N . The indexes k < N indicate past situations and J 0 is the 0 order Bessel function of first kind. The concept of path memory dynamics has been introduced in its discrete form by Eddi et al. [START_REF] Eddi | [END_REF]. It has been investigated further in its discrete form by Milewski et al. [25] and in a continuous form by Oza et al. [26] we shall discuss later. The path memory is an essential feature at the origin of the temporal non-local properties of walker systems. An iterative discrete model based on this path memory has been benchmarked and so far has been able to reproduce quantitatively the experimental results [17]. In this article, for the sake of tractability, we take the continuous limit, meaning the above summation yields

h h 0 t -∞ dT T F J 0 (k F r(t) -r (T ) ) e -(t-T )/τ . (3)
We shall discuss in Sec. III A 2 the consequences of reaching the continuous limit. This continuous form has been proposed by Moláček et Bush [27,28] and yields

r = -μṙ -ω 2 r -C [∇h] r . (4)
r(t) denotes the drop position at time t. μ stands for an averaged dissipation friction per unit mass. -C [∇h] r is the force guiding the drop with a coupling coefficient per unit mass C. The external force -m W ω 2 r is applied by means of a magnetic field as described in Sec. II A. The hydrodynamic origin of each coefficient has been investigated in [28].

We now reformulate the equations of motion into a normal form. The main source of complexity is embedded within the nonlinear integral term as it encodes long-term correlations. A part of the information is encoded into the wave field and the rest of information is stored in the drop position and momentum. The aim of this paragraph is to clarify what means this interplay of information storage. Which dimensions store wave information? How is information and momentum exchanged between the wave and the drop? And eventually what sets the route toward chaos? Under the integro-differential form (Eq. 4), the dynamics appears as non-local in time. We first show that this equation can be reformulated as a local normal form Ż = F (Z) where Z is the state vector of the dynamics. This normal form will allow us in section III A 2 to compute the phase space divergence thanks to the partial derivatives of F. M. Miskin proposed in [3] to decompose the wave terms into a Bessel wave basis. The Graf's decomposition theorem enables one to project each "J 0 " terms of Eq 3 into a Bessel wave basis centred at the origin imposed by the external harmonic potential. The derivation of a local normal form describing the walker dynamics from Eq. 4 is detailed in appendix. We recall here its final form, obtained from a central mode decomposition of the wave field into the radial basis (r, θ). It yields

                       ẍ = -μ ẋ -Ω 2 x + Ck F h 0 n≥0 C n T n,x + S n V n,x ÿ = -μ ẏ -Ω 2 y + Ck F h 0 n≥0 C n V n,y -S n,y T n,y Ċn = -C n /τ + J n k F x 2 + y 2 T n (x, y) ∀n ≥ 0 Ṡn = -S n /τ + J n k F x 2 + y 2 V n (x, y) ∀n ≥ 1 , (5)
where the functions T n are the first kind Tchebychev polynomials that can be calculated by recurrence, and V n , T n,x and V n,x can be calculated from the T n and J n functions (see appendix for details). The mode amplitude C n and its phase quadrature S n describes the amplitude and the phase of the mode J n in the Bessel decomposition of the wave field. They are given by

               C n = N k=-∞ J n (k F r k ) e -(t N -t k )/τ cos (nθ k ) ≈ t -∞ dT T F J n (k F r(T )) e -(t-T )/τ cos (nθ(T )) S n = N k=-∞ J n (k F r k ) e -(t N -t k )/τ sin (nθ k ) ≈ t -∞ dT T F J n (k F r(T )) e -(t-T )/τ sin (nθ(T )) . (6)
The true interest of the rewriting stands in the clear separation between the evolution of the particle trajectory and its associated wave fields. We can now explicitly write a state vector

Z =             x ẋ y ẏ C 0 C 1 S 1 . . .             (7)
and remark that the dynamical system of Eq. 5 formally writes

Ż = F (Z) (8)
The memory information stored in the wave can be stored into a set of modes supporting the wave dimensions of the dynamics. Conceptually, the first four dimensions (x, ẋ, y, ẏ) corresponds to the particle phase space while the other dimensions (C 0 , C 1 , S 1 . . .) correspond to the wave counterpart. Note that the wave is responsible for supporting an arbitrary large number of additional degrees of freedom in the same spirit as the theoretical investigations of Labousse et al. [17,[START_REF] Labousse | Investigation of a path-memory dynamics : a theoretical trial[END_REF] and of T. Gilet [START_REF] Gilet | [END_REF]. The equation A.6 indicates a linear relation between [∂ x h] r and [∂ y h] r and C i et S i but a strong non-linear relation with the position coordinates. So the wave dimensions couple linearly to the particle dimensions. In contrast, the particle dimensions couple non-linearly to the wave dimensions. In this Section III A 1, we have translated the path memory integral into a set of modes evolution. Eq. 8 is still very complex, but this form of equations enables i) to define a phase space and ii) to calculate in the next section some very general properties of the system and its phase space contraction.

Phase space contraction and noise-driven chaos

In this section we aim at investigating theoretically the phase contraction properties. We consider a small volume of phase space

δV d = dxd ẋdyδ ẏdC 0 dC 1 dS 1 . . . dC d dS d = (2d+1)+4 j=1 dZ i (9)
We first consider an arbitrary large but finite dimension cut-off d so that Eq. 9 and every equations that follow have an unambiguous mathematical meaning. We will consider the limit d → +∞ at the end of the section. The divergence rate of this small phase space volume δV d towards a given direction is formally related to the Jacobian matrix

J i,j = ∂F i dZ j (10)
Specifically J i,j is the phase space deformation dZ j towards the direction of phase space i. The eigenvalues of J are directly related to the Lyapunov spectrum of the system {λ i } i=1,...,2d+5 . This information depends in general on the given position in the phase space Z and is not an invariant quantity of the dynamics. However we found that it is possible to compute an invariant of the system. Indeed, we recall that the time evolution of the phase space volume is given by [START_REF] Manneville | Instabilités, Chaos et Turbulence[END_REF][START_REF] Castiglione | Physique statisque, Chaos et approches multiéchelle[END_REF]

d dt δV d = (divF) δV d (11)
The divergence rate divF is the trace of the Jacobian matrix J and is base invariant. Its value does not depend of the choice of representation but could depend on position in the phase space Z. We find that is not the case and that divF is an invariant of the dynamics itself, specifically

divF = -2μ - 2d + 1 τ . (12)
Integrating Eq. 11 gives

δV d (t) = δV d (t = 0)e divF t ( 13
)
First divF is constant and always negative, implying that any initial volume of phase space contracts. Note that this does not prevent the existence of diverging directions of phase space, but this means that this divergence will be overcompensated by a contraction in another direction of phase space. Secondly, a contraction of phase space represents a loss of information. This loss has two distinct origins: a mechanical loss of information through the dissipation with the bath and ambient air at rate -2μ and a wave loss of information at rate -(2d + 1)/τ . Note that the loss of information per mode 1/τ is governed by the sole memory parameter.

Finally we calculate the global Lyapunov exponent λ = i λ i of the system. From any arbitrary infinitesimal volume δV d we can construct an arbitrary infinitesimal distance in the phase space δL as δL ∝ δV

1/D d , (14)
with the D = (2d + 1) + 4 the phase space of dimension. Combining Eqs 13 and 14 and then taking the limit d → +∞ we get

δL(t) = δL(t = 0)e -t/τ (15)
We may now identify the global Lyapunov exponent λ of the system as

λ = - 1 τ (16)
By taking the limit d → +∞ the wave dimensions dominates the scaling law (Eq. 14) and signifies that the phase contraction is mainly driven by the wave loss of information. Thus the wave dimensions act as an infinite number of additional degrees of freedom. We stress again that λ < 0 does not prevent from chaotic behaviours. Interestingly, for large memory parameter τ , λ converges towards 0 with negative values. This convergence to an asymptotic wave neutrality has direct consequences if the dynamics is discrete. Indeed the coarse-graining of the dynamics implies finite amplitude perturbation that can break the stability of an asymptotically neutrally stable attractor. The loss of stability in a phase space of large dimension triggers a route towards chaos. The driven mechanism is the encoding of the finite amplitude perturbation within the infinite wave degrees of freedom. This mechanism will be at play in Sec. IV B. We now turn to the existence of a different fundamental limit triggering a distinct instability mechanism.

B. General results for wave modes for periodic orbits: a route towards low-dimensional chaotic scenario

The existence of periodic solutions with central symmetry implies mathematical constraints on the dynamics we shall now investigate. These constraints may originate either from the dynamics or the geometry. We first exhibit these different constraints and show their interplay. Finally we show that these constraints lead to a non-linear mechanism promoting unstable situations.

Mechanical and geometrical constraints

We consider an orbit with a central symmetry of period T o . We first derive dynamical constraints. By integrating the equations of motion (Eq. 4) over one period of time between t 0 and t 0 +T o , we obtain the first mechanical constrain

T0 0 [∇h] r dt = 0 ( 17
)
For symmetry reasons, all the other terms disappear. For periodic orbits, the origin of time t 0 is irrelevant and is accordingly set to zero. For periodic orbits, the path integral is a summation of identical but progressively damped sequences. This geometry series simplifies the computation of the path-memory integral to one single orbit through the relation

h = h 0 t -∞ dT T F J 0 (k F r(t) -r (T ) ) e -(t-T )/τ = h 0 1 -e -To/τ t t-To dT T F J 0 (k F r(t) -r (T ) ) e -(t-T )/τ (18)
So a periodic trajectory must satisfy the mechanical constrain

To 0 ∇ r t t-To dT T F J 0 (k F r(t) -r (T ) ) e -(t-T )/τ = 0 (19)
We now exhibit geometrical constraint arising from the periodicity of the paths. After a transient the existence of periodic trajectories imposes that

h(t) -h(t -T o ) = h 0 T F (t), (20)
with (t) → 0 after a typical period of time corresponding to the transient. This means that

t t-To dT T F J 0 (k F r(t) -r (T ) ) e -(t-T )/τ = (t) (21)
We find that the geometrical constrain (Eq. 21) implies that the mechanical constrain (Eq. 19) is always satisfied at long time. Conversely Eq 19 requires periodic motion and a radial force. As a contraposition, any non radial external force prevents periodic trajectories. Thus the geometrical constrain (Eq. 21) is more demanding than the mechanical one (Eq. 19) and will be a preferred starting point for pursuing the calculation.

Incompatible wave mode constraints and instability

Expanding the field in Eq. 21 into a Bessel basis leads to

n∈Z J n (k F r(t)) e inθ(t) t t-To dT J n (k F r (T )) e -inθ(T ) e -(t-T )/τ = (t) (22)
For long memory we get

n∈Z J n (k F r(t)) e inθ(t) t t-To dT J n (k F r (T )) e -inθ(T ) = (t) + O T 0 τ (23)
We denote T ) . We integrate all the terms between t -T o and t and by symmetry of the integration domains, we get

A n = t t-To dT J n (k F r (T )) e -inθ(
1 2 n∈Z A n A * n = ˜ (t), (24)
with ˜ (t) arising from the integration of the right-hand side of Eq. 24 and converging to 0 for large time and large memory. Now consider a path invariant under a rotation of angle 2π/N as observed experimentally. We have A n = 0 if n/N is not an integer. So the summation reduces to the modes symmetrically compatible specifically

1 2 p∈Z |A p×N | 2 = ˜ (t) (25)
As thoroughly investigated for circular paths [START_REF] Oza | [END_REF]20], Eq. 25 means that the mode 0 extinguishes, |A 0 | 0 at high memory. For higher order symmetries, say a lemniscate (N=2), Eq. 25 implies that |A 0 | 0 and {|A ±2k | 0} k simultaneously. The size quantisation of the higher order solutions of the dynamics originates in this simultaneous mode extinction. The rather limited number of constraints involved supports a low-dimensional mechanisms as described in Sec. IV A. However, at high memory, as ˜ → 0, these constraints are difficult to be simultaneously satisfied at all time, which leads to low-dimensional instabilities. We now turn back to the experimental characterisation of these different scenarios.

IV. EXPERIMENTAL CHARACTERISATION OF THE DIFFERENT CHAOTIC SCENARIO

In this section, we seek for deterministic behaviour within the intermittent trajectories obtained experimentally in a harmonic well. A unified picture of the transition to chaos in wave-memory dynamics in confined situation has not been achieved. Instead, two regimes qualitatively different have been observed. We will present the two archetypes of intermittency types that have been observed for various spatial extension of the trajectories and different symmetry of the modes involved. It concerns the intermittency within the single lemniscate mode (n, m) = (2, 0), the intermittency between circular orbits n = 1 and the lemniscate mode (n, m) = (2, 0), and eventually the intermittency between circular orbits at n = 2.

To look for order in disorder, one canonical method is the use of first return map. The main idea is to reduce the description of a dynamics of a continuous variable in time X(t) to a dynamics of a variable discrete in time X k [START_REF] Floquet | [END_REF]. We then look for an expression of X k+1 as a function of X k . For a quasi-periodic system in which the main temporal evolution is an oscillation at period T with additional slow variation in time one can naturally chose X k = X(kT ). From an experimental point of view, the period T may be ill-defined. In practice, the choice of the local extremum of X(t) for the X k is more suitable to identify deterministic behaviour [START_REF] Abarbanel | Analysis of Observed Chaotic Data[END_REF]. The time interval T between X k and X k+1 can therefore vary with time as it is self-determined by the duration between two successive extrema. We then seek to describe the time serie X k by the mean of a function G such that :

X k+1 = G(X k ) ( 26
)
If G exists it is called the application of first return. This iterative point of view has been widely used in the 70's and 80's to describe the temporal evolution of highly disordered temporal signals [START_REF] Eckmann | [END_REF]. The first experimental use of the first return map has been used in particular on a model of thermal convection by Lorenz who identified a chaotic behaviour in a system of low dimensionality [35]. For sake of simplicity we first analyse the low-dimension chaos type. We build a first return map from the temporal record of the distance r to the center. Figure 3 shows two realisations computed from a lemniscate (figures 3a and b) and from an intermittency between circular orbits n = 1 and lemniscate (n = 2) mode (figures 3c andd). The temporal record of the angular momentum has been added to evidence the transition between modes of different symmetry. We first focus on the case of a single mode shown in figures 3a and b. It was obtained with a confinement Λ = 0.57 so in the vicinity of the lemniscate stability range. The first return map has been computed from the successive maximum R k of the distance r to the center. Each black dot corresponds to a measure while the red circles were obtained from an average of all dots in the range R and R + δR with δR = 0.005λ F . As we can see on figure 3b, the experimental realisation of G is single valued : one can indeed describes the dynamics with a deterministic function G. The fixed points of the map G exp namely M 1,±1 , M 2,0 and M 3,±1 correspond to attractors of the dynamics. In the current case, none of them are stable. Most of the dynamics occurs in the vicinity of M 2,0 , but excursions near M 1,±1 and M 3,±1 have also been observed which correspond in real space of loops of variable extensions and orientations. As a consequence, the apparent signal is disordered even if there exists an underlying deterministic evolution.

2. Orbits n = 1 : deterministic chaos, type III.

The intermittency near the stability range of n = 1 orbits can also be studied using the first return map although it exhibits qualitative difference. Figure 3d extracted from [7] shows the first return map built from the successive maximum R k of the distance to the center for Λ = 0.49 and M = 63. Each grey dot corresponds to a measure of one period of motion. In a similar manner, it shows the existence of three fixed points M 1,±1 , M 2,0 and M 3,±1 along the line R k+1 = R k . The associated trajectory in the vicinity of these points corresponds respectively to the modes (1, ±1), (2, 0) and (3, ±1). The striking aspect of this first return map is the presence of two branches between M 1,±1 and M 2,0 located on each side of the line R k+1 = R k . These two branches can however be distinguished by adding the previous maximum radius R k-1 as a variable. The conditional average on R k > R k1 gives the red dots curve, while the conditional average on R k < R k-1 gives the blue dots curve. The red branch of our first return map is reminiscent of the Lorenz attractor first return map [35] or the Rössler attractor [36]. The high slope of G near its maximum makes this dynamics highly sensitive to the initial conditions. If one can predict the temporal evolution starting from any initial condition this predictability is in practice lost after a finite time.

The presence of two branches in this application breaks the deterministic evolution of a single variable R k . This predictability can be recovered by adding the new variable R k-1 . However, this choice introduces temporal correlation in the dynamics and thus does not correspond to a first return map. In order to build a two dimensional first return map without invoking temporal correlations, one has to look for a new relevant variable Y k describing the system state at instant k. Neither the angular position of the speed nor the drop position turn out to be relevant. In order to keep a description local in time of the instant k + 1 as a function of instant k, one has to add a variable which is not associated v or r. One way to rationalise the evolution is to choose for Y a variable associated to the wave field surrounding the drop. Following the Graf's theorem decomposition introduced in section III, a natural choice could be the amplitude A 0 of the centred mode J 0 .

Eventually, there exists two routes to rationalise this multivalued first return map. The first one includes variables associated to the wave field. Doing so, we introduce a field variable for the motion of a localised object. The dynamics is then no more equivalent to the motion of a localised particle in 2 dimensions of space. The second rationalising path adds a variable of the previous time R k-1 . It is still described solely by the motion of a localised particle in space, but the temporal locality is lost. In that sense, the multivalued first return map of figure 3 is a signature of the wave memory, that makes this dynamics either non local in time or in space.

B. Experimental characterisation of effective noise-driven chaos

The next parameter range that exhibits intermittency behaviour is located between the stability range of the lemniscate mode (n, m) = (2, 0) and the circular modes (n, m) = (2, ±2) for 0.65 < Λ < 0.8. Figure 4 shows four examples of trajectories obtained within this range for a memory parameter M = 50 and increasing values of Λ. The temporal record of the instantaneous angular momentum L z also shows evidence of a dynamics switching between values that correspond to eigenmodes (n, m) of the dynamics (dashed lines). The time spent in any of these modes can last hundreds of periods mediated by fast erratic transitions. Contrary to the previous cases around n = 1 and (n, m) = (2, 0), no slow growth of a wobbling amplitude have been observed. The first return map have been computed for the distance to the center R, the orientation θ, the angular moment L z and the amplitudes of the first modes of the wave field Bessel decomposition. None of them have revealed a deterministic function G. These maps only evidence an accumulation of points around unstable fixed points, without deterministic link between successive events. We did not succeed in revealing a global underlying deterministic behaviour.

The failure of the deterministic description can originate from an increase of the number of relevant degrees of freedom necessary to describe the dynamics. Indeed without any apriori symmetry consideration, the number of excited mode increases (as a square roots) with the orbit size. This increase of relevant dimensions is known to be an obstacle for the search of determinism in experimental realisations of a chaotic system [START_REF] Manneville | Instabilités, Chaos et Turbulence[END_REF]. However, even if the global deterministic description is apparently lost, the transition between modes still exhibit reproducible patterns. The two archetype transitions for (2, 0) → (2, 2) and (2, 2) → (2, 0) are shown in figure 5. From a lemniscate mode, an excursion of large amplitude generates one loop of larger radius of curvature that stay away from the center (red line of figure 5a). Once one large loop has been achieved, the drop keeps moving along a stadium shape characteristic of (n, m) = (2, ±2) modes. The transition back from stadium to lemniscate involves a transient two successive loops of small radius of curvature passing by the centre of the harmonic well (figure 5b). Once a loop of opposite angular momentum has been drawn, the drop keeps moving in a lemniscate mode until it will destabilise again. Due to the iterative characteristic of the walker's dynamics, the observation of these trajectories presents another interest. As can be observed from Eqs. 2&4, the drop trajectory encodes the entire information necessary to reconstruct the dynamics. The two transitions of figure 5 are therefore remarkable examples of representation of a bistable dynamics in real space.

A dynamics composed of an erratic evolution near the attractors with fast deterministic behaviour between at- tractors can be reminiscent of low-dimensional systems coupled to an external noise. It is in particular the case for magnetic field earth reversal in which the magnetic field is submitted to a multiplicative noise throughout the coupling with the turbulent velocity field [37,38]. Depending on the nature of stable and unstable fixed points, low-dimensional systems coupled to an external noise can also be associated to ON-OFF intermittency [39,40]. The low number of relevant dimensions, the reproducibility of transient regimes and the existence of a multiplicative noise are common features of these system. We will see that a multiplicative noise is also found in walker's dynamics.

In the walker dynamics, there is a source of multiplicative noise which amplitude may grow with the memory parameter. The temporal variability of the impact time and impact positions generate a fluctuation of the wave amplitude h 0 (1 + η k ) where η k 1 is the noise amplitude. The total wave force exerted on the drop writes :

F waves ∝ -∇ r N N k=-∞ h 0 (1 + η k )J 0 (k F ||r N -r k ||)e -(t N -t k )/τ (27)
If we consider e.g. the η k as independent random variables with η k = 0 and η 2 k = η 2 , the central limit theorem applies and the total fluctuating part of the wave force δF waves writes : Even if the noises η k on each individual source is of negligible amplitude, the temporal recording in the wave field accumulate noise effects. Numerical simulations performed with noise of increasing amplitude have previously revealed its key influence on the dynamics [41]. The temporal recording of noise effects in the wave field and its influence on the dynamics is still an open question that would need further investigations. The accumulation of noises and its multiplicative nature, as well as the temporal record of angular momentum shows strong evidences of an analogy with multi-stability cases driven by multiplicative noises. The existence of two modes of same symmetry (n, m) = (2, ±2) whose transitions are mediated by a more unstable mode of different symmetry is reminiscent of the magnetic earth reversal [42].

δF waves ≈ η √ M F waves (28)
When one observes more extended trajectories for increasing values of Λ, the complexity of the trajectories and the number of stable modes involved keeps increasing. As a result, no deterministic behaviour has been identified. A route to chaos as described numerically by Tambasco et al. [8] is more likely to happen. In this particular case, the absence of frequency doubling cascade in the presence of a harmonic well is an argument for the Ruelle-Takens-Newhouse scenario [43,44]. However, the existence of a multiplicative noise of increasing amplitude with memory may limit, in practice, the observation range of deterministic chaos in the experiments performed with walkers. The influence of noise effects on the dynamics have been directly pointed out [41]. Noise amplification will always eventually arise due to the amplification mechanism we present here, but the memory limit at which it dominates the dynamics depends on setups and initial noise level to be amplified. A direct study of noises encoding and its influence on the dynamics is still to explored.

As a summary to our experimental investigation, we show in figure 6 the diagram of all observed transitions between modes. The mean angular momentum L is represented as a function of the mean spatial extension R for all the stable trajectories observed experimentally in a harmonic well. Reminiscent long loops observed in chaotic regimes have also been added on this diagram as they mediate the transition between the other modes (n, m). The possible transitions are represented with left and right arrows. Solid lines are documented archetypes while dashed lined would need further experiments to be characterised. It is remarkable to note that all transitions are observed to be local in this representation. They involved either a shift of symmetry m → m ± 2 (red lines) or a shift of symmetry and spatial extension (m, n) → (n ± 1, m ± 1) (blue lines).

V. CONCLUSION

In this article we have presented the archetypes of transitions to chaos in walker dynamics under central force confinement, and we have discussed the different scenario leading to disordered trajectories. Section II has presented how to obtain a central force field confinement for a bouncing drop, and present the first sign of instabilities we observed experimentally. Then in Sec. III we have presented a theoretical investigation of walker dynamics stability in a central force field. We emphasize in particular two possible nonlinear mechanisms originating the transition to disordered trajectories. The first one is based on an analysis of the phase space contraction in the limit of large memory. The global Lyapunov exponent λ converges toward 0 from negative values at high memory. Moreover this exponent is dominated by the wave field that contains most of the relevant dimensions. In the high memory limit, any finite size perturbations lead to instability and may trigger chaotic behaviours. In that scenario the discreteness of the dynamics that is encoded within the waves degrees of freedom acts as a noise source and may play a crucial role. The second scenario is based on an analysis of the constraints imposed on the wave field by periodic trajectories. In the high memory limit, we show that the wave mode amplitude must all vanish for ensuring periodic motion. For orbits of small extension it implies only a limited number of relevant constraint. In that regime, the walker dynamics may be described by a small number of degrees of freedom and thus a low dimensional chaotic behaviour is expected.

In section IV, we present the experimental characterisation of disordered trajectories. Depending on spatial extension trajectories, two scenario have been identified. The first scenario is a low dimensional chaotic behaviour that can be rationalised using first return map. This map can be defined either from the drop position only (intermittency of lemniscate mode) or from a combination of drop position variables and degrees of freedom associated to the wave field (intermittency between n = 1 modes and lemniscate). It is in accordance with a low dimensional chaos intuited in section III. The second scenario is a loss of full determinism throughout noise amplification. The absence of experimental deterministic first return map for higher extension modes (Λ > 0.7) is in favour of this second scenario. In that scenario, the loss of stability is triggered by noise fluctuations. Once the motion has left a stable region of the phase space, a fast and deterministic transition occurs which brings the dynamics in the vicinity of another attractor of the phase space. The noise amplification mechanism detailed in Eqs. 27 and 28 is a key ingredient for understanding walker dynamics in the high memory limit. Finally, we present a diagram of all the observed transitions which shows that the symmetry change during a transition are local, i.e. only local shift of symmetry m → m ± 2 or local shift of symmetry and spatial extension (m, n) → (n ± 1, m ± 1) have been observed. The emergence of selection rules for confined walkers have been already evidenced in Coriolis force experiments [2,[START_REF] Oza | [END_REF] or in central force configuration [3,18,19]. We highlight here few mechanisms of transitions from stable orbits toward disordered states from an experimental and a theoretical point of view. It is remarkable that a macroscopic wave particle exhibits also selection rules for transitions between eigenstates.

The functions T n are the first kind Tchebychev polynomials that can be calculated by recurrence. Finally the integrodifferential equation can be rewritten into the desired local normal form

                       ẍ = -μ ẋ -Ω 2 x + Ck F h 0 n≥0 C n T n,x + S n V n,x ÿ = -μ ẏ -Ω 2 y + Ck F h 0 n≥0
C n V n,y -S n,y T n,y Ċn = -C n /τ + J n k F x 2 + y 2 T n (x, y) ∀n ≥ 0 Ṡn = -S n /τ + J n k F x 2 + y 2 V n (x, y) ∀n ≥ 1 (A.9)

FIG. 1 .

 1 FIG. 1. Experimental set up and region of interest. a) Sketch of the experimental set up. A drop of silicon oil encapsulating ferrofluid (zoomed in view) is generated on a vertically vibrating liquid bath. Immersed in a magnetic field B generated by Helmholtz coils (B0 field) and a cylindrical magnet (Bm field) the drop is submitted to a central attractive harmonic force. The drop exhibit b) Diagram of spatial extension R and angular momentum L computed on one period of motion for all observed periodic trajectories from[3]. c) Zoom in view on the mean spatial extension R of trajectories as a function of the confinement Λ for 50 < M < 100. Regions of stable periodic orbits (dark grey) and regions of intermittent behaviour (light grey). We focus here on the light grey regions. Figure1badapted from[3] with courtesy of Perrard et al.

Figure 5 . 2 -FIG. 2 .

 522 Figure 5.2 -Déstabilisation d'une orbite stable. Pour un paramètre de mémoire donné Me = 63, la valeur du confinement ⇤ est augmentée depuis une zone d'orbite circulaire n = 1 stable (⇤ = 0.46) vers une zone instable. a) ⇤ = 0.46. b) ⇤ = 0.49. c) ⇤ = 0.5. d) ⇤ = 0.62. Ces trajectoires ont été obtenues avec une goutte de vitesse V = 10±0.5 mm/s Les domaines correspondant à des modes différents ont été coloriées en jaune pour le mode (n, m) = (1, 1), en vert pour (n, m) = (1, 1) et en rouge pour (n, m) = (2, 0). Une apparition de comportements intermittents est observée.

FIG. 3 .

 3 FIG. 3. First return map a) Temporal record of dimensionless angular momentum L obtained for Λ = 0.57 and V = 8.6 mm/s for M = 100. b) Associated first return map (single valued case). The local maximum R k+1 is represented as a function of the previous maximum R k . Ensemble average on R k values for bins of width δR = 0.005 (•). Fix point M2,0 of the first return map defined by the intersection with R k+1 = R k (solid line). c) Temporal record of L for Λ = 0.49 and V = 9.7 mm/s for M = 63 using also the local maximum R k . d) Associated first return map (bi-valued case) for R k . Ensemble average on R k values for bins of width δR = 0.005 (•) for R k > R k-1 (•) and R k < R k-1 (•). The two branches are well separated between A and B, and mixed between B and C. Fix point M1,±1 and M2,0 of the first return map define by the intersection with R k+1 = R k (solid line). Fig 3c & 3c adapted from [7] with courtesy of Perrard [? ]

FIG. 4 .

 4 FIG. 4. Intermittency of circular orbits n = 2 Trajectory and angular momentum L for increasing values of Λ, M = 50 and V = 11.5±0.3 mm/s. a) Λ = 0.73. b) Λ = 0.75. c) Λ = 0.76. d) Λ = 0.77.

FIG. 5 .

 5 FIG. 5. Archetype of transition between modes n = 2. Trajectory and angular momentum near a spontaneous transition at the onset of intermittency (M = 25, Λ = 0.62 and V = 7.8 mm/s). Both transitions are extracted from the same recording.Trajectory before the transition (blue) during the transition (red) and after the transition (black). a) Transition (n, m) = (2, 0) → (2, -2). b) Transition (n, m) = (2, -2) → (2, 0).

FIG. 6 .

 6 FIG.6. Transitions between levels : an overview. Diagram in the representation ( R, L) of stable trajectories in the shape of circle (black circles), lemniscates (red squares) and trifolium (blue triangles) for M > 50. Unstable trajectories that mediate the transition between modes have been added. Arrows represent the documented spontaneous transitions between modes observed for M > 50.

VI. ACKNOWLEDGEMENTS

The authors thank M. Berhanu, Y. Couder, E. Fort and F. Petrelis for fruitful discussions. S.P. and M.L. acknowledge LASIPS program, 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-10-IDEX-0001-02 PSL.

Appendix: Derivation of the ODEs for pilot-wave dynamics

We detail here the derivation of Eq. 5 from the integro-differential form given by Eq. 4. The Graf's decomposition theorem enables one to project each "J 0 " terms of (Eq 3) into a Bessel wave basis centred at the origin imposed by the external harmonic potential. This central mode decomposition into the radial basis (r, θ) yields

dT TF Jn (kF r (T )) e -(t-T )/τ e in(θ(t)-θ(T )) , (A.1)

and the surface field rewrite

with δ n,0 the Kronecker symbol. The mode C n and its phase quadrature S n are given by

The true interest of the rewriting stands in the link between C n (or S n ) and its time derivative

This decomposition in a central wave basis is the first but crucial step toward writing the dynamics in a form local in time. This set of equation means that the modes C n and S n support all the waves information. The last step of reformulation consists in expressing the equation of motions into a Cartesian basis

The x and y component of the wave force [∂ x h] r and [∂ y h] r are given by