We introduce a dynamical model of a Bose-Einstein condensate based on the two-dimensional Gross-Pitaevskii equation, in which the nonlinear coefficient is a function of radius. The model describes a situation with spatial modulation of the negative atomic scattering length, via the Feshbach resonance controlled by a properly shaped magnetic of optical field. We focus on the configuration with the nonlinear coefficient different from zero in a circle or annulus, including the case of a narrow ring. Two-dimensional axisymmetric solitons are found in a numerical form, and also by means of a variational approximation; for an infinitely narrow ring, the soliton is found in an exact form (in the latter case, exact solitons are also found in a two-component model). A stability region for the solitons is identified by means of numerical and analytical methods. In particular, if the nonlinearity is supported on the annulus, the upper stability border is determined by azimuthal perturbations; the stability region disappears if the ratio of the inner and outer radii of the annulus exceeds a critical value . The model gives rise to bistability, as the stationary solitons coexist with stable axisymmetric breathers, whose stability region extends to higher values of the norm than that of the static solitons. The collapse threshold strongly increases with the radius of the inner hole of the annulus. Vortex solitons are found too, but they are unstable.
We introduce a dynamical model of a Bose-Einstein condensate based on the two-dimensional Gross-Pitaevskii equation, in which the nonlinear coefficient is a function of radius. The model describes a situation with spatial modulation of the negative atomic scattering length, via the Feshbach resonance controlled by a properly shaped magnetic of optical field. We focus on the configuration with the nonlinear coefficient different from zero in a circle or annulus, including the case of a narrow ring. Two-dimensional axisymmetric solitons are found in a numerical form, and also by means of a variational approximation; for an infinitely narrow ring, the soliton is found in an exact form (in the latter case, exact solitons are also found in a two-component model). A stability region for the solitons is identified by means of numerical and analytical methods. In particular, if the nonlinearity is supported on the annulus, the upper stability border is determined by azimuthal perturbations; the stability region disappears if the ratio of the inner and outer radii of the annulus exceeds a critical value . The model gives rise to bistability, as the stationary solitons coexist with stable axisymmetric breathers, whose stability region extends to higher values of the norm than that of the static solitons. The collapse threshold strongly increases with the radius of the inner hole of the annulus. Vortex solitons are found too, but they are unstable.
We introduce a dynamical model of a Bose-Einstein condensate based on the one-dimensional (1D) Gross-Pitaevskii equation (GPE) with a nonlinear optical lattice (NOL), which is represented by the cubic term whose coefficient is periodically modulated in the coordinate. The model describes a situation when the atomic scattering length is spatially modulated, via the optically controlled Feshbach resonance, in an optical lattice created by interference of two laser beams. Relatively narrow solitons supported by the NOL are predicted by means of the variational approximation (VA), and an averaging method is applied to broad solitons. A different feature is a minimum norm (number of atoms), N = N(min), necessary for the existence of solitons. The VA predicts N(min) very accurately. Numerical results are chiefly presented for the NOL with the zero spatial average value of the nonlinearity coefficient. Solitons with values of the amplitude A larger than at N = N(min) are stable. Unstable solitons with smaller, but not too small, A rearrange themselves into persistent breathers. For still smaller A, the soliton slowly decays into radiation without forming a breather. Broad solitons with very small A are practically stable, as their decay is extremely slow. These broad solitons may freely move across the lattice, featuring quasielastic collisions. Narrow solitons, which are strongly pinned to the NOL, can easily form stable complexes. Finally, the weakly unstable low-amplitude solitons are stabilized if a cubic term with a constant coefficient, corresponding to weak attraction, is included in the GPE.
Soliton is an exotic topological excitation, and it widely exists in various nonlinear systems, such as nonlinear optics, Bose-Einstein condensates, classical and quantum fluids, plasma, magnetic materials, etc. A stable soliton can propagate with constant amplitude and velocity, and maintain its shape. Two-dimensional and three-dimensional solitons are usually hard to stabilize, and how to realize stable two-dimensional or three-dimensional solitons has aroused the great interest of the researchers. Ring dark soliton is a kind of two-dimensional soliton, which was first theoretically predicted and experimentally realized in nonlinear optical systems. Compared with the usual two-dimensional solitons, ring dark solitons have good stability and rich dynamical behaviors. Owing to their highly controllable capability, Bose-Einstein condensates provide a new platform for studying the ring dark solitons. Based on the recent progress in Bose-Einstein condensates and solitons, this paper reviews the research on the analytic solutions, stability, as well as the decay dynamics of ring dark solitons in Bose-Einstein condensates. A transform method is introduced, which generalizes the analytic solutions of ring dark solitons from a homogeneous system with time-independent nonlinearity to a harmonically trapped inhomogeneous system with time-dependent nonlinearity. The stability phase diagram of the ring dark soliton under deformation perturbations is discussed by numerically solving the Gross-Pitaevskii equations in the mean-field theory. A method of enhancing the stability of ring dark solitons by periodically modulating the nonlinear coefficients is introduced. It is also shown that the periodically modulated nonlinear coefficient can be experimentally realized by the Feshbach resonance technology. In addition, we discuss the dynamics of the decay of ring dark solitons. It is found that the ring dark soliton can decay into various vortex clusters composed of vortices and antivortices. This opens a new avenue to the investigation of vortex dynamics and quantum turbulence. It is also found that the ring dark solitons combined with periodic modulated nonlinearity can give rise to the pattern formation, which is an interesting nonlinear phenomenon widely explored in all the fields of nature. Finally, some possible research subjects about ring dark solitons in future research are also discussed.
We consider, by means of the variational approximation (VA) and direct numerical simulations of the Gross-Pitaevskii (GP) equation, the dynamics of two-dimensional (2D) and 3D condensates with a scattering length containing constant and harmonically varying parts, which can be achieved with an ac magnetic field tuned to the Feshbach resonance. For a rapid time modulation, we develop an approach based on the direct averaging of the GP equation, without using the VA. In the 2D case, both VA and direct simulations, as well as the averaging method, reveal the existence of stable self-confined condensates without an external trap, in agreement with qualitatively similar results recently reported for spatial solitons in nonlinear optics. In the 3D case, the VA again predicts the existence of a stable self-confined condensate without a trap. In this case, direct simulations demonstrate that the stability is limited in time, eventually switching into collapse, even though the constant part of the scattering length is positive (but not too large). Thus a spatially uniform ac magnetic field, resonantly tuned to control the scattering length, may play the role of an effective trap confining the condensate, and sometimes causing its collapse.
Using variational and numerical solutions of the mean-field Gross-Pitaevskii equation we show that a bright soliton can be stabilized in a trapless three-dimensional attractive Bose-Einstein condensate (BEC) by a rapid periodic temporal modulation of scattering length alone by using a Feshbach resonance. This scheme also stabilizes a rotating vortex soliton in two dimensions. Apart from possible experimental application in BEC, the present study suggests that the spatiotemporal solitons of nonlinear optics in three dimensions can also be stabilized in a layered Kerr medium with sign-changing nonlinearity along the propagation direction.
We review our recent theoretical advances in the dynamics of Bose-Einstein condensates with tunable interactions using Feshbach resonance and external potential. A set of analytic and numerical methods for Gross-Pitaevskii equations are developed to study the nonlinear dynamics of Bose-Einstein condensates. Analytically, we present the integrable conditions for the Gross-Pitaevskii equations with tunable interactions and external potential, and obtain a family of exact analytical solutions for one- and two-component Bose-Einstein condensates in one and two-dimensional cases. Then we apply these models to investigate the dynamics of solitons and collisions between two solitons. Numerically, the stability of the analytic exact solutions are checked and the phenomena, such as the dynamics and modulation of the ring dark soliton and vector-soliton, soliton conversion via Feshbach resonance, quantized soliton and vortex in quasi-two-dimensional are also investigated. Both the exact and numerical solutions show that the dynamics of Bose-Einstein condensates can be effectively controlled by the Feshbach resonance and external potential, which offer a good opportunity for manipulation of atomic matter waves and nonlinear excitations in Bose-Einstein condensates.
We consider a model of Bose-Einstein condensates which combines a stationary optical lattice (OL) and periodic change of the sign of the scattering length (SL) due to the Feshbach-resonance management. Ordinary solitons and ones of the gap type being possible, respectively, in the model with constant negative and positive SL, an issue of interest is to find solitons alternating, in the case of the low-frequency modulation, between shapes of both types, across the zero-SL point. We find such \textit{alternate solitons} and identify their stability regions in the 2D and 1D models. Three types of the dynamical regimes are distinguished: stable, unstable, and semi-stable. In the latter case, the soliton sheds off a conspicuous part of its initial norm before relaxing to a stable regime. In the 2D case, the threshold (minimum number of atoms) necessary for the existence of the alternate solitons is essentially higher than its counterparts for the ordinary and gap solitons in the static model. In the 1D model, the alternate solitons are also found only above a certain threshold, while the static 1D models have no threshold. In the 1D case, stable antisymmetric alternate solitons are found too. Additionally, we consider a possibility to apply the variational approximation (VA) to the description of stationary gap solitons, in the case of constant positive SL. It predicts the solitons in the first finite bandgap very accurately, and does it reasonably well in the second gap too. In higher bands, the VA predicts a border between tightly and loosely bound solitons.