No Arabic abstract
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 systems main (semi-infinite) and annex (finite) bandgaps. Tilted (moving) solitons are unstable, spontaneously evolving into robust breathers. For broad solitons, diffraction, represented by second derivatives in the system, may be neglected, leading to a simplified model with a finite bandgap. It is populated by skew-antisymmetric gap solitons, which are nearly stable close to the gaps bottom.
We report results of the analysis for families of one-dimensional (1D) trapped solitons, created by competing self-focusing (SF) quintic and self-defocusing (SDF) cubic nonlinear terms. Two trapping potentials are considered, the harmonic-oscillator (HO) and delta-functional ones. The models apply to optical solitons in colloidal waveguides and other photonic media, and to matter-wave solitons in Bose-Einstein condensates (BEC) loaded into a quasi-1D trap. For the HO potential, the results are obtained in an approximate form, using the variational and Thomas-Fermi approximations (VA and TFA), and in a full numerical form, including the ground state and the first antisymmetric excited one. For the delta-functional attractive potential, the results are produced in a fully analytical form, and verified by means of numerical methods. Both exponentially localized solitons and weakly localized trapped modes are found for the delta-functional potential. The most essential conclusions concern the applicability of competing Vakhitov-Kolokolov (VK) and anti-VK criteria to the identification of the stability of solitons created under the action of the competing SF and SDF terms.
We consider the dynamical model of a binary bosonic gas trapped in a symmetric dual-core cigar-shaped potential. The setting is modeled by a system of linearly-coupled one-dimensional Gross-Pitaevskii equations with the cubic self-repulsive terms and quadratic attractive ones,which represent the Lee-Huang-Yang corrections to the mean-field theory in this geometry. The main subject is spontaneous symmetry breaking (SSB) of quantum droplets (QDs), followed by restoration of the symmetry, with respect to the symmetric parallel-coupled trapping cores, following the increase of the QDs total norm. The SSB transition and inverse symmetry-restoring one form a bifurcation loop, whose shape in concave at small values of the inter-core coupling constant, $kappa$, and convex at larger $kappa $. The loop does not exist above a critical value of $kappa $. At very large values of the norm, QDs do not break their symmetry, featuring a flat-top shape. Some results are obtained in an analytical form, including an exact front solution connecting constant zero and finite values of the wave function. Collisions between moving QDs are considered too, demonstrating a trend to merger into breathers.
We study coupled unstaggered-staggered soliton pairs emergent from a system of two coupled discrete nonlinear Schr{o}dinger (DNLS) equations with the self-attractive on-site self-phase-modulation nonlinearity, coupled by the repulsive cross-phase-modulation interaction, on 1D and 2D lattice domains. These mixed modes are of a symbiotic type, as each component in isolation may only carry ordinary unstaggered solitons. While most work on DNLS systems addressed symmetric on-site-centered fundamental solitons, these models give rise to a variety of other excited states, which may also be stable. The simplest among them are antisymmetric states in the form of discrete twisted solitons, which have no counterparts in the continuum limit. In the extension to 2D lattice domains, a natural counterpart of the twisted states are vortical solitons. We first introduce a variational approximation (VA) for the solitons, and then correct it numerically to construct exact stationary solutions, which are then used as initial conditions for simulations to check if the stationary states persist under time evolution. Two-component solutions obtained include (i) 1D fundamental-twisted and twisted-twisted soliton pairs, (ii) 2D fundamental-fundamental soliton pairs, and (iii) 2D vortical-vortical soliton pairs. We also highlight a variety of other transient dynamical regimes, such as breathers and amplitude death. The findings apply to modeling binary Bose-Einstein condensates, loaded in a deep lattice potential, with identical or different atomic masses of the two components, and arrays of bimodal optical waveguides.
We introduce a dual-core system with double symmetry, one between the cores, and one along each core, imposed by the spatial modulation of local nonlinearity in the form of two tightly localized spots, which may be approximated by a pair of ideal delta-functions. The analysis aims to investigate effects of spontaneous symmetry breaking in such systems. Stationary one-dimensional modes are constructed in an implicit analytical form. These solutions include symmetric ones, as well as modes with spontaneously broken inter-core and along-the-cores symmetries. Solutions featuring the simultaneous (double) breaking of both symmetries are produced too. In the model with the ideal delta-functions, all species of the asymmetric modes are found to be unstable. However, numerical consideration of a two dimensional extension of the system, which includes symmetric cores with a nonzero transverse thickness, and the nonlinearity-localization spots of a small finite size, produces stable asymmetric modes of all the types, realizing the separate breaking of each symmetry, and states featuring simultaneous (double) breaking of both symmetries.
Nonlinear periodic systems, such as photonic crystals and Bose-Einstein condensates (BECs) loaded into optical lattices, are often described by the nonlinear Schrodinger/Gross-Pitaevskii equation with a sinusoidal potential. Here, we consider a model based on such a periodic potential, with the nonlinearity (attractive or repulsive) concentrated either at a single point or at a symmetric set of two points, which are represented, respectively, by a single {delta}-function or a combination of two {delta}-functions. This model gives rise to ordinary solitons or gap solitons (GSs), which reside, respectively, in the semi-infinite or finite gaps of the systems linear spectrum, being pinned to the {delta}-functions. Physical realizations of these systems are possible in optics and BEC, using diverse variants of the nonlinearity management. First, we demonstrate that the single {delta}-function multiplying the nonlinear term supports families of stable regular solitons in the self-attractive case, while a family of solitons supported by the attractive {delta}-function in the absence of the periodic potential is completely unstable. We also show that the {delta}-function can support stable GSs in the first finite gap in both the self-attractive and repulsive models. The stability analysis for the GSs in the second finite gap is reported too, for both signs of the nonlinearity. Alongside the numerical analysis, analytical approximations are developed for the solitons in the semi-infinite and first two finite gaps, with the single {delta}-function positioned at a minimum or maximum of the periodic potential. In the model with the symmetric set of two {delta}-functions, we study the effect of the spontaneous symmetry breaking of the pinned solitons. Two configurations are considered, with the {delta}-functions set symmetrically with respect to the minimum or maximum of the potential.