🔢 mathematics

Sheets of Spectral Data of Stokes Waves in Weakly Nonlinear Models

This paper introduces a perturbation method to analytically derive sheets of spectral data for small-amplitude Stokes waves in weakly nonlinear unidirectional models, unifying the treatment of high-frequency and Benjamin–Feir instabilities (including those with discontinuous dispersion relations) and providing the first analytic approximation of the Benjamin–Feir spectrum for this class of models, which is validated against numerical computations.

Original authors: Benjamin Akers, Ryan P. Creedon

Published 2026-03-17

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Benjamin Akers, Ryan P. Creedon

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are standing by a calm lake. You toss a pebble in, creating a single, perfect ripple. Now, imagine you could keep adding energy to create a long, steady train of waves moving across the water. In physics, these are called Stokes waves.

For nearly 200 years, scientists have studied these waves. But here's the catch: even if a wave looks perfectly stable, it might actually be a "ticking time bomb." If you nudge it just right, it can suddenly break apart or grow chaotic. This paper is about figuring out exactly when and how these waves decide to explode into chaos.

Here is the breakdown of the research, explained without the heavy math jargon.

1. The Setup: The "Flat" State vs. The "Wave" State

Think of the water when it's perfectly still. In math, this is the "flat state." If you shine a light through it, the light passes straight through.

Now, imagine a wave starts moving. The researchers are looking at what happens when you take that flat state and add a tiny bit of energy (a small wave). They want to know: Does this tiny wave stay tiny, or does it grow into a monster?

To answer this, they use a "spectral map." Imagine a map where every point represents a possible way the wave could wiggle.

Safe Zone: Points on the map where the wave stays calm.
Danger Zone: Points where the wave grows unstable and breaks.

2. The Discovery: "Islands" of Chaos

The most exciting thing the authors found is that the "Danger Zones" aren't just random blobs. They are shaped like islands floating in a sea of safety.

The "Isola": In the math world, they call these islands isolated curves (or isolas). Imagine a tiny, elliptical island in the middle of a calm ocean. If your wave's energy lands on this island, it becomes unstable.
The "Figure-Eight": There is a special, famous type of instability called the Benjamin-Feir instability. On their map, this looks like a figure-eight shape (a lemniscate) right in the center. This is the classic way waves break up.

3. The New Tool: The "Sheet" of Data

Before this paper, scientists usually looked at these islands one slice at a time, like looking at a single slice of a loaf of bread. They would pick a specific wave size and see what the island looked like.

The authors developed a new method to look at the whole loaf at once.

The Analogy: Imagine the islands of instability are actually 3D sheets of paper floating in space.
The Magic: By changing a single "knob" (a parameter they call $\beta$ , which is like the ratio of two different wave frequencies mixing together), they can slide their view up and down this sheet.
The Result: This allows them to see the entire shape of the instability, not just a snapshot. They can see how the islands grow, shrink, and move as the wave gets bigger.

4. Two Types of "Explosions"

The paper identifies two main ways these waves can go unstable, which they call Triads and Quartets.

Triads (The Quick Spark):
- What it is: A collision between three wave components.
- The Analogy: Like a spark hitting a pile of dry leaves. It happens very fast and is very strong.
- The Shape: These create small, round islands that appear very close to the original wave's energy. They grow quickly (proportional to the wave size, $\epsilon$ ).
- Who wins? If these exist, they are the "bosses" of the instability. They dominate everything else.
Quartets (The Slow Burn):
- What it is: A collision between four wave components.
- The Analogy: Like a slow-burning fuse. It takes longer to kick in and is weaker than the Triad.
- The Shape: These islands are smaller and drift further away from the center. They grow slower (proportional to the square of the wave size, $\epsilon^2$ ).
The Benjamin-Feir (The Classic Figure-Eight):
- This is the old-school instability everyone knew about. The paper shows that if the "Triad" and "Quartet" islands aren't there, the Figure-Eight takes over. But if the Triads are present, they usually crush the Figure-Eight.

5. The "Broken" Wave Equation

Most math models assume the water behaves smoothly, like silk. But the authors also looked at a model called the Akers-Milewski equation, which has a "kink" or a "jump" in it (like a step in a staircase).

The Challenge: Standard math tools break when they hit a step. You can't take a smooth derivative of a step.
The Fix: The authors invented a new way to do the math that works even when the rules suddenly jump. They successfully mapped the instability for this "broken" model, proving their method is robust enough to handle rough, real-world scenarios that other theories can't touch.

6. Why Does This Matter?

You might ask, "Who cares about tiny islands on a math map?"

Predicting Tsunamis and Storms: Understanding exactly how waves break helps us predict when a calm ocean will suddenly turn violent.
Better Models: By comparing different mathematical models (like the Whitham equation vs. the Kawahara equation), they can tell engineers which model is best for predicting real water behavior.
The "Sheet" Concept: By seeing the whole 3D sheet of data, scientists can now predict how a wave will behave before it actually breaks, rather than just reacting to it.

Summary

In simple terms, this paper is like a weather forecast for water waves.

They found that waves don't just break randomly; they break on specific, predictable "islands" of chaos.
They built a 3D map (a "sheet") to see all these islands at once, rather than just looking at one slice.
They discovered that some instabilities are "fast and furious" (Triads) while others are "slow and steady" (Quartets), and they figured out which one wins.
They proved their math works even on "rough" models that other scientists couldn't solve.

It's a major step forward in understanding the hidden rules that govern the ocean's most dangerous waves.

1. Problem Statement

The paper investigates the spectral stability of small-amplitude Stokes waves within a general family of weakly nonlinear, unidirectional dispersive models governed by the equation:
$u_t + Lu + (u^2)_x = 0$
where $L$ is a Fourier multiplier with symbol $i\omega(k)$ , representing the linear dispersion relation.

The primary focus is on understanding how instabilities arise from collisions of eigenvalues in the linearized spectrum about the flat state ( $u=0$ ). Specifically, the authors analyze two distinct types of instabilities:

High-Frequency Instabilities: Arising from collisions of non-zero imaginary eigenvalues at non-zero Floquet exponents ( $p_0 \neq 0$ ). These manifest as "isolated" loops of instability (isolas) in the complex spectral plane.
Benjamin–Feir (Modulational) Instability: Arising from the collision of eigenvalues at the origin ( $\lambda=0, p_0=0$ ). This typically manifests as a "figure-eight" shaped instability region.

The authors aim to unify the analysis of these instabilities across various models (including those with discontinuous dispersion relations) and provide a direct comparison of their growth rates.

2. Methodology

The authors employ a perturbation method expanding the spectral data (eigenvalues $\lambda$ , eigenfunctions $v$ , and Floquet exponent $p$ ) in powers of the wave amplitude $\epsilon$ (where $0 < \epsilon \ll 1$ ).

Key Methodological Innovations:

Free Parameterization via Mode Ratio ( $\beta$ ): Unlike previous approaches that expanded the ratio of colliding eigenmode amplitudes as a power series in $\epsilon$ $ϵ$ , this paper treats the ratio $\beta$ $β$ of the amplitudes of the colliding flat-state eigenfunctions as a free parameter.
- By varying $\beta$ and $\epsilon$ , the authors generate 2D sheets of spectral data in the $(\text{Re}\lambda, \text{Im}\lambda, \epsilon)$ space.
- Slices of these sheets at fixed $\epsilon$ recover the known isolas and figure-eight curves.
Unified Framework: The same perturbation scheme is applied to both high-frequency collisions (non-zero $p_0$ ) and the Benjamin–Feir collision ( $p_0=0$ ).
Handling Discontinuous Dispersion: The method is adapted to handle models where the dispersion relation $\omega(k)$ is not differentiable (specifically the Akers–Milewski equation). This involves careful tracking of singular terms arising from jump discontinuities in $\omega(k)$ at $k=0$ , rather than relying on Taylor expansions of $\omega$ .
Validation: Asymptotic predictions are rigorously validated against numerical spectra computed via:
- Floquet–Fourier–Hill (FFH) method.
- Quasi-Newton continuation methods for individual eigenpairs.

3. Key Contributions and Results

A. Classification of High-Frequency Instabilities

The paper categorizes high-frequency instabilities based on the resonant interaction between the carrier wave ( $k=1$ ) and the colliding modes ( $k_1, k_2$ ). Let $m = k_1 - k_2$ .

Triad Resonances ( $|m|=1$ ):
- Mechanism: Interaction between modes where $|k_1 - k_2| = 1$ .
- Result: Produces elliptical isolas centered at the flat-state eigenvalue $\lambda_0$ .
- Scaling: The semi-major axes scale as $O(\epsilon)$ .
- Condition: Instability occurs if the Krein signature condition is met (necessary and sufficient). The authors show that the solvability conditions reduce to a linear system, simplifying the algebra compared to previous nonlinear approaches.
Quartet Resonances ( $|m|=2$ ):
- Mechanism: Interaction where $|k_1 - k_2| = 2$ .
- Result: Produces elliptical isolas, but the centers drift from $\lambda_0$ along the imaginary axis.
- Scaling: The drift and semi-major axes scale as $O(\epsilon^2)$ .
- Condition: Also governed by the Krein signature.
Higher-Order Resonances ( $|m| \geq 3$ ):
- Result: No instability is generated at $O(\epsilon)$ or $O(\epsilon^2)$ . Instability is expected to appear at $O(\epsilon^{|m|})$ .

B. Benjamin–Feir Instability and Discontinuous Dispersion

Lemniscate Approximation: The authors derive an analytic approximation for the Benjamin–Feir spectrum, showing it forms a lemniscate (figure-eight) shape.
Growth Rate: The instability grows as $O(\epsilon^2)$ .
Discontinuous Models: The method successfully derives the Benjamin–Feir spectrum for the Akers–Milewski equation, which has a dispersion relation $\omega(k) = \text{sgn}(k)(1+|k|)^2$ $ω (k) = sgn (k) (1 + ∣ k ∣)^{2}$ with a jump discontinuity at $k=0$ $k = 0$ . This is the first analytic approximation of the Benjamin–Feir spectrum for a model with discontinuous dispersion.
- For this specific model, the imaginary component of the instability vanishes at $O(\epsilon)$ and $O(\epsilon^2)$ , requiring an expansion to $O(\epsilon^3)$ to capture the leading-order behavior.

C. Comparative Analysis of Instability Strengths

By using a unified framework, the authors directly compare the growth rates of different instabilities:

Triad vs. Benjamin–Feir: Triad instabilities ( $O(\epsilon)$ ) always dominate Benjamin–Feir instabilities ( $O(\epsilon^2)$ ) at small amplitudes if triad resonances exist.
Quartet vs. Benjamin–Feir: In models without triad resonances, the competition depends on model-specific coefficients. The authors provide criteria to determine whether the quartet instability or the Benjamin–Feir instability is the dominant mechanism.
No High-Frequency Instabilities: If a model admits neither triad nor quartet instabilities, the Benjamin–Feir mechanism becomes the leading source of instability.

4. Significance

Unification: The paper provides a single, coherent perturbation framework that handles both high-frequency and modulational instabilities across a broad class of dispersive models.
Analytic Insight into Discontinuous Systems: It breaks the barrier of requiring smooth dispersion relations, extending rigorous spectral analysis to models like Akers–Milewski, which are crucial for modeling capillary effects and other physical phenomena where smoothness fails.
Predictive Power: The "sheet of spectral data" concept offers a comprehensive view of how instability regions evolve with amplitude and Floquet exponent, rather than just providing isolated snapshots.
Validation: The high degree of agreement between the asymptotic predictions and numerical computations (including error analysis showing $O(\epsilon^3)$ convergence) confirms the robustness of the derived formulas.
Future Directions: The bifurcation structure revealed serves as a foundation for future rigorous existence proofs of unstable spectra and guides the study of nonlinear dynamics near the onset of instability.

In summary, this work advances the theoretical understanding of water wave stability by providing a generalized, parameterized method to predict and compare the onset of various instabilities, successfully bridging the gap between smooth and discontinuous dispersive models.