Law of Large Numbers and Central Limit Theorem for random sets of solitons of the focusing nonlinear Schrödinger equation
Girotti, Grava, McLaughlin et al.
We study a random configuration of $N$ soliton solutions $Ï_N(x,t;\boldsymbolλ)$ of the cubic focusing Nonlinear Schrödinger (fNLS) equation in one space dimension. The $N$ soliton solutions are parametrized by $2N$ complex numbers $(\boldsymbolλ, \boldsymbol{c})$ where $\boldsymbolλ\in\mathbb{C}_+^N$ are the eigenvalues of the Zakharov-Shabat linear operator, and $ \boldsymbol{c}\in\mathbb{C}^N\backslash \{0\}$ are the norming constants of the corresponding eigenfunctions. The randomness is obtained by choosing the complex eigenvalues to be i.i.d. random variables sampled from a probability distribution with compact support in the complex plane. The corresponding norming constants are interpolated by a smooth function of the eigenvalues. Then we consider the expectation of the random measure associated to this random spectral data. Such expectation uniquely identifies, via the Zakharov-Shabat inverse spectral problem, a solution $Ï_\infty(x,t)$ of the fNLS equation. This solution can be interpreted as a soliton gas solution.
We prove a Law of Large Numbers and a Central Limit Theorem for the differences $Ï_N(x,t;\boldsymbolλ)-Ï_\infty(x,t)$ and $|Ï_N(x,t;\boldsymbolλ)|^2-|Ï_\infty(x,t)|^2$ when $(x,t)$ are in a compact set of $\mathbb R\times\mathbb R^+$; we additionally compute the correlation functions.
academic
Law of Large Numbers and Central Limit Theorem for random sets of solitons of the focusing nonlinear Schrödinger equation
This paper investigates random configurations ψN(x,t;λ) of N-soliton solutions to the one-dimensional cubic focusing nonlinear Schrödinger equation (fNLS). The N-soliton solutions are parametrized by 2N complex numbers (λ,c), where λ∈C+N are the eigenvalues of the Zakharov-Shabat linear operator and c∈CN∖{0} are normalization constants of the corresponding eigenfunctions. Randomness is introduced by choosing complex eigenvalues as independent identically distributed random variables sampled from a compactly supported probability distribution on the complex plane. The corresponding normalization constants are interpolated through smooth functions of the eigenvalues. The authors then consider the expectation value of random measures associated with this random spectral data, which uniquely determines through the Zakharov-Shabat inverse spectral problem a solution ψ∞(x,t) to the fNLS equation, interpretable as a soliton gas solution. The authors prove the law of large numbers and central limit theorem for the differences ψN(x,t;λ)−ψ∞(x,t) and ∣ψN(x,t;λ)∣2−∣ψ∞(x,t)∣2 when (x,t) belongs to compact subsets of R×R+, and compute the associated correlation functions.
Core Problem to be Addressed: The paper addresses the establishment of a rigorous statistical theory for random N-soliton solutions. Specifically, it seeks to understand the relationship between random N-soliton solutions ψN(x,t;λ) and the deterministic limiting solution ψ∞(x,t) as the number of solitons N approaches infinity.
Importance of the Problem:
Probability distributions of nonlinear waves undergo substantial deformation during time evolution, in stark contrast to the linear case
For integrable nonlinear partial differential equations, a predictive statistical theory for large-amplitude waves must be established
Soliton dynamics theory and generalized hydrodynamic equations lack rigorous mathematical derivations
Limitations of Existing Approaches:
Existing soliton dynamics theories are primarily qualitative
Lack of methods to rigorously derive dynamical equations through analysis of solutions to underlying nonlinear PDEs
Rigorous analysis of random N-soliton solutions has not been established
Research Motivation: By introducing randomness in the linear setting of scattering data, the paper studies random large-amplitude nonlinear dynamics, providing a predictive statistical theory for large-amplitude waves in the fNLS equation.
Established Law of Large Numbers for Random N-Soliton Solutions: Proved that as N→∞, ψN(x,t;λ) and ∣ψN(x,t;λ)∣2 converge in mean to the deterministic solutions ψ∞(x,t) and ∣ψ∞(x,t)∣2
Proved Central Limit Theorem: Established that N(ψN(x,t;λ)−ψ∞(x,t)) and N(∣ψN(x,t;λ)∣2−∣ψ∞(x,t)∣2) converge to Gaussian random variables
Computed Correlation Functions: Provided explicit expressions for two-point correlation functions
Provided Rigorous Mathematical Framework: First established rigorous analytical methods for random N-soliton solutions to nonlinear dispersive PDEs
Random N-Soliton RH Problem: Find a 2×2 matrix-valued function MN(z;x,t,λ) satisfying:
Analyticity in C∖γ
Boundary conditions: (MN)+(z)=(MN)−(z)JN(z;x,t,λ), where
\begin{pmatrix} 1 & 0 \\ -e^{\theta(z;x,t)}\sum_{k=1}^N \frac{c_k}{z-\lambda_k} & 1 \end{pmatrix}, & z \in \gamma_+ \\
\begin{pmatrix} 1 & e^{-\theta(z;x,t)}\sum_{k=1}^N \frac{c_k}{z-\lambda_k} \\ 0 & 1 \end{pmatrix}, & z \in \gamma_-
\end{cases}$$
**Average RH Problem**: Corresponding deterministic problem with jump matrix $J(z;x,t) = \mathbb{E}[J_N(z;x,t,\lambda)]$
#### 2. Error Analysis Framework
Define error matrix $E(z) = M_N(z)M(z)^{-1}$ with jump matrix:
$$J_E(z) = M_-(z)J_N(z)J(z)^{-1}M_-(z)^{-1}$$
Introduce linear statistics:
$$X_N^f(z) = \sum_{k=1}^N \frac{r(\lambda_k)}{z-\lambda_k} - N\iint_{D_+} \frac{r(w)}{z-w}d\mu(w)$$
### Technical Innovations
1. **Small Norm RH Theory with High Probability**:
- Introduce set $B_\delta = \{|\lambda_1,\ldots,\lambda_N| : \sup_{z\in\gamma_+} |X_N^f(z)|/N < \delta\}$
- Prove $P((B_\delta)^c) \leq \frac{c_1}{\delta^{2p}N^{p}} + \frac{c_2}{\delta^{2p+1}N^{p+1/2}}$
2. **Neumann Series Expansion**: When configurations lie in $B_\delta$ and $c_0\delta < 1$, the error matrix can be expressed as a convergent Neumann series
3. **Probabilistic Estimation Techniques**: Combine central limit theorem with uniform estimates to control the behavior of linear statistics
## Experimental Setup
### Theoretical Verification Framework
This is primarily a theoretical work, with results verified through rigorous mathematical proofs rather than numerical experiments. Main verifications include:
1. **Existence Proofs**: Use Zhou's vanishing lemma to prove existence and uniqueness of solutions to the average RH problem
2. **Uniform Estimates**: Establish uniform estimates through boundedness of Cauchy operators
3. **Probabilistic Estimates**: Use Markov inequality and gridding techniques to control probabilities
### Key Estimates
- **Upper Bound for N-Soliton Solutions**: $|\psi_N(x,t;\lambda)| \leq 4\sum_{k=1}^N \text{Im}(\lambda_k)$
- **Operator Norm Estimates**: $\|C_{W_N}\| \leq c_0\delta$
- **Probability Decay Rate**: $P((B_\delta)^c) = O(N^{-p})$ for appropriately chosen $p$
## Main Results
### Law of Large Numbers (Theorem 2.6)
For compact sets $(x,t) \in K \subset \mathbb{R}\times\mathbb{R}_+$:
$$\lim_{N\to\infty} \mathbb{E}[|\psi_N(x,t;\lambda) - \psi_\infty(x,t)|] = 0$$
$$\lim_{N\to\infty} \mathbb{E}[||\psi_N(x,t;\lambda)|^2 - |\psi_\infty(x,t)|^2|] = 0$$
### Central Limit Theorem (Theorem 2.7)
The random variables $\sqrt{N}(\psi_N(x,t;\lambda) - \psi_\infty(x,t))$ and $\sqrt{N}(|\psi_N(x,t;\lambda)|^2 - |\psi_\infty(x,t)|^2)$ converge respectively to complex and real Gaussian random variables $X_{G_1}$ and $X_{G_2}$, with covariances:
$$\mathbb{E}[X_{G_i}(x,t)^2] = \iint_{D_+} G_i(w;x,t)^2 d\mu(w) - \left(\iint_{D_+} G_i(w;x,t)d\mu(w)\right)^2$$
where functions $G_1$ and $G_2$ are explicitly given by solutions to the average RH problem.
### Correlation Functions (Theorem 2.8)
The limit of two-point correlation functions is:
$$\lim_{N\to\infty} R_{2,N}(x_1,t_1,x_2,t_2) = \iint_{D_+} G_1(s;x_1,t_1)G_1(s;x_2,t_2)d\mu(s) - \iint_{D_+} G_1(s;x_1,t_1)d\mu(s)\iint_{D_+} G_1(s';x_2,t_2)d\mu(s')$$
## Related Work
1. **Integrable Systems Theory**: Built upon Zakharov-Shabat inverse scattering theory and RH problem methods
2. **Soliton Dynamics Theory**: Extends soliton dynamics theory proposed by Zakharov, providing rigorous mathematical foundations
3. **Random Integrable Systems**: Related to work by Deift-Trogdon on random matrix spectral data
4. **Generalized Hydrodynamics**: Provides rigorous derivation framework for GHD equations
## Conclusions and Discussion
### Main Conclusions
1. First established rigorous statistical theory for random N-soliton solutions
2. Proved universal behavior of random soliton configurations in the large N limit
3. Provided rigorous mathematical foundations for soliton gas theory
### Limitations
1. Results restricted to compact sets $(x,t) \in K \subset \mathbb{R}\times\mathbb{R}_+$
2. Requires eigenvalue distributions to have compact support
3. Boundary behavior of initial profiles remains an open problem
### Future Directions
1. Extension to other integrable equations (KdV, sine-Gordon, etc.)
2. Investigation of non-compactly supported distributions
3. Establishment of connections with experimental observations
## In-Depth Evaluation
### Strengths
1. **Theoretical Rigor**: Provides complete and rigorous mathematical proof framework
2. **Methodological Innovation**: Cleverly combines RH problem theory, probability theory, and asymptotic analysis
3. **Depth of Results**: First establishment of CLT for random soliton systems, filling theoretical gaps
4. **Technical Sophistication**: Combination of small norm RH theory with probabilistic methods is highly skillful
### Weaknesses
1. **Limited Application Scope**: Results apply only to specific parameter regions and spacetime domains
2. **Computational Complexity**: Explicit computation of correlation functions may be difficult in practical applications
3. **Absence of Numerical Verification**: Lacks numerical experiments to validate theoretical predictions
### Impact
1. **Academic Value**: Establishes foundations for random theory of integrable systems
2. **Methodological Contribution**: Method combining RH problems with probability theory has broad applicability
3. **Application Prospects**: Provides new theoretical tools for understanding nonlinear wave phenomena
### Applicable Scenarios
1. Soliton propagation in nonlinear optics
2. Wave phenomena in fluid mechanics
3. Integrable models in condensed matter physics
4. Research on random integrable systems in mathematical physics
## References
The paper contains 53 references covering important works in integrable systems theory, RH problem methods, probability theory, and related physical applications across multiple fields.
