Radiation damping of the soliton internal mode in 1D quadratic Klein-Gordon equation

This paper demonstrates that on a codimension-one manifold of fine-tuned initial data, the internal mode of a soliton in the 1D quadratic Klein-Gordon equation undergoes slow, irreversible decay into dispersive radiation due to radiation damping, a process accurately described by a cubic resonant approximation with a damping rate governed by a Fermi golden rule-type coefficient.

The Gist

Imagine you have a perfectly balanced, solitary wave—a "soliton"—traveling through a medium. Think of it like a perfect, self-contained wave in a pond that doesn't spread out or fade away on its own. In the world of physics, this is a stable, localized packet of energy.

However, in the specific equation studied in this paper (the 1D quadratic Klein-Gordon equation), this soliton has a secret: it's not just a calm, static object. It has a "heartbeat."

The Soliton's Heartbeat and Its Wobbly Legs

The researchers discovered that this soliton has two distinct internal behaviors:

1. The Unstable Leg (The Wobble): The soliton has a tendency to tip over. If you nudge it slightly, it wants to collapse or explode. It's like a pencil balanced on its tip; it can stay there, but the slightest push sends it falling.
2. The Internal Mode (The Heartbeat): Inside the soliton, there is a localized vibration, a rhythmic oscillation. Think of it like a tiny bell ringing inside the wave.

The problem is that in nature, things rarely stay perfectly balanced. If you try to set this soliton in motion, the "Unstable Leg" usually wins, and the whole thing falls apart. But, the researchers asked: What if we could perfectly balance the pencil so it doesn't fall, but just keeps ringing?

The "Fine-Tuned" Balancing Act

To stop the soliton from collapsing, the researchers had to be incredibly precise. They had to set up the initial conditions (the starting push) with perfect accuracy. They call this a "codimension-one manifold," which is a fancy way of saying: "There is a very specific, thin slice of starting conditions where the soliton doesn't explode, but instead survives."

Once they found this perfect balance, they watched what happened to the "Heartbeat" (the internal mode).

The Leaky Antenna

Here is the core discovery: Even though the soliton is stable enough not to collapse, the heartbeat doesn't last forever. It slowly fades away.

Why? Because the soliton acts like a leaky antenna.

• The Metaphor: Imagine the soliton is a radio tower (the core) with a vibrating antenna (the internal mode). The antenna is vibrating at a specific frequency.
• The Interaction: Because the system is non-linear (the waves interact with themselves), the vibration of the antenna doesn't just stay inside the tower. It starts to "talk" to the surrounding ocean of waves (the continuum).
• The Result: The internal vibration transfers its energy into the surrounding ocean, creating ripples that travel away. This is called radiation damping.

The paper explains exactly how this energy leaks out. It's not a sudden explosion; it's a slow, steady drain. The amplitude of the heartbeat gets smaller and smaller over time, following a very specific mathematical rule (it decays like $1/\sqrt{t}$).

The "Fermi Golden Rule" and the Magic Number

The researchers used advanced math (called "normal form methods") to simplify this complex interaction. They found that the rate at which the energy leaks out is determined by a specific coefficient, which they call the Fermi Golden Rule.

Think of this rule as a "leakage coefficient." It's a number that tells you exactly how efficient the soliton is at turning its internal vibration into outgoing ripples.

• If the number is high, the soliton loses energy fast.
• If the number is low, it holds onto its energy longer.

In this specific equation, they calculated this number precisely. They also found that as the energy leaks out, the "pitch" of the heartbeat changes slightly (a frequency shift), just like a spinning ice skater changes speed when they pull their arms in.

Why Does This Matter?

You might ask, "Who cares about a mathematical wave in a 1D equation?"

This mechanism is everywhere in physics:

• Optical Fibers: Light pulses traveling through fiber optic cables can have internal modes that leak energy, affecting how far the signal travels.
• Bose-Einstein Condensates: These are super-cold clouds of atoms that act like a single giant wave. They have internal vibrations that can decay.
• Cosmology: Some theories suggest the universe was once filled with "solitons" (like cosmic strings or domain walls). Understanding how they lose energy helps us understand the evolution of the universe.

The Bottom Line

The paper is a story about metastability. It shows that even when you perfectly balance a system to prevent it from collapsing, it still has a slow, inevitable decay. The internal energy of the soliton is slowly siphoned off into the surrounding universe as radiation.

The authors didn't just guess this; they derived the exact formula for the decay rate and the frequency shift, and then they built a computer simulation to prove it. The computer numbers matched their math perfectly, confirming that the "leaky antenna" model is the correct way to understand how these solitary waves relax over time.

In short: They found the exact recipe for how a stable, vibrating wave slowly "sings" its energy away into the void, turning a complex, chaotic problem into a simple, predictable rhythm.

Technical Summary

1. Problem Statement

The paper investigates the long-time dynamics of small, even perturbations of a soliton solution in the one-dimensional quadratic nonlinear Klein-Gordon (NLKG) equation:
$$\phi_{tt} - \phi_{xx} + \phi = \phi^2$$
The system possesses a static soliton solution $S(x) = \frac{3}{2} \text{sech}^2(x/2)$. Linearizing around this soliton reveals a spectral structure containing:

• A continuous spectrum $[1, \infty)$ representing dispersive radiation.
• An unstable mode (eigenvalue $-\lambda^2 = -5/4$) leading to exponential growth.
• A zero mode (eigenvalue 0) associated with translation invariance (odd parity, thus not excited by even initial data).
• An internal mode (eigenvalue $\omega^2 = 3/4$), a localized oscillatory excitation.

The Core Challenge: The coexistence of an internal mode and an unstable mode makes the system sensitive. Generic perturbations lead to finite-time blowup due to the unstable mode. However, on a codimension-one manifold of fine-tuned initial data, the unstable direction is suppressed. The authors aim to characterize the long-time behavior of the internal mode on this stable manifold. Specifically, they seek to quantify how the internal mode decays by transferring energy to the continuous spectrum (radiation damping) and to determine the associated nonlinear frequency shift.

2. Methodology

The authors employ a combination of normal form theory, spectral analysis, and numerical validation.

A. Spectral Decomposition

The perturbation $u(t,x)$ is decomposed into orthogonal components:
$$u(t, x) = a(t)\psi(x) + b(t)\xi(x) + \eta(t, x)$$
where $\psi$ is the internal mode, $\xi$ is the unstable mode, and $\eta$ is the radiation field (projected onto the continuous spectrum). This yields a coupled system of ODEs for the amplitudes $a(t)$ and $b(t)$, and a PDE for $\eta$.

B. Normal Form Reduction

To isolate the dynamics of the internal mode, the authors perform a near-identity transformation to eliminate non-resonant quadratic terms from the equations of motion.

1. Elimination of Quadratic Terms: They transform variables to remove $O(\epsilon^2)$ terms, reducing the system to a cubic order approximation.
2. Handling the Unstable Mode: By restricting to the stable manifold, the unstable mode amplitude $b(t)$ is shown to be $O(\epsilon^2)$ and decays exponentially or remains bounded, allowing it to be treated as a higher-order correction.
3. Radiation Coupling: The radiation field $\eta$ is eliminated via a transformation involving the solution to inhomogeneous linear equations. Crucially, the frequency of the quadratic nonlinearity ($2\omega = \sqrt{3}$) lies within the continuous spectrum ($[1, \infty)$), enabling resonant energy transfer.

C. Derivation of the Effective Equation

The authors derive a cubic resonant approximation (Birkhoff normal form) for the complex amplitude $A(t)$ of the internal mode. The derivation involves calculating specific coefficients arising from the interaction between the internal mode, the unstable mode, and the radiation field.

• The coefficient of the cubic term $A|A|^2$ is split into a real part (damping) and an imaginary part (frequency shift).
• The damping term is identified with a Fermi Golden Rule coefficient, $\Gamma$, which quantifies the resonant coupling strength between the discrete internal mode and the continuum.

3. Key Contributions and Results

A. Analytical Derivation of Decay Rate

The authors derive an effective ordinary differential equation for the internal mode amplitude $A(t)$:
$$\dot{A} = \left( i\gamma - \frac{\Gamma}{2} \right) A|A|^2$$
where:

• $\gamma \approx 0.045938$ is the nonlinear frequency shift coefficient.
• $\Gamma \approx 0.008966$ is the damping rate coefficient.

The value of $\Gamma$ is explicitly calculated using the overlap of the internal mode squared with the continuum eigenfunctions (Jost functions):
$$\Gamma = \frac{1}{2p\omega} |\langle j_+(\cdot, p), \psi^2 \rangle|^2$$
This confirms the Fermi Golden Rule condition ($\Gamma > 0$), ensuring irreversible energy transfer.

B. Asymptotic Behavior

Solving the effective equation yields the long-time asymptotic behavior for the amplitude $R(t) = |A(t)|$ and phase $\theta(t)$:

• Amplitude Decay: $R(t) \sim \frac{1}{\sqrt{\Gamma t}}$ for $t \gg \epsilon^{-2}$.
  • This indicates a slow, algebraic decay ($t^{-1/2}$) rather than exponential decay.
  • The decay rate becomes independent of the initial amplitude $\epsilon$ at late times (universality).
• Frequency Shift: The phase evolves as $\theta(t) \sim \frac{\gamma}{\Gamma} \ln t$, leading to a logarithmic frequency shift.

C. Radiation Profile

The paper analyzes the "backreaction" of the internal mode on the radiation field. The radiated wave $\eta(t,x)$ is shown to consist of outgoing distorted plane waves with wavenumber $k=\sqrt{2}$ and frequency $2\omega=\sqrt{3}$.

• The radiation amplitude decays as $t^{-1}$ (slower than the standard $t^{-3/2}$ for linear dispersive waves) due to the near-resonance between the driving frequency and the continuum.
• A logarithmic phase modulation is predicted to detune the exact resonance, preventing unbounded growth.

D. Numerical Validation

The authors perform direct numerical simulations of the full NLKG equation with fine-tuned initial data (using a "shooting-to-threshold" method to suppress the unstable mode).

• They track the internal mode amplitude and instantaneous frequency over long times ($t \sim 5000$).
• The numerical fits for $\Gamma$ and $\gamma$ agree with the analytical predictions to within 1%, confirming the accuracy of the cubic resonant approximation.

4. Significance

• Quantitative Description of Metastability: The paper provides a rigorous, quantitative description of how a soliton's internal mode relaxes via radiation damping in the presence of an unstable direction. It moves beyond qualitative stability proofs to provide explicit decay rates.
• Fermi Golden Rule in Soliton Dynamics: It explicitly demonstrates the application of the Fermi Golden Rule to a soliton system with quadratic nonlinearity, linking the microscopic spectral properties to macroscopic dissipation.
• Universality: The result that the late-time decay profile is independent of the initial perturbation size ($\epsilon$) highlights a universal mechanism for energy transfer in nonlinear wave systems.
• Physical Relevance: The mechanism described (internal mode acting as an antenna leaking energy) is relevant to various physical systems, including optical fibers, Bose-Einstein condensates, and field-theoretic solitons, where similar internal modes and radiation continua exist.

In summary, the paper successfully bridges the gap between infinite-dimensional Hamiltonian dynamics and effective finite-dimensional dissipative models, offering a precise mathematical framework for understanding the slow decay of solitons in quadratic nonlinear media.