[论文解读] Numerical verification of a gap condition for linearized NLS

Laurent Demanet, Wilhelm Schlag|ArXiv.org|Aug 13, 2005

Advanced Mathematical Physics Problems参考文献 16被引用 7

一句话总结

该论文通过Birman-Schwinger方法和针对1/|x−y|核的谱精度求积规则，数值验证了三维线性化三次非线性薛定谔方程的谱隙条件。结果表明，当$\beta > \beta_* \approx 0.913958905$时，算子$L_+$和$L_-$在间隙$(0, \alpha^2]$内均无本征值或共振态，从而确立了$\beta > \beta_*$时的间隙性质，这对于孤立波的条件渐近稳定性至关重要。

ABSTRACT

We make a detailed numerical study of the spectrum of two Schroedinger operators L_- and L_+ arising in the linearization of the supercritical nonlinear Schroedinger equation (NLS) about the standing wave, in three dimensions. This study was motivated by a recent result of the second author on conditional asymptotic stability of solitary waves in the case of a cubic nonlinearity. Underlying the validity of this result is a spectral condition on the operators L_- and L_+, namely that they have no eigenvalues nor resonances in the gap (a region of the positive real axis between zero and the continuous spectrum,) which we call the gap property. The present numerical study verifies this spectral condition, and further shows that the gap property holds for NLS exponents of the form 2*beta + 1, as long as beta_* < beta <= 1, where beta_* = 0.913958905 +- 1e-8. Our strategy consists of rewriting the original eigenvalue problem via the Birman-Schwinger method. From a numerical analysis viewpoint, our main contribution is an efficient quadrature rule for the kernel 1/|x-y| in R^3, i.e., provably spectrally accurate. As a result, we are able to give similar accuracy estimates for all our eigenvalue computations. We also propose an improvement of the Petviashvili's iteration for the computation of standing wave profiles which automatically chooses the radial solution.

研究动机与目标

• 验证三维线性化NLS的谱隙条件，该条件是孤立波条件渐近稳定性的必要前提。
• 确定临界值$\beta_*$，即当$\beta < \beta_*$时，$L_+$在间隙$(0, \alpha^2]$内获得本征值，导致间隙性质失效。
• 开发一种数值上稳健且可复现的方法，用于计算由NLS线性化产生的Schrödinger算子的本征值。
• 改进Petviashvili迭代方法，以增强其收敛性并保持径向对称性，用于计算驻波解形貌。
• 提供一个公开可获取、可复现的数值框架实现。

提出的方法

• 利用Birman-Schwinger方法重写$L_\pm$的本征值问题，将其转化为紧算子的本征值问题。
• 实现一种可证明具有谱精度的求积规则，用于$\mathbb{R}^3$中的$1/|x-y|$核，确保高精度的本征值计算。
• 采用径向谱方法结合傅里叶变换，通过Petviashvili迭代计算孤立子形貌$\phi$，并自动保证径向对称性。
• 通过Petviashvili常数$M_n$和Euler-Lagrange残差$\| -\Delta\phi_n + \phi_n - \phi_n^{2\beta+1} \|_{L^2}$监控收敛性。
• 计算Birman-Schwinger算子$K_\pm$的最大本征值，以检测间隙中的本征值：当$\beta_*$处$\lambda_2(K_-) < 1$且$\lambda_5(K_+)$穿过1时，表明$K_+$在$\beta_*$处出现新本征值。
• 采用高分辨率计算，设置$L=25$，$N=200$，并使用三次插值，将$\beta_*$估计至8位有效数字精度。

实验结果

研究问题

• RQ1 当$\beta \leq 1$时，线性化NLS算子$L_+$在间隙$(0, \alpha^2]$内是否存在本征值？
• RQ2 临界值$\beta_*$为何值，使得当$\beta < \beta_*$时间隙条件失效？
• RQ3 能否在$\mathbb{R}^3$中为$1/|x-y|$核构造一种谱精度求积规则，以实现高精度本征值计算？
• RQ4 如何改进Petviashvili迭代方法，以确保径向对称性并收敛至基态？
• RQ5 用于计算间隙条件的数值框架是否可复现且公开可用？

主要发现

• 当$\beta > \beta_* \approx 0.913958905 \pm 10^{-8}$时，间隙条件成立，确认了超临界NLS中孤立波条件渐近稳定性的谱要求。
• Birman-Schwinger算子$K_+$的第五个本征值$\lambda_5(K_+)$在$\beta_*$处穿过1，标志$L_+$的间隙中开始出现离散本征值。
• 该数值方法对$\mathbb{R}^3$中的$1/|x-y|$核实现了谱精度，支持可靠本征值计算并可提供误差估计。
• 改进后的Petviashvili迭代收敛至径向基态，精度达机器精度，且径向对称性保持在16位有效数字精度。
• 数值实验完全可复现，相关MATLAB代码已公开发布于http://www.acm.caltech.edu/~demanet/NLS/。
• 本研究证实，$L_-$在$(0, \alpha^2]$内无本征值或共振态，而$L_+$仅在$\beta < \beta_*$时在间隙中存在一个本征值。