1 Introduction
Pilotwave hydrodynamics as a field of study was initiated in 2005 with the discovery by Couder et al. [Reference Couder, Protière, Fort and Boudaoud3] that a droplet may selfpropel along the surface of a vibrating bath of fluid, guided by the waves it generates on each impact. Moláček and Bush [Reference Moláček and Bush12] developed a hydrodynamically consistent equation of motion for this walking droplet (or ‘walker’) of mass m given by nondimensional position $\mathbf {x}_p$ governed by its real and analytic wave field $h(\mathbf {x},t)$ . By time averaging over the bouncing period, the vertical dynamics are eliminated from consideration which produces the trajectory equation
such that $\mathcal {F}=\mathcal {F}(\mathbf {\dot x}_p)$ is an external force and
where $F=mgAk_F$ , D is the drag coefficient, $T_F$ is the Faraday period, $T_d$ is the decay time of waves without forcing, $k_F$ is the Faraday wavenumber, $\gamma $ is forcing acceleration, $\gamma _F$ is the Faraday instability threshold, A is the amplitude of a single surface wave and t is nondimensional time (see [Reference Oza, Harris, Rosales and Bush14]). The guiding potential h is approximated by an integral over the particle’s prior trajectory:
such that W is real and entire. Writing $\mathbf {x}_p(t)$ as $x(t)\mathbf {j}$ , for an arbitrary horizontal unit vector $\mathbf {j}$ , we arrive at a dimensionless integrodifferential equation:
For vibration amplitudes close to but below the Faraday threshold, waves created by the walker on each bounce extend far in space and decay slowly in time. In this regime, predicting the walker’s future dynamics not only requires knowledge of its present state but also of its past, creating memory in the system and nonlocality in time. Durey [Reference Durey4] and Valani et al. [Reference Valani, Slim, Paganin, Samula and Vo16] investigated the highmemory regime using an idealised pilotwave model that implements a sinusoidal waveform for the waves generated by the droplet and discovered an equivalence in orbital dynamics to the Lorenz system.
For a twodimensional pilotwave system in a rotating frame, Oza et al. [Reference Oza, Harris, Rosales and Bush14] demonstrated quantisation of orbitally unstable and chaotic regions. The differential Galois criterion for partial integrability can be interpreted as a measure of orbital complexity (in the sense that complexity is determined by whether a solution to a given dynamical system $\mathcal {S}$ is expressible in terms of first integrals). So, it is natural that this criterion is connected to classifying whether or not the dynamics of $\mathcal {S}$ are chaotic in its parameter space.
The main methods used in the present paper are due to Ruiz’s applications of differential Galois theory [Reference MoralesRuiz13] and Bogoyavlensky’s extension of Liouvillian integrability to dynamical systems [Reference Bogoyavlenskij2]. In Section 2, the equivalence of (1.1) to the Lorenz system and its linear stability is provided. In Section 3, we shall formulate a criterion for partial integrability (specifically, ‘Bintegrability’ as defined in [Reference Bogoyavlenskij2]), which is the main result.
2 The trajectory equation as a nonlinear differential system in $\mathbb {R}^3$
The following theorem expresses (1.1) as an infinite dimensional system of nonlinear ODEs. Referring to [Reference Valani15, Reference Valani, Slim, Paganin, Samula and Vo16], this form of (1.1) is more tractable and, in special cases, such as when considering a waveform in the absence of spatial decay, the infinite system reduces to a finite dimensional system.
Theorem 2.1. Let
for all $n\in \mathbb {N}$ . Equation (1.1) holds if and only if
Proof. The first two linear ODEs are easily verified, since (1.1) may be expressed as
Given that the integrand of $M_n$ is analytic, we may apply the Leibniz integral rule. Observe that for all $n\geq 0$ ,
as claimed.
2.1 The linear stability of a sinusoidal waveform system
According to Durey [Reference Durey4], we may define $W(s)=\cos (s)$ for the high memory regime, such that letting $\sigma =1/\kappa $ , $r=\beta $ and $b=1$ , (2.1) is reduced to the vector field $V:\mathbb {R}^3\rightarrow \mathbb {R}^3$ , defined by
where $X=v$ , $Y=\beta M_0$ , $Z=\beta (1M_1)$ and $\mathbf {P}=(X,Y,Z)$ . For all $\mathcal {F}$ such that $\mathcal {F}(0)=0$ , an equilibrium point occurs at $\mathbf {P}^*=(0,0,0)$ and the Jacobian matrix of V at $\mathbf {P}^*$ , letting $\mathcal {F}'(X)=\partial \mathcal {F}/\partial X$ , is
So, we now consider the linearised system to first order in $\epsilon \in \mathbb {R}$ , where $\mathbf {P}=\mathbf {P}^*+\epsilon \mathbf {u}$ :
Setting $\det (\lambda IJV_{\mathbf {P}^*})=0$ for $\lambda \in \mathbb {R}$ yields the characteristic equation
Hence, we have $\lambda _0=b$ and may derive the remaining eigenvalues explicitly:
This explicit formula for eigenvalues is useful when calculating the linearised approximation $\mathbf {P}_t\approx \sum _{i=0}^2 C_i e^{\lambda _i t} \mathbf {v}_i$ for the stable region given in Theorem 2.3, where $\mathbf {v}_i$ is the corresponding eigenvector to $\lambda _i$ and $C_i\in \mathbb {R}$ . Note that the accuracy of the approximation is determined by how close the initial conditions are to $\mathbf {P}^*$ , which follows from Lemma 2.2.
Lemma 2.2 (Hartman–Grobman theorem, [Reference Hartman6])
If all the eigenvalues $\lambda _j$ of $JV_{\mathbf {P}^*}$ have nonzero real part, then the nonlinear flow is topologically conjugate to the flow of the linearised system in a neighbourhood of $\mathbf {P}^*$ . If $\Re (\lambda _j)>0$ for at least one j, the nonlinear flow is asymptotically unstable. If $\Re (\lambda _j)<0$ for all j, the nonlinear flow is asymptotically stable.
Theorem 2.3. The asymptotically stable region S of (2.2) is
such that $\sigma (\mathcal {F}'(0)1)<1$ and $r,\sigma \geq 0$ . Letting $\zeta $ denote the stability boundary, the remaining asymptotically unstable region $U:=\mathbb {R}^3\setminus (S\cup \zeta )$ satisfies $r>1\mathcal {F}'(0)$ .
Proof. If $\lambda =i\omega $ is a root with $\omega \in \mathbb {R}$ , then according to (2.4),
Separately considering real and imaginary parts, the stability boundary satisfies
Suppose that $\sigma (\mathcal {F}'(0)1)<1$ . For $\omega \lessgtr \pm \sqrt {\sigma (1\mathcal {F}'(0)r)},$ we must have each real part of $\lambda _{1,2}$ less than zero, which is the criterion for asymptotic stability according to Lemma 2.2 (the converse applies for asymptotic instability). Hence, simplifying the imaginary part of (2.5), we have
3 Nonlinear differential Galois theory and an integrability criterion
We will present an argument for algebraic Bintegrability introduced by MoralesRuiz [Reference MoralesRuiz13] and subsequently developed by Huang et al. [Reference Huang, Shi and Li7] for its application to the Lorenz system. From (2.2), we consider the system
such that $\mathbf {P}\in \mathcal {M}$ is an $(n=3)$ dimensional complex analytic manifold with $t\in \mathbb {C}$ .
Definition 3.1. System (3.1) is completely integrable if it possesses $n1$ functionally independent first integrals $\Phi _1,\Phi _2$ .
Definition 3.2. System (3.1) is Bintegrable if it possesses k functionally independent first integrals $\Phi _1,\ldots ,\Phi _k$ and $(nk)$ vector fields $w_1=V,\ldots ,w_{nk}$ such that
with $1\leq i\leq k\leq n$ and $1\leq j \leq nk$ . Note that $[\cdot ,\cdot ]:\mathfrak {g}\times \mathfrak {g}\rightarrow \mathfrak {g}$ denotes the Lie bracket for a vector space $\mathfrak {g}$ .
Bintegrability was introduced by Bogoyavlenskij [Reference Bogoyavlenskij2] and it was shown that if a system is Bintegrable, then it is integrable by quadrature. Intuitively, Bintegrability is a generalisation of Liouvillian integrability from Hamiltonian systems to dynamical systems (see [Reference Bogoyavlenskij2]). Referring to Llibre et al. [Reference Llibre, Valls and Zhang11], a completely integrable system is orbitally equivalent to a linear differential system.
Let $\hat {\mathbf {P}}=\hat {\mathbf {P}}(t)$ denote a nonequilibrium solution of (3.1). A form of the variational equation was briefly introduced as (2.3). For the known particular solution, the variational equation along the phase curve $\Gamma $ of $\hat {\mathbf {P}}$ is
such that $\nu \in T_{\Gamma }\mathcal {M}$ and $T_{\Gamma }\mathcal {M}$ is a vector bundle of $T \mathcal {M}$ restricted on $\Gamma $ . With the normal bundle $N=T_{\Gamma } \mathcal {M}/\Gamma $ , a natural projection $\pi : T_{\Gamma }\mathcal {M}\rightarrow N$ and $\eta \in N$ , the variational equations may be reduced to
The differential Galois group of (3.2) can be defined as a matrix group G, with $G\subset \mathrm{GL} (n1,\mathbb {C})$ acting on the fundamental solutions of (3.2) such that it does not change polynomial and differential relations between them.
The following theorem is a powerful result which extends Galoisian obstructions to mermorphic integrability from Hamiltonian systems, provided in [Reference MoralesRuiz13], to the more general nonHamiltonian case.
Theorem 3.3 [Reference Ayoul and Zung1]
System (3.1) is Bintegrable, in the meromorphic category, in a neighbourhood of $\Gamma $ if and only if the identity element $G^0$ of the differential Galois group of the normal variational equations along $\Gamma $ is abelian.
The following theorem provides necessary conditions for (3.1) to possess a certain number of first integrals. In the proof of Theorem 3.7, we will show that (3.1) violates condition (1) and provide a criterion for when condition (3) is violated.
Theorem 3.4 [Reference Li and Shi10]
If (3.1) has m functionally independent meromorphic first integrals in a neighbourhood of $\Gamma $ , then the Lie algebra $\mathcal {G}$ of the differential Galois group G of (3.2) has m meromorphic invariants and $G^0$ has at most $(nm1)(n1)$ generators. Hence:

(1) if $m=2$ , then (3.1) is completely integrable, $\mathcal {G}=0$ and $G^0=\{\mathbf {1}\}$ , which is the identity element;

(2) if $m=1$ , then $\mathcal {G}$ and $G^0$ have at most two generators;

(3) if $m=1$ , then $\mathcal {G}$ and $G^0$ are solvable.
The point $\hat {\mathbf {P}}=(0,0,e^{bt})$ is a nonequilibrium solution to (3.1). Hence, the variational equations along $\Gamma $ are
such that $\mathbf {P}=(\nu ,\eta ,\vartheta +e^{bt})$ is in (3.1). Trivially, this may be reduced to the closed subsystem
Applying the change of variable $\tau =e^{bt}$ , the above equation is expressed with rational coefficients:
Hence, letting $\nu (\tau )=\chi (\tau )\tau ^{(1\sigma (\mathcal {F}'(0)1)b)/2b}$ ,
Using the convention of Kovačič [Reference Kovačič9], since the coefficient of $\chi $ is a rational function, when we refer to the poles of the coefficient, we mean the poles in $\mathbb {C}$ . If $r=z_1/z_2$ with $z_1,z_2\in \mathbb {C}[\tau ]$ relatively prime, then the poles of r are the zeros of $z_2$ and the order of the pole is the multiplicity of the zero of $z_2$ . By the order of the coefficient at $\infty $ , we mean the order of $\infty $ as a zero of the coefficient so that the order of the coefficient of $\chi $ at $\infty $ is $\deg z_2\deg z_1$ .
Lemma 3.5. The differential Galois group G of (3.4) is infinite.
Proof. The coefficient of $\chi $ has two poles at $0$ and $\infty $ with orders $2$ and $1$ . According to the theorem provided in Section 2.1 of [Reference Kovačič9], the conditions for cases 1 and 3 do not hold and hence, G is not finite.
Lemma 3.6. The identity element $G^0$ of G is solvable if and only if
Proof. Letting $s=2\sqrt {\sigma \tau }/b$ and $\hat {\chi }=\chi /s$ , (3.4) becomes the Bessel equation
such that
Since $\hat {\chi }'+s^{1}\hat {\chi }=0$ has the nontrivial solution $s^{1}\in \mathbb {C}$ , we have $G\subset \text {SL}(2,\mathbb {C})$ . Referring to [Reference Kolchin8], when $2\rho \in \mathbb {Z}_{\text {odd}}$ , we may (replacing $\rho$ by $\rho$ if necessary) suppose that $\mu=\rho\frac{1}{2}\geq 0$ so that
is a fundamental system of solutions to (3.5), such that each solution is exponential over $\mathbb {C}$ and their product is in $\mathbb {C}$ .
Hence, the differential Galois group G is a diagonal group
which has a solvable (abelian) identity element $G^0$ if $2\rho \in \mathbb {Z}_{\text {odd}}$ . Otherwise, when $2\rho \notin \mathbb {Z}_{\text {odd}}$ , [Reference Kolchin8] demonstrates that $G=\text {SL}(2,\mathbb {C})$ . Kovačič’s algorithm may also be applied to derive the Galois group of (3.5) (see [Reference Duval and LodayRichaud5, Reference Kovačič9]).
Theorem 3.7. If (3.1) is Bintegrable in the meromorphic category, then
such that $\sigma ,b\neq 0$ and $n\in \mathbb {Z}$ .
Proof. Since the differential Galois group of (3.3) is a normal subgroup of the differential Galois group of (3.4), the differential Galois group of (3.3) is also infinite by Lemma 3.5. The identity element $G^0$ is a normal subgroup of G with finite index. Hence, $G^0$ is a trivial subgroup if and only if G is finite [Reference MoralesRuiz13], which implies that $G^0\neq \{\mathbf {1}\}$ and (3.1) is not completely integrable with meromorphic first integrals by Theorem 3.4.
According to Lemma 3.6, if
the identity elements of the differential Galois groups of (3.5) and (3.3) are not solvable. Therefore, the identity element $G^0$ of (3.3) is not abelian. After rearranging (3.7) and considering the contrapositive, the criterion in (3.6) follows trivially in view of Theorems 3.3 and 3.4.
Acknowledgements
I am grateful to Rahil Valani for his kind discussions and invaluable suggestions. The anonymous referee’s suggestions to improve the presentation of the paper are also very much appreciated.