Researcher profile

Boris A. Malomed

Boris A. Malomed contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
50works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

50 published item(s)

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&#39;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 &#34;debris&#34; 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&#39;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&#39;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&#39;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&#39;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&#39;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&#39;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&#39;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 &#34;bare&#34; 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 (&#34;bare charge&#34;) 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&#39;s main (semi-infinite) and annex (finite) bandgaps. Tilted (&#34;moving&#34;) 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&#39;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&#39;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 &#34;symbiotic&#34; 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&#39;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 &#39;minimum loss principle,&#39; 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&#39;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&#39; dynamics. In particular, dependence of stability domains in the model&#39;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.