In this paper we study the solutions to a nonlinear Schrödinger equation with harmonic and periodic potentials, motivated from the recent interest in these models as mean-field descriptions of Bose–Einstein condensates. We use a two-mode Galerkin approximation to study the dynamics of the full model. The phase plane and stability analysis of the reduced model yield very good agreement with the findings of the full partial differential equation. A particularly interesting finding of the stability analysis is a spontaneous symmetry breaking through a branching bifurcation, resulting in the stabilization of asymmetric states and the destabilization of symmetric or anti-symmetric ones. We also highlight the important differences between the cases of symmetric potentials and those of weakly asymmetric potentials.
In this paper we study the solutions to a nonlinear Schrödinger equation with harmonic and periodic potentials, motivated from the recent interest in these models as mean-field descriptions of Bose–Einstein condensates. We use a two-mode Galerkin approximation to study the dynamics of the full model. The phase plane and stability analysis of the reduced model yield very good agreement with the findings of the full partial differential equation. A particularly interesting finding of the stability analysis is a spontaneous symmetry breaking through a branching bifurcation, resulting in the stabilization of asymmetric states and the destabilization of symmetric or anti-symmetric ones. We also highlight the important differences between the cases of symmetric potentials and those of weakly asymmetric potentials.
We study the existence and stability of solutions of the two-dimensional nonlinear Schrodinger equation in the combined presence of a parabolic and a periodic potential. The motivating physical example consists of Bose-Einstein condensates confined in an harmonic (e.g., magnetic) trap and an optical lattice. By connecting the nonlinear problem with the underlying linear spectrum, we examine the bifurcation of nonlinear modes out of the linear ones for both focusing and defocusing nonlinearities. In particular, we find real-valued solutions (such as multipoles) and complex-valued ones (such as vortices). A primary motivation of the present work is to develop "rules of thumb" about what waveforms to expect emerging in the nonlinear problem and about the stability of those modes. As a case example of the latter, we find that among the real-valued solutions, the one with larger norm for a fixed value of the chemical potential is expected to be unstable.
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.
In this paper we consider solutions of a nonlinear Schrodinger equation with a parabolic and a periodic potential motivated from the dynamics of Bose-Einstein condensates. Our starting point is the corresponding linear problem which we analyze through regular perturbation and homogenization techniques. We then use Lyapunov-Schmidt theory to establish the persistence and bifurcation of the linear states in the presence of attractive and repulsive nonlinear inter-particle interactions. Stability of such solutions is also examined and a count is given of the potential real, complex and imaginary eigenvalues with negative Krein signature that such solutions may possess. The results are corroborated with numerical computations.
We consider the statics and dynamics of F = 1 spinor Bose–Einstein condensates (BECs) confined in double-well potentials. We use a two-mode Galerkin-type quasi-analytical approximation to describe the stationary states of the system. This way, we are able to obtain not only earlier results based on the single-mode approximation (SMA) frequently used in studies of spinor BECs, but also additional modes that involve either two or all three spinor components of the F = 1 spinor BEC. The results based on this Galerkin-type decomposition are in good agreement with the analysis of the full system. We subsequently analyze the stability of these multi-component states, as well as their dynamics when we find them to be unstable. The instabilities of the symmetric or anti-symmetric states exhibit symmetry-breaking and recurrent asymmetric patterns. Our results yield qualitatively similar bifurcation diagrams both for polar (such as 23 Na) and ferromagnetic (such as 87 Rb) spinor BECs.
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 treat the behavior of Bose-Einstein condensates in double square well potentials of both equal and different depths. For even depth, symmetry preserving solutions to the relevant nonlinear Schrödinger equation are known, just as in the linear limit. When the nonlinearity is strong enough, symmetry breaking solutions also exist, side by side with the symmetric one. Interestingly, solutions almost entirely localized in one of the wells are known as an extreme case. Here we outline a method for obtaining all these solutions for repulsive interactions. The bifurcation point at which, for critical nonlinearity, the asymmetric solutions branch off from the symmetry preserving ones is found analytically. We also find this bifurcation point and treat the solutions generally via a Josephson junction model. When the confining potential is in the form of two wells of different depth, interesting phenomena appear. This is true of both the occurrence of the bifurcation point for the static solutions and also of the dynamics of phase and amplitude varying solutions. Again a generalization of the Josephson model proves useful. The stability of solutions is treated briefly.