Deviation from one-dimensionality in stationary properties and collisional dynamics of matter-wave solitons

By means of analytical and numerical methods, we study how the residual three-dimensionality affects dynamics of solitons in an attractive Bose-Einstein condensate loaded into a cigar-shaped trap. Based on an effective 1D Gross-Pitaevskii equation that includes an additional quintic self-focusing term, generated by the tight transverse confinement, we find a family of exact one-soliton solutions and demonstrate stability of the entire family, despite the possibility of collapse in the 1D equation with the quintic self-focusing nonlinearity. Simulating collisions between two solitons in the same setting, we find a critical velocity, $V_{c}$, below which merger of identical in-phase solitons is observed. Dependence of $V_{c}$ on the strength of the transverse confinement and number of atoms in the solitons is predicted by means of the perturbation theory and investigated in direct simulations. Symmetry breaking in collisions of identical solitons with a nonzero phase difference is also shown in simulations and qualitatively explained by means of an analytical approximation.Comment: 10 pages, 7 figure

Symbiotic gap and semi-gap solitons in Bose-Einstein condensates

Using the variational approximation and numerical simulations, we study one-dimensional gap solitons in a binary Bose-Einstein condensate trapped in an optical-lattice potential. We consider the case of inter-species repulsion, while the intra-species interaction may be either repulsive or attractive. Several types of gap solitons are found: symmetric or asymmetric; unsplit or split, if centers of the components coincide or separate; intra-gap (with both chemical potentials falling into a single bandgap) or inter-gap, otherwise. In the case of the intra-species attraction, a smooth transition takes place between solitons in the semi-infinite gap, the ones in the first finite bandgap, and semi-gap solitons (with one component in a bandgap and the other in the semi-infinite gap).Comment: 5 pages, 9 figure

Matter-wave 2D solitons in crossed linear and nonlinear optical lattices

It is demonstrated the existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with linear OL in the $x-$direction and nonlinear OL (NOL) in the $y-$direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance. In particular, we show that such crossed linear and nonlinear OL allows to stabilize two-dimensional (2D) solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach (VA), with the Vakhitov-Kolokolov (VK) necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation (GPE). Very good agreement of the results corresponding to both treatments is observed.Comment: 8 pages (two-column format), with 16 eps-files of 4 figure

Universal fractal structures in the weak interaction of solitary waves in generalized nonlinear Schr\"{o}dinger equations

Weak interactions of solitary waves in the generalized nonlinear Schr\"{o}dinger equations are studied. It is first shown that these interactions exhibit similar fractal dependence on initial conditions for different nonlinearities. Then by using the Karpman-Solov'ev method, a universal system of dynamical equations is derived for the velocities, amplitudes, positions and phases of interacting solitary waves. These dynamical equations contain a single parameter, which accounts for the different forms of nonlinearity. When this parameter is zero, these dynamical equations are integrable, and the exact analytical solutions are derived. When this parameter is non-zero, the dynamical equations exhibit fractal structures which match those in the original wave equations both qualitatively and quantitatively. Thus the universal nature of fractal structures in the weak interaction of solitary waves is analytically established. The origin of these fractal structures is also explored. It is shown that these structures bifurcate from the initial conditions where the solutions of the integrable dynamical equations develop finite-time singularities. Based on this observation, an analytical criterion for the existence and locations of fractal structures is obtained. Lastly, these analytical results are applied to the generalized nonlinear Schr\"{o}dinger equations with various nonlinearities such as the saturable nonlinearity, and predictions on their weak interactions of solitary waves are made.Comment: 22pages, 15 figure

Matter-wave vortices in cigar-shaped and toroidal waveguides

We study vortical states in a Bose-Einstein condensate (BEC) filling a cigar-shaped trap. An effective one-dimensional (1D) nonpolynomial Schroedinger equation (NPSE) is derived in this setting, for the models with both repulsive and attractive inter-atomic 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 (GPE). 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.Comment: 6 pages, Physical Review A, in pres

Spatiotemporal discrete multicolor solitons

We have found various families of two-dimensional spatiotemporal solitons in quadratically nonlinear waveguide arrays. The families of unstaggered odd, even and twisted stationary solutions are thoroughly characterized and their stability against perturbations is investigated. We show that the twisted and even solutions display instability, while most of the odd solitons show remarkable stability upon evolution.Comment: 18 pages,7 figures. To appear in Physical Review

Vector solitons in nearly-one-dimensional Bose-Einstein condensates

We derive a system of nonpolynomial Schroedinger equations (NPSEs) for one-dimensional wave functions of two components in a binary self-attractive Bose-Einstein condensate loaded in a cigar-shaped trap. The system is obtained by means of the variational approximation, starting from the coupled 3D Gross-Pitaevskii equations and assuming, as usual, the factorization of 3D wave functions. The system can be obtained in a tractable form under a natural condition of symmetry between the two species. A family of vector (two-component) soliton solutions is constructed. Collisions between orthogonal solitons (ones belonging to the different components) are investigated by means of simulations. The collisions are essentially inelastic. They result in strong excitation of intrinsic vibrations in the solitons, and create a small orthogonal component ("shadow") in each colliding soliton. The collision may initiate collapse, which depends on the mass and velocities of the solitons.Comment: 7 pages, 6 figures; Physical Review A, in pres

Two-dimensional solitons with hidden and explicit vorticity in bimodal cubic-quintic media

We demonstrate that two-dimensional two-component bright solitons of an annular shape, carrying vorticities $(m,\pm m)$ in the components, may be stable in media with the cubic-quintic nonlinearity, including the \textit{hidden-vorticity} (HV) solitons of the type $(m,-m)$, whose net vorticity is zero. Stability regions for the vortices of both $(m,\pm m)$ types are identified for $m=1$, 2, and 3, by dint of the calculation of stability eigenvalues, and in direct simulations. A novel feature found in the study of the HV solitons is that their stability intervals never reach the (cutoff) point at which the bright vortex carries over into a dark one, hence dark HV solitons can never be stable, contrarily to the bright ones. In addition to the well-known symmetry-breaking (\textit{external}) instability, which splits the ring soliton into a set of fragments flying away in tangential directions, we report two new scenarios of the development of weak instabilities specific to the HV solitons. One features \textit{charge flipping}, with the two components exchanging the angular momentum and periodically reversing the sign of their spins. The composite soliton does not split in this case, therefore we identify such instability as an \textit{intrinsic} one. Eventually, the soliton splits, as weak radiation loss drives it across the border of the ordinary strong (external) instability. Another scenario proceeds through separation of the vortex cores in the two components, each individual core moving toward the outer edge of the annular soliton. After expulsion of the cores, there remains a zero-vorticity breather with persistent internal vibrations.Comment: 10 pages, 11 figure

Discrete Solitons and Vortices on Anisotropic Lattices

We consider effects of anisotropy on solitons of various types in two-dimensional nonlinear lattices, using the discrete nonlinear Schr{\"{o}}dinger equation as a paradigm model. For fundamental solitons, we develop a variational approximation, which predicts that broad quasi-continuum solitons are unstable, while their strongly anisotropic counterparts are stable. By means of numerical methods, it is found that, in the general case, the fundamental solitons and simplest on-site-centered vortex solitons ("vortex crosses") feature enhanced or reduced stability areas, depending on the strength of the anisotropy. More surprising is the effect of anisotropy on the so-called "super-symmetric" intersite-centered vortices ("vortex squares"), with the topological charge $S$ equal to the square's size $M$: we predict in an analytical form by means of the Lyapunov-Schmidt theory, and confirm by numerical results, that arbitrarily weak anisotropy results in dramatic changes in the stability and dynamics in comparison with the \emph{degenerate}, in this case, isotropic limit.Comment: 10 pages + 7 figure

Stability of vortex solitons in a photorefractive optical lattice

Stability of off-site vortex solitons in a photorefractive optical lattice is analyzed. It is shown that such solitons are linearly unstable in both the high and low intensity limits. In the high-intensity limit, the vortex looks like a familiar ring vortex, and it suffers oscillatory instabilities. In the low-intensity limit, the vortex suffers both oscillatory and Vakhitov-Kolokolov instabilities. However, in the moderate-intensity regime, the vortex becomes stable if the lattice intensity or the applied voltage is above a certain threshold value. Stability regions of vortices are also determined at typical experimental parameters.Comment: 3 pages, 5 figure