
Nonlinear Schr¨odinger equations with amplitude-dependent Wadati potentials
Dmitry A. Zezyulin∗
Department of Physics and Engineering, ITMO University, St. Petersburg 197101, Russia
(Dated: October 28, 2022)
Complex Wadati-type potentials of the form V(x) = −w2(x)+iwx(x), where w(x) is a real-valued
function, are known to possess a number of intriguing features, unusual for generic non-Hermitian
potentials. In the present work, we introduce a class of nonlinear Schr¨odinger-type problems which
generalize the Wadati potentials by assuming that the base function w(x) depends not only on the
transverse spatial coordinate but also on the amplitude of the field. Several examples of prospective
physical relevance are discussed, including models with the nonlinear dispersion or with the deriva-
tive nonlinearity. The numerical study indicates that the generalized model inherits the remarkable
features of standard Wadati potentials, such as the existence of continuous soliton families, the
possibility of symmetry-breaking bifurcations when the model obeys the parity-time symmetry, the
existence of constant-amplitude waves, and the eigenvalue quartets in the linear-instability spec-
tra. Our results deepen the current understanding of the interplay between nonlinearity and non-
Hermiticity and expand the class of systems which enjoy the exceptional combination of properties
unusual for generic dissipative nonlinear models.
I. INTRODUCTION
The combination of nonlinearity and non-Hermiticity
can have a dramatic impact on the properties of a waveg-
uiding system. In particular, the presence of a non-
Hermitian complex potential (which in physical terms
takes into account the energy exchange with the environ-
ment) can heavily alter the structure of stationary local-
ized nonlinear modes or solitons propagating in the sys-
tem. While in conservative waveguides the solitons exist
as continuous families parameterized by an ‘internal’ pa-
rameter, such as the soliton frequency or amplitude, the
dissipative solitons most usually exist as isolated points
[1, 2] which, from the dynamical point of view, behave as
attractors (provided that the soliton is stable). A proto-
typical model that illustrates this dissimilarity between
conservative and dissipative systems is the generalized
nonlinear Schr¨odinger equation (GNLSE)
iΨt=−Ψxx +V(x)Ψ + W(x)F(|Ψ|2)Ψ,(1)
where complex-valued functions V(x) and W(x) can be
referred to as a linear and nonlinear potential, respec-
tively, and real-valued function F(·), such that F(0) = 0,
specifies the nonlinearity. Equation (1) is a canonical
model describing evolution of nonlinear waves in various
physical settings. In particular, in optical applications, Ψ
corresponds to the dimensionless amplitude of the elec-
tric field. Its dependence on coordinate tdescribes evo-
lution of the pulse along the propagation direction, and
xis the transverse coordinate. Effective complex optical
potentials correspond to the presence of spatial regions
with gain and loss [3, 4], which, in particular, can be
implemented using coherent multilevel atoms driven by
external laser fields [5, 6].
∗email: d.zezyulin@gmail.com
When the imaginary parts of both potentials vanish,
the model becomes conservative, and its stationary non-
linear modes Ψ(x, t) = e−iµtψ(x) can be parameterized
by the continuous change of the real frequency µ. How-
ever, if at least one of the potentials is complex, then
the continuous families of solitons typically disappear.
At the same time, there exist at least two types of non-
Hermitian potentials that do support continuous fami-
lies of nonlinear modes (see a more detailed discussion in
[7]). Systems of the first type obey the parity-time (PT )
symmetry [3, 8]: in this case real and imaginary parts of
functions V(x) and W(x) are even and odd functions of
x, respectively [9, 10]. The second type corresponds to
the so-called Wadati-type potentials [11] for which
V(x) = −w2(x) + iwx(x), W (x)≡1,(2)
where w(x) is a real-valued differentiable function which
is not required to have any particular symmetry (at the
same time, if w(x) is even, then the resulting Wadati po-
tential V(x) becomes PT symmetric, i.e., the two types
of potentials are partially overlapping). In what follows,
it will be convenient to say that the GNLSE (1) with
potentials (2) corresponds to linear Wadati potentials.
Apart from the existence of continuous families of non-
linear modes [12], linear Wadati potentials are known
to support several other remarkable features. When the
nonlinearity is absent, i.e., F≡0, linear Wadati poten-
tials can have all-real spectrum of eigenvalues [13] and
undergo distinctive phase transition from all-real to com-
plex spectrum through an exceptional point or a self-dual
spectral singularity [14, 15]. Returning to the nonlinear
setup, Wadati-type potentials support the existence of
constant-amplitude nonlinear waves [16], feature unusual
dynamical behavior near the phase-transition threshold
[17], and eigenvalue quartets in the linear-stability spec-
trum [18]. Additionally, when a Wadati potential is PT
symmetric, it allows for bifurcations of families of non-
PT -symmetric modes [19], which is impossible for PT -
symmetric potentials of general form. While the the-
arXiv:2210.15552v1 [nlin.PS] 27 Oct 2022