Two-dimensional solitons in nonlinear lattices

Dynamical self-trapping of two-dimensional binary solitons in cross-combined linear and nonlinear optical lattices

Dynamical and self-trapping properties of two-dimensional (2D) binary mixtures of Bose-Einstein condensates in cross-combined lattices, consisting of a one-dimensional (1D) linear optical lattice (LOL) in the x direction for the first component and a 1D nonlinear optical lattice (NOL) in the y direction for the second component, are analytically and numerically investigated. The existence and stability of 2D binary matter wave solitons in these settings are demonstrated both by variational analysis and by direct numerical integration of the coupled Gross-Pitaevskii equations. We find that in the absence of the NOL, binary solitons, stabilized by the action of the 1D LOL and by the attractive intercomponent interaction, can freely move in the y direction. In the presence of the NOL, we find, quite remarkably, the existence of threshold curves in the parameter space separating regions where solitons can move from regions where the solitons become dynamically self-trapped. The mechanism underlying the dynamical self-trapping phenomenon (DSTP) is qualitatively understood in terms of a dynamical barrier induced by the NOL, similar to the Peirls-Nabarro barrier of solitons in discrete lattices. DSTP is numerically demonstrated for binary solitons that are put in motion both by phase imprinting and by the action of external potentials applied in the y direction. In the latter case, we show that the trapping action of the NOL allows one to maintain a 2D binary soliton at rest in a nonequilibrium position of a parabolic trap or to prevent it from falling under the action of gravity. Possible applications of the results are also briefly discussed.

Surface gap solitons in the Schrödinger equation with quintic nonlinearity and a lattice potential

We demonstrate the existence of surface gap solitons, a special type of asymmetric solitons, in the one-dimensional nonlinear Schrödinger equation with quintic nonlinearity and a periodic linear potential. The nonlinearity is suddenly switched in a step-like fashion in the middle of the transverse spatial region, while the periodic linear potential is chosen in the form of a simple sin 2 lattice. The asymmetric nonlinearities in this work can be realized by the Feshbach resonance in Bose-Einstein condensates or by the photorefractive effect in optics. The major peaks in the gap soliton families are asymmetric and they are located at the position of the jump in nonlinearity (at x = 0). In addition, the major peaks of the two-peak and multi-peak solitons at the position x = 0 are higher than those after that position, at x > 0. And such phenomena are more obvious when the value of chemical potential is large, or when the difference of nonlinearity values across the jump is big. Along the way, linear stability analysis of the surface gap solitons is performed and the stability domains are identified. It is found that in this model, the solitons in the first band gap are mostly stable (excepting narrow domains of instability at the edges of the gap), while those in the second band gap are mostly unstable (excepting extremely narrow domains of stability for fundamental solitons). These findings are also corroborated by direct numerical simulations.

Triangular bright solitons in nonlinear optics and Bose-Einstein condensates

We demonstrate what we believe to be novel triangular bright solitons that can be supported by the nonlinear Schrödinger equation with inhomogeneous Kerr-like nonlinearity and external harmonic potential, which can be realized in nonlinear optics and Bose-Einstein condensates. The profiles of these solitons are quite different from the common Gaussian or sech envelope beams, as their tops and bottoms are similar to the triangle and inverted triangle functions, respectively. The self-defocusing nonlinearity gives rise to the triangle-up solitons, while the self-focusing nonlinearity supports the triangle-down solitons. Here, we restrict our attention only to the lowest-order fundamental triangular solitons. All such solitons are stable, which is demonstrated by the linear stability analysis and also clarified by direct numerical simulations. In addition, the modulated propagation of both types of triangular solitons, with the modulated parameter being the strength of nonlinearity, is also presented. We find that such propagation is strongly affected by the form of the modulation of the nonlinearity. For example, the sudden change of the modulated parameter causes instabilities in the solitons, whereas the gradual variation generates stable solitons. Also, a periodic variation of the parameter causes the regular oscillation of solitons, with the same period. Interestingly, the triangle-up and triangle-down solitons can change into each other, when the parameter changes the sign.

Dark soliton families in quintic nonlinear lattices

We prove that the dark solitons can be stable in the purely quintic nonlinear lattices, including the fundamental, tripole and five-pole solitons. These dark soliton families are generated on the periodic nonlinear backgrounds. The propagation constant affects the forms of these solitons, while the number of poles does not lead to the variation of the backgrounds. The dark solitons are stable only when the propagation constant is moderately large.

Localized modes and dark solitons sustained by nonlinear defects

Dark solitons and localized defect modes against periodic backgrounds are considered in arrays of waveguides with defocusing Kerr nonlinearity, constituting a nonlinear lattice. Bright defect modes are supported by a local increase in nonlinearity, while dark defect modes are supported by a local decrease in nonlinearity. Dark solitons exist for both types of defects, although in the case of weak nonlinearity, they feature side bright humps, making the total energy propagating through the system larger than the energy transferred by the constant background. All considered defect modes are found stable. Dark solitons are characterized by relatively narrow windows of stability. Interactions of unstable dark solitons with bright and dark modes are described.

Suppression of the critical collapse for one-dimensional solitons by saturable quintic nonlinear lattices

The stabilization of one-dimensional solitons by a nonlinear lattice against the critical collapse in the focusing quintic medium is a challenging issue. We demonstrate that this purpose can be achieved by combining a nonlinear lattice and saturation of the quintic nonlinearity. The system supports three species of solitons, namely, fundamental (even-parity) ones and dipole (odd-parity) modes of on- and off-site-centered types. Very narrow fundamental solitons are found in an approximate analytical form, and systematic results for very broad unstable and moderately broad partly stable solitons, including their existence and stability areas, are produced by means of numerical methods. Stability regions of the solitons are identified by means of systematic simulations. The stability of all the soliton species obeys the Vakhitov-Kolokolov criterion.

Aberrated surface soliton formation in a nonlinear 1D and 2D photonic crystal

We discuss a novel type of surface soliton-aberrated surface soliton-appearance in a nonlinear one dimensional photonic crystal and a possibility of this surface soliton formation in two dimensional photonic crystal. An aberrated surface soliton possesses a nonlinear distribution of the wavefront. We show that, in one dimensional photonic crystal, the surface soliton is formed at the photonic crystal boundary with the ambient medium. Essentially, that it occupies several layers at the photonic crystal boundary and penetrates into the ambient medium at a distance also equal to several layers, so that one can infer about light energy localization at the lateral surface of the photonic crystal. In the one dimensional case, the surface soliton is formed from an earlier formed soliton that falls along the photonic crystal layers at an angle which differs slightly from the normal to the photonic crystal face. In the two dimensional case, the soliton can appear if an incident Gaussian beam falls on the photonic crystal face. The influence of laser radiation parameters, optical properties of photonic crystal layers and ambient medium on the one dimensional surface soliton formation is investigated. We also discuss the influence of two dimensional photonic crystal configuration on light energy localization near the photonic crystal surface. It is important that aberrated surface solitons can be created at relatively low laser pulse intensity and for close values of alternating layers dielectric permittivity which allows their experimental observation.

Stabilization of spatiotemporal solitons in Kerr media by dispersive coupling

We introduce a mechanism to stabilize spatiotemporal solitons in Kerr nonlinear media, based on the dispersion of linear coupling between the field components forming the soliton states. Specifically, we consider solitons in a two-core guiding structure with inter-core coupling dispersion (CD). We show that CD profoundly affects properties of the solitons, causing the complete stabilization of the otherwise highly unstable spatiotemporal solitons in Kerr media with focusing nonlinearity. We also find that the presence of CD stimulates the formation of bound states, which, however, are unstable.

Twisted toroidal vortex solitons in inhomogeneous media with repulsive nonlinearity

Toroidal modes in the form of so-called Hopfions, with two independent winding numbers, a hidden one (twist s), which characterizes a circular vortex thread embedded into a three-dimensional soliton, and the vorticity around the vertical axis (m), appear in many fields, including field theory, ferromagnetics, and semi- and superconductors. Such topological states are normally generated in multicomponent systems, or as trapped quasilinear modes in toroidal potentials. We uncover that stable solitons with this structure can be created, without any linear potential, in the single-component setting with the strength of repulsive nonlinearity growing fast enough from the center to the periphery, for both steep and smooth modulation profiles. Toroidal modes with s=1 and vorticity m=0, 1, 2 are produced. They are stable for m≤1, and do not exist for s>1. An approximate analytical solution is obtained for the twisted ring with s=1, m=0. Under the application of an external torque, it rotates like a solid ring. The setting can be implemented in a Bose-Einstein condensate (BEC) by means of the Feshbach resonance controlled by inhomogeneous magnetic fields.

Stabilization of multipole-mode solitons in mixed linear-nonlinear lattices with a PT symmetry

We report the evolution of higher-order nonlinear states in a focusing cubic medium, where both the linear refractive index and the nonlinearity are spatially modulated by a complex optical lattice exhibiting a parity-time (PT) symmetry. We reveal that introduction of out-of-phase nonlinearity modulation makes possible the stabilization of higher-order solitons with number of poles up to 7, which are highly unstable in linear PT lattices. Under appropriate conditions, multipole-mode solitons with out-of-phase components in the neighboring lattice sites are completely stable provided that their power or propagation constant exceeds a critical value. Thus, our findings suggest an effective way for the realization of stable multipole-mode solitons in periodic potentials with gain-loss components.

Bright solitons in defocusing media with spatial modulation of the quintic nonlinearity

It has been recently demonstrated that self-defocusing (SDF) media with cubic nonlinearity, whose local coefficient grows from the center to the periphery fast enough, support stable bright solitons without the use of any linear potential. Our objective is to test the genericity of this mechanism for other nonlinearities, by applying it to one- and two-dimensional (1D and 2D) quintic SDF media. The models may be implemented in optics (in particular, in colloidal suspensions of nanoparticles), and the 1D model may be applied to the description of the Tonks-Girardeau gas of ultracold bosons. In 1D, the nonlinearity-modulation function is taken as g0+sinh2(βx). This model admits a subfamily of exact solutions for fundamental solitons. Generic soliton solutions are constructed in a numerical form and also by means of the Thomas-Fermi and variational approximations (TFA and VA). In particular, a new ansatz for the VA is proposed, in the form of "raised sech," which provides for an essentially better accuracy than the usual Gaussian ansatz. The stability of all the fundamental (nodeless) 1D solitons is established through the computation of the corresponding eigenvalues for small perturbations and also verified by direct simulations. Higher-order 1D solitons with two nodes have a limited stability region, all the modes with more than two nodes being unstable. It is concluded that the recently proposed inverted Vakhitov-Kolokolov stability criterion for fundamental bright solitons in systems with SDF nonlinearities holds here too. Particular exact solutions for 2D solitons are produced as well.

Spatial solitons in optofluidic waveguide arrays with focusing ultrafast Kerr nonlinearity

We present an optofluidic nonlinear waveguide array that is fabricated by selectively filling several strands of a photonic crystal fiber with the liquid CCl(4), which exhibits a large focusing ultrafast Kerr nonlinearity. We demonstrate a power dependent formation of a spatial soliton in this novel optofluidic device. The large thermo-optical effect of liquids enables us to control the characteristics of the spatial soliton formation in these nonlinear structures. This opens the road toward flexible designs and the realization of a new class of optofluidic devices with complex nonlinear landscapes and novel effects.

Stable bright and vortex solitons in photonic crystal fibers with inhomogeneous defocusing nonlinearity

We predict that a photonic crystal fiber whose strands are filled with a defocusing nonlinear medium can support stable bright solitons and also vortex solitons if the strength of the defocusing nonlinearity grows toward the periphery of the fiber. The domains of soliton existence depend on the transverse growth rate of the filling nonlinearity and nonlinearity of the core. Remarkably, solitons exist even when the core material is linear.

Stable two-dimensional solitons supported by radially inhomogeneous self-focusing nonlinearity

We demonstrate that modulation of the local strength of the cubic self-focusing (SF) nonlinearity in the two-dimensional geometry, in the form of a circle with contrast Δg of the SF coefficient relative to the ambient medium with a weaker nonlinearity, stabilizes a family of fundamental solitons against the critical collapse. The result is obtained in an analytical form, using the variational approximation and Vakhitov-Kolokolov stability criterion, and corroborated by numerical computations. For the small contrast, the stability interval of the soliton's norm scales as ΔN~Δg (the replacement of the circle by an annulus leads to a reduction of the stability region by perturbations breaking the axial symmetry). To further illustrate this mechanism, we demonstrate, in an exact form, the stabilization of one-dimensional solitons against the critical collapse under the action of a locally enhanced quintic SF nonlinearity.

Transfer and scattering of wave packets by a nonlinear trap

In the framework of a one-dimensional model with a tightly localized self-attractive nonlinearity, we study the formation and transfer (dragging) of a trapped mode by "nonlinear tweezers," as well as the scattering of coherent linear wave packets on the stationary localized nonlinearity. The use of a nonlinear trap for dragging allows one to pick up and transfer the relevant structures without grabbing surrounding "radiation." A stability border for the dragged modes is identified by means of analytical estimates and systematic simulations. In the framework of the scattering problem, the shares of trapped, reflected, and transmitted wave fields are found. Quasi-Airy stationary modes with a divergent norm, which may be dragged by a nonlinear trap moving at a constant acceleration, are briefly considered too.

Bright solitons from defocusing nonlinearities

We report that defocusing cubic media with spatially inhomogeneous nonlinearity, whose strength increases rapidly enough toward the periphery, can support stable bright localized modes. Such nonlinearity landscapes give rise to a variety of stable solitons in all three dimensions, including one-dimensional fundamental and multihump states, two-dimensional vortex solitons with arbitrarily high topological charges, and fundamental solitons in three dimensions. Solitons maintain their coherence in the state of motion, oscillating in the nonlinear potential as robust quasiparticles and colliding elastically. In addition to numerically found soliton families, particular solutions are found in an exact analytical form, and accurate approximations are developed for the entire families, including moving solitons.

Algebraic bright and vortex solitons in defocusing media

We demonstrate that spatially inhomogeneous defocusing nonlinear landscapes with the nonlinearity coefficient growing toward the periphery as (1+|r|(α)) support one- and two-dimensional fundamental and higher-order bright solitons, as well as vortex solitons, with algebraically decaying tails. The energy flow of the solitons converges as long as nonlinearity growth rate exceeds the dimensionality, i.e., α>D. Fundamental solitons are always stable, while multipoles and vortices are stable if the nonlinearity growth rate is large enough.

Two-dimensional solitons in media with stripe-shaped nonlinearity modulation

We introduce a model of media with the cubic attractive nonlinearity concentrated along a single or double stripe in the two-dimensional (2D) plane. The model can be realized in terms of nonlinear optics (in the spatial and temporal domains alike) and BEC. It is known from recent works that search for stable 2D solitons in models with a spatially localized self-attractive nonlinearity is a challenging problem. We make use of the variational approximation (VA) and numerical methods to investigate conditions for the existence and stability of solitons in the present setting. The result crucially depends on the transverse shape of the stripe: while the rectangular profile supports stable 2D solitons, its smooth Gaussian-shaped counterpart makes all the solitons unstable. This difference is explained, in particular, by the VA. The double stripe with the rectangular profile admits stable solitons of three distinct types: symmetric and asymmetric ones with a single-peak, and double-peak symmetric solitons. The shape and stability of the single-peak solitons of either type are accurately predicted by the VA. Collisions between identical stable solitons are briefly considered too, by means of direct simulations. Depending on the collision velocity, we observe excitation of intrinsic oscillations of the solitons, or their decay, or the collapse (catastrophic self-focusing).

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.