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.

Magneto-optical control of light collapse in bulk Kerr media

The Cotton-Mouton (Voigt) and Faraday effects induce adjustable linear and circular birefringence in optical media with external magnetic fields. We consider these effects as a technique for magneto-optical control of the transmission of bimodal light beams through Kerr-nonlinear crystals. Numerical analysis suggests that a properly applied magnetic field may accelerate, delay, or arrest the collapse of (2+1)D beams. Experimentally, the magnetic collapse acceleration is demonstrated in a bulk yttrium iron garnet (YIG) crystal.

Collisions between counter-rotating solitary vortices in the three-dimensional Ginzburg-Landau equation

We report results of collisions between coaxial vortex solitons with topological charges +/-S in the complex cubic-quintic Ginzburg-Landau equation. With the increase of the collision momentum, merger of the vortices into one or two dipole or quadrupole clusters of fundamental solitons (for S=1 and 2, respectively) is followed by the appearance of pairs of counter-rotating "unfinished vortices," in combination with a soliton cluster or without it. Finally, the collisions become elastic. The clusters generated by the collisions are very robust, while the "unfinished vortices," eventually split into soliton pairs.

Surface solitons in three dimensions

We study localized modes on the surface of a three-dimensional dynamical lattice. The stability of these structures on the surface is investigated and compared to that in the bulk of the lattice. Typically, the surface makes the stability region larger, an extreme example of that being the three-site "horseshoe"-shaped structure, which is always unstable in the bulk, while at the surface it is stable near the anticontinuum limit. We also examine effects of the surface on lattice vortices. For the vortex placed parallel to the surface, the increased stability-region feature is also observed, while the vortex cannot exist in a state normal to the surface. More sophisticated localized dynamical structures, such as five-site horseshoes and pyramids, are also considered.

Anisotropic solitons in dipolar bose-einstein condensates

Starting with a Gaussian variational ansatz, we predict anisotropic bright solitons in quasi-2D Bose-Einstein condensates consisting of atoms with dipole moments polarized perpendicular to the confinement direction. Unlike isotropic solitons predicted for the moments aligned with the confinement axis [Phys. Rev. Lett. 95, 200404 (2005)10.1103/PhysRevLett.95.200404], no sign reversal of the dipole-dipole interaction is necessary to support the solitons. Direct 3D simulations confirm their stability.

Symmetry breaking in linearly coupled dynamical lattices

We examine one- and two-dimensional models of linearly coupled lattices of the discrete-nonlinear-Schrödinger type. Analyzing ground states of the system with equal powers (norms) in the two components, we find a symmetry-breaking phenomenon beyond a critical value of the total power. Asymmetric states, with unequal powers in their components, emerge through a subcritical pitchfork bifurcation, which, for very weakly coupled lattices, changes into a supercritical one. We identify the stability of various solution branches. Dynamical manifestations of the symmetry breaking are studied by simulating the evolution of the unstable branches. The results present the first example of spontaneous symmetry breaking in two-dimensional lattice solitons. This feature has no counterpart in the continuum limit because of the collapse instability in the latter case.

Mobility of discrete solitons in quadratically nonlinear media

We study the mobility of solitons in lattices with quadratic (chi(2), alias second-harmonic-generating) nonlinearity. Using the notion of the Peierls-Nabarro potential and systematic numerical simulations, we demonstrate that, in contrast with their cubic (chi(3)) counterparts, the discrete quadratic solitons are mobile not only in the one-dimensional (1D) setting, but also in two dimensions (2D), in any direction. We identify parametric regions where an initial kick applied to a soliton leads to three possible outcomes: staying put, persistent motion, or destruction. On the 2D lattice, the solitons survive the largest kick and attain the largest speed along the diagonal direction.

Discrete surface solitons in two dimensions

We investigate fundamental localized modes in two-dimensional lattices with an edge (surface). The interaction with the edge expands the stability area for fundamental solitons, and induces a difference between dipoles oriented perpendicular and parallel to the surface. On the contrary, lattice vortex solitons cannot exist too close to the border. We also show, analytically and numerically, that the edge supports a species of localized patterns, which exists too but is unstable in the uniform lattice, namely, a horseshoe-shaped soliton, whose "skeleton" consists of three lattice sites. Unstable horseshoes transform themselves into a pair of ordinary solitons.

Mismatch management for optical and matter-wave quadratic solitons

We propose a way to control solitons in chi(2) (quadratically nonlinear) systems by means of periodic modulation imposed on the phase-mismatch parameter ("mismatch management," MM). It may be realized in the cotransmission of fundamental-frequency (FF) and second-harmonic (SH) waves in a planar optical waveguide via a long-period modulation of the usual quasi-phase-matching pattern of ferroelectric domains. In an altogether different physical setting, the MM may also be implemented by dint of the Feshbach resonance in a harmonically modulated magnetic field in a hybrid atomic-molecular Bose-Einstein condensate (BEC), with the atomic and molecular mean fields (MFs) playing the roles of the FF and SH, respectively. Accordingly, the problem is analyzed in two different ways. First, in the optical model, we identify stability regions for spatial solitons in the MM system, in terms of the MM amplitude and period, using the MF equations for spatially inhomogeneous configurations. In particular, an instability enclave is found inside the stability area. The robustness of the solitons is also tested against variation of the shape of the input pulse, and a threshold for the formation of stable solitons is found in terms of the power. Interactions between stable solitons are virtually unaffected by the MM. The second method (parametric approximation), going beyond the MF description, is developed for spatially homogeneous states in the BEC model. It demonstrates that the MF description is valid for large modulation periods, while, at smaller periods, non-MF components acquire gain, which implies destruction of the MF under the action of the high-frequency MM.

Skyrmion-like states in two- and three-dimensional dynamical lattices

We construct, in discrete two-component systems with cubic nonlinearity, stable states emulating Skyrmions of the classical field theory. In the two-dimensional case, an analog of the baby Skyrmion is built on the square lattice as a discrete vortex soliton of a complex field [whose vorticity plays the role of the Skyrmion's winding number (WN)], coupled to a radial "bubble" in a real lattice field. The most compact quasi-Skyrmion on the cubic lattice is composed of a nearly planar complex-field discrete vortex and a three-dimensional real-field bubble; unlike its continuum counterpart which must have WN=2, this stable discrete state exists with WN=1. Analogs of Skyrmions in the one-dimensional lattice are also constructed. Stability regions for all these states are found in an analytical approximation and verified numerically. The dynamics of unstable discrete Skyrmions (which leads to the onset of lattice turbulence) and their partial stabilization by external potentials are explored too.

X , Y , and Z waves: extended structures in nonlinear lattices

We propose a new type of waveforms in two-dimensional (2D) and three-dimensional (3D) discrete media-multilegged extended nonlinear structures (ENSs), built as arrays of lattice solitons (tiles and stones, in the 2D and 3D cases, respectively). We study the stability of the tiles and stones analytically, and then extend them numerically to complete ENS forms for both 2D and 3D lattices, aiming to single out stable ENSs. The predicted patterns can be realized in Bose-Einstein condensates trapped in deep optical lattices, crystals built of microresonators, and 2D photonic crystals. In the latter case, the patterns provide for a technique for writing reconfigurable virtual partitions in multipurpose photonic devices.

Stable three-dimensional optical solitons supported by competing quadratic and self-focusing cubic nonlinearities

We show that the quadratic (chi(2)) interaction of fundamental and second harmonics in a bulk dispersive medium, combined with self-focusing cubic (chi(3)) nonlinearity, give rise to stable three-dimensional spatiotemporal solitons (STSs), despite the possibility of the supercritical collapse, induced by the chi(3) nonlinearity. At exact phase matching (beta = 0) , the STSs are stable for energies from zero up to a certain maximum value, while for beta not equal 0 the solitons are stable in energy intervals between finite limits.

Stable solitons of even and odd parities supported by competing nonlocal nonlinearities

We introduce a one-dimensional phenomenological model of a nonlocal medium featuring focusing cubic and defocusing quintic nonlocal optical nonlinearities. By means of numerical methods, we find families of solitons of two types, even-parity (fundamental) and dipole-mode (odd-parity) ones. Stability of the solitons is explored by means of computation of eigenvalues associated with modes of small perturbations, and tested in direct simulations. We find that the stability of the fundamental solitons strictly follows the Vakhitov-Kolokolov criterion, whereas the dipole solitons can be destabilized through a Hamiltonian-Hopf bifurcation. The solitons of both types may be stable in the nonlocal model with only quintic self-attractive nonlinearity, in contrast with the instability of all solitons in the local version of the quintic model.

Discrete solitons and vortices in the two-dimensional Salerno model with competing nonlinearities

An anisotropic lattice model in two spatial dimensions, with on-site and intersite cubic nonlinearities (the Salerno model), is introduced, with emphasis on the case in which the intersite nonlinearity is self-defocusing, competing with on-site self-focusing. The model applies, for example, to a dipolar Bose-Einstein condensate trapped in a deep two-dimensional (2D) optical lattice. Soliton families of two kinds are found in the model: ordinary ones and cuspons, with peakons at the border between them. Stability borders for the ordinary solitons are found, while all cuspons (and peakons) are stable. The Vakhitov-Kolokolov criterion does not apply to cuspons, but for the ordinary solitons it correctly identifies the stability limits. In direct simulations, unstable solitons evolve into localized pulsons. Varying the anisotropy parameter, we trace a transition between the solitons in 1D and 2D versions of the model. In the isotropic model, we also construct discrete vortices of two types, on-site and intersite centered (vortex crosses and squares, respectively), and identify their stability regions. In simulations, unstable vortices in the noncompeting model transform into regular solitons, while in the model with the competing nonlinearities they evolve into localized vortical pulsons, which maintain their topological character. Bound states of regular solitons and vortices are constructed too, and their stability is identified.

Stable vortex tori in the three-dimensional cubic-quintic Ginzburg-Landau equation

We demonstrate the existence of stable toroidal dissipative solitons with the inner phase field in the form of rotating spirals, corresponding to vorticity S=0, 1, and 2, in the complex Ginzburg-Landau equation with the cubic-quintic nonlinearity. The stable solitons easily self-trap from pulses with embedded vorticity. The stability is corroborated by accurate computation of growth rates for perturbation eigenmodes. The results provide the first example of stable vortex tori in a 3D dissipative medium, as well as the first example of higher-order tori (with S=2) in any nonlinear medium. It is found that all stable vortical solitons coexist in a large domain of the parameter space; in smaller regions, there coexist stable solitons with either S=0 and S=1, or S=1 and S=2.

Discrete vector solitons in one-dimensional lattices in photorefractive media

We construct families of two-component spatial solitons in a one-dimensional lattice with saturable on-site nonlinearity (focusing or defocusing) in a photorefractive crystal. We identify 14 species of vector solitons, depending on their type (bright/dark), phase (in-phase/staggered), and location on the lattice (on/off-site). Two species of the bright/bright type form entirely stable soliton families, four species are partially stable (depending on the value of the propagation constant), while the remaining eight species are completely unstable. "Symbiotic" soliton pairs (of the bright/dark type), which contain components that cannot exist in isolation in the same model, are found as well.

Moving embedded lattice solitons

It was recently proved that solitons embedded in the spectrum of linear waves may exist in discrete systems, and explicit solutions for isolated unstable embedded lattice solitons (ELS) of a differential-difference version of a higher-order nonlinear Schrodinger equation were found [Gonzalez-Perez-Sandi, Fujioka, and Malomed, Physica D 197, 86 (2004)]. The discovery of these ELS gives rise to relevant questions such as the following: (1) Are there continuous families of ELS? (2) Can ELS be stable? (3) Is it possible for ELS to move along the lattice? (4) How do ELS interact? The present work addresses these questions by showing that a novel equation (a discrete version of a complex modified Korteweg-de Vries equation that includes next-nearest-neighbor couplings) has a two-parameter continuous family of exact ELS. These solitons can move with arbitrary velocities across the lattice, and the numerical simulations demonstrate that these ELS are completely stable. Moreover, the numerical tests show that these ELS are robust enough to withstand collisions, and the result of a collision is only a shift in the positions of the solitons. The model may apply to the description of a Bose-Einstein condensate with dipole-dipole interactions between the atoms, trapped in a deep optical-lattice potential.

Solitons in the Salerno model with competing nonlinearities

We consider a lattice equation (Salerno model) combining onsite self-focusing and intersite self-defocusing cubic terms, which may describe a Bose-Einstein condensate of dipolar atoms trapped in a strong periodic potential. In the continuum approximation, the model gives rise to solitons in a finite band of frequencies, with sechlike solitons near one edge, and an exact peakon solution at the other. A similar family of solitons is found in the discrete system, including a peakon; beyond the peakon, the family continues in the form of cuspons. Stability of the lattice solitons is explored through computation of eigenvalues for small perturbations, and by direct simulations. A small part of the family is unstable (in that case, the discrete solitons transform into robust pulsonic excitations); both peakons and cuspons are stable. The Vakhitov-Kolokolov criterion precisely explains the stability of regular solitons and peakons, but does not apply to cuspons. In-phase and out-of-phase bound states of solitons are also constructed. They exchange their stability at a point where the bound solitons are peakons. Mobile solitons, composed of a moving core and background, exist up to a critical value of the strength of the self-defocusing intersite nonlinearity. Colliding solitons always merge into a single pulse.

Three-dimensional spatiotemporal optical solitons in nonlocal nonlinear media

We demonstrate the existence of stable three-dimensional spatiotemporal solitons (STSs) in media with a nonlocal cubic nonlinearity. Fundamental (nonspinning) STSs forming one-parameter families are stable if their propagation constant exceeds a certain critical value that is inversely proportional to the range of nonlocality of nonlinear response. All spinning three-dimensional STSs are found to be unstable.

Vector solitons with an embedded domain wall

We present a class of soliton solutions to a system of two coupled nonlinear Schrödinger equations, with an intrinsic domain wall (DW) which separates regions occupied by two different fields. The model describes a binary mixture of two Bose-Einstein condensates (BECs) with interspecies repulsion. For the attractive or repulsive interactions inside each species, we find solutions which are bright or dark solitons in each component, while for the opposite signs of the intraspecies interaction, a bright-dark soliton pair is found (each time, with the intrinsic DW). These solutions can arise in the context of discrete lattices, and most of them can be supported in continuum settings by an external parabolic trap. The stability of the solitons with intrinsic DWs is examined, and the evolution of unstable ones is analyzed. We also briefly discuss the possibility of generating such families of solutions in the presence of linear coupling between the components, and an application of the model to bimodal light propagation in nonlinear optics.

Discrete solitons and vortices on anisotropic lattices

We consider the effects of anisotropy on solitons of various types in two-dimensional nonlinear lattices, using the discrete nonlinear Schrödinger equation as a paradigm model. For fundamental solitons, we develop a variational approximation that predicts that broad quasicontinuum solitons are unstable, while their strongly anisotropic counterparts are stable. By means of numerical methods, it is found that, in the general case, the fundamental solitons and simplest on-site-centered vortex solitons ("vortex crosses") feature enhanced or reduced stability areas, depending on the strength of the anisotropy. More surprising is the effect of anisotropy on the so-called "super-symmetric" intersite-centered vortices ("vortex squares"), with the topological charge equal to the square's size : we predict in an analytical form by means of the Lyapunov-Schmidt theory, and confirm by numerical results, that arbitrarily weak anisotropy results in dramatic changes in the stability and dynamics in comparison with the degenerate, in this case, isotropic, limit.

Accumulation of embedded solitons in systems with quadratic nonlinearity

Previous numerical studies have revealed the existence of embedded solitons (ESs) in a class of multiwave systems with quadratic nonlinearity, families of which seem to emerge from a critical point in the parameter space, where the zero solution has a fourfold zero eigenvalue. In this paper, the existence of such solutions is studied in a three-wave model. An appropriate rescaling casts the system in a normal form, which is universal for models supporting ESs through quadratic nonlinearities. The normal-form system contains a single irreducible parameter epsilon, and is tantamount to the basic model of type-I second-harmonic generation. An analytical approximation of Wentzel-Kramers-Brillouin type yields an asymptotic formula for the distribution of discrete values of epsilon at which the ESs exist. Comparison with numerical results shows that the asymptotic formula yields an exact value of the scaling index, -65, and a fairly good approximation for the numerical factor in front of the scaling term.

Fully three dimensional breather solitons can be created using feshbach resonances

We investigate the stability properties of breather solitons in a three-dimensional Bose-Einstein condensate with Feshbach resonance management of the scattering length and confined only by a one-dimensional optical lattice. We compare regions of stability in parameter space obtained from a fully 3D analysis with those from a quasi-two-dimensional treatment. For moderate confinement we discover a new island of stability in the 3D case, not present in the quasi-2D treatment. Stable solutions from this region have non-trivial dynamics in the lattice direction; hence, they describe fully 3D breather solitons. We demonstrate these solutions in direct numerical simulations and outline a possible way of creating robust 3D solitons in experiments in a Bose-Einstein condensate in a one-dimensional lattice. We point out other possible applications.

Stable spatiotemporal solitons in bessel optical lattices

We investigate the existence and stability of three-dimensional solitons supported by cylindrical Bessel lattices in self-focusing media. If the lattice strength exceeds a threshold value, we show numerically, and using the variational approximation, that the solitons are stable within one or two intervals of values of their norm. In the latter case, the Hamiltonian versus norm diagram has a swallowtail shape with three cuspidal points. The model applies to Bose-Einstein condensates and to optical media with saturable nonlinearity, suggesting new ways of making stable three-dimensional solitons and "light bullets" of an arbitrary size.

Three-dimensional nonlinear lattices: from oblique vortices and octupoles to discrete diamonds and vortex cubes

We construct a variety of novel localized topological structures in the 3D discrete nonlinear Schrödinger equation. The states can be created in Bose-Einstein condensates trapped in strong optical lattices and crystals built of microresonators. These new structures, most of which have no counterparts in lower dimensions, range from multipole patterns and diagonal vortices to vortex "cubes" (stack of two quasiplanar vortices) and "diamonds" (formed by two orthogonal vortices).