Spontaneous symmetry breaking of fundamental states, vortices, and dipoles in two- and one-dimensional linearly coupled traps with cubic self-attraction
Zhaopin Chen, Yongyao Li, Boris A. Malomed, Luca Salasnich
We introduce two- and one-dimensional (2D and 1D) systems of two linearly coupled Gross-Pitaevskii equations (GPEs) with the cubic self-attraction and harmonic-oscillator (HO) trapping potential in each GPE. The system models a Bose-Einstein condensate with a negative scattering length, loaded in a double-pancake trap, combined with the in-plane HO potential. In addition to that, the 1D version applies to the light transmission in a dual-core waveguide with the Kerr nonlinearity and in-core confinement represented by the HO potential. The subject of the analysis is spontaneous symmetry breaking in 2D and 1D ground-state (GS, alias fundamental) modes, as well as in 2D vortices and 1D dipole modes. (The latter ones do not exist without the HO potential.) By means of the variational approximation and numerical analysis, it is found that both the 2D and 1D systems give rise to a symmetry-breaking bifurcation (SBB) of the supercritical type. The stability of symmetric and asymmetric states, produced by the SBB, is analyzed through the computation of eigenvalues for perturbation modes and verified by direct simulations. The asymmetric GSs are always stable, while the stability region for vortices shrinks and eventually disappears with the increase of the linear-coupling constant, $\ensuremath{\kappa}$. The SBB in the 2D system does not occur if $\ensuremath{\kappa}$ is too large (at $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{max}$); in that case, the two-component system behaves, essentially, as its single-component counterpart. In the 1D system, both asymmetric and symmetric dipole modes feature an additional oscillatory instability, unrelated to the symmetry breaking. This instability occurs in several regions which expand with the increase of $\ensuremath{\kappa}$.
We introduce two- and one-dimensional (2D and 1D) systems of two linearly coupled Gross-Pitaevskii equations (GPEs) with the cubic self-attraction and harmonic-oscillator (HO) trapping potential in each GPE. The system models a Bose-Einstein condensate with a negative scattering length, loaded in a double-pancake trap, combined with the in-plane HO potential. In addition to that, the 1D version applies to the light transmission in a dual-core waveguide with the Kerr nonlinearity and in-core confinement represented by the HO potential. The subject of the analysis is spontaneous symmetry breaking in 2D and 1D ground-state (GS, alias fundamental) modes, as well as in 2D vortices and 1D dipole modes. (The latter ones do not exist without the HO potential.) By means of the variational approximation and numerical analysis, it is found that both the 2D and 1D systems give rise to a symmetry-breaking bifurcation (SBB) of the supercritical type. The stability of symmetric and asymmetric states, produced by the SBB, is analyzed through the computation of eigenvalues for perturbation modes and verified by direct simulations. The asymmetric GSs are always stable, while the stability region for vortices shrinks and eventually disappears with the increase of the linear-coupling constant, \ensuremathκ. The SBB in the 2D system does not occur if \ensuremathκ is too large (at \ensuremathκ>\ensuremathκmax); in that case, the two-component system behaves, essentially, as its single-component counterpart. In the 1D system, both asymmetric and symmetric dipole modes feature an additional oscillatory instability, unrelated to the symmetry breaking. This instability occurs in several regions which expand with the increase of \ensuremathκ.
Models of two-dimensional (2D) traps, with double-well structure in the third direction, for Bose-Einstein condensates are introduced with attractive or repulsive interactions between atoms. The models are based on systems of linearly coupled 2D Gross-Pitaevskii equations, where the coupling accounts for tunneling between the wells. Each well carries an optical lattice (OL) (stable 2D solitons cannot exist without OLs). The linear coupling splits each finite band gap in the spectrum of the single-component model into two subgaps. The main subject of the work is spontaneous symmetry breaking (SSB) in two-component 2D solitons and localized vortices (SSB was not considered before in 2D settings). Using variational approximation (VA) and numerical methods, we demonstrate that, in a system with attraction or repulsion, SSB occurs in families of symmetric or antisymmetric solitons (or vortices), respectively. The corresponding bifurcation destabilizes the original solution branch and gives rise to a stable branch of asymmetric solitons or vortices. The VA provides for an accurate description of the emerging branch of asymmetric solitons. In the model with attraction, all stable branches eventually terminate due to the onset of collapse. Stable asymmetric solitons in higher finite band gaps and vortices with a multiple topological charge are found too. The models also give rise to first examples of embedded solitons and embedded vortices (the states located inside Bloch bands) in two dimensions. In the linearly coupled system with opposite signs of the nonlinearity in the two cores, two distinct types of stable solitons and vortices are found, dominated by either the self-attractive component or the self-repulsive one. In the system with a mismatch between the two OLs, a pseudobifurcation is found: when the mismatch attains its largest value $(\ensuremath{\pi})$, the bifurcation does not happen, as branches of different solutions asymptotically approach each other, but fail to merge.
We consider a two-component linearly coupled system with the intrinsic cubic nonlinearity and the harmonic-oscillator (HO) confining potential. The system models binary settings in BEC and optics. In the symmetric system, with the HO trap acting in both components, we consider Josephson oscillations (JO) initiated by an input in the form of the HO’s ground state (GS) or dipole mode (DM), placed in one component. With the increase of the strength of the self-focusing nonlinearity, spontaneous symmetry breaking (SSB) between the components takes place in the dynamical JO state. Under still stronger nonlinearity, the regular JO initiated by the GS input carries 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.
We consider a two-component linearly-coupled system with the intrinsic cubic nonlinearity and the harmonic-oscillator (HO) confining potential. The system models binary settings in BEC and optics. In the symmetric system, with the HO trap acting in both components, we consider Josephson oscillations (JO) initiated by an input in the form of the HO's ground state (GS) or dipole mode (DM), placed in one component. With the increase of the strength of the self-focusing nonlinearity, spontaneous symmetry breaking (SSB) between the components takes place in the dynamical JO state. Under still stronger nonlinearity, the regular JO initiated by the GS input carry over into a chaotic dynamical state. For the DM input, the chaotization happens at smaller powers than for the GS, which is followed by SSB at a slightly stronger nonlinearity. In the system with the defocusing nonlinearity, SSB does not take place, and dynamical chaos occurs in a small area of the parameter space. In the asymmetric half-trapped system, with the HO potential applied to a single component, we first focus on the spectrum of confined binary modes in the linearized system. The spectrum is found analytically in the limits of weak and strong inter-component coupling, and numerically in the general case. Under the action of the coupling, the existence region of the confined modes shrinks for GSs and expands for DMs. In the full nonlinear system, the existence region for confined modes is identified in the numerical form. They are constructed too by means of the Thomas-Fermi approximation, in the case of the defocusing nonlinearity. Lastly, particular (non-generic) exact analytical solutions for confined modes, including vortices, in one- and two-dimensional asymmetric linearized systems are found. They represent bound states in the continuum.
We study spontaneous symmetry breaking in a system of two parallel quasi-one-dimensional traps (cores), equipped with optical lattices (OLs) and filled with a Bose-Einstein condensate (BEC). The cores are linearly coupled by tunneling (the model may also be interpreted in terms of spatial solitons in parallel planar optical waveguides with a periodic modulation of the refractive index). Analysis of the corresponding system of linearly coupled Gross-Pitaevskii equations (GPEs) reveals that spectral band gaps of the single GPE split into subgaps. Symmetry breaking in two-component BEC solitons is studied in cases of the attractive (AA) and repulsive (RR) nonlinearity in both traps; the mixed situation, with repulsion in one trap and attraction in the other (RA), is considered too. In all the cases, stable asymmetric solitons are found, bifurcating from symmetric or antisymmetric ones (and destabilizing them), in the AA and RR systems, respectively. In either case, bistability is predicted, with a nonbifurcating stable branch, either antisymmetric or symmetric, coexisting with asymmetric ones. Solitons destabilized by the bifurcation tend to rearrange themselves into their stable asymmetric counterparts. In addition to the fundamental solitons, branches of twisted (odd) solitons in the AA system, and twisted bound states of fundamental solitons in both AA and RR systems, are found too. The impact of a phase mismatch, $\ensuremath{\Delta}$, between the OLs in the two cores is also studied. It is concluded that $\ensuremath{\Delta}=\ensuremath{\pi}∕2$ only mildly deforms the picture, while $\ensuremath{\Delta}=\ensuremath{\pi}$ changes it drastically, replacing the symmetry-breaking bifurcations by pseudobifurcations, with the branch of asymmetric solutions asymptotically approaching its symmetric or antisymmetric counterpart (in the AA and RR system, respectively), rather than splitting off from it. Also considered is a related model, for a binary BEC in a single-core trap with the OL, assuming that the two species (representing different spin states of the same atom) are coupled by linear interconversion. In that case, the symmetry-breaking bifurcations in the AA and RR models switch their character, if the interspecies nonlinear interaction becomes stronger than the intraspecies nonlinearity.
We report results of an analysis of the spontaneous symmetry breaking (SSB) in a basic (actually, simplest) model that is capable of producing the SSB phenomenology in a one-dimensional setting. It is based on the Gross-Pitaevskii-nonlinear Schrödinger equation with the cubic self-attractive term and a double-well potential built as an infinitely deep potential box split by a narrow (δ functional) barrier. The barrier's strength ɛ is the single free parameter of the scaled form of the model. It may be implemented in atomic Bose-Einstein condensates and nonlinear optics. The SSB bifurcation of the symmetric ground state (g.s.) is predicted analytically in two limit cases, viz., for deep or weak splitting of the potential box by the barrier (ɛ≫1 or ɛ≪1, respectively). For the generic case, a variational approximation (VA) is elaborated. The analytical findings are presented along with systematic numerical results. The stability of stationary states is studied through the calculation of eigenvalues for small perturbations and by means of direct simulations. The g.s. always undergoes the SSB bifurcation of the supercritical type, as predicted by the VA at moderate values of ɛ, although the VA fails at small ɛ, due to inapplicability of the underlying ansatz in that case. However, the latter case is correctly treated by the approximation based on a soliton ansatz. On top of the g.s., the first and second excited states are studied too. The antisymmetric mode (the first excited state) is destabilized at a critical value of its norm. The second excited state undergoes SSB bifurcation, like the g.s., but, unlike it, the bifurcation produces an unstable asymmetric mode. All unstable modes tend to spontaneously reshape into the asymmetric g.s.
We study stationary states of a two-dimensional (2D) Bose-Einstein condensate with both attractive and repulsive nonlinearities in a combination of a double-square-well potential in one direction and a perpendicular optical lattice. We look for dual-core solitons in this configuration, focusing on their symmetry-breaking bifurcations. For attractive interactions, without the lattice, a similar analysis was performed [M. Matuszewski et al., Phys. Rev. A 75, 063621 (2007)], where subcritical bifurcation transforming antisymmetric gap solitons into asymmetric ones was found. Here we focus on the effect of an optical lattice and the so created gap solitons. We discover that a phase transition occurs when the lattice depth increases and the additional dimension becomes strongly suppressed. The bifurcation type changes from subcritical (typical for a 2D system, with hysteresis) to supercritical (typical for a one-dimensional system). An additional advantage of the lattice is that gap solitons exist even for repulsive interactions. In this case we also discover bifurcation of a supercritical type. The analysis is based on a variational approximation, which is surprisingly well verified by numerical results.