PaperlayerView exact revision on arXiv

Source-linked AI summary

Robustness of mKdV Breathers: A Numerical Perspective

Chandler Haight, Svetlana Roudenko, Diana Son, Kai Yang

arXiv:2609.27080v1nlin.PScs.NAmath.APmath.NA

TL;DR

The paper asks how mKdV breathers behave under perturbations and numerically studies perturbed breathers, interactions, and changes to the equation’s nonlinearity. It finds stability under small perturbations, resolution into coherent structures under larger ones, and breather-like persistence in some non-integrable models.

  • Problem

    The paper asks how perturbations of the nonlinear term affect the existence and stability of breather-type solutions.

    Source: §5.1 · p. 25

  • Method

    The authors numerically perturb and evolve mKdV breathers, extracting asymptotic parameters and local masses to characterize coherent states and radiation.

    Source: §6 · p. 29

  • Results

    Small perturbations leave mKdV breathers stable, while larger ones yield breather or soliton structures with leftward radiation; breather-like dynamics also persist in non-integrable mKdV-type models.

    Source: §1 · p. 2, §6 · p. 29

  • Takeaways & Limitations

    The numerical observations motivate rigorous analysis of breather stability and long-time asymptotics.

    Source: §6 · p. 30

    The L∞ norm may not reveal breather periodicity, especially when α is much larger than β.

    Source: §2.2 · p. 7

Abstract

from arXiv · show

We systematically investigate breather solutions in the modified Korteweg-de Vries (mKdV) equation via numerical simulations. We show that the breather solutions are stable under a variety of perturbations, including amplitude changes, modifications of internal parameters, and perturbations of the nonlinearity of the equation. Our results are not only consistent with the known analytical studies on the stability of breathers in integrable systems, extending them further, but also provide numerical evidence that such breather-type structures can persist and remain stable over the simulated time scales in some non-integrable settings with symmetric potentials.

1. Introduction

The paper numerically examines mKdV breathers, their interactions, and perturbations of both the solutions and the equation’s nonlinearity. It reports stability under small perturbations, resolution into coherent structures under larger perturbations, and stable breather-like behavior in a non-integrable setting.

1. Introduction

  • A breather is a spatially localized, time-periodic oscillation under a solitary-wave-shaped envelope, moving in one direction at constant speed.Its oscillatory profile distinguishes it from a soliton, while both are localized and travel in a specific direction.Source: §1 · p. 1
  • Wadati wrote the studied mKdV breather family using integrable-systems methods, following early discussions of breathing solitons and non-purely imaginary poles.The inverse scattering method was introduced in.Source: §1 · p. 1
  • Breathers have received less study than solitons, although recent research addresses their properties, stability, interactions, applications, and behavior in non-integrable systems.Examples of applications include fluid dynamics, oceanic models, dispersion-managed systems, and nonlinear wave equations.Source: §1 · p. 2
  • Two perturbed forms add lower-order nonlinearities: a generalized-Gardner variant and an absolute-value variant.The paper considers a small added nonlinear term in each form.Source: §1 · p. 2
  • The paper reviews mKdV breathers and develops numerical identification by internal parameters for their long-term evolution, including interactions and multi-breathers.It also studies breather-soliton interactions and perturbed mKdV models.Source: §1 · p. 2
  • Small perturbations leave mKdV breathers stable, while larger ones yield breather or soliton structures with leftward radiation, supporting the soliton-breather resolution conjecture.The study also reports breather-like stability in a non-integrable mKdV-type model with a symmetric potential.Source: §1 · p. 2

2. Review of solitons and breathers in the mKdV equation

The mKdV equation supports traveling solitons and localized, time-periodic breather solutions. This review relates breather parameters to their motion and oscillations, verifies numerical profiles against exact solutions, and shows why the L∞ norm alone can obscure breather periodicity.

2. Review of solitons and breathers in the mKdV equation

  • Breathers are localized in space and periodic in time; γ determines envelope velocity, while δ sets the internal oscillation rate.The physical velocity is −γ; γ > 0 corresponds to leftward motion.Source: §2 · p. 3

2.1. Pokhozhaev identities.

  • Pokhozhaev identities and the mass–energy connection are useful for numerical testing and checking scheme accuracy.The review identifies the breather parameters from mass and energy: mass gives β, then energy gives α.Source: §2.1 · p. 4, §2 · p. 4

2.2. Concise review of mKdV breathers.

  • Increasing α from 1 to 2 at β = 1 increases the internal oscillation rate and the distance traveled over the same time.For B1,1, δ = −2; for B2,1, δ = 1.Source: §2.2 · p. 6
  • Reducing β from 1 to 0.5 at α = 1 makes the breather wider and its envelope oscillations slower.For B1,0.5, δ = 1.Source: §2.2 · p. 6
  • The simulated breather profiles nearly match the exact profiles, supporting the numerical scheme’s accuracy.The examples include matching the B1,1 profile at t = 1.58 with an exact shifted breather.Source: Figure 1 · p. 5, §2.2 · p. 6
  • For B3,0.2, the L∞ norm looks nearly constant despite small oscillations, so that norm alone cannot distinguish a breather from a soliton.The oscillations become visible on zooming in to changes in the third decimal place.Source: §2.2 · p. 7, Figure 5 · p. 7

3. Interactions of multiple mKdV coherent structures

The section tests forced mKdV interactions involving breathers traveling left, including breather pairs, triples, and soliton combinations. Across the simulations, interacting structures recover their original profiles, with no radiation or dispersion observed at the numerical resolution used.

3. Interactions of multiple mKdV coherent structures

  • Across the simulated interactions, coherent structures recovered their original shapes and speeds without radiation or dispersion at the numerical resolution used.Source: §3 · p. 8, §3.2 · p. 10
  • After temporary distortion in the soliton-breather collision, both structures returned to their original size and shape and continued in opposite directions.No visible radiation appeared at the numerical resolution used.Source: Figure 6 · p. 8, §3.1 · p. 8

3.2. Breather-Breather Interaction.

  • For two left-traveling breathers with equal β, differing α values produced different speeds and internal oscillations, while their widths and amplitudes matched.The breathers separated after collision without observed radiation or dispersion.Source: §3.2 · p. 9, Figure 7 · p. 9
  • In the oppositely traveling breather pair, the profiles temporarily distorted during collision, then recovered their initial shapes and sizes without observed radiation.The β = 3 breather traveled right and the β = 0.5 breather left.Source: §3.2 · pp. 9–10, Figure 8 · p. 10

3.3. Triple Breathers Interaction.

  • Three left-traveling breathers with different speeds collided and distorted, then separated with their original profiles and no observed radiation.The two faster breathers separated from the slowest before all three were apart by t = 1.6.Source: §3.3 · p. 10, Figure 9 · p. 11

3.4. Triple Breathers-Soliton Interaction.

  • In the triple-breather–soliton collision, all four structures separated with their original profiles, and no radiation or dispersion was observed.The soliton emerged to the right of the collision and the fastest breather to the left.Source: §3.4 · p. 11

4. Perturbations of mKdV breathers

This section uses numerical simulations to examine how amplitude perturbations affect mKdV breathers, tracking their parameters and separating radiation from localized structures. It characterizes outcomes ranging from radiation and approach to a breather profile to loss of breather structure and formation of opposite-sign solitons.

4.1. The computational approach.

  • The computational approach tracks β and γ from mass and energy, then solves for α and δ; localized mass and energy better identify the final breather while quantifying shed radiation.Localization separates the breather from faster left-moving radiation; global parameter tracking can obscure parameter changes after initial jumps.Source: §4.1 · pp. 12–13, §4.2 · p. 14
  • For A = 0.9, the breather sheds radiation leftward and stabilizes; localized tracking gives α′ = 1.4703 and β′ = 1.6032, and the dispersed mass reaches 0.4030, about 1.04% of initial mass.At t = 4.97, the evolved profile matches an exact breather using those tracked parameters; its L∞ peaks oscillate between 5.5 and 8.Source: §4.2 · pp. 14–15, Figure 13 · p. 15
  • Across A < 1 cases, amplitude reduction increases internal oscillations, with initial α below α(0) and β above β(0); simulations below A = 0.6 may require β > 4 to distinguish breathers from radiation.The authors report A = 0.6 with β = 4 as an approximate threshold for numerical distinction.Source: §4.2 · p. 15
  • For A = 1.1, perturbed breathers shed radiation and match nearby asymptotic breathers, with β determining whether they travel left, remain near the origin, or travel right.The tracked α′ decreases and β′ increases as initial β increases; the profile matches support asymptotic stability in these cases.Source: §4.3 · pp. 16, 17, 19, Figure 16 · p. 17
  • For A = 1.2, β = 1.15 and 1.16 evolve into breather profiles, whereas β = 1.3 breaks into positive and negative solitons with radiation.The β = 1.15 profile completes a period by t = 8.01; the β = 1.16 profile is matched at its half-period around t = 9.52.Source: Figure 19 · p. 19, §4.3 · p. 20
  • For A = 1.2, initial β below β* leads to radiation shedding and convergence to an asymptotic breather, whereas β above β* yields a positive and a negative soliton.At β = 1.3, periodic breathing disappears; convergence to a double pole at the numerically observed threshold is expected but not rigorously proved.Source: §4.3 · pp. 21–22
  • For A = 1.3, the numerically estimated critical α* is approximately 1.08: below it the breather breaks into two solitons, while above it the breather stabilizes.The α = 1.08 case is near the boundary and has a slow period; α = 1.12 stabilizes with faster periodicity.Source: Figure 23 · p. 22, §4.3 · pp. 23–24
  • As initial α increases for A = 1.3, the initial mass stays fixed at ∥u0∥2 L2 = 40.56, while local mass and dispersion mass both decrease.The passage reports these mass trends across the tested initial α values.Source: §4.3 · p. 24
  • For larger amplitude perturbations, the approximate critical α* increases while β* decreases, with higher α or lower β favoring breather stability.The authors attribute this to profile stretching and increased height; the critical-value estimates are accurate only to order O(h^2).Source: §4.3 · pp. 24–25
  • At the numerically observed critical threshold, the breather sheds radiation and is suggested to converge to a double pole solution, but this behavior is not proved.The authors relate this suggestion to the known result that a breather converges to a double pole solution as α converges to 0.Source: §4.3 · pp. 22, 25

5. Perturbations of the nonlinear term

The section studies mKdV breathers under combined nonlinearities, including an integrable Gardner-type case and a non-integrable case with a symmetric potential. Numerical profile matching and local-mass tracking support stable breather-like behavior in both cases over the simulated intervals, alongside radiative shedding.

5. Perturbations of the nonlinear term

  • With p=2, the section compares a Gardner-variant integrable perturbation with a non-integrable absolute-value perturbation; mass and energy remain conserved in the latter.The absolute-value potential is symmetric, and the paper reports positive numerical evidence for the associated breather-existence conjecture.Source: §5.1 · pp. 25–26, §5.1.2 · p. 27

5.1. The mKdV with combined nonlinearities.

  • In the Gardner-type case, small leftward radiation precedes a negative breather profile by t=2.55, while matching profiles at t=100 and t=300 supports periodicity and stability.The initial condition is an mKdV breather; the study does not track Gardner breather parameters.Source: §5.1.1 · p. 26, Figure 27 · p. 26
  • In the symmetric-potential, non-integrable case, profile matching and local-mass stabilization around t=200 support a stable breather-like state, with radiative mass about 0.0037.The matched profiles are from t=200.4 and t=399.9; the passage reports the radiative mass for this simulation.Source: §5.1.2 · pp. 27–28, Figure 29 · p. 28, Figure 30 · p. 29
  • Local-mass tracking captures perturbation changes that the consistent L∞ norm misses, and the local mass stabilizes around t=200 as radiation plateaus.The L∞ norm is consistent throughout the simulation but does not capture the small changes caused by perturbing N(u).Source: §5.1.2 · p. 28
  • For the non-integrable simulation, the final energy error is 1.3417 × 10^-9; increasing spatial resolution improves mass and energy errors but increases computational time.The energy error is reported for the domain size L = 500π used in Figures 29 and 30.Source: §5.1.2 · p. 28, Figure 30 · p. 29

6. Conclusions

The conclusions describe a numerical study of mKdV breather structure, perturbations, and interactions, using localized-structure and mass diagnostics to characterize long-time behavior. The simulations provide evidence that breather-like dynamics can persist in non-integrable mKdV-type models, while motivating further analytical work on stability and long-time asymptotics.

6. Conclusions

  • Local-mass tracking diagnoses convergence and persistence of breather-like structures, while estimating left-propagating radiation as a fraction of the initial breather mass.The study also applies these diagnostics to multi-breather and breather-soliton interactions.Source: §6 · p. 29
  • The simulations provide numerical evidence that breather behavior can persist in non-integrable mKdV-type models, with both coherent breather components and accompanying radiation identified and quantified.The authors suggest adapting these diagnostics to related equations in the sine-Gordon/mKdV hierarchy.Source: §6 · p. 29
  • The numerical findings motivate future analytical work to rigorously justify the observed stability and long-time asymptotics.Source: §6 · p. 30
Loading 2609.27080v1…