Prerpints.org logo
Preprint
Review

Discrete and Semi-discrete Multidimensional Solitons and Vortices: Established Results and Novel Findings

Altmetrics

Downloads

134

Views

40

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

08 January 2024

Posted:

09 January 2024

You are already at the latest version

Alerts
Abstract
This article presents a concise survey of basic discrete and semi-discrete nonlinear models which produce two- and three-dimensional (2D and 3D) solitons, and a summary of main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ⁽²⁾) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of basic species of multidimensional discrete modes, including fundamental and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.
Keywords: 
Subject: Physical Sciences  -   Theoretical Physics

1. Introduction

1.1. Discrete Nonlinear Schrödinger (DNLS) Equations

1.1.1. The Basic Equation

Commonly adopted models of physical media are based on linear and nonlinear partial differential equations, such as the Gross-Pitaevskii (GP) equations for the mean-field wave function ψ x , y , z ; t in atomic Bose-Einstein condensates (BECs) [1], and nonlinear Schrödinger (NLS) equations for the amplitude of the electromagnetic field in optical waveguides [2,3]. The scaled form of the three-dimensional (3D) GP/NLS equation is
i ψ t = − ( 1 / 2 ) ∇ 2 ψ + σ | ψ | 2 ψ + U x , y , z ψ ,
where σ = + 1 and − 1 correspond to the self-defocusing and focusing signs of the cubic nonlinearity, and U x , y , z is an external potential. In the application to optics, time t is replaced, as the evolution variable, by the propagation distance z, while the original coordinate z is then replaced by the temporal one, τ = t − z / V gr , where V gr is the group velocity of the carrier wave [2]. In optics, the effective potential may be two-dimensional (2D), − U ( x , y ) , which represents a local variation of the refractive index in the transverse plane.
In many cases the potential is spatially periodic, such as one induced by optical lattices (OLs) in BEC [4,5], or by photonic crystals which steer the propagation of light in optics [6],
U latt x , y , z = − ε cos ( 2 π L / x ) + cos ( 2 π L / y ) + cos ( 2 π L / z ) ,
as well as its 2D and 1D reductions. A deep lattice potential, which corresponds to large amplitude ε in Eq. (2), splits the continuous wave function into a set of “droplets” trapped in local potential wells, which are linearly coupled by tunneling. Accordingly, in the framework of the tight-binding approximation, the NLS equation is replaced by a discrete NLS (DNLS) equation with the linear coupling between adjacent sites of the discrete lattice (nearest neighbors). This equation was derived, in the 1D form, for arrays of optical fibers [7,8,9,10] and plasmonic nanowires [11], as well as for BEC loaded in a deep OL potential [12]:
i ψ ˙ l , m , n = − ( 1 / 2 ) u l + 1 , m , n + ψ l − 1 , m , n − 2 ψ l , m , n + ψ l , m + 1 , n + ψ l , m − 1 , n − 2 ψ l , m , n + ψ l , m , n + 1 + ψ l , m , n − 1 − 2 ψ l , m , n + σ ψ l , m , n 2 ψ l , m . n + V l , m , n ψ l , m , n ,
where the overdot stands for d / d t , and the set of integer indices, l , m , n , replaces continuous coordinates x , y , z in Eq. (1). Here, potential V l , m , n is a possible smooth addition to the deep lattice potential which imposes the discretization. The full DNLS equation (3) is often reduced to 2D and 1D forms. 1D lattices are sometimes built in the form of zigzag chains, making it relevant to add couplings between the next-nearest neighbors [13,14]. 2D lattices with similar additional couplings were elaborated too [15].
Unlike the continuous NLS equation (1), the nonlinearity signs σ = ± 1 in Eq. (3) are actually equivalent to each other. Indeed, the staggering transformation of the discrete wave function,
ψ l , m , n ( t ) ≡ ( − 1 ) l + m + n exp − 6 i t ψ ˜ l , m , n * ( t )
where * stands for the complex-conjugate, transforms Eq. (3) with σ = − 1 into the same equation for ψ ˜ l , m , n with σ = + 1 . For the 2D and 1D DNLS equations, exp − 6 i t in the corresponding versions of staggering transformation (4) is replaced, respectively, by exp − 4 i t and exp − 2 i t .

1.1.2. Extended Equations

The DNLS equation and its extensions, such as systems of coupled DNLS equations [16,17], constitute a class of models with a large number of physical realizations [10,18,20,21,22]. They have also drawn much interest as subjects of mathematical studies [23]. One of incentives for this interest is the fact that the discreteness arrests the development of the critical and supercritical collapse, which is driven by the self-focusing nonlinear term with σ = − 1 in the 2D and 3D continuous NLS equations (1), respectively. The collapse leads to the emergence of singular solutions in the form of infinitely tall peaks, after a finite evolution time [3]. Naturally, the discreteness causes the arrest of the collapse, replacing it by a quasi-collapse [24], when the width of the shrinking peak becomes comparable to the spacing of the DNLS lattice.
The possibility of the collapse destabilizes formal 2D and 3D soliton solutions which are produced by Eq. (1), therefore a challenging problem is prediction and experimental realization of physical settings than make it possible to produce stable multidimensional solitons [25,26]. Thus, the discreteness provides a general method for the stabilization of 2D and 3D solitons.

The Gross-Pitaevskii (GP) Equations Amended by Effects of Quantum Fluctuations

Another promising possibility for the suppression of the collapse is offered by binary BEC, in which the cubic inter-component attraction creates 3D soliton-like states in the form of “quantum droplets" (QDs), while the development of the supercritical collapse is arrested by the self-repulsive quartic term that takes into account a correction to the mean-field nonlinearity produced by quantum fluctuations (known as the celebratedLee-Huang-Yang effect [28]). For the symmetric binary condensate, with identical wave functions ψ of its components, the accordingly amended scaled GP equations (in the absence of the trapping potential) was derived by Petrov [29]
i ψ t = − ( 1 / 2 ) ∇ 2 ψ − | ψ | 2 ψ + | ψ | 3 ψ .
Surprisingly quickly, the QD modes predicted by Eq. (5) have been created experimentally, in the quasi-2D [30,31] and full 3D [32] forms. The reduction of the spatial dimension to 2D and 1D replaces Eq. (5) by the following GP equations, respectively [33]:
2 D : i ψ t = − ( 1 / 2 ) ∇ 2 ψ + | ψ | 2 ln ( | ψ | 2 ) · ψ ,
1 D : i ψ t = − ( 1 / 2 ) ψ x x + | ψ | 2 ψ − | ψ | ψ .
Note, in particular, that in the 1D equation (7) the quantum correction is represented by the attractive term − | ψ | ψ , on the contrary to its repulsive counterpart + | ψ | 3 ψ in the 3D equation (5). For this reason, the usual mean-field cubic term is taken in Eq. (7) with the self-repulsion sign, to make it possible to study effects of competition of the quadratic self-attraction and cubic repulsion [34]. A semi-discrete version of Eq. (7) is considered below, see Eq. (66).

The Ablowitz-Ladik (AL) and Salerno-Model (SM) Equations

The 1D continuous NLS equation without an external potential is integrable by means of the inverse-scattering transform, with either sign of the nonlinearity, σ = ± 1 [35,36,37,38], although 2D and 3D extensions of the NLS equation are nonintegrable. The straightforward discretization destroys the integrability of the 1D NLS equation [39,40]. Nevertheless, the NLS equation admits a specially designed 1D discretization, which leads to an integrable discrete model, viz., the Ablowitz-Ladik (AL) equation [41]:
i ψ ˙ n = − ψ n + 1 + ψ n − 1 1 + μ ψ n 2 ,
where positive and negative values of the nonlinearity coefficient, μ , correspond to the self-focusing and defocusing, respectively. Integrable discrete equations, such as the AL one, are exceptional models which provide exact solutions for discrete solitons [42].
Equation (8) gives rise to an exact solution for solitons in the case of μ > 0 . Setting μ ≡ + 1 by means of rescaling, the solution is
ψ n ( t ) = sinh β sech β ( n − ξ ( t ) ) exp i α n − ξ ( t ) − i φ ( t ) ,
where β and α are arbitrary real parameters that determine the soliton’s amplitude, A ≡ sinh β , its velocity, V ≡ ξ ˙ = 2 β − 1 sinh β sin α , and overall frequency Ω ≡ φ ˙ = − 2 cosh β cos α + ( α / β ) sinh β sin α . The existence of exact solutions for traveling solitons in the discrete system is a highly nontrivial property of the AL equation, which follows from its integrability. Generically, motion of a discrete soliton through a lattice is braked by emission of radiation, even if this effect may seem very weak in direct simulations [65]. Another integrable discrete model which admits exact solutions for moving solitons is the Toda-lattice equation for real coordinates x n ( t ) of particles with unit mass and exponential potential of interaction between adjacent ones [66]:
x ¨ n + exp − x n + 1 − x n − exp − x n − x n − 1 = 0 .
Considerable interest was also drawn to the nonintegrable combination of the AL and DNLS equations, in the form of the Salerno model (SM) [43], with an additional onsite cubic term, different from the intersite one in Eq. (8):
i ψ ˙ n = − ψ n + 1 + ψ n − 1 1 + μ ψ n 2 − 2 ψ n 2 ψ n .
Here, the signs and magnitude of the onsite nonlinearity coefficient are fixed by means of the staggering transformation (4) and rescaling. The SM finds a physical realization in the context of the Bose-Hubbard model, which represents the BEC loaded in a deep OL, taking into regard the nonlinearity of the intersite coupling [44].
The AL and SM equations (8) and (11) conserve the total norm N, whose definition is different from the straightforward one for the DNLS equation, given below by Eq. (27); namely,
N AL , SM = ( 1 / μ ) ∑ n ln 1 + μ | ψ n | 2
[41,107]. The Hamiltonians of the AL and SM equations, which are dynamical invariants too, are
H AL = − ∑ n ψ n ψ n + 1 * + ψ n + 1 ψ n * ,
H SM = − ∑ n ψ n ψ n + 1 * + ψ n + 1 ψ n * + ( 2 / μ ) | ψ n | 2 + ( 2 / μ ) N SM .
In particular, the ostensible “simplicity” of Hamiltonian (13) is related to the complexity of the respective Poisson brackets (symplectic structure), which determine the evolution equations for ψ n as d ψ n / d t = H , ψ n . For the AL and SM models, the Poisson brackets, written for a pair of arbitrary functions B ψ n , ψ n * and C ψ n , ψ n * , are [43,44]
B , C = i ∑ n ∂ B ∂ ψ n ∂ C ∂ ψ n * − ∂ B ∂ ψ n * ∂ C ∂ ψ n 1 + μ ψ n 2 .
It is also relevant to consider the continuum limit of the SM, which is introduced by approximating the intersite combination of the discrete fields by a truncated Taylor’s expansion,
ψ n ( t ) ≡ e 2 i t Ψ ( x , t ) , Ψ x = n ± 1 , t ≈ Ψ x = n ± Ψ x | x = n + ( 1 / 2 ) Ψ x x | x = n ,
where Ψ ( x ) is considered as a function of continuous coordinate x, whose integer values coincide with the discrete lattice coordinate n. The substitution of approximation (16) in Eq. (11) leads to a generalized (nonintegrable) NLS equation [55]
i Ψ t = − 1 + μ Ψ 2 Ψ x x − 2 1 + μ Ψ 2 Ψ .
Equation (17) conserves the total norm and Hamiltonian, which are continuum counterparts of expressions (12) and (14):
N AL cont = 1 μ ∫ − ∞ + ∞ d x ln 1 + μ | Ψ | 2 ,
H SM cont = ∫ − ∞ + ∞ d x Ψ x 2 − 2 1 μ + 1 | Ψ | 2 + 2 μ N AL cont
It is relevant to mention that the general approximation opposite to the continuum limit is the anti-continuum limit [18,56]. This approach starts with the limit form of the DNLS equation, in which the linear couplings between adjacent sites are dropped. Then, one can try to construct various states, including solitons, by introducing an input composed of simple solutions of the single-site equations corresponding to Eq. (3), viz., i ψ ˙ l , m , n = σ ψ l , m , n 2 ψ l , m . n + V l , m , n ψ l , m , n , at a finite set of sites, and keeping the zero solution at all others. The single-site “simple solutions" are
ψ l , m , n = a l , m , n exp − i a l , m , n 2 t − i V l , m , n t ,
where a l , m , n is an arbitrary set of complex amplitudes. Next, one reintroduces weak linear intersite couplings and attempts to identify nontrivial solutions that may thus appear from the finite-set input composed of the single-site solutions (20).

Self-Trapping in Lattices with the Self-Repulsion Strength Growing from the Center to Periphery

DNLS equations with the onsite self-repulsive nonlinearity, corresponding to σ > 0 in Eq. (3), may support discrete-soliton (self-trapped) states without the resort to the staggering transform (4) if the local self-repulsion strength is made a function of the lattice coordinates, growing fast enough from the center to periphery. Originally, this option was elaborated in the framework of the continuum NLS and GP equations in the space of dimension D, with the local self-defocusing (repulsion) coefficients growing at r → ∞ (r is the radial coordinate) faster than r D [45]. In terms of the 1D and 2D DNLS equations, similar settings were introduced in Refs. [46] and [47], with the site-dependent self-attraction coefficients, σ n 1 D = σ 0 exp α | n | and σ m , n 2 D = σ 0 exp α | m | + | n | , respectively, with positive constants σ 0 and α . In the 2D model, solutions were constructed, and their stability analyzed, for fundamental, dipole, quadrupole, and vortical discrete solitons [45].

DNLS Equations with Long-Range Dipole-Dipole and Quadrupole-Quadrupole Intersite Interactions

It is well known that atomic BEC formed of ultracold atoms carrying permanent magnetic moments feature specific dynamical effects due to the long-range interactions between atomic moments [48]. This fact suggests to combine the dipole-dipole interactions and a deep OL potential, thus introducing DNLS equations with the nonlocal (long-range) coupling between lattice sites. In the 2D setting, this model gives rise to different forms of the DNLS equations. The simplest setup corresponds to the case when the atomic moments are polarized by external dc magnetic field perpendicular to the system’s plane. In this case, the dipole-dipole interactions amount to the isotropic nonlocal repulsion, accounted for by the respective interaction coefficient Γ > 0 [49]
i ∂ ψ m , n ∂ t = − 1 2 ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n + σ ψ m , n 2 + Γ ∑ m ′ , n ′ ≠ m , n m − m ′ 2 + n − n ′ 2 − 3 / 2 ψ m ′ , n ′ 2 ψ m , n ,
where σ is the same coefficient of the onsite self-interaction as in Eq. (3). A more sophisticated setup corresponds to the atomic magnetic moments polarized parallel to the system’s plane. In the former case, the nonlocal term term in the respective DNLS equation is anisotropic, being attractive in one in-plane direction and repulsive in the other, cf. Ref. [51]:
i ∂ ψ m , n ∂ t = − 1 2 ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n + σ ψ m , n 2 + Γ ∑ m ′ , n ′ ≠ m , n n − n ′ 2 − 2 m − m ′ 2 m − m ′ 2 + n − n ′ 2 5 / 2 ψ m ′ , n ′ 2 ψ m , n .
The analysis reported in Ref. [49] demonstrates that the nonlocal repulsion, added to Eq. (21), helps to stabilize discrete solitons with embedded vorticity. Solutions of Eq. (22) for anisotropic vortex solitons can be found too, but they are completely unstable [49].
A 2D DNLS model which combines the local onsite nonlinearity and long-range interaction between particles carrying permanent quadrupole electric moments was elaborated in Ref. [52]:
i ∂ ψ m , n ∂ t = − 1 2 ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m n + σ | ψ m n | 2 + Γ ∑ m ′ , n ′ ≠ m , n ( n − n ′ ) 2 − 4 ( m − m ′ ) 2 [ ( m − m ′ ) 2 + ( n − n ′ ) 2 ] 7 / 2 | ψ m ′ n ′ | 2 ψ m , n ,
cf. Eq. (21). This model also gives rise to families of stable 2D discrete solitons [52].

The 2D Discrete Second-Harmonic-Generating ( χ ( 2 ) ) system

The quadratic (alias χ ( 2 ) ) nonlinearity is a fundamentally important effect which gives rise to coherent generation of the second harmonic in optics. In terms of the 2D spatial-domain propagation in a continuum material, the standard χ ( 2 ) system for amplitudes ψ x , y , z and ϕ x , y , z of the fundamental-frequency (FF) and second-harmonic (SH) waves is [53]
i ψ z = − ( 1 / 2 ) ∇ 2 ψ − ψ * ϕ , 2 i ϕ z = − ( 1 / 2 ) ∇ 2 ϕ − Q ϕ − 1 2 ψ 2 ,
where z is the propagation distance, the paraxial-diffraction operator ( 1 / 2 ) ∇ 2 acts on the transverse coordinates x , y , Q is a real mismatch parameter, and * stands for the complex conjugate. The discretized version of Eqs. (24), which represents, in the tightly-binding approximation, light propagation in a photonic crystal made of the χ ( 2 ) material, is [54]
i d ψ m , n d z = − C 1 ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m n − ψ m , n * ϕ m , n , 2 i d ϕ m , n d z = − C 2 2 ϕ m + 1 , n + ϕ m − 1 , n + ϕ m , n + 1 + ϕ m , n − 1 − 4 ϕ m n − Q ϕ m , n − 1 2 ψ m , n 2 ,
where C 1 and C 2 are effective lattice-coupling constants for the FF and SH waves. The role of the conserved norm of the discrete χ ( 2 ) system is played by the Manley-Rowe invariant, i.e., the total optical power, I = ∑ m , n ψ m , n 2 + 4 ϕ m , n 2 .
An essential property of 2D discrete solitons produced by Eqs. (25) is their mobility [54]. In this connection, it is relevant to mention that, while the development of the quasi-collapse in the 2D discrete NLS equation with the cubic self-attraction is arrested by the underlying lattice structure, the quasi-collapse strongly pins the 2D solitons to the same structure, and thus makes them immobile. On the other hand, the χ ( 2 ) nonlinearity does not give rise to the collapse in the 2D (and 3D) space, therefore 2D χ ( 2 ) solitons do not demonstrate a trend for strong pinning, remaining effectively mobile robust localized modes [54].

1.2. One-Dimensional DNLS Solitons

1.2.1. Fundamental Solitons

In the 1D setting, the model of basic interest is the DNLS equation with self-attraction, which corresponds to σ = − 1 in the 1D version of Eq. (3), without the external potential ( V l , m , n = 0 ):
i ψ ˙ n = − ( 1 / 2 ) ψ n + 1 + ψ n − 1 − 2 ψ n − ψ n 2 ψ n .
This equation conserves two dynamical invariants, viz., the total norm,
N = ∑ n = − ∞ + ∞ ψ n 2 ,
and Hamiltonian (energy),
H = ∑ n = − ∞ + ∞ ( 1 / 2 ) ψ n − ψ n − 1 2 − ( 1 / 4 ) ψ n 4 .
Stationary solutions to Eq. (26) with real frequency ω are looked for as
ψ n ( t ) = e − i ω t u n ,
with real amplitudes u n satisfying the discrete equation,
ω u n = − ( 1 / 2 ) u n + 1 + u n − 1 − 2 u n − u n 3 .
Note that Eq. (30) can be derived by varying its Lagrangian,
L = ∑ n = − ∞ + ∞ ( 1 / 4 ) u n − u n − 1 2 − u n 4 − ω u n 2 ,
with respect to the discrete real field u n .
A fundamental property of the DNLS equation (26) with the self-attractive onsite nonlinearity is the modulational instability (MI) of its spatially homogeneous continuous-wave (CW) state [57], ψ n = a exp i a 2 t , with an arbitrary amplitude a [cf. Eq. (29)]. MI breaks the CW state into a chain of discrete solitons [10]. Analytical solutions for these solitons are not available, as the DNLS equation is not integrable. The solitons can be readily found in a numerical form, and studied in the framework of the variational approximation (VA) [23]. The VA is based on a particular ansatz, i.e. an analytical expression which aims to approximate the solution [58]. The only discrete ansatz for which analytical calculations are feasible is represented by the exponential function [59,60,61,62], namely,
u n onsite = A exp − a | n | ,
with a > 0 . The corresponding norm, calculated as per Eq. (27), is
N ansatz = A 2 coth a .
Note that ansatz (32) is appropriate for strongly and moderately discrete solitons, as shown in Figure 1, but not for broad (quasi-continuum) ones, which may be approximated by the commonly known soliton solution of the NLS equation (the 1D version of (1) with U = 0 ),
ψ x , t = η sech η x − ξ exp i η 2 t ,
with width η − 1 which must be large in comparison with the discreteness spacing, η − 1 ≫ 1 , and central coordinate ξ .
The substitution of ansatz (32) in Lagrangian (31) produces the corresponding VA Lagrangian:
L VA = ( A 2 / 2 ) coth ( a / 2 ) − ( A 4 / 4 ) coth 2 a − ω A 2 coth a .
Then, for given ω < 0 (solitons do not exist for ω > 0 ), the squared amplitude, A 2 , and inverse width, a, of the discrete soliton are predicted by the Euler-Lagrange equations,
∂ L eff ∂ A 2 = ∂ L eff ∂ a = 0 .
This corresponding system of algebraic equations for A 2 and a can be easily solved numerically. The VA produces an accurate predictions for the solitons, as shown in Figure 1 and Ref. [63].
Rigorous justification of the VA was elaborated in Ref. [64]. Furthermore, the VA and a full numerical solution of Eq. (26) demonstrate that the entire family of the discrete solitons is stable [23].

1.2.2. Higher-Order One-Dimensional Modes: Twisted Discrete Solitons and Bound States

In addition to the fundamental (single-peak) solitons outlined above, Eq. (30) admits stable second-order states in the form of twisted modes, which are subject to the antisymmetry condition, u n = − u 1 − n [67]. Such states exist and are stable only in a strongly discrete form, vanishing in the continuum limit. In particular, the above-mentioned anti-continuum approximation is appropriate for the construction of the twisted modes.
Stable 1D DNLS solitons may form bound states, which also represent higher-order modes of the DNLS equation. They are stable in the out-of-phase form, i.e., with opposite signs of the constituent solitons [68,69], which resembles the structure of the twisted modes (however, the tight antisymmetric structure of the twisted modes cannot be considered as a bound state of fundamental solitons), the same being true for 2D DNLS solitons [70]. Note that stationary bound states of fundamental solitons do not exist in the continuum limit.

1.2.3. 1D Solitons in the Salerno Model (SM)

The substitution of the usual ansatz (29) in Eq. (11) produces a stationary discrete equation for real u n :
ω u n = − u n + 1 + u n − 1 1 + μ u n 2 − 2 u n 3 ,
cf. (30). Discrete solitons produced by the SM equation (11) with μ > 0 , i.e. with noncompeting intersite and onsite self-focusing nonlinearities, were investigated by means of numerical methods [107,108,109]. It was found that the SM gives rise to the 1D solitons at all positive values of μ .
Another option is to consider the SM with μ < 0 , which features competing nonlinearities, as the intersite cubic term, with coefficient μ < 0 in Eq. (11), which accounts for the nonlinear coupling between adjacent sites of the lattice, and the onsite cubic term in Eq. (11) represent, respectively, repulsive and attractive nonlinear interactions. This version of the SM gives rise to families of discrete solitons, in the usual form (29), with ω < 0 and real u n , of two different types. One family represents ordinary discrete solitons, similar to those generated by the DNLS equation. Another family represents cuspons, featuring higher curvature of their profile at the center. A small subfamily of ordinary solitons produced by the SM with the competing nonlinearities is unstable, while all cuspons are stable.
As mentioned above, antisymmetric bound states of DNLS solitons are stable, while symmetric bound states are unstable [68,69]. The same is true for bound states of ordinary discrete solitons in the SM [55]. However, in the framework of the SM with the competing nonlinearities, the situation is exactly opposite for the cuspons: their symmetric and antisymmetric bound states are stable and unstable, respectively [55].

1.3. The Subject and Structure of the Present Article

The above-mentioned reviews [10,18,20,21,22] and [23] produce a comprehensive survey of theoretical and experimental results for discrete solitons in various 1D systems. The objective of this article is to produce a relatively brief summary of results for multidimensional (chiefly, two-dimensional) discrete and semi-discrete solitons, which were considered in less detail in previous reviews and, on the other hand, draw growing interest in the context of the current work with 2D and 3D solitons in diverse physical contexts [25,26]. In this context, the presence of the two or three coordinates makes it also possible to define semi-discrete states as ones which are discrete in one direction and continuous in the perpendicular one [71,72,73,74,77]. The article chiefly represents theoretical results, but some experimental findings for quasi-discrete 2D solitons in photonic lattices [79,80] are included too.
The review presented below does not claim to be comprehensive. The comprises results that are produced by conservative models of the DNLS types (including, in particular, the 2D SM). Discrete models of other types – in particular, those similar to the Toda lattice, see Eq. (10), Fermi-Pasta-Ulam-Tsingou lattices [81], and Frenkel-Kontorova systems [82] – are not considered here. Dissipative systems are not considered either, except for a 2D model with the parity-time ( PT ) symmetry [83], see Eq. (70) below.
The rest of the article is arranged as follows. Basic results for fundamental (zero-vorticity) and vortex solitons, as well as bound states of such solitons, produced by the 2D DNLS equation and its generalizations, are summarized in Section II, which is followed, in Section III, by brief consideration of fundamental and vortex solitons in the 2D SM (Salerno model). Section IV addresses discrete solitons of the semi-vortex and mixed-mode types in the 2D spin-orbit-coupled (SOC) system of GP equations for a two-component BEC. Basic results for discrete self-trapped modes produced by 3D DNLS equations, including fundamental and vortex solitons, along with skyrmions, are presented in Section V. The findings for 2D semi-discrete systems, again including fundamental and vortex solitons, supported by combined quadratic-cubic and cubic-quintic nonlinearities (that are relevant for BEC and optics, respectively), are summarized in Section VI. This section also addresses transverse mobility of confined spatiotemporal modes in an array of optical fibers with the intrinsic cubic self-focusing (Kerr nonlinearity). Fundamental and vortex solitons produced by a PT -symmetric discrete 2D coupler with the cubic nonlinearity are considered in Section VII. The article is concluded by Section VIII, which, in particular, suggests directions for the further work in this area, and mentions particular topics which are not included in the present review.

2. Two-Dimensional (2D) Nonlinear-Schrödinger Lattices: Fundamental and Vortex Solitons, and Their Bound States

2.1. Vortex Solitons: Theoretical and Experimental Results

The basic 2D cubic DNLS equation is the 2D version of Eq. (3) with the self-attraction ( σ = − 1 ) and without the external potential:
i ψ ˙ m , n = − ( 1 / 2 ) ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n − ψ m , n 2 ψ m . n ,
cf. Eq. (26). The substitution of ψ m , n = exp − i ω t u m , n in Eq. (37) produces the stationary equation, where the stationary discrete wave function, u m , n , may be complex:
ω u m , n = − ( 1 / 2 ) u m + 1 , n + u m − 1 , n + u m , n + 1 + u m , n − 1 − 4 u m , n − u m , n 2 u m , n ,
cf. (30). Fundamental-soliton solutions to Eq. (38) can be predicted by means of VA [84,85], using an exponential ansatz, see (47) below (cf. Eq. (32) for the 1D soliton). More interesting in the 2D setting are discrete solitons with embedded vorticity, which were introduced in Ref. [86] (see also Ref. [87]). Vorticity, alias the topological charge, or winding number, is defined as Δ φ / 2 π , where Δ φ is a total change of the phase of the complex discrete function u m , n along a contour surrounding the vortex’ pivot. Stability is an important issue for 2D discrete solitons, because, in the continuum limit, the 2D NLS equation gives rise to the well-known Townes solitons [88], which are unstable against the onset of the critical collapse [3]. In the same limit, the Townes solitons with embedded vorticity (vortex rings [89]) are subject to much stronger instability against spontaneous splitting of the ring in fragments [90].
The lattice structure of the DNLS equation provides for stabilization of both fundamental (zero-vorticity) and vortex solitons [86]. A typical example of a stable 2D vortex soliton with topological charge S = 1 is displayed in Figure 2. 2D fundamental and vortex solitons, with topological charges S = 0 and 1, are stable in regions − ω > | ω cr ( S = 0 ) | ≈ 0.50 and − ω > | ω cr ( S = 1 ) | ≈ 1.23 , respectively [86], while the higher-order discrete vortices with charges S = 2 and 4 are unstable, being replaced by stable modes in the form of quadrupoles and octupoles [91]. The vortex solitons with S = 3 may be stable, but only in an extremely discrete form, viz., at − ω > | ω cr ( S = 2 ) | ≈ 4.94 . In agreement with what is said above, these results imply that all the solitons are unstable in the continuum limit, corresponding, in the present notation, to ω → 0 .
The experimentally relevant lattice structure may be anisotropic, with the linear combination ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n in Eq. (37) replaced by
α ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 2 1 + α ψ m , n , with anisotropy parameter α ≠ 1 . Effects of the anisotropy on the structure and stability of the fundamental and vortical discrete solitons were explored in Ref. [102].
The theoretically predicted 2D discrete solitons with vorticity S = 1 were experimentally created in Refs. [79] and [80], using a photorefractive crystal. Unlike uniform media of this type, where delocalized (“dark") optical vortices were originally demonstrated [92,93], those works made use of a deep virtual photonic lattice as a quasi-discrete structure supporting the self-trapping of nonlinear modes in the optical field with extraordinary polarization (while the photonic lattice was induced as the interference pattern of quasi-linear beams in the ordinary polarization). Intensity distributions observed in vortex solitons of the onsite- and intersite-centered types (i.e., with the vortex’ pivot coinciding with a site of the underlying lattice, or set between the sites, respectively), are displayed in Figure 3.
Another interesting result demonstrated (and theoretically explained) in deep virtual photonic lattices is a possibility of periodic flipping of the topological charge of a vortex soliton initially created with topological charge S = 2 [94]. Stable vortex solitons with S = 2 were created using a hexagonal virtual photonic lattice [95].

2.2. Bound States of 2D Discrete Solitons and Solitary Vortices

As mentioned above, stable 2D discrete solitons may form stable bound states, composed of two or several items. Vortex solitons may form bound states as well. This possibility and stability of the resulting bound states are determined by an effective potential of interaction between two identcial discrete solitons with intrinsic vorticity S, which are separated by large distance L. The potential can be derived from an asymptotic expression for exponentially decaying tails of the soliton. In the quasi-continuum approximation, it is u m , n ∼ m 2 + n 2 − 1 / 4 exp − − 2 ω m 2 + n 2 (recall soliton solutions to Eq. (38) exist for ω < 0 ). Then, the overlap of the tail of each soliton with the central body of the other one gives rise to the following interaction potential:
U int ( L ) ≈ const · cos δ · ( − 1 ) S L − 1 / 2 exp − 2 ω L ,
with const > 0 , where δ is the phase shift between the solitons. Thus, for the fundamental solitons with S = 0 , Eq. (39) predicts the attractive interaction between the in-phase solitons ( δ = 0 ), and repulsion between out-of-phase ones ( δ = π ). Accordingly, the interplay of the repulsive interaction with the effective Peierls-Nabarro potential, which is pinning the soliton to the underlying lattice [23], produces stable bound states of two or several mutually out-of-phase solitons, while the in-phase bound states are unstable. These predictions were confirmed by numerical results [70]. As an example, Figure 4 displays a numerically found stable bound state in the form of a string of three fundamental solitons with alternating phases.
For the pair of identical vortex solitons with S = 1 , Eq. (39) predicts the opposite result, viz., the repulsive interaction and stability of the ensuing bound states for in-phase vortices ( δ = 0 ), and the attraction leading to instability of the bound state in the case of δ = π . These predictions were also corroborated by numerical findings [70], see an example of a stable bound state of two identical vortex solitons in Figure 5.

2.3. 2D Discrete Solitons in Mini-Gaps of a Spatially Modulated Lattice

A specific class of self-trapped modes are gap solitons which may populate finite bandgaps in linear spectra of various nonlinear systems originating in optics and BEC [96,97,98]. While in most cases gap solitons are predicted theoretically and created experimentally in the context of matter waves [99] and optical pulses [100] in the continuum, they may naturally appear as discrete modes in mini-gaps, which are induced in the linear spectrum of lattice media by superimposed periodic spatial modulations (superlattice).
Such a 2D lattice model was introduced in Ref. [101], based on the following DNLS equation:
i ψ ˙ m , n + C m , n ψ m + 1 , n + ψ m − 1 , n + K m , n ψ m , n + 1 + ψ m , n − 1 + ψ m , n 2 ψ m , n = 0 ,
where the horizontal and vertical coupling constants are modulated as follows:
C m , n = 1 + Δ cos ( Q m ) , K m , n = 1 + Δ cos ( Q n ) ,
cf. Eq. (3). The superlattice represented by Eqs. (40) and (41) can be created by means of the technique used for making OLs in experiments with BEC [101].
Looking for solutions in the usual form, ψ m , n ( t ) = exp − i ω t u m , n , with real frequency (chemical potential) ω , the numerical analysis produces the linear spectrum of this system, including the usual semi-infinite bandgap and a pair of additional narrow mini-gaps. Further, a family of fundamental 2D discrete solitons populating the mini-gaps was furnished by the numerical solution of the full nonlinear system, being stable in a small section of the mini-gap, as shown in Figure 5. The stable 2D soliton displayed in panel Figure 5(a) features a typical shape of gaps solitons, with a number of small satellite peaks surrounding the tall central one [97]. Bound states of two and four fundamental solitons were found too, featuring weak instability [101].

2.4. 2D Discrete Solitons in a Rotating Lattice

Dynamics of BEC loaded in OLs rotating at angular velocity Ω , as well as the propagation of light in a twisted nonlinear photonic crystal with pitch Ω , is modeled by the 2D version of Eq. (1) including the lattice potential, with depth ε and period 2 π / k , written in the rotating reference frame:
i ψ t = − ( 1 / 2 ) ∇ 2 + Ω L ^ z ψ − ε cos k x + cos ( k y ) ψ + σ | ψ | 2 ψ ,
where L ^ z = i ( x ∂ y − y ∂ x ) ≡ i ∂ θ is the operator of the z-component of the orbital momentum, θ being the angular coordinate in the x , y plane. In the tight-binding approximation, Eq. (42) is replaced by the following variant of the DNLS equation [103]:
i ψ ˙ m , n = − ( C / 2 ) ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n − i Ω m ψ m , n + 1 − ψ m , n − 1 − n ψ m + 1 , n − ψ m − 1 , n + σ | ψ m , n | 2 ψ m , n ,
where C is the intersite coupling constant. In Ref. [103], stationary solutions to Eq. (43) were looked for in the usual form (29), fixing ω ≡ − 1 and varying C in (43) as a control parameter. Two species of localized states were thus constructed: off-axis fundamental discrete solitons, placed at distance R from the origin, and on-axis ( R = 0 ) vortex solitons, with topological numbers S = 1 and 2. At a fixed value of rotation frequency Ω , a stability interval for the fundamental soliton, 0 < C < C max ( fund ) ( R ) , monotonously shrinks with the increase of R, i.e., most stable are the discrete fundamental solitons with the center placed at the rotation pivot. Vortices with S = 1 are gradually destabilized with the increase of Ω (i.e., their stability interval, 0 < C < C max ( vort ) ( Ω ) , shrinks). On the contrary, a remarkable finding is that vortex solitons with S = 2 , which, as said above, are completely unstable in the usual DNLS equation with Ω = 0 , are stabilized by the rotation, in an interval 0 < C < C cr ( S = 2 ) ( Ω ) , with C cr ( S = 2 ) ( Ω ) growing as a function of Ω . In particular, C cr ( S = 2 ) ( Ω ) ≈ 2.5 Ω at small Ω [103].

2.5. Spontaneous Symmetry Breaking of the 2D Discrete Solitons in Linearly-Coupled Lattices

A characteristic feature of many nonlinear dual-core systems, built of two identical linearly-coupled waveguides with intrinsic self-attractive nonlinearity, is a spontaneous-symmetry-breaking (SSB) bifurcation, which destabilizes the symmetric ground state, with equal components of the wave function in the coupled cores, and creates stable asymmetric states. The SSB bifurcation takes place at a critical value of the nonlinearity strength, the asymmetric state existing above this value [105]. A system of linearly-coupled DNLS equations is a basic model for SSB in discrete settings. Its 2D form is [17]
i ϕ ˙ n = − ( 1 / 4 ) ϕ m + 1 , n + ϕ m − 1 , n + ϕ m , n + 1 + ϕ m , n − 1 − 4 ϕ m , n − ϕ m , n 2 ϕ m , n − K ψ m , n , i ψ ˙ n = − ( 1 / 4 ) ψ m + 1 , n + ψ m − 1 , n + ψ m , n + 1 + ψ m , n − 1 − 4 ψ m , n − ψ m , n 2 ψ m , n − K ϕ m , n ,
where ϕ m , n and ψ m , n are wave functions of the discrete coordinates m and n, and K > 0 represents the linear coupling between the cores. Stationary states with frequency ω are looked for as ϕ m , n , ψ m , n = exp − i ω t u m , n , v m . n . Real stationary fields in the two components are characterized by their norms,
E u , v = ∑ m , n = − ∞ + ∞ u m , n 2 , v m , n 2 ,
and the asymmetry of the symmetry-broken state is determined by parameter
r = E u − E v / E u + E v .
The system under the consideration can be analyzed by means of the VA, based on the 2D ansatz
u m , n , v m , n = A , B exp − a | m | + | n | ,
with inverse width a and amplitudes, A and B, of the two components (cf. the 1D ansatz (32)). The SSB is represented by solutions with A ≠ B . An example of a stable 2D discrete soliton is displayed in Figure 7(a), which corroborates accuracy of the VA. In Figure 7(b), the families of symmetric and asymmetric 2D discrete solitons is characterized by the dependence of asymmetry parameter r, defined as per Eq. (46), on the total norm, E ≡ E u + E v , see Eq. (45). Figure 7(b) demonstrates the SSB bifurcation of the subcritical type [106], with the two branches of broken-symmetry states originally going backward in the E direction, as unstable ones; they become stable after passing the turning point. Accordingly, Figure 7(b) demonstrates a bistability area, where symmetric and asymmetric states coexist as stable ones.

3. 2D Discrete Solitons in the Salerno Model (SM)

The 2D version of the SM was introduced in Ref. [111]:
i ψ ˙ n , m = − ψ n + 1 , m + ψ n − 1 , m + ψ n , m + 1 + ψ n , m − 1 × 1 + μ ψ n , m 2 − 2 ψ n , m 2 ψ n , m ,
Similar to its 1D version (11), Eq. (48) conserves the norm and Hamiltonian, cf. Eqs. (12) and (14),
N SM 2 D = ( 1 / μ ) ∑ m , n ln 1 + μ | ψ n , m | 2 ,
H SM 2 D = − ∑ n , m ψ n , m ψ n + 1 , m * + ψ n + 1 , m ψ n , m * + ψ n , m ψ n , m + 1 * + ψ n , m + 1 ψ n , m * + ( 2 / μ ) | ψ n , m | 2 + ( 2 / μ ) N SM 2 D .
The continuum limit of this model is the 2D equation which is an extension of its 1D counterpart (17):
i Ψ t + 1 + μ Ψ 2 Ψ x x + Ψ y y + 2 2 μ + 1 | Ψ | 2 Ψ = 0 .
Note that term in μ Ψ 2 Ψ x x + Ψ y y prevents the onset of the collapse in Eq. (51), making it possible to create stable 2D solitons.
2D discrete solitons are looked for as solutions to Eq. (51) in the usual form, ψ m , n ( t ) = e − i ω t Φ m , n . In the most interesting case of the competing nonlinearities, μ < 0 , the situation is similar to that outlined above for the one-dimensional SM: there are ordinary discrete solitons, which have their stability and instability regions, and 2D cuspons, which are entirely stable in their existence region. Typical 2D solitons of both types are displayed in Figure 8. Antisymmetric bound states of ordinary 2D discrete solitons, and symmetric complexes built of 2D cuspons, are stable, while the bound states of cuspons with opposite parities are unstable, also like in the 1D model.
Along with the fundamental solitons, the 2D SM with the competing nonlinearities gives rise to vortex-soliton modes which may be stable in narrow parameter regions [111]. Examples of onsite- and intersite-centered vortex solitons (alias vortex cross and vortex square, respectively) are presented in Figure 9. In the 2D SM with non-competing nonlinearities ( μ > 0 in Eq. (51)), vortex solitons are unstable, spontaneously transforming into fundamental ones and losing their vorticity. This transition is possible because the angular momentum is not conserved in the lattice system. The situation is different in the 2D SM with competing nonlinearities ( μ < 0 ), where unstable vortex modes transform into vortical breathers, i.e., persistently oscillating localized modes that keep the original vorticity.

4. Solitons of the Semi-Vortex (SV) and Mixed-Mode (MM) Types in the Discrete 2D Spin-Orbit-Coupling (SOC) System

Recently, much interest has been drawn to emulation of the spin-orbit-coupling (SOC) effect, which is well known in physics of semiconductors, in binary BEC [112,113,114,115]. While SOC is a linear effect, its interplay with the intrinsic mean-field nonlinearity of atomic BEC gives rise to predictions of new species of 1D, 2D, and 3D solitons [116]. In particular, the effectively 2D binary BEC with SOC of the Rashba type is modeled by the following system of coupled GP equations for two components ϕ ( ± ) of the pseudo-spinor wave function [117],
i ∂ ϕ ( + ) ∂ t = − 1 2 ∇ 2 ϕ ( + ) − ( | ϕ ( + ) | 2 + γ | ϕ ( − ) | 2 ) ϕ ( + ) + 2 λ ∂ ϕ ( − ) ∂ x − i ∂ ϕ ( − ) ∂ y , i ∂ ϕ ( − ) ∂ t = − 1 2 ∇ − 2 ϕ ( − ) − ( | ϕ ( − ) | 2 + γ | ϕ ( + ) | 2 ) ϕ ( − ) + 2 λ − ∂ ϕ ( + ) ∂ x − i ∂ ϕ ( + ) ∂ y .
In this system, SOC of the Rashba type is represented in by terms with coefficient λ , which couple the two equations through the first-order spatial derivatives. The system also includes the self- and cross-attractive nonlinearities, with scaled coefficients 1 and γ , respectively.
The system of coupled GP equations (52) maintains 2D solitons of two different types, namely, semi-vortices (SVs) and mixed modes (MMs) [117], The SV solitons, written in polar coordinates r , θ , have zero vorticity in one component, and vorticity + 1 or − 1 in the other:
ϕ 1 ( + ) = e − i ω t f 1 ( r 2 ) , ϕ 1 ( − ) = e − i ω t + i θ r f 2 ( r 2 ) ,
ϕ 2 ( + ) = e − i ω t − i θ r f 2 ( r 2 ) , ϕ 2 ( − ) = e − i ω t f 1 ( r 2 ) ,
where ω is the chemical potential, and f 1 , 2 r 2 are real functions which take finite values at r = 0 and exponentially decay ∼ sin 2 λ r , cos 2 λ r exp − 2 − ω − 2 λ 2 r at r → ∞ . The two SV solutions (53) and (), which are mirror images of each other, exist in the semi-infinite bandgap, ω < − 2 λ 2 .
The combination of zero and nonzero vorticities in the SV solutions (53) and () is exactly compatible with the structure of the coupled GP equations (52). On the contrary to this, MM solitons cannot be represented by an exact ansatz similar to Eqs. (53) and (), but they may be approximated by a linear combination of both types of the SVs, ϕ 1 ( + ) + ϕ 2 ( + ) , ϕ 1 ( − ) ( 1 ) + ϕ 2 ( + ) . An essential result is that the SVs and MMs are stable and represent the system’s ground state in the cases of γ < 1 and γ > 1 , respectively, i.e., when the self-attraction is stronger or weaker than the cross-attraction in Eqs. (52) [117]. On the other hand, the SVs and MMs are unstable, as excited states, in the opposite cases, i.e., γ > 1 and γ < 1 , respectively.
The discretized version of the SOC GP system (52), which corresponds to the spin-orbit-coupled binary BEC trapped in a deep OL potential, with discrete coordinates ( m , n ) , was introduced in Ref. [118]:
i ∂ ϕ m , n ( + ) ∂ t = − 1 2 ϕ m + 1 , n ( + ) + ϕ m − 1 , n ( + ) + ϕ m , n + 1 ( + ) + ϕ m , n − 1 ( + ) − 4 ϕ m , n ( + ) − ( | ϕ m , n ( + ) | 2 + γ | ϕ − | 2 ) ϕ m , n ( + ) + λ ϕ m + 1 , n ( − ) − ϕ m − 1 , n ( − ) − i ϕ m , n + 1 ( − ) − ϕ m , n − 1 ( − ) , i ∂ ϕ m , n ( − ) ∂ t = − 1 2 ϕ m + 1 , n ( − ) + ϕ m − 1 , n ( − ) + ϕ m , n + 1 ( − ) + ϕ m , n − 1 ( − ) − 4 ϕ m , n ( − ) − ( | ϕ m , n ( − ) | 2 + γ | ϕ m , n ( + ) | 2 ) ϕ − + λ − ϕ m + 1 , n ( + ) − ϕ m − 1 , n ( + ) − i ϕ m , n + 1 ( + ) − ϕ m , n − 1 ( + ) .
The linearized version of this system gives rise to the following dispersion relation for the plane waves, ϕ m , n ( ± ) ∼ exp − i ω t + i p x + i q y , with wavenumbers taking values in the first Brillouin zone, 0 < p , q < 2 π :
ω = 2 sin 2 p 2 + sin 2 q 2 ± 2 λ sin 2 p + sin 2 q .
The numerical solution of Eq. (55) has produced 2D modes which are discrete counterparts of the SV and MM solitons of the continuum system (52), see examples in Figure 10. As concerns the stability, the discreteness extends the stability of the SV and MM solitons towards γ > 1 and γ < 1 , respectively.
A drastic difference of the discrete solitons of both the SV and MM types from their counterparts in the continuum is that they suddenly suffer delocalization (decay) when the SOC strength λ in Eq. (55) exceeds a certain critical value, λ cr . The dependence of λ cr on the soliton’s norm,
N = ∑ m , n ϕ m , n ( + ) 2 + ϕ m , n ( − ) 2 ,
for the SV and MM solitons is displayed in Figure 11. The onset of the delocalization may be explained as a transition of the solution from the spectral bandgap to the band populated by the small-amplitude plane-wave states in the system’s linear spectrum, which is produced by Eq. (56).

5. Stable Soliton Species in the 3D DNLS Equation

5.1. The 3D Setting

A natural development of the analysis of the solitons and solitary vortices, and their bound states, produced by the 2D discrete DNLS equation and its extensions, which is outlined above in Section 2 – 4, is to construct self-trapped states (solitons) in the framework of the 3D equation:
i ϕ ˙ l , m , n + C ϕ l + 1 , m , n + ϕ l , m + 1 , n + ϕ l , m , n + 1 + ϕ l − 1 , m , n + ϕ l , m − 1 , n − 1 + ϕ l , m , n − 1 − 6 ϕ l , m , n + ϕ l , m , n 2 ϕ l , m , n = 0 ,
where, as above, the overdot stands for the time derivative, l , m , n is the set of the 3D discrete coordinates, and C > 0 is the coefficient of the intersite coupling. The 3D DNLS equation cannot be realized in optics, but it admits natural implementation for BEC loaded in a deep 3D OL potential [97,98]. In that case, ϕ l , m , n ( t ) is the the respective effectively discretized BEC wave function.
As above, stationary soliton solutions to Eq. (58) with chemical potential ω are looked for as
ϕ l , m , n = exp ( − i ω t ) u l , m , n ,
where the stationary discrete wave function u l , m , n obeys the corresponding equation,
ω u l , m , n + C u l + 1 , m , n + u l , m + 1 , n + u l , m , n + 1 + u l − 1 , m , n + u l , m − 1 , n − 1 + u l , m , n − 1 − 6 u l , m , n + u l , m , n 2 u l , m , n = 0 .
In particular, numerical solutions of Eq. (60) for 3D discrete solitons with embedded vorticity S = 0 , 1 , 2 , . . . ( S = 0 corresponds to the fundamental solitons, for which the wave function u l , m , n is real) can be obtained, starting from the natural input
u l , m , n ( init ) = A l + i m S sech η l 2 + m 2 exp − | n | ,
where η is a real scale parameter, and it is implied that the vorticity axis is directed along coordinate n [119].
It is also relevant to consider a two-component system of nonlinearly-coupled 3D DNLS equations, for wave functions ϕ l , m , n ( t ) and ψ l . m . n ( t ) of two interacting BEC species (most typically, these are different hyperfine states of the same atom) [119]:
i ϕ ˙ l , m , n + C ϕ l + 1 , m , n + ϕ l , m + 1 , n + ϕ l , m , n + 1 + ϕ l − 1 , m , n + ϕ l , m − 1 , n − 1 + ϕ l , m , n − 1 − 6 ϕ l , m , n + ϕ l , m , n 2 + β | ψ l , m , n | 2 ϕ l , m , n = 0 , i ψ ˙ l , m , n + C ψ l + 1 , m , n + ψ l , m + 1 , n + ψ l , m , n + 1 + ψ l − 1 , m , n + ψ l , m − 1 , n − 1 + ψ l , m , n − 1 − 6 ψ l , m , n + ψ l , m , n 2 + β | ϕ l , m , n | 2 ψ l , m , n = 0 ,
Here β > 0 is the relative strength of the inter-component attractive interaction with respect to the intra-component self-attraction.

5.2. Results

5.2.1. Single-Component 3D Solitons

The numerical analysis, starting from input (61), has provided families of fundamental and vortex solitons. Here, following Ref. [119], the results are displayed for a fixed value of the chemical potential, ω = − 2 in Eq. (59), while varying coupling constant C in Eqs. (58) and (60). In particular, the discrete fundamental solitons with S = 0 are stable at C < C cr ( 0 ) ≈ 2 , and the vortex modes with S = 1 are stable at C < C cr ( 1 ) ≈ 0.65 . Note that the limit of C → ∞ corresponds to the 3D NLS equation in the continuum, in which all solitons are definitely unstable; therefore, all discrete solitons become unstable at sufficiently large values of C. At C > C cr ( 1 ) the simulations demonstrate that the development of the instability destroys the vortical structure and, eventually, transforms the soliton into a fundamental one, with S = 0 (not shown here in detail). The change of the topological charge is possible, as the angular momentum is not a dynamical invariant of the lattice dynamics.
The vortex solitons with S = 2 are completely unstable, but an unusual feature of these states is that, at sufficiently small values of C (in particular, at C = 0.01 ), the instability spontaneously rebuilds them into stable discrete solitons with a larger vorticity, S = 3 [119]. An example of a stable soliton with S = 3 is displayed in Figure 12.
In addition to the fundamental and vortex solitons, Eqs. (58) and (60) produce diverse multimode species of stable discrete 3D solitons in the form of dipoles, quadrupoles and octupoles [120]. Examples of such states are presented in Figure 13 for C = 0.1 . This figure displays tightly-bound dipoles with different orientations with respect to the lattice, v i z ., straight, 2D-diagonal, and 3D-diagonal ones (panels (a,b,c)), quadrupole (panel (d)), and octupole (panel (f)), in which the field fills adjacent sites of the lattice (with the lattice distance between them d = 1 ). Also displayed are loosely-bound quadrupole and octupole (in panels (e) and (g), respectively), with distance d = 2 between the filled sites. Similar multimode states with still larger separations d between the filled sites were found too. The results are summarized in Figure 13(h), which shows stability boundaries C c r ( 3 D , d ) for dipoles, quadrupoles, and octupoles vs. d. Naturally, the stability region, C < C c r ( 3 D , d ) , increases with the increase of d, as the interaction between the filled sites, which leads to the possible dynamical instability of the multipole states, is weaker for larger d.
In addition to the above-mentioned states, Eq. (60) admits more sophisticated stable composite states, such as “vortex cubes", built as a pair of identical parallel quasi-planar vortices with topological numbers S 1 = S 2 = 1 , with opposite signs (phase shift π ). set in parallel planes, as shown in Figure 14(a). Stationary solutions representing vortex-antivortex cubes, in the form of parallel quasi-planar vortices with opposite topological charges, S 1 = − S 2 = 1 , can be found too, as shown in Figure 14(b), but they are completely unstable.
The same Eq. (58) gives rise to other stable self-trapped modes, such as vortex solitons with the axis directed along the 2D diagonal, cf. Figure 13(b). Vortex modes with the axis parallel to the 3D diagonal exist too, but they are unstable, see further details in Ref. [120].

Two-Component 3D Solitons (Including Skyrmions)

The system of coupled 3D DNLS equations (62) produces stable soliton states which are specific to the two-component nonlinear lattice medium. A noteworthy example is a composite mode built as a bound state of vortex solitons in the two components with mutually orthogonal orientations, see an example in Figure 15. These bound states are stable for sufficiently small values of the coupling constant, such as C = 0.01 in Figure 15, and for β < 1 in Eq. (62), i.e., under the condition that the self-attraction in each component is stronger than the inter-component attraction.
The system of coupled GP equations with the repulsive sign of the nonlinearity may be used as the simplest model producing skyrmions in the binary BEC [121,122,123]. The discretization of the GP system leads to Eqs. (62) with the opposite sign in front of the nonlinear terms [124]. Then, these equations are reduced to stationary ones by the usual substitution with chemical potential ω , ϕ , ψ = exp − i ω t u l , m , n , v l , m , n :
C u l + 1 , m , n + u l , m + 1 , n + u l , m , n + 1 + u l − 1 , m , n + u l , m − 1 , n − 1 + u l , m , n − 1 − 6 u l , m , n − u l , m , n 2 + β | v l , m , n | 2 u l , m , n = − ω u l , m , n , C v l + 1 , m , n + v l , m + 1 , n + v l , m , n + 1 + v l − 1 , m , n + v l , m − 1 , n − 1 + v l , m , n − 1 − 6 v l , m , n − v l , m , n 2 + β | u l , m , n | 2 v l , m , n = − ω v l , m , n ,
where the relative strength β of the inter-component repulsion with respect to the self-repulsion remain a positive coefficient. For ω > 0 , skyrmions can be constructed by choosing field u l , m , n as a complex one, representing a quasi-flat vortex soliton with topological charge S = 1 , and real field  v l , m , n as a bubble into which the vortex soliton is embedded, with a nonzero background value at | l | , | m | , | n | → ∞ , viz., v background 2 = ω [124]. An example of a numerically found skyrmion solution of this type is displayed in Figure 16.
The same work [124] reported solutions for 2D discrete “baby skyrmions", which are lattice counterparts of the modes produced by the 2D reduction of the Skyrme model [125,126]. The have a simple structure similar to its 3D counterpart displayed in Figure 16, i.e., a complex 2D vortex soliton in one component, embedded into a bubble of the delocalized field in the other real component.

6. 2D solitons and Solitary Vortices in Semi-Discrete Systems

6.1. Spatiotemporal Optical Solitons in Arrayed Waveguides

The consideration of 2D and 3D settings suggests a natural option to introduce 2D semi-discrete systems, with a continuous coordinate in one direction and a discrete coordinate in another, as well as 3D systems, where one or two coordinates are continuous, while the remaining two or one coordinates are discrete. In optics, a well-known 2D setting belonging to this class represents spatiotemporal propagation of light in an array of optical fibers [127]. It is modeled by the system of linearly-coupled NLS equations for complex amplitudes u n z , τ of optical fields in individual fibers:
i ∂ z u n + ( 1 / 2 ) D ∂ τ 2 u n + ( κ / 2 ) u n + 1 + u n − 1 − 2 u n + u n 2 u n = 0 ,
where z is the propagation distance, τ ≡ t − x / V gr (with time t and carrier group velocity V gr ) is the usual temporal variable, real D is the group-velocity-dispersion coefficient in each fiber, κ > 0 is the coefficient of coupling between adjacent fibers in the array, and the nonlinearity coefficient is normalized to be 1. It is commonly known that optical solitons (semi-discrete ones, in the present case) can be supported in the case of anomalous dispersion, i.e., D > 0 .
A remarkable counter-intuitive property of semi-discrete localized modes generated by Eq. (64) is their ability to stably move across the array, under the action of a kick applied at z = 0 [71]:
u n ( τ ) → exp i a n u n ( τ ) ,
with real a, in spite of the presence of the respective quasi-1D Peierls-Nabarro potential, An example of such a moving mode is displayed in Figure 17. This property may be compared to the above-mentioned mobility of 1D discrete solitons in the DNLS equation [65], and of 2D discrete solitons in the framework of the χ ( 2 ) system (25).
Similarly, quasi-discrete settings modeled by an extension of (64) with two transverse spatial coordinates were used for the creation for spatiotemporal optical solitons (“light bullets") [128], as well as soliton-like transient modes with embedded vorticity [129]. Waveguides employed in those experiments feature a transverse hexagonal-lattice structure, written in bulk silica by means of an optical technology. A spatiotemporal vortex state observed in the bundle-like structure (in the experiment, it is actually a transient one) is represented by Figure 18, which displays both numerically predicted and experimentally observed distributions of intensity of light in the transverse plane, together with a phase plate used in the experiment to embed the vorticity into the incident spatiotemporal pulse which was used to create the mode.

6.2. Semi-Discrete Quantum and Photonic Droplets

A new type of semi-discrete solitons was recently elaborated in Ref. [73], in the framework of an array of linearly coupled 1D GP equations, including the above-mentioned Lee-Huang-Yang correction, which represents an effect of quantum fluctuations around the mean-field states of a binary BEC [29,33]. The system is
i ∂ t ψ j = − ( 1 / 2 ) ∂ z z ψ j − C / 2 ψ j + 1 − 2 ψ j + ψ j − 1 + g | ψ j | 2 ψ j − | ψ j | ψ j ,
where ψ j ( z ) is the mean-field wave function in the j-th core with coordinate z, C is the effective inter-core coupling constant, the self-attractive quadratic term represents the Lee-Hung-Yang correction in the 1D limit, cf. Eq. (7), and g > 0 accounts for the mean-field self-repulsion.
A semi-discrete system similar to the one modeled by Eq. (66), but with the cubic-quintic nonlinearity instead of the combination of the quadratic and cubic terms in Eq. (66), was derived in the context of nonlinear optics [74]:
i ∂ z u n = − ( 1 / 2 ) ∂ x x u n − ( C / 2 ) u n + 1 − 2 u n + u n − 1 − | u n | 2 u n + | u n | 4 u n .
It corresponds to the array of parallel-coupled planar waveguides, as shown in Figure 19. In this case, u n ( x , z ) is the complex local amplitude of the optical wave in the n-th waveguide, z is the propagation distance, and x is the transverse coordinate in each waveguide, while C is the effective coupling constant, similar to Eq. (66). By analogy with the quantum droplets, semi-discrete solitons produced by Eq. (67) may be called “photonic droplets".
The droplets produced by Eqs. (66) and (67) are characterized by the total norm, which is proportional to the number of atoms in BEC,
N = ∑ j ∫ − ∞ + ∞ ψ j ( z ) 2 d z ,
or the total power of the photonic droplet,
P = ∑ n ∫ − ∞ + ∞ u n ( x ) 2 d x .
For solitons produced by Eqs. (66) and (67) sets of control parameters are, respectively, C , g for fixed N, or C , P .
The models based on Eqs. (66) and (67) give rise to many families of semi-discrete solitons, including a novel species of semi-discrete vortex solitons. Typical examples of the onsite- and intersite vortices with topological charge S = 1 , produced by Eq. (66), are displayed in Figure 20.
An example of a stable semi-discrete vortex soliton produced by Eq. (67) with S = 2 is displayed in Figure 21.
Getting back to the semi-discrete system (66), a chart in the plane of C , g which displays stability areas for the semi-discrete vortex solitons with topological charges S = 2 , 3, 4, and 5, is plotted in Figure 22. The chart demonstrates abundant multistability – for instance, the stable solitons with S = 2 coexist with all higher-order ones (with S = 3 , 4, and 5).
Self-trapped solutions of a continuum model, which are similar to semi-discrete vortex solitons outlined above, were recently reported for a photonic crystal built in a χ ( 2 ) material with a checkerboard structure representing quasi-phase matching [75].
Semi-discreteness of another type is possible in two-component systems, where one component is governed by a discrete equation, and the other one by a continuous equation. This type of two-component systems was proposed in [76]. It introduced a χ ( 2 ) model, assuming that the continuous second-harmonic wave propagates in a slab with a continuous transverse coordinate, while the fundamental-harmonic field is trapped in a discrete waveguiding array built on top of the slab.

7. 2D Fundamental and Vortical Discrete Solitons in a Two-Component PT (parity-time) symmetric lattice

While the above presentation deals solely with conservative discrete systems, many properties of conservative settings are shared by a very special type of dissipative ones, viz., systems with the parity-time ( PT ) symmetry. They include mutually symmetric spatially separated elements carrying linear gain and loss [138,139,140]. The experimental realization of such systems in optics [140] suggests one to include the Kerr nonlinearity, thus opening the way to the prediction an creation of PT -symmetric solitons [141,142]. In particular, exact solutions for 1D PT -symmetric solitons and exact results for their stability boundaries were found in the model of the nonlinear PT -symmetric coupler (dual-core waveguide), with mutually symmetric linear gain and loss carried by the linearly coupled cores [143,144]. Stability limits for 2D fundamental solitons in the 2D PT -symmetric coupler with the cubic-quintic nonlinearity in each core (essentially the same as in Eqs. (67), chosen to prevent the critical-collapse instability) were identified in Ref. [145].
The definition of the PT symmetry makes it also natural to consider discrete PT -symmetric systems. Various species of stable discrete solitons were predicted in chains of PT -symmetric elements [83,146,147,148,149,150], and the existence of such solitons was demonstrated experimentally [151].
A natural model for the creation of PT -symmetric discrete 2D solitons is a generalization of the 2D discrete nonlinear coupler, based on Eqs. (44), by adding the linear gain and loss terms with strength γ > 0 to the coupled equations [83]:
i d ψ m , n d z = − C 2 ( ψ m , n + 1 + ψ m , n − 1 + ψ m − 1 , n + ψ m + 1 , n − 4 ψ m , n ) − | ψ m , n | 2 ψ m , n − φ m , n + i γ ψ m , n , i d φ m , n d z = − C 2 ( φ m , n + 1 + φ m , n − 1 + φ m − 1 , n + φ m + 1 , n − 4 φ m , n ) − | φ m , n | 2 φ m , n − ψ m , n − i γ φ m , n .
Here, in terms of the optical realization, the evolution variable z is the propagation distance, the inter-core coupling coefficient is scaled to be 1, and C > 0 is constant of the intra-core coupling between adjacent sites of the lattice. The dispersion relations for plane-wave solutions to the linearized version of Eqs. (70), ψ m , n ( z ) , φ m , n ( z ) ∼ exp i p m + i q n + i k z , is
k = − 2 C sin 2 p 2 + sin 2 q 2 ± 1 − γ 2 .
As it follows from Eq. (71), the PT symmetry holds under condition γ < γ max ≡ 1 , i.e., the gain-loss strength γ must be smaller than the linear-coupling coefficient, that is 1 in the present notation, which is a generic property of PT -symmetric couplers [143,144].
Stationary modes with real propagation constant k are looked for as solutions to the full nonlinear system of Eqs. (70) in the usual form, ψ m , n ( z ) , φ m , n ( z ) = e i k z u m , n , v m , n , with stationary amplitudes obeying equations
k + i γ u m , n = C 2 ( u m , n + 1 + u m , n − 1 + u m − 1 , n + u m + 1 , n − 4 u m , n ) + | u m , n | 2 u m , n + v m , n , k − i γ v m , n = C 2 ( v m , n + 1 + v m , n − 1 + v m − 1 , n + v m + 1 , n − 4 v m , n ) + | v m , n | 2 v m , n + u m , n .
Localized states produced by Eqs. (72) are characterized, as above, by the total power,
P = ∑ m , n | u m , n | 2 + | v m , n | 2 .
Straightforward analysis of Eqs. (72) demonstrates that the system produces PT -symmetric fundamental-soliton solutions, which must be subject to the relation v m . n = u m . n * (with * standing for the complex conjugate), in the form of
u m . n , v m n = w m , n exp ± ( i / 2 ) arcsin γ ,
where real discrete distribution w m , n should be found as a solution of the usual stationary equation for 2D discrete solitons,
− k + 1 − γ 2 w m , n + C 2 w m , n + 1 + w m , n − 1 + w m − 1 , n + w m + 1 , n − 4 w m , n + w m , n 3 = 0 ,
cf. Eq. (38). In agreement with the linear spectrum (71), Eq. (75) may produce soliton solutions for k > 1 − γ 2 . An example of a stable fundamental PT -symmetric soliton is displayed, by means of its cross-section shapes, in Figure 23.
The existence and stability of the PT -symmetric fundamental discrete solitons is summarized in the plane of γ , P for the fundamental solitons in Figure 24. It is seen that, naturally, the stability area shrinks as the gain-loss coefficient γ is approaching its limit value, γ max = 1 (cf. the 1D situation considered in Refs. [143] and [144]). The existence boundary, i.e., the minimum value of P, below which no solitons are found (in the white area), corresponds to the limit of very broad small-amplitude solitons. In this limit, the discrete soliton may be approximated by its counterpart in the continuum NLS equation, i.e., the above-mentioned Townes soliton, whose power takes the unique value, which thus determines the existence boundary in Figure 24.
The stability boundary in Figure 24 may be understood as the one at which the symmetric soliton is destabilized by the spontaneous symmetry breaking (as described in detail above for 2D solitons produced by the linearly-coupled conservative DNLS equations (44), see also Ref. [17]), which is here modified by the presence of the linear gain and loss. Because asymmetric solitons cannot exist in the system with the balanced gain and loss, the symmetry breaking always leads to either blowup or decay of the soliton [83]. In their stability region, the PT -symmetric fundamental discrete solitons actually represent the system’s ground state [83].
Alongside the fundamental discrete PT -symmetric solitons, the same system of Eqs. (72) produces PT -symmetric vortex solitons, which also have their stability area, see details in Ref. [83]. An example of a stable PT -symmetric vortex soliton is presented in Figure 25.
In addition to the PT -symmetric solitons, Eqs. (72) give rise to anti- PT -symmetric ones, defined by relation v m . n = − u m . n * . They, as well as anti- PT -symmetric vortex solitons, are stable in some parameter areas (see details in Ref. [83]), but those areas are essentially smaller than their counterparts for the PT -symmetric modes. The reduced stability area for the anti- PT -symmetric fundamental solitons is explained by the fact that they cannot be the system’s ground state.

8. Conclusion

The interplay of the discreteness and intrinsic nonlinearity in various physical media – chiefly, in nonlinear optics and BEC – gives rise to a great variety of self-trapped localized states, in the form of discrete solitons. This article aims to produce a concise review, starting from the brief survey of basic theoretical models combining the discreteness in 1D, 2D, and 3D geometries and various nonlinearities, such as cubic, quadratic, and quintic. The main subject addressed in the article is a summary of basic results for 2D and 3D discrete solitons produced by such models. Unlike the topic of 1D discrete solitons, the multidimensional ones were not previously reviewed in a systematic form. Along with the fundamental solitons, topologically organized ones, in the form of solitary vortices and discrete skyrmions, are considered too. Some experimental findings are also included, such as the observation of 2D discrete optical solitons with embedded vorticity.
In many cases, the discreteness helps to produce states which either do not exist or are definitely unstable in continuum analogs of the discrete settings. In particular, these are 2D fundamental and vortex solitons, which may be stable in the discrete form, while their continuum counterparts are completely unstable in the free space. On the other hand, mobility of solitons, which is their obvious property in the continuum, is a nontrivial issue for the lattice (discrete) solitons.
The work in this area remains a subject of ongoing theoretical and experimental work, promising new findings. A perspective direction is to produce 2D and 3D self-trapped states with intrinsic topological structures. Some results obtained in this direction are presented in this article, such as discrete solitons in the system with spin-orbit coupling [118] (see also Ref. [130]), sophisticated 3D discrete modes with embedded vorticity [119,120], and discrete skyrmions [124]. A challenging task is experimental realization of these states which, thus far, were only predicted in the theoretical form.
It is relevant to mention some topics which may be relevant in the present context but are not included here, to keep a reasonable size of the review. In particular, these are interactions of discrete solitons with local defects in the underlying lattice, as well as with interfaces and edges. It is known that defects and surfaces may often help to create and stabilize localized modes which do not exist or are unstable in uniform lattices, such as Tamm [131] and topological-insulator [132,133] states. Another vast area of studies, which is not considered here, deals with dissipative discrete nonlinear systems. In this article, only the very special case of PT -symmetric systems is addressed. Basic nonlinear dissipative models are represented by discrete complex Ginzburg-Landau equations, i.e., DNLS equations with complex coefficients in front of the onsite linear and nonlinear terms, which account for losses and compensating gain [134]. Unlike conservative and PT -symmetric models, the dissipative ones may only give rise to stable discrete solitons which do not exist in continuous families, but rather as isolated attractors [135,136,137].

Funding

This work was supported, in part, by the Israel Science Foundation, through grant No. 1695/22.

Acknowledgments

I appreciate the invitation of Editors of the special issue of Entropy, Profs. Lars English and Faustino Palmero, to submit this article as a contribution to the special issue on the topic of “Recent Advances in the Theory of Nonlinear Lattices". I would like to thank colleagues in collaboration with whom I have been working on various topics related to the review: G. E. Astrakharchik, B. B. Baizakov, P. Beličev, A. R. Bishop, R. Blit, L. L. Bonilla, O. V. Borovkova, R. Carretero-González, Zhaopin Chen, Zhigang Chen, C. Chong, J. Cuevas-Maraver, J. D’Ambroise, F. K. Diakonos, S. V. Dmitriev, R. Driben, N. Dror, O. Dutta, A. Eckhard, L. M. Floría, D. J. Frantzeskakis, S. Fu, G. Gligorić, J. Gómez-Gardeñes, L. Hadžievsli, D. Herring, J. Hietarinta, P. Hauke, T. Kapitula, Y. V. Kartashov, D. J. Kaup (deceased). P. G. Kevrekidis, V. V. Konotop, M. Lewenstein, Ben Li, Yongyao Li, A. Maluckov, T. Meier, D. Mihalache, N. C. Panoiu, I. E. Papachalarampus, J. Petrović, M. A. Porter, K. . Rasmussen, H. Sakaguchi, M. Salerno, B. Sánchez-Rey, H. Susanto, L. Torner, M. Trippenbach, A. V. Ustinov, R. A. Van Gorder, A. Vardi, M. I. Weinstein, X. Xu, A. V. Yulin, and K. Zegadlo.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. L. P. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford University Press, Oxford, 2003).
  2. Yu. S. Kivshar, G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, San Diego, 2003).
  3. G. Fibich, The Nonlinear Schrödinger Equation: Singular Solutions and Optical Collapse (Springer, Heidelberg, 2015).
  4. O. Morsch and M. Oberthaler, Dynamics of Bose-Einstein condensates in optical lattice, Rev. Mod. Phys. 78, 179-215 (2006). [CrossRef]
  5. M. A. Porter, R. Carretero-González, P. G. Kevrekidis, and B.A. Malomed, Nonlinear lattice dynamics of Bose-Einstein condensates, Chaos 15, 015115 (2005). [CrossRef]
  6. M. Skorobogatiy, J. Yang, Fundamentals of Photonic Crystal Guiding (Cambridge University Press, Cambridge, 2008).
  7. D. N. Christodoulides and R. I. Joseph, Discrete self-focusing in nonlinear arrays of coupled waveguides, Opt. Lett. 13, 794-796 (1988). [CrossRef]
  8. H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, Discrete spatial optical solitons in waveguide arrays, Phys. Rev. Lett. 81, 3383-3386 (1998). [CrossRef]
  9. D. N. Christodoulides, F. Lederer, and Y. Silberberg, Discretizing light behaviour in linear and nonlinear waveguide lattices, Nature 424, 817-823 (2003). [CrossRef]
  10. F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberg, Discrete solitons in optics, Phys. Rep. 463, 1-126 (2008). [CrossRef]
  11. F. Ye, D. Mihalache, B. Hu, N. C. Panoiu, Subwavelength plasmonic lattice solitons in arrays of metallic nanowires, Phys. Rev. Lett. 104, 106802 (2010). [CrossRef]
  12. A. Smerzi and A. Trombettoni, Nonlinear tight-binding approximation for Bose-Einstein condensates in a lattice, Phys. Rev. A 68, 023613 (2003). [CrossRef]
  13. A. Szameit, R. Keil, F. Dreisow, M. Heinrich, T. Pertsch, S. Nolte, and A. Tünnermann, Observation of discrete solitons in lattices with second-order interaction, Opt. Lett. 34, 2838-2840 (2009). [CrossRef]
  14. C. Chong, R. Carretero-González, B. A. Malomed, and P. G. Kevrekidis, Variational approximations in discrete nonlinear Schrödinger equations with next-nearest-neighbor couplings, Physica D 240, 1205-1212 (2011). [CrossRef]
  15. A. Szameit, T. Pertsch, S. Nolte, A. Tünnermann, F. Lederer, Long-range interaction in waveguide lattices, Phys. Rev. A 77, 043804 (2008). [CrossRef]
  16. A. Locatelli, D. Modotto, D. Paloschi, and C. De Angelis, All optical switching in ultrashort photonic crystal couplers, Opt. Commun. 237, 97-102 (2004). [CrossRef]
  17. G. Herring, P. G. Kevrekidis, B. A. Malomed, R. Carretero-González, D. J. Frantzeskakis, Symmetry breaking in linearly coupled dynamical lattices, Phys. Rev. E 76, 066606 (2007). [CrossRef]
  18. S. Aubry, Breathers in nonlinear lattices: Existence, linear stability and quantization, Physica D 103, 201-250 (1997). [CrossRef]
  19. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, Soliton shape and mobility control in optical lattices, Prog. in Optics 52, 63-148 (2009). [CrossRef]
  20. V. M. Rothos, Nonlinear wave propagation in discrete and continuous systems, Eur. Phys. J. Special Topics 225, 943–958 (2016). [CrossRef]
  21. E. N. Tsoy and B. A. Umarov, Introduction to nonlinear discrete systems: theory and modelling, Eur. J. Phys. 39, 055803 (2018). [CrossRef]
  22. B. A. Malomed, Nonlinearity and discreteness: Solitons in lattices, in Emerging Frontiers in Nonlinear Science, pp. 81-110, ed. by P. G. Kevrekidis, J. Cuevas-Maraver, and A. Saxena (Springer Nature Switzerland AG: Cham, 2020).
  23. P. G. Kevrekidis, The Discrete Nonlinear Schrödinger Equation: Mathematical Analysis, Numerical Computations, and Physical Perspectives (Springer, Berlin Heidelberg, 2009).
  24. E. W. Laedke, K. H. Spatschek, and S. K. Turitsyn, Stability of discrete solitons and quasicollapse to intrinsically localized modes, Phys. Rev. Lett. 73, 1055-1059 (1994). [CrossRef]
  25. Y. Kartashov, G. Astrakharchik, B. Malomed, and L. Torner, Frontiers in multidimensional self-trapping of nonlinear fields and matter, Nature Reviews Physics 1, 185-197 (2019). [CrossRef]
  26. B. A. Malomed, Multidimensional Solitons (AIP Publishing, Melville, NY, 2022).
  27. B. A. Malomed, Multidimensional Soliton Systems, Advances in Physics: X, in press.
  28. T. D. Lee, K. Huang, and C. N. Yang, Eigenvalues and eigenfunctions of a Bose system of hard spheres and its low temperature properties, Phys. Rev. 106, 1135-1145 (1957). [CrossRef]
  29. D. S. Petrov, Quantum mechanical stabilization of a collapsing Bose-Bose mixture, Phys. Rev. Lett. 115, 155302 (2015). [CrossRef]
  30. C. Cabrera, L. Tanzi, J. Sanz, B. Naylor, P. Thomas, P. Cheiney, and L. Tarruell, Quantum liquid droplets in a mixture of Bose-Einstein condensates, Science 359, 301-304 (2018). [CrossRef]
  31. P. Cheiney, C. R. Cabrera, J. Sanz, B. Naylor, L. Tanzi, and L. Tarruell, Bright soliton to quantum droplet transition in a mixture of Bose-Einstein condensates, Phys. Rev. Lett. 120, 135301 (2018). [CrossRef]
  32. G. Semeghini, G. Ferioli, L. Masi, C. Mazzinghi, L. Wolswijk, F. Minardi, M. Modugno, G. Modugno, M. Inguscio, and M. Fattori, Self-bound quantum droplets of atomic mixtures in free space?, Phys. Rev. Lett. 120, 235301 (2018). [CrossRef]
  33. D. S. Petrov and G. E. Astrakharchik, Ultradilute low-dimensional liquids, Phys. Rev. Lett. 117, 100401 (2016). [CrossRef]
  34. G. E. Astrakharchik and B. A. Malomed, Dynamics of one-dimensional quantum droplets, Phys. Rev. A 98, 01363 (2018). [CrossRef]
  35. V. E. Zakharov, S. V. Manakov, S. P. Novikov, and L. P. Pitaevskii, Solitons: The Inverse Scattering Method (Nauka Publishers, Moscow, 1980) [English translation: Consultants Bureau, New York, 1984].
  36. M. J. Ablowitz and H. Segur, Solitons and Inverse Scattering Method (SIAM, Philadelphia, 1981).
  37. F. Calogero and A. Degasperis, Spectral Transform and Solitons: Tools to Solve and Investigate Nonlinear Evolution Equations (North-Holland, New York, 1982).
  38. A. C. Newell, Solitons in Mathematics and Physics (SIAM: Philadelphia, 1985).
  39. M. J. Ablowitz and B. M. Herbst, On homoclinic structure and numerically induced chaos for the nonlinear Schrödinger equation, SIAM J. Appl. Math. 50, 339-351 (1990). [CrossRef]
  40. D. Levi, M. Petrera, and C. Scimiterna, On the integrability of the discrete nonlinear Schrödinger equation, Europhys. Lett. 84, 10003 (2008). [CrossRef]
  41. M. J. Ablowitz and J. F. Ladik, Nonlinear differential–difference equations and Fourier analysis, J. Math. Phys. 17, 1011-1018 (1976). [CrossRef]
  42. Yu. B. Suris, The problem of integrable discretization: Hamiltonian approach (Birkhauser, Basel, 2003).
  43. M. Salerno, A new method to solve the quantum Ablowitz-Ladik system, Phys. Lett. A 162, 381-384 (1992). [CrossRef]
  44. O. Dutta, M. Gajda, P. Hauke, M. Lewenstein, D.-S. Luhmann, B.A. Malomed, T. Sowinski, J. Zakrzewski, Non-standard Hubbard models in optical lattices: a review, Rep. Prog. Phys. 78, 066001 (2015). [CrossRef]
  45. O. V. Borovkova, Y. V. Kartashov, L. Torner, and B. A. Malomed, Bright solitons from defocusing nonlinearities, Phys. Rev. E 84, 035602 (R) (2011). [CrossRef]
  46. G. Gligorić, A. Maluckov, L. Hadzievski, and B. A. Malomed, Discrete localized modes supported by an inhomogeneous defocusing nonlinearity, Phys. Rev. E 88, 032905 (2013). [CrossRef]
  47. P. G. Kevrekidis, B. A. Malomed, A. Saxena, A. R. Bishop, and D. J. Frantzeskakis, Solitons and vortices in two-dimensional discrete nonlinear Schrödinger systems with spatially modulated nonlinearity, Phys. Rev. E 91, 043201 (2015). [CrossRef]
  48. T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. Pfau, The physics of dipolar bosonic quantum gases, Rep. Prog. Phys. 72, 126401 (2009). [CrossRef]
  49. G. Gligorić, A. Maluckov, M. Stepić, L. Hadžievski, and B. A. Malomed, Discrete vortex solitons in dipolar Bose-Einstein condensates, J. Phys. B: At. Mol. Opt. Phys. 43, 055303 (2010). [CrossRef]
  50. P. Pedri and L. Santos, Two-dimensional bright solitons in dipolar Bose-Einstein condensates, Phys. Rev. Lett. 95, 200404 (2005). [CrossRef]
  51. I. Tikhonenkov, B. A. Malomed, and A. Vardi, Anisotropic solitons in dipolar Bose-Einstein condensates,. Phys. Rev. Lett. 100, 090406 (2008). [CrossRef]
  52. Y. Li, J. Liu, W. Pang, and B. A. Malomed, Lattice solitons with quadrupolar intersite interactions, Phys. Rev. A 88, 063635 (2013). [CrossRef]
  53. V. Buryak, P. Di Trapani, D. V. Skryabin, and S. Trillo, Optical solitons due to quadratic nonlinearities: from basic physics to futuristic applications, Phys. Rep. 370, 63-235 (2002).
  54. H. Susanto, P. G. Kevrekidis, R. Carretero-González, B. A. Malomed, and D. J. Frantzeskakis, Mobility of discrete solitons in quadratically nonlinear media, Phys. Rev. Lett. 99, 214103 (2007). [CrossRef]
  55. J. Gómez-Gardeñes, B. A. Malomed, L. M. Floría, A. R. Bishop, Solitons in the Salerno model with competing nonlinearities, Phys. Rev. E 73, 036608 (2006). [CrossRef]
  56. D. Chen, S. Aubry, G. P. Tsironis, Breather mobility in discrete φ4 nonlinear lattices, Phys. Rev. Lett. 77, 4776-4779 (1996). [CrossRef]
  57. Yu. S. Kivshar and M. Peyrard, Modulational instabilities in discrete lattices, Phys. Rev. A 46, 3198-3205 (1992). [CrossRef]
  58. B. A. Malomed, Variational methods in nonlinear fiber optics and related fields, Prog. Optics 43, 71-193 (2002). [CrossRef]
  59. B.A. Malomed and M. I. Weinstein, Soliton dynamics in the discrete nonlinear Schrödinger equation, Phys. Lett. A 220, 91-96 (1996). [CrossRef]
  60. I. E. Papacharalampous, P. G. Kevrekidis, B. A. Malomed, D. J. Frantzeskakis, Soliton collisions in the discrete nonlinear Schrödinger equation, Phys. Rev. E 68, 046604 (2003). [CrossRef]
  61. D. J. Kaup, Variational solutions for the discrete nonlinear Schrödinger equation, Math. Comput. Simulat. 69, 322-333 (2005). [CrossRef]
  62. B. A. Malomed, D. J. Kaup, and R. A. Van Gorder, Unstaggered-staggered solitons in two-component discrete nonlinear Schrödinger lattices, Phys. Rev. E 85, 026604 (2012). [CrossRef]
  63. J. Cuevas, G. James, P. G. Kevrekidis, B. A. Malomed, B. Sánchez-Rey, Approximation of solitons in the discrete NLS equation, J. Nonlinear Math. Phys. 15 (Suppl. 3), 124-136 (2008). [CrossRef]
  64. C. Chong, D. E. Pelinovsky, and G. Schneider, On the validity of the variational approximation in discrete nonlinear Schrödinger equations, Physica D 241, 115-124 (2012). [CrossRef]
  65. D. B. Duncan, J. C. Eilbeck, H. Feddersen, and A. D. Wattis, Solitons on lattices, Physica D 68, 1-11 (1993). [CrossRef]
  66. M. Toda, Vibration of a chain with a non-linear interaction, J. Phys. Soc. Jpn. 22, 431-436 (1967). [CrossRef]
  67. S. Darmanyan, A. Kobyakov, and F. Lederer, Stability of strongly localized excitations in discrete media with cubic nonlinearity, J. Exp. Theor. Phys. 86, 682-686 (1998). [CrossRef]
  68. T. Kapitula, P. G. Kevrekidis, and B. A. Malomed, Stability of multiple pulses in discrete systems, Phys. Rev. E 63, 036604 (2001). [CrossRef]
  69. D. E. Pelinovsky, P. G. Kevrekidis, and D. J. Frantzeskakis, Stability of discrete solitons in nonlinear Schrödinger lattices, Physica D 212, 1-19 (2005). [CrossRef]
  70. P. G. Kevrekidis, B. A. Malomed, and A. R. Bishop, Bound states of two-dimensional solitons in the discrete nonlinear Schrödinger equation, J. Phys. A: Math. Gen. 34, 9615-9629 (2001). [CrossRef]
  71. R. Blit and B. A. Malomed, Propagation and collisions of semi-discrete solitons in arrayed and stacked waveguides, Phys. Rev. A 86, 043841 (2012). [CrossRef]
  72. R. Driben, V. V. Konotop, B. A. Malomed, T. Meier, and A. V. Yulin, Nonlinearity-induced localization in a periodically driven semidiscrete system, Phys. Rev. E 97, 062210 (2018). [CrossRef]
  73. X. Zhang, X. Xu, Y. Zheng, Z. Chen, B. Liu, C. Huang, B. A. Malomed, and Y. Li, Semidiscrete quantum droplets and vortices, Phys. Rev. Lett. 123, 133901 (2019). [CrossRef]
  74. X. Xu, G. Ou, Z. Chen, B. Liu, B. A. Malomed, and Y. Li, Semidiscrete vortex solitons, Advanced Photonics Research 2, 2000082 (2021). [CrossRef]
  75. X. Xu, F. Zhao, J. Huang, H. Xiang, L. Zhang, Z. Chen, Z. Nie, B. A Malomed, and Y. Li, Semidiscrete optical vortex droplets in quasi-phase-matched photonic crystals, Opt. Exp. 31, 38343-38354 (2023). [CrossRef]
  76. N. C. Panoiu, R. M. Osgood, and B. A. Malomed, Semi-discrete composite solitons in arrays of quadratically nonlinear waveguides, Opt. Lett. 31, 1097-1099 (2006). [CrossRef]
  77. X. Xu, G. Ou, Z. Chen, B. Liu, B. A. Malomed, and Y. Li, Semidiscrete vortex solitons, Advanced Photonics Research 2, 2000082 (2021). [CrossRef]
  78. X. Xu, F. Zhao, J. Huang, H. Xiang, L. Zhang, Z. Chen, Z. Nie, B. A Malomed, and Y. Li, Semidiscrete optical vortex droplets in quasi-phase-matched photonic crystals, Opt. Exp. 31, 38343-38354 (2023). [CrossRef]
  79. D. N. Neshev, T. J. Alexander, E. A. Ostrovskaya, Yu. S. Kivshar, H. Martin, I. Makasyuk, and Z. G. Chen, Observation of discrete vortex solitons in optically induced photonic lattices, Phys. Rev. Lett. 92, 123903 (2004). [CrossRef]
  80. J. W. Fleischer, G. Bartal, O. Cohen, O. Manela, M. Segev, J. Hudock, and D. N. Christodoulides, Observation of vortex-ring “discrete” solitons in 2D photonic lattices, Phys. Rev. Lett. 92, 123904 (2004). [CrossRef]
  81. J. Ford, The Fermi-Pasta-Ulam problem – paradox turns discovery, Phys. Rep. 213, 271-310 (1992). [CrossRef]
  82. O. M. Braun and Yu. S. Kivshar, The Frenkel-Kontorova Model: Concepts, Methods, and Applications (Springer, Berlin, 2004).
  83. Z. Chen, J. Liu, S. Fu, Y. Li, and B. A. Malomed, Discrete solitons and vortices on two-dimensional lattices of PT-symmetric couplers, Opt. Exp. 22, 29679-29692 (2014). [CrossRef]
  84. M. I. Weinstein, Excitation thresholds for nonlinear localized modes on lattices, Nonlinearity 12, 673-621 (1999). [CrossRef]
  85. C. Chong, R. Carretero-González, B. A. Malomed, and P. G. Kevrekidis, Multistable solitons in higher-dimensional cubic-quintic nonlinear Schrödinger lattices, Physica D 238, 126-136 (2009). [CrossRef]
  86. B. A. Malomed and P. G. Kevrekidis, Discrete vortex solitons, Phys. Rev. E 64, 026601 (2001). [CrossRef]
  87. J. Cuevas, G. James, P. G. Kevrekidis, and K. J. H. Law, Vortex solutions of the discrete Gross-Pitaevskii equation starting from the anti-continuum limit, Physica D 238, 1422-1431 (2009). [CrossRef]
  88. R. Y. Chiao, E. Garmire, and C. H. Townes, Self-Trapping of Optical Beams, Phys. Rev. Lett. 13, 479-482 (1964). [CrossRef]
  89. V. I. Kruglov, Yu. A. Logvin, and V. M. Volkov, The theory of spiral laser beams in nonlinear media, J. Mod. Opt. 39, 2277-2291 (1992). [CrossRef]
  90. B. A. Malomed, (INVITED) Vortex solitons: Old results and new perspectives, Physica D 399, 108-137 (2019). [CrossRef]
  91. P. G. Kevrekidis, B. A. Malomed, Z. Chen, and D. J. Frantzeskakis, Stable higher-order vortices and quasivortices in the discrete nonlinear Schrödinger equation, Phys. Rev. E 70, 056612 (2004). [CrossRef]
  92. Z. Chen, M. Segev, D. W. Wilson, R. E. Muller, P. D. Maker, Self-trapping of an optical vortex by use of the bulk photovoltaic effect, Phys. Rev. Lett. 78, 2948-2951 (1997). [CrossRef]
  93. Z. Chen, M.-F. Shih, M. Segev, D. W. Wilson, R. E. Muller, and P. D. Maker, Steady-state vortex-screening solitons formed in biased photorefractive media, Opt. Lett. 22, 1751-1753 (1997). [CrossRef]
  94. A. Bezryadina, E. Eugenieva, and Z. Chen, Self-trapping and flipping of double-charged vortices in optically induced photonic lattices, Opt. Lett. 31, 2456-2458 (2006). [CrossRef]
  95. B. Terhalle, T. Richter, K. J. H. Law, D. Göries, P. Rose, T. J. Alexander, P. G. Kevrekidis, Anton S. Desyatnikov, W. Krolikowski, F. Kaiser, C. Denz, and Y. S. Kivshar, Observation of double-charge discrete vortex solitons in hexagonal photonic lattices, Phys. Rev. A 79, 043821 (2009). [CrossRef]
  96. C. M. de Sterke and J. E. Sipe, Gap solitons, Prog. Opt., vol. XXXIII, pp. 203-260 (1994).
  97. V. A. Brazhnyi and V. V. Konotop, Theory of nonlinear matter waves in optical lattices, Mod. Phys. Lett. B 18, 627-651 (2004). [CrossRef]
  98. O. Morsch and M. Oberthaler, Dynamics of Bose-Einstein condensates in optical lattices, Rev. Mod. Phys. 78, 179-212 (2006). [CrossRef]
  99. B. Eiermann., Th. Anker, M. Albiez, M. Taglieber, P. Treutlein, K.-P. Marzlin, and M. K. Oberthaler, Bright Bose-Einstein gap solitons of atoms with repulsive interaction, Phys. Rev. Lett. 92, 230401 (2004). [CrossRef]
  100. J. T. Mok, C. M. de Sterke, I. C. M. Litte, and B. J. Eggleton, Dispersionless slow light using gap solitons, Nature Phys. 2, 775-780 (2006).
  101. G. Gligorić, A. Maluckov, L. Hadzievski, and B. A. Malomed, Localized modes in mini-gaps opened by periodically modulated intersite coupling in two-dimensional nonlinear lattices, Chaos 24, 023124 (2014). [CrossRef]
  102. P. G. Kevrekidis, D. J. Frantzeskakis, R. Carretero-González, B. A. Malomed, and A. R. Bishop. Discrete solitons and vortices on anisotropic lattices, Phys. Rev. E 72, 046613 (2005). [CrossRef]
  103. J. Cuevas, B. A. Malomed, and P. G. Kevrekidis, Two-dimensional discrete solitons in rotating lattices, Phys. Rev. E 76, 046608 (2007). [CrossRef]
  104. B. A. Malomed, Potential of interaction between two- and three-dimensional solitons, Phys. Rev. E 58, 7928-7933 (1998). [CrossRef]
  105. B. A. Malomed, Spontaneous symmetry breaking in nonlinear systems: An overview and a simple model, in: Nonlinear Dynamics: Materials, Theory and Experiments, ed. by M. Tlidi, M. Clerc, Springer Proceedings in Physics, vol. 173 (Springer, Cham, 2016), p. 97.
  106. G. Iooss, D. D. Joseph, Elementary Stability Bifurcation Theory (Springer, New York, 1980).
  107. D. Cai, A. R. Bishop, and N. Grønbech-Jensen, Perturbation theories of a discrete, integrable nonlinear Schrödinger equation, Phys. Rev. E 53, 4131-4136 (1996). [CrossRef]
  108. D. Cai, A. R. Bishop, and N. Grønbech-Jensen, Resonance in the collision of two discrete intrinsic localized excitations, Phys. Rev. E 56, 7246-7252 (1997). [CrossRef]
  109. S. V. Dmitriev, P. G. Kevrekidis, B. A. Malomed, and D. J. Frantzeskakis, Two-soliton collisions in a near-integrable lattice system, Phys. Rev. E 68, 056603 (2003). [CrossRef]
  110. N. G. Vakhitov and A. A. Kolokolov, Stationary solutions of the wave equation in a medium with nonlinearity saturation, Radiophys. Quantum Electron. 16, 783-789 (1973). [CrossRef]
  111. J. Gómez-Gardeñes, B. A. Malomed, L. M. Floría, and A. R. Bishop, Discrete solitons and vortices in the two-dimensional Salerno model with competing nonlinearities, Phys. Rev. E 74, 036607 (2006). [CrossRef]
  112. Y.-J. Lin, K. Jiménez-García, and I. B. Spielman, Spin-orbit-coupled Bose–Einstein condensates, Nature 471, 83–86 (2011).
  113. V. Galitski and I. B. Spielman, Spin-orbit coupling in quantum gases, Nature 494, 49-54 (2013).
  114. N. Goldman, G. Juzeliunas, P. Öhberg, and I. B. Spielman, Light-induced gauge fields for ultracold atoms, Rep. Prog. Phys. 77, 126401 (2014). [CrossRef]
  115. H. Zhai, Degenerate quantum gases with spin-orbit coupling, Rep. Prog. Phys. 78, 026001 (2015). [CrossRef]
  116. B. A. Malomed, Creating solitons by means of spin-orbit coupling, EPL 122, 36001 (2018). [CrossRef]
  117. H. Sakaguchi, B. Li, and B. A. Malomed, Creation of two-dimensional composite solitons in spin-orbit-coupled self-attractive Bose-Einstein condensates in free space, Phys. Rev. E 89, 032920 (2014). [CrossRef]
  118. H. Sakaguchi and B. A. Malomed, Discrete and continuum composite solitons in Bose-Einstein condensates with the Rashba spin-orbit coupling in one and two dimensions, Phys. Rev. E 90, 062922 (2014). [CrossRef]
  119. P. G. Kevrekidis, B. A. Malomed, D. J. Frantzeskakis and R. Carretero-González, Three-dimensional solitary waves and vortices in a discrete nonlinear Schrödinger lattice, Phys. Rev. Lett. 93, 080403 (2004). [CrossRef]
  120. R. Carretero-González, P. G. Kevrekidis, B. A. Malomed and D. J. Frantzeskakis, Three-dimensional nonlinear lattices: From oblique vortices and octupoles to discrete diamonds and vortex cubes, Phys. Rev. Lett. 94, 203901 (2005). [CrossRef]
  121. J. Ruostekoski and J. R. Anglin, Creating vortex rings and three-dimensional skyrmions in Bose-Einstein condensates, Phys. Rev. Lett. 86, 3934-3937 (2001). [CrossRef]
  122. U. Al Khawaja and H. Stoof, Skyrmions in a ferromagnetic Bose-Einstein condensate, Nature 411, 918–920 (2001). [CrossRef]
  123. R. A. Battye, N. R. Cooper, and P. M. Sutcliffe, Stable skyrmions in two-component Bose-Einstein condensates, Phys. Rev. Lett. 88, 080401 (2002). [CrossRef]
  124. P. G. Kevrekidis, R. Carretero-González, D. J. Frantzeskakis, B. A. Malomed, and F. K. Diakonos. Skyrmion-like states in two- and three-dimensional dynamical lattices. Phys. Rev. E 75, 026603 (2007). [CrossRef]
  125. A. Kudryavtsevy, B. Piette, and W. J. Zakrzewski, Skyrmions and domain walls in (2+1) dimensions, Nonlinearity 11, 783-795 (1998). [CrossRef]
  126. T. Weidig, The baby Skyrme models and their multi-skyrmions, Nonlinearity 12, 1489-1503 (1999). [CrossRef]
  127. A. B. Aceves, C. De Angelis, A. M. Rubenchik, and S. K. Turitsyn, Multidimensional solitons in fiber arrays, Opt. Lett. 19, 329 (1994). [CrossRef]
  128. S. Minardi, F. Eilenberger, Y. V. Kartashov, A. Szameit, U. Röpke, J. Kobelke, K. Schuster, H. Bartelt, S. Nolte, L. Torner, F. Lederer, A. Tünnermann, and T. Pertsch, Three-dimensional light bullets in arrays of waveguides, Phys. Rev. Lett. 105, 263901 (2010). [CrossRef]
  129. F. Eilenberger, K. Prater, S. Minardi, R. Geiss, U. Röpke, J. Kobelke, K. Schuster, H. Bartelt, S. Nolte, A. Tünnermann, and T. Pertsch, Observation of discrete, vortex light bullets, Phys. Rev. X 3, 041031 (2013). [CrossRef]
  130. P. Beličev, G. Gligorić, J. Petrović, A. Maluckov, L. Hadžievski, B.A. Malomed, Composite localized modes in discretized spin-orbit-coupled Bose-Einstein condensates, J. Phys. B At. Mol. Opt. Phys. 48, 065301 (2015). [CrossRef]
  131. Yu. S. Kivshar, Nonlinear Tamm states and surface effects in periodic photonic structures, Laser Phys. Lett. 5, 703-713 (2008). [CrossRef]
  132. D. R. Gulevich, D. Yudin, D. V. Skryabin, I. V. Iorsh, and I. A. Shelykh, Edge solitons in kagome lattice, Sci. Rep. 7, 1780 (2017). [CrossRef]
  133. Y. V. Kartashov and D. V. Skryabin, Modulational instability and solitary waves in polariton topological insulators, Optica 3, 1228 (2016). [CrossRef]
  134. V. Hakim and W. J. Rappel, Dynamics of the globally coupled complex Ginzburg-Landau equation, Phys. Rev. A 46, 7347-7350(R) (1992). [CrossRef]
  135. N. K. Efremidis and D. N. Christodoulides, Discrete Ginzburg-Landau solitons, Phys. Rev. E 67, 026606 (2003). [CrossRef]
  136. K. Maruno, A. Ankiewicz, and N. Akhmediev, Exact localized and periodic solutions of the discrete complex Ginzburg-Landau equation, Opt. Commun. 221, 199-209 (2003). [CrossRef]
  137. N. K. Efremidis, D. N. Christodoulides, and K. Hizanidis, Two-dimensional discrete Ginzburg-Landau solitons, Phys. Rev. A 76, 043839 (2007). [CrossRef]
  138. C. M. Bender, Making sense of non-Hermitian Hamiltonians, Rep. Prog. Phys. 70, 947-1018 (2007). [CrossRef]
  139. K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and Z. H, Musslimani, Beam dynamics in PT symmetric optical lattices, Phys. Rev. Lett. 100, 103904 (2008). [CrossRef]
  140. C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, Observation of parity–time symmetry in optics, Nature Phys. 6, 192-195 (2010).
  141. V. V. Konotop, J. Yang, and D. Zezyulin, Nonlinear waves in PT -symmetric systems, Rev. Mod. Phys. 88, 035002 (2016).
  142. S. V. Suchkov, A. A. Sukhorukov, J. Huang, S. V. Dmitriev, C. Lee, and Yu. S. Kivshar, Nonlinear switching and solitons in PT-symmetric photonic systems, Laser Phot. Rev. 10, 177-213 (2016). [CrossRef]
  143. R. Driben and B. A. Malomed, Stability of solitons in parity-time-symmetric couplers, Opt. Lett. 36, 4323-4325 (2011). [CrossRef]
  144. N. V. Alexeeva, I. V. Barashenkov, A. A. Sukhorukov, and Y. S. Kivshar, Optical solitons in PT -symmetric nonlinear couplers with gain and loss, Phys. Rev. A 85, 063837 (2012). [CrossRef]
  145. G. Burlak and B. A. Malomed, Stability boundary and collisions of two-dimensional solitons in PT-symmetric couplers with the cubic-quintic nonlinearity, Phys. Rev. E 88, 062904 (2013). [CrossRef]
  146. V. V. Konotop, D. E. Pelinovsky, and D. A. Zezyulin, Discrete solitons in PT-symmetric lattices, Europhys. Lett. 100, 56006 (2012).
  147. C. Huang, C. Li, and L. Dong, Stabilization of multipole-mode solitons in mixed linear-nonlinear lattices with a PT symmetry, Opt. Exp. 21, 3917-3925 (2013). [CrossRef]
  148. D. Leykam, V. V. Konotop, and A. S. Desyatnikov, Discrete vortex solitons and parity time symmetry, Opt. Lett. 38, 371-373 (2013). [CrossRef]
  149. D. E. Pelinovsky, D. A. Zezyulin, and V. V. Konotop, Nonlinear modes in a generalized PT-symmetric discrete nonlinear Schrödinger equation, J. Phys. A: Math. Gen. 47, 085204 (2014). [CrossRef]
  150. J. D’Ambroise, P. G. Kevrekidis, and B. A. Malomed, Staggered parity-time-symmetric ladders with cubic nonlinearity, Phys. Rev. E 91, 033207 (2015). [CrossRef]
  151. M. Wimmer, A. Regensburger, M. A. Miri, C. Bersch, D. N. Christodoulides, U. Peschel, Observation of optical solitons in PT-symmetric lattices, Nature Commun. 6, 7782 (2015).
Preprints 95761 g001
Figure 1. Comparison of a typical discrete soliton, obtained as the numerical solution of Eq. (30), shown by chains of blue dots, and its counterpart produced by the VA (shown by red open circles). In this figure, ω = − 1 , see Eq. (29). The figure is borrowed from Ref. [22].
Figure 1. Comparison of a typical discrete soliton, obtained as the numerical solution of Eq. (30), shown by chains of blue dots, and its counterpart produced by the VA (shown by red open circles). In this figure, ω = − 1 , see Eq. (29). The figure is borrowed from Ref. [22].
Preprints 95761 g001
Preprints 95761 g002
Figure 2. A stable discrete vortex soliton with topological charge S = 1 , produced by Eq. (38) with ω = − 3.2 . The left and right panels show, respectively, distributions of the absolute value and phase of the complex wave function u m , n in the plane with coordinates m , n . The figure is borrowed from Ref. [86].
Figure 2. A stable discrete vortex soliton with topological charge S = 1 , produced by Eq. (38) with ω = − 3.2 . The left and right panels show, respectively, distributions of the absolute value and phase of the complex wave function u m , n in the plane with coordinates m , n . The figure is borrowed from Ref. [86].
Preprints 95761 g002
Preprints 95761 g003
Figure 3. Quasi-discrete optical solitons with vorticity S = 1 , created in the photorefractive crystal with an induced deep photonic lattice. The left and right panels display, respectively, the onsite- and intersite-centered vortex solitons. The figure is borrowed from Ref. [80].
Figure 3. Quasi-discrete optical solitons with vorticity S = 1 , created in the photorefractive crystal with an induced deep photonic lattice. The left and right panels display, respectively, the onsite- and intersite-centered vortex solitons. The figure is borrowed from Ref. [80].
Preprints 95761 g003
Preprints 95761 g004
Figure 4. An example of a stable bound state of three 2D fundamental (zero-vorticity) solitons, which are mutually out of phase. The solution is produced by the numerical solution of Eq. (38) with ω = − 1 . The figure is borrowed from Ref. [70].
Figure 4. An example of a stable bound state of three 2D fundamental (zero-vorticity) solitons, which are mutually out of phase. The solution is produced by the numerical solution of Eq. (38) with ω = − 1 . The figure is borrowed from Ref. [70].
Preprints 95761 g004
Preprints 95761 g005
Figure 5. An example of a stable bound state of two in-phase discrete vortex solitons, with topological charge S = 1 . The result is produced by the numerical solution of Eq. (38) with ω = − 1 . The figure is borrowed from Ref. [70].
Figure 5. An example of a stable bound state of two in-phase discrete vortex solitons, with topological charge S = 1 . The result is produced by the numerical solution of Eq. (38) with ω = − 1 . The figure is borrowed from Ref. [70].
Preprints 95761 g005
Preprints 95761 g006
Figure 6. (a) An example of a stable 2D discrete soliton with chemical potential μ = 4.05 , which corresponds to the right vertical line in (b), found in the mini-gap of the system based on Eqs. (40) and (41) with Δ = 0.5 and Q = π / 3 . (b) The dependence P ( μ ) of the norm of the solitons populating the mini-gap, which is identical to the interval of values of μ represented in the panel. The solitons are stable in the narrow shaded interval. The figure is borrowed from Ref. [101].
Figure 6. (a) An example of a stable 2D discrete soliton with chemical potential μ = 4.05 , which corresponds to the right vertical line in (b), found in the mini-gap of the system based on Eqs. (40) and (41) with Δ = 0.5 and Q = π / 3 . (b) The dependence P ( μ ) of the norm of the solitons populating the mini-gap, which is identical to the interval of values of μ represented in the panel. The solitons are stable in the narrow shaded interval. The figure is borrowed from Ref. [101].
Preprints 95761 g006
Preprints 95761 g007
Figure 7. Left: A stable 2D two-component discrete soliton with spontaneously broken symmetry between the components, generated by system (44). The 2D soliton, with total norm E ≡ E u + E v = 1.435 , is displayed by means of its cross section. Symbols labelled U N , V N and U A , V A stand, respectively, for the components of the numerically constructed soliton and its analytical counterpart predicted by the VA based on ansatz (47). Right: Families of 2D onsite-centered discrete solitons, generated by system (44), are shown by means of curves r ( E ) , where r is the asymmetry parameter (46). The dashed-dotted curve shows the VA prediction, while the solid and dashed ones represent stable and unstable solitons produced by the numerical solution. The figure is borrowed from Ref. [17].
Figure 7. Left: A stable 2D two-component discrete soliton with spontaneously broken symmetry between the components, generated by system (44). The 2D soliton, with total norm E ≡ E u + E v = 1.435 , is displayed by means of its cross section. Symbols labelled U N , V N and U A , V A stand, respectively, for the components of the numerically constructed soliton and its analytical counterpart predicted by the VA based on ansatz (47). Right: Families of 2D onsite-centered discrete solitons, generated by system (44), are shown by means of curves r ( E ) , where r is the asymmetry parameter (46). The dashed-dotted curve shows the VA prediction, while the solid and dashed ones represent stable and unstable solitons produced by the numerical solution. The figure is borrowed from Ref. [17].
Preprints 95761 g007
Preprints 95761 g008
Figure 8. Stable discrete solitons in the 2D Salerno model with competing nonlinearities [ μ < 0 in Eq. (48)], obtained for ω = − 4.22 : (a) a regular soliton at μ = − 0.2 ; (b) a cuspon at μ = − 0.88 . The figure is borrowed from Ref. [111].
Figure 8. Stable discrete solitons in the 2D Salerno model with competing nonlinearities [ μ < 0 in Eq. (48)], obtained for ω = − 4.22 : (a) a regular soliton at μ = − 0.2 ; (b) a cuspon at μ = − 0.88 . The figure is borrowed from Ref. [111].
Preprints 95761 g008
Preprints 95761 g009
Figure 9. Examples of discrete vortex solitons with topological charge S = 1 , produced by the 2D SM, based on Eq. (48). Profiles of the real part of the stationary wave function Φ m , n for the vortices of the onsite-centered (stable vortex cross) and intersite-centered (unstable vortex square) types are displayed in the top and bottom panels, respectively. Both solutions are obtained for μ = − 0.4 and ω = 7.0 . The figure is borrowed from Ref. [111].
Figure 9. Examples of discrete vortex solitons with topological charge S = 1 , produced by the 2D SM, based on Eq. (48). Profiles of the real part of the stationary wave function Φ m , n for the vortices of the onsite-centered (stable vortex cross) and intersite-centered (unstable vortex square) types are displayed in the top and bottom panels, respectively. Both solutions are obtained for μ = − 0.4 and ω = 7.0 . The figure is borrowed from Ref. [111].
Preprints 95761 g009
Preprints 95761 g010
Figure 10. Left: Juxtaposed profiles of ϕ m , n ( + ) and ϕ m , n ( − ) (solid and dashed lines, respectively) of a stable 2D discrete soliton of the semi-vortex type, in the central cross section, produced by a numerical solution of Eq. (55) with λ = 0.53 and γ = 0 . The soliton’s norm (see Eq. (57)) is N = 3.5 . Right: The same for a stable discrete soliton of the mixed-mode type, with λ = 0.58 , γ = 2 , and N = 2 . Values of the discrete fileds at lattice sites are connected by lines, for better visualization. The figure is borrowed from Ref. [118].
Figure 10. Left: Juxtaposed profiles of ϕ m , n ( + ) and ϕ m , n ( − ) (solid and dashed lines, respectively) of a stable 2D discrete soliton of the semi-vortex type, in the central cross section, produced by a numerical solution of Eq. (55) with λ = 0.53 and γ = 0 . The soliton’s norm (see Eq. (57)) is N = 3.5 . Right: The same for a stable discrete soliton of the mixed-mode type, with λ = 0.58 , γ = 2 , and N = 2 . Values of the discrete fileds at lattice sites are connected by lines, for better visualization. The figure is borrowed from Ref. [118].
Preprints 95761 g010
Preprints 95761 g011
Figure 11. Left and right: The dependence of the critical value of the SOC strength, λ cr , above which the 2D discrete solitons of the SV and MM types, produced by the numerical solution of Eq. (55) with γ = 0 and 2, respectively, suffer the delocalization, on the total soliton’s norm (57). The figure is borrowed from Ref. [118].
Figure 11. Left and right: The dependence of the critical value of the SOC strength, λ cr , above which the 2D discrete solitons of the SV and MM types, produced by the numerical solution of Eq. (55) with γ = 0 and 2, respectively, suffer the delocalization, on the total soliton’s norm (57). The figure is borrowed from Ref. [118].
Preprints 95761 g011
Preprints 95761 g012
Figure 12. A stable discrete vortex soliton with topological charge S = 3 , produced by the numerical solution of Eq. (60) with C = 0.01 and ω = − 2 . Panels (a) and (d) display level contours corresponding to Re u l , m , n = ± 0.25 and Im u l , m , n = ± 0.25 , respectively.The figure is borrowed from Ref. [119].
Figure 12. A stable discrete vortex soliton with topological charge S = 3 , produced by the numerical solution of Eq. (60) with C = 0.01 and ω = − 2 . Panels (a) and (d) display level contours corresponding to Re u l , m , n = ± 0.25 and Im u l , m , n = ± 0.25 , respectively.The figure is borrowed from Ref. [119].
Preprints 95761 g012
Preprints 95761 g013
Figure 13. Stable 3D multipole solutions of Eq. (60) with C = 0.1 and ω = − 2 . The top row depicts stable tightly-bound dipoles (with intersite separation d = 1 ): (a) straight, (b) 2D-diagonal, and (c) 3D-diagonal ones. (d) and (e): Quadrupoles set in the n = 0 plane, with intrinsic separation d = 1 and d = 2 , respectively. (f) and (g): Octupoles with d = 1 and d = 2 . Panel (h) displays the stability boundary C c r ( 3 D , d ) as a function of the intrinsic separation d for diagonal, oblique, and straight dipoles, octupoles, and quadrupoles, from top to bottom. The horizontal dashed line designates the stability threshold for the fundamental discrete soliton. Note that, for the quadrupoles (the bottom boundary), C quad ( 3 D , d ) is a linear function of d at d ≤ 4 (see the dashed straight line with slope 0.325 , included for the guidance). In panels (a)-(g), level contour corresponding to Re ( u l , m , n ) = ± 0.5 are shown by blue and red (colors, respectively. The figure is borrowed from Ref. [120].
Figure 13. Stable 3D multipole solutions of Eq. (60) with C = 0.1 and ω = − 2 . The top row depicts stable tightly-bound dipoles (with intersite separation d = 1 ): (a) straight, (b) 2D-diagonal, and (c) 3D-diagonal ones. (d) and (e): Quadrupoles set in the n = 0 plane, with intrinsic separation d = 1 and d = 2 , respectively. (f) and (g): Octupoles with d = 1 and d = 2 . Panel (h) displays the stability boundary C c r ( 3 D , d ) as a function of the intrinsic separation d for diagonal, oblique, and straight dipoles, octupoles, and quadrupoles, from top to bottom. The horizontal dashed line designates the stability threshold for the fundamental discrete soliton. Note that, for the quadrupoles (the bottom boundary), C quad ( 3 D , d ) is a linear function of d at d ≤ 4 (see the dashed straight line with slope 0.325 , included for the guidance). In panels (a)-(g), level contour corresponding to Re ( u l , m , n ) = ± 0.5 are shown by blue and red (colors, respectively. The figure is borrowed from Ref. [120].
Preprints 95761 g013
Preprints 95761 g014
Figure 14. Vortex cubes produced by the numerical solution of Eq. (60) with Λ = 2 and C = 0.1 . Panel (a) shows a stable composite mode, built of two parallel identical quasi-planar vortices with topological numbers S 1 = S 2 = 1 and a phase shift of π . Panel (b) shows an unstable vortex-antivortex cube, formed by vortices with opposite topological charges, S 1 = − S 2 . Level contours corresponding to Re ( u l , m , n ) = ± 0.5 are shown by blue and red colors, and the contours corresponding to Im ( u l , m , n ) = ± 0.5 are shown by green and yellow colors, respectively. The figure is borrowed from Ref. [120].
Figure 14. Vortex cubes produced by the numerical solution of Eq. (60) with Λ = 2 and C = 0.1 . Panel (a) shows a stable composite mode, built of two parallel identical quasi-planar vortices with topological numbers S 1 = S 2 = 1 and a phase shift of π . Panel (b) shows an unstable vortex-antivortex cube, formed by vortices with opposite topological charges, S 1 = − S 2 . Level contours corresponding to Re ( u l , m , n ) = ± 0.5 are shown by blue and red colors, and the contours corresponding to Im ( u l , m , n ) = ± 0.5 are shown by green and yellow colors, respectively. The figure is borrowed from Ref. [120].
Preprints 95761 g014
Preprints 95761 g015
Figure 15. A stable bound state of two mutually orthogonal vortex solitons with topological charges S = − 1 . The bound state is obtained as a numerical solution of Eqs. (62) with β = 0.5 , C = 0.01 , and chemical potentials of both components ω = − 2 . Top and bottom panels represent, respectively, to the first and second components of the system. The panels show contours of Re ( u l , m , n ) = ± 0.5 and Im ( u l , m , n ) = ± 0.5 . The bluish and reddish colors designate the positive and negative values, respectively. The figure is borrowed from Ref. [119].
Figure 15. A stable bound state of two mutually orthogonal vortex solitons with topological charges S = − 1 . The bound state is obtained as a numerical solution of Eqs. (62) with β = 0.5 , C = 0.01 , and chemical potentials of both components ω = − 2 . Top and bottom panels represent, respectively, to the first and second components of the system. The panels show contours of Re ( u l , m , n ) = ± 0.5 and Im ( u l , m , n ) = ± 0.5 . The bluish and reddish colors designate the positive and negative values, respectively. The figure is borrowed from Ref. [119].
Preprints 95761 g015
Preprints 95761 g016
Figure 16. An example of a 3D discrete skyrmion, produced by the numerical solution of Eqs. (63) with C = 0.05 , β = 0.25 , and ω = 2 . The left and right panels show, respectively, contours of the complex field, corresponding to Re ( u ) = + 1 / − 1 (blue/red colors) and Im ( u ) = + 1 / − 1 (green/yellow colors), and of the real field, corresponding to v = ( 1 , 0 , − 1 ) (blue, green and red colors, designating the contours from the inside out). The figure is borrowed from Ref. [124].
Figure 16. An example of a 3D discrete skyrmion, produced by the numerical solution of Eqs. (63) with C = 0.05 , β = 0.25 , and ω = 2 . The left and right panels show, respectively, contours of the complex field, corresponding to Re ( u ) = + 1 / − 1 (blue/red colors) and Im ( u ) = + 1 / − 1 (green/yellow colors), and of the real field, corresponding to v = ( 1 , 0 , − 1 ) (blue, green and red colors, designating the contours from the inside out). The figure is borrowed from Ref. [124].
Preprints 95761 g016
Preprints 95761 g017
Figure 17. An example of a semi-discrete localized spatiotemporal mode, generated by Eq. (64), which performs stable transverse motion under the action of the kick, defined according to (65), with a = 1.5 . The cross section of the plot at any fixed z shows the distribution of power u n ( τ ) 2 for each n. The figure is borrowed from Ref. [71].
Figure 17. An example of a semi-discrete localized spatiotemporal mode, generated by Eq. (64), which performs stable transverse motion under the action of the kick, defined according to (65), with a = 1.5 . The cross section of the plot at any fixed z shows the distribution of power u n ( τ ) 2 for each n. The figure is borrowed from Ref. [71].
Preprints 95761 g017
Preprints 95761 g018
Figure 18. A semi-discrete vortex soliton in a hexagonal quasi-discrete array of waveguides made in bulk silica. (a) The phase plate used for imprinting the vortex structure into the input beam. (b,c) Numerically simulated and experimentally observed (transient) intensity distributions in the transverse plane, with phase shifts 2 π / 3 between adjacent peaks, which represents the simplest vortical phase structure with topological charge S = 1 . The figure is borrowed from Ref. [129].
Figure 18. A semi-discrete vortex soliton in a hexagonal quasi-discrete array of waveguides made in bulk silica. (a) The phase plate used for imprinting the vortex structure into the input beam. (b,c) Numerically simulated and experimentally observed (transient) intensity distributions in the transverse plane, with phase shifts 2 π / 3 between adjacent peaks, which represents the simplest vortical phase structure with topological charge S = 1 . The figure is borrowed from Ref. [129].
Preprints 95761 g018
Preprints 95761 g019
Figure 19. The realization of the semi-discrete system 67: the array of planar optical waveguides (blue slabs), separated by gray isolating layers, with the continuous transverse coordinate, x, and the discrete one, n. As shown by the arrow, light is coupled into the array along the z direction. The figure is borrowed from Ref. [74].
Figure 19. The realization of the semi-discrete system 67: the array of planar optical waveguides (blue slabs), separated by gray isolating layers, with the continuous transverse coordinate, x, and the discrete one, n. As shown by the arrow, light is coupled into the array along the z direction. The figure is borrowed from Ref. [74].
Preprints 95761 g019
Preprints 95761 g020
Figure 20. Left and right panels display amplitude and phase profiles of stable onsite- and intersite-centered semi-discrete vortex solitons with topological charge S = 1 , produced by the numerical solution of Eq. (66) with parameters g , C = 0.48 , 0.1 and 0.77 , 0.15 , respectively. The norm (68) of both solutions is N = 100 . The figure is borrowed from Ref. [73].
Figure 20. Left and right panels display amplitude and phase profiles of stable onsite- and intersite-centered semi-discrete vortex solitons with topological charge S = 1 , produced by the numerical solution of Eq. (66) with parameters g , C = 0.48 , 0.1 and 0.77 , 0.15 , respectively. The norm (68) of both solutions is N = 100 . The figure is borrowed from Ref. [73].
Preprints 95761 g020
Preprints 95761 g021
Figure 21. Amplitude and phase profiles of a stable intersite-centered semi-discrete vortex soliton with topological charge S = 2 , produced by the numerical solution of Eq. (67) with C = 0.018 . The total power (69) of this soliton is P = 250 . The figure is borrowed from Ref. [74].
Figure 21. Amplitude and phase profiles of a stable intersite-centered semi-discrete vortex soliton with topological charge S = 2 , produced by the numerical solution of Eq. (67) with C = 0.018 . The total power (69) of this soliton is P = 250 . The figure is borrowed from Ref. [74].
Preprints 95761 g021
Preprints 95761 g022
Figure 22. Stability areas in the parameter plane C , g , produced by the numerical solution of Eq. (66) for onsite-centered semi-discrete vortex solitons with S = 2 (all colored regions), 3 (orange + brown + green), 4 (brown + green + blue), and 5 (green + blue + dark gray). For the convenience of plotting, the normalizations for S = 2 , 3, 4, and 5 are fixed as N = 400 , 900, 2500, and 4500, respectively. The figure is borrowed from Ref. [73].
Figure 22. Stability areas in the parameter plane C , g , produced by the numerical solution of Eq. (66) for onsite-centered semi-discrete vortex solitons with S = 2 (all colored regions), 3 (orange + brown + green), 4 (brown + green + blue), and 5 (green + blue + dark gray). For the convenience of plotting, the normalizations for S = 2 , 3, 4, and 5 are fixed as N = 400 , 900, 2500, and 4500, respectively. The figure is borrowed from Ref. [73].
Preprints 95761 g022
Preprints 95761 g023
Figure 23. Cross-sections, along m = 0 , of real and imaginary parts of components u m . n and v m . n of a stable PT -symmetric fundamental soliton produced by the numerical solution of Eq. (75), with ( C , γ ) = ( 0.3 , 0.4 ) , taking relation (74) into regard. The total power (73) of the soliton is P = 2 . The figure is borrowed from Ref. [83].
Figure 23. Cross-sections, along m = 0 , of real and imaginary parts of components u m . n and v m . n of a stable PT -symmetric fundamental soliton produced by the numerical solution of Eq. (75), with ( C , γ ) = ( 0.3 , 0.4 ) , taking relation (74) into regard. The total power (73) of the soliton is P = 2 . The figure is borrowed from Ref. [83].
Preprints 95761 g023
Preprints 95761 g024
Figure 24. Red and gray colors designate, respectively, stability and instability areas for PT -symmetric fundamental discrete solitons, produced by Eqs. (75) and (74), in the plane of the gain-loss coefficient, γ , and total power P, which is defined as per Eq. (73). The soliton solutions do not exist in the white area. The figure is borrowed from Ref. [83].
Figure 24. Red and gray colors designate, respectively, stability and instability areas for PT -symmetric fundamental discrete solitons, produced by Eqs. (75) and (74), in the plane of the gain-loss coefficient, γ , and total power P, which is defined as per Eq. (73). The soliton solutions do not exist in the white area. The figure is borrowed from Ref. [83].
Preprints 95761 g024
Preprints 95761 g025
Figure 25. Real (a) and imaginary (b) parts, and the phase structure (c), of field u m , n of a stable discrete PT -symmetric vortex soliton, produced by Eqs. (72) for ( C , γ ) = ( 0.06 , 0.4 ) , with propagation constant k = 1 and total power P = 1.65 , defined as per Eq. (73). The figure is borrowed from Ref. [83].
Figure 25. Real (a) and imaginary (b) parts, and the phase structure (c), of field u m , n of a stable discrete PT -symmetric vortex soliton, produced by Eqs. (72) for ( C , γ ) = ( 0.06 , 0.4 ) , with propagation constant k = 1 and total power P = 1.65 , defined as per Eq. (73). The figure is borrowed from Ref. [83].
Preprints 95761 g025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated