Bright Fractional Single and Multi-Solitons in a Prototypical Nonlinear Schr{ö}dinger Paradigm: Existence, Stability and Dynamics

In the present work we explore features of single and pairs of solitary waves in a fractional variant of the nonlinear Schr{ö}dinger equation. Motivated by the recent experimental realization of arbitrary fractional exponents, upon quantifying the tail properties of such coherent structures, we detail their destabilization when the fractional exponent $α$ acquires values $α<1$ and showcase how the relevant destabilization is associated with collapse type phenomena. We then turn to in- and out-of-phase pairs of such waveforms and illustrate how they generically exist for arbitrary $α$ when we cross the harmonic limit, i.e., for $α>2$. Importantly, we use the parameter $α$ as a ``bifurcation parameter’’ in order to connect the harmonic ($α=2$) and biharmonic ($α=4$) limits. Remarkably, not only do we retrieve the instability of all solitonic pairs in the biharmonic case, but showcase a stabilization feature of particular branches of such multipulses that is {\it unique} to the fractional case and does not arise – to our knowledge – for integer multi-pulse settings. We explain systematically this stabilization via spectral analysis and expand upon the implications of our results for the potential observability of fractional multipulse solitary waves.

💡 Research Summary

The manuscript investigates bright solitary waves—both single‑pulse and multi‑pulse configurations—in a fractional‑order nonlinear Schrödinger (FF‑NLS) model. The governing equation is
( i,u_t = D_\alpha u + 2|u|^2 u,\qquad 0<\alpha\le 4, )
where (D_\alpha) denotes the Riesz fractional derivative defined in Fourier space by multiplication with (-|k|^\alpha). When (\alpha=2) the operator reduces to the standard Laplacian, while (\alpha=4) corresponds to a biharmonic (fourth‑order) dispersion. The authors treat the fractional exponent (\alpha) as a continuous bifurcation parameter that interpolates between the harmonic (second‑order) and biharmonic (fourth‑order) limits.

Key Findings

1. Existence and Tail Structure of Single Solitons

   • Numerical continuation of the stationary equation ( D_\alpha\phi + 2\phi^3 - \phi =0) (with (\omega=1)) yields real, localized profiles for every (\alpha) in the interval (0,4].
   • For (0<\alpha\le2) the solutions are strictly positive and possess monotonic, exponentially decaying tails.
   • When (2<\alpha<4) the tails become oscillatory: they can be accurately fitted by an exponentially damped sinusoid (the “near tail”) together with a far‑field algebraic decay (\phi\sim|x|^{-(\alpha+1)}). As (\alpha\to4) the number of zero‑crossings grows without bound, reproducing the infinite‑oscillation pattern known for pure quartic solitons.

2. Linear Stability of the Single Pulse

   • Linearization around the stationary state leads to the eigenvalue problem (\lambda\begin{pmatrix}P\Q\end{pmatrix}= \begin{pmatrix}0&L_-\-L_+&0\end{pmatrix}\begin{pmatrix}P\Q\end{pmatrix}) with (L_\pm = D_\alpha + (2\phi^2-1)) or ((6\phi^2-1)).
   • For (\alpha<1) a pair of real eigenvalues (\pm\lambda) appears, signalling an exponential growth mode and a collapse‑type instability (the critical dimension for the 1‑D fractional NLS is (\alpha=d=1)).
   • At (\alpha=1) the real pair collides at the origin and moves onto the imaginary axis for (\alpha>1). Consequently the single pulse is spectrally stable for all (\alpha\ge1), including the entire super‑harmonic regime (2<\alpha\le4).

3. Two‑Pulse (Multi‑Soliton) Branches

   • Starting from the exact (\alpha=2) sech solution, the authors construct in‑phase and out‑of‑phase two‑pulse states by superposing two well‑separated copies with separation (\delta). Continuation in (\alpha) from 4 down to 2 reveals a finite set of admissible separations for each (\alpha).
   • For (\alpha>2) both in‑phase and out‑of‑phase pairs exist; as (\alpha\to2^+) the separation diverges, indicating a bifurcation from the harmonic limit where no multi‑pulse steady state exists.
   • In the biharmonic limit ((\alpha=4)) all multi‑pulse configurations are unstable, reproducing earlier results. However, in the intermediate fractional regime (2<\alpha<4) the authors discover stability windows for even‑pulse branches (e.g., two‑pulse, four‑pulse). These windows are absent in the integer‑order cases and are attributed to the complex eigenvalues of the linearized operator that remain purely imaginary for certain (\alpha).

4. Nonlinear Dynamics and Collapse

   • Direct time‑integration of the FF‑NLS confirms the linear predictions. For (\alpha<1) the real eigenvalue triggers a rapid blow‑up of the amplitude, a hallmark of fractional collapse.
   • In the super‑harmonic regime, unstable two‑pulse states exhibit two distinct mechanisms: (i) a symmetry‑breaking drift that changes the inter‑pulse distance, leading to either separation or merger, and (ii) a phase‑flip instability for out‑of‑phase pairs.
   • Stable even‑pulse branches display only internal breathing oscillations; perturbations do not destroy the multi‑pulse structure over long propagation distances.

5. Experimental Relevance

   • Recent optical experiments have demonstrated the ability to engineer arbitrary fractional dispersion by shaping the Fourier‑space response of a waveguide or using programmable spatial light modulators. The paper cites works where (\alpha) has been tuned continuously between 1 and 4.
   • The authors argue that the newly identified stability windows for even‑pulse families could be observed in such platforms, provided the nonlinearity is sufficiently strong and loss is minimized. They suggest a practical protocol: (a) generate two (or more) sech‑like input beams, (b) set the fractional exponent to a value in the range (2.5\lesssim\alpha\lesssim3.5), and (c) monitor the evolution with high‑resolution imaging to detect either persistent multi‑pulse propagation or the predicted drift/phase‑flip instabilities.

Overall Assessment

The work makes several substantive contributions to the theory of fractional nonlinear wave equations:

• It treats the fractional exponent (\alpha) as a genuine continuation parameter, thereby constructing a unified bifurcation diagram that links the well‑studied harmonic NLS ((\alpha=2)) with the less‑explored biharmonic case ((\alpha=4)).
• It provides a clear numerical and analytical characterization of tail behavior across the full (\alpha) range, confirming theoretical predictions of power‑law decay and oscillatory near‑tails.
• The identification of a critical threshold at (\alpha=1) for collapse aligns with dimensional analysis for fractional operators and is supported by both eigenvalue spectra and dynamical simulations.
• Most notably, the discovery of stability islands for even‑pulse multi‑solitons in the fractional regime is novel; such stabilization does not occur for integer‑order dispersion and highlights the richer spectral landscape introduced by fractional derivatives.
• The paper connects these theoretical insights to recent experimental capabilities, outlining realistic conditions under which the predicted phenomena could be observed.

Potential avenues for future research include extending the analysis to higher spatial dimensions (where the critical exponent changes), incorporating external potentials or lattice structures, and exploring other nonlinearities (e.g., saturable or PT‑symmetric terms). Overall, the manuscript offers a comprehensive and technically solid exploration of fractional soliton physics, with clear implications for both mathematical theory and experimental optics.