# Drifting Faraday patterns under localised driving

Juan Marín, Rafael Riveros-Ávila, Saliya Coulibaly, Majid Taki, Leonardo

Gordillo, Mónica García-Ñustes

► To cite this version:

Juan Marín, Rafael Riveros-Ávila, Saliya Coulibaly, Majid Taki, Leonardo Gordillo, et al. Drifting Faraday patterns under localised driving. Communications Physics, 2023, 6 (1), pp.63. 10.1038/s42005-023-01170-8. https://hal-04505240

HAL Id: hal-04505240

https://hal.science/hal-04505240v1

Submitted on 16 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Drifting Faraday patterns under localised driving

Juan F. Marín ID <fnref n="1" />✉, Rafael Riveros-Ávila ID ², Saliya Coulibaly ID ³, Majid Taki³, Leonardo Gordillo ID ⁴✉ &
Mónica A. García-Ñustes ID ²

Physical systems undergoing spontaneous pattern formation are governed by intrinsic length scales that may compete with extrinsic ones, resulting in exceptional spatiotemporal behaviour. In this work, we report experimental and theoretical evidence that spatial non-uniformity sets Faraday-wave patterns in motion, which are noticeable in the zigzag and drift dynamics exhibited by their wave crests. We provide a minimal theoretical model that succeeds in characterising the growth of localised patterns under nonuniform parametric driving. Furthermore, the derivation accounts for symmetry-breaking nonlinear gradients that we show are the source of the drift mechanism, which comes into play right after the system has crossed a secondary bifurcation point. Numerical solutions of the governing equations match our experimental findings and theoretical predictions. Our results advance the understanding of pattern behaviour induced by nonuniformity in generic nonlinear extended systems far from equilibrium.

<footnote n="1">¹Departamento de Física, Facultad de Ciencias Naturales, Matemática y del Medio Ambiente, Universidad Tecnológica Metropolitana, Las Palmeras 3360, Nuñoa 780-0003, Santiago, Chile. ²Instituto de Física, Pontificia Universidad Católica de Valparaíso, Casilla 4059, Chile. ³Université de Lille, CNRS, UMR 8523, Laboratoire de Physique des Lasers Atomes et Molécules, F-59000 Lille, France. ⁴Departamento de Física, Universidad de Santiago de Chile, Av. Victor Jara 3493, Estación Central, Santiago, Chile. ✉email: j.marinm@utem.cl; leonardo.gordillo@usach.cl</footnote>
Highly ordered replication of individual structures in patterns is generally related to the natural selection of an intrinsic length scale. This length is deeply rooted in the properties of the physical system supporting them, which can drastically differ. Patterns are frequently observed in nonequilibrium systems¹ and in very different contexts from the inanimate world, such as fluids, chemical reactions, magnetic textures, nonlinear optics, and granular layers²⁻⁵. Often, most of the geometric and dynamic features of the patterns can be described by finding the link between the intrinsic length scale and the physical quantities in play. However, there is little understanding of how patterns may be affected if a second extrinsic length scale comes into play. This is a relevant question considering that systems are often driven into out-of-equilibrium states in limited regions of space, and the driving extent may compare to the intrinsic length scale of the system.

Standing gravity waves emerging on the free surface of a vertically vibrated fluid are one of the simplest examples of patterns in out-of-equilibrium systems. The vertical oscillation of the container has the effect of modulating the effective gravity. When the amplitude of such oscillations is above a certain threshold, one observes the spontaneous formation of a regular pattern of standing waves with a well-defined wavelength oscillating at half the driving frequency (parametric instability). These waves are known as Faraday patterns and were first reported in Faraday's seminal experimental work on vertically vibrated fluid and granular layers⁶,⁷. When the container is uniformly vibrated as a whole, Faraday patterns and localised structures support a rich assortment of spatiotemporal phenomena, such as pattern selection⁸, quasi-patterns⁹, alternating patterns¹⁰, supersquare patterns¹¹, oscillons¹², underlying streaming flows¹³, wavelength changing instabilities (Eckhaus instability)¹⁴,¹⁵, turbulence generation¹⁶,¹⁷, and complex dynamics of localised structures¹⁸⁻²⁰. However, the influence of nonuniform vertical vibration on the pattern dynamics²¹ remains a relevant open question that is puzzling due to the lack of suitable models.

In this work, we report the existence of drift instabilities in localised Faraday patterns induced at the surface of a fluid under a spatially nonuniform parametric drive. The counterintuitive effect of inducing a propagating phase in the system, even when the vertical vibration modulation is fixed in space, is an intriguing attribute and differs from other experimental realisations featuring special types of boundary conditions, such as annular channels²²⁻²⁴. We show that such drift is entirely induced by the nonuniform nature of the vertical drive at the bottom of the fluid container and thus cannot be observed under uniform driving. We use the normal form theory to explain the observed drift through an amplitude equation for Faraday patterns under localised driving. We demonstrate that the evolution of the instabilities at the first order of nonlinearity is described by a quintic Complex Ginzburg-Landau equation with Weber-like and self-phase modulation terms in the form of nonlinear gradients. Such nonlinear gradients trigger drift instabilities through a spontaneous nonlinear symmetry breaking above a secondary bifurcation threshold.

## Results

**Drift instabilities in a vertically vibrated fluid.** The experimental setup in Fig. 1a was suitably devised to study the spatiotemporal interactions of patterns under localised driving. Here, Faraday patterns, i.e., subharmonic standing waves induced by vertical vibrations, are sustained on the surface of a quasi-one-dimensional fluid of depth $h = 20$ mm in a container with a soft elastic bed. The free surface of the fluid is at $y = 0$ when at rest, and displacements of the surface from the equilibrium condition

are described by the time-space dependent field $\eta(x, t)$ as shown in the Fig. 1d, e. The novelty of the setup is the way in which the system is forced. Instead of uniform vibrations, an array of linearly spaced pistons mechanically actuates on the soft bed (at $y = -h$), thus introducing localised vertical vibrations into the system. As detailed in the Methods, pistons are driven by a brushless motor through an array of rotating cams, which control the frequency and amplitude of oscillations of the set of pistons. For the bottom motion, we programmed a controlled oscillating Gaussian profile, $\exp(-x^2/2\sigma_t^2)$, with a typical extrinsic length of injection $L \sim 20$ cm. The wavelength of the standing wave is related to $L$ in a nontrivial way, but it is close to the uniform-forcing case, $\lambda_n = L/n\pi$, where $n$ is the number of visible anti-nodes on the stationary surface wave. The time-dependent part of the oscillatory drive is controlled in the laboratory by

$$ \Gamma(t) = \Gamma_0 \cos(2\pi ft), \quad (1) $$

where $f$ is the forcing frequency. The acceleration amplitude of the forcing $\Gamma_0$ can also be written as a function of the maximum bottom displacement $a_0$ and $f$.

To characterise the spatiotemporal patterns in our experiment, we recorded the fluid surface with a high-velocity camera (400 frames per second), as shown in Fig. 1b. The spatiotemporal data of the water level were obtained after several image-processing steps, including brightness and gamma corrections, resulting in a clear view of the liquid interface on the glass wall, as shown in Fig. 1c and d. Afterwards, a Gaussian filter with kernel size $3 \times 3$ pixels was applied. A Canny edge detector was used to find continuous edges in the image. For further analysis, the final data were stored as a discrete matrix $\eta_{ij}$ representing the continuous water level $\eta(x, t)$. A typical spatiotemporal diagram is shown in Fig. 1e. The temporal envelope of the pattern oscillations can be extracted by applying a Hilbert transform over time to the spatiotemporal data. Figure 1e shows a locally driven Faraday pattern for $f = 14.8$ Hz and $\Gamma_0 = 0.613$, and Fig. 1f is a spatiotemporal representation of the modulus of the Hilbert transform. The real part of the Hilbert transform relates to a stroboscopic view of the phenomenon at the frequency of the wave oscillation $f/2$.

At a fixed forcing frequency $f = 14.8$ Hz, a locally driven standing Faraday pattern with mode $n = 6$ forms above the threshold $\Gamma_0^f = 0.331$, as shown in Fig. 1e and f. By increasing the drive amplitude above $\Gamma_0^f = 0.613$, the pattern develops a zigzag motion, noticeable in the oscillations exhibited by the wave crests in the spatiotemporal diagrams of Fig. 2a. By further increasing the driving amplitude above $\Gamma_0 = 0.691$, we discovered drifting patterns and nonzero mean velocity in the system, as depicted in Fig. 2a. The drift superimposes with the zigzag motion and features the continuous annihilation and creation of modes at the opposite sides of the pattern. Such a drift of localised Faraday patterns is a nontrivial effect that has not been reported in the literature.

Both the drift velocity and zigzag frequency display a strong dependence on the drive amplitude. To characterise the dependence, we analysed the mean trajectory of the crests to compute the drift mean velocity as a function of $\Gamma_0$. Figure 2b shows the experimental measurements of the drift average velocity $\langle v \rangle$ for different values of $\Gamma_0$. Additionally, the frequencies of the zigzag motion were calculated and fitted to a power law, as shown in the inset of Fig. 2b. The curve fits the experimental observations for $\Gamma_0 \gtrsim 0.676$ well, especially when we consider the nonnegligible effect of unavoidable small additive noise, which is typical in experiments. Additive noise in bifurcations of emerging patterns²⁵ affects the most likely value
<figure><img src="image_3.png" /><figcaption>Fig. 1 Experimental generation of localised Faraday patterns. a Experimental setup for generating Faraday patterns in the surface of fluid under localised vertical vibrations. b Image received from the high-velocity camera. c Processed image with gamma manipulation, brightness manipulation and Gaussian filter. d Overlap of the processed image and detected free water surface η(x) at an arbitrary time t*. e Spatiotemporal diagram of the free surface corresponding to a Faraday pattern η(x, t). f Modulus of η̃(x, t) corresponding to the time envelope of the free surface oscillations.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>Fig. 2 Experiments of zigzag and drifting patterns. a Typical behaviour of the drifting patterns in experiments with different drive amplitudes. b Measurements of the mean velocity 〈v〉 of the drifting patterns as a function of the drive amplitude Γ₀. The fit using Eq. (8) (without including noise effects) is shown as dashed lines, while the fit including noise is shown as a solid line. The inlet figure shows a power-law fit of the crest oscillation frequency as a function of the drive amplitude Γ₀. Error bars represent the mean of the standard deviation of measurements performed for the correspondent amplitude Γ₀ at the same frequency.</figcaption></figure>

of the mean drift velocity as follows:

$$ \langle v \rangle = A \sqrt{\frac{1}{2} \left( \Delta_D + \sqrt{\Delta_D^2 + 2\xi} \right)}, \quad (2) $$

where Δ_D := Γ₀ − Γ_c, Γ_c is the critical drive amplitude for the emergence of the drift velocity, and ξ is the noise intensity. Figure 2b shows the fitting curve using Eq. (2) in a solid red line, where the fitted parameter values are A = 2.463, Γ_c = 0.683, and ξ = 0.01. We observed agreement between the experimental measurements and the fit using Eq. (2).

Figure 3 summarises our experimental results. In Fig. 3a, Γ₀ < Γ₀ < Γ_c, and the pattern is absolutely steady. Note that although the Gaussian driving profile in the bottom is centred around the origin, the maximum of the pattern is slightly shifted towards the negative values of x. Figure 3b, Γ₀ > Γ_c, illustrates a combination of a unidirectional drift with zigzag motion. A further increase in Γ₀ beyond the threshold of drift instability creates a complex convective state, as shown in Fig. 3c. In this latter case, the pattern dynamics exhibit convection towards both possible directions. We also noticed the emergence of branching phenomena, where a drifting maximum of the pattern breaks up into two crests momentarily travelling in opposite directions. We

also observed intermittence in this regime: drift and branching ceased to occur for a short period of time, as evidenced in Fig. 3c for t ≳ 5 × 10⁵T. As explained below, the interplay of nonlinearities and the localised vertical drive plays a fundamental role in the underlying bifurcations behind these phenomena.

**Modelling drift in localised Faraday patterns.** It has been shown that pattern formation in the surface of a weakly viscous fluid in a trough with a large aspect ratio can be modelled by the parametrically driven nonlinear Schrödinger (pdlS) equation³,

$$ \partial_t \psi = -(\mu + i\nu)\psi - i\alpha\partial_{xx}\psi - i|\psi|^2\psi + \gamma(x)\psi^*, \quad (3) $$

written here in dimensionless form. The pdlS model is simple yet can provide profound insights into the generic phenomenon of pattern formation in fluids and many extended nonequilibrium systems²⁶⁻³¹. In Eq. (3), pattern dynamics are described through a complex order parameter ψ that depends on space and time and is related to the displacement η of the free surface by η(x, t) = Reψ(x, t) exp(iπft). The parameter μ is the effective dissipation of the system, which accounts for both the fluid bulk viscosity and wave damping due to the friction with the walls³². The parameter α is related to a spatial coupling length, and the
<figure><img src="image_5.png" /><figcaption>Fig. 3 Comparison between experiments and simulations of drifting patterns. a-c Experimental evidence of the drift instability in a fluid subject to vertical and localised drive. The amplitude of the vertical drive is (a) $\Gamma_0 = 0.56$, (b) $\Gamma_0 = 0.70$, and (c) $\Gamma_0 = 0.84$. The detected level of the surface of the fluid is indicated by the colour scale. d-f Drift instability in the numerical simulations of the pdnIS equation (3) with $\mu = 0.45$, $\alpha = 1$, $\nu = 1.0$, and $\sigma_i = 16$ for (d) $\gamma_0 = 0.55$, (e) $\gamma_0 = 0.8375$, and (f) $\gamma_0 = 0.9000$.</figcaption></figure>

parameter $\nu = [(f/f_1)^2 - 1]/2$ is a small frequency detuning between the driving frequency and the natural frequency of the first transverse mode, $f_1 = (2\pi)^{-1}\sqrt{g\kappa\tanh\kappa}$, with $\kappa = 2\pi/\lambda$ being the wavenumber of the fundamental mode$^{33}$.

Under uniform driving, the system destabilises when$^3 \gamma > \gamma_c^h = \mu$ through the spontaneous growth of a pattern for positive values of $\nu$ via a supercritical bifurcation$^{31,34}$ with leading eigenmodes of wavenumber $\kappa_c = \pm\sqrt{\nu/\alpha}$. The effect of a localised profile in the vertical drive has been analysed$^{21}$ by introducing the space-dependent Gaussian profile $\gamma(x) = \gamma_0 \exp(-x^2/2\sigma_i^2)$, where $\sigma_i$ gives the extension of the injection zone and $\gamma_0 = \Gamma_0/4g$ is the dimensionless amplitude of the vertical drive. The linear stability analysis under a slightly localised Gaussian profile, i.e., $\sigma_i \ll \sqrt{\alpha}$, shows that the envelope $C_0$ of the pattern has a discrete response given by the Gauss-Hermite polynomials$^{21}$. Urra et al.$^{21}$ showed that the threshold of the instability of the m-th Gauss-Hermite mode is given by

$$ \gamma_0^{(m)} = \mu + (2m + 1) \frac{1}{\sigma_i} \sqrt{\frac{\nu}{\alpha}} \quad (4) $$

and therefore, localised Faraday patterns appear at $\gamma_0 > \gamma_0^{(0)}$ for $m=0$.

In the lower panel of Fig. 3, we show the results from the numerical simulation of the pdnIS equation (3) for the given values of the parameters. Near the fundamental threshold of instability, the Faraday patterns are rather regular and highly ordered, as shown in Fig. 3d. The amplitude of the pattern increases with $\gamma_0$. However, at some value $\gamma_D \in (\gamma_0^{(2)}, \gamma_0^{(3)})$, a secondary instability emerges, and the pattern dynamics become complex for both the amplitude and phase. In Fig. 3d, we fixed $\gamma_0$ slightly below the secondary instability onset $\gamma_D$ so the localised pattern is absolutely stable. Note that the maximum of the wave amplitude is shifted towards the left of the origin, as we observe in Fig. 3a. If $\gamma_0$ is slightly above the threshold $\gamma_D$, the unidirectional zigzag drift becomes evident, as shown in Fig. 3e. Experiments and simulations suggest the existence of a nonlinear symmetry-breaking mechanism in the system. Further increasing the drive

amplitude $\gamma_0$, the system evolves towards a more complex state with a superposition of oscillations (suggesting an underlying Hopf instability), branching and drifting in both directions, as shown in Fig. 3f.

To provide a theoretical explanation of the observed behaviour, we characterised the dynamics of localised Faraday patterns near the fundamental threshold of instability through Eq. (3) using the method of normal forms (see Methods). Under this widely used formalism, a nonlinear system takes the “simplest”—or so-called normal—form preserving its essential features, i.e., dynamic behaviour, near a bifurcation point$^{35}$. To capture the dynamics close to the spatial instability, we introduced a bifurcation parameter $\delta = \gamma_0 - \gamma_0^{(0)}$ and a small dimensionless parameter $\epsilon = \sqrt{\alpha}/\sigma_i \ll 1$ characterising the length of the driving zone. With the condition $\epsilon \ll 1$ we assumed that the injection zone is large compared to the spatial coupling length, i.e., the system behaves as a weakly nonuniform system. We studied the behaviour of solutions around a Faraday pattern with critical wavenumber $\kappa_c$ introducing a modulating complex amplitude function $C(x, t)$—an order parameter—, i.e.

$$ \psi(x, t) = C(x, t)e^{ik_c x} + c.c. + h.o.h., \quad (5) $$

where h.o.h denotes the higher-order harmonics. Below the bifurcation point, $\gamma_0 < \gamma_0^{(0)}$, one observes a stable uniform solution $C=0$ (disordered phase). Such a uniform solution becomes unstable just above the bifurcation point, and $C(x, t)$ has its own spatiotemporal dynamics modulating the underlying Faraday pattern (self-organised phase). At the leading order of non-linearity, one obtains the following governing equation for the amplitude:

$$ \partial_t C = \frac{2\nu}{\alpha\mu} \partial_x C + \left( \delta - \frac{\mu}{2\sigma_i^2} x^2 \right) C - i \frac{2}{\mu} \sqrt{\frac{\nu}{\alpha}} [3\partial_x(|C|^2C) - 2C\partial_x(|C|^2)] - \frac{9}{2\mu} |C|^4 C. \quad (6) $$

For nonextended systems, the homogeneous limit of Eq. (6) turns into the well-known normal form derived by Coullet et al.$^{34}$. The first two terms on the right-hand side of Eq. (6) are Weber-like terms that give the known linear limit of the nonuniform system$^{21}$. At the modulational instability, the
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Comparison between theory and numerical simulations. **a, b** Normalised amplitudes $|\tilde{\psi}|$ and $|C|$ obtained numerically from Eq. (3) and Eq. (6), respectively, with $\mu = 0.45$, $\nu = 1$, $\alpha = 1$, $\sigma_1 = 16$, and (a)$\gamma_0 = 0.67$, (b)$\gamma_0 = 0.8$. **c** Normalised maximum amplitude of the saturated Faraday pattern, $|C_{\text{max}}|$, as a function of $\gamma_0$, according to Eq. (7) (solid line) and the numerical solutions of pdnIS Eq. (3) (stars). Drift is observed above a critical value $\gamma_D \approx 0.837$ indicated with a vertical dashed line.</figcaption></figure>

amplitude of the unstable Gauss-Hermite modes is saturated by the quintic nonlinearity given by the last term in Eq. (6).

Surprisingly, Eq. (6) displays unexpected terms proportional to $\partial_x(|C|^2C)$ and $C\partial_x(|C|^2)$. These are gradient-driven terms featuring self-phase modulation (SPM) effects in the slowly evolving modes of the localised Faraday patterns. Such terms are often introduced in the nonlinear Schrödinger equation to account for self-steepening effects and self-frequency shifting via stimulated Raman scattering$^{36-40}$. To our knowledge, this is the first time a complex quintic Ginzburg-Landau equation with Weber-like terms has been derived with these SPM terms appearing naturally from the weakly nonlinear analysis. We show below that such terms explain both the absolute and drift instabilities observed in our experiments of localised Faraday patterns.

In Fig. 4a, b, we compare the predictions of our normal form (6) and the pdnIS equation (3). In both cases, we numerically solved Eqs. (3) and (6) with no-flux boundary conditions. The initial conditions were given by small-amplitude random distributions of values in space for the real and imaginary parts of $\psi$ and $C$. There is agreement between both models near the fundamental threshold of instability. Figure 4a and b confirm that the extension $\sigma$ and the envelope profile of the Faraday pattern are well described by the normal form. Beyond the fundamental threshold, as one would expect, predictions from the normal form show some mismatch from the outcome of the pdnIS equation. Figure 4c shows the maximum amplitude of the Faraday pattern as a function of $\gamma_0$, and demonstrates that the normal form (6) describes the nature of the bifurcation of the pattern well, as explained in detail below.

The SPM terms in Eq. (6) are induced by nonlinear gradients, which reveal underlying complex phase dynamics. The gradient terms break the transverse reflection symmetry ($x \to -x$) and dramatically affect the pattern selection in the system$^{41,42}$. Thus,

the nonlinear gradients in Eq. (6) yield a nonlinear symmetry breaking where convection is characterised by an amplitude-dependent group velocity. In Fig. 5, we demonstrate examples of drift instabilities in the system obtained from numerical simulation of Eq. (6) for some given parameter values.

In Fig. 5a, we show data for a run in which the system exhibited absolute stability with a pattern amplitude dominating over drifting and zigzag-like motion. In the absence of SPM terms, the localised Faraday pattern spontaneously grows centred around the Gaussian driving with width$^{21} \sigma_w = (\sqrt{\nu/\alpha} \sigma_i / \mu)^{1/2}$. In Fig. 5a, we indicate with dashed lines the $\sigma_w$ centred at the Gaussian driving for reference, evidencing the clear shift in Re(C) due to the SPM terms. The first frames of the simulation in Fig. 5b show that there is a transient drift towards the left as the pattern grows. However, the convection is too weak to overcome saturation, and the system rapidly becomes absolutely stable. Although the imaginary part of $C$ displays a similar shift, the absolute value $|C|$ exhibits no shift, and the phase lacks any spatiotemporal dynamics in this case, as shown in Fig. 5c.

Increasing the value of $\gamma_0$, we found for some value $\gamma_D \in (\gamma_0^{(2)}, \gamma_0^{(3)})$ that the propagation overcomes amplification and the system exhibits a nonlinear transition to drift instability. Figure 5d shows this scenario, in which a Faraday pattern displays a drift in the real part of its envelope. The real and imaginary parts of $C$ exhibit the same drift, except for a phase factor. However, the absolute value $|C|$ has no drift, as we show in Fig. 5e. Thus, the pattern drift is entirely due to the phase dynamics, as evidenced in Fig. 5f.

A further increase in $\gamma_0$ results in acceleration of the drift, as shown in Fig. 5g. The drifting wave also undergoes changes in both the frequency and wavelength. These observations confirm the nonlinear nature of the drift instability. Note that the drifting pattern becomes nonpropagating near the tails of the Gaussian
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Drift instabilities in the normal form. Numerical results obtained from Eq. (6) with $\mu = 0.45$, $\nu = 1$, $\alpha = 1$, and $\sigma_t = 16$. **a–c** A wave exhibiting absolute stability for $\gamma_0 = 0.80$. **d–f** Drift instability for $\gamma_0 = 1.2$. **g–i** A strongly nonlinear drift instability for $\gamma_0 = 1.62$.</figcaption></figure>

drive, as highlighted in the inset of Fig. 5h. This pinning effect stems from the spatial nonuniformity of the vertical drive, which imposes smooth amplitude variations at the tails of the Gaussian drive and induces a boundary-like condition. When spatial variations in the pattern envelope become comparable to modulation amplitude, they could result in coupling between the slow scale of the envelope to the fast scale of the modulation of the underlying Faraday pattern, leading to pinning of the drifting wave$^{43}$. Moreover, the phase dynamics exhibit travelling waves inside the injection zone defined by the Gaussian drive, as evidenced in Fig. 5f and i. The mean velocity of such travelling waves decreases as $\gamma_0$ decreases.

**Drift velocity characterisation.** To gain further insight into the dynamic behaviour of the drifting localised patterns, we looked for an explicit expression for the average pattern velocity, $\langle v \rangle$. The key idea is that at the tails of the Gaussian driving, the normal form of Eq. (6) has an extra nonresonant term arising from the coupling between the slow-scale dynamics of $C(x, t)$ and

the fast scale of the modulation of the Faraday pattern. This effect has been observed in the pinning of patterns by a boundary-like condition$^{43}$. By calculating the normal form assuming that scale separation is no longer valid at the tails of the localised drive, we obtain the full amplitude equation,

$$ \partial_t C = \frac{2\nu}{\alpha\mu}\partial_{xx}C + \left(\delta - \frac{\mu}{2\sigma_t^2}x^2\right)C - i\frac{2}{\mu}\sqrt{\frac{\nu}{\alpha}}[3\partial_x(|C|^2C) - 2C\partial_x(|C|^2)] \\ - \frac{9}{2\mu}|C|^4C - \frac{3i}{2\sqrt{\nu}}e^{-i2\kappa_x x}C^2\partial_x C, \qquad (7) $$

where the last term proportional to $\exp(-i2\kappa_c x)$ is nonresonant and has the form of a nonlinear gradient. Further calculations lead to a formula for the average velocity (see Methods)

$$ \langle v \rangle = \sqrt{h-h_c}, \qquad (8) $$

where $h$ is a drift bifurcation parameter and $h_c$ is a critical value. Thus, we conclude that for $h$ close to $h_c$, the system exhibits a
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Zigzag drifting and branching phenomena. Phenomena captured by the normal form with nonresonant terms. **a–c** Numerical simulations of Eq. (7) for $\nu = 0.25$, $\mu = 0.15$, and (a)$\gamma_0 = 0.3$, (b)$\gamma_0 = 0.37$, and (c)$\gamma_0 = 0.38$. **d–f** Real part of the complex field $\psi$ reconstructed from (a), (b) and (c) through Eq. (5).</figcaption></figure>

saddle-node bifurcation to drift instability with the pattern velocity increasing with the square root of $h$. In Fig. 2b, dashed lines indicate the fit of the velocity experimental measurements using Eq. (8). Note that neglecting noise effects on the bifurcation, we have $\xi \to 0$ in Eq. (2), and we recover our theoretical prediction on Eq. (8), which is valid for the deterministic case.

Figure 6 shows the results from numerical simulations of the full normal form (7) for the given parameter values. Figure 6 shows that our full normal form captures the zigzag drifting and branching phenomena observed in the experiments. Figure 6a–c shows the real part of the order parameter C. The pattern in Fig. 6a exhibits a unidirectional drift with a zigzag-like motion for $\gamma_0 = 0.3$. Increasing the dimensionless drive amplitude to $\gamma_0 = 0.37$, we observe that the frequency of the zigzag-like motion increases, as shown in Fig. 6b, which is in agreement with our experimental observations in Fig. 2. The drift velocity also increases from Fig. 6a to Fig. 6b, as expected. Further increasing the value of $\gamma_0$, we observe the branching phenomenon shown in Fig. 6c. Figure 6d–f shows the corresponding dynamics of the real part of the complex field $\psi$, reconstructed from the numerical solutions of Eq. (7) through Eq. (5). The similarity between Fig. 6d, e and f with Fig. 3a–c, respectively, is noticeable.

**Self-phase modulation and nonlinear saturation.** Here, we address the question of how drift instabilities impact our current understanding of the nonlinear saturation mechanism of localised patterns. This problem was first assessed in ref. 21 but without taking into account the SPM terms. Following previous calculations, we performed a multiple-scale analysis on Eq. (6) to investigate the role of the SPM terms. We found that SPM has no effect on the saturation of the fundamental mode close to the instability threshold. Assuming that the spatial and temporal scales of perturbations on the pattern amplitude are large, one obtains that the evolution of the amplitude $D$ of the fundamental mode for $\gamma_0 \sim \gamma_0^{(0)}$ is governed by

$$ \frac{d}{dt}D = \delta D - \frac{9}{2\sqrt{3}\mu}D^5, \qquad (9) $$

in full agreement with previous results derived without the SPM terms. Equation (9) predicts the scaling law $D \propto \delta^{1/4}$, which has

been previously verified in experiments$^{21}$ and gives the signature of a supercritical bifurcation. Figure 4c shows the maximum stationary amplitude $|C_{max}|$ of the Faraday wave as a function of $\gamma_0$: solid lines are predicted from Eq. (9), and symbols depict the results from the numerical simulations of the pdnls Eq. (3). We conclude that the results of the multiple-scale analysis on the normal form (6) are in good agreement with both the numerical simulations of the pdnls equation (3) and the Faraday-wave experiments.

## Discussion

In ref. 21, we derived an amplitude equation describing the growth of nonlinear instabilities around the fundamental threshold $\gamma_0^{(0)}$, performing an appropriate expansion on the dispersion relation obtained from the WKBJ analysis. The resulting amplitude equation provides insight into how the nonlinear instability of the Gauss-Hermite modes gives rise to a saturated pattern. However, such an amplitude equation was not formally derived from the weakly nonlinear analysis.

Our findings unveil a connection between the drift instability of coherent structures and heterogeneity in general nonlinear systems far from equilibrium. In our system, such a connection is provided by the interplay between the wavelength of the Faraday pattern (an intrinsic length scale which depends mainly on the intrinsic and drive parameters of the system) and the spatial width $\sigma_i$ of the Gaussian drive (an extrinsic length scale which is independent on the intrinsic parameters of the system and tunable). Indeed, the drift instabilities reported in this work are spontaneously induced by the spatial nonuniformity of the drive, i.e. a spatially featured non-drifting drive can generate drifting patterns. Under uniform drive, the amplitude of the Faraday pattern will be space-independent, and hence, the SPM terms causing the drift in Eq. (6) vanish. Note that only the nonuniform Weber-like term in Eq. (6) vanishes in the limit $\sigma_i \to \infty$ (i.e., $\gamma(x) \to \gamma_0$). However, although the SPM terms are not directly coupled to the Gaussian driving width $\sigma_i$, they are indirectly linked through the nonuniform profile of $C(x, t)$ induced by the Gaussian drive. We believe this is the reason why such drift has not been reported in previous studies, which have been focused on uniformly accelerated systems. Recently, the role of spatial nonuniformity has often been identified as an important factor in
nonlinear systems<sup>21,44–47</sup>. Our theoretical insights on the underlying instabilities of localised patterns may be relevant in understanding other wandering localised structures observed in one and two-dimensional heterogeneous nonequilibrium systems. Moreover, drift instabilities in two dimensions open up possibilities for coherent structures arising from convection and multi-directed stream. It also remains the open question of how a spatial distribution of driving forces could interact with other self-driven drifting phenomena, such as those observed in walking droplets<sup>48</sup>. Our findings open avenues of research, validating the use of SPM terms such as those naturally emerging in Eq. (6) to capture drift instabilities induced by nonuniform drives in different physical, chemical, and biological systems.

To summarise, we have shown that nonuniformity can trigger drift instabilities in the evolving amplitude of localised Faraday patterns. We performed an experimental study of Faraday waves on the free surface of a fluid under localised vertical driving. We show that above a secondary threshold of instability, the localised and nondrifting drive triggers drift instabilities in the surface patterns. We provide a theoretical explanation of the phenomenon based on a nonuniform version of the pdnLS equation as a simple model for parametric driven patterns. Under the assumption of a vertical drive modulated by a spatially spread Gaussian profile featured by an extrinsic length scale, we derived an evolution equation governing the slow dynamics of the pattern amplitude. The resulting normal form includes two unreported SPM terms that are linked to complex dynamics in the underlying patterns. We show that such terms trigger drift via nonlinear symmetry breaking induced by a nonlinear gradient. Numerical simulations of the pdnLS equation are also in agreement with both the theory and experimental data. These results extend our general understanding of nontrivial dynamical effects of non-uniformity in nonlinear extended parametric oscillators and how pattern selection can be dramatically affected by the nonuniform nature of the input of energy.

## Methods

**Experimental setup.** We used a rectangular water basin, 15 mm long, 490 mm wide, and 100 mm deep, whose bottom centre consists of a rectangular soft elastomer made out of high-performance platinum cure liquid silicone Dragon Skin 10 medium 10A. The Gaussian drive, directly linked to $y(x) = y_0 \exp(-x^2/2\sigma_t^2)$, is mechanically achieved in this waterproof soft zone through bottom vibrations. To control both the driving amplitude and frequency, we used a brushless servomotor with feedback (Model No. BLM-N23-50-1000-B) controlled with a DMC-40x0 motion controller from Galil Motion Control. The servomotor was coupled to a shaft to transmit the rotational motion to an evenly spaced array of 13 pistons ($\Delta x = 16$ mm), which were placed underneath the soft zone. Each piston holds linear bearings to maintain vertical motion and has a tiny roller in contact with a rotary cam fixed to the shaft, making the transmission of rotational oscillatory motion to linear oscillatory motion possible as shown in Fig. 7a. The key of the mechanism is the special design of the cams (see Fig. 7b) according to the specific formulae in Fig. 7c. The pistons are pushed with springs towards the cams to avoid detachment under peak accelerations. The setup resembles the mechanical transmission system of a music box, as shown in Fig. 1a, and allows to deform the bottom of the channel with a programmed spatial distribution.

For the experiment, the basin was filled up to 20 mm with a mixture of water and white ink for image detection purposes. For data acquisition, we used a Phantom VEO 440L camera controlled by PCC 3.5 software. The resolution of the images was 1280 × 252 pixels, and the acquisition time was approximately 94 s at 400 frames per second for every amplitude. Images were saved in greyscale to simplify the contour detection process.

**Numerical methods.** Numerical simulations of the pdnLS equation (3) and our normal form (6) were performed with no-flux boundary conditions (rigid impermeable wall) using a 401-point spatial grid with resolution $dx = 0.25$. We used finite differences of the second-order accuracy for the spatial derivatives, and time integration was achieved using a fourth-order Runge-Kutta scheme with time-step $dt = 0.0001$. To obtain the amplitude of the drifting localised Faraday patterns shown in Fig. 3, we first numerically obtained the pattern profile $\psi(x, t)$ from the numerical simulation of Eq. (3). Then, the spatial complex amplitude $C(x, t)$ was obtained from the Hilbert transform in space,

<figure><img src="image_9.png" /><figcaption>**Fig. 7** The cam-roller-spring-piston mechanism. a Angular oscillatory motion displaces the piston to generate vertical oscillations. The special design of the cam (b) according to the formulae shown in the table (c), allows control of vertical displacement amplitude through angular amplitude in a range suitable for experiments (shaded region). The functions f and g smoothly close the cam edge, shown as the red curve in b.</figcaption></figure>

$\hat{\psi}(x, t) = (1/\pi) \int_{-L/2}^{L/2} d\chi \, \psi(\chi, t)/(x - \chi)$, which was computed at each time step of the numerical integration. When needed, the temporal complex amplitude is obtained via $\hat{\psi}(x, t) = (1/\pi) \int_{-L/2}^{L/2} d\tau \, \psi(\tau, t)/(x - \tau)$. Both were obtained numerically with Matlab®.

**Normal forms.** The theory of normal forms is a powerful method for the systematic construction of local, near-identity, and nonlinear transformations to simplify the equations describing the dynamics of complex nonlinear problems near a bifurcation point<sup>35</sup>. We started by writing Eq. (3) in the form $\partial_t \Psi = L\Psi + N$, where $\Psi := (\psi_R, \psi_I)$, $\psi := (\psi_R + i\psi_I, L)$ is a linear operator, and N is a nonlinear vector. We performed our transformations around a Faraday-pattern solution with critical wavenumber $\kappa_c$ according to

$$ \Psi = C(\eta, t) \begin{pmatrix} 1 \\ 0 \end{pmatrix} e^{i\kappa_x t} + \sum_{n=1}^{\infty} \mathbf{W}^{[n]} + \text{c.c.}, \quad (10) $$

where $C(\eta, t)$ is a complex amplitude that varies as a function of $\eta := \epsilon^{1/2}x$, as described in the Results, and the vector functions $\mathbf{W}^{[n]} \in \mathbb{C}^2$, $n \ge 1$ are higher-order corrections to be calculated. Assuming in Eq. (3) the scaling laws $C \sim \epsilon^{1/4}$, $\delta \sim \epsilon$, $\partial_t \sim \epsilon$ and $\partial_x \sim \epsilon^{1/2}$, one can expand the pdnLS equation at different orders of $\epsilon^{1/4}$ and obtain a hierarchy of linear problems of the form

$$ \mathcal{L} \mathbf{W}^{[n]} = \mathbf{b}_n, \quad (11) $$

which can be solved for $\mathbf{W}^{[n]}$ only if $\mathbf{b}_n$ is in the image of the linear operator $\mathcal{L}$. According to the Fredholm alternative<sup>2</sup>, at least one solution for $\mathbf{W}^{[n]}$ exists if $\exists!|\mathbf{v}| \in \ker \mathcal{L}^l$ such that $\langle \mathbf{v}|\mathbf{b}_n \rangle = 0$. Defining the inner product in the vector function space as

$$ (\mathbf{f}|\mathbf{g}) = \kappa_c \int_{x_o}^{x_o+1/\kappa_c} dx (\mathbf{f}^*)^T \mathbf{g}, \quad |\mathbf{f}|, |\mathbf{g}| \in \mathbb{C}^2, \quad (12) $$

the Fredholm alternative yields a solvability condition to calculate $\mathbf{W}^{[n]}$ at each order of $\epsilon^{1/4}$. Our normal form (6) is the solvability condition obtained at order $\epsilon^{5/4}$.

**Calculation of the drift velocity.** We wrote our full normal form Eq. (7) in polar coordinates by setting $C = \rho \exp(i\theta)$. After removing the phase factor, the imaginary part reads

$$ \partial_t \theta = 4\partial_{xx}\theta + \frac{8}{\rho}\partial_{x\rho}\rho\partial_{x}\theta - 4\rho^2(3\partial_{x}\rho + 2|\partial_{x}\rho|) - \frac{3\mu}{\beta}\rho(\rho\sin[2(\beta X - \theta)]\partial_{x}\theta + \cos[2(\beta X - \theta)]\partial_{x}\rho), \quad (13) $$

with $\beta = \nu/\sqrt{\alpha}$, $X = x\sqrt{\nu/\alpha}$ and $\tau = 2\mu t$. As depicted in Fig. 5f and i, the phase $\theta(x, t) = \arg(C)$ consists of two different superimposed dynamic behaviours: (i) a
travelling wave with some phase velocity $c = \omega/k$ and (ii) a monotonically increasing one modulo $2\pi$, such that $\arg(C) \in [-\pi, \pi]$ (following the conventional restriction of the inverse trigonometric functions to the principal branch). Based on this observation, we propose the ansatz
$$ \theta(X, \tau) = [\theta_c(X - c\tau) + \chi(\tau)] \mod (2\pi) - \pi, \text{ where } \theta_c \text{ and } \chi \text{ account for the tra-} $$
velling and increasing phase dynamics, respectively.

Close to the drift instability threshold, $d \langle\theta\rangle_X/dt$ slowly evolves over time, and therefore, $\chi(\tau) \simeq \langle\chi(\tau')\rangle_\tau$. Averaging Eq. (13) in both space and time, with $\langle\theta(X, \tau)\rangle_X = \sigma_w^{-1} \int_{-\sigma_w/2}^{\sigma_w/2} dX \theta(X, \tau)$ and $\langle\theta(X, \tau')\rangle_\tau = (\omega/T) \int_{\tau}^{\tau+T} d\tau' \theta(X, \tau')$, one obtains

$$ \frac{d\chi}{d\tau} = \frac{h}{2} - \frac{3\mu}{\beta} [A_0 \cos(2\chi) - B_0 \sin(2\chi)], \quad (14) $$

where $A_0 = \langle\rho(\rho-1)\partial_X\theta_c\sin[2(\beta X - \theta_c)]\rangle_{X,\tau}$, $B_0 = \langle\rho(\rho-1)\partial_X\theta_c\cos[2(\beta X - \theta_c)]\rangle_{X,\tau}$, and $h = -2(\partial_X\theta_c)_{X,\tau}$. Finally, analytically solving Eq. (14), we obtain Eq. (8) with $h_c = (6\mu/\beta)^2(A_0^2 + B_0^2)$.

**Method of multiple scales.** The method of multiple scales to study the effects of SPM in the Faraday pattern near the fundamental threshold of instability separates the time-scale dynamics of the different Gauss-Hermite modes controlling divergences in the perturbative developments⁴⁹. Just above the fundamental threshold of instability, the fundamental Gauss-Hermite mode, $\mathcal{G}_0$, dominates the amplitude dynamics of the Faraday pattern as higher-order modes decay²¹. As the bifurcation parameter is further increased, more Gauss-Hermite modes are excited. Each mode grows at its own growth rate or characteristic time. Assuming that the spatial and temporal scales of perturbations on the wave amplitude are large, the spatial derivative is rescaled by a small parameter $\zeta: \partial_x \rightarrow \zeta^{1/2}\partial_\sigma$, and the time scales slowly evolve as $T_i = \zeta^i t$ (with $i = 1, 2, ...$). Writing $C$ as a perturbation series of the different time scales,

$$ C(x, t) := \sum_{k=1}^{\infty} \zeta^{k/4} A_k(x, T_1, T_2, \dots), \quad (15) $$

and using a similar expansion in the bifurcation parameter,

$\delta = \delta_0 + \zeta\delta_1 + \zeta^2\delta_2 + \dots$, results in a hierarchy of problems at each order of $\zeta^{1/4}$, similar to Eq. (11) for $A_k$, $k \ge 1$.

At order $\zeta^{1/4}$, one obtains that $A_1$ is governed by the (time-independent) Weber equation²¹, whose general solution is given by a linear combination of the Gauss-Hermite modes with time-dependent coefficients $D_n(T_1, \dots), n \ge 0$. In particular, the saturation of the fundamental Gauss-Hermite mode will be given by the time evolution of the first of such time-dependent coefficients, $D_o(T_1, \dots)$.

At order $\zeta^{5/4}$, in the analogous hierarchical equation (11), one obtains a vector $\mathbf{b}_5$ that contains SPM terms. Given that the set of all Gauss-Hermite modes $\{\mathcal{G}_n\}_{n=0}^\infty$ is a basis of the space, we write the functions $A_2$ and $A_3$ as the linear expansion

$$ A_k = D_o^{(k)}(T_1)\mathcal{G}_0 + \sum_{n=1}^{\infty} D_n^{(k)}(T_1)\mathcal{G}_n, \quad k = 2, 3. \quad (16) $$

However, near the threshold of instability of the fundamental mode, the modes $\mathcal{G}_n$ with $n \ge 2$ are stable. Thus, the coefficients $D_n^{(k)}(T_1)$ for $n \ge 2$ decay over time. Once the Faraday pattern has completely evolved, the maximum amplitude of the pattern is $C_{\max} \simeq D_o := D$, and the functions $A_k$ ($k = 1, 2, 3$) are even in space. Thus, using the Fredholm alternative, from symmetry arguments and the inner product (12), one obtains that the contribution of the self-phase modulation terms vanishes, and the solvability conditions lead to Eq. (7).

## Data availability

The data recorded during the localised Faraday experiments are available from the corresponding author upon reasonable request. The remaining data that support the findings of this work are generated from computational simulations (see code availability), and are available from the corresponding author upon reasonable request.

## Code availability

The code for image processing and computing the fields $\eta(x, t)$ and $|\psi(x, t)|$ are available from our GitHub repository.

Received: 14 September 2022; Accepted: 10 March 2023; Published online: 07 April 2023

## References

1. Cross, M. & Greenside, H. *Pattern Formation and Dynamics in Nonequilibrium Systems* (Cambridge University Press, 2009).

2. Pismen, L. M. *Patterns and Interfaces in Dissipative Dynamics* (Springer Science & Business Media, 2006).

3. Miles, J. W. *Parametrically excited solitary waves*. *J. Fluid Mech.* **148**, 451–460 (1984).

4. Kenig, E., Malomed, B. A., Cross, M. C. & Lifshitz, R. Intrinsic localized modes in parametrically driven arrays of nonlinear resonators. *Phys. Rev. E* **80**, 046202 (2009).

5. Clerc, M. G., Coulibaly, S. & Laroze, D. Localized states and non-variational Ising-Bloch transition of a parametrically driven easy-plane ferromagnetic wire. *Physica D* **239**, 72–86 (2010).

6. Faraday, M. On the forms and states of fluids on vibrating elastic surfaces. *Philos. Trans. R. Soc. A* https://doi.org/10.1098/rstl.1831.0018 (1831).

7. Faraday, M. XVII. On a peculiar class of acoustical figures; and on certain forms assumed by groups of particles upon vibrating elastic surfaces. *Philos. Trans. R. Soc. A* https://doi.org/10.1098/rstl.1831.0018 (1831).

8. Douady, S. & Fauve, S. Pattern selection in Faraday instability. *EPL* **6**, 221–226 (1988).

9. Edwards, W. S. & Fauve, S. Patterns and quasi-patterns in the Faraday experiment. *J. Fluid Mech.* **278**, 123–148 (1994).

10. Péринet, N., Juric, D. & Tuckerman, L. S. Alternating hexagonal and striped patterns in Faraday surface waves. *Phys. Rev. Lett.* **109**, 164501 (2012).

11. Kahouadi, L. et al. Numerical simulation of supersquare patterns in Faraday waves. *J. Fluid Mech.* **772**, R2.1–R2.12 (2015).

12. Shats, M., Xia, H. & Punzmann, H. Parametrically excited water surface ripples as ensembles of oscillons. *Phys. Rev. Lett.* **108**, 034502 (2012).

13. Péринet, N., Gutiérrez, P., Urra, H., Mujica, N. & Gordillo, L. Streaming patterns in Faraday waves. *J. Fluid Mech.* **819**, 285–310 (2017).

14. Fauve, S., Douady, S. & Thual, O. Drift instabilities of cellular patterns. *J. Phys. II* **1**, 311–322 (1991).

15. Pétrélis, F., Laroche, C., Gallet, B. & Fauve, S. Drifting patterns as field reversals. *Europhys. Lett.* **112**, 54007 (2015).

16. Francois, N., Xia, H., Punzmann, H., Ramsden, S. & Shats, M. Three-dimensional fluid motion in Faraday waves: Creation of vorticity and generation of two-dimensional turbulence. *Phys. Rev. X* **4**, 021021 (2014).

17. von Kameke, A., Huhn, F., Fernández-García, G., Muñuzuri, A. P. & Pérez-Muñuzuri, V. Double cascade turbulence and richardson dispersion in a horizontal fluid flow induced by faraday waves. *Phys. Rev. Lett.* **107**, 074502 (2011).

18. Gordillo, L. et al. Can non-propagating hydrodynamic solitons be forced to move? *Euro. Phys. J. D* **62**, 39–49 (2011).

19. Gordillo, L. & García-Ñustes, M. A. Dissipation-driven behavior of nonpropagating hydrodynamic solitons under confinement. *Phys. Rev. Lett.* **112**, 164101 (2014).

20. Gordillo, L. & Mujica, N. Measurement of the velocity field in parametrically excited solitary waves. *J. Fluid Mech.* **754**, 590–604 (2014).

21. Urra, H. et al. Localized Faraday patterns under heterogeneous parametric excitation. *Phys. Rev. E* **99**, 033115 (2019).

22. Douady, S., Fauve, S. & Thual, O. Oscillatory phase modulation of parametrically forced surface waves. *Europhys. Lett.* **10**, 309–315 (1989).

23. Martin, E., Martel, C. & Vega, J. M. Drift instability of standing Faraday waves. *J. Fluid Mech.* **467**, 57–79 (2002).

24. Martin, E. & Vega, J. M. The effect of surface contamination on the drift instability of standing Faraday waves. *J. Fluid Mech.* **546**, 203–225 (2006).

25. Agez, G., Clerc, M. G., Louvergneaux, E. & Rojas, R. G. Bifurcations of emerging patterns in the presence of additive noise. *Phys. Rev. E* **87**, 042919 (2013).

26. Barashenkov, I. V., Bogdan, M. M. & Korobov, V. I. Stability diagram of the phase-locked solitons in the parametrically driven, damped nonlinear Schrödinger equation. *Europhys. Lett.* **15**, 113–118 (1991).

27. Denardo, B. et al. Observations of localized structures in nonlinear lattices: domain walls and kinks. *Phys. Rev. Lett.* **68**, 1730–1733 (1992).

28. Kutz, J. N., Kath, W. L., Li, R.-D. & Kumar, P. Long-distance pulse propagation in nonlinear optical fibers by using periodically spaced parametric amplifiers. *Opt. Lett.* **18**, 802–804 (1993).

29. Longhi, S. Stable multipulse states in a nonlinear dispersive cavity with parametric gain. *Phys. Rev. E* **53**, 5520–5522 (1996).

30. Alexeeva, N. V., Barashenkov, I. V. & Tsironis, G. P. Impurity-induced stabilization of solitons in arrays of parametrically driven nonlinear oscillators. *Phys. Rev. Lett.* **84**, 3053–3056 (2000).

31. Clerc, M. G., Coulibaly, S. & Laroze, D. Parametrically driven instability in quasi-reversal systems. *Int. J. Bifurcation Chaos* **19**, 3525–3532 (2009).

32. Miles, J. W. & Benjamin, T. B. Surface-wave damping in closed basins. *Proc. R. Soc. Lond.* **297**, 459–475 (1967).

33. Gordillo, L. *J. Non-Propagating Hydrodynamic Solitons in a Quasi-One Dimensional Free Surface Subject to Vertical Vibrations*. Ph.D. Thesis (Universidad de Chile, 2012).

34. Coullet, P., Frisch, T. & Sonnino, G. Dispersion-induced patterns. *Phys. Rev. E* **49**, 2087–2090 (1994).
35. Haragus, M. & Iooss, G. *Local Bifurcations, Center Manifolds, and Normal Forms in Infinite-Dimensional Dynamical Systems* (Springer Science & Business Media, 2010).

36. Kumar, A. Soliton dynamics in a monomode optical fibre. *Phys. Rep.* **187**, 63–108 (1990).

37. Porsezian, K. & Nakkeeran, K. Optical solitons in presence of Kerr dispersion and self-frequency shift. *Phys. Rev. Lett.* **76**, 3955–3958 (1996).

38. Gedalin, M., Scott, T. C. & Band, Y. B. Optical solitary waves in the higher order nonlinear Schrödinger equation. *Phys. Rev. Lett.* **78**, 448–451 (1997).

39. Nakkeeran, K., Porsezian, K., Sundaram, P. S. & Mahalingam, A. Optical solitons in N-coupled higher order nonlinear Schrödinger equations. *Phys. Rev. Lett.* **80**, 1425–1428 (1998).

40. Liu, W.-J., Tian, B., Zhang, H.-Q., Xu, T. & Li, H. Solitary wave pulses in optical fibers with normal dispersion and higher-order effects. *Phys. Rev. A* **79**, 063810 (2009).

41. Louvergneaux, E., Szwaj, C., Agez, G., Glorieux, P. & Taki, M. Experimental evidence of absolute and convective instabilities in optics. *Phys. Rev. Lett.* **92**, 043901 (2004).

42. Coulibaly, S., Durniak, C. & Taki, M. *Spatial Dissipative Solitons Under Convective and Absolute Instabilities in Optical Parametric Oscillators* 1–27 (Springer Berlin Heidelberg, Berlin, Heidelberg, 2008).

43. Clerc, M. G., Fernandez-Oto, C., Garcia-Ñustes, M. A. & Louvergneaux, E. Origin of the pinning of drifting monostable patterns. *Phys. Rev. Lett.* **109**, 104101 (2012).

44. Edri, Y., Meron, E. & Yochelis, A. Spatial asymmetries of resonant oscillations in periodically forced heterogeneous media. *Physica D* **410**, 132501 (2020).

45. Edri, Y., Meron, E. & Yochelis, A. Spatial heterogeneity may form an inverse camel shaped Arnol'd tongue in parametrically forced oscillations. *Chaos* **30**, 023120 (2020).

46. García-Ñustes, M. A., Humire, F. R. & Leon, A. O. Self-organization in the one-dimensional Landau-Lifshitz-Gilbert-Slonczewski equation with non-uniform anisotropy fields. *Commun. Nonlinear Sci. Numer. Simul.* **96**, 105674 (2021).

47. Nicolaou, Z. G., Case, D. J., Wee, E. B. V. D., Driscoll, M. M. & Motter, A. E. Heterogeneity-stabilized homogeneous states in driven media. *Nat. Commun.* **12**, 1–9 (2021).

48. Bush, J. W. Pilot-wave hydrodynamics. *Ann. Rev. Fluid Mech.* **47**, 269–292 (2015).

49. Peyrard, M. & Dauxois, T. *Physique des Solitons* (Savoirs Actuels, 2004).

## Acknowledgements

The authors thank Milena Páez-Silva for improvements of the experimental setup. J.F.M. thanks Laurette Tuckerman for fruitful discussions. J.F.M. was funded by ANID through the FONDECYT POSTDOCTORADO grant No. 3200499. R.R.-A. was funded by PUCV and UTFSM scholarship N°065/2022. R.R.A. and M.A.G.-Ñ. acknowledge the support of PROGRAMA DE COOPERACIÓN CIENTÍFICA ECOS-ANID ECOS200006, and S.C., the Program ECOS-Sud C20E07. M.A.G.-Ñ. and L.G. were funded by ANID FONDECYT

grants No. 1201434 and 1221103. J.F.M. also acknowledges ECOS-Sud project No. C15E06 for supporting a collaboration visit at Université de Lille, France, where part of this work was conceived. Finally, the authors are grateful to Dicyt-USACH for providing the funds required for publication.

## Author contributions

Conceptualisation, J.F.M., S.C. and M.A.G.-Ñ.; formal analysis, J.F.M., S.C., M.T. and M.A.G.-Ñ.; experimental investigation, R.R.-A., L.G. and M.A.G.-Ñ.; methodology, J.F.M., S.C., L.G. and M.A.G.-Ñ.; software, J.F.M., R.R.-A. and S.C.; supervision, S.C., M.T., L.G. and M.A.G.-Ñ.; validation, S.C., M.T., L.G. and M.A.G.-Ñ.; writing—original draft, J.F.M.; writing—review and editing, R.R.-A., L.G. and M.A.G.-Ñ. All authors have read and agreed to the published version of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42005-023-01170-8.

**Correspondence** and requests for materials should be addressed to Juan F. Marín or Leonardo Gordillo.

**Peer review information** *Communications Physcis* thanks Lyes Kahouadji and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**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/.

© The Author(s) 2023