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.

Robust [Formula: see text] symmetry of two-dimensional fundamental and vortex solitons supported by spatially modulated nonlinearity

The real spectrum of bound states produced by [Formula: see text]-symmetric Hamiltonians usually suffers breakup at a critical value of the strength of gain-loss terms, i.e., imaginary part of the complex potential. The breakup essentially impedes the use of [Formula: see text]-symmetric systems for various applications. On the other hand, it is known that the [Formula: see text] symmetry can be made unbreakable in a one-dimensional (1D) model with self-defocusing nonlinearity whose strength grows fast enough from the center to periphery. The model is nonlinearizable, i.e., it does not have a linear spectrum, while the (unbreakable) [Formula: see text] symmetry in it is defined by spectra of continuous families of nonlinear self-trapped states (solitons). Here we report results for a 2D nonlinearizable model whose [Formula: see text] symmetry remains unbroken for arbitrarily large values of the gain-loss coefficient. Further, we introduce an extended 2D model with the imaginary part of potential ~xy in the Cartesian coordinates. The latter model is not a [Formula: see text]-symmetric one, but it also supports continuous families of self-trapped states, thus suggesting an extension of the concept of the [Formula: see text] symmetry. For both models, universal analytical forms are found for nonlinearizable tails of the 2D modes, and full exact solutions are produced for particular solitons, including ones with the unbreakable [Formula: see text] symmetry, while generic soliton families are found in a numerical form. The [Formula: see text]-symmetric system gives rise to generic families of stable single- and double-peak 2D solitons (including higher-order radial states of the single-peak solitons), as well as families of stable vortex solitons with m = 1, 2, and 3. In the model with imaginary potential ~xy, families of single- and multi-peak solitons and vortices are stable if the imaginary potential is subject to spatial confinement. In an elliptically deformed version of the latter model, an exact solution is found for vortex solitons with m = 1.

Purely Kerr nonlinear model admitting flat-top solitons

We elaborate one- and two-dimensional (1D and 2D) models of media with self-repulsive cubic nonlinearity, whose local strength is subject to spatial modulation that admits the existence of flat-top solitons of various types, including fundamental ones, 1D multipoles, and 2D vortices. Previously, solitons of this type were only produced by models with competing nonlinearities. The present setting may be implemented in optics and Bose-Einstein condensates. The 1D version gives rise to an exact analytical solution for stable flat-top solitons, and generic families may be predicted by means of the Thomas-Fermi approximation. Stability of the obtained flat-top solitons is analyzed by means of the linear-stability analysis and direct simulations. Fundamental solitons and 1D multipoles with k=1 and 2 nodes, as well as vortices with winding number m=1, are completely stable. For multipoles with k≥3 and vortices with m≥2, alternating stripes of stability and instability are identified in their parameter spaces.

Optical solitons in media with focusing and defocusing saturable nonlinearity and a parity-time-symmetric external potential

We report results for solitons in models of waveguides with focusing or defocusing saturable nonlinearity and a parity-time ([Formula: see text])-symmetric complex-valued external potential of the Scarf-II type. The model applies to the nonlinear wave propagation in graded-index optical waveguides with balanced gain and loss. We find both fundamental and multipole solitons for both focusing and defocusing signs of the saturable nonlinearity in such [Formula: see text]-symmetric waveguides. The dependence of the propagation constant on the soliton's power is presented for different strengths of the nonlinearity saturation, S The stability of fundamental, dipole, tripole and quadrupole solitons is investigated by means of the linear-stability analysis and direct numerical simulations of the corresponding (1+1)-dimensional nonlinear Schrödinger-type equation. The results show that the instability of the stationary solutions can be mitigated or completely suppressed, increasing the value of SThis article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 1)'.

𝒫𝒯-symmetric and antisymmetric nonlinear states in a split potential box

We introduce a one-dimensional [Formula: see text]-symmetric system, which includes the cubic self-focusing, a double-well potential in the form of an infinitely deep potential box split in the middle by a delta-functional barrier of an effective height ε, and constant linear gain and loss, γ, in each half-box. The system may be readily realized in microwave photonics. Using numerical methods, we construct [Formula: see text]-symmetric and antisymmetric modes, which represent, respectively, the system's ground state and first excited state, and identify their stability. Their instability mainly leads to blowup, except for the case of ε=0, when an unstable symmetric mode transforms into a weakly oscillating breather, and an unstable antisymmetric mode relaxes into a stable symmetric one. At ε>0, the stability area is much larger for the [Formula: see text]-antisymmetric state than for its symmetric counterpart. The stability areas shrink with increase of the total power, P In the linear limit, which corresponds to [Formula: see text], the stability boundary is found in an analytical form. The stability area of the antisymmetric state originally expands with the growth of γ, and then disappears at a critical value of γThis article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 1)'.