We derive an effective nonpolynomial Schr\"odinger equation (NPSE) for self-repulsive or attractive BEC in the nearly one-dimensional cigar-shaped trap, with the transverse confining frequency periodically modulated along the axial direction. In addition to the usual linear cigar-shaped trap, where the periodic modulation emulates the action of an optical lattice (OL), the model may be also relevant to toroidal traps, where an ordinary OL cannot be created. For either sign of the nonlinearity, extended and localized states are found, in the numerical form [using both the effective NPSE and the full three-dimensional (3D) Gross-Pitaevskii equation] and by means of the variational approximation (VA). The latter is applied to construct ground-state solitons and predict the collapse threshold in the case of self-attraction. It is shown that numerical solutions provided by the one-dimensional NPSE are always very close to full 3D solutions, and the VA yields quite reasonable results too. The transition from delocalized states to gap solitons, in the first finite bandgap of the linear spectrum, is examined in detail, for the repulsive and attractive nonlinearities alike.
We derive an effective nonpolynomial Schr"odinger equation (NPSE) for self-repulsive or attractive BEC in the nearly one-dimensional cigar-shaped trap, with the transverse confining frequency periodically modulated along the axial direction. In addition to the usual linear cigar-shaped trap, where the periodic modulation emulates the action of an optical lattice (OL), the model may be also relevant to toroidal traps, where an ordinary OL cannot be created. For either sign of the nonlinearity, extended and localized states are found, in the numerical form [using both the effective NPSE and the full three-dimensional (3D) Gross-Pitaevskii equation] and by means of the variational approximation (VA). The latter is applied to construct ground-state solitons and predict the collapse threshold in the case of self-attraction. It is shown that numerical solutions provided by the one-dimensional NPSE are always very close to full 3D solutions, and the VA yields quite reasonable results too. The transition from delocalized states to gap solitons, in the first finite bandgap of the linear spectrum, is examined in detail, for the repulsive and attractive nonlinearities alike.
We study vortical states in a Bose-Einstein condensate (BEC) filling a cigar-shaped trap. An effective one-dimensional (1D) nonpolynomial Schr\"odinger equation (NPSE) is derived in this setting, for the models with both repulsive and attractive interatomic interactions. Analytical formulas for the density profiles are obtained from the NPSE in the case of self-repulsion within the Thomas-Fermi approximation, and in the case of the self-attraction as exact solutions (bright solitons). A crucially important ingredient of the analysis is the comparison of these predictions with direct numerical solutions for the vortex states in the underlying 3D Gross-Pitaevskii equation. The comparison demonstrates that the NPSE provides for a very accurate approximation, in all the cases, including the prediction of the stability of the bright solitons and collapse threshold for them. In addition to the straight cigar-shaped trap, we also consider a torus-shaped configuration. In that case, we find a threshold for the transition from the axially uniform state, with the transverse intrinsic vorticity, to a symmetry-breaking pattern, due to the instability in the self-attractive BEC filling the circular trap.
Abstract Starting from the standard three-dimensional (3D) Gross–Pitaevskii equation (GPE) and using a variational approximation, we derive an effective one-dimensional nonpolynomial Schrödinger equation (1D-NPSE) governing the axial dynamics of atomic Bose–Einstein condensates (BECs) under the action of a singular but physically relevant funnel-shaped transverse trap, i.e. an attractive 2D potential ∼−1/ r (where r is the radial coordinate in the transverse plane), in combination with the repulsive self-interaction. Wave functions of the trapped BEC are regular, in spite of the potential’s singularity. The model applies to a condensate of particles (small molecules) carrying a permanent electric dipole moment in the field of a uniformly charged axial thread, as well as to a quantum gas of magnetic atoms pulled by an axial electric current. By means of numerical simulations, we verify that the effective 1D-NPSE provides accurate static and dynamical results, in comparison to the full 3D GPE, for both repulsive and attractive signs of the intrinsic nonlinearity.
Within the framework of a mean-field description, we investigate atomic Bose-Einstein condensates, with attraction between atoms, under the action of a strong transverse confinement and periodic [optical-lattice (OL)] axial potential. Using a combination of the variational approximation, one-dimensional (1D) nonpolynomial Schr\"odinger equation, and direct numerical solutions of the underlying 3D Gross-Pitaevskii equation, we show that the ground state of the condensate is a soliton belonging to the semi-infinite band gap of the periodic potential. The soliton may be confined to a single cell of the lattice or extended to several cells, depending on the effective self-attraction strength g (which is proportional to the number of atoms bound in the soliton) and depth of the potential, V0, the increase of V0 leading to strong compression of the soliton. We demonstrate that the OL is an effective tool to control the soliton's shape. It is found that, due to the 3D character of the underlying setting, the ground-state soliton collapses at a critical value of the strength, g=gc, which gradually decreases with the increase of V0; under typical experimental conditions, the corresponding maximum number of 7Li atoms in the soliton, Nmax, ranges between 8000 and 4000. Examples of stable multipeaked solitons are also found in the first finite band gap of the lattice spectrum. The respective critical value gc again slowly decreases with the increase of V0, corresponding to ${N}_{\mathrm{max}}\ensuremath{\simeq}5000$.
Starting from the three-dimensional (3D) Gross-Pitaevskii equation and using a variational approach, we derive an effective 1D wave equation that describes the axial dynamics of a Bose condensate confined in an external potential with cylindrical symmetry. The trapping potential is harmonic in the transverse direction and generic in the axial one. Our equation, that is a time-dependent nonpolynomial nonlinear Schr\"odinger equation (1D NPSE), can be used to model cigar-shaped condensates, whose dynamics is essentially 1D. We show that 1D NPSE gives much more accurate results than all other effective equations recently proposed. By using 1D NPSE we find analytical solutions for bright and dark solitons, which generalize the ones known in the literature. We deduce also an effective 2D nonpolynomial Schr\"odinger equation (2D NPSE) that models disk-shaped Bose condensates confined in an external trap that is harmonic along the axial direction and generic in the transverse direction. In the limiting cases of weak and strong interaction, our approach gives rise to Schr\"odinger-like equations with different polynomial nonlinearities.
The Bose-Einstein condensate (BEC), confined in a combination of the cigar-shaped trap and axial optical lattice, is studied in the framework of two models described by two versions of the one-dimensional (1D) discrete nonpolynomial Schrodinger equation (NPSE). Both models are derived from the three-dimensional Gross-Pitaevskii equation (3D GPE). To produce "model 1" (which was derived in recent works), the 3D GPE is first reduced to the 1D continual NPSE, which is subsequently discretized. "Model 2," which was not considered before, is derived by first discretizing the 3D GPE, which is followed by the reduction in the dimension. The two models seem very different; in particular, model 1 is represented by a single discrete equation for the 1D wave function, while model 2 includes an additional equation for the transverse width. Nevertheless, numerical analyses show similar behaviors of fundamental unstaggered solitons in both systems, as concerns their existence region and stability limits. Both models admit the collapse of the localized modes, reproducing the fundamental property of the self-attractive BEC confined in tight traps. Thus, we conclude that the fundamental properties of discrete solitons predicted for the strongly trapped self-attracting BEC are reliable, as the two distinct models produce them in a nearly identical form. However, a difference between the models is found too, as strongly pinned (very narrow) discrete solitons, which were previously found in model 1, are not generated by model 2-in fact, in agreement with the continual 1D NPSE, which does not have such solutions either. In that respect, the newly derived model provides for a more accurate approximation for the trapped BEC.
We report results of a systematic analysis of matter-wave gap solitons (GSs) in three-dimensional self-repulsive Bose-Einstein condensates (BECs) loaded into a combination of a cigar-shaped trap and axial optical-lattice (OL) potential. Basic cases of the strong, intermediate, and weak radial (transverse) confinement are considered, as well as settings with shallow and deep OL potentials. Only in the case of the shallow lattice combined with tight radial confinement, which actually has little relevance to realistic experimental conditions, does the usual one-dimensional (1D) cubic Gross-Pitaevskii equation (GPE) furnish a sufficiently accurate description of GSs. However, the effective 1D equation with the nonpolynomial nonlinearity, derived in Ref. [Phys. Rev. A 77, 013617 (2008)], provides for quite an accurate approximation for the GSs in all cases, including the situation with weak transverse confinement, when the soliton's shape includes a considerable contribution from higher-order transverse modes, in addition to the usual ground-state wave function of the respective harmonic oscillator. Both fundamental GSs and their multipeak bound states are considered. The stability is analyzed by means of systematic simulations. It is concluded that almost all the fundamental GSs are stable, while their bound states may be stable if the underlying OL potential is deep enough.