Source author record

Boris A. Malomed

Boris A. Malomed appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

nlin.PS physics.optics cond-mat.quant-gas quant-ph math-ph math.MP cond-mat.other nlin.SI hep-th math.AP physics.atom-ph cond-mat cond-mat.mes-hall cond-mat.soft nlin.CD nlin.CG physics.acc-ph

Catalog footprint

What is connected

224works

17topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Collisions and fusion of one- and two-dimensional solitons driven by potential troughs in the cubic-quintic nonlinear Schrödinger equations

We study the formation and collision of one- and two-dimensional (1D and 2D) Gaussian-shaped and flat-top (FT) solitons in the framework of the nonlinear Schrödinger equation with the cubic-quintic nonlinearity and two intersecting potential troughs. We find that Gaussian-Gaussian and Gaussian-FT collisions between the solitons, steered by the troughs, are quasi-elastic, while the collisions between FT solitons may be either quasi-elastic or inelastic, in the form of merger into a single FT soliton, thus spontaneously breaking the symmetry between the steering troughs. The Gaussian-FT collisions, being overall quasi-elastic, generate weak radiation fields.

preprint2026arXiv

Dynamical correlations and chimera-like states of nanoemitters coupled to plasmon-polaritons in a lattice of conducting nanorings

We systematically investigate semiclassical dynamics of the optical field produced by quantum nanoemitters (NEs) embedded in a periodic lattice of conducting nanorings (NRs), in which plasmon polaritons (PPs) are excited. The coupling between PPs and NEs through the radiated optical field leads to establishment of a significant cross-correlation between NEs, so that their internal dynamics (photocurrent affected by the laser irradiation) depends on the NR's plasma frequency $ω_{p}$. The transition to this regime,combined with the nonlinearity of the system, leads to a steep increase of the photocurrent in the NEs, as well as to non-smooth (chimera-like or chaotic) behavior in the critical (transition) region, where small variations of $ω_{p}$ lead to significant changes in the level of the NE pairwise cross-correlations. The chimera-like state is realized as coexistence of locally synchronized and desynchronized NE dynamical states. A fit of the dependence of the critical current on $ω_{p}$ is found, being in agreement with results of numerical simulations. The critical effect may help to design new optical devices, using dispersive nanolattices which are made available by modern nanoelectronics.

preprint2022arXiv

Experimental realisations of the fractional Schrödinger equation in the temporal domain

The fractional Schrödinger equation (FSE) -- a natural extension of the standard Schrödinger equation -- is the basis of fractional quantum mechanics. It can be obtained by replacing the kinetic-energy operator with a fractional derivative. Here, we report the experimental realisation of an optical FSE for femtosecond laser pulses in the temporal domain. Programmable holograms and the single-shot measurement technique are respectively used to emulate a \textit{Lévy waveguide} and to reconstruct the amplitude and phase of the pulses. Varying the Lévy index of the FSE and the initial pulse, the temporal dynamics is observed in diverse forms, including solitary, splitting and merging pulses, double Airy modes, and ``rain-like'' multi-pulse patterns. Furthermore, the transmission of input pulses carrying a fractional phase exhibits a ``fractional-phase protection'' effect through a regular (non-fractional) material. The experimentally generated fractional time-domain pulses offer the potential for designing optical signal-processing schemes.

preprint2022arXiv

Fundamental and vortex dissipative quadratic solitons supported by spatially localized gain

We consider settings providing the existence of stable two-dimensional (2D) dissipative solitons with zero and nonzero vorticity in optical media with the quadratic ($χ^{(2)}$) nonlinearity. To compensate the spatially uniform loss in both the fundamental-frequency (FF) and second-harmonic (SH) components of the system, a strongly localized amplifying region ("hot spot",HS), carrying the linear gain, is included, acting onto either the FF component or SH one. In both cases, the Gaussian radial gain profile supports stable fundamental dissipative solitons pinned to the HS. The structure of existence and stability domains for the 2D solitons is rather complex. They demonstrate noteworthy features, such as bistability and spontaneous symmetry breaking. A ring-shaped gain profile acting onto the FF component supports stable vortex solitons, with the winding number up to 5, and multipoles. Nontrivial transformation of vortex-soliton profiles upon either growth of the gain value or stretching the gain profile is observed.

preprint2022arXiv

General higher-order breathers and rogue waves in the two-component long-wave--short-wave resonance-interaction model

General higher-order breather and rogue wave (RW) solutions to the two-component long wave--short wave resonance interaction (2-LSRI) model are derived via the bilinear Kadomtsev-Petviashvili hierarchy reduction method and are given in terms of determinants. Under particular parametric conditions, the breather solutions can reduce to homoclinic orbits, or a mixture of breathers and homoclinic orbits. There are three families of RW solutions, which correspond to a simple root, two simple roots, and a double root of an algebraic equation related to the dimension reduction procedure. The first family of RW solutions consists of $\frac{N(N+1)}{2}$ bounded fundamental RWs, the second family is composed of $\frac{N_1(N_1+1)}{2}$ bounded fundamental RWs coexisting with another $\frac{N_2(N_2+1)}{2}$ fundamental RWs of different bounded state ($N,N_1,N_2$ being positive integers), while the third one have ${[\widehat{N}_1^2+\widehat{N}_2^2-\widehat{N}_1(\widehat{N}_2-1)]}$ fundamental bounded RWs ($\widehat{N}_1,\widehat{N}_2$ being non-negative integers). The second family can be regarded as the superpositions of the first family, while the third family can be the degenerate case of the first family under particular parameter choices. These diverse RW patterns are illustrated graphically.

preprint2022arXiv

Localized modes in nonlinear fractional systems with deep lattices

Solitons in the fractional space, supported by lattice potentials, have recently attracted much interest. We consider the limit of deep one- and two-dimensional (1D and 2D) lattices in this system, featuring finite bandgaps separated by nearly flat Bloch bands. Such spectra are also a subject of great interest in current studies. The existence, shapes, and stability of various localized modes, including fundamental gap and vortex solitons, are investigated by means of numerical methods; some results are also obtained with the help of analytical approximations. In particular, the 1D and 2D gap solitons, belonging to the first and second finite bandgaps, are tightly confined around a single cell of the deep lattice. Vortex gap solitons are constructed as four-peak \textquotedblleft squares" and \textquotedblleft rhombuses" with imprinted winding number $S=1$. Stability of the solitons is explored by means of the linearization and verified by direct simulations.

preprint2022arXiv

Multidimensional dissipative solitons and solitary vortices

This article offers a review of results for solitons in 2D and 3D models of nonlinear dissipative media. The existence of such solitons requires to maintain two balances: between nonlinear self-focusing and linear diffraction and/or dispersion, and between loss and gain. Due to these conditions, dissipative solitons (DSs) exist not in families, but as isolated solutions. The main issue is stability of 2D and 3D DSs, especially vortical ones. First, stable 2D DSs are presented in the framework of the complex Ginzburg-Landau equation with the cubic-quintic (CQ) nonlinearity, which combines linear and quintic loss with cubic gain. In addition to fundamental (zero-vorticity) DSs, stable spiral DSs are presented too, with vorticities 1 and 2. Stable 2D solitons were also found in a system of two linearly-coupled fields, with linear gain acting in one and linear loss in the other. In this case, the cubic loss (without quintic terms) is sufficient for the stability of fundamental and vortex DSs. In addition to truly localized states, weakly localized ones are presented too, in a model with nonlinear loss without explicit gain, the losses being compensated by influx of power from infinity. Other classes of 2D models which are considered here use spatially modulated loss or gain to predict many species of robust DSs, including ones featuring complex periodically recurring metamorphoses. Stable fundamental and vortex solitons are also produced by models with a trapping or spatially periodic potential. In the latter case, 2D gap DSs are considered. Further, 2D dissipative models with spin-orbit coupling give rise to stable semi-vortices, with vorticity carried by one component. Along with the 2D solitons, the review includes 3D fundamental and vortex DSs, stabilized by the CQ nonlinearity and/or external potentials, and collisions between them.

preprint2022arXiv

Observation of SQUID-like behavior in fiber laser with intra-cavity epsilon-near-zero effect

Establishing relations between fundamental effects in far-flung areas of physics is a subject of great interest in the current research. We here report realization of a novel photonic system akin to the radio-frequency superconducting quantum interference device (RF-SQUID), in a fiber laser cavity with epsilon-near-zero (ENZ) nanolayers as intra-cavity components. Emulating the RF-SQUID scheme, the photonic counterpart of the supercurrent, represented by the optical wave, circulates in the cavity, passing through effective optical potential barriers. Different ENZ wavelengths translate into distinct spectral outputs through the variation of cavity resonances, emulating the situation with a frequency-varying tank circuit in the RF-SQUID. Due to the presence of the ENZ element, the optical potential barrier is far lower for selected frequency components, granting them advantage in the gain-resource competition. The findings reported in this work provide a deeper insight into the ultrafast ENZ photonics, revealing a new path towards the design of nanophotonic on-chip devices with various operational functions, and offer a new approach to study superconducting and quantum-mechanical systems.

preprint2022arXiv

On-demand harnessing of photonic soliton molecules

Soliton molecules (SMs) are fundamentally important modes in nonlinear optical systems. It is a challenge to experimentally produce SMs with a required temporal separation in mode-locked fiber lasers. Here, we propose and realize an experimental scenario for harnessing SM dynamics in a laser setup. In particular, we tailor SMs in a mode-locked laser controlled by second-order group-velocity dispersion and dispersion losses: the real part of dispersion maintains the balance between the dispersion and nonlinearity, while the dispersion loss determines the balance of gain and losses. The experimental results demonstrate that the dispersion loss makes it possible to select desired values of the temporal separation (TS) in bound pairs of SMs in the system. Tunability of the SM's central wavelength and the corresponding hysteresis are addressed too. The demonstrated regime allows us to create multiple SMs with preselected values of the TS and central wavelength, which shows the potential of our setup for the design of optical data-processing schemes.

preprint2022arXiv

One- and two-dimensional solitons in spin-orbit-coupled Bose-Einstein condensates with fractional kinetic energy

We address effects of spin-orbit coupling (SOC), phenomenologically added to a two-component Bose-Einstein condensate composed of particles moving by Levy flights, in one- and two-dimensional (1D and 2D) settings. The corresponding system of coupled Gross-Pitaevskii equations includes fractional kinetic-energy operators, characterized by the Levy index, α< 2 (the normal kinetic energy corresponds to α= 2). The SOC terms, with strength λ, produce strong effects in the 2D case: they create families of stable solitons of the semi-vortex (SV) and mixed-mode (MM) types in the interval of 1 < α< 2, where the supercritical collapse does not admit the existence of stable solitons in the absence of the SOC. At λ--> 0, amplitudes of these solitons vanish as (λ)^{1/(α- 1)}.

preprint2022arXiv

Solitons in $\mathcal{PT}$-symmetric systems with spin-orbit coupling and critical nonlinearity

We construct families of one-dimensional (1D) stable solitons in two-component $\mathcal{PT}$-symmetric systems with spin-orbit coupling (SOC) and quintic nonlinearity, which plays the critical role in 1D setups. The system models light propagation in a dual-core waveguide with skewed coupling between the cores. Stability regions for the solitons are identified in the system's parameter space. They include the main semi-infinite gap, and an additional finite $\textit{annex gap}$. Stability boundaries are identified by means of simulations of the perturbed evolution, which agree with results produced by the linear-stability analysis for small perturbations. Distinct evolution scenarios are identified for unstable solitons. Generally, they suffer blowup or decay, while weakly unstable solitons transform into breathers. Due to a regularizing effect of SOC, stationary solitons are also found beyond the exceptional point, at which the $\mathcal{PT}$ symmetry breaks down, but they are unstable. Interactions between adjacent solitons are explored too, featuring rebound or merger followed by blowup. Slowly moving (tilted) solitons develop weak oscillations, while fast ones are completely unstable. Also considered is the reduced diffractionless system, which creates only unstable solitons.

preprint2022arXiv

Stability limits for modes held in alternating trapping-expulsive potentials

We elaborate a scheme of trapping-expulsion management (TEM), in the form of the quadratic potential periodically switching between confinement and expulsion, as a means of stabilization of two-dimensional dynamical states against the backdrop of the critical collapse driven by the cubic self-attraction with strength g. The TEM scheme may be implemented, as spatially or temporally periodic modulations, in optics or BEC, respectively. The consideration is carried out by dint of numerical simulations and variational approximation (VA). In terms of the VA, the dynamics amounts to a nonlinear Ermakov equation, which, in turn, is tantamount to a linear Mathieu equation. Stability boundaries are found as functions of g and parameters of the periodic modulation of the trapping potential. Below the usual collapse threshold, which is known, in the numerical form, as g < 5.85 (in the standard notation), the stability is limited by the onset of the parametric resonance. This stability limit, including the setup with the self-repulsive sign of the cubic term (g < 0), is accurately predicted by the VA. At g > 5.85, the collapse threshold is found with the help of full numerical simulations. The relative increase of the critical value of g above 5.85 is ~ 1.5%, which is a meaningful result, even if its size is small, because the collapse threshold is a universal constant, which is difficult to change.

preprint2022arXiv

Stable high-dimensional weak-light soliton molecules and their active control

Bound states of solitons, alias soliton molecules (SMs), are well known in one-dimensional (1D) systems, while making stable bound states of multidimensional solitons is a challenging problem because of the underlying instabilities. Here we propose a scheme for the creation of stable (2+1)D and (3+1)D optical SMs in a gas of cold Rydberg atoms, in which electromagnetically induced transparency (EIT) is induced by a control laser field. We show that, through the interplay of the EIT and the strong long-range interaction between the Rydberg atoms, the system gives rise to giant nonlocal Kerr nonlinearity, which in turn supports stable (2+1)D spatial optical SMs, as well as ring-shaped soliton necklaces, including rotating ones. They feature a large size, low generation power, and can be efficiently manipulated by tuning the nonlocality degree of the Kerr nonlinearity. Stable (3+1)D spatiotemporal optical SMs, composed of fundamental or vortex solitons, with low power and ultraslow propagation velocity, can also be generated in the system. These SMs can be stored and retrieved through the switching off and on of the control laser field. The findings reported here suggest applications to data processing and transmission in optical systems.

preprint2022arXiv

Trapping wave fields in an expulsive potential by means of linear coupling

We demonstrate the existence of confined states in one- and two-dimensional (1D and 2D) systems of two linearly-coupled components, with the confining harmonic-oscillator (HO) potential acting upon one component, and an expulsive anti-HO potential acting upon the other. The systems ca be implemented in optical and BEC dual-core waveguides. In the 1D linear system, codimension-one solutions are found in an exact form for the ground state (GS) and dipole mode (the first excited state). Generic solutions are produced by means of the variational approximation, and are found in a numerical form. Exact codimension-one solutions and generic numerical ones are also obtained for the GS and vortex states in the 2D system (the exact solutions are found for all values of the vorticity). Both the trapped and anti-trapped components of the bound states may be dominant ones, in terms of the norm. The localized modes may be categorized as bound states in continuum, as they coexist with delocalized ones. The 1D states, as well as the GS in 2D, are weakly affected and remain stable if the self-attractive or repulsive nonlinearity is added to the system. The self-attraction makes the vortex states unstable against splitting, while they remain stable under the action of the self-repulsion.

preprint2022arXiv

Two-dimensional Airy waves and three-wave solitons in quadratic media

We address the dynamics of two-dimensional (2D) truncated Airy waves and three-component solitons in the system of two fundamental-frequency and second-harmonic fields, coupled by quadratic (chi^(2)) terms. The system models second-harmonic-generating optical media and atomic-molecular mixtures in Bose-Einstein condensates. In addition to stable solitons, the system maintains truncated-Airy-waves states in either one of the fundamental-frequency components, represented by exact solutions, which are stable, unlike the Airy waves in the degenerate (two-component) chi^(2) system. It is also possible to imprint vorticity onto the 2D Airy modes. By means of systematic simulations, we examine interactions between truncated Airy waves originally carried by different fundamental-frequency components which are bending in opposite directions, through the second-harmonic field. The interaction leads to fusion of the input into a pair of narrow solitons. This is opposed to what happens in the 1D system, where the interacting Airy waves split into a large number of solitons. The interaction of truncated Airy waves carrying identical imprinted vorticities creates an additional pair of solitons, while opposite vorticities create a set of small-amplitude "debris" in the output. Slowly moving solitons colliding with a heavy truncated Airy wave bounce back, faster ones are absorbed by it, and collisions are quasi-elastic for fast solitons. Soliton-soliton collisions lead to merger into a single mode, or elastic passage, for lower and higher velocities, respectively.

preprint2022arXiv

Two-dimensional solitons in nonlocal media: a brief review

This is a review addressing soliton-like states in systems with nonlocal nonlinearity. The work on this topic has long history. Some findings, such as optical solitons supported by thermal nonlinearity, and by the orientational nonlinearity in liquid crystals, have been reviewed in the literature, therefore they are outlined in the present review in a brief form. Some other studies, such as those addressing models with fractional diffraction, which is represented by a linear nonlocal operator, have started recently, therefore it will be relevant to review them in detail when more results are accumulated; the present article provides a short outline of the latter topic. The main part of the article is a summary of results obtained for two-dimensional (2D) solitons in specific models originating in studies of Bose-Einstein condensates (BECs), which are sufficiently mature but have not yet been reviewed. These are, in particular, anisotropic quasi-2D solitons supported by long-range dipole-dipole interactions in a condensate of magnetic atoms, and giant vortex solitons, which are stable for high values of the winding number, as well as 2D vortex solitons of the latter type moving with self-acceleration. The vortex solitons are states of a hybrid type, which include matter-wave and electromagnetic-wave components. They are supported, in a binary BEC composed of two different atomic states, by the resonant interaction of the two-component matter waves with a microwave field which couples the two atomic states. Some other topics are included in the review in a brief form.

preprint2022arXiv

Vortex-ring quantum droplets in a radially-periodic potential

We establish stability and characteristics of two-dimensional (2D) vortex ring-shaped quantum droplets (QDs) formed by binary Bose-Einstein condensates (BECs). The system is modeled by the Gross-Pitaevskii (GP) equation with the cubic term multiplied by a logarithmic factor (as produced by the Lee-Huang-Yang correction to the mean-field theory) and a potential which is a periodic function of the radial coordinate. Narrow vortex rings with high values of the topological charge, trapped in particular circular troughs of the radial potential, are produced. These results suggest an experimentally relevant method for the creation of vortical QDs (thus far, only zero-vorticity ones have been reported). The 2D GP equation for the narrow rings is approximately reduced to the 1D form, which makes it possible to study the modulational stability of the rings against azimuthal perturbations. Full stability areas are delineated for these modes. The trapping capacity of the circular troughs is identified for the vortex rings with different winding numbers (WNs). Stable compound states in the form of mutually nested concentric multiple rings are constructed too, including ones with opposite signs of the WNs. Other robust compound states combine a modulationally stable narrow ring in one circular potential trough and an azimuthal soliton performing orbital motion in an adjacent one. The results may be used to design a device employing coexisting ring-shaped modes with different WNs for data storage.

preprint2021arXiv

Dragging spin-orbit-coupled solitons by a moving optical lattice

It is known that the interplay of the spin-orbit-coupling (SOC) and mean-field self-attraction creates stable two-dimensional (2D) solitons (ground states) in spinor Bose-Einstein condensates. However, SOC destroys the system's Galilean invariance, therefore moving solitons exist only in a narrow interval of velocities, outside of which the solitons suffer delocalization. We demonstrate that the application of a relatively weak moving optical lattice (OL), with the 2D or quasi-1D structure, makes it possible to greatly expand the velocity interval for stable motion of the solitons. The stability domain in the system's parameter space is identified by means of numerical methods. In particular, the quasi-1D OL produces a stronger stabilizing effect than its full 2D counterpart. Some features of the domain are explained analytically.

preprint2021arXiv

Embedded solitons in second-harmonic-generating lattices

Embedded solitons are exceptional modes in nonlinear-wave systems with the propagation constant falling in the system's propagation band. An especially challenging topic is seeking for such modes in nonlinear dynamical lattices (discrete systems). We address this problem for a system of coupled discrete equations modeling the light propagation in an array of tunnel-coupled waveguides with a combination of intrinsic quadratic (second-harmonic-generating) and cubic nonlinearities. Solutions for discrete embedded solitons (DESs) are constructed by means of two analytical approximations, adjusted, severally, for broad and narrow DESs, and in a systematic form with the help of numerical calculations. DESs of several types, including ones with identical and opposite signs of their fundamental-frequency and second-harmonic components, are produced. In the most relevant case of narrow DESs, the analytical approximation produces very accurate results, in comparison with the numerical findings. In this case, the DES branch extends from the propagation band into a semi-infinite gap as a family of regular discrete solitons. The study of stability of DESs demonstrates that, in addition to ones featuring the well-known property of semi-stability, there are linearly stable DESs which are genuinely robust modes.

preprint2021arXiv

External light control of three-dimensional ultrashort far-infrared pulses in an inhomogeneous array of carbon nanotubes

We present a study of the propagation of three-dimensional (3D) bipolar electromagnetic ultrashort pulses in an inhomogeneous array of semiconductor carbon nanotubes (CNTs) in the presence of a control high-frequency (HF) electric field. The inhomogeneity is present in the form of a layer with an increased concentration of conduction electrons, which acts as a barrier for the propagation of ultrashort electromagnetic pulses through the CNT array. The dynamics of the pulse is described by a nonlinear equation for the vector potential of the electromagnetic field (it takes the form of a 3D generalization of the sine-Gordon equation), derived from the Maxwell's equations and averaged over the period of the HF control field. By means of systematic simulations, we demonstrate that, depending on the amplitude and frequency of the HF control, the ultrashort pulse approaching the barrier layer either passes it or bounces back. The layer's transmissivity for the incident pulse is significantly affected by the amplitude and frequency of the HF control field, with the reflection coefficient nearly vanishing in intervals that make up a discrete set of transparency windows, which resembles the effect of the electromagnetically-induced transparency. Having passed the barrier, the ultrashort pulse continues to propagate, keeping its spatiotemporal integrity. The results may be used for the design of soliton valves, with the transmissivity of the soliton stream accurately controlled by the HF field.

preprint2021arXiv

Nonlinear dynamics of wave packets in tunnel-coupled harmonic-oscillator traps

We consider a two-component linearly-coupled system with the intrinsic cubic nonlinearity and the harmonic-oscillator (HO) confining potential. The system models binary settings in BEC and optics. In the symmetric system, with the HO trap acting in both components, we consider Josephson oscillations (JO) initiated by an input in the form of the HO's ground state (GS) or dipole mode (DM), placed in one component. With the increase of the strength of the self-focusing nonlinearity, spontaneous symmetry breaking (SSB) between the components takes place in the dynamical JO state. Under still stronger nonlinearity, the regular JO initiated by the GS input carry over into a chaotic dynamical state. For the DM input, the chaotization happens at smaller powers than for the GS, which is followed by SSB at a slightly stronger nonlinearity. In the system with the defocusing nonlinearity, SSB does not take place, and dynamical chaos occurs in a small area of the parameter space. In the asymmetric half-trapped system, with the HO potential applied to a single component, we first focus on the spectrum of confined binary modes in the linearized system. The spectrum is found analytically in the limits of weak and strong inter-component coupling, and numerically in the general case. Under the action of the coupling, the existence region of the confined modes shrinks for GSs and expands for DMs. In the full nonlinear system, the existence region for confined modes is identified in the numerical form. They are constructed too by means of the Thomas-Fermi approximation, in the case of the defocusing nonlinearity. Lastly, particular (non-generic) exact analytical solutions for confined modes, including vortices, in one- and two-dimensional asymmetric linearized systems are found. They represent bound states in the continuum.

preprint2021arXiv

Quadratic fractional solitons

We introduce a system combining the quadratic self-attractive or composite quadratic-cubic nonlinearity, acting in the combination with the fractional diffraction, which is characterized by its Lévy index $α$. The model applies to a gas of quantum particles moving by Lévy flights, with the quadratic term representing the Lee-Huang-Yang correction to the mean-field interactions. A family of fundamental solitons is constructed in a numerical form, while the dependence of its norm on the chemical potential characteristic is obtained in an exact analytical form. The family of \textit{quasi-Townes solitons}, appearing in the limit case of $α=1/2$, is investigated by means of a variational approximation. A nonlinear lattice, represented by spatially periodical modulation of the quadratic term, is briefly addressed too. The consideration of the interplay of competing quadratic (attractive) and cubic (repulsive) terms with a lattice potential reveals families of single-, double-, and triple-peak gap solitons (GSs) in two finite bandgaps. The competing nonlinearity gives rise to alternating regions of stability and instability of the GS, the stability intervals shrinking with the increase of the number of peaks in the GS.

preprint2021arXiv

Semidiscrete vortex solitons

We demonstrate a possibility of the creation of stable optical solitons combining one continuous and one discrete coordinate, with embedded vorticity, in an array of planar waveguides with intrinsic cubic-quintic nonlinearity. The same system may be realized in terms of the spatiotemporal light propagation in an array of tunnel-coupled optical fibers with the cubic-quintic nonlinearity. In contrast with zero-vorticity states, semidiscrete vortex solitons do not exist without the quintic term in the nonlinearity. Two types of the solitons, \emph{viz.}, intersite- and onsite-centered ones (IC and OC, respectively), with even and odd numbers $N$ of actually excited sites in the discrete direction, are identified. We consider the modes carrying the embedded vorticity $S=1$ and $2$. In accordance with their symmetry, the vortex solitons of the OC type exhibit an intrinsic core, while the IC solitons with a small $N$ may have a coreless structure. Facilitating their creation in the experiment, the modes reported in the present work may be much more compact states than their counterparts considered in other systems, and they feature strong anisotropy. They can be set in motion in the discrete direction, provided that the coupling constant exceeds a certain minimum value. Collisions between moving vortex solitons are considered too.

preprint2020arXiv

Collective Excitations of a One-Dimensional Quantum Droplet

We calculate the excitation spectrum of a one-dimensional self-bound quantum droplet in a two-component bosonic mixture described by the Gross-Pitaevskii equation (GPE) with cubic and quadratic nonlinearities. The cubic term originates from the mean-field energy of the mixture proportional to the effective coupling constant $δg$, whereas the quadratic nonlinearity corresponds to the attractive beyond-mean-field contribution. The droplet properties are governed by a control parameter $γ\propto δg N^{2/3}$, where $N$ is the particle number. For large $γ>0$ the droplet features the flat-top shape with the discrete part of its spectrum consisting of plane-wave Bogoliubov phonons propagating through the flat-density bulk and reflected by edges of the droplet. With decreasing $γ$ these modes cross into the continuum, sequentially crossing the particle-emission threshold at specific critical values. A notable exception is the breathing mode which we find to be always bound. The balance point $γ= 0$ provides implementation of a system governed by the GPE with an unusual quadratic nonlinearity. This case is characterized by the ratio of the breathing-mode frequency to the particle-emission threshold equal to 0.8904. As $γ$ tends to $-\infty$ this ratio tends to 1 and the droplet transforms into the soliton solution of the integrable cubic GPE.

preprint2020arXiv

Holding and transferring matter-wave solitons against gravity by spin-orbit-coupling tweezers

We consider possibilities to grasp and drag one-dimensional solitons in two-component Bose- Einstein condensates (BECs), under the action of gravity, by tweezers induced by spatially confined spin-orbit (SO) coupling applied to the BEC, with the help of focused laser illumination. Solitons of two types are considered, semi-dipoles and mixed modes. We find critical values of the gravity force, up to which the solitons may be held or transferred by the tweezers. The dependence of the critical force on the magnitude and spatial extension of the localized SO interaction, as well as on the solitons norm and speed (in the transfer regime), are systematically studied by means of numerical methods, and analytically with the help of a quasi-particle approximation for the soliton. In particular, a noteworthy finding is that the critical gravity force increases with the increase of the transfer speed (i.e., moving solitons are more robust than quiescent ones). Nonstationary regimes are addressed too, by considering abrupt application of gravity to solitons created in the weightless setting. In that case, solitons feature damped shuttle motion, provided that the gravity force does not exceed a dynamical critical value, which is smaller than its static counterpart. The results may help to design gravimeters based on ultracold atoms.

preprint2020arXiv

Josephson oscillations of chirality and identity in two-dimensional solitons in spin-orbit-coupled condensates

We investigate dynamics of two-dimensional chiral solitons of semi-vortex (SV) and mixed-mode (MM) types in spin-orbit-coupled Bose-Einstein condensates with the Manakov nonlinearity, loaded in a dual-core (double-layer) trap. The system supports two novel manifestations of Josephson phenomenology: one in the form of persistent oscillations between SVs or MMs with opposite chiralities in the two cores, and another one demonstrating robust periodic switching (identity oscillations) between SV in one core and MM in the other, provided that the strength of the inter-core coupling exceeds a threshold value. Below the threshold, the system creates composite states, which are asymmetric with respect to the two cores, or suffer the collapse. Robustness of the chirality and identity oscillations against deviations from the Manakov nonlinearity is investigated too. These dynamical regimes are possible only in the nonlinear system. In the linear one, exact stationary and dynamical solutions for SVs and MMs of the Bessel type are found. They sustain Josephson self-oscillations in different modes, with no interconversion between them.

preprint2020arXiv

Modulational instability, inter-component asymmetry and formation of quantum droplets in one-dimensional binary Bose gases

Quantum droplets are ultradilute liquid states which emerge from the competitive interplay of two Hamiltonian terms, the mean-field energy and beyond-mean-field correction, in a weakly interacting binary Bose gas. We relate the formation of droplets in symmetric and asymmetric two-component one-dimensional boson systems to the modulational instability of a spatially uniform state driven by the beyond-mean-field term. Asymmetry between the components may be caused by their unequal populations or unequal intra-component interaction strengths. Stability of both symmetric and asymmetric droplets is investigated. Robustness of the symmetric solutions against symmetry-breaking perturbations is confirmed.

preprint2020arXiv

Nonlinearity and discreteness: Solitons in lattices

An overview is given of basic models combining discreteness in their linear parts (i.e. the models are built as dynamical lattices) and nonlinearity acting at sites of the lattices or between the sites. The considered systems include the Toda and Frenkel-Kontorova lattices (including their dissipative versions), as well as equations of the discrete nonlinear Schroedinger (DNLS) and Ablowitz-Ladik (AL) types, and DNLS-AL combination in the form of the Salerno model. The interplay of discreteness and nonlinearity gives rise to a variety of states, most important ones being self-trapped discrete solitons. Basic results for one- and two-dimensional (1D and 2D) discrete solitons are collected in the review, including 2D solitons with embedded vorticity, and some results concerning mobility of discrete solitons. Main experimental findings are overviewed too. Models of the semi-discrete type, and basic results for solitons supported by them, are also considered, in a brief form. Perspectives for the development of topics covered the review are discussed throughout the text.

preprint2020arXiv

Reversible ultrafast soliton switching in dual-core highly nonlinear optical fibers

We experimentally investigate a nonlinear switching mechanism in a dual-core highly nonlinear optical fiber. We focus the input beam of femtosecond pulses on one core only, to identify transitions between inter-core oscillations, self-trapping in the cross core, and self-trapping of the pulse in the straight core. A model based in the system of coupled nonlinear Schrodinger equations provides surprisingly good agreement with the experimental findings.

preprint2020arXiv

Shortcuts to adiabaticity for an interacting Bose-Einstein condensate via exact solutions of the generalized Ermakov equation

Shortcuts to adiabatic expansion of the effectively one-dimensional Bose-Einstein condensate (BEC) loaded in the harmonic-oscillator (HO) trap is investigated by combining techniques of the variational approximation and inverse engineering. Piecewise-constant (discontinuous) intermediate trap frequencies, similar to the known bang-bang forms in the optimal-control theory, are derived from an exact solution of a generalized Ermakov equation. Control schemes considered in the paper include imaginary trap frequencies at short time scales, i.e., the HO potential replaced by the quadratic repulsive one. Taking into regard the BEC's intrinsic nonlinearity, results are reported for the minimal transfer time, excitation energy (which measures deviation from the effective adiabaticity), and stability for the shortcut-to-adiabaticity protocols. These results are not only useful for the realization of fast frictionless cooling, but also help to address fundamental problems of the quantum speed limit and thermodynamics.

preprint2020arXiv

Singular mean-field states: A brief review of recent results

This article provides a focused review of recent findings which demonstrate, in some cases quite counter-intuitively, the existence of bound states with a singularity of the density pattern at the center, while the states are physically meaningful because their total norm converges. One model of this type is based on the 2D Gross-Pitaevskii equation (GPE) which combines the attractive potential ~ 1/r^2 and the quartic self-repulsive nonlinearity, induced by the Lee-Huang-Yang effect (quantum fluctuations around the mean-field state). The GPE demonstrates suppression of the 2D quantum collapse, driven by the attractive potential, and emergence of a stable ground state (GS), whose density features an integrable singularity ~1/r^{4/3} at r --> 0. Modes with embedded angular momentum exist too, and they have their stability regions. A counter-intuitive peculiarity of the model is that the GS exists even if the sign of the potential is reversed from attraction to repulsion, provided that its strength is small enough. This peculiarity finds a relevant explanation. The other model outlined in the review includes 1D, 2D, and 3D GPEs, with the septimal (seventh-order), quintic, and cubic self-repulsive terms, respectively. These equations give rise to stable singular solitons, which represent the GS for each dimension D, with the density singularity ~1/r^{2/(4-D). Such states may be considered as a result of screening of a "bare" delta-functional attractive potential by the respective nonlinearity.

preprint2020arXiv

Singular solitons

We demonstrate that the commonly known concept, which treats solitons as nonsingular solutions produced by the interplay of nonlinear self-attraction and linear dispersion, may be extended to include modes with a relatively weak singularity at the central point, which keeps their integral norm convergent. Such states are generated by self-repulsion, which should be strong enough, namely, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schroedinger equations (NLSEs), respectively. Although such solutions seem counterintuitive, we demonstrate that they admit a straightforward interpretation as a result of screening of an additionally introduced attractive delta-functional potential by the defocusing nonlinearity. The strength ("bare charge") of the attractive potential is infinite in 1D, finite in 2D, and vanishingly small in 3D. Analytical asymptotics of the singular solitons at small and large distances are found, entire shapes of the solitons being produced in a numerical form. Complete stability of the singular modes is accurately predicted by the anti-Vakhitov-Kolokolov criterion (under the assumption that it applies to the model), as verified by means of numerical methods. In 2D, the NLSE with a quintic self-focusing term admits singular soliton solutions with intrinsic vorticity too, but they are fully unstable. We also mention that dissipative singular solitons can be produced by the model with a complex coefficient in front of the nonlinear term.

preprint2020arXiv

Splitting of two-component solitary waves from collisions with narrow potential barriers

We consider the interaction of two-component bright-bright solitons with a narrow potential barrier (splitter) in the framework of a system of two Gross-Pitaevskii (nonlinear-Schrödinger) equations modeling a binary Bose-Einstein condensate, with self-attraction in each component and cross-attraction between them. The objective is to study splitting of composite solitons, which may be used in the design of two-component soliton-interferometer schemes. We produce approximate analytic results, assuming a weak barrier and applying the perturbation theory in the limit in which the system is integrable and the solitary waves may be considered as exact solitons. We do this in the case of negligible interspecies interactions, and also when the nonlinearities are strongly asymmetric, allowing one to neglect the self-interaction in one of the species. Then, we use systematic simulations to study the transmissions of both components in regions outside these approximations and, in particular, to compare numerical results with their analytical countertparts. We concluded that there is an appreciable parameter range where one component is almost entirely transmitted through the barrier, while the other one is reflected. Excitation of internal vibrations in the passing and rebounding solitons is explored too, with a conclusion that it is weak in the regime of high-quality splitting.

preprint2020arXiv

Stabilization of one-dimensional Townes solitons by spin-orbit coupling in a dual-core system

It was recently demonstrated that two-dimensional Townes solitons (TSs) in two-component systems with cubic self-focusing, which are normally made unstable by the critical collapse, can be stabilized by linear spin-orbit coupling (SOC), in Bose-Einstein condensates and optics alike. We demonstrate that one-dimensional TSs, realized as optical spatial solitons in a planar dual-core waveguide with dominant quintic self-focusing, may be stabilized by SOC-like terms emulated by obliquity of the coupling between cores of the waveguide. Thus, SOC offers a universal mechanism for the stabilization of the TSs. A combination of systematic numerical considerations and analytical approximations identifies a vast stability area for skew-symmetric solitons in the system's main (semi-infinite) and annex (finite) bandgaps. Tilted ("moving") solitons are unstable, spontaneously evolving into robust breathers. For broad solitons, diffraction, represented by second derivatives in the system, may be neglected, leading to a simplified model with a finite bandgap. It is populated by skew-antisymmetric gap solitons, which are nearly stable close to the gap's bottom.

preprint2020arXiv

Stabilization of single- and multi-peak solitons in the fractional nonlinear Schroedinger equation with a trapping potential

We address the existence and stability of localized modes in the framework of the fractional nonlinear Schroedinger equation (FNSE) with the focusing cubic or focusing-defocusing cubic-quintic nonlinearity and a confining harmonic-oscillator (HO) potential. Approximate analytical solutions are obtained in the form of Hermite-Gauss modes. The linear stability analysis and direct simulations reveal that, under the action of the cubic self-focusing, the single-peak ground state and dipole mode are stabilized by the HO potential at values of the Levy index (the fractionality degree) alpha = 1 and alpha < 1, which lead, respectively, to the critical or supercritical collapse in free space. In addition to that, the inclusion of the quintic self-defocusing provides stabilization of higher-order modes, with the number of local peaks up to seven, at least.

preprint2020arXiv

Suppression of the quasi-two-dimensional quantum collapse in the attraction field by the Lee-Huang-Yang effect

Quantum collapse in three and two dimensions (3D and 2D) is induced by attractive potential ~ -1/r^2. It was demonstrated that the mean-field (MF) cubic self-repulsion in the 3D bosonic gas suppresses the collapse and creates the missing ground state (GS). However, the cubic nonlinearity is not strong enough to suppress the 2D collapse. We demonstrate that the Lee-Hung-Yang (LHY) quartic term, induced by quantum fluctuations around the MF state, is sufficient for the stabilization of the 2D gas against the collapse. By means of numerical solution of the Gross-Pitaevskii equation including the LHY term, as well as with the help of analytical methods, such as expansions of the wave function at small and large distances from the center and the Thomas-Fermi approximation, we construct stable GS, with a singular density, ~ 1/r^{4/3}, but convergent integral norm. Counter-intuitively, the stable GS exists even if the external potential is repulsive, with the strength falling below a certain critical value. An explanation to this finding is given. Along with the GS, singular vortex states are produced too, and their stability boundary is found analytically. Unstable vortices spontaneously transform into the stable GS, expelling the vorticity to periphery.

preprint2020arXiv

Symmetry breaking in a two-component system with repulsive interactions and linear coupling

We extend the well-known theoretical treatment of the spontaneous symmetry breaking (SSB) in two-component systems, combining linear coupling and self-attractive nonlinearity, to a system in which the linear coupling competes with repulsive interactions. First, we address one- and two-dimensional (1D and 2D) ground-state (GS) solutions and 2D vortex states with topological charges S=1 and 2, maintained by a confining harmonic-oscillator potential. The system can be implemented in BEC and optics. By means of the Thomas-Fermi approximation and numerical solution of the underlying coupled Gross-Pitaevskii equations, we demonstrate that SSB takes place, in the GSs and vortices alike, when the cross-component repulsion is stronger that the self-repulsion in each component. The SSB transition is categorized as a supercritical bifurcation, which gives rise to states featuring broken symmetry in an inner area, and intact symmetry in a surrounding layer. Unlike stable GSs and vortices with S=1, the states with S=2 are unstable against splitting. We also address SSB for 1D gap solitons in the system including a lattice potential. In this case, SSB takes place under the opposite condition, i.e., the cross-component repulsion must be weaker than the self-repulsion, and SSB is exhibited by antisymmetric solitons.

preprint2020arXiv

Symmetry breaking of spatial Kerr solitons in fractional dimension

We study symmetry breaking of solitons in the framework of a nonlinear fractional Schrödinger equation (NLFSE), characterized by its Lévy index, with cubic nonlinearity and a symmetric double-well potential. Asymmetric, symmetric, and antisymmetric soliton solutions are found, with stable asymmetric soliton solutions emerging from unstable symmetric and antisymmetric ones by way of symmetry-breaking bifurcations. Two different bifurcation scenarios are possible. First, symmetric soliton solutions undergo a symmetry-breaking bifurcation of the pitchfork type, which gives rise to a branch of asymmetric solitons, under the action of the self-focusing nonlinearity. Second, a family of asymmetric solutions branches off from antisymmetric states in the case of self-defocusing nonlinearity through a bifurcation of an inverted-pitchfork type. Systematic numerical analysis demonstrates that increase of the Lévy index leads to shrinkage or expansion of the symmetry-breaking region, depending on parameters of the double-well potential. Stability of the soliton solutions is explored following the variation of the Lévy index, and the results are confirmed by direct numerical simulations.

preprint2020arXiv

The Gross-Pitaevskii-Poisson model for an ultracold plasma: density waves and solitons

We introduce 1D and 2D models of a degenerate bosonic gas composed of ions with positive and negative charges (cations and anions). The system may exist in the mean-field condensate state, enabling the competition of the Coulomb coupling, contact repulsion, and kinetic energy of the particles, provided that their effective mass is reduced by means of a lattice potential. The model combines the Gross-Pitaevskii (GP) equations for the two-component wave function of the cations and anions, coupled to the Poisson equation for the electrostatic potential mediating the Coulomb coupling. The contact interaction in the GP system can be derived, in the Thomas-Fermi approximation, from a system of three GP equations, which includes the wave function of heavy neutral atoms. In the system with fully repulsive contact interactions, we construct stable spatially periodic patterns (density waves, DWs). The transition to DWs is identified by analysis of the modulational instability of a uniformly mixed neutral state. The DW pattern, which represents the system's ground state (GS), is predicted by a variational approximation. In 2D, a stable pattern is produced too, with a quasi-1D shape. The 1D system with contact self-attraction in each component produces bright solitons of three types: neutral ones, with fully mixed components; dipoles, with the components separated by the inter-species contact repulsion; and quadrupoles, with a layer of one component sandwiched between side lobes formed by the other. The transition from the neutral solitons to dipoles is accurately modeled analytically. A chart of the GSs of the different types (neutral solitons, dipoles, or quadrupoles) is produced. Different soliton species do not coexist as stable states. Collisions between traveling solitons are elastic for dipole-dipole pairs, while dipole-antidipole ones merge into stable quadrupoles via multiple collisions.

preprint2020arXiv

Tunable nonlinear spectra of anti-directional couplers

We produce transmission and reflection spectra of the anti-directional coupler (ADC) composed of linearly-coupled positive- and negative-refractive-index arms, with intrinsic Kerr nonlinearity. Both reflection and transmission feature two highly amplified peaks at two distinct wavelengths in a certain range of values of the gain, making it possible to design a wavelength-selective mode-amplification system. We also predict that a blend of gain and loss in suitable proportions can robustly enhance reflection spectra which are detrimentally affected by the attenuation, in addition to causing red and blue shifts owing to the Kerr effect. In particular, ADC with equal gain and loss coefficients, is considered in necessary detail.

preprint2020arXiv

Unstaggered-staggered solitons on one- and two-dimensional two-component discrete nonlinear Schrödinger lattices

We study coupled unstaggered-staggered soliton pairs emergent from a system of two coupled discrete nonlinear Schrödinger (DNLS) equations with the self-attractive on-site self-phase-modulation nonlinearity, coupled by the repulsive cross-phase-modulation interaction, on 1D and 2D lattice domains. These mixed modes are of a "symbiotic" type, as each component in isolation may only carry ordinary unstaggered solitons. While most work on DNLS systems addressed symmetric on-site-centered fundamental solitons, these models give rise to a variety of other excited states, which may also be stable. The simplest among them are antisymmetric states in the form of discrete twisted solitons, which have no counterparts in the continuum limit. In the extension to 2D lattice domains, a natural counterpart of the twisted states are vortical solitons. We first introduce a variational approximation (VA) for the solitons, and then correct it numerically to construct exact stationary solutions, which are then used as initial conditions for simulations to check if the stationary states persist under time evolution. Two-component solutions obtained include (i) 1D fundamental-twisted and twisted-twisted soliton pairs, (ii) 2D fundamental-fundamental soliton pairs, and (iii) 2D vortical-vortical soliton pairs. We also highlight a variety of other transient dynamical regimes, such as breathers and amplitude death. The findings apply to modeling binary Bose-Einstein condensates, loaded in a deep lattice potential, with identical or different atomic masses of the two components, and arrays of bimodal optical waveguides.

preprint2020arXiv

Vortex solitons in fractional nonlinear Schrödinger equation with the cubic-quintic nonlinearity

We address the existence and stability of vortex-soliton (VS) solutions of the fractional nonlinear Schrödinger equation (NLSE) with competing cubic-quintic nonlinearities and the Lévy index (fractionality) taking values 1 \leqα\leq2. Families of ring-shaped VSs with vorticities s = 1,2, and 3 are constructed in a numerical form. Unlike the usual two-dimensional NLSE (which corresponds to α = 2), in the fractional model VSs exist above a finite threshold value of the total power,P. Stability of the VS solutions is investigated for small perturbations governed by the linearized equation, and corroborated by direct simulations. Unstable VSs are broken up by azimuthal perturbations into several fragments, whose number is determined by the fastest growing eigenmode of small perturbations. The stability region, defined in terms of P, expands with the increase of α from 1 up to 2 for all s = 1, 2, and 3, except for steep shrinkage for s = 2 in the interval of 1\leqα\leq1.3.

preprint2019arXiv

Asymptotic dynamics of three-dimensional bipolar ultrashort electromagnetic pulses in an array of semiconductor carbon nanotubes

We study the propagation of three-dimensional bipolar ultrashort electromagnetic pulses in an array of semiconductor carbon nanotubes at times much longer than the pulse duration, yet still shorter than the relaxation time in the system. The interaction of the electromagnetic field with the electronic subsystem of the medium is described by means of Maxwell's equations, taking into account the field inhomogeneity along the nanotube axis beyond the approximation of slowly varying amplitudes and phases. A model is proposed for the analysis of the dynamics of an electromagnetic pulse in the form of an effective equation for the vector potential of the field. Our numerical analysis demonstrates the possibility of a satisfactory description of the evolution of the pulse field at large times by means of a three-dimensional generalization of the sine-Gordon and double sine-Gordon equations.

preprint2019arXiv

Coupled density-spin Bose-Einstein condensates dynamics and collapse in systems with quintic nonlinearity

We investigate the effects of spin-orbit coupling and Zeeman splitting on the coupled density-spin dynamics and collapse of the Bose-Einstein condensate driven by the quintic self-attraction in the same- and cross-spin channels. The characteristic feature of the collapse is the decrease in the width as given by the participation ratio of the density rather than by the expectation values of the coordinate. Qualitative arguments and numerical simulations reveal the existence of a critical spin-orbit coupling strength which either prohibits or leads to the collapse, and its dependence on other parameters, such as the condensates norm, spin-dependent nonlinear coupling, and the Zeeman splitting. The entire nonlinear dynamics critically depend on the initial spin sate.

preprint2019arXiv

Creating Oscillons and Oscillating Kinks in Two Scalar Field Theories

Oscillons are time-dependent, localized in space, extremely long-lived states in nonlinear scalar-field models, while kinks are topological solitons in one spatial dimension. In the present work, we show new classes of oscillons and oscillating kinks in a system of two nonlinearly coupled scalar fields in $1 + 1$ spatiotemporal dimensions. The solutions contain a control parameter, the variation of which produces oscillons and kinks with a flat-top shape. The model finds applications to condensed matter, cosmology, and high-energy physics.

preprint2019arXiv

Gap and embedded solitons in microwave-coupled binary condensates

It was recently found that, under the action of the spin-orbit coupling (SOC) and Zeeman splitting (ZS), binary BEC with intrinsic cubic nonlinearity supports families of gap solitons, provided that the kinetic energy is negligible in comparison with the SOC and ZS terms. We demonstrate that, also under the action of SOC and ZS, a similar setting may be introduced for BEC with two components representing different atomic states, resonantly coupled by microwave radiation, while the Poisson equation accounts for the feedback of the two-component atomic wave function onto the radiation. The microwave-mediated interaction induces an effective nonlinear trapping potential, which strongly affects the purport of the linear spectrum in this system. As a result, families of both gap and embedded solitons (those overlapping with the continuous spectrum) are found, being chiefly stable. The shape of the solitons features exact or broken skew symmetry. In addition to fundamental solitons (whose shape may or may not include a node), a family of dipole solitons is constructed too, which are even more stable than their fundamental counterparts. A nontrivial stability area is identified for moving solitons in the present system, which lacks Galilean invariance. Colliding solitons merge into a single one.

preprint2019arXiv

Mechanisms of Spatiotemporal Mode-Locking

Mode-locking is a process in which different modes of an optical resonator establish, through nonlinear interactions, stable synchronization. This self-organization underlies light sources that enable many modern scientific applications, such as ultrafast and high-field optics and frequency combs. Despite this, mode-locking has almost exclusively referred to self-organization of light in a single dimension - time. Here we present a theoretical approach, attractor dissection, for understanding three-dimensional (3D) spatiotemporal mode-locking (STML). The key idea is to find, for each distinct type of 3D pulse, a specific, minimal reduced model, and thus to identify the important intracavity effects responsible for its formation and stability. An intuition for the results follows from the 'minimum loss principle,' the idea that a laser strives to find the configuration of intracavity light that minimizes loss (maximizes gain extraction). Through this approach, we identify and explain several distinct forms of STML. These novel phases of coherent laser light have no analogues in 1D and are supported by experimental measurements of the three-dimensional field, revealing STML states comprising more than $10^7$ cavity modes. Our results should facilitate the discovery and understanding of new higher-dimensional forms of coherent light which, in turn, may enable new applications.

preprint2019arXiv

Quasi-one-dimensional approximation for Bose-Einstein condensates transversely trapped by a funnel potential

Starting from the standard three-dimensional (3D) Gross-Pitaevskii equation (GPE) and using a variational approximation, we derive an effective one-dimensional nonpolynomial Schrödinger equation (1D-NPSE) governing the axial dynamics of atomic Bose-Einstein condensates (BECs) under the action of a singular but physically relevant funnel-shaped transverse trap, i.e., an attractive 2D potential $\sim-1/r$ (where $r$ is the radial coordinate in the transverse plane), in combination with the repulsive self-interaction. Wave functions of the trapped BEC are regular, in spite of the potential's singularity. The model applies to a condensate of particles (small molecules) carrying a permanent electric dipole moment in the field of a uniformly charged axial thread, as well as to a quantum gas of magnetic atoms pulled by an axial electric current. By means of numerical simulations, we verify that the effective 1D-NPSE provides accurate static and dynamical results, in comparison to the full 3D GPE, for both repulsive and attractive signs of the intrinsic nonlinearity.

preprint2019arXiv

Soliton dynamics in a fractional complex Ginzburg-Landau model

The general objective of the work is to study dynamics of dissipative solitons in the framework of a one-dimensional complex Ginzburg-Landau equation (CGLE) of a fractional order. To estimate the shape of solitons in fractional models, we first develop the variational approximation for solitons of the fractional nonlinear Schrodinger equation (NLSE), and an analytical approximation for exponentially decaying tails of the solitons. Proceeding to numerical consideration of solitons in fractional CGLE, we study, in necessary detail, effects of the respective Levy index (LI) on the solitons' dynamics. In particular, dependence of stability domains in the model's parameter space on the LI is identified. Pairs of in-phase dissipative solitons merge into single pulses, with the respective merger distance also determined by LI.

preprint2019arXiv

Spatiotemporal dissipative solitons and vortices in a multi-transverse-mode fiber laser

We introduce a model for spatiotemporal modelocking in multimode fiber lasers, which is based on the (3+1)-dimensional cubic-quintic complex Ginzburg-Landau equation (cGLE) with conservative and dissipative nonlinearities and a 2-dimensional transverse trapping potential. Systematic numerical analysis reveals a variety of stable nonlinear modes, including stable fundamental solitons and breathers, as well as solitary vortices with winding number $n=1$, while vortices with $n=2$ are unstable, splitting into persistently rotating bound states of two unitary vortices. A characteristic feature of the system is bistability between the fundamental and vortex spatiotemporal solitons.

preprint2016arXiv

Analytical Solutions for Beams Passing Apertures with Sharp Boundaries

An approximation is elaborated for the paraxial propagation of diffracted beams, with both one- and two-dimensional cross sections, which are released from apertures with sharp boundaries. The approximation applies to any beam under the condition that the thickness of its edges is much smaller than any other length scale in the beam's initial profile. The approximation can be easily generalized for any beam whose initial profile has several sharp features. Therefore, this method can be used as a tool to investigate the diffraction of beams on complex obstacles. The analytical results are compared to numerical solutions and experimental findings, which demonstrates high accuracy of the approximation. For an initially uniform field confined by sharp boundaries, this solution becomes exact for any propagation distance and any sharpness of the edges. Thus, it can be used as an efficient tool to represent the beams, produced by series of slits with a complex structure, by a simple but exact analytical solution.

preprint2016arXiv

Commentary (invited update) to article "Spatiotemporal optical solitons", by B. A. Malomed, D. Mihalache, F. Wise, and L. Torner, J. Opt. B: Quantum Semiclass. Opt. 7, R53-R72 (2005)

This is an invited short update of the topic covered by the review article, which aims to briefly survey progress made in theoretical and experimental studies of multidimensional solitons since the publication of the review. The Commentary was invited as an addition to the original review article that will be reprinted in an e-book celebrating the 50th anniversary of Journal of Physics B (incorporating Journal of Optics B: Quantum and Semiclassical Optics).

preprint2016arXiv

Matter-wave soliton interferometer based on a nonlinear splitter

We elaborate a model of the interferometer which, unlike previously studied ones, uses a local (delta-functional) nonlinear repulsive potential, embedded into a harmonic-oscillator trapping potential, as the splitter for the incident soliton. An estimate demonstrates that this setting may be implemented by means of the localized Feshbach resonance controlled by a focused laser beam. The same system may be realized as a nonlinear waveguide in optics. Subsequent analysis produces an exact solution for scattering of a plane wave in the linear medium on the delta-functional nonlinear repulsive potential, and an approximate solution for splitting of the incident soliton when the ambient medium is nonlinear. The most essential result, obtained by means of systematic simulations, is that the use of the nonlinear splitter provides the sensitivity of the soliton-based interferometer to the target, inserted into one of its arms, which is much higher than the sensitivity provided by the usual linear splitter.

preprint2016arXiv

Motion of solitons in one-dimensional spin-orbit-coupled Bose-Einstein condensates

Solitons play a fundamental role in dynamics of nonlinear excitations. Here we explore the motion of solitons in one-dimensional uniform Bose-Einstein condensates subjected to a spin-orbit coupling (SOC). We demonstrate that the spin dynamics of solitons is governed by a nonlinear Bloch equation. The spin dynamics influences the orbital motion of the solitons leading to the spin-orbit effects in the dynamics of the macroscopic quantum objects (mean-field solitons). The latter perform oscillations with a frequency determined by the SOC, Raman coupling, and intrinsic nonlinearity. These findings reveal unique features of solitons affected by the SOC, which is confirmed by analytical considerations and numerical simulations of the underlying Gross-Pitaevskii equations.

preprint2016arXiv

Multidimensional solitons: Well-established results and novel findings

A brief review is given of some well-known and some very recent results obtained in studies of two- and three-dimensional (2D and 3D) solitons. Both zero-vorticity (fundamental) solitons and ones carrying vorticity S = 1 are considered. Physical realizations of multidimensional solitons in atomic Bose-Einstein condensates (BECs) and nonlinear optics are briefly discussed too. Unlike 1D solitons, which are typically stable, 2D and 3D ones are vulnerable to instabilities induced by the occurrence of the critical and supercritical collapse, respectively, in the same 2D and 3D models that give rise to the solitons. Vortex solitons are subject to a still stronger splitting instability. For this reason, a central problem is search for physical settings in which 2D and 3D solitons may be stabilized. The brief review addresses one well-established topic, viz., the stabilization of the 3D and 2D states, with S = 0 and 1, trapped in harmonic-oscillator (HO) potentials, and another topic which was developed very recently: the stabilization of 2D and 3D free-space solitons, which juxtapose components with S = 0 and (+/-)1 (semi-vortices and mixed modes), in a binary system with the spin-orbit coupling (SOC) between its components. The former model is based on the single cubic nonlinear Schroedinger/Gross-Pitaevskii equation (NLSE/GPE), while the latter one is represented by a system of two coupled GPEs. In both cases, generic situations are drastically different in the 2D and 3D geometries. In the 2D settings, the stabilization mechanism creates a stable ground state (GS, which was absent without the stabilization), whose norm falls below the threshold value at which the critical collapse sets in. In the 3D geometry, the supercritical collapse does not allow to create a GS, but metastable solitons can be constructed.

preprint2016arXiv

One- and two-dimensional solitons in PT-symmetric systems emulating the spin-orbit coupling

We introduce a two-dimensional (2D) system, which can be implemented in dual-core planar optical couplers with the Kerr nonlinearity in its cores, making it possible to blend effects of the PT symmetry, represented by the balanced linear gain and loss in the two cores, and spin-orbit coupling (SOC), emulated by a spatially biased coupling between the cores. Families of 1D and 2D solitons and their stability boundaries are identified. In the 1D setting, the addition of the SOC terms leads, at first, to shrinkage of the stability area for PT-symmetric solitons, which is followed by its rapid expansion. 2D solitons have their stability region too, in spite of the simultaneous action of two major destabilizing factors, viz., the collapse driven by the Kerr nonlinearity, and a trend towards spontaneous breakup of the gain-loss balance. In the limit of the SOC terms dominating over the intrinsic diffraction, the 1D system gives rise to a new model for gap solitons, which admits exact analytical solutions.

preprint2016arXiv

Resonant Meta-atoms with Nonlinearities on Demand

Nonlinear light-matter interactions and their applications are constrained by properties of available materials. The use of metamaterials opens the way to achieve precise control over electromagnetic properties at a microscopic level, providing new tools for experimental studies of complex nonlinear phenomena in photonics. Here a doubly resonant nonlinear meta-atom is proposed, analyzed and characterized in the GHz spectral range. The underlying structure is composed of a pair of split rings, resonant at both fundamental and nonlinear frequencies. The rings share a varactor diode, which serves as a microscopic source of nonlinearity. Flexible control over the coupling and near- and far-field patterns are reported, favoring the doubly resonant structure over other realizations. Relative efficiencies of the second and third harmonics, generated by the diode, are tailored by dint of the double-ring geometry, providing a guideline for selecting one frequency against another, using the design of the auxiliary structures. The on-demand control over the microscopic nonlinear properties enables developing a toolbox for experimental emulation of complex nonlinear phenomena.

preprint2016arXiv

Robust self-trapping of vortex beams in a saturable optical medium

We report the first observation of robust self-trapping of vortex beams propagating in a uniform condensed medium featuring local saturable self-focusing nonlinearity. Optical vortices with topological charge m=1, that remain self-trapped over ~ 5 Rayleigh lengths, are excited in carbon disulfide using a helical light beam at 532 nm and intensities from 8 to 10 GW/cm^2. At larger intensities, the vortex beams lose their stability, spontaneously breaking into two fragments. Numerical simulations based on the nonlinear Schrödinger equation including the three-photon absorption and nonpolynomial saturation of the refractive nonlinearity demonstrate close agreement with the experimental findings.

preprint2016arXiv

Spatiotemporal accessible solitons in fractional dimensions

We report solutions for solitons of the "accessible" type in globally nonlocal nonlinear media of fractional dimension (FD), viz., for self-trapped modes in the space of effective dimension $2<D\le3$ with harmonic-oscillator potential whose strength is proportional to the total power of the wave field. The solutions are categorized by a combination of radial, orbital and azimuthal quantum numbers $(n,l,m)$. They feature coaxial sets of vortical and necklace-shaped rings of different orders, and can be exactly written in terms of special functions that include Gegenbauer polynomials, associated Laguerre polynomials, and associated Legendre functions. The validity of these solutions is verified by direct simulation. The model can be realized in various physical settings emulated by FD spaces; in particular, it applies to excitons trapped in quantum wells.

preprint2016arXiv

Spontaneous symmetry breaking in a split potential box

We report results of the analysis of the spontaneous symmetry breaking (SSB) in the basic (actually, simplest) model which is capable to produce the SSB phenomenology in the one-dimensional setting. It is based on the Gross-Pitaevskii - nonlinear Schroedinger equation with the cubic self-attractive term and a double-well-potential built as an infinitely deep potential box split by a narrow (delta-functional) barrier. The barrier's strength, epsilon, is the single free parameter of the scaled form of the model. It may be implemented in atomic Bose-Einstein condensates and nonlinear optics. The SSB bifurcation of the symmetric ground state (GS) is predicted analytically in two limit cases, viz., for deep or weak splitting of the potential box by the barrier. For the generic case, a variational approximation (VA) is elaborated. The analytical findings are presented along with systematic numerical results. Stability of stationary states is studied through the calculation of eigenvalues for small perturbations, and by means of direct simulations. The GS always undergoes the SSB bifurcation of the supercritical type, as predicted by the VA at moderate values of epsilon, although the VA fails at small epsilon, due to inapplicability of the underlying ansatz in that case. However, the latter case is correctly treated by the approximation based on a soliton ansatz. On top of the GS, the first and second excited states are studied too. The antisymmetric mode (the first excited state) is destabilized at a critical value of its norm. The second excited state undergoes the SSB bifurcation, like the GS, but, unlike it, the bifurcation produces an unstable asymmetric mode. All unstable modes tend to spontaneously reshape into the asymmetric GS.

preprint2016arXiv

Spontaneous symmetry breaking of self-trapped and leaky modes in quasi-double-well potentials

We investigate competition between two phase transitions of the second kind induced by the self-attractive nonlinearity, viz., self-trapping of the leaky modes, and spontaneous symmetry breaking (SSB) of both fully trapped and leaky states. We use a one-dimensional mean-field model, which combines the cubic nonlinearity and a double-well-potential (DWP) structure with an elevated floor, which supports leaky modes (quasi-bound states) in the linear limit. The setting can be implemented in nonlinear optics and BEC. The order in which the SSB and self-trapping transitions take place with the growth of the nonlinearity strength depends on the height of the central barrier of the DWP: the SSB happens first if the barrier is relatively high, while self-trapping comes first if the barrier is lower. The SSB of the leaky modes is characterized by specific asymmetry of their radiation tails, which, in addition, feature a resonant dependence on the relation between the total size of the system and radiation wavelength. As a result of the SSB, the instability of symmetric modes initiates spontaneous Josephson oscillations. Collisions of freely moving solitons with the DWP structure admit trapping of an incident soliton into a state of persistent shuttle motion, due to emission of radiation. The study is carried out numerically, and basic results are explained by means of analytical considerations.

preprint2016arXiv

Stabilization of vortex beams in Kerr media by nonlinear absorption

We elaborate a new solution for the problem of stable propagation of transversely localized vortex beams in homogeneous optical media with self-focusing Kerr nonlinearity. Stationary nonlinear Bessel-vortex states are stabilized against azimuthal breakup and collapse by multiphoton absorption, while the respective power loss is offset by the radial influx of the power from an intrinsic reservoir. A linear stability analysis and direct numerical simulations reveal a region of stability of these vortices. Beams with multiple vorticities have their stability regions too. These beams can then form robust tubular filaments in transparent dielectrics as common as air, water and optical glasses at sufficiently high intensities. We also show that the tubular, rotating and speckle-like filamentation regimes, previously observed in experiments with axicon-generated Bessel beams, can be explained as manifestations of the stability or instability of a specific nonlinear Bessel-vortex state, which is fully identified.

preprint2016arXiv

Stable dipole solitons and soliton complexes in the nonlinear Schrodinger equation with periodically modulated nonlinearity

We develop a general classification of the infinite number of families of solitons and soliton complexes in the one-dimensional Gross-Pitaevskii/nonlinear Schrodinger equation with a nonlinear lattice pseudopotential, i.e., periodically modulated coefficient in front of the cubic term, which takes both positive and negative local values. This model finds direct implementations in atomic Bose-Einstein condensates and nonlinear optics. The most essential finding is the existence of two branches of dipole solitons (DSs), which feature an antisymmetric shape, essentially squeezed into a single cell of the nonlinear lattice. This soliton species was not previously considered in nonlinear lattices. We demonstrate that one branch of the DS family (namely, the one which obeys the Vakhitov-Kolokolov criterion) is stable, while unstable DSs spontaneously transform into stable fundamental solitons (FSs). The results are obtained in numerical and approximate analytical forms, the latter based on the variational approximation. Some stable bound states of FSs are found too.

preprint2016arXiv

Two-dimensional $χ^{2}$ solitons generated by the downconversion of Airy waves

Conversion of truncated Airy waves (AWs) carried by the second-harmonic (SH) component into axisymmetric $χ^{2}$ solitons is considered in the 2D system with the quadratic nonlinearity. The spontaneous conversion is driven by the parametric instability of the SH wave. The input in the form of the AW vortex is considered too. As a result, one, two, or three stable solitons emerge in a well-defined form, unlike the recently studied 1D setting, where the picture is obscured by radiation jets. Shares of the total power captured by the emerging solitons and conversion efficiency are found as functions of parameters of the AW input.

preprint2016arXiv

Two-dimensional solitons in dipolar Bose-Einstein condensates with spin-orbit coupling

We report families of two-dimensional (2D) composite solitons in spinor dipolar Bose-Einstein condensates, with two localized components linearly mixed by the spin-orbit coupling (SOC), and the intrinsic nonlinearity represented by the dipole-dipole interaction (DDI) between atomic magnetic moments polarized in-plane by an external magnetic field. Recently, stable solitons were predicted in the form of \textit{semi-vortices} (composites built of coupled fundamental and vortical components) in the 2D system combining the SOC and contact attractive interactions. Replacing the latter by the anisotropic long-range DDI, we demonstrate that, for a fixed norm of the soliton, the system supports a \emph{continuous family} of stable spatially asymmetric vortex solitons (AVSs), parameterized by an offset of the pivot of the vortical component relative to its fundamental counterpart. The offset is limited by a certain maximum value, while the energy of the AVS practically does not depend on the offset. At small values of the norm, the vortex solitons are subject to a weak oscillatory instability. In the present system, with the Galilean invariance broken by the SOC, the composite solitons are set in motion by a kick whose strength exceeds a certain depinning value. The kicked solitons feature a negative effective mass, drifting along a spiral trajectory opposite to the direction of the kick. A critical angular velocity, up to which the semi-vortices may follow rotation of the polarizing magnetic field, is found too.

preprint2016arXiv

Unidirectional transport of wave packets through tilted discrete breathers in nonlinear lattices with asymmetric defects

We consider the transfer of lattice wave packets through a tilted discrete breather (TDB) in opposite directions in the discrete nonlinear Schrödinger model with asymmetric defects, which may be realized as a Bose-Einstein condensate trapped in a deep optical lattice, or as optical beams in a waveguide array. A unidirectional transport mode is found, in which the incident wave packets, whose energy belongs to a certain interval between full reflection and full passage regions, pass the TDB only in one direction, while, in the absence of the TDB, the same lattice admits bi-directional propagation. The operation of this mode is accurately explained by an analytical consideration of the respective energy barriers. The results suggest that the TDB may emulate the unidirectional propagation of atomic and optical beams in various settings. In the case of the passage of the incident wave packet, the scattering TDB typically shifts by one lattice unit in the direction from which the wave packet arrives, which is an example of the tractor-beam effect, provided by the same system, in addition to the rectification of incident waves.

preprint2016arXiv

Vortex solitons in two-dimensional spin-orbit coupled Bose-Einstein condensates: effects of the Rashba-Dresselhaus coupling and the Zeeman splitting

We present an analysis of two-dimensional (2D) matter-wave solitons, governed by the pseudo-spinor system of Gross-Pitaevskii equations with self- and cross-attraction, which includes the spin-orbit coupling (SOC) in the general Rashba-Dresselhaus form, and, separately, the Rashba coupling and the Zeeman splitting. Families of semi-vortex (SV) and mixed-mode (MM) solitons are constructed, which exist and are stable in free space, as the SOC terms prevent the onset of the critical collapse and create the otherwise missing ground states in the form of the solitons. The Dresselhaus SOC produces a destructive effect on the vortex solitons, while the Zeeman term tends to convert the MM states into the SV ones, which eventually suffer delocalization. Existence domains and stability boundaries are identified for the soliton families. For physically relevant parameters of the SOC system, the number of atoms in the 2D solitons is limited by $\sim 1.5\times 10^{4}$. The results are obtained by means of combined analytical and numerical methods.

preprint2015arXiv

Composite localized modes in discretized spin-orbit-coupled Bose-Einstein condensates

We introduce a discrete model for binary spin-orbit-coupled (SOC) Bose-Einstein condensates (BEC) trapped in a deep one-dimensional optical lattice. Two different types of the couplings are considered, with spatial derivatives acting inside each species, or between the species. The discrete system with inter-site couplings dominated by the SOC, while the usual hopping is negligible, \textit{emulates} condensates composed of extremely heavy atoms, as well as those with opposite signs of the effective atomic masses in the two components.\ Stable localized composite states of miscible and immiscible types are constructed. The effect of the SOC on the immiscibility-miscibility transition in the localized complexes, which emulates the phase transition between insulating and conducting states in semiconductors, is studied.

preprint2015arXiv

Correlations and synchronization in a Bose-Fermi mixture

We study a Bose-Fermi mixture within the framework of the mean-field theory, including three possible regimes for the fermionic species: fully polarized, BCS, and unitarity. Starting from the 3D description and using the variational approximation (VA), we derive 1D and 2D systems of equations, under the corresponding confining potentials. This method produces a pair of nonlinear Schrödinger (NLS) equations coupled to algebraic equations for the transverse widths of the confined state. The equations incorporate interactions between atoms of the same species and between the species, assuming that the latter can be manipulated by means of the Feshbach resonance (FR). As an application, we explore spatial density correlations in the ground state (GS) between the species, concluding that they strongly depend on the sign and strength of the inter-species interaction. Also studied are the dynamics of the mixture in a vicinity of the GS and the corresponding spatiotemporal inter-species correlation. The correlations are strongly affected by the fermionic component, featuring the greatest variation in the unitary regime. Results produced by the VA are verified by comparison with full numerical solutions.

preprint2015arXiv

CP symmetry in optical systems

We introduce a model of a dual-core optical waveguide with opposite signs of the group-velocity-dispersion (GVD) in the two cores, and a phase-velocity mismatch between them. The coupler is embedded into an active host medium, which provides for the linear coupling of a gain-loss type between the two cores. The same system can be derived, without phenomenological assumptions, by considering the three-wave propagation in a medium with the quadratic nonlinearity, provided that the depletion of the second-harmonic pump is negligible. This linear system offers an optical realization of the charge-parity ($\mathcal{CP}$) symmetry, while the addition of the intra-core cubic nonlinearity breaks the symmetry. By means of direct simulations and analytical approximations, it is demonstrated that the linear system generates expanding Gaussian states, while the nonlinear one gives rise to broad oscillating solitons, as well as a general family of stable stationary gap solitons.

preprint2015arXiv

Dark solitons in dual-core waveguides with dispersive coupling

We report on new types of two-component one-dimensional dark solitons (DSs) in a model of a dual-core waveguide with normal group-velocity dispersion and Kerr nonlinearity in both cores, the coupling between which is dispersive too. In the presence of the dispersive coupling, quiescent DSs supported by the zero-frequency background are always gray, being stable with the out-of-phase background, i.e., for opposite signs of the fields in the cores. On the contrary, the background with a nonzero frequency supports quiescent black solitons which may be stable for both out- and in-phase backgrounds, if the dispersive coupling is sufficiently strong. Only DSs supported by the out-of-phase background admit an extension to the case of nonzero phase mismatch between the cores.

preprint2015arXiv

Discrete solitons in self-defocusing systems with $\mathcal{PT}$-symmetric defects

We construct families of discrete solitons (DSs) in an array of self-defocusing waveguides with an embedded $\mathcal{PT}$ (parity-time)-symmetric dimer, which is represented by a pair of waveguides carrying mutually balanced gain and loss. Four types of states attached to the embedded defect are found, namely, staggered and unstaggered bright localized modes and gray or anti-gray DSs. Their existence and stability regions expand with the increase of the strength of the coupling between the dimer-forming sites. The existence of the gray and staggered bright DSs is qualitatively explained by dint of the continuum limit. All the gray and anti-gray DSs are stable (some of them are unstable if the dimer carries the nonlinear $\mathcal{PT}$ symmetry, represented by balanced nonlinear gain and loss; in that case, the instability does not lead to a blowup, but rather creates oscillatory dynamical states). The boundary between the gray and anti-gray DSs is predicted in an approximate analytical form.

preprint2015arXiv

Engineering Bright Solitons to Enhance the Stability of Two-Component Bose-Einstein Condensates

We consider a system of coupled Gross-Pitaevskii (GP) equations describing a binary quasi-one-dimensional Bose-Einstein condensate (BEC) with intrinsic time-dependent attractive interactions, placed in a time-dependent expulsive parabolic potential, in a special case when the system is integrable (a deformed Manakov's system). Since the nonlinearity in the integrable system which represents binary attractive interactions exponentially decays with time, solitons are also subject to decay. Nevertheless, it is shown that the robustness of bright solitons can be enhanced in this system, making their respective lifetime longer, by matching the time dependence of the interaction strength (adjusted with the help of the Feshbach-resonance management) to the time modulation of the strength of the parabolic potential. The analytical results, and their stability, are corroborated by numerical simulations. In particular, we demonstrate that the addition of random noise does not impact the stability of the solitons.

preprint2015arXiv

Gap solitons attached to a gapless layer

We consider linear and nonlinear modes pinned to a grating-free (gapless) layer placed between two symmetric or asymmetric semi-infinite Bragg gratings (BGs), with a possible phase shift between them, in a medium with the uniform Kerr nonlinearity. The asymmetry is defined by a difference between bandgap widths in the two BGs. In the linear system, exact defect modes (DMs) are found. Composite gap solitons pinned to the central layer are found too, in analytical and numerical forms, in the nonlinear model. In the asymmetric system, existence boundaries for the DMs and gap solitons, due to the competition between attraction to the gapless layer and repulsion from the reflectivity step, are obtained analytically. Stability boundaries for solitons in the asymmetric system are identified by means of direct simulations. Collisions of moving BG solitons with the gapless layer are studied too.

preprint2015arXiv

Hybrid matter-wave-microwave solitons produced by the local-field effect

It was recently found that the electric local-field effect (LFE) can lead to strong coupling of atomic Bose-Einstein condensates (BECs) to off-resonant optical fields. We demonstrate that the magnetic LFE gives rise to a previously unexplored mechanism for coupling a (pseudo) spinor BEC or fermion gas to microwaves (MWs). We present a theory for the magnetic LFE, and find that it gives rise to a short-range attractive interaction between two components of the (pseudo) spinor, and a long-range interaction between them. The latter interaction, resulting from deformation of the magnetic field, is locally repulsive but globally attractive, in sharp contrast with its counterpart for the optical LFE, produced by phase modulation of the electric field. Our analytical results, confirmed by the numerical computations, show that the long-range interaction gives rise to modulational instability of the spatially uniform state, and creates stable ground states in the form of hybrid matter-wave-microwave solitons (which seem like one-dimensional magnetic monopoles), with a size much smaller than the MW wavelength, even in the presence of arbitrarily strong contact inter-component repulsion. The setting is somewhat similar to exciton-polaritonic condensates in semiconductor microcavities. The release of matter waves from the soliton may be used for the realization of an atom laser. The analysis also applies to molecular BECs with rotational states coupled by the electric MW field.

preprint2015arXiv

Interfaces between Bose-Einstein and Tonks-Girardeau atomic gases

We consider one-dimensional mixtures of an atomic Bose-Einstein condensate (BEC) and Tonks- Giradeau (TG) gas. The mixture is modeled by a coupled system of the Gross-Pitaevskii equation for the BEC and the quintic nonlinear Schroedinger equation for the TG component. An immiscibility condition for the binary system is derived in a general form. Under this condition, three types of BEC-TG interfaces are considered: domain walls (DWs) separating the two components; bubble-drops (BDs), in the form of a drop of one component immersed into the other (BDs may be considered as bound states of two DWs); and bound states of bright and dark solitons (BDSs). The same model applies to the copropagation of two optical waves in a colloidal medium. The results are obtained by means of systematic numerical analysis, in combination with analytical Thomas-Fermi approximations (TFAs). Using both methods, families of DW states are produced in a generic form. BD complexes exist solely in the form of a TG drop embedded into the BEC background. On the contrary, BDSs exist as bound states of TG bright and BEC dark components, and vice versa.

preprint2015arXiv

Nonlinear light propagation in cholesteric liquid crystals with a helical Bragg microstructure

Nonlinear optical propagation in cholesteric liquid crystals (CLC) with a spatially periodic helical molecular structure is studied experimentally and modeled numerically. This periodic structure can be seen as a Bragg grating with a propagation stopband for circularly polarized light. The CLC nonlinearity can be strengthened by adding absorption dye, thus reducing the nonlinear intensity threshold and the necessary propagation length. As the input power increases, a blue shift of the stopband is induced by the self-defocusing nonlinearity, leading to a substantial enhancement of the transmission and spreading of the beam. With further increase of the input power, the self-defocusing nonlinearity saturates, and the beam propagates as in the linear-diffraction regime. A system of nonlinear couple-mode equations is used to describe the propagation of the beam. Numerical results agree well with the experiment findings, suggesting that modulation of intensity and spatial profile of the beam can be achieved simultaneously under low input intensities in a compact CLC-based micro-device.

preprint2015arXiv

One- and two-dimensional solitons supported by singular modulation of quadratic nonlinearity

We introduce a model of one- and two-dimensional (1D and 2D) optical media with the $χ^{(2)}$ nonlinearity whose local strength is subject to cusp-shaped spatial modulation, $χ^{(2)}\sim r^{-α}$, with $α>0$, which can be induced by spatially nonuniform poling. Using analytical and numerical methods, we demonstrate that this setting supports 1D and 2D fundamental solitons, at $α<1$ and $α<2$, respectively. The 1D solitons have a small instability region, while the 2D solitons have a stability region at $α<0.5$ and are unstable at $α>0.5$. 2D solitary vortices are found too. They are unstable, splitting into a set of fragments, which eventually merge into a single fundamental soliton pinned to the cusp. Spontaneous symmetry breaking of solitons is studied in the 1D system with a symmetric pair of the cusp-modulation peaks.

preprint2015arXiv

Optically and Electrically Tunable Dirac Points and Zitterbewegung in Graphene-Based Photonic Superlattices

We demonstrate that graphene-based photonic superlattices provide a versatile platform for electrical and all-optical control of photonic beams with deep-subwavelength accuracy. Specifically, by inserting graphene sheets into periodic metallo-dielectric structures one can design optical superlattices that posses photonic Dirac points (DPs) at frequencies at which the spatial average of the permittivity of the superlattice, $\bar{ \varepsilon}$, vanishes. Similar to the well-known zero-$\bar{n}$ bandgaps, we show that these zero-$\bar{\varepsilon}$ DPs are highly robust against structural disorder. We also show that, by tuning the graphene permittivity via the optical Kerr effect or electrical doping, one can induce a spectral variation of the DP exceeding \SI{30}{\nano\meter}, at mid-IR and THz frequencies. The implications of this wide tunability for the photonic Zitterbewegung effect in a vicinity of the DP are explored too.

preprint2015arXiv

Rogue waves in a two-component Manakov system with variable coefficients and an external potential

We construct rogue waves (RWs) in a coupled two-mode system with the self-focusing nonlinearity of the Manakov type (equal SPM and XPM coefficients), spatially modulated coefficients, and a specially designed external potential. The system may be realized in nonlinear optics and Bose-Einstein condensates. By means of a similarity transformation, we establish a connection between solutions of the coupled Manakov system with spatially-variable coefficients and the basic Manakov model with constant coefficients. Exact solutions in the form of two-component Peregrine and dromion waves are obtained. The RW dynamics is analyzed for different choices of parameters in the underlying parameter space. Different classes of RW solutions are categorized by means of a naturally introduced control parameter which takes integer values.

preprint2015arXiv

Solitons and vortices in nonlinear potential wells

We consider self-trapping of topological modes governed by the one- and two-dimensional (1D and 2D) nonlinear-Schrodinger/Gross-Pitaevskii equation with effective single- and double-well (DW) nonlinear potentials induced by spatial modulation of the local strength of the self-defocusing nonlinearity. This setting, which may be implemented in optics and Bose-Einstein condensates, aims to extend previous studies, which dealt with single-well nonlinear potentials. In the 1D setting, we find several types of symmetric, asymmetric and antisymmetric states, focusing on scenarios of the spontaneous symmetry breaking. The single-well model is extended by including rocking motion of the well, which gives rise to Rabi oscillations between the fundamental and dipole modes. Analysis of the 2D single-well setting gives rise to stable modes in the form of ordinary dipoles, vortex-antivortex dipoles (VADs), and vortex triangles (VTs), which may be considered as produced by spontaneous breaking of the axial symmetry. The consideration of the DW configuration in 2D reveals diverse types of modes built of components trapped in the two wells, which may be fundamental states and vortices with topological charges m = 1 and 2, as well as VADs (with m = 0) and VTs (with m = 2).

preprint2015arXiv

Solitons in a forced nonlinear Schrödinger equation with the pseudo-Raman effect

Dynamics of solitons is considered in the framework of an extended nonlinear Schrödinger equation (NLSE), which is derived from a Zakharov-type model for wind-driven high-frequency (HF) surface waves in the ocean, coupled to damped low-frequency (LF) internal waves. The drive gives rise to a convective (but not absolute) instability in the system. The resulting NLSE includes a pseudo-stimulated-Raman-scattering (pseudo-SRS) term, which is a spatial-domain counterpart of the SRS term, a well-known ingredient of the temporal-domain NLSE in optics. Analysis of the field-momentum balance and direct simulations demonstrate that wavenumber downshift by the pseudo-SRS may be compensated by the upshift induced by the wind traction, thus maintaining robust bright solitons in both stationary and oscillatory forms; in particular, they are not destroyed by the underlying convective instability. Analytical soliton solutions are found in an approximate form and verified by numerical simulations. Solutions for soliton pairs are obtained in the numerical form.

preprint2015arXiv

Spontaneous symmetry breaking in nonlinear systems: An overview and a simple model

The paper combines two topics belonging to the general theme of the spontaneous symmetry breaking (SSB) in systems including two basic competing ingredients: the self-focusing cubic nonlinearity and a double-well-potential (DWP) structure. Such systems find diverse physical realizations, chiefly in optical waveguides, made of a nonlinear material and featuring a transverse DWP structure, and in models of atomic BEC with attractive inter-atomic interactions, loaded into a pair of symmetric potential wells coupled by tunneling across the separating barrier. With the increase of the nonlinearity strength, the SSB occurs at a critical value of the strength. The first part of the paper offers a brief overview of the topic. The second part presents a model which is designed as the simplest one capable to produce the SSB phenomenology in the one-dimensional geometry. The model is based on the DWP built as an infinitely deep potential box, which is split into two wells by a delta-functional barrier at the central point. Approximate analytical predictions for the SSB are produced for two limit cases: strong (deep) or weak (shallow) splitting of the potential box by the central barrier. Critical values of the strength of the nonlinearity at the SSB point, represented by the norm of the stationary wave field, are found in both cases (the critical strength is small in the former case, and large in the latter one). For the intermediate case, a less accurate variational approximation (VA) is developed.

preprint2015arXiv

Stability conditions for one-dimensional optical solitons in cubic-quintic-septimal media

Conditions for stable propagation of one-dimensional bright spatial solitons in media exhibiting optical nonlinearities up to the seventh-order are investigated. The results show well-defined stability regions even when all the nonlinear terms are focusing. Conditions for onset of the supercritical collapse of the optical beam are identified too. A variational approximation is used to predict dependence of the soliton propagation constant on the norm, and respective stability regions are identified using the Vakhitov-Kolokolov criterion. Analytical results obtained by means of the variational approximation are corroborated by numerical simulations of the cubic-quintic-septimal nonlinear Schroedinger equation.

preprint2015arXiv

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.

preprint2015arXiv

Stabilization of three-wave vortex beams in the waveguide

We consider two-dimensional (2D) localized vortical modes in the three-wave system with the quadratic ($χ^{(2)}$) nonlinearity, alias nondegenerate second-harmonic-generating system, guided by the isotropic harmonic-oscillator (HO) (alias parabolic) confining potential. In addition to the straightforward realization in optics, the system models mixed atomic-molecular Bose-Einstein condensates (BECs). The main issue is stability of the vortex modes, which is investigated through computation of instability growth rates for eigenmodes of small perturbations, and by means of direct simulations. The threshold of parametric instability for single-color beams, represented solely by the second harmonic (SH) with zero vorticity, is found in an analytical form with the help of the variational approximation (VA). Trapped states with vorticities $\left( +1,-1,0\right) $ in the two fundamental-frequency (FF) components and the SH one [the so-called \textit{hidden-vorticity} (HV) modes] are completely unstable. Also unstable are \textit{semi-vortices} (SVs), with component vorticities $% \left( 1,0,1\right) $. However, full vortices, with charges $\left( 1,1,2\right) $, have a well-defined stability region. Unstable full vortices feature regions of robust dynamical behavior, where they periodically split and recombine, keeping their vortical content.

preprint2015arXiv

Stable multiple vortices in collisionally inhomogeneous attractive Bose-Einstein condensates

We study stability of solitary vortices in the two-dimensional trapped Bose-Einstein condensate (BEC) with a spatially localized region of self-attraction. Solving the respective Bogoliubov-de Gennes equations and running direct simulations of the underlying Gross-Pitaevskii equation reveals that vortices with topological charge up to S = 6 (at least) are stable above a critical value of the chemical potential (i.e., below a critical number of atoms, which sharply increases with S). The largest nonlinearity-localization radius admitting the stabilization of the higher-order vortices is estimated analytically and accurately identified in a numerical form. To the best of our knowledge, this is the first example of a setting which gives rise to stable higher-order vortices, S > 1, in a trapped self-attractive BEC. The same setting may be realized in nonlinear optics too.

preprint2015arXiv

Stable NLS solitons in a cubic-quintic medium with a delta-function potential

We study the one-dimensional nonlinear Schrödinger equation with the cubic-quintic combination of attractive and repulsive nonlinearities, and a trapping potential represented by a delta-function. We determine all bound states with a positive soliton profile through explicit formulas and, using bifurcation theory, we describe their behavior with respect to the propagation constant. This information is used to prove their stability by means of the rigorous theory of orbital stability of Hamiltonian systems. The presence of the trapping potential gives rise to a regime where two stable bound states coexist, with different powers and same propagation constant.

preprint2015arXiv

Staggered $\mathcal{PT}$-symmetric ladders with cubic nonlinearity

We introduce a ladder-shaped chain with each rung carrying a $\mathcal{PT}$ -symmetric gain-loss dipole. The polarity of the dipoles is staggered along the chain, meaning that a rung bearing gain-loss is followed by one bearing loss-gain. This renders the system $\mathcal{PT}$-symmetric in both horizontal and vertical directions. The system is governed by a pair of linearly coupled discrete nonlinear Schrödinger (DNLS) equations with self-focusing or defocusing cubic onsite nonlinearity. Starting from the analytically tractable anti-continuum limit of uncoupled rungs and using the Newton's method for identifying solutions and parametric continuation in the inter-rung coupling for following the associated branches, we construct families of $\mathcal{PT}$-symmetric discrete solitons and identify their stability regions. Waveforms stemming from a single excited rung, as well as ones from multiple rungs are identified. Dynamics of unstable solitons is presented too.

preprint2015arXiv

Symmetry breaking in laser cavities

A brief introduction to the topic of spontaneous symmetry breaking (SSB) in conservative and dissipative nonlinear systems with an underlying double-well-potential structure is given. The reason is a discussion of a recent observation of the SSB a dual-core nanolaser cavity [5]. The effect is illustrated by means of a simple semi-analytically-tractable model (Fig. 1).

preprint2015arXiv

Transverse Anderson localization of light near Dirac points of photonic nanostructures

We perform a comparative study of the Anderson localization of light beams in disordered layered photonic nanostructures that, in the limit of periodic layer distribution, possess either a Dirac point or a Bragg gap in the spectrum of the wavevectors. In particular, we demonstrate that the localization length of the Anderson modes increases when the width of the Bragg gap decreases, such that in the vanishingly small bandgap limit, namely when a Dirac point is formed, even extremely high levels of disorders are unable to localize the optical modes located near the Dirac point. A comparative analysis of the key features of the propagation of Anderson modes formed in the Bragg gap or near the Dirac point is also presented. Our findings could provide valuable guidelines in assessing the influence of structural disorder on the functionality of a broad array of optical nanodevices.

preprint2015arXiv

Tunable storage of optical pulses in a tailored Bragg-grating structure

Scenarios for controllable creation, trapping and holding of single and multiple solitons in a specially designed nonlinear Bragg grating (BG) are proposed. The setting includes a chirped BG segment, which is linked via a local defect to a uniform BG with a built-in array of defects. A parabolic relation between the trapping position of the incident soliton and its power is obtained. Simultaneous trapping of two and three solitons at different locations is demonstrated too.

preprint2014arXiv

$\mathcal{PT}$-symmetric ladders with a scattering core

We consider a $\mathcal{PT}$-symmetric chain (ladder-shaped) system governed by the discrete nonlinear Schrödinger equation where the cubic nonlinearity is carried solely by two central "rungs" of the ladder. Two branches of scattering solutions for incident plane waves are found. We systematically construct these solutions, analyze their stability, and discuss non-reciprocity of the transmission associated with them. To relate these results to finite-size wavepacket dynamics, we also perform direct simulations of the evolution of the wavepackets, which confirm that the transmission is indeed asymmetric in this nonlinear system with the mutually balanced gain and loss.

preprint2014arXiv

Anderson Localization at the Subwavelength Scale and Loss Compensation for Surface-Plasmon Polaritons in Disordered Arrays of Metallic Nanowires

Using a random array of coupled metallic nanowires as a generic example of disordered plasmonic systems, we demonstrate that the structural disorder induces localization of light in these nanostructures at a deep-subwavelength scale. The ab initio analysis is based on solving the complete set of 3D Maxwell equations. We find that random variations of the radius of coupled plasmonic nanowires are sufficient to induce the Anderson localization (AL) of surface-plasmon polaritons (SPPs), the size of these trapped modes being significantly smaller than the wavelength. Remarkably, the optical-gain coefficient, needed to compensate losses in the plasmonic components of the system, is much smaller than the loss coefficient of the metal, which is obviously beneficial for the realization of the AL in plasmonic nanostructures. The dynamics of excitation and propagation of the Anderson-localized SPPs are addressed too.

preprint2014arXiv

Breather-like solitons extracted from the Peregrine rogue wave

Based on the Peregrine solution (PS) of the nonlinear Schrödinger (NLS) equation, the evolution of rational fraction pulses surrounded by zero background is investigated. These pulses display the behavior of a breather-like solitons. We study the generation and evolution of such solitons extracted, by means of the spectral-filtering method, from the PS in the model of the optical fiber with realistic values of coefficients accounting for the anomalous dispersion, Kerr nonlinearity, and higher-order effects. The results demonstrate that the breathing solitons stably propagate in the fibers. Their robustness against small random perturbations applied to the initial background is demonstrated too.

preprint2014arXiv

Breathing solitary-pulse pairs in a linearly coupled system

It is shown that pairs of solitary pulses (SPs) in a linearly-coupled system with opposite group-velocity dispersions form robust breathing bound states. The system can be realized by temporal-modulation coupling of SPs with different carrier frequencies propagating in the same medium, or by coupling of SPs in a dual-core waveguide. Broad SP pairs are produced in a virtually exact form by means of the variational approximation. Strong nonlinearity tends to destroy the periodic evolution of the SP pairs.

preprint2014arXiv

Building patterns by traveling vortices and dipoles in periodic dissipative media

We analyze pattern-formation scenarios in the two-dimensional (2D) complex Ginzburg-Landau (CGL) equation with the cubic-quintic (CQ) nonlinearity and a cellular potential. The equation models laser cavities with built-in gratings, which are used to stabilize 2D patterns. The pattern-building process is initiated by kicking a localized compound mode, in the form of a dipole, quadrupole, or vortex which is composed of four local peaks. The hopping motion of the kicked mode through the cellular structure leads to the generation of various extended patterns pinned by the structure. In the ring-shaped system, the persisting freely moving dipole hits the stationary pattern from the opposite side, giving rise to several dynamical regimes, with the pinned multi-soliton chain playing the role of the Newton's cradle (NC).

preprint2014arXiv

Bulk vortices and half-vortex surface modes in parity-time-symmetric media

We demonstrate that in-bulk vortex localized modes, and their surface half-vortex ("horseshoe") counterparts self-trap in two-dimensional (2D) nonlinear optical systems with $\mathcal{PT}$-symmetric photonic lattices (PLs). The respective stability regions are identified in the underlying parameter space. The in-bulk states are related to truncated nonlinear Bloch waves in gaps of the PL-induced spectrum. The basic vortex and horseshoe modes are built, severally, of four and three beams with appropriate phase shifts between them. Their stable complex counterparts, built of up to 12 beams, are reported too.

preprint2014arXiv

Collisions of matter-wave solitons

Solitons are localised wave disturbances that propagate without changing shape, a result of a nonlinear interaction which compensates for wave packet dispersion. Individual solitons may collide, but a defining feature is that they pass through one another and emerge from the collision unaltered in shape, amplitude, or velocity. This remarkable property is mathematically a consequence of the underlying integrability of the one-dimensional (1D) equations, such as the nonlinear Schrödinger equation, that describe solitons in a variety of wave contexts, including matter-waves$^{1,2}$. Here we explore the nature of soliton collisions using Bose-Einstein condensates of atoms with attractive interactions confined to a quasi-one-dimensional waveguide. We show by real-time imaging that a collision between solitons is a complex event that differs markedly depending on the relative phase between the solitons. Yet, they emerge from the collision unaltered in shape or amplitude, but with a new trajectory reflecting a discontinuous jump. By controlling the strength of the nonlinearity we shed new light on these fundamental features of soliton collisional dynamics, and explore the implications of collisions that bring the wave packets out of the realm of integrability, where they may undergo catastrophic collapse. 1. Zabusky, N.J. & Kruskal, M.D. Interaction of "solitons" in a collisionless plasma and the recurrence of initial states. Phys. Rev. Lett. 15, 240 (1965). 2. Zakharov, V.E. & Shabat, A.B. Exact theory of two-dimensional self-focusing and one-dimensional self-moduation of waves in nonlinear media. Sov. Phys. JEPT. 34, 62 (1972).

preprint2014arXiv

Discrete solitons and scattering of lattice waves in guiding arrays with a nonlinear $\mathcal{PT}$-symmetric defect

Discrete fundamental and dipole solitons are constructed, in an exact analytical form, in an array of linear waveguides with an embedded $\mathcal{PT}$-symmetric dimer, which is composed of two nonlinear waveguides carrying equal gain and loss. Fundamental solitons in tightly knit lattices, as well as all dipole modes, exist above a finite threshold value of the total power. However, the threshold vanishes for fundamental solitons in loosely knit lattices. The stability of the discrete solitons is investigated analytically by means of the Vakhitov-Kolokolov (VK) criterion, and, in the full form, via the computation of eigenvalues for perturbation modes. Fundamental and dipole solitons tend to be stable at smaller and larger values of the total power (norm), respectively. The increase of the strength of the coupling between the two defect-forming sites leads to strong expansion of the stability areas. The scattering problem for linear lattice waves impinging upon the defect is considered too.

preprint2014arXiv

Discrete solitons and vortices on two-dimensional lattices of $\mathcal{PT}$-symmetric couplers

We introduce a 2D network built of $\mathcal{PT}$-symmetric dimers with on-site cubic nonlinearity, the gain and loss elements of the dimers being linked by parallel square-shaped lattices. The system may be realized as a set of $\mathcal{PT}$-symmetric dual-core waveguides embedded into a photonic crystal. The system supports $\mathcal{PT}$-symmetric and antisymmetric fundamental solitons (FSs) and on-site-centered solitary vortices (OnVs). Stability of these discrete solitons is the central topic of the consideration. Their stability regions in the underlying parameter space are identified through the computation of stability eigenvalues, and verified by direct simulations. Symmetric FSs represent the system's ground state, being stable at lowest values of the power, while anti-symmetric FSs and OnVs are stable at higher powers. Symmetric OnVs, which are also stable at lower powers, are remarkably robust modes: on the contrary to other PT-symmetric states, unstable OnVs do not blow up, but spontaneously rebuild themselves into stable FSs.

preprint2014arXiv

Dissipative quadratic solitons supported by a localized gain

We propose two models for the creation of stable dissipative solitons in optical media with the $χ^{(2)}$ (quadratic) nonlinearity. To compensate spatially uniform loss in both the fundamental-frequency (FF) and second-harmonic (SH) components of the system, a strongly localized "hot spot", carrying the linear gain, is added, acting either on the FF component, or on the SH one. In both systems, we use numerical methods to find families of dissipative $χ^{(2)}$ solitons pinned to the "hot spot". The shape of the existence and stability domains may be rather complex. An existence boundary for the solitons, which corresponds to the guided mode in the linearized version of the systems, is obtained in an analytical form. The solitons demonstrate noteworthy features, such as spontaneous symmetry breaking (of spatially symmetric solitons) and bistability.

preprint2014arXiv

Emulation of lossless exciton-polariton condensates by dual-core optical waveguides: Stability, collective modes, and dark solitons

We propose a possibility to simulate the exciton-polariton (EP) system in the lossless limit, which is not currently available in semiconductor microcavities, by means of a simple optical dual-core waveguide, with one core carrying the nonlinearity and operating close to the zero-group-velocity-dispersion (GVD) point, and the other core being linear and dispersive. Both 2D and 1D EP systems may be emulated by means of this optical setting. In the framework of this system, we find that, while the uniform state corresponding to the lower branch of the nonlinear dispersion relation is stable against small perturbations, the upper branch is always subject to the modulational instability (MI). The stability and instability are verified by direct simulations too. We analyze collective excitations on top of the stable lower-branch state, which include a Bogoliubov-like gapless mode and a gapped one. Analytical results are obtained for the corresponding sound velocity and energy gap. The effect of a uniform phase gradient (superflow) on the stability is considered too, with a conclusion that the lower-branch state becomes unstable above a critical wavenumber of the flux. Finally, we demonstrate that the stable 1D state may carry robust dark solitons.

preprint2014arXiv

Interactions of spatial solitons with fused couplers

We study dynamical and stationary states of solitons in dual-core waveguides which are locally coupled (fused) at one or several short segments. The model applies to planar optical waveguides, and to BEC in dual traps. Collisions of an incident soliton with single and double locally-fused couplers are investigated by means of systematic simulations and several analytical methods (quasi- linear, fast-soliton, and adiabatic approximations). Excitation dynamics of a soliton trapped by a local coupler is studied by means of the variational approximation, and verified by simulations. Shuttle motion of a soliton trapped in a cavity between two local couplers, and in a finite array of couplers, is studied too.

preprint2014arXiv

Localized modes in dissipative lattice media: An overview

We overview recent theoretical studies of the dynamics of one- and two-dimensional spatial dissipative solitons in models based on the complex Ginzburg-Landau equations with the cubic-quintic combination of loss and gain terms, which include imaginary, real, or complex spatially periodic potentials. The imaginary potential represents periodic modulation of the local loss and gain. It is shown that the effective gradient force, induced by the inhomogeneous loss distribution, gives rise to three generic propagation scenarios for one-dimensional (1D) dissipative solitons: transverse drift, persistent swing motion, and damped oscillations. When the lattice-average loss/gain value is zero, and the real potential has spatial parity opposite to that of the imaginary component, the respective complex potential is a realization of the parity-time symmetry. Under the action of lattice potentials of the latter type, 1D solitons feature unique motion regimes in the form of transverse drift and persistent swing. In the 2D geometry, three types of axisymmetric radial lattices are considered, viz., ones based solely on the refractive-index modulation, or solely on the linear-loss modulation, or on a combination of both. The rotary motion of solitons in such axisymmetric potentials can be effectively controlled by varying the strength of the initial tangential kick.

preprint2014arXiv

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.

preprint2014arXiv

Localized modes in quasi-2D Bose-Einstein condensates with spin-orbit and Rabi couplings

We consider a two-component pancake-shaped, i.e., effectively two-dimensional (2D), Bose-Einstein condensate (BEC) coupled by the spin-orbit (SO) and Rabi terms. The SO coupling adopted here is of the mixed Rashba-Dresselhaus type. For this configuration, we derive a system of two 2D nonpolynomial Schrödinger equations (NPSEs), for both attractive and repulsive interatomic interactions. In the low- and high-density limits, the system amounts to previously known models, namely, the usual 2D Gross-Pitaevskii equation, or the Schrödinger equation with the nonlinearity of power 7/3. We present simple approximate localized solutions, obtained by treating the SO and Rabi terms as perturbations. Localized solutions of the full NPSE system are obtained in a numerical form. Remarkably, in the case of the attractive nonlinearity acting in free space (i.e., without any 2D trapping potential), we find parameter regions where the SO and Rabi couplings make 2D fundamental solitons dynamically stable.

preprint2014arXiv

Localized modes in two-dimensional Schroedinger lattices with a pair of nonlinear sites

We address the existence and stability of localized modes in the two-dimensional (2D) linear Schroedinger lattice with two symmetric nonlinear sites embedded into it, and a generalization for moderately localized nonlinearity featuring two distinct symmetric maxima. The latter setting admits a much greater variety of localized modes. Symmetric, antisymmetric, and asymmetric discrete solitons are found, and a subcritical bifurcation, accounting for the spontaneous symmetry breaking (SSB) of the symmetric modes and transition to asymmetric ones, is identified. Existence and stability of more complex 2D solutions in the form of discrete symmetric and asymmetric vortices are also discussed.

preprint2014arXiv

Mini-review: Spatial solitons supported by localized gain

The creation of stable 1D and 2D localized modes in lossy nonlinear media is a fundamental problem in optics and plasmonics. This article gives a short review of theoretical methods elaborated for this purpose, using localized gain applied at one or several "hot spots" (HSs). The introduction surveys a broad class of models for which this approach was developed. Other sections focus in some detail on basic 1D continuous and discrete systems, where the results can be obtained, partly or fully, in an analytical form (verified by comparison with numerical results), which provides a deeper insight into the nonlinear dynamics of optical beams in dissipative nonlinear media. In particular, considered are the single and double HSs in the usual waveguide with the self-focusing (SF) or self-defocusing (SDF) Kerr nonlinearity, which give rise to rather sophisticated results, in spite of apparent simplicity of the model; solitons attached to a PT-symmetric dipole embedded into the SF or SDF medium; gap solitons pinned to an HS in a Bragg grating; and discrete solitons in a 1D lattice with a "hot site".

preprint2014arXiv

Multisoliton Newton's cradles and supersolitons in regular and PT-symmetric nonlinear couplers

We demonstrate the existence of stable collective excitation in the form of "supersolitons" propagating through chains of solitons with alternating signs (i.e., Newton's cradles built of solitons) in nonlinear optical couplers, including the PT-symmetric version thereof. In the regular coupler, stable supersolitons are created in the cradles composed of both symmetric solitons and asymmetric ones with alternating polarities. Collisions between moving supersolitons are investigated too, by the means of direct simulations in both the regular and PT-symmetric couplers.

preprint2014arXiv

Non-standard Hubbard models in optical lattices: a review

Originally, the Hubbard model has been derived for describing the behaviour of strongly-correlated electrons in solids. However, since over a decade now, variations of it are also routinely being implemented with ultracold atoms in optical lattices. We review some of the rich literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so called extended Bose-Hubbard models, that include interactions between neighboring sites, next-neighboring sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation on the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, also the effects related to higher Bloch bands become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher $p$-orbitals, as well as models that couple dynamically spin and orbital degrees of freedom. Finally, we review mean-field nonlinear-Schrödinger models of the Salerno type that share with the non-standard Hubbard models the nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of the consideration. We conclude by listing some future open problems.

preprint2014arXiv

Persistence of the Thomas-Fermi approximation for ground states supported by the nonlinear confinement

We justify the Thomas--Fermi approximation for the elliptic problem with the repulsive nonlinear confinement used in the recent physical literature. The method is based on the resolvent estimates and the fixed-point iterations.

preprint2014arXiv

Pinned modes in two-dimensional lossy lattices with local gain and nonlinearity

We introduce a system with one or two amplified nonlinear sites ("hot spots", HSs) embedded into a two-dimensional linear lossy lattice. The system describes an array of evanescently coupled optical or plasmonic waveguides, with gain applied at selected HS cores. The subject of the analysis is discrete solitons pinned to the HSs. The shape of the localized modes is found in quasi-analytical and numerical forms, using a truncated lattice for the analytical consideration. Stability eigenvalues are computed numerically, and the results are supplemented by direct numerical simulations. In the case of self-focusing nonlinearity, the modes pinned to a single HS are stable or unstable when the nonlinearity includes the cubic loss or gain, respectively. If the nonlinearity is self-defocusing, the unsaturated cubic gain acting at the HS supports stable modes in a small parametric area, while weak cubic loss gives rise to a bistability of the discrete solitons. Symmetric and antisymmetric modes pinned to a symmetric set of two HSs are considered too.

preprint2014arXiv

Quasi-energies, parametric resonances, and stability limits in ac-driven $\mathcal{PT}$-symmetric systems

We introduce a simple model for implementing the concepts of quasi-energy and parametric resonances (PRs) in systems with the $\mathcal{PT}$ symmetry, i.e., a pair of coupled and mutually balanced gain and loss elements. The parametric (ac) forcing is applied through periodic modulation of the coefficient accounting for the coupling of the two degrees of freedom. The system may be realized in optics as a dual-core waveguide with the gain and loss applied to different cores, and the thickness of the gap between them subject to a periodic modulation. The onset and development of the parametric instability for a small forcing amplitude ($V_{1}$) is studied in an analytical form. The full dynamical chart of the system is generated by systematic simulations. At sufficiently large values of the forcing frequency, $ω$, tongues of the parametric instability originate, with the increase of $V_{1}$, as predicted by the analysis. However, the tongues following further increase of $V_{1}$ feature a pattern drastically different from that in usual (non-$\mathcal{PT}$)\ parametrically driven systems: instead of bending down to larger values of the dc coupling constant, $V_{0}$, they maintain a direction parallel to the $V_{1}$ axis. The system of the parallel tongues gets dense with the decrease of $ω$, merging into a complex small-scale structure of alternating regions of stability and instability. The cases of $ω\rightarrow 0$ and $ω\rightarrow \infty $ are studied analytically by means of the adiabatic and averaging approximation, respectively. The cubic nonlinearity, if added to the system, alters the picture, destabilizing many originally robust dynamical regimes, and stabilizing some which were unstable.

preprint2014arXiv

Spatial control of the competition between self-focusing and self-defocusing nonlinearities in one- and two-dimensional systems

We introduce a system with competing self-focusing (SF) and self-defocusing (SDF) terms, which have the same scaling dimension. In the one-dimensional (1D) system, this setting is provided by a combination of the SF cubic term multiplied by the delta-function, $δ(x)$, and a spatially uniform SDF quintic term. This system gives rise to the most general family of 1D-Townes solitons, the entire family being unstable. However, it is completely stabilized by a finite-width regularization of the $δ$-function. The results are produced by means of numerical and analytical methods. We also consider the system with a symmetric pair of regularized $δ$-functions, which gives rise to a wealth of symmetric, antisymmetric, and asymmetric solitons, linked by a bifurcation loop, that accounts for the breaking and restoration of the symmetry. Soliton families in 2D versions of both the single- and double-delta-functional systems are also studied. The 1D and 2D settings may be realized for spatial solitons in optics, and in Bose-Einstein condensates.

preprint2014arXiv

Spontaneous symmetry breaking in dual-core baby-Skyrmion systems

We introduce a system composed of two (2+1)-dimensional baby-skyrmion models (BSMs) set on parallel planes and linearly coupled by tunneling of fields. This system can be realized in a dual-layer ferromagnetic medium. Unlike dual-core models previously studied\ in nonlinear optics and BEC, here the symmetry-breaking bifurcation (SBB) in solitons (baby skyrmions) occurs with the increase of the inter-core coupling ($κ$), rather than with its decrease, due to the fact that even in the uncoupled system neither core may be empty. Prior to the onset of the symmetry breaking between the two components of the solitons, they gradually separate in the lateral direction, due to the increase of $κ$, which is explained in an analytical form by means of an effective interaction potential. Such evolution scenario are produced for the originally symmetric states with topological charges in the two cores $\left(B^{(1)},B^{(2)}\right)$, with $% B^{(1)}=B^{(2)}=1,2,3,4$. The evolution of mixed states, of the $\left(1,2\right)$ and $\left(2,4\right)$ types, with the variation of $κ$ is studied too.

preprint2014arXiv

The creation of two-dimensional composite solitons in spin-orbit-coupled self-attractive Bose-Einstein condensates in free space

It is commonly known that two-dimensional mean-field models of optical and matter waves with the cubic self-attraction cannot produce stable solitons in free space because of the occurrence of the collapse in the same setting. By means of the numerical analysis and variational approximation, we demonstrate that the two-component model of the Bose-Einstein condensate, with the spin-orbit Rashba coupling and cubic attractive interactions, gives rise to solitary-vortex complexes of two types: semi-vortices (SVs, with a vortex in one component and a fundamental soliton in the other), and mixed modes (MMs, with topological charges 0 and +1/-1 mixed in both components). These two- dimensional composite modes can be created using the trapping harmonic-oscillator (HO) potential, but remain stable in the free space, if the trap is gradually removed. The SVs and MMs realize the ground state of the system, provided that the self-attraction in the two components is, respectively, stronger or weaker than the cross-attraction between them. The SVs and MMs which are not the ground state are subject to a drift instability. In the free space (in the absence of the HO trap), modes of both types degenerate into unstable Townes solitons when their norms attain the respective critical values, while there is no lower existence threshold for the stable modes. Moving free-space stable solitons are also found in the present non-Galilean-invariant system, up to a critical velocity. Collisions between two moving solitons lead to their merger into a single one.

preprint2014arXiv

Three-dimensional hybrid vortex solitons

We show, by means of numerical and analytical methods, that media with a repulsive nonlinearity which grows from the center to the periphery support a remarkable variety of previously unknown complex stationary and dynamical three-dimensional solitary-wave states. Peanut-shaped modulation profiles give rise to vertically symmetric and antisymmetric vortex states, and novel stationary hybrid states, built of top and bottom vortices with opposite topological charges, as well as robust dynamical hybrids, which feature stable precession of a vortex on top of a zero-vorticity base. The analysis reveals stability regions for symmetric, antisymmetric, and hybrid states. In addition, bead-shaped modulation profiles give rise to the first example of exact analytical solutions for stable three-dimensional vortex solitons. The predicted states may be realized in media with a controllable cubic nonlinearity, such as Bose-Einstein condensates.

preprint2014arXiv

Unbreakable PT-symmetry of solitons supported by inhomogeneous defocusing nonlinearity

We consider bright solitons supported by a symmetric inhomogeneous defocusing nonlinearity growing rapidly enough toward the periphery of the medium, combined with an antisymmetric gain-loss profile. Despite the absence of any symmetric modulation of the linear refractive index, which is usually required to establish a PT-symmetry in the form of a purely real spectrum of modes, we show that the PT-symmetry is never broken in the present system, and that the system always supports stable bright solitons, fundamental and multi-pole ones. Such phenomenon is connected to non-linearizability of the underlying evolution equation. The increase of the gain-losses strength results, in lieu of the PT-symmetry breaking, in merger of pairs of different soliton branches, such as fundamental and dipole, or tripole and quadrupole ones. The fundamental and dipole solitons remain stable for all values of the gain-loss coefficient.

preprint2014arXiv

Unveiling the orbital angular momentum and acceleration of electron beams

New forms of electron beams have been intensively investigated recently, including vortex beams carrying orbital angular momentum, as well as Airy beams propagating along a parabolic trajectory. Their traits may be harnessed for applications in materials science, electron microscopy and interferometry, and so it is important to measure their properties with ease. Here we show how one may immediately quantify these beams' parameters without need for additional fabrication or non-standard microscopic tools. Our experimental results are backed by numerical simulations and analytic derivation.

preprint2013arXiv

A solvable model for solitons pinned to a PT-symmetric dipole

We introduce the simplest one-dimensional nonlinear model with the parity-time (PT) symmetry, which makes it possible to find exact analytical solutions for localized modes ("solitons"). The PT-symmetric element is represented by a point-like (delta-functional) gain-loss dipole δ^{\prime}(x), combined with the usual attractive potential δ(x). The nonlinearity is represented by self-focusing (SF) or self-defocusing (SDF) Kerr terms, both spatially uniform and localized ones. The system can be implemented in planar optical waveguides. For the sake of comparison, also introduced is a model with separated δ-functional gain and loss, embedded into the linear medium and combined with the δ-localized Kerr nonlinearity and attractive potential. Full analytical solutions for pinned modes are found in both models. The exact solutions are compared with numerical counterparts, which are obtained in the gain-loss-dipole model with the δ^{\prime}- and δ- functions replaced by their Lorentzian regularization. With the increase of the dipole's strength, γ, the single-peak shape of the numerically found mode, supported by the uniform SF nonlinearity, transforms into a double-peak one. This transition coincides with the onset of the escape instability of the pinned soliton. In the case of the SDF uniform nonlinearity, the pinned modes are stable, keeping the single-peak shape.

preprint2013arXiv

Bright solitons in quasi-one dimensional dipolar condensates with spatially modulated interactions

We introduce a model for the condensate of dipolar atoms or molecules, in which the dipole-dipole interaction (DDI) is periodically modulated in space, due to a periodic change of the local orientation of the permanent dipoles, imposed by the corresponding structure of an external field (the necessary field can be created, in particular, by means of magnetic lattices, which are available to the experiment). The system represents a realization of a nonlocal nonlinear lattice, which has a potential to support various spatial modes. By means of numerical methods and variational approximation (VA), we construct bright one-dimensional solitons in this system, and study their stability. In most cases, the VA provides good accuracy, and correctly predicts the stability by means of the Vakhitov-Kolokolov (VK)\ criterion. It is found that the periodic modulation may destroy some solitons, which exist in the usual setting with unmodulated DDI, and can create stable solitons in other cases, not verified in the absence of modulations. Unstable solitons typically transform into persistent localized breathers. The solitons are often mobile, with inelastic collisions between them leading to oscillating localized modes.

preprint2013arXiv

Buffering and Trapping Ultrashort Optical Pulses in Concatenated Bragg Gratings

Strong retardation of ultrashort optical pulses, including their deceleration and stoppage in the form of Bragg solitons in a cascaded Bragg-grating (BG) structure, is proposed. The manipulations of the pulses are carried out, using nonlinear effects, in a chirped BG segment which is linked, via a defect, to a uniform grating. The storage of the ultrashort pulses is shown to be very robust with respect to variations of the input field intensity, suggesting the feasibility of storing ultrafast optical pulses in such a structure. Physical estimates are produced for the BGs written in silicon.

preprint2013arXiv

Discrete solitons in an array of quantum dots

We develop a theory for the interaction of classical light fields with an a chain of coupled quantum dots (QDs), in the strong-coupling regime, taking into account the local-field effects. The QD chain is modeled by a one-dimensional (1D) periodic array of two-level quantum particles with tunnel coupling between adjacent ones. The local-field effect is taken into regard as QD depolarization in the Hartree-Fock-Bogoliubov approximation. The dynamics of the chain is described by a system of two discrete nonlinear Schrödinger (DNLS) equations for local amplitudes of the probabilities of the ground and first excited states. The two equations are coupled by a cross-phase-modulation cubic terms, produced by the local-field action, and by linear terms too. In comparison with previously studied DNLS systems, an essentially new feature is a phase shift between the intersite-hopping constants in the two equations. By means of numerical solutions, we demonstrate that, in this QD chain, Rabi oscillations (RO) self-trap into stable bright\textit{\ Rabi solitons} or \textit{Rabi breathers}. Mobility of the solitons is considered too. The related behavior of observable quantities, such as energy, inversion, and electric-current density, is given a physical interpretation. The results apply to a realistic region of physical parameters.

preprint2013arXiv

Generation of tightly compressed solitons with a tunable frequency shift in Raman-free fibers

Optimization of the compression of input N-solitons into robust ultra-narrow fundamental solitons, with a tunable up- or downshifted frequency, is proposed in photonic crystal fibers (PCF) free of the Raman effect. Due to the absence of the Raman self-frequency shift, these fundamental solitons continue propagation maintaining the acquired frequency, once separated from the input N-soliton's temporal slot. A universal optimal value of the relative strength of the third-order-dispersion (TOD) is found, providing the strongest compression of the fundamental soliton, is found. It depends only the order of the injected N-soliton. The largest compression degree significantly exceeds the analytical prediction supplied by the Satsuma-Yajima (SY) formula. The mechanism behind this effect, which remains valid in the presence of the self-steepening, is explained.

preprint2013arXiv

High- and low-frequency phonon modes in dipolar quantum gases trapped in deep lattices

We study normal modes propagating on top of the stable uniform background in arrays of dipolar Bose-Einstein condensate (BEC) droplets trapped in a deep optical lattice. Both the on-site mean-field dynamics of the droplets and their displacement due to the repulsive dipole-dipole interactions (DDIs) are taken into account. Dispersion relations for two modes, \textit{viz}., high- and low- frequency counterparts of optical and acoustic phonon modes in condensed matter, are derived analytically and verified by direct simulations, for both cases of the repulsive and attractive contact interactions. The (counterpart of the) optical-phonon branch does not exist without the DDIs. These results are relevant in the connection to emerging experimental techniques enabling real-time imaging of the condensate dynamics and direct experimental measurement of phonon dispersion relations in BECs.

preprint2013arXiv

Holding spatial solitons in a pumped cavity with the help of nonlinear potentials

We introduce a one-dimensional model of a cavity with the Kerr nonlinearity and saturated gain, designed so as to keep solitons in the state of shuttle motion. The solitons are always unstable in the cavity bounded by the usual potential barriers, due to accumulation of noise generated by the linear gain. Complete stabilization of the shuttling soliton is achieved if the linear barrier potentials are replaced by nonlinear ones, which trap the soliton, being transparent to the radiation. The removal of the noise from the cavity is additionally facilitated by an external ramp potential. The stable dynamical regimes are found numerically, and their basic properties are explained analytically.

preprint2013arXiv

Inversion and tight focusing of Airy pulses under the action of third-order dispersion

By means of direct simulations and theoretical analysis, we study the nonlinear propagation of truncated Airy pulses in an optical fiber exhibiting both anomalous second-order and strong positive third-order dispersions. It is found that the Airy pulse first reaches a finite-size focal area as determined by the relative strength of the two dispersion terms, and then undergoes an inversion transformation such that it continues to travel with an opposite acceleration. The system notably features tight focusing if the third-order dispersion is a dominant factor. These effects are partially reduced by a Kerr nonlinearity.

preprint2013arXiv

Lattice solitons with quadrupolar intersite interactions

We study two-dimensional (2D) solitons in the mean-field models of ultracold gases with long-range quadrupole-quadrupole interaction (QQI) between particles. The condensate is loaded into a deep optical-lattice (OL) potential, therefore the model is based on the 2D discrete nonlinear Schrödinger equation with contact onsite and long-range intersite interactions, which represent the QQI. The quadrupoles are built as pairs of electric dipoles and anti-dipoles orientated perpendicular to the 2D plane to which the gas is confined. Because the quadrupoles interact with the local gradient of the external field, they are polarized by inhomogeneous dc electric field that may be supplied by a tapered capacitor. Shapes, stability, mobility, and collisions of fundamental discrete solitons are studied by means of systematic simulations. In particular, threshold values of the norm, necessary for the existence of the solitons, are found, and anisotropy of their static and dynamical properties is explored. As concerns the mobility and collisions, it is the first analysis of such properties for discrete solitons on 2D lattices with long-range intersite interactions of any type. Estimates demonstrate that the setting can be realized under experimentally available conditions, predicting solitons built of $\sim$ 10,000 particles.

preprint2013arXiv

Localized modes in dense repulsive and attractive Bose-Einstein condensates with spin-orbit and Rabi couplings

We consider a binary Bose-Einstein condensate with linear and nonlinear interactions between its components, which emulate the spinor system with spin-orbit (SO) and Rabi couplings. For a relatively dense condensate, 1D coupled equations with the nonpolynomial nonlinearity of both repulsive and attractive signs are derived from the 3D Gross-Pitaevskii equations. Profiles of modes confined in an external potential under the action of the self-repulsion, and self-trapped solitons in the case of the self-attraction, are found in a numerical form and by means of analytical approximations. In the former case, the interplay of the SO and Rabi couplings with the repulsive nonlinearity strongly distorts shapes of the trapped modes, adding conspicuous side lobes to them. In the case of the attractive nonlinearity, the most essential result is reduction of the collapse threshold under the action of the SO and Rabi couplings.

preprint2013arXiv

Matter-wave solitons supported by field-induced dipole-dipole repulsion with a spatially modulated strength

We demonstrate the existence of one and two-dimensional bright solitons in the Bose-Einstein condensate with repulsive dipole-dipole interactions induced by a combination of dc and ac polarizing fields, oriented perpendicular to the plane in which the BEC is trapped, assuming that the strength of the fields grows in the radial ($r$) direction faster than $r^{3}$. Stable tightly confined 1D and 2D fundamental solitons, twisted solitons in 1D, and solitary vortices in 2D are found in a numerical form. The fundamental solitons remain robust under the action of an expulsive potential, which is induced by the interaction of the dipoles with the polarizing field. The confinement and scaling properties of the soliton families are explained analytically. The Thomas-Fermi approximation is elaborated for fundamental solitons. The mobility of the fundamental solitons is limited to the central area. Stable 1D even and odd solitons are also found in the setting with a double-well modulation function, along with a regime of Josephson oscillations.

preprint2013arXiv

Mobility of solitons in one-dimensional lattices with the cubic-quintic nonlinearity

We investigate mobility regimes for localized modes in the discrete nonlinear Schrödinger (DNLS) equation with the cubic-quintic onsite terms. Using the variational approximation (VA), the largest soliton's total power admitting progressive motion of kicked discrete solitons is predicted, by comparing the effective kinetic energy with the respective Peierls-Nabarro (PN) potential barrier. The prediction is novel for the DNLS model with the cubic-only nonlinearity too, demonstrating a reasonable agreement with numerical findings. Small self-focusing quintic term quickly suppresses the mobility. In the case of the competition between the cubic self-focusing and quintic self-defocusing terms, we identify parameter regions where odd and even fundamental modes exchange their stability, involving intermediate asymmetric modes. In this case, stable solitons can be set in motion by kicking, so as to let them pass the PN barrier. Unstable solitons spontaneously start oscillatory or progressive motion, if they are located, respectively, below or above a mobility threshold. Collisions between moving discrete solitons, at the competing nonlinearities frame, are studied too.

preprint2013arXiv

Photon-atomic solitons in a Bose-Einstein condensate trapped in a soft optical lattice

We investigate the ground state (GS) of a collisionless Bose-Einstein condensate (BEC) trapped in a soft one-dimensional optical lattice (OL), which is formed by two counterpropagating optical beams perturbed by the BEC density profile through the local-field effect (LFE). We show that LFE gives rise to an envelope-deformation potential, a nonlocal potential resulting from the phase deformation, and an effective self-interaction of the condensate. As a result, stable photon-atomic lattice solitons, including an optical component, in the form of the deformation of the soft OL, in a combination with a localized matter-wave component, are generated in the blue-detuned setting, without any direct interaction between atoms. These self-trapped modes, which realize the system's GS, are essentially different from the gap solitons supported by the interplay of the OL potential and collisional interactions between atoms. A transition to tightly bound modes from loosely bound ones occurs with the increase of the number of atoms in the BEC.

preprint2013arXiv

Self-trapping of Fermi and Bose gases under spatially modulated repulsive nonlinearity and transverse confinement

We show that self-localized ground states can be created in the spin-balanced gas of fermions with repulsion between the spin components, whose strength grows from the center to periphery, in combination with the harmonic-oscillator (HO) trapping potential acting in one or two transverse directions. We also consider the ground state in the non-interacting Fermi gas under the action of the spatially growing tightness of the one- or two-dimensional (1D or 2D) HO confinement. These settings are considered in the framework of the Thomas-Fermi-von Weizsacker (TF-vW) density functional. It is found that the vW correction to the simple TF approximation (the gradient term) is nearly negligible in all situations. The properties of the ground state under the action of the 2D and 1D HO confinement with the tightness growing in the transverse directions is investigated too for the Bose-Einstein condensate (BEC) with the self-repulsive nonlinearity.

preprint2013arXiv

Solitary modes in nonlocal media with inhomogeneous self-repulsive nonlinearity

We demonstrate the existence of two species of stable bright solitons, fundamental and dipole ones, in one-dimensional self-defocusing nonlocal media, with the local value of nonlinearity coefficient having one or several minima and growing at any rate faster than |x| at large values of coordinate x. The model can be derived for a slab optical waveguide with thermal nonlinearity. The most essential difference from the local counterpart of this system is the competition between two different spatial scales, the one determining the modulation pattern of the nonlinearity coefficient, and the correlation length of the nonlocality. The competition is explicitly exhibited by analytically obtained asymptotic form of generic solutions. Particular exact solutions are found analytically, and full soliton families are constructed in a numerical form. The multi-channel settings, with two or three local minima of the nonlinearity coefficient, are considered here for the first time, for both local and nonlocal models of the present type. States with multiple solitons launched into different channels are stable if the spacing between them exceeds a certain minimum value. A regime of stable Josephson oscillations of solitons between parallel channels is reported too.

preprint2013arXiv

Solitary vortices supported by localized parametric gain

We demonstrate the existence and stability of bright vortex solitons sustained by a ring-shaped parametric gain, for both focusing and defocusing Kerr nonlinearities in lossy optical media. With the defocusing nonlinearity, the vortices are stable at all values of the detuning parameter, while under the focusing nonlinearity their stability region is limited to some positive values of the detuning. Unstable vortices in the focusing medium transform into stable rotating azimuthons.

preprint2013arXiv

Soliton gyroscopes in media with spatially growing repulsive nonlinearity

We find that the recently introduced model of self-trapping supported by a spatially growing strength of a repulsive nonlinearity gives rise to robust vortex-soliton tori, i.e., three-dimensional vortex solitons, with topological charges S. The family with S=1 is completely stable, while the one with S=2 has alternating regions of stability and instability. The families are nearly exactly reproduced in an analytical form by the Thomas-Fermi approximation (TFA). Unstable states with S=2 and 3 split into persistently rotating pairs or triangles of unitary vortices. Application of a moderate torque to the vortex torus initiates a persistent precession mode, with the torus' axle moving along a conical surface. A strong torque heavily deforms the vortex solitons, but, nonetheless, they restore themselves with the axle oriented according to the vectorial addition of angular momenta.

preprint2013arXiv

Solitons supported by localized parametric gain

We address the existence and properties of one-dimensional solitons maintained by localized parameter gain in focusing and defocusing lossy nonlinear media. Localized parametric gain supports both fundamental and multipole solitons. We found that the family of fundamental solitons is partly stable in focusing nonlinear medium, and completely stable in defocusing medium, while all higher-order solitons are unstable. In addition to numerical results, the existence threshold for the solitons, and a particular stable exact solution are obtained in an exact analytical form.

preprint2013arXiv

Stability boundaries and collisions of two-dimensional solitons in PT-symmetric couplers with the cubic-quintic nonlinearity

We introduce one- and two-dimensional (1D and 2D) models of parity-time ($% \mathcal{PT}$) -symmetric couplers with the mutually balanced linear gain and loss applied to the two cores, and cubic-quintic (CQ) nonlinearity acting in each one. The 2D and 1D models may be realized in dual-core optical waveguides, in the spatiotemporal and spatial domains, respectively. Stationary solutions for $\mathcal{PT}$-symmetric solitons in these systems reduce to their counterparts in the usual coupler. The most essential problem is the stability of the solitons, which become unstable against symmetry breaking with the increase of the energy (norm), and retrieve the stability at still larger energies. The boundary value of the intercore-coupling constant, above which the solitons are completely stable, is found by means of an analytical approximation, based on the CW (zero-dimensional) counterpart of the system. The approximation demonstrates good agreement with numerical findings for the 1D and 2D solitons. Numerical results for the stability limits of the 2D solitons are obtained by means of the computation of eigenvalues for small perturbations, and verified in direct simulations. Although large parts of the solitons families are unstable, the instability is quite weak. Collisions between 2D solitons in the $\mathcal{PT}$-symmetric coupler are studied by means of simulations. Outcomes of the collisions are inelastic but not destructive, as they do not break the $\mathcal{PT}$ symmetry.

preprint2013arXiv

Suppression of the quantum collapse in binary bosonic gases

Attraction of the quantum particle to the center in the 3D space with potential V/r^2 gives rise to the quantum collapse, i.e., nonexistence of the ground state (GS) when the attraction strength exceeds a critical value (V = 1/8, in the present notation). Recently, we have demonstrated that the quantum collapse is suppressed, and the GS is restored, if repulsive interactions between particles in the quantum gas are taken into account, in the mean-field approximation. This setting can be realized in a gas of dipolar molecules attracted to the central charge, with dipole-dipole interactions taken into regard too. Here we analyze this problem for a binary gas. GSs supported by the repulsive interactions are constructed in a numerical form, as well as by means of analytical approximations for both miscible and immiscible binary systems. In particular, the Thomas-Fermi (TF) approximation is relevant if V is large enough. It is found that the GS of the miscible binary gas, both balanced and imbalanced, features a weak phase transition at another critical value, V = 1/2. The transition is characterized by an analyticity-breaking change in the structure of the wave functions at small r. To illustrate the generic character of the present phenomenology, we also consider the binary system with the attraction between the species (rather than repulsion), in the case when the central potential pulls a single component only.

preprint2013arXiv

Symmetric and antisymmetric nonlinear modes supported by dual local gain in lossy lattices

We introduce a discrete lossy system, into which a double hot spot (HS) is inserted, i.e., two mutually symmetric sites carrying linear gain and cubic nonlinearity. The system can be implemented as an array of optical or plasmonic waveguides, with a pair of amplified nonlinear cores embedded into it. We focus on the case of the self-defocusing nonlinearity and cubic losses acting at the HSs. Symmetric localized modes pinned to the double HS are constructed in an implicit analytical form, which is done separately for the cases of odd and even numbers of intermediate sites between the HSs. In the former case, some stationary solutions feature a W-like shape, with a low peak at the central site, added to tall peaks at the positions of the embedded HSs. The special case of two adjacent HSs is considered too. Stability of the solution families against small perturbations is investigated in a numerical form, which reveals stable and unstable subfamilies. The instability of symmetric modes accounting for by an isolated positive eigenvalue leads to their spontaneous transformation into co-existing stable antisymmetric modes, while the instability represented by a pair of complex-conjugate eigenvalues gives rise to persistent breathers.

preprint2013arXiv

Symmetric and asymmetric solitons in dual-core couplers with competing quadratic and cubic nonlinearities

We consider the model of a dual-core spatial-domain coupler with chi^(2) and chi^(3) nonlinearities acting in two parallel cores. We construct families of symmetric and asymmetric solitons in the system with self-defocusing chi^(3) terms, and test their stability. The transition from symmetric to asymmetric soliton branches, and back to the symmetric ones proceeds via a bifurcation loop. A pair of stable asymmetric branches emerge from the symmetric family via a supercritical bifurcation; eventually, the asymmetric branches merge back into the symmetric one through a reverse bifurcation. The existence of the loop is explained by means of an extended version of the cascading approximation for the chi^(2) interaction, which takes into regard the XPM part of the chi(3) interaction. When the inter-core coupling is weak, the bifurcation loop features a concave shape, with the asymmetric branches losing their stability at the turning points. In addition to the two-color solitons, which are built of the fundamental-frequency (FF) and second-harmonic (SH) components, in the case of the self-focusing chi^(3) nonlinearity we also consider single-color solitons, which contain only the SH component, but may be subject to the instability against FF perturbations. Asymmetric single-color solitons are always unstable, whereas the symmetric ones are stable, provided that they do not coexist with two-color counterparts. Collisions between tilted solitons are studied too.

preprint2013arXiv

Symmetry breaking and restoring wave transmission in diode-antidiode double chains

We introduce a system of two parallel-coupled discrete nonlinear Schroedinger (DNLS) inhomogeneous chains. Each one favors the unidirectional transmission of incident packets, in the opposite directions with respect to each other. Two different configurations of the diode-antidiode pair are considered, a ladder and a plaquette. They feature, respectively, the uniform transverse linear coupling, or the coupling limited to the central nonlinear segment of the system. In the case of weak linear coupling, the symmetry breaking is observed (i.e., the pair still features the diode behavior), while the moderately strong coupling restores the symmetry, making the transmission effectively bidirectional. In the case of the ladder, an oscillatory dependence of the transmission on the strength of the coupling is observed and qualitatively explained.

preprint2013arXiv

Two-dimensional intra-band solitons in lattice potentials with local defects and self-focusing nonlinearity

It is commonly known that stable bright solitons in periodic potentials, which represent gratings in photonics/plasmonics, or optical lattices in quantum gases, exist either in the spectral semi-infinite gap (SIG) or in finite bandgaps. Using numerical methods, we demonstrate that, under the action of the cubic self-focusing nonlinearity, defects in the form of "holes" in two-dimensional (2D) lattices support continuous families of 2D solitons \textit{embedded} into the first two Bloch bands of the respective linear spectrum, where solitons normally do not exist. The two families of the \textit{embedded defect solitons} (EDSs) are found to be continuously linked by the branch of \textit{gap defect solitons} (GDSs) populating the first finite bandgap. Further, the EDS branch traversing the first band links the GDS family with the branch of regular defect-supported solitons populating the SIG. Thus, we construct a continuous chain of regular, embedded, and gap-mode solitons ("superfamily") threading the entire bandgap structure considered here. The EDSs are stable in the first Bloch band, and partly stable in the second. They exist with the norm exceeding a minimum value, hence they do not originate from linear defect modes. Further, we demonstrate that double, triple and quadruple lattice defects support stable dipole-mode solitons and vortices.

preprint2013arXiv

Vortical light bullets in second-harmonic-generating media supported by a trapping potential

We introduce a three-dimensional (3D) model of optical media with the quadratic ($χ^{(2)}$) nonlinearity and an effective 2D isotropic harmonic-oscillator (HO) potential. While it is well known that 3D χ^2 solitons with embedded vorticity ("vortical light bullets") are unstable in the free space, we demonstrate that they have a broad stability region in the present model, being supported by the HO potential against the splitting instability. The shape of the vortical solitons may be accurately predicted by the variational approximation (VA). They exist above a threshold value of the total energy (norm) and below another critical value, which determines a stability boundary. The existence threshold vanishes is a part of the parameter space, depending on the mismatch parameter, which is explained by means of the comparison with the 2D counterpart of the system. Above the stability boundary, the vortex features shape oscillations, periodically breaking its axisymmetric form and restoring it. Collisions between vortices moving in the longitudinal direction are studied too. The collision is strongly inelastic at relatively small values of the velocities, breaking the two colliding vortices into three, with the same vorticity. The results suggest a possibility of the creation of stable 3D optical solitons with the intrinsic vorticity.

preprint2012arXiv

Asymmetric solitons and domain walls supported by inhomogeneous defocusing nonlinearity

We show that an inhomogeneous defocusing nonlinearity that grows toward the periphery in the positive and negative transverse directions at different rates can support strongly asymmetric fundamental and multipole bright solitons, which are stable in wide parameter regions. In the limiting case when nonlinearity is uniform in one direction, solitons transform into stable domain walls (fronts), with constant or oscillating intensity in the homogeneous region, attached to a tail rapidly decaying in the direction of growing nonlinearity.

preprint2012arXiv

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

It has been recently demonstrated that self-defocusing (SDF) media with the cubic nonlinearity, whose local coefficient grows from the center to 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 $% g_{0}+\sinh ^{2}(β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 $\mathrm{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 "anti-VK" stability criterion for fundamental bright solitons in systems with SDF nonlinearities holds here too. Particular exact solutions for 2D solitons are produced as well.

preprint2012arXiv

Families of solitons in Bragg supergratings

We study fundamental optical gap solitons in the model of a fiber Bragg grating (BG), which is subjected to a periodic modulation of the local reflectivity, giving rise to a supergrating. In addition, the local refractive index is also periodically modulated with the same period. It is known that the supergrating opens an infinite system of new bandgaps in the BG's spectrum. We use a combination of analytical and computational methods to show that each emerging bandgap is filled with gap solitons (GSs), including asymmetric ones and bound states of the GSs. In particular, bifurcations of the GSs created by the supergrating are studied in terms of a geometric analysis.

preprint2012arXiv

Matter-wave solitons with the minimum number of particles in two-dimensional quasiperiodic potentials

We report results of systematic numerical studies of 2D matter-wave soliton families supported by an external potential, in a vicinity of the junction between stable and unstable branches of the families, where the norm of the solution attains a minimum, facilitating the creation of the soliton. The model is based on the Gross-Pitaevskii equation for the self-attractive condensate loaded into a quasiperiodic (QP) optical lattice (OL). The same model applies to spatial optical solitons in QP photonic crystals. Dynamical properties and stability of the solitons are analyzed with respect to variations of the depth and wavenumber of the OL. In particular, it is found that the single-peak solitons are stable or not in exact accordance with the Vakhitov-Kolokolov (VK) criterion, while double-peak solitons, which are found if the OL wavenumber is small enough, are always unstable against splitting.

preprint2012arXiv

Nonlinear modes and symmetry breaking in rotating double-well potentials

We study modes trapped in a rotating ring carrying the self-focusing (SF) or defocusing (SDF) cubic nonlinearity and double-well potential $\cos^{2}θ$, where $θ$ is the angular coordinate. The model, based on the nonlinear Schrödinger (NLS) equation in the rotating reference frame, describes the light propagation in a twisted pipe waveguide, as well as in other optical settings, and also a Bose-Einstein condensate (BEC)trapped in a torus and dragged by the rotating potential. In the SF and SDF regimes, five and four trapped modes of different symmetries are found, respectively. The shapes and stability of the modes, and transitions between them are studied in the first rotational Brillouin zone. In the SF regime, two symmetry-breaking transitions are found, of subcritical and supercritical types. In the SDF regime, an antisymmetry-breaking transition occurs. Ground-states are identified in both the SF and SDF systems.

preprint2012arXiv

Nonlinear PT-symmetric plaquettes

We introduce four basic two-dimensional (2D) plaquette configurations with onsite cubic nonlinearities, which may be used as building blocks for 2D PT -symmetric lattices. For each configuration, we develop a dynamical model and examine its PT symmetry. The corresponding nonlinear modes are analyzed starting from the Hamiltonian limit, with zero value of the gain-loss coefficient. Once the relevant waveforms have been identified (chiefly, in an analytical form), their stability is examined by means of linearization in the vicinity of stationary points. This reveals diverse and, occasionally, fairly complex bifurcations. The evolution of unstable modes is explored by means of direct simulations. In particular, stable localized modes are found in these systems, although the majority of identified solutions is unstable.

preprint2012arXiv

Pattern formation by kicked solitons in the two-dimensionnal Ginzburg-Landau medium with a transverse grating

We consider the kick-induced mobility of two-dimensional (2D) fundamental dissipative solitons in models of lasing media based on the 2D complex Ginzburg-Landau (CGL) equation including a spatially periodic potential (transverse grating). The depinning threshold is identified by means of systematic simulations, and described by means of an analytical approximation, depending on the orientation of the kick. Various pattern-formation scenarios are found above the threshold. Most typically, the soliton, hopping between potential cells, leaves arrayed patterns of different sizes in its wake. In the laser cavity, this effect may be used as a mechanism for selective pattern formation controlled by the tilt of the seed beam. Freely moving solitons feature two distinct values of the established velocity. Elastic and inelastic collisions between free solitons and pinned arrayed patterns are studied too.

preprint2012arXiv

Pinned modes in lossy lattices with local gain and nonlinearity

We introduce a discrete linear lossy system with an embedded "hot spot" (HS), i.e., a site carrying linear gain and complex cubic nonlinearity. The system can be used to model an array of optical or plasmonic waveguides, where selective excitation of particular cores is possible. Localized modes pinned to the HS are constructed in an implicit analytical form, and their stability is investigated numerically. Stability regions for the modes are obtained in the parameter space of the linear gain and cubic gain/loss. An essential result is that the interaction of the unsaturated cubic gain and self-defocusing nonlinearity can produce stable modes, although they may be destabilized by finite amplitude perturbations. On the other hand, the interplay of the cubic loss and self-defocusing gives rise to a bistability.

preprint2012arXiv

Propagation and collisions of semi-discrete solitons in arrayed and stacked waveguides

We consider shapes and dynamics of semi-discrete solitons (SDSs) in the known model of the set of linearly-coupled waveguides with the intrinsic cubic nonlinearity. The model applies to the description of a planar array of optical fibers, or of a stack of parallel planar waveguides, in the temporal and spatial domains, respectively, as well as to the self-attractive Bose-Einstein condensate (BEC) loaded into an array of parallel tunnel-coupled cigar-shaped traps. It was found previously that the interplay of the group-velocity dispersion, discrete diffraction (in the longitudinal and transverse directions, respectively) and intrinsic self-focusing gives rise to SDSs in the array. We here develop the variational approximation (VA) and additional analytical methods for the description of the SDSs, and study their mobility and collisions by means of systematic simulations. The VA and an exact solution of the linearized equation in the cores adjacent to the central one produce an accurate description for the family of stable fundamental onsite-centered SDS solutions, as well of surface SDSs in the semi-infinite array. The VA is also presented for transversely unstable intersite-centered solitons. In simulations, the solitons are not mobile in the discrete direction (non-soliton semi-discrete modes may be mobile across the array). Collisions between SDSs traveling in the longitudinal direction feature a threshold separating the passage and merger or destruction. The exact shape of the threshold, considered as a function of the solitons' energy, features irregularities, while its average form is explained analytically. "Shifted" collisions between SDSs centered at adjacent cores are studied too.

preprint2012arXiv

Quenched dynamics of two-dimensional solitons and vortices in the Gross-Pitaevskii equation

We consider a two-dimensional (2D) counterpart of the experiment that led to the creation of quasi-1D bright solitons in Bose-Einstein condensates (BECs) [Nature 417, 150--153 (2002)]. We start by identifying the ground state of the 2D Gross-Pitaevskii equation for repulsive interactions, with a harmonic-oscillator (HO) trap, and with or without an optical lattice (OL). Subsequently, we switch the sign of the interaction to induce interatomic attraction and monitor the ensuing dynamics. Regions of the stable self-trapping and catastrophic collapse of 2D fundamental solitons are identified in the parameter plane of the OL strength and BEC norm. The increase of the OL strength expands the persistence domain for the solitons to larger norms. For single-charged solitary vortices, in addition to the survival and collapse regimes, an intermediate one is identified, where the vortex resists the collapse but loses its structure, transforming into a fundamental soliton. The same setting may also be implemented in the context of optical solitons and vortices, using photonic-crystal fibers.

preprint2012arXiv

Solitons supported by singular spatial modulation of the Kerr nonlinearity

We introduce a setting based on the one-dimensional (1D) nonlinear Schroedinger equation (NLSE) with the self-focusing (SF) cubic term modulated by a singular function of the coordinate, |x|^{-a}. It may be additionally combined with the uniform self-defocusing (SDF) nonlinear background, and with a similar singular repulsive linear potential. The setting, which can be implemented in optics and BEC, aims to extend the general analysis of the existence and stability of solitons in NLSEs. Results for fundamental solitons are obtained analytically and verified numerically. The solitons feature a quasi-cuspon shape, with the second derivative diverging at the center, and are stable in the entire existence range, which is 0 < a < 1. Dipole (odd) solitons are found too. They are unstable in the infinite domain, but stable in the semi-infinite one. In the presence of the SDF background, there are two subfamilies of fundamental solitons, one stable and one unstable, which exist together above a threshold value of the norm (total power of the soliton). The system which additionally includes the singular repulsive linear potential emulates solitons in a uniform space of the fractional dimension, 0 < D < 1. A two-dimensional extension of the system, based on the quadratic nonlinearity, is formulated too.

preprint2012arXiv

Solitons supported by spatially inhomogeneous nonlinear losses

We uncover that, in contrast to the common belief, stable dissipative solitons exist in media with uniform gain in the presence of nonuniform cubic losses, whose local strength grows with coordinate x (in one dimension) faster than |x|. The spatially-inhomogeneous absorption also supports new types of solitons, that do not exist in uniform dissipative media. In particular, single-well absorption profiles give rise to spontaneous symmetry breaking of fundamental solitons in the presence of uniform focusing nonlinearity, while stable dipoles are supported by double-well absorption landscapes. Dipole solitons also feature symmetry breaking, but under defocusing nonlinearity.

preprint2012arXiv

Spatial solitons under competing linear and nonlinear diffractions

We introduce a general model which augments the one-dimensional nonlinear Schrödinger (NLS) equation by nonlinear-diffraction terms competing with the linear diffraction. The new terms contain two irreducible parameters and admit a Hamiltonian representation in a form natural for optical media. The equation serves as a model for spatial solitons near the supercollimation point in nonlinear photonic crystals. In the framework of this model, a detailed analysis of the fundamental solitary waves is reported, including the variational approximation (VA), exact analytical results, and systematic numerical computations. The Vakhitov-Kolokolov (VK) criterion is used to precisely predict the stability border for the solitons, which is found in an exact analytical form, along with the largest total power (norm) that the waves may possess. Past a critical point, collapse effects are observed, caused by suitable perturbations. Interactions between two identical parallel solitary beams are explored by dint of direct numerical simulations. It is found that in-phase solitons merge into robust or collapsing pulsons, depending on the strength of the nonlinear diffraction.

preprint2012arXiv

Spontaneous symmetry breaking of binary fields in a nonlinear double-well structure

We introduce a one-dimensional two-component system with the self-focusing cubic nonlinearity concentrated at a symmetric set of two spots. Effects of the spontaneous symmetry breaking (SSB) of localized modes were previously studied in the single-component version of this system. In this work, we study the evolution (in the configuration space of the system) and SSB scenarios for two-component modes of three generic types, as concerns the spatial symmetry of each component: symmetric-symmetric (Sm-Sm), antisymmetric-antisymmetric (AS-AS), and symmetric-antisymmetric (S-AS) ones. In the limit case of the nonlinear potential represented by two $% δ$-functions, solutions are obtained in a semi-analytical form. They feature novel properties, in comparison with the previously studied single-component model. In particular, the SSB of antisymmetric modes is possible solely in the two-component system, and, obviously, S-AS states exist only in the two-component system too. In the general case of the symmetric pair of finite-width nonlinear potential wells, evolution scenarios are very complex. In this case, new results are reported, first, for the single-component model. These are pairs of broken-antisymmetry modes, and of twin-peak symmetric ones, which are generated by saddle-mode bifurcations separated from the transformations previously studied in the the single-component setting. With regard to these findings, complex scenarios of the evolution of the two-component solution families are realized in terms of links connecting pairs of modes of three simplest types: (A) two-component ones with unbroken symmetries; (B) single-component modes featuring density peaks in both potential wells; (C) single-component modes which are trapped, essentially, in a single well.

preprint2012arXiv

Stabilization of one-dimensional solitons against the critical collapse by quintic nonlinear lattices

It has been recently discovered that stabilization of two-dimensional (2D) solitons against the critical collapse in media with the cubic nonlinearity by means of nonlinear lattices (NLs) is a challenging problem. We address the 1D version of the problem, i.e., the nonlinear-Schrödinger equation (NLSE) with the quintic or cubic-quintic (CQ) terms, the coefficient in front of which is periodically modulated in space. The models may be realized in optics and Bose-Einstein condensates (BECs). Stability diagrams for the solitons are produced by means of numerical methods and analytical approximations. It is found that the sinusoidal NL stabilzes solitons supported by the quintic-only nonlinearity in a narrow stripe in the respective parameter plane, on the contrary to the case of the cubic nonlinearity in 2D, where the stabilization of solitons by smooth spatial modulations is not possible at all. The stability region is much broader in the 1D CQ model, where higher-order solitons may be stable too.

preprint2012arXiv

Stabilizing single- and two-color vortex beams in quadratic media by a trapping potential

We consider two-dimensional (2D) localized modes in the second-harmonic-generating χ^{(2)} system with the harmonic-oscillator (HO) trapping potential. In addition to its realization in optics, the system describes the mean-field dynamics of mixed atomic-molecular Bose-Einstein condensates (BECs). The existence and stability of various modes is determined by their total power, N, topological charge, m/2 [m is the intrinsic vorticity of the second-harmonic (SH) field], and χ^{(2)} mismatch, q. The analysis is carried out in a numerical form and, in parallel, by means of the variational approximation (VA), which produces results that agree well with numerical findings. Below a certain power threshold, N\leq N_{c}^{(m)}(q), all trapped modes are of the single-color type, represented by the SH component only, while the fundamental-frequency (FF) one is absent. In contrast with the usual situation, where such modes are always unstable, we demonstrate\ that they are stable, for m=0,1,2 (the mode with m=1 may be formally considered as a semi-vortex with topological charge m/2=1/2), at N\leq N_{c}^{(m)}(q), and unstable above this threshold. On the other hand, N_{c}^{(m)}(q)\equiv 0 at q\geq q_{\max} (in our notation, q_{\max}=1$), hence the single-color modes are unstable in the latter case. At N=N_{c}^{(m)}, the modes with m=0 and m=2 undergo a pitchfork bifurcation, which gives rise to two-color states, which remain completely stable for m=0. The two-color vortices with m=2 (topological charge 1) have an upper stability border, N=N_{c2}(q). Above the border, they exhibit periodic splittings and recombinations, while keeping their vorticity. The semi-vortex does not bifurcate; at N=N_{c}^{(m=1)}, it exhibits quasi-chaotic oscillations and a rotating "groove" resembling a screw-edge dislocation induced by the semi-integer vorticity.

preprint2012arXiv

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 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.

preprint2012arXiv

Stable dark solitons in PT-symmetric dual-core waveguides

We construct dark solitons in the recently introduced model of the nonlinear dual-core coupler with the mutually balanced gain and loss applied to the two cores, which is a realization of parity-time symmetry in nonlinear optics. The main issue is stability of the dark solitons. The modulational stability of the CW (continuous-wave) background, which supports the dark solitons, is studied analytically, and the full stability is investigated in a numerical form, via computation of eigenvalues for modes of small perturbations. Stability regions are thus identified in the parameter space of the system, and verified in direct simulations. Collisions between stable dark solitons are briefly considered too.

preprint2012arXiv

Stable nonlinear amplification of solitons without gain saturation

We demonstrate that the cubic gain applied in a localized region, which is embedded into a bulk waveguide with the cubic-quintic nonlinearity and uniform linear losses, supports stable spatial solitons in the absence of the quintic dissipation. The system, featuring the bistability between the solitons and zero state (which are separated by a family of unstable solitons), may be used as a nonlinear amplifier for optical and plasmonic solitons, which, on the contrary to previously known settings, does not require gain saturation. The results are obtained in an analytical form and corroborated by the numerical analysis.

preprint2012arXiv

Stable periodic density waves in dipolar Bose-Einstein condensates trapped in optical lattices

Density-wave patterns in (quasi-) discrete media with local interactions are known to be unstable. We demonstrate that \emph{stable} double- and triple- period patterns (DPPs and TPPs), with respect to the period of the underlying lattice, exist in media with nonlocal nonlinearity. This is shown in detail for dipolar Bose-Einstein condensates (BECs), loaded into a deep one-dimensional (1D) optical lattice (OL), by means of analytical and numerical methods in the tight-binding limit. The patterns featuring multiple periodicities are generated by the modulational instability of the continuous-wave (CW) state, whose period is identical to that of the OL. The DPP and TPP emerge via phase transitions of the second and first kind, respectively. The emerging patterns may be stable provided that the dipole-dipole (DD) interactions are repulsive and sufficiently strong, in comparison with the local repulsive nonlinearity. Within the set of the considered states, the TPPs realize a minimum of the free energy. Accordingly, a vast stability region for the TPPs is found in the parameter space, while the DPP\ stability region is relatively narrow. The same mechanism may create stable density-wave patterns in other physical media featuring nonlocal interactions, such as arrayed optical waveguides with thermal nonlinearity.

preprint2012arXiv

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 (2D) 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 (VA) and Vakhitov-Kolokolov (VK) stability criterion, and corroborated by numerical computations. For the small contrast, the stability interval of the soliton's norm scales as $ΔN\sim Δ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 1D solitons against the critical collapse under the action of a locally enhanced quintic SF nonlinearity.

preprint2012arXiv

Symbiotic two-component gap solitons

We consider a two-component one-dimensional model of gap solitons (GSs), which is based on two nonlinear Schrödinger equations, coupled by repulsive XPM (cross-phase-modulation) terms, in the absence of the SPM (self-phase-modulation) nonlinearity. The equations include a periodic potential acting on both components, thus giving rise to GSs of the "symbiotic" type, which exist solely due to the repulsive interaction between the two components. The model may be implemented for "holographic solitons" in optics, and in binary bosonic or fermionic gases trapped in the optical lattice. Fundamental symbiotic GSs are constructed, and their stability is investigated, in the first two finite bandgaps of the underlying spectrum. Symmetric solitons are destabilized, including their entire family in the second bandgap, by symmetry-breaking perturbations above a critical value of the total power. Asymmetric solitons of intra-gap and inter-gap types are studied too, with the propagation constants of the two components falling into the same or different bandgaps, respectively. The increase of the asymmetry between the components leads to shrinkage of the stability areas of the GSs. Inter-gap GSs are stable only in a strongly asymmetric form, in which the first-bandgap component is a dominating one. Intra-gap solitons are unstable in the second bandgap. Unstable two-component GSs are transformed into persistent breathers. In addition to systematic numerical considerations, analytical results are obtained by means of an extended ("tailed") Thomas-Fermi approximation (TFA).

preprint2012arXiv

Symmetric and asymmetric solitons in a nonlocal nonlinear coupler

We study effects of nonlocality of the cubic self-focusing nonlinearity on the stability and symmetry-breaking bifurcation (SBB) of solitons in the model of a planar dual-core optical waveg- uide with nonlocal (thermal) nonlinearity. In comparison with the well-known coupled systems with the local nonlinearity, the present setting is affected by the competition of different spatial scales, viz., the coupling length and correlation radius of the nonlocality, $\sqrt{d}$. By means of numer- ical methods and variational approximation (VA, which is relevant for small d), we find that, with the increase of the correlation radius, the SBB changes from subcritical into supercritical, which makes all the asymmetric solitons stable. On the other hand, the nonlocality has little influence on the stability of antisymmetric solitons. Analytical results for the SBB are also obtained (actually, for antisymmetric "accessible solitons") in the opposite limit of the ultra-nonlocal nonlinearity, using a coupler based on the Snyder-Mitchell model. The results help to grasp the general picture of the symmetry breaking in nonlocal couplers.

preprint2012arXiv

Symmetry breaking in dipolar matter-wave solitons in dual-core couplers

We study effects of the spontaneous symmetry-breaking (SSB) in solitons built of the dipolar Bose-Einstein condensate (BEC), trapped in a dual-core system with the dipole-dipole interactions (DDIs) and hopping between the cores. Two realizations of such a matter-wave coupler are introduced, weakly- and strongly-coupled. The former one in based on two parallel pipe-shaped traps, while the latter one is represented by a single pipe sliced by an external field into parallel layers. The dipoles are oriented along axes of the pipes. In these systems, the dual-core solitons feature the SSB of the supercritical type and subcritical types, respectively. Stability regions are identified for symmetric and asymmetric solitons, and, in addition, for non-bifurcating antisymmetric ones, as well as for symmetric flat states, which may also be stable in the strongly-coupled system, due to competition between the attractive and repulsive intra- and inter-core DDIs. Effects of the contact interactions are considered too. Collisions between moving asymmetric solitons in the weakly-symmetric system feature elastic rebound, merger into a single breather, and passage accompanied by excitation of intrinsic vibrations of the solitons, for small, intermediate, and large collision velocities, respectively. A PT-symmetric version of the weakly-coupled system is briefly considered too, which may be relevant for matter-wave lasers. Stability boundaries for PT-symmetric and antisymmetric solitons are identified.

preprint2012arXiv

Two-dimensional solitons and vortices in media with incommensurate linear and nonlinear lattice potentials

We construct families of ordinary and gap solitons (GSs), including solitary vortices, in the two-dimensional (2D) system based on the nonlinear-SchrÄodinger/Gross-Pitaevskii equation with the 2D or quasi-1D (Q1D) periodic linear potential, combined with the periodic modulation of the cubic nonlinearity (also in the 2D or Q1D form), which is, generally, incommensurate with the linear potential, thus forming a \nonlinear quasicrystal". Stable vortices are built as complexes of four peaks with the separation between them equal to the double period of the linear potential. The system may be realized in photonic crystals or Bose-Einstein condensates (BECs). The variational approximation (VA) is applied to ordinary solitons (residing in the semi-infinite gap), and numerical methods are used to construct solitons of all the types. Stability regions are identified for soliton families in all the versions of the model.

preprint2012arXiv

Unstaggered-staggered solitons in two-component discrete nonlinear Schrödinger lattices

We present stable bright solitons built of coupled unstaggered and staggered components in a symmetric system of two discrete nonlinear Schrödinger (DNLS) equations with the attractive self-phase-modulation (SPM) nonlinearity, coupled by the repulsive cross-phase-modulation (XPM) interaction. These mixed modes are of a "symbiotic" type, as each component in isolation may only carry ordinary unstaggered solitons. The results are obtained in an analytical form, using the variational and Thomas-Fermi approximations (VA and TFA), and the generalized Vakhitov-Kolokolov (VK) criterion for the evaluation of the stability. The analytical predictions are verified against numerical results. Almost all the symbiotic solitons are predicted by the VA quite accurately, and are stable. Close to a boundary of the existence region of the solitons (which may feature several connected branches), there are broad solitons which are not well approximated by the VA, and are unstable.

preprint2012arXiv

Wave modes trapped in rotating nonlinear potentials

We study modes trapped in a rotating ring with the local strength of the nonlinearity modulated as cos(2θ), where θis the azimuthal angle. This modulation pattern may be of three different types: self-focusing (SF), self-defocusing (SDF), and alternating SF-SDF. The model, based on the nonlinear Schrodinger (NLS) equation with periodic boundary conditions, applies to the light propagation in a twisted pipe waveguide, and to a Bose-Einstein condensate (BEC) loaded into a toroidal trap, under the action of the rotating nonlinear pseudopotential induced by means of the Feshbach resonance in an inhomogeneous external field. In the SF, SDF, and alternating regimes, four, three, and five different types of stable trapped modes are identified, respectively: even, odd, second-harmonic (2H), symmetry-breaking, and 2H-breaking ones. The shapes and stability of these modes, together with transitions between them, are investigated in the first rotational Brillouin zone. Ground-state modes are identified in each regime. Boundaries between symmetric and asymmetric modes are also obtained in an analytical form, by of a two-mode approximation.

preprint2012arXiv

Wave scattering on a domain wall in a chain of PT-symmetric couplers

We study wave propagation in linear arrays composed of pairs of conjugate waveguides with balanced gain and loss, i.e. arrays of the PT-symmetric couplers, where the linear spectrum is known to feature high-frequency and low-frequency branches. We introduce a domain wall by switching the gain and loss in a half of the array, and analyze the scattering of linear waves on this defect. The analysis reveals two major effects: amplification of both reflected and transmitted waves, and excitation of the reflected and transmitted low-frequency and high-frequency waves by the incident high-frequency and low-frequency waves, respectively.

preprint2011arXiv

A nonpolynomial Schroedinger equation for resonantly absorbing gratings

We derive a nonlinear Schroedinger equation with a radical term, in the form of the square root of (1-|V|^2), as an asymptotic model of the optical medium built as a periodic set of thin layers of two-level atoms, resonantly interacting with the electromagnetic field and inducing the Bragg reflection. A family of bright solitons is found, which splits into stable and unstable parts, exactly obeying the Vakhitov-Kolokolov criterion. The soliton with the largest amplitude, which is |V| = 1, is found in an explicit analytical form. It is a "quasi-peakon", with a discontinuity of the third derivative at the center. Families of exact cnoidal waves, built as periodic chains of quasi-peakons, are found too. The ultimate solution belonging to the family of dark solitons, with the background level |V| = 1, is a dark compacton, also obtained in an explicit analytical form. Those bright solitons which are unstable destroy themselves (if perturbed) attaining the critical amplitude, |V| = 1. The dynamics of the wave field around this critical point is studied analytically, revealing a switch of the system into an unstable phase. Collisions between bright solitons are investigated too. The collisions between fast solitons are quasi-elastic, while slowly moving ones merge into breathers, which may persist or perish (in the latter case, also by attaining |V| = 1).

preprint2011arXiv

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+abs(r)]**a 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., a>D. Fundamental solitons are always stable, while multipoles and vortices are stable if the nonlinearity growth rate is large enough.

preprint2011arXiv

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 1D fundamental and multihump states, 2D vortex solitons with arbitrarily high topological charges, and fundamental solitons in 3D. Solitons maintain their coherence in the state of motion, oscillating in the nonlinear potential as robust quasi-particles 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.

preprint2011arXiv

Domain walls and vortices in linearly coupled systems

We investigate 1D and 2D radial domain-wall (DW) states in the system of two nonlinear-Schrödinger/Gross-Pitaevskii equations, which are coupled by the linear mixing and by the nonlinear XPM (cross-phase-modulation). The system has straightforward applications to two-component Bose-Einstein condensates, and to the bimodal light propagation in nonlinear optics. In the former case, the two components represent different hyperfine atomic states, while in the latter setting they correspond to orthogonal polarizations of light. Conditions guaranteeing the stability of flat continuous wave (CW) asymmetric bimodal states are established, followed by the study of families of the corresponding DW patterns. Approximate analytical solutions for the DWs are found near the point of the symmetry-breaking bifurcation of the CW states. An exact DW solution is produced for ratio 3:1 of the XPM and SPM coefficients. The DWs between flat asymmetric states, which are mirror images to each other, are completely stable, and all other species of the DWs, with zero crossings in one or two components, are fully unstable. Interactions between two DWs are considered too, and an effective potential accounting for the attraction between them is derived analytically. Direct simulations demonstrate merger and annihilation of the interacting DWs. The analysis is extended for the system including single- and double-peak external potentials. Generic solutions for trapped DWs are obtained in a numerical form, and their stability is investigated. An exact stable solution is found for the DW trapped by a single-peak potential. In the 2D geometry, stable two-component vortices are found, with topological charges s=1,2,3. Radial oscillations of annular DW-shaped pulsons, with s=0,1,2, are studied too. A linear relation between the period of the oscillations and the mean radius of the DW ring is derived analytically.

preprint2011arXiv

Double symmetry breaking of solitons in one-dimensional virtual photonic crystals

We demonstrate that spatial solitons undergo two consecutive spontaneous symmetry breakings (SSBs), with the increase of the total power, in nonlinear photonic crystals (PhCs) built as arrays of alternating linear and nonlinear stripes, in the case when maxima of the effective refractive index coincide with minima of the self-focusing coefficient, and vice versa, i.e.,the corresponding linear and nonlinear periodic potentials are in competition. This setting may be induced, as a virtual PhC, by means of the EIT (electromagnetically-induced-transparency) technique, in a uniform optical medium. It may also be realized as a Bose-Einstein condensate (BEC) subject to the action of combined periodic optical potential and periodically modulated Feshbach resonance. The first SSB happens at the center of a linear stripe, pushing a broad low-power soliton into an adjacent nonlinear stripe and gradually suppressing side peaks in the soliton's shape. Then, the soliton restores its symmetry, being pinned to the midpoint of the nonlinear stripe. The second SSB occurs at higher powers, pushing the narrow soliton off the center of the nonlinear channel,while the soliton keeps its internal symmetry. The results are obtained by means of numerical and analytical methods. They may be employed to control switching of light beams by means of the varying power.

preprint2011arXiv

Dynamics of bright solitons and soliton arrays in the nonlinear Schrodinger equation with a combination of random and harmonic potentials

We report results of systematic simulations of the dynamics of solitons in the framework of the one-dimensional nonlinear Schrödinger equation (NLSE), which includes the harmonic-oscillator (HO) potential and a random potential. The equation models experimentally relevant spatially disordered settings in Bose-Einstein condensates (BECs) and nonlinear optics. First, the generation of soliton arrays from a broad initial quasi-uniform state by the modulational instability (MI) is considered, following a sudden switch of the nonlinearity from repulsive to attractive. Then, we study oscillations of a single soliton in this setting, which models a recently conducted experiment in BEC. Basic characteristics of the MI-generated array, such as the number of solitons and their mobility, are reported as functions of the strength and correlation length of the disorder, and of the total norm. For the single oscillating soliton, its survival rate is found. Main features of these dependences are explained qualitatively.

preprint2011arXiv

Gap solitons under competing local and nonlocal nonlinearities

We analyze the existence, bifurcations, and shape transformations of one-dimensional gap solitons (GSs) in the first finite bandgap induced by a periodic potential built into materials with local self-focusing and nonlocal self-defocusing nonlinearities. Originally stable on-site GS modes become unstable near the upper edge of the bandgap with the introduction of the nonlocal self-defocusing nonlinearity with a small nonlocality radius. Unstable off-site GSs bifurcate into a new branch featuring single-humped, double-humped, and flat-top modes due to the competition between local and nonlocal nonlinearities. The mechanism underlying the complex bifurcation pattern and cutoff effects (termination of some bifurcation branches) is illustrated in terms of the shape transformation under the action of the varying degree of the nonlocality. The results of this work suggest a possibility of optical-signal processing by means of the competing nonlocal and local nonlinearities.

preprint2011arXiv

Localization and delocalization of two-dimensional discrete solitons pinned to linear and nonlinear defects

We study the dynamics of two-dimensional (2D) localized modes in the nonlinear lattice described by the discrete nonlinear Schrödinger (DNLS) equation, including a local linear or nonlinear defect. Discrete solitons pinned to the defects are investigated by means of the numerical continuation from the anti-continuum limit and also using the variational approximation (VA), which features a good agreement for strongly localized modes. The models with the time-modulated strengths of the linear or nonlinear defect are considered too. In that case, one can temporarily shift the critical norm, below which localized 2D modes cannot exists, to a level above the norm of the given soliton, which triggers the irreversible delocalization transition.

preprint2011arXiv

Multi-stable dissipative structures pinned to dual hot spots

We analyze the formation of one-dimensional localized patterns in a nonlinear dissipative medium including a set of two narrow "hot spots" (HSs), which carry the linear gain, local potential, cubic self-interaction, and cubic loss, while the linear loss acts in the host medium. This system can be realized, as a spatial-domain one, in optics, and also in Bose-Einstein condensates of quasi-particles in solid-state settings. Recently, exact solutions were found for localized modes pinned to the single HS represented by the delta-function. The present paper reports analytical and numerical solutions for coexisting two- and multi-peak modes, which may be symmetric or antisymmetric with respect to the underlying HS pair. Stability of the modes is explored through simulations of their perturbed evolution. The sign of the cubic nonlinearity plays a crucial role: in the case of the self-focusing, only the fundamental symmetric and antisymmetric modes, with two local peaks tacked to the HSs, and no additional peaks between them, may be stable. In this case, all the higher-order multi-peak modes, being unstable, evolve into the fundamental ones. Stability regions for the fundamental modes are reported. A more interesting situation is found in the case of the self-defocusing cubic nonlinearity, with the HS pair giving rise to a multi-stability, with up to eight coexisting stable multi-peak patterns, symmetric and antisymmetric ones. The system without the self-interaction, the nonlinearity 2 being represented only by the local cubic loss, is investigated too. This case is similar to those with the self-focusing or defocusing nonlinearity, if the linear potential of the HS is, respectively, attractive or repulsive. An additional feature of the former setting is the coexistence of the stable fundamental modes with robust breathers.

preprint2011arXiv

Nonlinearly-PT-symmetric systems: spontaneous symmetry breaking and transmission resonances

We introduce a class of PT-symmetric systems which include mutually matched nonlinear loss and gain (inother words, a class of PT-invariant Hamiltonians in which both the harmonic and anharmonic parts are non-Hermitian). For a basic system in the form of a dimer, symmetric and asymmetric eigenstates, including multistable ones, are found analytically. We demonstrate that, if coupled to a linear chain, such a nonlinear PT-symmetric dimer generates new types of nonlinear resonances, with the completely suppressed or greatly amplified transmission, as well as a regime similar to the electromagnetically-induced transparency (EIT). The implementation of the systems is possible in various media admitting controllable linear and nonlinear amplification of waves.

preprint2011arXiv

Quasi-compactons in inverted nonlinear photonic crystals

We study large-amplitude one-dimensional solitary waves in photonic crystals featuring competition between linear and nonlinear lattices, with minima of the linear potential coinciding with maxima of the nonlinear pseudopotential, and vice versa (inverted nonlinear photonic crystals, INPhCs), in the case of the saturable self-focusing nonlinearity. Such crystals were recently fabricated using a mixture of SU-8 and Rhodamine-B optical materials. By means of numerical methods and analytical approximations, we find that large-amplitude solitons are broad sharply localized stable pulses (quasi-compactons, QCs). With the increase of the totalpower, P, the QC's centroid performs multiple switchings between minima and maxima of the linear potential. Unlike cubic INPhCs, the large-amplitude solitons are mobile in the medium with the saturable nonlinearity. The threshold value of the kick necessary to set the soliton in motion is found as a function of P. Collisions between moving QCs are considered too.

preprint2011arXiv

Self-trapping and splitting of bright vector solitons under inhomogeneous defocusing nonlinearities

We show that bimodal systems with a spatially nonuniform defocusing cubic nonlinearity, whose strength grows toward the periphery, can support stable two-component solitons. For a sufficiently strong XPM interaction, vector solitons with overlapping components become unstable, while stable families of solitons with spatially separated components emerge. Stable complexes with separated components may be composed not only of fundamental solitons, but of multipoles too.

preprint2011arXiv

Solitons in a chain of PT-invariant dimers

Dynamics of a chain of interacting parity-time invariant nonlinear dimers is investigated. A dimer is built as a pair of coupled elements with equal gain and loss. A relation between stationary soliton solutions of the model and solitons of the discrete nonlinear Schrodinger (DNLS) equation is demonstrated. Approximate solutions for solitons whose width is large in comparison to the lattice spacing are derived, using a continuum counterpart of the discrete equations. These solitons are mobile, featuring nearly elastic collisions. Stationary solutions for narrow solitons, which are immobile due to the pinning by the effective Peierls-Nabarro potential, are constructed numerically, starting from the anti-continuum limit. The solitons with the amplitude exceeding a certain critical value suffer an instability leading to blowup, which is a specific feature of the nonlinear PT-symmetric chain, making it dynamically different from DNLS lattices. A qualitative explanation of this feature is proposed. The instability threshold drops with the increase of the gain-loss coefficient, but it does not depend on the lattice coupling constant, nor on the soliton's velocity.

preprint2011arXiv

Solitons supported by localized nonlinearities in periodic media

Nonlinear periodic systems, such as photonic crystals and Bose-Einstein condensates (BECs) loaded into optical lattices, are often described by the nonlinear Schrödinger/Gross-Pitaevskii equation with a sinusoidal potential. Here, we consider a model based on such a periodic potential, with the nonlinearity (attractive or repulsive) concentrated either at a single point or at a symmetric set of two points, which are represented, respectively, by a single δ-function or a combination of two δ-functions. This model gives rise to ordinary solitons or gap solitons (GSs), which reside, respectively, in the semi-infinite or finite gaps of the system's linear spectrum, being pinned to the δ-functions. Physical realizations of these systems are possible in optics and BEC, using diverse variants of the nonlinearity management. First, we demonstrate that the single δ-function multiplying the nonlinear term supports families of stable regular solitons in the self-attractive case, while a family of solitons supported by the attractive δ-function in the absence of the periodic potential is completely unstable. We also show that the δ-function can support stable GSs in the first finite gap in both the self-attractive and repulsive models. The stability analysis for the GSs in the second finite gap is reported too, for both signs of the nonlinearity. Alongside the numerical analysis, analytical approximations are developed for the solitons in the semi-infinite and first two finite gaps, with the single δ-function positioned at a minimum or maximum of the periodic potential. In the model with the symmetric set of two δ-functions, we study the effect of the spontaneous symmetry breaking of the pinned solitons. Two configurations are considered, with the δ-functions set symmetrically with respect to the minimum or maximum of the potential.

preprint2011arXiv

Spatial quadratic solitons guided by narrow layers of a nonlinear material

We report analytical solutions for spatial solitons supported by layers of a quadratically nonlinear material embedded into a linear planar waveguide. A full set of symmetric, asymmetric, and antisymmetric modes pinned to a symmetric pair of the nonlinear layers is obtained. The solutions describe a bifurcation of the subcritical type, which accounts for the transition from the symmetric to asymmetric modes. The antisymmetric states (which do not undergo the bifurcation) are completely stable (the stability of the solitons pinned to the embedded layers is tested by means of numerical simulations). Exact solutions are also found for nonlinear layers embedded into a nonlinear waveguide, including the case when the uniform and localized nonlinearities have opposite signs (competing nonlinearities). For the layers embedded into the nonlinear medium, stability properties are explained by comparison to the respective cascading limit.

preprint2011arXiv

Spontaneous symmetry breaking in linearly coupled disk-shaped Bose-Einstein condensates

We study effects of tunnel coupling on a pair of parallel disk-shaped Bose-Einstein condensates with the self-attractive intrinsic nonlinearity. Each condensate is trapped in a combination of in-plane and transverse harmonic-oscillator potentials. It is shown that, depending on the self-interaction strength and tunneling coupling, the ground state of the system exhibits a phase transition which links three configurations: a symmetric one with equal numbers of atoms in the coupled condensates, an asymmetric configuration with a population imbalance (a manifestation of the macroscopic quantum self-trapping), and the collapsing state. A modification of the phase diagram of the system in the presence of vortices in the disk-shaped condensates is reported too. The study of dynamics around the stationary configurations reveals properties which strongly depend on the symmetry of the configuration.

preprint2011arXiv

Spontaneous symmetry breaking in Schrödinger lattices with two nonlinear sites

We introduce discrete systems in the form of straight (infinite) and ring-shaped chains, with two symmetrically placed nonlinear sites. The systems can be implemented in nonlinear optics (as waveguiding arrays) and BEC (by means of an optical lattice). A full set of exact analytical solutions for symmetric, asymmetric, and antisymmetric localized modes is found, and their stability is investigated in a numerical form. The symmetry-breaking bifurcation (SBB), through which the asymmetric modes emerge from the symmetric ones, is found to be of the subcritical type. It is transformed into a supercritical bifurcation if the nonlinearity is localized in relatively broad domains around two central sites, and also in the ring of a small size, i.e., in effectively nonlocal settings. The family of antisymmetric modes does not undergo bifurcations, and features both stable and unstable portions. The evolution of unstable localized modes is investigated by means of direct simulations. In particular, unstable asymmetric states, which exist in the case of the subcritical bifurcation, give rise to breathers oscillating between the nonlinear sites, thus restoring an effective dynamical symmetry between them.

preprint2011arXiv

Spontaneous symmetry breaking of photonic and matter waves in two-dimensional pseudopotentials

We introduce the two-dimensional Gross-Pitaevskii/nonlinear-Schrodinger (GP/NLS) equation with the self-focusing nonlinearity confined to two identical circles, separated or overlapped. The model can be realized in terms of Bose-Einstein condensates (BECs) and photonic-crystal fibers. Following the recent analysis of the spontaneous symmetry breaking (SSB) of localized modes trapped in 1D and 2D double-well nonlinear potentials (also known as pseudopotentials), we aim to find 2D solitons in the two-circle setting, using numerical methods and the variational approximation (VA). Well-separated circles support stable symmetric and antisymmetric solitons. The decrease of separation L between the circles leads to destabilization of the solitons. The symmetric modes undergo two SSB transitions. First, they are transformed into weakly asymmetric breathers, which is followed by a transition to single-peak modes trapped in one circle. The antisymmetric solitons perform a direct transition to the single-peak mode. The symmetric solitons are described reasonably well by the VA. For touching (L=0) and overlapping (L<0) circles, single-peak solitons are found-asymmetric ones, trapped in either circle, and symmetric solitons centered at the midpoint of the bi-circle configuration. If the overlap is weak, the symmetric soliton is unstable. It may spontaneously leap into either circle and perform shuttle motion in it. A region of stability of the symmetric solitons appears with the increase of the overlap degree. In the case of a moderately strong overlap, another SSB effect is found, in the form of a pair of symmetry-breaking and restoring bifurcations which link families of the symmetric and asymmetric solitons.

preprint2011arXiv

Stability of solitons in PT-symmetric couplers

Families of analytical solutions are found for symmetric and antisymmetric solitons in the dual-core system with the Kerr nonlinearity and PT-balanced gain and loss. The crucial issue is stability of the solitons. A stability region is obtained in an analytical form, and verified by simulations, for the PT-symmetric solitons. For the antisymmetric ones, the stability border is found in a numerical form. Moving solitons of both types collide elastically. The two soliton species merge into one in the "supersymmetric" case, with equal coefficients of the gain, loss and inter-core coupling. These solitons feature a subexponential instability, which may be suppressed by periodic switching ("management").

preprint2011arXiv

Symmetric and asymmetric solitons trapped in H-shaped potentials

We report results of numerical and analytical studies of the spontaneous symmetry breaking in solitons, both two- and one-dimensional, which are trapped in H-shaped potential profiles, built of two parallel potential troughs linked by a narrow rung in the transverse direction. This system can be implemented in self-attractive Bose-Einstein condensates (BECs), as well as in a nonlinear bulk optical waveguide.We demonstrate that the introduction of the transverse link changes the character of the symmetry-breaking bifurcation (SBB) in the system from subcritical to supercritical (in terms of the corresponding phase transition, it is a change between the first and second kinds). A noteworthy feature of the SBB in this setting is a non-monotonous dependence of the soliton's norm at the bifurcation point on the strength of the transverse link. In the full 2D system, the results are obtained in a numerical form. An exact analytical solution is found for the bifurcation in the 1D version of the model, with the transverse rung modeled by the local linear coupling between the parallel troughs with the Delta-functional longitudinal profile. Replacing the Delta-function by its finite-width Gaussian counterpart, similar results are obtained by means of the variational approximation (VA). The VA is also applied to the 1D system with a mixed linear and nonlinear transverse localized coupling. Comparison of the results produced by the different varieties of the system clearly reveals basic features of the symmetry-breaking transition in it.

preprint2011arXiv

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 the nonlinear trap for the dragging allows one to pick up and transfer the relevant structures without grabbing surrounding "garbage". A stability border for the dragged modes is identified by means of 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, that may be dragged by the nonlinear trap moving at a constant acceleration, are briefly considered too.

preprint2010arXiv

Arrayed and checkerboard optical waveguides controlled by the electromagnetically-induced transparency

We introduce two models of quasi-discrete optical systems: an array of waveguides doped by four-level N-type atoms, and a nonlinear checkerboard pattern, formed by doping with three-level atoms of the Λ-type. The dopant atoms are driven by external fields, to induce the effect of the electromagnetically-induced transparency (EIT). These active systems offer advantages and addition degrees of freedom, in comparison with ordinary passive waveguiding systems. In the array of active waveguides, the driving field may adjust linear and nonlinear propagation regimes for a probe signal. The nonlinear checkerboard system supports the transmission of stable spatial solitons and their "fuzzy" counterparts, straight or oblique.

preprint2010arXiv

Collapse suppression and stabilization of dipole solitons in two-dimensional media with anisotropic semi-local nonlinearity

We consider the impact of anisotropic nonlocality on the arrest of the collapse and stabilization of dipole-mode (DM) solitons in two-dimensional (2D) models of optical media with the diffusive nonlinearity. The nonlocal nonlinearity is made anisotropic through elliptic diffusivity. The medium becomes semi-local in the limit case of 1D diffusivity. Families of fundamental and DM solitons are found by means of the variational approximation (VA) and in a numerical form. We demonstrate that the collapse of 2D beams is arrested even in the semi-local system. The anisotropic nonlocality readily stabilizes the DM solitons, which are completely unstable in the isotropic medium.

preprint2010arXiv

Exact exact solutions of the Gross-Pitaevskii equation for stable vortex modes

We construct exact solutions of the Gross-Pitaevskii equation for solitary vortices, and approximate ones for fundamental solitons, in 2D models of Bose-Einstein condensates with a spatially modulated nonlinearity of either sign and a harmonic trapping potential. The number of vortex-soliton (VS) modes is determined by the discrete energy spectrum of a related linear Schrödinger equation. The VS families in the system with the attractive and repulsive nonlinearity are mutually complementary. \emph{% Stable} VSs with vorticity $S\geq 2$ and those corresponding to higher-order radial states are reported for the first time, in the case of the attraction and repulsion, respectively.

preprint2010arXiv

Exact solutions for periodic and solitary matter waves in nonlinear lattices

We produce three vast classes of exact periodic and soliton solutions to the one-dimensional Gross-Pitaevskii equation (GPE) with the pseudopotential in the form of a nonlinear lattice (NL), induced by a spatially periodic modulation of the local nonlinearity. It is well known that NLs in Bose-Einstein condensates (BECs) may be created by means of the Feshbach-resonance technique. The model may also include linear potentials with the same periodicity. The NL modulation function, the linear potential (if any), and the corresponding exact solutions are expressed in terms of the Jacobi's elliptic functions of three types, cn, dn, and sn, which give rise to the three different classes of the solutions. The potentials and associated solutions are parameterized by two free constants and an additional sign parameter in the absence of the linear potential. In the presence of the latter, the solution families feature two additional free parameters. The families include both sign-constant and sign-changing NLs. Density maxima of the solutions may coincide with either minima or maxima of the periodic pseudopotential. The solutions reduce to solitons in the limit of the infinite period. The stability of the solutions is tested via direct simulations of the GPE. As a result, stability regions are identified for the periodic solutions and solitons. The periodic patterns of the cn type, and the respective limit-form solutions in the form of bright solitons, may be stable both in the absence and presence of the linear potential. The stability of the two other solution classes, of the dn and sn types, is only possible with the linear potential.

preprint2010arXiv

Gap solitons in optical lattices embedded into nonlocal media

We analyze the existence, stability, and mobility of gap solitons (GSs) in a periodic photonic structure built into a nonlocal self-defocusing medium. Counter-intuitively, the GSs are supported even by a highly nonlocal nonlinearity, which makes the system quasi-linear. Unlike local models, the variational approximation (VA) predicts the GSs in a good agreement with numerical findings, due to the suppression of undulating tails of the solitons.

preprint2010arXiv

Matter-wave solitons and finite-amplitude Bloch waves in optical lattices with a spatially modulated nonlinearity

We investigate solitons and nonlinear Bloch waves in Bose-Einstein condensates trapped in optical lattices. By introducing specially designed localized profiles of the spatial modulation of the attractive nonlinearity, we construct an infinite number of exact soliton solutions in terms of the Mathieu and elliptic functions, with the chemical potential belonging to the semi-infinite bandgap of the optical-lattice-induced spectrum. Starting from the exact solutions, we employ the relaxation method to construct generic families of soliton solutions in a numerical form. The stability of the solitons is investigated through the computation of the eigenvalues for small perturbations, and also by direct simulations. Finally, we demonstrate a virtually exact (in the numerical sense) composition relation between nonlinear Bloch waves and solitons.

preprint2010arXiv

Solitons in combined linear and nonlinear lattice potentials

We study ordinary solitons and gap solitons (GSs) in the effectively one-dimensional Gross-Pitaevskii equation, with a combination of linear and nonlinear lattice potentials. The main points of the analysis are effects of the (in)commensurability between the lattices, the development of analytical methods, viz., the variational approximation (VA) for narrow ordinary solitons, and various forms of the averaging method for broad solitons of both types, and also the study of mobility of the solitons. Under the direct commensurability (equal periods of the lattices, the family of ordinary solitons is similar to its counterpart in the free space. The situation is different in the case of the subharmonic commensurability, with L_{lin}=(1/2)L_{nonlin}, or incommensurability. In those cases, there is an existence threshold for the solitons, and the scaling relation between their amplitude and width is different from that in the free space. GS families demonstrate a bistability, unless the direct commensurability takes place. Specific scaling relations are found for them too. Ordinary solitons can be readily set in motion by kicking. GSs are mobile too, featuring inelastic collisions. The analytical approximations are shown to be quite accurate, predicting correct scaling relations for the soliton families in different cases. The stability of the ordinary solitons is fully determined by the VK (Vakhitov-Kolokolov) criterion, while the stability of GS families follows an inverted ("anti-VK") criterion, which is explained by means of the averaging approximation.

preprint2010arXiv

Solitons in nonlinear lattices

This article offers a comprehensive survey of results obtained for solitons and complex nonlinear wave patterns supported by purely nonlinear lattices (NLs), which represent a spatially periodic modulation of the local strength and sign of the nonlinearity, and their combinations with linear lattices. A majority of the results obtained, thus far, in this field and reviewed in this article are theoretical. Nevertheless, relevant experimental settings are surveyed too, with emphasis on perspectives for implementation of the theoretical predictions in the experiment. Physical systems discussed in the review belong to the realms of nonlinear optics (including artificial optical media, such as photonic crystals, and plasmonics) and Bose-Einstein condensation (BEC). The solitons are considered in one, two, and three dimensions (1D, 2D, and 3D). Basic properties of the solitons presented in the review are their existence, stability, and mobility. Although the field is still far from completion, general conclusions can be drawn. In particular, a novel fundamental property of 1D solitons, which does not occur in the absence of NLs, is a finite threshold value of the soliton norm, necessary for their existence. In multidimensional settings, the stability of solitons supported by the spatial modulation of the nonlinearity is a truly challenging problem, for the theoretical and experimental studies alike. In both the 1D and 2D cases, the mechanism which creates solitons in NLs is principally different from its counterpart in linear lattices, as the solitons are created directly, rather than bifurcating from Bloch modes of linear lattices.

preprint2010arXiv

Suppression of the quantum-mechanical collapse by repulsive interactions in a quantum gas

The quantum-mechanical collapse (alias fall onto the center of particles attracted by potential -1/r^2), or "quantum anomaly", is a well-known issue in the quantum theory. We demonstrate that the mean-field repulsive nonlinearity prevents the collapse and thus puts forward a solution to the quantum-anomaly problem different from that previously developed in the framework of the linear quantum-field theory. This solution may be realized in the 3D or 2D gas of dipolar bosons attracted by a central charge, and in the 2D gas of magnetic dipoles attracted by a current filament. In the 3D setting, the dipole-dipole interactions are also taken into regard, in the mean-field approximation. In lieu of the collapse, the cubic nonlinearity creates a 3D ground state (GS), which does not exist in the respective linear Schroedinger equation (SE). The addition of the harmonic trap gives rise to a tristability, in the case when the SE still does not lead to the collapse. In the 2D setting, the cubic nonlinearity is not strong enough to prevent the collapse; however, the quintic term does it, creating the GS, as well as its counterparts carrying the angular momentum (vorticity). Counter-intuitively, such self-trapped 2D modes exist even in the case of a weakly repulsive potential 1/r^2. In the presence of the harmonic trap, the 2D quintic model with a weakly repulsive central potential 1/r^2 gives rise to three confined modes, the middle one being unstable, spontaneously developing into a breather. In both the 3D and 2D cases, the GS wave functions are found in a numerical form, and also in the form of an analytical approximation, which is asymptotically exact in the limit of the large norm.

preprint2010arXiv

Symmetric and asymmetric solitons and vortices in linearly coupled two-dimensional waveguides with the cubic-quintic nonlinearity

It is well known that the two-dimensional (2D) nonlinear Schrödinger equation (NLSE) with the cubic-quintic (CQ) nonlinearity supports a family of stable fundamental solitons, as well as solitary vortices (alias vortex rings), which are stable for sufficiently large values of the norm. We study stationary localized modes in a symmetric linearly coupled system of two such equations, focusing on asymmetric states. The model may describe "optical bullets" in dual-core nonlinear optical waveguides (including spatiotemporal vortices that were not discussed before), or a Bose-Einstein condensate (BEC) loaded into a "dual-pancake" trap. Each family of solutions in the single-component model has two different counterparts in the coupled system, one symmetric and one asymmetric. Similarly to the earlier studied coupled 1D system with the CQ nonlinearity, the present model features bifurcation loops, for fundamental and vortex solitons alike: with the increase of the total energy (norm), the symmetric solitons become unstable at a point of the direct bifurcation, which is followed, at larger values of the energy, by the reverse bifurcation restabilizing the symmetric solitons. However, on the contrary to the 1D system, the system may demonstrate a double bistability for the fundamental solitons. The stability of the solitons is investigated via the computation of instability growth rates for small perturbations. Vortex rings, which we study for two values of the "spin", s = 1 and 2, may be subject to the azimuthal instability, like in the single-component model. We also develop a quasi-analytical approach to the description of the bifurcations diagrams, based on the variational approximation. Splitting of asymmetric vortices, induced by the azimuthal instability, is studied by means of direct simulations. Interactions between initially quiescent solitons of different types are studied too.

preprint2010arXiv

Symmetry breaking, coupling management, and localized modes in dual-core discrete nonlinear-Schrödinger lattices

We introduce a system of two linearly coupled discrete nonlinear Schrödinger equations (DNLSEs), with the coupling constant subject to a rapid temporal modulation. The model can be realized in bimodal Bose-Einstein condensates (BEC). Using an averaging procedure based on the multiscale method, we derive a system of averaged (autonomous) equations, which take the form of coupled DNLSEs with additional nonlinear coupling terms of the four-wave-mixing type. We identify stability regions for fundamental onsite discrete symmetric solitons (single-site modes with equal norms in both components), as well as for two-site in-phase and twisted modes, the in-phase ones being completely unstable. The symmetry-breaking bifurcation, which destabilizes the fundamental symmetric solitons and gives rise to their asymmetric counterparts, is investigated too. It is demonstrated that the averaged equations provide a good approximation in all the cases. In particular, the symmetry-breaking bifurcation, which is of the pitchfork type in the framework of the averaged equations, corresponds to a Hopf bifurcation in terms of the original system.

preprint2010arXiv

The ground state of binary systems with a periodic modulation of the linear coupling

We consider a quasi-one-dimensional two-component systm, described by a pair of Nonlinear Schrödinger/Gross-Pitaevskii Equations (NLSEs/GPEs), which are coupled by the linear mixing, with local strength $Ω$, and by the nonlinear incoherent interaction. We assume the self-repulsive nonlinearity in both components, and include effects of a harmonic trapping potential. The model may be realized in terms of periodically modulated slab waveguides in nonlinear optics, and in Bose-Einstein condensates too. Depending on the strengths of the linear and nonlinear couplings between the components, the ground states (GSs) in such binary systems may be symmetric or asymmetric. In this work, we introduce a periodic spatial modulation of the linear coupling, making $Ω$ an odd, or even function of the coordinate. The sign flips of $Ω(x)$ strongly modify the structure of the GS in the binary system, as the relative sign of its components tends to lock to the local sign of $Ω$. Using a systematic numerical analysis, and an analytical approximation, we demonstrate that the GS of the trapped system contains one or several kinks (dark solitons) in one component, while the other component does not change its sign. Final results are presented in the form of maps showing the number of kinks in the GS as a function of the system's parameters, with the odd/even modulation function giving rise to the odd/even number of the kinks. The modulation of $Ω(x)$ also produces a strong effect on the transition between states with nearly equal and strongly unequal amplitudes of the two components.

preprint2010arXiv

The inverse problem for the Gross - Pitaevskii equation

Two different methods are proposed for the generation of wide classes of exact solutions to the stationary Gross - Pitaevskii equation (GPE). The first method, suggested by the work by Kondrat'ev and Miller (1966), applies to one-dimensional (1D) GPE. It is based on the similarity between the GPE and the integrable Gardner equation, all solutions of the latter equation (both stationary and nonstationary ones) generating exact solutions to the GPE, with the potential function proportional to the corresponding solutions. The second method is based on the "inverse problem" for the GPE, i.e. construction of a potential function which provides a desirable solution to the equation. Systematic results are presented for 1D and 2D cases. Both methods are illustrated by a variety of localized solutions, including solitary vortices, for both attractive and repulsive nonlinearity in the GPE. The stability of the 1D solutions is tested by direct simulations of the time-dependent GPE.

preprint2010arXiv

Two - dimensional solitons in media with the 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. In recent works, it was concluded 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. The double stripe with the rectangular profile admit stable solitons of three distinct types: symmetric and asymmetric ones with a single peak, and double-peak symmetric solitons. The shape and stability of single-peak solitons of either type are accurately predicted by the VA. Collisions between stable solitons are briefly considered too, by means of direct simulations. Depending on the relative velocity we observe excitation, decay or catastrophic self focusing.

preprint2010arXiv

Two-dimensional gap solitons in elliptic-lattice potentials

We study two-dimensional (2D) matter-wave gap solitons trapped in an elliptically deformed concentric lattice potential, within the framework of the Gross-Pitaevskii equation (GPE) with self-attraction or self-repulsion. For a fixed eccentricity of the lattice, soliton families are found in both the repulsive and attractive models. In the former case, the analysis reveals two kinds of gap solitons trapped in the first oval trough (the ring-shaped potential minimum closest to the center): elliptic annular solitons (EASs), and double solitons (DSs), which are formed by two tightly localized density peaks located at diametrically opposite points of the trough, with zero phase difference between them. With the decrease of the norm, the density distribution in the EAS along the azimuthal direction changes from nearly-uniform to double-peaked and, eventually, to the DS. In the attractive model, there exist only DSs in the oval trough, while EASs are not found. All such solitons without the angular momentum (l = 0) are fully stable. For l is not equal to 0, vortical solitons - both EASs with a sufficiently large norm (in the repulsive model) and DSs (in models with both signs of the nonlinearity) - are quasi-stable, exhibiting rocking motion in the elliptic trough (we consider the cases of l=1 and l=2). At smaller values of the norm, the vortical annular solitons (in the repulsive model) are unstable. Stable fundamental solitons trapped in the central potential well are investigated too, in both the attractive and repulsive models, by means of the variational approximation and numerical methods.

preprint2010arXiv

Two-dimensional vector solitons stabilized by a linear or nonlinear lattice acting in one component

The subject of the work is the stabilization of two-dimensional (2D) two-component (vector) solitons, in media with the attractive cubic nonlinearity, against the collapse by a linear lattice (LL, which is induced by a periodic modulation of the refractive index in optics, or created as an optical lattice in BEC), or by a nonlinear lattice (NL, induced by a periodic modulation of the nonlinearity coefficient). We demonstrate that, due to the XPM (cross-phase-modulation) coupling between the components, the LL or NL acting on a single component is sufficient for the stabilization of vector solitons, that include a component for which the self-focusing medium is uniform. In the case of the LL, the vector solitons are stable almost in their entire existence domain, while the NL can only stabilize the solitons in which the component affected by the lattice carries the norm which is comparable to, or larger than the norm of the component in the uniform medium.

preprint2010arXiv

Vortices and ring dark solitons in nonlinear amplifying waveguides

We consider the generation and propagation of (2+1)-dimensional beams in a nonlinear waveguide with the linear gain. Simple self-similar evolution of the beams is achieved at the asymptotic stage, if the input beams represent the fundamental mode. On the contrary, if they carry vorticity or amplitude nodes (or phase slips), vortex tori and ring dark solitons (RDSs) are generated, featuring another type of the self-similar evolution, with an exponentially shrinking vortex core or notch of the RDS. Numerical and analytical considerations reveal that these self-similar structures are robust entities in amplifying waveguides, being \emph{stable} against azimuthal perturbations.

preprint2009arXiv

Application of the Feshbach-resonance management to a tightly confined Bose-Einstein condensate

We study suppression of the collapse and stabilization of matter-wave solitons by means of time-periodic modulation of the effective nonlinearity, using the nonpolynomial Schroedinger equation (NPSE) for BEC trapped in a tight cigar-shaped potential. By means of systematic simulations, a stability region is identified in the plane of the modulation amplitude and frequency. In the low-frequency regime, solitons feature chaotic evolution, although they remain robust objects.

preprint2009arXiv

Interactions between two-dimensional solitons in the diffractive-diffusive Ginzburg-Landau equation with the cubic-quintic nonlinearity

We report results of systematic numerical analysis of collisions between two and three stable dissipative solitons in the two-dimensional (2D) complex Ginzburg- Landau equation (CGLE) with the cubic-quintic (CQ) combination of gain and loss terms. The equation may be realized as a model of a laser cavity which includes the spatial diffraction, together with the anomalous group-velocity dispersion (GVD) and spectral filtering acting in the temporal direction. Collisions between solitons are possible due to the Galilean invariance along the spatial axis. Outcomes of the collisions are identified by varying the GVD coefficient, D, and the collision "velocity" (actually, it is the spatial slope of the soliton's trajectory). At small velocities, two or three in-phase solitons merge into a single standing one. At larger velocities, both in-phase soliton pairs and pairs of solitons with opposite signs suffer a transition into a delocalized chaotic state. At still larger velocities, all collisions become quasi-elastic. A new outcome is revealed by collisions between slow solitons with opposite signs: they self-trap into persistent wobbling dipoles, which are found in two modifications (horizontal at smaller D, and vertical if D is larger (the horizontal ones resemble "zigzag" bound states of two solitons known in the 1D CGL equation of the CQ type). Collisions between solitons with a finite mismatch between their trajectories are studied too.

preprint2009arXiv

Single- and multi-peak solitons in two-component models of metamaterials and photonic crystals

We report results of the study of solitons in a system of two nonlinear-Schrodinger (NLS) equations coupled by the XPM interaction, which models the co-propagation of two waves in metamaterials(MMs). The same model applies to photonic crystals (PCs), as well as to ordinary optical fibers, close to the zero-dispersion point. A peculiarity of the system is a small positive or negative value of the relative group-velocity dispersion (GVD) coefficient in one equation, assuming that the dispersion is anomalous in the other. In contrast to earlier studied systems of nonlinearly coupled NLS equations with equal GVD coefficients, which generate only simple single-peak solitons, the present model gives rise to families of solitons with complex shapes, which feature extended oscillatory tails and/or a double-peak structure at the center. Regions of existence are identified for single- and double-peak bimodal solitons, demonstrating a broad bistability in the system. Behind the existence border, they degenerate into single-component solutions. Direct simulations demonstrate stability of the solitons in the entire existence regions. Effects of the group-velocity mismatch (GVM) and optical loss are considered too. It is demonstrated that the solitons can be stabilized against the GVM by means of the respective "management" scheme. Under the action of the loss, complex shapes of the solitons degenerate into simple ones, but periodic compensation of the loss supports the complexity.

preprint2009arXiv

Solitons and solitary vortices in "pancake"-shaped Bose-Einstein condensates

We study fundamental and vortical solitons in disk-morphed Bose-Einstein condensates (BECs) subject to strong confinement along the axial direction. Starting from the three-dimensional (3D) Gross-Pitaevskii equation (GPE), we proceed to an effective 2D nonpolynomial Schroeodinger equation (NPSE) derived by means of the integration over the axial coordinate. Results produced by the latter equation are in very good agreement with those obtained from the full 3D GPE, including cases when the formal 2D equation with the cubic nonlinearity is unreliable. The 2D NPSE is used to predict density profiles and dynamical stability of repulsive and attractive BECs with zero and finite topological charge in various planar trapping configurations, including the axisymmetric harmonic confinement and 1D periodic potential. In particular, we find a stable dynamical regime that was not reported before, viz., periodic splitting and recombination of trapped vortices with topological charge 2 or 3 in the self-attractive BEC.

preprint2009arXiv

Solitons in quasi one-dimensional Bose-Einstein condensates with competing dipolar and local interactions

We study families of one-dimensional matter-wave bright solitons supported by the competition of contact and dipole-dipole (DD) interactions of opposite signs. Soliton families are found, and their stability is investigated in the free space, and in the presence of an optical lattice (OL). Free-space solitons may exist with an arbitrarily weak local attraction if the strength of the DD repulsion is fixed. In the case of the DD attraction, solitons do not exist beyond a maximum value of the local-repulsion strength. In the system which includes the OL, a stability region for \textit{subfundamental solitons} (SFSs) is found in the second finite bandgap. For the existence of gap solitons (GSs) under the attractive DD interaction, the contact repulsion must be strong enough. In the opposite case of the DD repulsion, GSs exist if the contact attraction is not too strong. Collisions between solitons in the free space are studied too. In the case of the local attraction, they merge or pass through each other at small and large velocities, respectively. In the presence of the local repulsion, slowly moving solitons bounce from each other.

preprint2009arXiv

Stabilization of multi-beam necklace solitons in circular arrays with spatially modulated nonlinearity

We address necklace solitons supported by circular waveguide arrays with out-of-phase modulation of nonlinearity and linear refractive index. Such two-dimensional necklace solitons appear as rings of multiple out-of-phase bright spots. We show that necklace solitons are stable if their peak intensity is low or moderate and that the domain of stability shrinks with increasing nonlinearity modulation depth. However, we find that the stability domains expand with increasing number of necklace spots.

preprint2009arXiv

Two-dimensional discrete solitons in dipolar Bose-Einstein condensates

We analyze the formation and dynamics of bright unstaggered solitons in the disk-shaped dipolar Bose-Einstein condensate, which features the interplay of contact (collisional) and long-range dipole-dipole (DD) interactions between atoms. The condensate is assumed to be trapped in a strong optical-lattice potential in the disk's plane, hence it may be approximated by a two-dimensional (2D) discrete model, which includes the on-site nonlinearity and cubic long-range (DD) interactions between sites of the lattice. We consider two such models, that differ by the form of the on-site nonlinearity, represented by the usual cubic term, or more accurate nonpolynomial one, derived from the underlying 3D Gross-Pitaevskii equation. Similar results are obtained for both models. The analysis is focused on effects of the DD interaction on fundamental localized modes in the lattice (2D discrete solitons). The repulsive isotropic DD nonlinearity extends the existence and stability regions of the fundamental solitons. New families of on-site, inter-site and hybrid solitons, built on top of a finite background, are found as a result of the interplay of the isotropic repulsive DD interaction and attractive contact nonlinearity. By themselves, these solutions are unstable, but they evolve into robust breathers which exist on an oscillating background. In the presence of the repulsive contact interactions, fundamental localized modes exist if the DD interaction (attractive isotropic or anisotropic) is strong enough. They are stable in narrow regions close to the anticontinuum limit, while unstable solitons evolve into breathers. In the latter case, the presence of the background is immaterial.

preprint2009arXiv

Vector solitons in nonlinear lattices

We consider two-component solitons in a medium with a periodic modulation of the nonlinear coefficient. The modulation enables the existence of complex multihump vector states. In particular, vector solitons composed of dipole and fundamental, or dipole and even double-hump components exist and may be stable. Families of unstable scalar solitons can be stabilized in the vectorial form, due to the coupling to a stable second component.

preprint2006arXiv

Mobility of Discrete Solitons in Quadratically Nonlinear Media

We study the mobility of solitons in second-harmonic-generating lattices. Contrary to what is known for their cubic counterparts, discrete quadratic solitons are mobile not only in the one-dimensional (1D) setting, but also in two dimensions (2D). We identify parametric regions where an initial kick applied to a soliton leads to three possible outcomes, namely, staying put, persistent motion, or destruction. For the 2D lattice, it is found that, for the solitary waves, the direction along which they can sustain the largest kick and can attain the largest speed is the diagonal. Basic dynamical properties of the discrete solitons are also discussed in the context of an analytical approximation, in terms of an effective Peierls-Nabarro potential in the lattice setting.

preprint2004arXiv

Modulational Instabilities and Domain Walls in Coupled Discrete Nonlinear Schrödinger Equations

We consider a system of two discrete nonlinear Schrödinger equations, coupled by nonlinear and linear terms. For various physically relevant cases, we derive a modulational instability criterion for plane-wave solutions. We also find and examine domain-wall solutions in the model with the linear coupling.

preprint2004arXiv

Motion of discrete solitons assisted by nonlinearity management

We demonstrate that periodic modulation of the nonlinearity coefficient in the discrete nonlinear Schrödinger (DNLS) equation can strongly facilitate creation of traveling solitons in the lattice. We predict this possibility in an analytical form, and test it in direct simulations. Systematic simulations reveal several generic dynamical regimes, depending on the amplitude and frequency of the time modulation, and on initial thrust which sets the soliton in motion. These regimes include irregular motion, regular motion of a decaying soliton, and regular motion of a stable one. The motion may occur in both the straight and reverse directions, relative to the initial thrust. In the case of stable motion, extremely long simulations in a lattice with periodic boundary conditions demonstrate that the soliton keeps moving as long as we can monitor without any visible loss. Velocities of moving stable solitons are in good agreement with the analytical prediction, which is based on requiring a resonance between the ac drive and motion of the soliton through the periodic potential. All the generic dynamical regimes are mapped in the model's parameter space. Collisions between moving stable solitons are briefly investigated too, with a conclusion that two different outcomes are possible: elastic bounce, or bounce with mass transfer from one soliton to the other. The model can be realized experimentally in a Bose-Einstein condensate trapped in a deep optical lattice.

preprint2003arXiv

Families of Matter-Waves for Two-Component Bose-Einstein Condensates

We produce several families of solutions for two-component nonlinear Schrödinger/Gross-Pitaevskii equations. These include domain walls and the first example of an antidark or gray soliton in the one component, bound to a bright or dark soliton in the other. Most of these solutions are linearly stable in their entire domain of existence. Some of them are relevant to nonlinear optics, and all to Bose-Einstein condensates (BECs). In the latter context, we demonstrate robustness of the structures in the presence of parabolic and periodic potentials (corresponding, respectively, to the magnetic trap and optical lattices in BECs).

preprint1999arXiv

Embedded Solitons in a Three-Wave System

We report a rich spectrum of isolated solitons residing inside ({\it embedded } into) the continuous radiation spectrum in a simple model of three-wave spatial interaction in a second-harmonic-generating planar optical waveguide equipped with a quasi-one-dimensional Bragg grating. An infinite sequence of fundamental embedded solitons are found, which differ by the number of internal oscillations. Branches of these zero-walkoff spatial solitons give rise, through bifurcations, to several secondary branches of walking solitons. The structure of the bifurcating branches suggests a multistable configuration of spatial optical solitons, which may find straightforward applications for all-optical switching.

Boris A. Malomed

What is connected

Connect this record

See the researcher in context

Building this map preview

224 published item(s)

Collisions and fusion of one- and two-dimensional solitons driven by potential troughs in the cubic-quintic nonlinear Schrödinger equations

Dynamical correlations and chimera-like states of nanoemitters coupled to plasmon-polaritons in a lattice of conducting nanorings

Experimental realisations of the fractional Schrödinger equation in the temporal domain

Fundamental and vortex dissipative quadratic solitons supported by spatially localized gain

General higher-order breathers and rogue waves in the two-component long-wave--short-wave resonance-interaction model

Localized modes in nonlinear fractional systems with deep lattices

Multidimensional dissipative solitons and solitary vortices

Observation of SQUID-like behavior in fiber laser with intra-cavity epsilon-near-zero effect

On-demand harnessing of photonic soliton molecules

One- and two-dimensional solitons in spin-orbit-coupled Bose-Einstein condensates with fractional kinetic energy

Solitons in $\mathcal{PT}$-symmetric systems with spin-orbit coupling and critical nonlinearity

Stability limits for modes held in alternating trapping-expulsive potentials

Stable high-dimensional weak-light soliton molecules and their active control

Trapping wave fields in an expulsive potential by means of linear coupling

Two-dimensional Airy waves and three-wave solitons in quadratic media

Two-dimensional solitons in nonlocal media: a brief review

Vortex-ring quantum droplets in a radially-periodic potential

Dragging spin-orbit-coupled solitons by a moving optical lattice

Embedded solitons in second-harmonic-generating lattices

External light control of three-dimensional ultrashort far-infrared pulses in an inhomogeneous array of carbon nanotubes

Nonlinear dynamics of wave packets in tunnel-coupled harmonic-oscillator traps

Quadratic fractional solitons

Semidiscrete vortex solitons

Collective Excitations of a One-Dimensional Quantum Droplet

Holding and transferring matter-wave solitons against gravity by spin-orbit-coupling tweezers

Josephson oscillations of chirality and identity in two-dimensional solitons in spin-orbit-coupled condensates

Modulational instability, inter-component asymmetry and formation of quantum droplets in one-dimensional binary Bose gases

Nonlinearity and discreteness: Solitons in lattices

Reversible ultrafast soliton switching in dual-core highly nonlinear optical fibers

Shortcuts to adiabaticity for an interacting Bose-Einstein condensate via exact solutions of the generalized Ermakov equation

Singular mean-field states: A brief review of recent results

Singular solitons

Splitting of two-component solitary waves from collisions with narrow potential barriers

Stabilization of one-dimensional Townes solitons by spin-orbit coupling in a dual-core system

Stabilization of single- and multi-peak solitons in the fractional nonlinear Schroedinger equation with a trapping potential

Suppression of the quasi-two-dimensional quantum collapse in the attraction field by the Lee-Huang-Yang effect

Symmetry breaking in a two-component system with repulsive interactions and linear coupling

Symmetry breaking of spatial Kerr solitons in fractional dimension

The Gross-Pitaevskii-Poisson model for an ultracold plasma: density waves and solitons

Tunable nonlinear spectra of anti-directional couplers

Unstaggered-staggered solitons on one- and two-dimensional two-component discrete nonlinear Schrödinger lattices

Vortex solitons in fractional nonlinear Schrödinger equation with the cubic-quintic nonlinearity

Asymptotic dynamics of three-dimensional bipolar ultrashort electromagnetic pulses in an array of semiconductor carbon nanotubes

Coupled density-spin Bose-Einstein condensates dynamics and collapse in systems with quintic nonlinearity

Creating Oscillons and Oscillating Kinks in Two Scalar Field Theories

Gap and embedded solitons in microwave-coupled binary condensates

Mechanisms of Spatiotemporal Mode-Locking

Quasi-one-dimensional approximation for Bose-Einstein condensates transversely trapped by a funnel potential

Soliton dynamics in a fractional complex Ginzburg-Landau model

Spatiotemporal dissipative solitons and vortices in a multi-transverse-mode fiber laser

Analytical Solutions for Beams Passing Apertures with Sharp Boundaries

Commentary (invited update) to article "Spatiotemporal optical solitons", by B. A. Malomed, D. Mihalache, F. Wise, and L. Torner, J. Opt. B: Quantum Semiclass. Opt. 7, R53-R72 (2005)

Matter-wave soliton interferometer based on a nonlinear splitter

Motion of solitons in one-dimensional spin-orbit-coupled Bose-Einstein condensates

Multidimensional solitons: Well-established results and novel findings

One- and two-dimensional solitons in PT-symmetric systems emulating the spin-orbit coupling

Resonant Meta-atoms with Nonlinearities on Demand

Robust self-trapping of vortex beams in a saturable optical medium

Spatiotemporal accessible solitons in fractional dimensions

Spontaneous symmetry breaking in a split potential box

Spontaneous symmetry breaking of self-trapped and leaky modes in quasi-double-well potentials

Stabilization of vortex beams in Kerr media by nonlinear absorption

Stable dipole solitons and soliton complexes in the nonlinear Schrodinger equation with periodically modulated nonlinearity

Two-dimensional $χ^{2}$ solitons generated by the downconversion of Airy waves

Two-dimensional solitons in dipolar Bose-Einstein condensates with spin-orbit coupling

Unidirectional transport of wave packets through tilted discrete breathers in nonlinear lattices with asymmetric defects

Vortex solitons in two-dimensional spin-orbit coupled Bose-Einstein condensates: effects of the Rashba-Dresselhaus coupling and the Zeeman splitting

Composite localized modes in discretized spin-orbit-coupled Bose-Einstein condensates

Correlations and synchronization in a Bose-Fermi mixture

CP symmetry in optical systems

Dark solitons in dual-core waveguides with dispersive coupling

Discrete solitons in self-defocusing systems with $\mathcal{PT}$-symmetric defects

Engineering Bright Solitons to Enhance the Stability of Two-Component Bose-Einstein Condensates

Gap solitons attached to a gapless layer