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<g_{c}^{(num)}≈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>g_{c}^{(num)}, the collapse threshold is found with the help of full numerical simulations. The relative increase of g_{c} above g_{c}^{(num)} is ≈1.5%. It is a meaningful result, even if its size is small, because the collapse threshold is a universal constant which is difficult to change.

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

Fermi-Pasta-Ulam phenomena and persistent breathers in the harmonic trap

We consider the long-term weakly nonlinear evolution governed by the two-dimensional nonlinear Schrödinger (NLS) equation with an isotropic harmonic oscillator potential. The dynamics in this regime is dominated by resonant interactions between quartets of linear normal modes, accurately captured by the corresponding resonant approximation. Within this approximation, we identify Fermi-Pasta-Ulam-like recurrence phenomena, whereby the normal-mode spectrum passes in close proximity of the initial configuration, and two-mode states with time-independent mode amplitude spectra that translate into long-lived breathers of the original NLS equation. We comment on possible implications of these findings for nonlinear optics and matter-wave dynamics in Bose-Einstein condensates.

Symmetry-breaking bifurcations and ghost states in the fractional nonlinear Schrödinger equation with a PT-symmetric potential

We report symmetry-breaking and restoring bifurcations of solitons in a fractional Schrödinger equation with cubic or cubic-quintic (CQ) nonlinearity and a parity-time-symmetric potential, which may be realized in optical cavities. Solitons are destabilized at the bifurcation point, and, in the case of CQ nonlinearity, the stability is restored by an inverse bifurcation. Two mutually conjugate branches of ghost states (GSs), with complex propagation constants, are created by the bifurcation, solely in the case of fractional diffraction. While GSs are not true solutions, direct simulations confirm that their shapes and results of their stability analysis provide a "blueprint" for the evolution of genuine localized modes in the system.

Modulational instability and soliton generation in chiral Bose-Einstein condensates with zero-energy nonlinearity

By means of analytical and numerical methods, we address the modulational instability (MI) in chiral condensates governed by the Gross-Pitaevskii equation including the current nonlinearity. The analysis shows that this nonlinearity partly suppresses the MI driven by the cubic self-focusing, although the current nonlinearity is not represented in the system's energy (although it modifies the momentum), hence it may be considered as zero-energy nonlinearity. Direct simulations demonstrate generation of trains of stochastically interacting chiral solitons by MI. In the ring-shaped setup, the MI creates a single traveling solitary wave. The sign of the current nonlinearity determines the direction of propagation of the emerging solitons.

Temporal cavity solitons in a laser-based microcomb: a path to a self-starting pulsed laser without saturable absorption

We theoretically present a design of self-starting operation of microcombs based on laser-cavity solitons in a system composed of a micro-resonator nested in and coupled to an amplifying laser cavity. We demonstrate that it is possible to engineer the modulational-instability gain of the system's zero state to allow the start-up with a well-defined number of robust solitons. The approach can be implemented by using the system parameters, such as the cavity length mismatch and the gain shape, to control the number and repetition rate of the generated solitons. Because the setting does not require saturation of the gain, the results offer an alternative to standard techniques that provide laser mode-locking.

Metastable soliton necklaces supported by fractional diffraction and competing nonlinearities

We demonstrate that the fractional cubic-quintic nonlinear Schrödinger equation, characterized by its Lévy index, maintains ring-shaped soliton clusters ("necklaces") carrying orbital angular momentum. They can be built, in the respective optical setting, as circular chains of fundamental solitons linked by a vortical phase field. We predict semi-analytically that the metastable necklace-shaped clusters persist, corresponding to a local minimum of an effective potential of interaction between adjacent solitons in the cluster. Systematic simulations corroborate that the clusters stay robust over extremely large propagation distances, even in the presence of strong random perturbations.

Creation and Characterization of Matter-Wave Breathers

We report the creation of quasi-1D excited matter-wave solitons, "breathers," by quenching the strength of the interactions in a Bose-Einstein condensate with attractive interactions. We characterize the resulting breathing dynamics and quantify the effects of the aspect ratio of the confining potential, the strength of the quench, and the proximity of the 1D-3D crossover for the two-soliton breather. Furthermore, we demonstrate the complex dynamics of a three-soliton breather created by a stronger interaction quench. Our experimental results, which compare well with numerical simulations, provide a pathway for utilizing matter-wave breathers to explore quantum effects in large many-body systems.

Double-layer Bose-Einstein condensates: A quantum phase transition in the transverse direction, and reduction to two dimensions

We revisit the problem of the reduction of the three-dimensional (3D) dynamics of Bose-Einstein condensates, under the action of strong confinement in one direction (z), to a 2D mean-field equation. We address this problem for the confining potential with a singular term, viz., V_{z}(z)=2z^{2}+ζ^{2}/z^{2}, with constant ζ. A quantum phase transition is induced by the latter term, between the ground state (GS) of the harmonic oscillator and the 3D condensate split in two parallel noninteracting layers, which is a manifestation of the "superselection" effect. A realization of the respective physical setting is proposed, making use of resonant coupling to an optical field, with the resonance detuning modulated along z. The reduction of the full 3D Gross-Pitaevskii equation (GPE) to the 2D nonpolynomial Schrödinger equation (NPSE) is based on the factorized ansatz, with the z -dependent multiplier represented by an exact GS solution of the 1D Schrödinger equation with potential V_{z}(z). For both repulsive and attractive signs of the nonlinearity, the 2D NPSE produces GS and vortex states, that are virtually indistinguishable from the respective numerical solutions provided by full 3D GPE. In the case of the self-attraction, the threshold for the onset of the collapse, predicted by the 2D NPSE, is also virtually identical to its counterpart obtained from the 3D equation. In the same case, stability and instability of vortices with topological charge S=1, 2, and 3 are considered in detail. Thus, the procedure of the spatial-dimension reduction, 3D → 2D, produces very accurate results, and it may be used in other settings.

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 stream 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 on the system of coupled nonlinear Schrödinger equations provides surprisingly good agreement with the experimental findings.

Quantum Fluctuations of the Center of Mass and Relative Parameters of Nonlinear Schrödinger Breathers

We study quantum fluctuations of macroscopic parameters of a nonlinear Schrödinger breather-a nonlinear superposition of two solitons, which can be created by the application of a fourfold quench of the scattering length to the fundamental soliton in a self-attractive quasi-one-dimensional Bose gas. The fluctuations are analyzed in the framework of the Bogoliubov approach in the limit of a large number of atoms N, using two models of the vacuum state: white noise and correlated noise. The latter model, closer to the ab initio setting by construction, leads to a reasonable agreement, within 20% accuracy, with fluctuations of the relative velocity of constituent solitons obtained from the exact Bethe-ansatz results [Phys. Rev. Lett. 119, 220401 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.220401] in the opposite low-N limit (for N≤23). We thus confirm, for macroscopic N, the breather dissociation time to be within the limits of current cold-atom experiments. Fluctuations of soliton masses, phases, and positions are also evaluated and may have experimental implications.

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 are investigated by combining techniques of 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 us to address fundamental problems of the quantum speed limit and thermodynamics.

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

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, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schrödinger 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 δ-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.

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.

Self-interaction of ultrashort pulses in an epsilon-near-zero nonlinear material at the telecom wavelength

Dynamics of femtosecond pulses with the telecom carrier wavelength is investigated numerically in a subwavelength layer of an indium tin oxide (ITO) epsilon-near-zero (ENZ) material with high dispersion and high nonlinearity. Due to the subwavelength thickness of the ITO ENZ material, and the fact that the pulse's propagation time is shorter than its temporal width, multiple reflections give rise to self-interaction in both spectral and temporal domains, especially at wavelengths longer than at the ENZ point, at which the reflections are significantly stronger. A larger absolute value of the pulse's chirp strongly affects the self-interaction by redistributing energy between wavelengths, while the sign of the chirp affects the interaction in the temporal domain. It is also found that, when two identical pulses are launched simultaneously from both ends, a subwavelength counterpart of a standing-wave state can be established. It shows robust energy localization in the middle of the sample, in terms of both the spectral and temporal intensity distributions.

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 in condensed matter, cosmology, and high-energy physics.

Nonlinear anti-directional couplers with gain and loss

Following the concept of PT-symmetric couplers, we propose a linearly coupled system of nonlinear waveguides, made of positive- and negative-index materials, which carry, respectively, gain and loss. We report novel bi- and multistability states pertaining to transmitted and reflective intensities, which are controlled by the ratio of the gain and loss coefficients, and phase mismatch between the waveguides. These states offer transmission regimes with extremely low threshold intensities for transitions between coexisting states, and very large amplification ratios between the input and output intensities leading to an efficient way of controlling light with light.

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.

Semidiscrete Quantum Droplets and Vortices

We consider a binary bosonic condensate with weak mean-field (MF) residual repulsion, loaded in an array of nearly one-dimensional traps coupled by transverse hopping. With the MF force balanced by the effectively one-dimensional attraction, induced in each trap by the Lee-Hung-Yang correction (produced by quantum fluctuations around the MF state), stable on-site- and intersite-centered semidiscrete quantum droplets (QDs) emerge in the array, as fundamental ones and self-trapped vortices, with winding numbers, at least, up to five, in both tightly bound and quasicontinuum forms. The application of a relatively strong trapping potential leads to squeezing transitions, which increase the number of sites in fundamental QDs and eventually replace vortex modes by fundamental or dipole ones. The results provide the first realization of stable semidiscrete vortex QDs, including ones with multiple vorticity.

Observation of incoherently coupled dark-bright vector solitons in single-mode fibers

We report experimental observation of incoherently coupled dark-bright vector solitons in single-mode fibers. Properties of the vector solitons accord well with those predicted by the respective systems of incoherently coupled nonlinear Schrödinger equations. To our knowledge, this is the first experimental observation of temporal incoherently coupled dark-bright solitons in single-mode fibers.

Attraction centers and parity-time-symmetric delta-functional dipoles in critical and supercritical self-focusing media

We introduce a model based on the one-dimensional nonlinear Schrödinger equation with critical (quintic) or supercritical self-focusing nonlinearity. We demonstrate that a family of solitons, which are unstable in this setting against the critical or supercritical collapse, is stabilized by pinning to an attractive defect, that may also include a parity-time (PT)-symmetric gain-loss component. The model can be realized as a planar waveguide in nonlinear optics, and in a super-Tonks-Girardeau bosonic gas. For the attractive defect with the delta-functional profile, a full family of the pinned solitons is found in an exact analytical form. In the absence of the gain-loss term, the solitons' stability is investigated in an analytical form too, by means of the Vakhitov-Kolokolov criterion; in the presence of the PT-balanced gain and loss, the stability is explored by means of numerical methods. In particular, the entire family of pinned solitons is stable in the quintic (critical) medium if the gain-loss term is absent. A stability region for the pinned solitons persists in the model with an arbitrarily high power of the self-focusing nonlinearity. A weak gain-loss component gives rise to intricate alternations of stability and instability in the system's parameter plane. Those solitons which are unstable under the action of the supercritical self-attraction are destroyed by the collapse. On the other hand, if the self-attraction-driven instability is weak and the gain-loss term is present, unstable solitons spontaneously transform into localized breathers, while the collapse does not occur. The same outcome may be caused by a combination of the critical nonlinearity with the gain and loss. Instability of the solitons is also possible when the PT-symmetric gain-loss term is added to the subcritical nonlinearity. The system with self-repulsive nonlinearity is briefly considered too, producing completely stable families of pinned localized states.

Metastability of Quantum Droplet Clusters

We show that metastable ring-shaped clusters can be constructed from two-dimensional quantum droplets in systems described by the Gross-Pitaevskii equations augmented with Lee-Huang-Yang quantum corrections. The clusters exhibit dynamical behavior ranging from contraction to rotation with simultaneous periodic pulsations, or expansion, depending on the initial radius of the necklace pattern and phase shift between adjacent quantum droplets. We show that, using an energy-minimization analysis, one can predict equilibrium values of the cluster radius that correspond to rotation without radial pulsations. In such a regime, the clusters evolve as metastable states, withstanding abrupt variations in the underlying scattering lengths and keeping their azimuthal symmetry in the course of evolution, even in the presence of considerable perturbations.

Interactions of solitons with positive and negative masses: Shuttle motion and coacceleration

We consider a possibility to realize self-accelerating motion of interacting states with effective positive and negative masses in the form of pairs of solitons in two-component BEC loaded in an optical-lattice (OL) potential. A crucial role is played by the fact that gap solitons may feature a negative dynamical mass, keeping their mobility in the OL. First, the respective system of coupled Gross-Pitaevskii equations (GPE) is reduced to a system of equations for envelopes of the lattice wave functions. Two generic dynamical regimes are revealed by simulations of the reduced system, viz., shuttle oscillations of pairs of solitons with positive and negative masses, and splitting of the pair. The coaccelerating motion of the interacting solitons, which keeps constant separation between them, occurs at the boundary between the shuttle motion and splitting. The position of the coacceleration regime in the system's parameter space can be adjusted with the help of an additional gravity potential, which induces its own acceleration, that may offset the relative acceleration of the two solitons, while gravity masses of both solitons remain positive. The numerical findings are accurately reproduced by a variational approximation. Collisions between shuttling or coaccelerating soliton pairs do not alter the character of the dynamical regime. Finally, regimes of the shuttle motion, coacceleration, and splitting are corroborated by simulations of the original GPE system, with the explicitly present OL potential.