We analyze gap solitons in trapped Bose-Einstein condensates (BECs) in optical lattice potentials under Feshbach resonance management. Starting with an averaged Gross-Pitaevsky (GP) equation with a periodic potential, we employ an envelope wave approximation to derive coupled-mode equations describing the slow BEC dynamics in the first spectral gap of the optical lattice. We construct exact analytical formulas describing gap soliton solutions and examine their spectral stability using the Chebyshev interpolation method. We show that these gap solitons are unstable far from the threshold of local bifurcation and that the instability results in the distortion of their shape. We also predict the threshold of the power of gap solitons near the local bifurcation limit.
INTRODUCTIONAt sufficiently low temperatures, particles in a dilute boson gas can condense in the ground state, forming a Bose-Einstein condensate (BEC) [1]. Under the typical confining conditions of experimental settings, BECs are inhomogeneous and the number of condensed atoms (N ) ranges from several thousand (or less) to several million (or more). The magnetic traps that confine them are usually approximated by harmonic potentials. There are two characteristic length scales: the harmonic oscillator length a ho = /(mω ho ) [which is on the order of a few microns], where ω ho = (ω x ω y ω z ) 1/3 is the geometric mean of the trapping frequencies, and the mean healing length χ = 1/ 8π|a|n [which is on the order of a micron], wheren is the mean particle density and a, the (two-body) s-wave scattering length, is determined by the atomic species of the condensate. Interactions between atoms are repulsive when a > 0 and attractive when a < 0. For a dilute ideal gas, a ≈ 0.If considering only two-body, mean-field interactions, a dilute Bose-Einstein gas can be modeled using a cubic nonlinear Schrödinger equation (NLS) with an external potential; this is also known as the Gross-Pitaevsky (GP) equation. BECs are modeled in the quasi-onedimensional (quasi-1D) regime when the transverse dimensions of the condensate are on the order of its healing length and its longitudinal dimension is much larger than its transverse ones [2]. The GP equation for the condensate wavefunction ψ(x, t) takes the formwhere |ψ| 2 is the number density, V (x) is the external trapping potential,] is proportional to the two-body scattering length, and ζ = |ψ| 2 |a| 3 is the dilute gas parameter [1,2].Experimentally realizable potentials V (x) include harmonic traps, quartic double-well traps, optical lattices and superlattices, and superpositions of lattices or superlattices with harmonic traps. The existence of quasi-1D ("cigar-shaped") BECs motivates the study of lower dimensional models such as Eq. (1). We focus here on the case of spatially periodic potentials without a confining trap along the dimension of the lattice, as that is of particular theoretical and experimental interest. For example, such potentials have been used to study Josephson effects [3], squeezed states [4], L...