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.
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π), the bifurcation does not happen, as branches of different solutions asymptotically approach each other, but fail to merge.
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 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 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.
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 elaborate a mechanism for the formation of stable solitons of the semivortex type (with vorticities 0 and 1 in their two components), populating a finite band gap in the spectrum of the spin-orbit-coupled binary Bose-Einstein condensate with the Zeeman splitting, in the two-dimensional (2D) free space, under conditions which make the kinetic-energy terms in the respective coupled Gross-Pitaevskii equations negligible. Unlike a recent work which used long-range dipole-dipole interactions to construct stable gap solitons in a similar setting, we here demonstrate that stable solitons are supported by generic local interactions of both attractive and repulsive signs, provided that the relative strength of the cross- and self-interactions in the two-component system does not exceed a critical value $\ensuremath{\approx}0.77$. A boundary between stable and unstable fundamental 2D gap solitons is precisely predicted by the Vakhitov-Kolokolov criterion, while all excited states of the 2D solitons, with vorticities (m,1+m) in the two components, m=1,2,..., are unstable. The analysis of the one-dimensional (1D) reduction of the system produces an exact analytical solution for the family of gap solitons which populate the entire band gap, the family being fully stable. Motion of the 1D solitons in the trapping potential is considered too, showing that their effective mass is positive or negative if the cubic nonlinearity is attractive or repulsive, respectively.
We consider one-dimensional solitons in a binary Bose-Einstein condensate with linear coupling between the components, trapped in an optical-lattice potential. The inter- and intraspecies interactions may be either repulsive or attractive. The main effects considered here are spontaneous breaking of the symmetry between components in symmetric and antisymmetric solitons, and spatial splitting between the components. These effects are studied by means of a variational approximation and numerical simulations.