[论文解读] Number of solitons produced from a large initial pulse in the generalized NLS dispersive hydrodynamics theory

L. F. Calazans de Brito, A. M. Kamchatnov|arXiv (Cornell University)|Aug 19, 2021

Nonlinear Waves and Solitons参考文献 31被引用 8

一句话总结

本文推导出在广义非线性薛定谔方程（gNLS）中，即使该方程并非完全可积，大初始脉冲激发的孤子数量的解析公式。通过基于色散冲击波中波峰演化特性的直接方法，作者获得了一个通用公式 N = (1/2π)∫k(c(x))dx，其中 k(c) 是由局域声速 c(x) 导出的波数，该公式推广了对可积系统（如标准 NLS 方程）的先前结果。

ABSTRACT

We show that the number of solitons produced from an arbitrary initial pulse of the simple wave type can be calculated analytically if its evolution is governed by a generalized nonlinear Schr\"{o}dinger equation provided this number is large enough. The final result generalizes the asymptotic formula derived for completely integrable nonlinear wave equations like the standard NLS equation with the use of the inverse scattering transform method.

研究动机与目标

• 推导广义 NLS 方程中由初始脉冲产生的孤子数量的一般解析公式。
• 将孤子计数方法扩展至完全可积系统之外，如标准 NLS 方程。
• 通过与不同非线性度（p ≠ 1）的数值解比较，验证该方法。
• 在不依赖于光滑区域中波数分布的非物理假设的前提下，提供孤子计数公式的物理解释推导。

提出的方法

• 应用色散冲击波（DSW）理论中的波峰计数方法，追踪 DSW 区域内的波峰。
• 利用 DSW 小振幅边缘处波峰数变化率 dN/dt = (1/2π)k(vg - vph)。
• 依赖 El 的方法来确定沿 DSW 边缘的波数 k(u) 和时间演化 t(u)。
• 将时间积分的波峰数转换为关于初始局域声速 c(x) 的空间积分。
• 通过将 k 表示为初始声速分布和非线性参数 p 的函数，推导出最终公式 N = (1/2π)∫k(c(x))dx。
• 通过不同 p 值（p = 1/2, 1, 2）的数值模拟验证该公式，且初始脉冲形状各异。

实验结果

研究问题

• RQ1 当广义 NLS 方程并非完全可积时，是否能从大初始脉冲中解析预测产生的孤子数量？
• RQ2 色散冲击波中的波峰计数方法是否能导出一个适用于孤子计数的通用公式，并在可积极限下与已知结果一致？
• RQ3 在 gNLS 方程中，孤子数量如何依赖于初始脉冲形状和非线性强度（p）？
• RQ4 该推导公式在 gNLS 方程的不同物理参数区域（包括非可积情形）中是否具有鲁棒性和准确性？
• RQ5 对初始平滑状态中波数分布的物理假设，能否通过直接推导加以证明？

主要发现

• 孤子数量由 N = (1/2π)∫k(c(x))dx 给出，其中 k(c) 依赖于非线性参数 p 和初始声速 c(x)。
• 当 p = 1（标准 NLS）时，公式简化为 N = (2/π)∫√[c0(c0 - c(x))]dx，与亮孤子的已知结果一致。
• 该公式在 p = 1/2, 1, 2 时通过数值模拟得到验证，与 gNLS 方程的直接模拟结果高度吻合。
• 推导过程避免了对光滑区域中波数分布的非物理假设，使该方法比以往方法更具物理透明性。
• 该方法为实验中非可积非线性波系统（如超冷费米气体或玻色-爱因斯坦凝聚体）中的孤子计数提供了可靠的预测工具。
• 该方法证实了孤子形成过程在可积与非可积非线性波方程中的普适性。