Five-Wave Resonances in Deep Water Gravity Waves: Integrability, Numerical Simulations and Experiments

In this work we consider the problem of finding the simplest arrangement of resonant deep-water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wavevectors K1+K2=K3 (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wavepackets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction in a symmetric configuration. Numerical simulations of the governing equations in natural variables using pseudospectral methods require the inclusion of up to 6-wave interactions, which imposes a strong dealiasing cut-off in order to properly resolve the evolving waves. We study the resonance numerically by looking at a target mode in the base triad and showing that the energy transfer to this mode is more efficient when the system is close to satisfying the resonant conditions. We first look at encountering plane waves with base frequencies in the range 1.32–2.35 Hz and steepnesses below 0.1, and show that the time evolution of the target mode’s energy is dramatically changed at the resonance. We then look at a scenario that is closer to experiments: Encountering wavepackets in a 400-m long numerical tank, where the interaction time is reduced with respect to the plane-wave case but the resonance is still observed; by mimicking a probe measurement of surface elevation we obtain efficiencies of up to 10% in frequency space after including near-resonant contributions. Finally, we perform preliminary experiments of encountering wavepackets in a 35-m long tank, which seem to show that the resonance exists physically. The measured efficiencies via probe measurements of surface elevation are relatively small, indicating that a finer search is needed along with longer wave flumes with much larger amplitudes and lower frequency waves. A further analysis of phases generated from probe data via the analytic signal approach (using the Hilbert transform) shows a strong triad phase synchronisation at the resonance, thus providing independent experimental evidence of the resonance.


Introduction
Historically, nonlinear resonant interactions in surface water waves have focused mainly on the so-called exact resonances, defined by the equations where N denotes the number of interacting waves, k j denote the wave-vectors, and ω j denote the frequencies where ω j = ω(k j ) is provided by a dispersion relation. Early theoretical developments for these interactions gave rise to resonance interaction theory, originated by Phillips (1960) and developed by several authors (see the reviews by [1,2], and references therein). The key working hypothesis of this theory is the smallness of the wave steepness, namely the product between wave amplitude and magnitude of the wavevector. A formal perturbation theory based on multiple-scale methods is thus developed which allows one to obtain the exact resonances as necessary conditions to avoid secular behaviour. Solutions to these exact resonance equations can be found in many cases. Case N = 3, called triad resonances, applies to capillary-gravity waves (and also to a variety of wave systems such as Charney-Hasegawa-Mima, quasi-geostrophic equations, internal waves, inertial waves, etc.); case N = 4, known as quartet resonances, applies to gravity waves (including deep water). The case N = 5 is relevant in the important context of gravity waves in one-dimensional propagation, which is the topic to be discussed in this paper. In this context, it was shown in [3] that the interaction coefficients for the 4-wave resonances vanish identically. Later on, in [4] it was demonstrated that 5-wave resonances can exist and their interaction coefficients are nonzero. Then, in [5] (see also [6]) the explicit calculation of all 5-wave resonances along with their interaction coefficients was performed. Once solutions to the exact resonance equations are found, the dynamics of the wave amplitudes (at dominant orders) is dictated by coupled nonlinear PDE systems. In some simple scenarios the equations are integrable and can be solved analytically via the inverse scattering transform, leading to recurrent behaviour, showing periodic exchanges of energy between modes. In some scenarios, however, particularly in cases where multiple resonant interactions are coupled, the ODE and PDE systems obtained are known to display chaotic behaviour [7][8][9][10].
For practical reasons, the multiple resonant interaction dynamics is discussed using stochastic approaches, quite successfully. Pioneered by Hasselmann [11][12][13] and Zakharov [14,15], the now established theory of wave turbulence describes the energy exchanges across scales due to resonant quartet interactions in the case of deep water surface gravity waves. A key assumption of this theory, in addition to the smallness of steepness and the limit of large box size, is an asymptotic closure of the hierarchy of cumulants [16], which leads to a set of evolution equations for the so-called spectrum variables, namely the individual quadratic energies of the spatial Fourier transforms of the original field variables. This approach assumes that the phases of the Fourier transforms do not play an important role in the dynamics of energy transfers. Interesting modern departures from these theories for water waves and other wave systems consider quasi-resonant scenarios [17,18] and finite-amplitude regimes in discrete wave turbulence for other systems (such as the barotropic vorticity equation and more recently in wave-mean flow models of solar cycle modulations), where the phases interact with the spectrum variables producing interesting effects, such as precession resonance [19][20][21].
Up to and including the papers by [22][23][24], and the review by [1], only two experiments studied resonant interactions between wavetrains with comparable initial amplitudes and only one experiment studied the relative roles of nonlinearity and randomness for a broad spectrum of waves. Since then, experiments on triad resonances in gravity-capillary waves were performed by [25,26]. In the case of gravity waves, experiments on quartet resonances based on earlier work on modulational instability [27] were conducted by [28,29]. Resonant interactions in the presence of an underwater current were investigated by [30]. Experiments on persistent wave patterns and steady-state resonant waves were performed by [31,32]. Finally, experiments on degenerate quartets of oblique waves not influenced by modulational instability were conducted by [33]. All these experiments have confirmed the resonant interaction theory in the case of small steepness (ka < 0.1 roughly). The latter experiments also confirmed the expected nonlinear corrections due to marginally off-resonance terms.
In this paper we consider the simplest experimental/numerical/theoretical setup to search for a resonance in deep water gravity waves in one-dimensional propagation. Namely, we look for a 5-wave resonance using the minimum possible number of wavevectors (along with their negatives, representing an experiment of encountering wavepackets). Such a minimal arrangement is obtained in terms of a triad of positive wavevectors: K 1 + K 2 = K 3 , K j > 0, where the wavevectors satisfy quite specific ratios. Of course, triads are not resonant in water gravity waves because the dispersion relation is subadditive. The quintet resonances obtained represent 3 → 2 processes (transforming three waves into two waves), and the interacting wavenumbers must contain positive and negative elements from the set {±K 1 , ±K 2 , ±K 3 }. Notably, the fact that minimal 5-wave resonances can be based on non-resonant triads was demonstrated by one of us in another one-dimensional resonant system (also with subadditive dispersion relation) of historical importance: The Fermi-Pasta-Ulam-Tsingou system [34,35].
In Section 2 we review the Hamiltonian theory in natural variables up to and including 6-wave interactions, which will be useful for numerical simulations. We then review the approach in normal-form variables, and briefly describe all possible resonant quintets obtained in the simplified scenario of encountering wavevectors based on non-resonant triads. We then consider the dynamical system for such a resonant configuration and demonstrate that it is not integrable in general, although it has several integrals of motion and in particular the amplitudes are bounded; by reducing to a symmetric scenario of encountering wavevectors with the same amplitudes one obtains an integrable system based on a homoclinic orbit.
In Section 3 we provide detailed results of our water-wave pseudospectral numerical simulations, in two main scenarios: Encountering plane waves and encountering wave packets. For each scenario only the Fourier modes with wavevectors ±K 2 and ±K 3 have initial energy in them. The energy transfer to the 'target' Fourier modes with wavevectors ±K 1 is therefore measured after a time period of nonlinear interaction, by defining an 'energy efficiency' as the proportion of the total energy that goes to these target wavevectors. Several numerical experiments are conducted, off resonant and at resonance, in order to quantify the peak of efficiency at the resonant case. As an alternative way to measure efficiency, closer to the experimental way, we consider a 'numerical probe measurement' corresponding to what the surface elevation at a given point along the numerical tank would look like. The resulting time series is analysed in terms of frequency and time using a wavelet transform. The resulting scaleograms are shown and an alternative efficiency is defined, this time in terms of the wavelet frequency spectrum. This provides an independent quantification of the resonance. We also show that the resonance can be diagnosed by looking at the coherence of the quintet phases over their time evolution.
In Section 4 we provide a detailed description of the linear tank experimental setup and calibration of the two wavemakers and the two measuring probes that measure surface elevation at two symmetrically situated points. We perform 9 experiments: In each experiment the wavemakers are oscillated with two frequencies each: f 2 and f 3 , corresponding to wavevectors ±K 2 and ±K 3 , respectively. We analyse the probes' time series (surface elevation) and produce frequency-time scaleograms as in the previous section, obtaining preliminary evidence that the resonance can be captured in this type of experiment.
Finally, in Section 5 we summarise the results and discuss the challenges posed, as well as possible improvements to experiments and extensions to other scenarios. In addition, we present independent evidence of the resonance by showing that phase synchronisation is strong at the resonance, in all scenarios studied: Numerical plane waves, numerical wave packets and experimental waves.

Fifth-Order Hamiltonian Theory in Natural Variables
We consider the evolution of the free surface of a water column (in one dimension) via the boundary conditions for the surface elevation ζ(x, t) and the velocity potential Φ(x, t). A widely used form of these conditions in the study of nonlinear phenomena on water waves is derived by expanding the vertical velocity at the free surface to fifth order [15,[36][37][38][39][40], leading to the system where Q n and R n are given by where g is the acceleration of gravity and the Fourier transform of the linear operator L is defined by: and h (> 0) is the water depth. Notice that in the above definitions for Q n and R n we introduced a useful convention on the use of parenthesis: Terms of order 1, 2, 3, 4, 5 in the field amplitudes are enclosed, respectively, by curly {}, square [], round () brackets, and so on, repeating the cycle. From here we restrict our analysis to the deep-water limit, h → ∞, leaving the finite-depth case for a subsequent work. In the deep-water limit the L-operator's action reduces to LΦ k (t) = −|k|Φ k (t).
Computationally we will assume periodic boundary conditions in the longitudinal direction with computational box size L, so that for all x, t ∈ R.
The Hamiltonian up to and including 6-wave interactions is given by Explicit expressions for this, order by order, can be found in [40]. The evolution Equations (1) are canonical Hamiltonian equations: We note in passing that there is a relation between the above formulation and the Hamiltonian formulation using Dirichlet-Neumann operators by [41].
Defining the Fourier components of surface elevation and velocity potential as we consider now the canonical transformation from natural variables (ζ k ,Φ k ) to the socalled normal variables a k . For this the quadratic part of the Hamiltonian must be diagonal in the normal variables and the Poisson bracket be preserved. The resulting transformation isζ Notice that whileζ k andΦ k satisfy reality conditions (ζ −k =ζ * k etc.) the normal variables a k do not. In fact, we may interpret a k and a −k as independent waves propagating in direction k and −k, respectively. Thus, the number of degrees of freedom is preserved under this linear transformation. The variables a k have dimensions L 3/2 T −1/2 .

The Simplest Arrangement of Five-Wave Resonances
Expressions for the Hamiltonian in normal variables a k are usually cumbersome. The usual approach is to further transform the system to the so-called normal form via near-identity canonical transformations relating a k to normal-form variables b k : where the constant coefficients A are provided explicitly in [6]. In these new variables, only resonant interactions are present by definition. For deepwater surface waves in one-dimensional propagation, it is known that the lowest-order resonances appearing in the normal-form Hamiltonian are five-wave resonances [4,42]. More specifically: Three-wave resonances are absent; four-wave resonances are present but the corresponding interaction coefficients vanish at the resonant manifold, with only Birkhoff four-wave resonances surviving; and five-wave resonances are present and most of them have nonzero interaction coefficients, with the essential requirement that the interacting wavenumbers have opposite signs [5].
The calculation of the normal form up to and including resonant five-wave interactions was performed by [4][5][6], with two important results:

1.
Explicit formulae for the five-wave resonant manifold, leading to five different types of resonance, numbered from (i) to (v).

2.
Closed-form expressions for the five-wave interaction coefficients evaluated at the resonant manifold. In particular, it was shown that these interaction coefficients are equal to zero in resonance types (iii) and (iv), and nonzero in the resonance types (i), (ii) and (v).
In recent results by one of us regarding the Fermi-Pasta-Ulam-Tsingou (FPUT) system with periodic boundary conditions [34,35], it has been demonstrated that five-wave resonances are effectively constructed out of three-wave "frequency-only" resonances, namely where the momentum condition for the three interacting waves is not satisfied. It is to be noted that in the FPUT system the dispersion relation is strictly subadditive, just like in deep water waves, so there are no exact three-wave resonances in either system. Therefore it comes as a surprise that a five-wave resonance can be based on a three-wave "mismatched" resonance, and especially with the momentum conditions involving only three different wavenumbers along with their negatives.
This motivated the question of this paper: Can a five-wave resonance be found in terms of three positive wavenumbers satisfying the momentum condition, K + K = K , along with their negatives? This question has an experimental appeal too, as its answer would provide a minimal scenario involving two opposing wave trains, each consisting of just three central wavenumbers.
• Second resonant quintet: along with its negative version (−K 2 ) + (−K 2 ) + K 3 = (−K 3 ) + K 4 . In terms of linear frequencies via the dispersion relation f (k) = g|k|/2π, this resonance reads f 2 : f 3 : f 4 = 1 : There is a third resonant quintet but it is of type (iv) in the notation of [5], where the interaction coefficient is equal to zero. This quintet could appear in numerical simulations using the third-order equations such as [36], so we list it anyway: • Third resonant quintet (vanishing interaction coefficient): along with its negative version (−K 2 ) + (−K 2 ) + (−K 2 ) = (−K 4 ) + K 2 . Notice that one of the wavenumbers (K 3 ), along with its negative, do not appear explicitly in this resonant quintet.
The proof that the five-wave resonances listed above are the only possible ones based on three wavenumbers and their negatives is elementary but long: See Appendix A for details. It is based on a direct implementation of the inequalities in the classification given in [5], except that we allow some of the inequalities to saturate by virtue of the degeneracy of the solutions (as clearly some frequencies must coincide since we use three basic frequencies only).

Non-Integrability of the Five-Wave Resonance in the Case of Encountering Waves
In this subsection we consider the normal form Hamiltonian up to five-wave resonances.
If we consider just one of the quintets from the first and second resonant quintets of the previous subsection, and assume all other modes in the system are initially zero, it is then possible to show that the normal form Hamiltonian up to five-wave interactions is integrable, as it will consist of Birkhoff-only four-wave terms plus one five-wave term only. The system has four degrees of freedom, and four constants of motion in involution: The Hamiltonian and three quadratic invariants.
In contrast, as soon as we allow for the corresponding "negative" quintet (and assume all other modes in the system are initially zero), we obtain two interacting quintets whose Hamiltonian is not integrable. The system has six degrees of freedom, but only four constants of motion: The Hamiltonian and three quadratic invariants. To see this, consider for example the first resonant quintet in the previous subsection along with its negative version. That is, we have the quintets It is useful to define the set of relevant wavevectors as C = {K 1 , K 2 , K 3 , −K 1 , −K 2 , −K 3 }, so that the normal form amplitudes satisfy b k := 0 ⇐⇒ k / ∈ C. Then, the Hamiltonian in the normal form variables can be obtained by gathering the analyses of the interaction coefficients done in [4,5]. We get (symmetric on the upper indices and the lower indices separately) is known explicitly on the resonant manifold. For our resonant quintet and its negative version we get (case (v)(a) in [5]) Introducing the notation b n = b K n , b −n = b −K n , for n = 1, 2, 3, we get the following explicit expression for the Hamiltonian in normal form: The equations of motion will not be shown explicitly; they are obtained simply by following the canonical formalism: Is this dynamical system integrable? Notice that the quadratic and quartic terms in the Hamiltonian depend on the moduli squared |b n | 2 of the normal form variables only, so these terms commute with any of the moduli squared (i.e., these terms do not change the energy of the evolving modes). In contrast, the quintic term depends on the quintet phases φ 1 as well as the moduli, and is thus responsible for the energy exchanges across modes. We can therefore state that the system has a reduced number of degrees of freedom: 6 moduli + 2 quintic phases = 8 degrees of freedom, the 4 remaining phases being obtainable by quadrature. This result and the structure of the normal form imply that there are 4 Manley-Rowe-type quadratic invariants: Thus, along with the Hamiltonian, there are 5 independent invariants, leading to an effective 3-dimensional first-order autonomous dynamical system. If we could find one more independent invariant then we could conclude that the system is integrable. It turns out that one can simplify the search as follows. A rigorous result, motivated by the question of integrability (and clearly unrelated to the actual physical situation of low steepnesses), considers the limit of large amplitudes |b n | 1, where we can neglect the quadratic and quartic terms in the Hamiltonian (18), leading to a homogeneous system. Then, if the original system has an extra independent invariant then this new homogeneous system has an extra independent invariant (other than the 5 invariants already found). So we focus on the homogeneous system for this question of integrability.
We can apply two ideas to study the homogeneous system. First, we can use an idea that was first applied in [43], whereby the Kovalevskaya exponents of a homogeneous dynamical system are computed, and then a theorem by Yoshida [44] is used in order to obtain the degree (in terms of the amplitudes) that a potentially new invariant of such a system would have [45]. We refer the reader to [43] for details of the method, as it is not relevant in our case because we could not establish in this way whether our system is integrable or not. On applying this method we found a Kovalevskaya exponent of 7/9, which would correspond to an invariant of degree 7/3 in the amplitudes. Unfortunately this method neither guarantees nor precludes the existence of an invariant: The result only states that if a new invariant existed, then it should be of degree 7/3 in the amplitudes.
The second idea is to take advantage of the fact that the homogeneous system is effectively 3-dimensional, so Poincaré cuts can be produced. If these clearly show the existence of chaotic regions then we can safely conclude that the system is chaotic. We can use the quadratic invariants to reduce the analysis to the plane spanned by the variables |b 1 | 2 and |b 3 | 2 . The remaining degree of freedom can be chosen to be the phase . Using a highly accurate numerical scheme that conserves the Hamiltonian up to O(10 −11 ) we integrate the system for several choices of initial conditions and generate a Poincaré cut on the (|b 1 | 2 , |b 3 | 2 ) plane at the instances when φ + ∈ 2πZ. The resulting cuts are plotted in Figure 1. The evidence for chaos is clear: We can visually identify sets of islands (representing quasi-periodic motion) which are surrounded by the usual chaotic sea.
We conclude that the homogeneous system is chaotic, and therefore the original system is not integrable.

Analysis of the Scenario of Encountering Waves in the Small-Steepness Case
In preparation for the numerical and experimental studies of the next sections, we consider a situation where the amplitudes are initially small so that wave steepnesses |k| 5/4 |b k |/(4g) 1/4 are small (typically of order 0.1 or less). In this situation, the truncated Hamiltonian (18) approximates well the system, as terms in it appear in a hierarchical way, with quadratic terms about 100 times larger than quartic terms, and with quartic terms about 10 times larger than quintic terms, and so on if we had continued to higher orders. Despite the smallness of the quintic terms in the Hamiltonian, the energy-exchange dynamics is still controlled solely by that term so it cannot be neglected. In fact, at resonance, which is the case when this Hamiltonian describes the system, it is possible to show that the quadratic part of the Hamiltonian is conserved separately, as it is a linear combination of the four Manley-Rowe invariants (19) and (20). As a result, the sum of the quartic part and the quintic part of the Hamiltonian is conserved separately, and these two parts differ by a factor of order 10 only. As these terms are not sign-definite it is possible to have significant transfers.
We are interested in a scenario where modes b ±2 and b ±3 are non-zero initially, with b ±1 = 0 initially. We expect that the resonant case will display active energy transfers towards the "target" modes b ±1 , in sharp contrast with the non-resonant cases, where there is no quintic term in the Hamiltonian and therefore no energy transfer amongst modes at this level of nonlinearity. Energy exchanges can be studied using the Manley-Rowe invariants along with Fjørtoft-like arguments, as follows. The Manley-Rowe invariants I ± from Equation (19) do not provide a lot of insight except that any exchanges between the already energetic modes b ±2 and b ±3 must be balanced. On the other hand, the Manley-Rowe invariants J ± from Equation (20) provide a clear scenario of active transfer, whereby energy in modes b ±2 can be efficiently transferred to modes b ±1 . While it is not possible to determine from this argument alone how much energy can actually be transferred, a rigorous analysis of the reduced case where initially |b 2 | = |b −2 |, |b 3 | = |b −3 |, |b 1 | = |b −1 | = 0 shows that it is possible to transfer all the energy from |b 2 |(= |b −2 |) to |b 1 |(= |b −1 |) by virtue of the conservation of J + (= J − ), leading to a homoclinic orbit with lim t→±∞ |b 1 | = |b 2 | t=0 / √ 3 while lim t→±∞ |b 2 | = 0.

Numerical Simulations of the Fifth-Order Governing Equations in Natural Variables
Here we present results of the numerical simulations of the full PDE (1) up to and including six-wave interactions in the natural variables, in the scenario of waves with central wavenumbers K 1 , K 2 , K 3 encountering their negative versions, namely the waves with central wavenumbers −K 1 , −K 2 , −K 3 . We will focus on demonstrating numerically the quintet resonance given by Equation (16).
The evolution Equation (1) are solved using a standard pseudospectral method, dealiasing such that n max = N/6 for resolution N and fourth order Runge-Kutta timestepping. In all the results that follow N = 2 14 and the timestep was chosen as ∆t = 0.0005. The code was validated against an existing third order code using the equations of [37], experimental data from [38,39] and monitoring conservation of the Hamiltonian.
We consider initial conditions which are as close as possible to the experiments described in Section 4. Numerically the permitted wavenumbers K in the domain are of the form K = n2π/L, n ∈ Z. At first it might seem sensible to choose L close to that of the real experiments (35 m; see Section 4), but in reality the pertinent and comparable behaviours should be captured using the same 'physical' wavenumbers, or frequencies, within a similar interaction time interval in any sufficiently large numerical domain. It is of more importance to account for the discreteness of this wavenumber spectrum and the fact that, at significant amplitudes, nonlinear corrections will adjust the precise location of the resonance in k. For this reason it is advantageous to use a large L which has the consequence of allowing smaller δk, and therefore a relatively higher fidelity grid in k space, but naturally necessitating higher resolutions to account for a certain spatial scale, and hence temporal frequency. To this end we choose L = 400 m and consider the scaled resonance condition n 1 = 16p, n 2 = 9p, with p = 50 which results in resonant wavenumbers corresponding to linear frequencies via the dispersion relation f (k) = g|k|/2π: To demonstrate excitation of f 1 or a ±K 1 at the point of resonance we conduct a series of numerical experiments which fix a ±K 2 and vary a ±K 3 as indicated in the subsequent plots, keeping the initial amplitudeŝ where β is an overall rescaling of the initial amplitudes and all other a k = 0 initially. Note that Equation (21) indicates that the phases of the waves are set initially to zero. Experiments, not shown, demonstrate that varying the phases has little qualitative effect on the results to follow, only adding some corrections, and maintaining the a K = a * −K symmetry leaves the results quantitatively identical.

Encountering Plane Waves
In order to diagnose growth of energy in the target mode we define a measure of efficiency using the normal variables, themselves computed from the natural ones by inverting (14) and (15): where H is the full Hamiltonian, defined in Equation (12). We show max t∈[0,500] E K 1 (t) against K 3 in Figure 2 for β = 1 and β = 0.75. A distinct, but narrow peak is observed near the predicted resonant K 3 (marked with vertical line). The efficiencies are around 1.5% at the peak for β = 1, with the time series of E K (t) shown in Figure 3 demonstrating sustained growth in mode a K 1 over approximately 250 s in the resonant case, and no growth off resonance. Of note is the width of the resonant peak in K 3 ; coarser sampling of the initial wavenumber (or frequency) space will likely miss the optimal value for the resonance in the simulations. Due to nonlinear corrections the peak does not appear precisely at K 3 = 1250 2π L but slightly displaced at K 3 = 1248 2π L ; namely, a difference in wavelength of δλ ≈ 0.5 mm or in frequency δ f ≈ 0.0018 Hz. Figure 3 also shows the time series for the efficiency (or amplitude) of a K 3 . This shows a decrease of energy at late time which is generic across all of the calculations and only moderated slightly by the resonance. To determine which other nonlinear interactions are important in this system we examine the wavenumber spectrum as a function of time, plotted in Figure 4. Aside from K 1 , K 2 , K 3 we can observe 2K 2 having energy from early time, and side-bands of K 3 namely K 3 ± δK with δK = 2K 2 − K 1 growing at late times due to modulational instability. This is the reason for the decrease observed in a K 3 . Importantly this is a background effect, present in all cases, and does not hinder the five-wave resonance under consideration.
To mimic the diagnostics available in a physical experiment we also seek evidence of resonant transfer in a time series of a numerical 'probe' measurement of the surface elevation, namely ζ(x p , t), where x p = L/2 + 2 m is the fixed location in the domain for the probe. The choice of x p is arbitrary in the plane wave case, however this choice is made to be in keeping with the experiments in Section 4. The resonance should be observed by an increase in the temporal linear frequency f 1 in this measure. In order to diagnose such growth we create a frequency-time figure using a wavelet transform in the form of a 'scaleogram'. With a duration of t = 200 s, sampling 10,000 points in time and performing a wavelet transform using Gabor wavelets with 48 oscillations, with a scale resolution of 11 frequency scales along with 64 voices per scale, we obtain Figure 5. Three cases are shown (at the same K 3 as the time series of Figure 3) which demonstrates, for f 3 ≈ 2.206 Hz, the growth in time of the signal with frequency f 1 , marked with a horizontal blue dashed line. A similar efficiency can be defined using this measure defined as the sum of the squares of the scaleogram signals in Figure 5, on a strip of width 0.05 f 1 about f 1 and between t ∈ [10,190], divided by the sum of squares over the whole range of frequencies, shown in Figure 6. From these data emerges a qualitatively similar result as for the spatial Fourier efficiency shown in Figure 2, with a distinct peak around K 3 = 1248 2π L .
(a) (b)      To finalise the plane-wave case we provide a piece of analysis that connects with the theoretical derivation of the Hamiltonian system in the resonant case. Recall that for the symmetric initial conditions we took, the resonant case is expected to be integrable. A useful integrability test is to plot variables that would normally behave independently, except when the system is near integrable. To this end we provide joint probability density functions over the simulation time t ∈ [0, 200 s] of the quintet phase φ Figure 7. The resonant case, f 1 = 2.206 Hz, shows a clear reduction of the system's dimension, evidencing the onset of a coherent, near-integrable regime which we attribute to the resonance. A careful investigation of the source of this coherence shows that it really arises due to the increase in amplitude felt by a K 1 at resonance. Off resonance the phase rotates with a greater deal of variability, since the complex amplitude is closer to the origin, than in the resonant case where a K 1 has larger absolute value and therefore rotates at a more consistent rate, given by the linear frequency ω K . Notice also that in the symmetrical "encountering waves" setup studied in our numerical simulations, the symmetry φ −j = −φ j is valid at all times. This leads to the identity φ + = φ 1 + φ 2 − φ 3 , indicating that the 5-wave resonance is directly related to a phase locking between the three basic waves.

Encountering Wave Packets
In order to determine how reproducible the above result using plane waves may be in a physical experiment, we adapt the initial condition of the previous section into two types of spatially localised, propagating initial data. We wish to examine the joint effects of the dispersion of waves from a source (wavemaker), the limited interaction time and the limited spatial interaction region, due to the finiteness of the domain. The first such initial condition we consider is a localised wave packet at the centre of the domain containing the same wave content as in the planar case of the previous section, i.e., with x 0 = L/2 and σ = L/4 giving a packet occupying half of the domain, composed of both left propagating waves (i.e., propagating towards −x direction) and right propagating waves (i.e., propagating towards +x direction). This situation is close to the plane-wave case: Setting progressively larger σ would interpolate towards plane waves. It is also analogous to starting measurements in an experiment once the waves have encountered and overlapped, but neglects any interactions which occur as they start to encounter. As the waves propagate away from the centre of the domain the region nonlinearly interacting shrinks until eventually the required waves are no longer coincident. This will be referred to as the 'A' configuration. The 'B' configuration has similar half-domain localised wave packets, but now two of them, which contain either left or right propagating waves in such a way as to have encountering waves at the centre, i.e., where the subscripts r and l are the right and left propagating packets, respectively, and the centre of each packet is x 0 = L/4 from the centre. The full initial condition is ζ = ζ r + ζ l and φ = φ r + φ l , see upper panel of Figure 8. This configuration will therefore capture more dispersion and temporal behaviour as the waves propagate across one another, but will have have a smaller overall region at any given time where the necessary waves are coincident. Between these two configurations we have convenient middle ground between the plane wave case and the full physical experiment, to help ascertain where differences are observed and perhaps show how to design better experiments in the future. We perform the same parameterised sweep of f 3 as in Section 3.1 to intersect the resonance condition and seek increased transfer efficiency. max E K 1 is shown in Figure 2 along side the plane wave case and time series are shown in Figure 9. The efficiency is now much lower but a broader peak persists, case A showing a stronger signal than B. Now the time series demonstrates, in both cases but less pronounced in case A, that growth into K 1 is more transient and less sustained than in the plane case, the peak in |a K 1 | occurring within the first 10-20 s.
The reduction in overall efficiency and reduction in sustained growth is, of course, to be expected since the interaction is now localised in space and is not globally instantaneous across the domain as in the plane wave case. In addition, although the initial local surface elevations are comparable the overall energy content in the modes is reduced by the localising kernel used to create the individual packets. Finally it should be noted that the initial spectrum is broadened by the localisation, this means that first there is a wider range of K 3 experiencing resonant (or quasi-resonant) transfer, but also it is possible for more than just the single target K 1 to experience resonant growth. We may take this into account by allowing more modes to contribute to the efficiency definition about the central K 1 . However a better approach, having established that some resonant interactions do persist is to return to the wavelet, scaleogram, wave probe analysis as described above and see if we can observe transfers using a local measure of the surface height.  Given the moderately stronger signal in the efficiency E K using the A initial condition, we analyse the probe measurements from this case. We now use two probe locations x p = L/2 ± 2m first to correspond to the physical experiments in the proceeding section, but also noticing that a symmetry breaking occurs, with a more definite peak in the wavelet based efficiency in the left hand probe at L/2 − 2m as shown in the scaleograms in Figure 10 and the 'sum of squares' efficiency in Figure 11. Figure 10 shows the same 3 cases as in Figure 9a and the initial frequency f 3 is in evidence. Naturally, due to the broadening of the spectrum and the shorter time window in which the waves will be encountering the probe, the scaleogram shows a far more scattered set of frequencies. As described the left hand probe shows a discernible brightness about the target f 1 line for f 3 = 2.204 Hz where the peak in E K 1 is also observed in Figure 2, however this is not easily seen in the right hand probe. This opens up some interesting questions about symmetry breaking and measurement for a physical experiment. The symmetry breaking presumably arises due to numerical error in preparing the localised initial condition and/or due to instability of the localised initial condition, including processes due to nonlinear interactions. There are other frequencies in these plots at which a coherent signal is observed; these are due to other nonlinear interactions, e.g., the side-band instability at K 3 ± δK described earlier, which we will not elaborate upon here.
From these results we may conclude that quite long interaction times and regions will be required to make definitive observations of a five-wave resonance of this type.

Experiments
We performed preliminary experiments in a water wave tank to search for this resonance. The experimental setup, shown in Figure 12, consists of a linear wave tank (length = 35 m) with identical right-angle wedges (wavemakers) on each end that are oscillated vertically. The water mass at rest has a depth of 0.55 m. Two probes are located symmetrically at a distance of 15.5 m from the corresponding nearest end. The distance between the probes is 4 m. See, in Appendix B, probe calibration curves and consistency tests.

Encountering Wavepacket Experiments
In our experiments the water depth is finite and therefore one should use the gravitywave dispersion relation 2π f (k) = gk tanh(kh), where g = 9.8 m s −1 is the acceleration of gravity and h = 0.55 m is the water depth. However, for the frequencies studied ( f ≥ 1.33 Hz) the relative error with respect to the infinite-depth approximation 2π f (k) = g|k| is less than 10 −5 so we can use this latter approximation. Equivalently, for the frequencies tested, deep water can be assumed. We perform 9 experiments. On each experiment, the two wavemakers operate with the same oscillation pattern: consisting of the sum of two monochromatic oscillations: The first one with a base fixed frequency f 2 = 1.33 Hz corresponding to wavevectors ±K 2 = ±7.131 m −1 and with a fixed amplitude of approximately A 2 ≈ 0.06/K 2 (so steepness is A 2 K 2 = 0.06), while the second one has a changeable frequency f 3 , depending on the experiment, and a changeable amplitude A 3 ≈ 0.12/K 3 again so that the steepness is fixed at 0.12. Table 1 summarises the range of values of wavemaker frequencies f 3 used in the 9 experiments, with the corresponding amplitudes A 3 used. Table 1. For each experiment, frequency f 3 along with the wavemaker oscillation amplitude A 3 , chosen so that A 3 K 3 ≈ 0.12, where K 3 = 4π 2 g −1 f 2 3 is the wavenumber. In the last column, the theoretical target frequency f 1 = gK 1 /2π that would be observed if wavenumber K 1 = K 3 − K 2 was produced. We highlight in boldface the resonant case f 1 : f 2 : f 3 = 4 : 3 : 5, corresponding to Equation (16 In each experiment, therefore, as a result of the wavemakers' activity, two encountering wave trains with central wavevectors K 2 , K 3 , −K 2 , −K 3 will approach the probes, and will eventually encounter each other, triggering nonlinear interactions. If the theory is correct, when K 3 /K 2 = 25/9 the five-wave resonance should produce waves with wavevectors ±K 1 = ±16K 2 /9, as explained in Equation (16). The frequencies are proportional to the square roots of the wavenumbers, and as explained in the text right after Equation (16), at resonance we would have f 1 : f 2 : f 3 = 4 : 3 : 5. Now, eventually these daughter waves would reach the probes, which would record a signal with frequency f 1,reso = 4 f 2 /3. Theoretically, this five-wave resonance should occur near Experiment 6, because in that experiment we have K 3 = 4π 2 g −1 f 2 3 = 19.85 m −1 and as K 2 = 7.131 m −1 , we would get K 3 /K 2 ≈ 2.78 which is very close to 25/9. Similarly, the corresponding frequency of the oscillations produced by the nonlinear resonant interactions, 78 Hz, would be very close to f 1,reso = 4 f 2 /3 ≈ 1.77 Hz. Probe measurements are started when the wavepackets reach the probes. For each experiment we perform measurements on each probe, at intervals of ∆t = 0.02 s and during a time of 81.9 s, corresponding to a total of N = 4096 data points per probe per experiment. Figure 13 shows the time series for experiment 7, which is the one that shows a stronger signal for the target frequency f 1 , as we will see when we analyse the scalograms. From the time series we can see that the two probes' measured oscillations are quite close and in phase: The dominant frequency in the plot is f 2 = 1.33 Hz corresponding to the wavevectors ±K 2 . Modulations about this main frequency correspond to a combination of frequencies f 3 and f 1 along with a spectrum of less important frequencies, including harmonics. The fact that the two probes' signals are not exactly the same is reminiscent of the symmetry breaking observed in the numerical case, and may be due to slight differences in the experimental conditions and the subsequent nonlinear amplification of these differences.
As the fastest waves are the ones with wavevectors ±K 2 (at the fixed frequency f 2 = 1.33 Hz in the experiment), one can assume these waves will reach the probes first. We recall that, in the experiments, the initial measurement time t = 0 corresponds precisely to the time when the waves reach the probes for the first time. As the group velocity of the waves with wavevectors ±K 3 is less than the group velocity of the waves with wavevectors ±K 2 , we expect to see no effect from the nonlinear interactions (and in particular no generation of ±K 1 wavevectors) in the early stages. To be more precise, a simple analysis of the five-wave resonance Equation (16) shows that nonlinearly-produced wavevectors ±K 1 could be measured by the probes from the moment the wavevectors ±K 3 reach the probes. This time can be calculated explicitly for each experiment, giving, for experiments 1 to 9, t 0 = {14. 44, 15.04, 15.44, 16.44, 17.14, 17.84, 18.42, 19.42, 20.42} s. The leftmost vertical marker line in Figure 13 corresponds to this "first measurement" time.
It is also evident from Figure 13 that wave amplitudes become larger near t ≈ 60 s. This is due to the reflected wavepackets with wavevectors ±K 2 , which after reflection from the walls of the tank arrive back at the probes. After this time we expect a stronger signal for the nonlinearly-produced wavevectors ±K 1 in the probes. The time this happens is t 1 = 59.38 s, independent from the experiment. The rightmost vertical marker line in Figure 13 corresponds to this "second measurement" time.
In order to understand the frequency content of these signals, we apply a wavelet transform on them. Wavelet scaleograms show the time evolution of the energy content over a range of frequencies. By trial and error we found that a good compromise is struck between frequency resolution and time resolution in the wavelet transform by choosing a Gabor wavelet with 96 oscillations and with the usual scale resolution of log 2 (4096/2) = 11 and a total of 256 voices per scale. Before analysing the real data, Figure 14 shows densityplot scaleograms corresponding to artificial signals containing only four frequencies: where f 3 is chosen from the list in Table 1 and f 1 , f 4 are two possible target frequencies defined in terms of the dispersion relation by 2π f 1 = gK 1 and 2π f 4 = gK 4 , where K 1 = K 3 − K 2 and K 4 = K 3 + K 2 are target wavevectors that could be generated by nonlinear interactions in the experiments (notice that 2π f 3 = gK 3 defines the wavevector K 3 ). Apart from f 2 = 1.33 Hz which is fixed, these frequencies depend on the experiment (but are completely determined by it) according to Table 1. The scaleograms show spurious boundary effects at early and late times for the frequency f 3 . For the frequencies f 1 and f 4 the boundary effects are minimal. The plots show that these two frequencies are resolved within a band of width ∆ f = 0.09 f . Thus, a real signal from the experiments at a given frequency will be expected to have a comparable width. This will be useful when estimating the relative energy in a frequency.
In Figures 15-17 we show the scaleograms obtained from probe 1 data (left panels) and probe 2 data (right panels). For each panel, vertical lines analogous to those in Figure 13 indicate the theoretical first time t 0 when nonlinearly produced target frequency f 1 due to the five-wave resonance can be measured by the probes (leftmost line), the theoretical time t 1 when the wavepackets arrive back at the probes after reflecting from the tank walls (middle line), and the time beyond which the scaleograms show spurious boundary effects, as evidenced in Figure 14 (rightmost line). The figures also show horizontal lines corresponding to the theoretical frequencies f 1 , f 3 , f 4 derived from the wavevectors ±K 1 , ±K 3 , ±K 4 , and a theoretical harmonic f 1 + f 2 . It is apparent that experiment 6 shows some activity in the f 1 line, as well as the line f 1 + f 2 . However, the next experiment (number 7) shows a more convincing signal at f 1 , apparently resolved, particularly at the late stages after t = t 1 = 59.38 s.
In order to provide a more quantitative measure of the signals shown by the scaleograms, we calculate the relative energy of the wavelet transforms that are plotted in the scaleograms, for the frequency f 1 . We do this by defining an early-time energy by the sum of squares between t = t 0 and t = t 1 = 59.38 s, and on a strip of width ∆ f = 0.09 f 1 about f 1 . We then divide this energy by the total energy over the time interval [t 0 , 59.38 s] and over the whole frequency range (including f 2 which is not shown). The result is an "early-time efficiency" which we plot in Figure 18, dashed lines, for each probe. Similarly, we define a late-time energy by the sum of squares between t = t 1 = 59.38 s and t = 77.92 s, and on a strip of width ∆ f = 0.09 f 1 about f 1 . We then divide this energy by the total energy over the time interval [59.38 s, 77.92 s] and over the whole frequency range (including f 2 , not shown). The result is a "late-time efficiency" which we plot in Figure 18, solid lines, for each probe. It is evident from these plots that there is an efficiency peak that stands out, and although it is modest (0.5% efficiency) it nevertheless shows that considering the short times involved there is scope for a measurement of the theoretical five-wave resonance, perhaps in a larger tank and with larger amplitudes so that the effect can be clearly obtained.

Conclusions and Discussion
In this work we proposed to study the minimal resonant configuration of water gravity waves in one-dimensional propagation. The resulting 5-wave resonances are based on a triad of wavevectors along with their negatives, leading to a scenario of encountering wavepackets, which has an experimental appeal. The problem poses interesting challenges in all fronts: • On the theoretical front, the 5-wave resonant manifold and the interaction coefficients on this resonant manifold had been obtained in a series of papers over two decades ago [4][5][6], triggered by the discovery that the interaction coefficients vanish identically on the 4-wave resonant manifold for one-dimensional propagation of water gravity waves [3]. We used these results to find the simplest 5-wave resonance that can be made out of a triad of wavevectors and their negatives, and calculated its normalform Hamiltonian. We proved that the system is not integrable, but it lacks just one constant of motion to become integrable, so symmetric scenarios can produce integrable systems. • On the front of numerical simulations of the governing partial differential equations, the equations had been obtained over five decades ago [15], as a power series in terms of the steepness effectively. The numerical implementation we needed to use in order to accurately resolve 5-wave interactions is the one that uses up to 6-wave interactions in a pseudo-spectral setting [40]. Such an implementation requires a higher-thanusual dealiasing and thus can get quite expensive in terms of the required spectral resolution. We managed to validate our implementation against some benchmarks, and were able to establish the existence of the resonance and quantify its effects in terms of the energy transferred to target modes. We considered encountering plane waves and also encountering wave packets, in a simulated tank that is 300 m long. The plane-wave case provided the most efficient energy transfer in terms of Hamiltonian energy, while the wavepackets case provided a higher efficiency in terms of the probe measurements of surface elevation at a point along the tank.
• On the experimental front the main difficulty is to fine tune the amplitudes and frequencies in order to capture the resonance, but we got this from hindsight. Our preliminary experiments seem to show that the resonance exists physically, although as can be seen in Figure 18, the efficiency is relatively small.
We learned from our more recent numerical experiments that the resonance peak has a narrow width. For this reason we propose to conduct experiments in much longer wave flumes with much larger amplitudes of lower frequency waves. In addition, the study of the second resonant quintet (17) obtained theoretically, and/or the exploration of optimal initial conditions, may make it easier to find a resonance in experiments as well as numerical simulations. In future work we plan to extend this type of minimal search to gravity-capillary waves and to gravity waves in two-dimensional propagation.
We would like to end on a positive note, by reporting on the following phase-locking analysis between the three scenarios studied: Numerical encountering plane waves, numerical encountering wave packets, and experimental encountering waves. We will establish that the triad phases get locked near our theoretical resonance, thus providing an independent piece of evidence that the resonance is indeed taking place. We explain briefly how this is done. Using the probe data for each scenario and each "experiment" within a given scenario, we perform pass-band filters around each of the three important frequencies f 1 , f 2 , f 3 , with width 0.09 f . Each one of these three filtered signals is then extended via the analytic signal approach using a Hilbert transform (see, for example, [46]) in order to get the "conjugate" signal and thus calculate the "analytic phases" denoted by the time series ϕ 1 (t), ϕ 2 (t), ϕ 3 (t). These phases contain the fast time dependence 2π f j t, respectively. After subtracting this fast time dependence, obtain the new phases ϕ j (t) = ϕ 1 (t) − 2π f j t, j = 1, 2, 3. We then calculate the triad phase The phase-locking value for this signal is defined by (see, e.g., [47]): where · t denotes time average over the whole time series, and | · | denotes modulus. A phase-locking value close to 1 denotes a situation where the corresponding phase is locked at a given value for most of the time. The resulting phase-locking values, for the different scenarios and the different experiments, are plotted in Figure 19, using similar notation as in the efficiency Figures 6, 11 and 18. It is evident that in each scenario the phase-locking value is maximum precisely at the resonant case (denoted in each panel with a vertical dashed line), thus confirming that the 5-wave resonance not only provides larger target amplitudes but also strong phase synchronisations amongst the interacting modes.
Notice that the effect of symmetry breaking is again patent in the numerical wave packets case, cf. Figure 19b, as well as in the experimental case, cf. Figure 19c. In each case, one of the probes shows a clear global maximum of phase-locking R at the resonance, while the other probe has a local maximum at the resonance and a global maximum at another place. This behaviour is consistent with the efficiency plots in Figures 11 and 18: The probes that show a higher efficiency peak in these figures, correspond precisely to the probes that show a clear and strong maximum of the phase-locking value R in Figure 19b,c.  Figure 19. Phase-locking values of analytic-signal-approach triad phases, Equations (30) and (31), obtained from probe data, for different experiments within the three scenarios studied in this paper: (a) Encountering numerical plane waves; (b) encountering numerical wave packets; (c) encountering experimental waves. In all three scenarios, a clear local maximum of the phase-locking value is attained at the predicted resonance, corresponding to the maximum in efficiency for each case, cf. Figures 6, 11 and 18, respectively. In all panels, the dashed grey vertical line corresponds to the theoretical resonant case Funding: R.I. is partially supported by the Bulgarian National Science Fund, grant KΠ -06H42/2 from 27.11.2020.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author (miguel.bustamante@ucd.ie).
(v.2.b) k 3 = −J 1 . Here the momentum condition reads 2J 2 + (−J 1 ) = J 1 + (−J 3 ), which simplifies to J 1 = 3J 2 , thus J 3 = 4J 2 . We turn to the resonance condition to find 2 √ J 2 + √ 3J 2 = √ 3J 2 + √ 4J 2 , which is satisfied. In summary this leads to a 5-wave resonance parameterised by K ∈ Z + as follows: K + K + (−3K) = 3K + (−4K), Figure A1 shows a picture of the probes arrangement in the experimental tank. Figure A1. Picture of the experimental tank with a close up on the installed probes. The wire of each capacitance probe is directly facing each wedge. Figure A2 shows the results of the calibration of the probes. A linear regression gives R 2 = 0.999 for the relation between surface elevation and voltage, over a wide range of surface elevations (between −9 and 9 centimetres).  Figure A3 shows the results of several cases of monochromatic wavemaker oscillations over a range of frequencies to be used in the interacting experiment. On each panel, solid black lines show the wave surfaces measured by probe 1 when the north wedge oscillates monochromatically, while red dashed lines correspond to probe 2 measurements when the south wedge is oscillated. Apart from slight differences between the two probes in panels (c) and (d) ( f = 2.22 Hz and f = 2.3 Hz, respectively), the waves are stable for the range of frequencies studied. It was found (figure not shown) that for frequencies above f = 2.35 Hz it was difficult to get regular wave trains, so we do not consider wavemaker frequencies higher than f = 2.35 Hz.