We investigate the spectral and dynamical properties of the fractional nonlinear Schrödinger equation with harmonic confinement. In this setting, the classical Laplacian is replaced by its fractional power (-∂x2)α/2 with α∈(1,2], introducing nonlocal, Lévy-type dispersion. This modification fundamentally alters the balance between nonlinearity, dispersion, and trapping, reshaping both the structure and stability of stationary states. Using a Fourier pseudo-spectral discretization, we compute stationary branches as functions of the temporal frequency Ω in focusing (σ=+1) and defocusing (σ=-1) regimes, and assess spectral stability via the linearized eigenvalue problem. Direct simulations, performed with split-step and exponential time-differencing integrators, confirm these predictions and reveal α-dependent transitions between coherent oscillations, bounded breathing dynamics, and decoherence or fragmentation. Our results show that decreasing α systematically shifts bifurcation curves, fragments stability windows for excited states, and amplifies instability in the focusing regime, while supporting robust coherence in the defocusing case. Beyond clarifying how harmonic confinement mediates the interplay between nonlinearity and fractional dispersion, the study also provides benchmarks for numerical treatments of fractional operators and points toward potential applications in nonlinear optics, Bose-Einstein condensates, and anomalous transport phenomena.
Kusdiantara et al. (2026) studied this question.