Hirota Bilinear Method for Soliton Solutions of the Qiao-Li Equation

Technical Summary: Hirota Bilinear Method for Soliton Solutions of the Qiao-Li Equation

Problem Statement
The Qiao-Li equation, a new integrable nonlinear wave equation proposed by Professors Zhijun Qiao and Shengtai Li, extends the integrable hierarchy of shallow water wave models (such as the Camassa-Holm and Degasperis-Procesi equations). While previous research has established the integrability of the Qiao-Li equation and investigated its parametric and stationary solutions, a systematic exploration of its multi-soliton solutions and interaction dynamics has been lacking. Specifically, explicit expressions for single, double, triple, and general $N$-soliton solutions, along with a detailed analysis of their collision properties, remained unexplored. This paper addresses this gap by applying the Hirota bilinear method to construct these exact solutions and analyze their nonlinear wave dynamics.

Methodology
The study employs the Hirota bilinear method, an algebraic approach for constructing exact solutions to nonlinear integrable equations. The methodology proceeds through the following steps:

1. Variable Transformation: To handle the complex terms involving fractional powers ($u^{-1/3}$ and $u^{-1/6}$) in the original Qiao-Li equation, a power transformation $u = v^{-6}$ is introduced. This converts the original equation into a more tractable form involving the new dependent variable $v(x,t)$.
2. Bilinearization: A further transformation is applied using a tau function $F(x,t)$, defined by $v(x,t) = [2\partial_x^2 \ln F(x,t)]^{1/7}$. This reduces the nonlinear partial differential equation to a bilinear form:
   $$(D_x^6 + D_x D_t) F \cdot F = 0$$
   where $D$ represents the Hirota bilinear derivative operator.
3. Solution Construction:
   • Single-Soliton: A tau function of the form $F = 1 + e^\eta$ is substituted to derive the dispersion relation $\omega = k^5$ and the explicit single-soliton solution.
   • Multi-Soliton: The tau function is extended to include multiple exponential terms and interaction coefficients ($A_{ij}$). For the two-soliton case, the coefficient $A_{12}$ is determined by ensuring the bilinear equation is satisfied, yielding $A_{12} = (k_1 - k_2)^4 / (k_1 + k_2)^4$. This logic is extended to derive the three-soliton solution and the general $N$-soliton solution.
4. Verification: The validity of the general $N$-soliton solution is confirmed through mathematical induction, demonstrating that the constructed tau functions satisfy the bilinear equation for any positive integer $N$.
5. Numerical Analysis: High-precision numerical simulations are conducted to visualize the solutions, specifically analyzing the propagation, amplitude, width, and collision dynamics of single, double, and triple solitons.

Key Contributions and Results

• Explicit Exact Solutions: The paper successfully derives closed-form expressions for the single, double, triple, and general $N$-soliton solutions of the Qiao-Li equation. The general $N$-soliton solution is presented in a compact tau function form involving a summation over binary combinations of phase variables and interaction coefficients.
• Dispersion and Interaction Laws: The study establishes the dispersion relation $\omega = k^5$ and derives the specific interaction coefficients $A_{ij} = (k_i - k_j)^4 / (k_i + k_j)^4$, which govern the nonlinear interaction between solitons.
• Collision Dynamics: Numerical simulations reveal that the Qiao-Li equation exhibits elastic collisions. During interactions (both double and triple soliton collisions), solitons undergo transient amplitude enhancements and phase shifts but fully recover their original shapes, amplitudes, and velocities after separation.
• Waveform Characteristics:
  • Amplitude Invariance: The maximum amplitude of the single soliton is found to be constant ($1/64$) regardless of the wavenumber $k$.
  • Velocity and Width: The propagation velocity follows a strong nonlinear dependence $v = k^5$. The soliton width is inversely proportional to $k$; as $k$ increases, the soliton becomes narrower and more localized.
  • Peakon-like Structure: The solutions exhibit sharp peaks and steep flanks resembling "peakons" (characteristic of the Camassa-Holm equation). However, the paper notes that these peaks remain mathematically smooth and differentiable, lacking the true cusps found in standard peakon solutions.
  • Valley-Peak Transformation: The analysis of two-soliton collisions characterizes a specific nonlinear phenomenon involving the mutual conversion between antisolitons and solitons, described as a "valley peak transformation."

Significance
The paper claims to fill a critical research gap by providing the first systematic construction of multi-soliton solutions for the Qiao-Li equation using the Hirota method. By confirming the equation's rich soliton structure and its membership in the family of integrable systems, the study validates the equation's theoretical robustness.

The work distinguishes the Qiao-Li equation from other integrable models (like the Camassa-Holm equation) by revealing its unique wave dynamics, specifically the smooth yet peakon-like nature of its solutions and the specific mechanics of its soliton collisions. The authors assert that these findings provide a theoretical foundation for understanding complex nonlinear wave phenomena, with potential relevance to fields such as optical fiber communications, fluid dynamics, and plasma physics, although the paper primarily focuses on the mathematical derivation and theoretical characterization of these solutions.