[논문 리뷰] Numerical verification of a gap condition for linearized NLS

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

Advanced Mathematical Physics Problems참고 문헌 16인용 수 7

한 줄 요약

이 논문은 3차원 선형화된 입자형 비선형 슈뢰딩거 방정식에 대해 스펙트럼 갭 조건을 수치적으로 검증한다. 비르만-슈윙거 방법과 1/|x−y| 커널에 대한 스펙트럼 정확도를 갖춘 구적법을 사용하여, $L_+$와 $L_-$가 갭 $(0, \alpha^2]$ 내에 고유값이나 공명이 없음을 확인한다. 이는 $\beta > \beta_* \approx 0.913958905$일 때 갭 성질이 성립함을 보여주며, 이는 고립파의 조건부 점근적 안정성에 필수적이다.

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.

연구 동기 및 목표

• 3차원 선형화된 NLS에 대한 스펙트럼 갭 조건을 검증하여 고립파의 조건부 점근적 안정성에 필수적인 조건을 확보한다.
• $L_+$가 갭 $(0, \alpha^2]$ 내에 고유값을 갖는 조건을 만족하지 못하는 임계값 $\beta_*$를 특정하여, 갭 성질이 붕괴되는 조건을 규명한다.
• NLS 선형화에서 유도된 슈뢰딩거 연산자의 고유값을 계산하기 위한 수치적으로 안정적이고 재현 가능한 방법을 개발한다.
• 정적 파동 프로파일을 계산하기 위한 페트비아슈빌리 반복의 수렴성과 원형 대칭성을 향상시킨다.
• 수치적 프레임워크를 공개하고 재현 가능한 형태로 제공한다.

제안 방법

• 비르만-슈윙거 방법을 사용하여 $L_\pm$의 고유값 문제를 컴act 연산자 고유값 문제로 재구성한다.
• $\mathbb{R}^3$에서 $1/|x-y|$ 커널에 대해 증명 가능하게 스펙트럼 정확도를 갖춘 구적법을 구현하여 고유값 계산의 고정밀도를 확보한다.
• 푸리에 변환을 활용한 원형 스펙트럼 방법을 사용하여 솔리톤 프로파일 $\phi$를 페트비아슈빌리 반복을 통해 계산하고, 자동으로 원형 대칭성을 유지한다.
• 페트비아슈빌리 상수 $M_n$과 오일러-라그랑주 잔여항 $\| -\Delta\phi_n + \phi_n - \phi_n^{2\beta+1} \|_{L^2}$을 통해 수렴성을 모니터링한다.
• 비르만-슈윙거 연산자 $K_\pm$의 최대 고유값을 계산하여 갭 내 고유값 존재 여부를 탐지한다: $\lambda_2(K_-) < 1$ 이고 $\lambda_5(K_+)$ 가 $\beta_*$에서 1을 초과한다.
• 고해상도 계산을 위해 $L=25$, $N=200$ 및 삼차 보간을 사용하여 $\beta_*$를 8자리 정밀도로 추정한다.

실험 결과

연구 질문

• RQ1 3차원 선형화된 NLS 연산자 $L_+$가 $\beta \leq 1$일 때 갭 $(0, \alpha^2]$ 내에 고유값을 갖는가?
• RQ2 고유값이 갭 내에 존재하게 되는 임계값 $\beta_*$는 무엇인가? 이 값 이하에서는 갭 조건이 실패한다.
• RQ3 $\mathbb{R}^3$에서 $1/|x-y|$ 커널에 대해 스펙트럼 정확도를 갖는 구적법을 구성할 수 있는가? 이를 통해 고정밀도 고유값 계산이 가능해지는가?
• RQ4 페트비아슈빌리 반복을 어떻게 개선하여 원형 대칭성을 유지하고 기저 상태로 수렴하게 할 수 있는가?
• RQ5 스펙트럼 갭 조건을 계산하기 위한 수치적 프레임워크는 재현 가능하며 공개 가능한가?

주요 결과

• $\beta > \beta_* \approx 0.913958905 \pm 10^{-8}$일 때 갭 조건이 성립함을 확인하였으며, 이는 초임계 NLS에서 고립파의 조건부 점근적 안정성에 필요한 스펙트럼 조건을 충족함을 의미한다.
• 비르만-슈윙거 연산자 $K_+$의 다섯 번째 고유값인 $\lambda_5(K_+)$가 $\beta_*$에서 1을 초과함으로써, $L_+$의 갭 내에 이산 고유값이 발생하기 시작하는 신호를 제공한다.
• 이 수치적 방법은 $\mathbb{R}^3$에서 $1/|x-y|$ 커널에 대해 스펙트럼 정확도를 확보하여 오차 추정이 가능한 신뢰할 수 있는 고유값 계산을 가능하게 한다.
• 개선된 페트비아슈빌리 반복은 기계 정밀도 수준에서 원형 기저 상태로 수렴하며, 원형 대칭성은 16자리 정밀도까지 유지된다.
• 수치 실험은 완전히 재현 가능하며, 관련 MATLAB 코드는 http://www.acm.caltech.edu/~demanet/NLS/ 에 공개되어 있다.
• 이 연구는 $L_-$가 $(0, \alpha^2]$ 내에 고유값이나 공명이 없으며, $L_+$는 $\beta < \beta_*$일 때만 갭 내에 단일 고유값을 갖는다를 확인하였다.