
Fractional Integrable and Related Discrete Nonlinear Schr¨odinger Equations
Mark J. Ablowitz,1Joel B. Been,2, 3, ∗and Lincoln D. Carr2, 3, 4
1Department of Applied Mathematics, University of Colorado, Boulder, Colorado 80309, U.S.A.
2Department of Applied Mathematics and Statistics,
Colorado School of Mines, Golden, Colorado 80401, U.S.A.
3Department of Physics, Colorado School of Mines, Golden, Colorado 80401, U.S.A.
4Quantum Engineering Program, Colorado School of Mines, Golden, Colorado 80401, U.S.A.
Integrable fractional equations such as the fractional Korteweg-deVries and nonlinear Schr¨odinger
equations are key to the intersection of nonlinear dynamics and fractional calculus. In this
manuscript, the first discrete/differential difference equation of this type is found, the fractional inte-
grable discrete nonlinear Schr¨odinger equation. This equation is linearized; special soliton solutions
are found whose peak velocities exhibit more complicated behavior than other previously obtained
fractional integrable equations. This equation is compared with the closely related fractional aver-
aged discrete nonlinear Schr¨odinger equation which has simpler structure than the integrable case.
For positive fractional parameter and small amplitude waves, the soliton solutions of the integrable
and averaged equations have similar behavior.
Keywords: Integrable equations; Fractional calculus; Discrete nonlinear Schr¨odinger equation; Fourier split
step
INTRODUCTION
Integrable systems play a central role in nonlinear dy-
namics because they provide exactly solvable models for
important physical systems. Notable examples of inte-
grable equations are the Korteweg-deVries (KdV), ap-
plicable to shallow water waves, plasma physics, and
lattice dynamics among others [1–3], and the nonlinear
Schr¨odinger (NLS) equation, which finds applications in
nonlinear optics, Bose-Einstein condensates, spin waves
in ferromagnetic films, plasma physics, water waves, etc.
[2–5]. These integrable nonlinear evolution equations
have an infinite number of conservation laws and soli-
ton solutions [6]. Solitons, the fundamental solutions
of such equations, are stable, localized nonlinear waves
which propagate without dispersing and interact elas-
tically with other solitons. Nonlinear integrable evolu-
tion equations have these surprising properties because
of their deep mathematical structure described by the
inverse scattering transform (IST).
IST is a method of solving nonlinear equations which
generalizes Fourier transforms. It solves these equations
in three steps: mapping the initial condition into scat-
tering space, evolving the intial data in scattering space
in time, and mapping the evolved scattering data back
to physical space; i.e., inverse scattering. This pro-
cess gives the solution to nonlinear equations solvable
by IST in terms of linear integral equations; such non-
linear equations are called integrable. Recently, we used
the mathematical structure of IST associated with for
the Korteweg-deVries (KdV) and nonlinear Schr¨odinger
(NLS) equations to develop a method of finding and solv-
ing the fractional KdV (fKdV) and fractional NLS (fNLS)
∗Corresponding author: joelbeen@mines.edu
equations [7]. We also showed that this method could
be applied to find fractional extensions of the modified
KdV, sine-Gordon, and sinh-Gordon equations [8]. These
equations represent the first known fractional integrable
nonlinear evolution equations with smooth (physical) so-
lutions and deeply connect the fields of nonlinear dynam-
ics and fractional calculus.
Fractional calculus is a mathematical structure origi-
nally designed to define non-integer derivatives and inte-
grals. It has sense become an effective way of modeling
many physical processes that exist in multi-scale media
[9, 10] or exhibit non-Gaussian statistics or power law
behavior [11–13]. A particularly important example is
anomalous diffusion, where the mean squared displace-
ment is proportional to tα,α > 0 [11, 14–16]. Trans-
port that follows this rule has been observed in biology
[17–20], amorphous materials [21–23], porous media [24–
27], climate science [28], and attenuation in materials [29]
amongst others. As we have shown, the merger of frac-
tional and nonlinear characteristics in integrable equa-
tions such as fKdV and fNLS predict anomalous disper-
sion, where the velocity and amplitude of solitonic solu-
tions are related by a power law [7].
In this article, we demonstrate how the method in-
troduced in Ref. [7] can be applied to discrete (or
differential-difference) systems to define integrable dis-
crete fractional nonlinear evolution equations by present-
ing a fractional generalization of the integrable discrete
nonlinear Schr¨odinger (IDNLS) equation. We do this by
demonstrating the three key mathematical ingredients of
our method — IST, power law dispersion relations, and
completeness relations — for the Ablowitz-Ladik (AL)
discrete scattering problem. In the linear limit, the frac-
tional IDNLS (fIDNLS) equation is a discretization of the
fractional Schr¨odinger equation which was derived with
Feynman path integrals over L´evy flights [30, 31].
arXiv:2210.01229v1 [nlin.SI] 3 Oct 2022