We present numerical solutions in the form of spatially localized nonrotating and rotating (azimuthon) multisolitons in the two-dimensional (``pancake-shaped configuration'') Bose-Einstein condensate (BEC) confined by a parabolic trap. The numerical relaxation procedure with a stabilizing factor is used.
We present exact numerical solutions in the form of spatially localized three-dimensional (3D) nonrotating and rotating (azimuthon) multipole solitons in the Bose-Einstein condensate (BEC) confined by a parabolic trap. We numerically show that the 3D azimuthon solutions exist as a continuous family parametrized by the angular velocity (or, equivalently, the modulational depth). By a linear stability analysis we show that 3D azimuthons with a sufficiently large phase modulational depth can be stable. The results are confirmed by direct numerical simulations of the Gross-Pitaevskii equation.
We have performed numerical analysis of the two-dimensional (2D) soliton solutions in Bose-Einstein condensates with nonlocal dipole-dipole interactions. For the modified 2D Gross-Pitaevski equation with nonlocal and attractive local terms, we have found numerically different types of nonlinear localized structures such as fundamental solitons, radially symmetric vortices, nonrotating multisolitons (dipoles and quadrupoles), and rotating multisolitons (azimuthons). By direct numerical simulations we show that these structures can be made stable.
We introduce matter-wave vector azimuthons, i.e., spatially localized vortex states with azimuthal modulations of density, in multicomponent Bose-Einstein condensates. These localized states generalize spatially modulated vortex solitons introduced earlier in nonlinear optics [A. S. Desyatnikov, A. A. Sukhorukov, and Yu. S. Kivshar, Phys. Rev. Lett. 95, 203904 (2005)] and Bose-Einstein condensates [V. M. Lashkin, Phys. Rev. A 77, 025602 (2008)]. We find, numerically, nonrotating and rotating two-component azimuthons in a Bose-Einstein condensate with a negative scattering length confined by a quasi-two-dimensional parabolic trap.
We numerically obtain stationary rings of vortices in pancake-shaped Bose-Einstein condensates confined in a three-dimensional nonrotating trap. For this purpose we use a static axisymmetric trapping potential that can sustain locally stable off-axis vortices. We analyze the maximum number of vortices the system can host as a function of the number of particles. We also show that this system provides a very suitable scenario for predicting vortex dynamics in inhomogeneous media. Specifically, on the one hand, we derive an exact and simple analytical expression for the velocity field in a particular vortex position due to the presence of the other vortices of the array. On the other hand, using the fact that in stationary conditions this field should balance all other contributions to the velocity, we investigate the applicability of approximated expressions for the precession velocity obtained in previous works for condensates with a single vortex.
We perform accurate investigation of stability of localized vortices in an effectively two-dimensional (``pancake-shaped'') trapped Bose-Einstein condensate with negative scattering length. The analysis combines computation of the stability eigenvalues and direct simulations. The states with vorticity S=1 are stable in a third of their existence region, 0<N<(1∕3)Nmax(S=1), where N is the number of atoms, and Nmax(S=1) is the corresponding collapse threshold. Stable vortices easily self-trap from arbitrary initial configurations with embedded vorticity. In an adjacent interval, (1∕3)Nmax(S=1)<N<0.43Nmax(S=1), the unstable vortex periodically splits in twofragments and recombines. At N>0.43Nmax(S=1), the fragments do not recombine, as each one collapses by itself. The results are compared with those in the full three-dimensional (3D) Gross-Pitaevskii equation. In a moderately anisotropic 3D configuration, with the aspect ratio 10, the stability interval of the S=1 vortices occupies $\ensuremath{\approx}40%$ of their existence region, hence the two-dimensional (2D) limit provides for a reasonable approximation in this case. For the isotropic 3D configuration, the stability interval expands to 65% of the existence domain. Overall, the vorticity heightens the actual collapse threshold by a factor of up to 2. All vortices with $S\ensuremath{\geqslant}2$ are unstable.