High-order-mode soliton structures in two-dimensional lattices with defocusing nonlinearity

Localized modes in a two-dimensional lattice with a pluslike geometry

We investigate analytically and numerically the existence and dynamical stability of different localized modes in a two-dimensional photonic lattice comprising a square plaquette inscribed in the dodecagon lattices. The eigenvalue spectrum of the underlying linear lattice is characterized by a net formed of one flat band and four dispersive bands. By tailoring the intersite coupling coefficient ratio, opening of gaps between two pairs of neighboring dispersive bands can be induced, while the fully degenerate flat band characterized by compact eigenmodes stays nested between two inner dispersive bands. The nonlinearity destabilizes the compact modes and gives rise to unique families of localized modes in the newly opened gaps, as well as in the semi-infinite gaps. The governing mechanism of mode localization in that case is the light energy self-trapping effect. We have shown the stability of a few families of nonlinear modes in gaps. The suggested lattice model may serve for probing various artificial flat-band systems such as ultracold atoms in optical lattices, periodic electronic networks, and polariton condensates.

Solitons and vortices in two-dimensional discrete nonlinear Schrödinger systems with spatially modulated nonlinearity

We consider a two-dimensional (2D) generalization of a recently proposed model [Gligorić et al., Phys. Rev. E 88, 032905 (2013)], which gives rise to bright discrete solitons supported by the defocusing nonlinearity whose local strength grows from the center to the periphery. We explore the 2D model starting from the anticontinuum (AC) limit of vanishing coupling. In this limit, we can construct a wide variety of solutions including not only single-site excitations, but also dipole and quadrupole ones. Additionally, two separate families of solutions are explored: the usual "extended" unstaggered bright solitons, in which all sites are excited in the AC limit, with the same sign across the lattice (they represent the most robust states supported by the lattice, their 1D counterparts being those considered as 1D bright solitons in the above-mentioned work), and the vortex cross, which is specific to the 2D setting. For all the existing states, we explore their stability (also analytically, when possible). Typical scenarios of instability development are exhibited through direct simulations.

Localized modes in mini-gaps opened by periodically modulated intersite coupling in two-dimensional nonlinear lattices

Spatially periodic modulation of the intersite coupling in two-dimensional (2D) nonlinear lattices modifies the eigenvalue spectrum by opening mini-gaps in it. This work aims to build stable localized modes in the new bandgaps. Numerical analysis shows that single-peak and composite two- and four-peak discrete static solitons and breathers emerge as such modes in certain parameter areas inside the mini-gaps of the 2D superlattice induced by the periodic modulation of the intersite coupling along both directions. The single-peak solitons and four-peak discrete solitons are stable in a part of their existence domain, while unstable stationary states (in particular, two-soliton complexes) may readily transform into robust localized breathers.

Two-dimensional paradigm for symmetry breaking: the nonlinear Schrödinger equation with a four-well potential

We study the existence and stability of localized modes in the two-dimensional (2D) nonlinear Schrödinger/Gross-Pitaevskii (NLS/GP) equation with a symmetric four-well potential. Using the corresponding four-mode approximation, we trace the parametric evolution of the trapped stationary modes, starting from the linear limit, and thus derive a complete bifurcation diagram for families of the stationary modes. This provides the picture of spontaneous symmetry breaking in the fundamental 2D setting. In a broad parameter region, the predictions based on the four-mode decomposition are found to be in good agreement with full numerical solutions of the NLS/GP equation. Stability properties of the stationary states coincide with those suggested by the corresponding discrete model in the large-amplitude limit. The dynamics of unstable modes is explored by means of direct simulations. Finally, in addition to the full analysis for the case of the self-attractive nonlinearity, the bifurcation diagram for the case of self-repulsion is briefly considered too.

Composition relation between gap solitons and Bloch waves in nonlinear periodic systems

We show with numerical computation and analysis that Bloch waves, at either the center or edge of the Brillouin zone, of a one dimensional nonlinear periodic system can be regarded as infinite chains composed of fundamental gap solitons (FGSs). This composition relation between Bloch waves and FGSs leads us to predict that there are n families of FGSs in the nth band gap of the corresponding linear periodic system, which is confirmed numerically. Furthermore, this composition relation can be extended to construct a class of solutions similar to Bloch waves but with multiple periods.

Two-dimensional discrete solitons in rotating lattices

We introduce a two-dimensional discrete nonlinear Schrödinger (DNLS) equation with self-attractive cubic nonlinearity in a rotating reference frame. The model applies to a Bose-Einstein condensate stirred by a rotating strong optical lattice, or light propagation in a twisted bundle of nonlinear fibers. Two types of localized states are constructed: off-axis fundamental solitons (FSs), placed at distance R from the rotation pivot, and on-axis (R=0) vortex solitons (VSs), with vorticities S=1 and 2 . At a fixed value of rotation frequency Omega , a stability interval for the FSs is found in terms of the lattice coupling constant C , 0<C<C_{cr}(R) , with monotonically decreasing C_{cr}(R) . VSs with S=1 have a stability interval, C[over ]_{cr};{(S=1)}(Omega)<C<C_{cr};{(S=1)}(Omega) , which exists for Omega below a certain critical value, Omega_{cr};{(S=1)} . This implies that the VSs with S=1 are destabilized in the weak-coupling limit by the rotation. On the contrary, VSs with S=2 , that are known to be unstable in the standard DNLS equation, with Omega=0 , are stabilized by the rotation in region 0<C<C_{cr};{(S=2)} , with C_{cr};{(S=2)} growing as a function of Omega . Quadrupole and octupole on-axis solitons are considered too, their stability regions being weakly affected by Omega not equal 0 .

Observation of higher-order solitons in defocusing waveguide arrays

We observe experimentally higher-order solitons in waveguide arrays with defocusing saturable nonlinearity. Such solitons can comprise several in-phase bright spots and are stable above a critical power threshold. We elucidate the impact of the nonlinearity saturation on the domains of existence and stability of the observed complex soliton states.