ar X iv :c on d- m at /0 30 84 15 v4 [ co nd -m at .s of t] 2 6 N ov 2 00 3 Few-Body Systems 0, 1–6 (2013) Few- Body Systems c© by Springer-Verlag 2013 Printed in Austria Bright vortex solitons in Bose Condensates Sadhan K. Adhikari∗ Instituto de F́isica Teórica, Universidade Estadual Paulista, 01405-900 São Paulo, São Paulo, Brazil Abstract. We suggest the possibility of observing and studying bright vortex solitons in attractive Bose-Einstein condensates in three dimensions with a ra- dial trap. Such systems lie on the verge of critical stability and we discuss the conditions of their stability. We study the interaction between two such soli- tons. Unlike the text-book solitons in one dimension, the interaction between two radially trapped and axially free three-dimensional solitons is inelastic in nature and involves exchange of particles and deformation in shape. The interaction remains repulsive for all phase δ between them except for δ ≈ 0. 1 Introduction Solitary waves or solitons are a consequence of nonlinear dynamics. A classic text-book example of soliton appears in the following one-dimensional nonlinear free Schrödinger equation in dimensionless units [ −i ∂ ∂t − ∂2 ∂y2 − |Ψ(y, t)|2 ] Ψ(y, t) = 0. (1) The solitons of this equation are localized solution due to the attractive nonlinear interaction −|Ψ(y, t)|2 with wave function at time t and position y: Ψ(y, t) = √ 2|Ω| exp(−iΩt)sech(y √ |Ω|), with Ω the energy [1]. Solitons have been noted in optics [2], high-energy physics and water waves [3], and more recently in Bose-Einstein condensates (BEC) [4, 5]. The Schrödinger equation with a nonlinear interaction −|Ψ |2 does not sustain a lo- calized solitonic solution in three dimensions. However, a radially trapped and axially free version of this equation in three dimensions does sustain such a bright solitonic solution [6] which has been observed experimentally [4, 5]. Here we study the dynamics of these bright solitons. We also suggest that such soli- tons can be generated in an axially rotating nonzero angular momentum state, which are called bright vortex solitons and can be observed in BEC. ∗E-mail address: adhikari@ift.unesp.br http://arxiv.org/abs/cond-mat/0308415v4 2 Bright vortex solitons in Bose-Einstein Condensates A number of bright solitons constituting a soliton train was observed in an experiment by Strecker et al. [4], where they turned a repulsive BEC of 7Li atoms attractive by manipulating the background magnetic field near a Feshbach resonance [7]. It was found [4] that solitons in such a train usually stay apart. Also, often a soliton was found to be missing from a train [4]. There have been theoretical attempts [8, 9, 10] to simulate essentials of these experiments [4, 5]. We use the explicit numerical solution of the axially-symmetric mean-field Gross-Pitaevskii (GP) equation [11] to study the dynamics of bright solitons in a soliton train [4]. Attractive BEC’s may not form vortices in a thermodynam- ically stable state. However, due to the conservation of angular momentum, a vortex soliton train could be generated by suddenly changing the inter-atomic interaction in an axially-symmetric rotating vortex condensate from repulsive to attractive near a Feshbach resonance [7] in the same fashion as in the experi- ment by Strecker et al. [4] for a non-rotating BEC. Alternatively, a single vortex soliton could be prepared and studied in the laboratory by forming a vortex in a small repulsive condensate and then making the interaction attractive via a Feshbach resonance and subsequently reducing the axial trap slowly. 2 Mean-field Model and Results Mean-field Model: In a quantized vortex state [12], with each atom having an- gular momentum L~ along the axial y axis, the axially-symmetric wave function can be written as Ψ(r, τ) = ϕ(r, y, τ) exp(iLθ) where θ is the azimuthal angle and r the radial direction. The dynamics of the BEC in an axially-symmetric trap can be described by the following GP equation [11, 12] [ −i ∂ ∂t − ∂2 ∂r2 + 1 r ∂ ∂r − ∂2 ∂y2 + 1 4 ( r2 + λ2y2 ) + L2 − 1 r2 + 8 √ 2πn ∣ ∣ ∣ ∣ ϕ(r, y; t) r ∣ ∣ ∣ ∣ 2] ϕ(r, y; t) = 0, (2) where the length, time, and wave function are expressed in units of √ ~/(2mω), ω−1, and [r √ l3/ √ 8]−1, respectively. Here radial and axial trap frequencies are ω and λω, respectively, m is the atomic mass, l ≡ √ ~/(mω) is the harmonic oscillator length, and n = Na/l the nonlinearity with a the interatomic scat- tering length. For solitonic states n is negative. In terms of the one-dimensional probability P (y, t) defined by P (y, t) = 2π ∫ ∞ 0 dr|ϕ(r, y, t)|2/r, (3) the normalization of the wave function is given by ∫ ∞ −∞ dyP (y, t) = 1 We solve the GP equation (2) numerically using split-step time-iteration method using the Crank-Nicholson discretization scheme described recently [13]. The details of the numerical scheme for this problem can be found in [14]. Results: For bright solitons the nonlinearity n is negative. Under the con- ditions n < 0 and λ = 0, a soliton-type BEC state can be generated only for S. K. Adhikari 3 n = -0.2 L = 0 (a) 0 2 4 6 r -30 -15 0 15 30 y 0 0.05 0.1 0.15 |ϕ(r,y)/r| n = -1 L = 1 (b) 0 2 4 6 r -30 -15 0 15 30 y 0 0.05 0.1 0.15 |ϕ(r,y)/r| Figure 1. Three-dimensional wave function |ϕ(r, y)/r| vs. r and y for a single soliton with λ = 0, and (a) L = 0, n = −0.2, and (b) L = 1, n = −1. n greater than a critical value (ncr): ncr < n < 0. For n < ncr, the system becomes too attractive and collapses and no stable soliton could be generated. The actual value of ncr is a function of the trap parameter λ. For the spherically symmetric case λ = 1, and ncr = −0.575 [11, 12]. Numerically solving the GP equation (2) for λ = 0, we find that the critical n for collapse of a single soliton is ncr = −0.67 for L = 0 and ncr = −2.10 for L = 1. For L = 0 ncr = −0.67 is in close agreement with with ncr = −0.676 obtained by other workers [6, 15]. For L = 1, ncr = −2.10 should be contrasted with ncr = −2.20 obrained by Salasnich [16]. For L = 0, ω = 2π×800 Hz and final scattering length −3a0 as in the experiment of Strecker et al. [4], ncr = −0.67 corresponds to about 6000 7Li atoms. One can have proportionately about three times more atoms in the L = 1 state. A L = 0 soliton with n = −0.2 is illustrated in Fig. 1 (a) where we plot the three-dimensional wave function |ϕ(r, y)/r| vs. r and y. For L = 1 we calculated the soliton for n = −1 and plot |ϕ(r, y)/r| in Fig. 1 (b). Because of the radial trap the soliton remains confined in the radial direction r, although free to move in the axial y direction. The nature of the two wave functions are different. For L = 0, the condensate has maximum density for r = 0. For L = 1, because of rotation a vortex has been generated along the axial direction corresponding to a zero density for r = 0. Next we consider the interaction between two solitons with a phase difference δ given by the following superposition of two solitons ϕ̄ at ±y0 at time t = 0 ϕ(r, y) = |ϕ̄(r, y +y0)|+eiδ |ϕ̄(r, y−y0)|. The time evolution of these two solitons is found using the solution of (2) for different δ. In the present simulation we consider two L = 1 equal vortex solitons each of n = −0.4 for λ = 0 initially at positions y0 = ±15 and observe them for an interval of time t = 400. We also 4 Bright vortex solitons in Bose-Einstein Condensates t = 0 (a) 2 4 6r -80 -40 0 40 80 y 0 0.04 0.08 |ϕ(r,y)/r| t = 100 (b) 2 4 6r -80 -40 0 40 80 y 0 0.04 0.08 |ϕ(r,y)/r| t = 200 (c) 2 4 6r -80 -40 0 40 80 y 0 0.04 0.08 |ϕ(r,y)/r| t = 300 (d) 2 4 6r -80 -40 0 40 80 y 0 0.04 0.08 |ϕ(r,y)/r| Figure 2. The wave function |ϕ(r, y)/r| vs. r and y of two L = 0 solitons with n = −0.2 each and δ = π/2 at times (a) t = 0, (b) 100, (c) 200 and (d) 300. consider the evolution of two L = 0 solitons each of n = −0.2 for y0 = ±15. The solitons interact by exchanging particles and after an interval of time two unequal solitons are generated from two equal solitons. The evolution of the two solitons in the L = 0 case for δ = π/2 is shown in Fig. 2. Similar evolution for the L = 1 case is reported in [14]. From Fig. 2 we find that at t = 0 the two solitons are equal and symmetri- cally located. However, this symmetry is broken for t > 0. The asymmetry and separation in the final position of the solitons y1 and y2 are best studied via |(|y1| − |y2|)| and |y1 − y2|, respectively, at time t = 400 for different phase δ between the solitons and in Fig. 3 (a) we plot the same for different δ for the L = 0 case. We find that the asymmetry is zero for δ = π and 0 and is largest for a δ in between. However the separation increases monotonically as δ increases from 0 to π. Hence the interaction is repulsive for almost all δ except for δ ≈ 0. Closely associated with the asymmetry and separation in the final position of the solitons is the number of exchanged atoms between the two solitons, which demonstrates the change in the sizes of the solitons. Actually, the smaller soliton travels faster and the larger one travels slower. This results in the asymmetry in S. K. Adhikari 5 0 20 40 60 80 0 0.2 0.4 0.6 0.8 1 0 20 40 60 80 s e p a ra ti o n | y 1 -y 2 | a s y m m e tr y | (| y 1 |- |y 2 |) | δ (π) (a) L = 0 asymmetry separation 0.3 0.4 0.5 0 100 200 300 400 500 600 700 N R / N t (b) L = 0 δ = 0, π δ = 3π/4 δ = π/8 δ = π/2 δ = π/4 0 20 40 60 80 0 0.2 0.4 0.6 0.8 1 0 20 40 60 80 s e p a ra ti o n | y 1 -y 2 | a s y m m e tr y | (| y 1 |- |y 2 |) | δ (π) (c) L = 1 asymmetry separation 0.3 0.4 0.5 0 100 200 300 400 500 600 700 N R / N t (d) L = 1 δ = 0, π δ = 3π/4 δ = π/2 δ = π/4 δ = π/8 Figure 3. (a) Final asymmetry and separation in the position of the two solitons of L = 0 and n = −0.2 each at t = 400 vs. phase δ. (b) Ratio NR/N vs. t for different phase δ for two L = 0 and n = −0.2 solitons initially at ±15. NR is the number of atoms in the right soliton and N the total number of atoms. (c) Same as (a) for two solitons of L = 1 and n = −0.4. (d) Same as (b) for two solitons of L = 1 and n = −0.4. the final positions. The change in the sizes of the solitons for L = 0 is demon- strated in the plot of NR/N vs. t for different δ in Fig. 3 (b), where NR is the number of atoms in the right soliton and N the total number of atoms in the two solitons. In Figs. 3 (c) and (d) we plot the same for L = 1. The variation of NR/N is qualitatively similar for L = 0 and 1 in Figs. 3 and also to that found in [10] for L = 0. However, there are quantitative differences, specially at large times. In the present simulation we find that, for the change δ → −δ, NR/N → NL/N , where NL ≡ (N − NR) is the number of atoms in the left soliton. If the phase difference δ between two neighboring solitons is not close to zero, they experience overall repulsion and stay apart. However, for δ close to zero they interact attractively and often a soliton could be lost as observed in the experiment of Strecker et al. [4]. Throughout this investigation in the interaction of two equal solitons we assumed that the nonlinearity |n| for each is less than |ncr|/2, so that a stable solitonic condensate with total |n| < |ncr| exists when the two coalesce. However, if two solitons each with |n| > |ncr|/2 encounter for δ = 0, the system is expected to coalesce, collapse and emit atoms via three- body recombination. It is possible that in this case only a smaller single soliton survives. This might also explain some missing soliton(s) in experiment. 6 Bright vortex solitons in Bose-Einstein Condensates 3 Conclusion We emphasize the possibility of creating and studying bright vortex solitons of an attractive BEC in laboratory under radial trapping. We determine the condition of critical stability of bright solitons. Employing a numerical solution of the GP equation with axial symmetry, we have performed a realistic mean-field study of the interaction among two bright solitons in a train and find the overall interaction to be repulsive except for phase δ between neighbors close to 0. There is an inelastic exchange of atoms between two solitons resulting in a change of size and shape. Except in the δ ≈ 0 case, the solitons in a train stay apart and never cross each other as observed in the experiment by Strecker et al. [4]. For δ ≈ 0 a single soliton can often disappear as a result of the attractive interaction among solitons, as observed experimentally by Strecker et al.. The L = 1 vortex solitons can accommodate a larger number of atoms and the present study may motivate future experiments with them. Acknowledgement. Work partially supported by the CNPq of Brazil. References 1. Kivshar Y. S., Alexander T. J. et al: Phys. Lett. A 278, 225 (2001) 2. Hasegawa A. and Tappert F.: Appl. Phys. Lett. 23, 171 (1973); 23, 142 (1973); Zakharov V. E. and Shabat A. B.: Sov. Phys. JETP 34, 62 (1972) 3. Taylor J. R.: Editor Optical Solitons - Theory and Experiment (Cambridge, Cambridge, 1992) 4. Strecker K. E., Partridge G. B. et al: Nature 417, 150 (2002) 5. Khaykovich L., Schreck F., Ferrari G. et al: Science 296, 1290 (2002) 6. Pérez-Garćia V. M., Michinel H. et al: Phys. Rev. A 57, 3837 (1998) 7. Inouye S., Andrews M. R., Stenger J. et al: Nature 392, 151 (1998) 8. Al Khawaja U., Stoof H. T. C. et al: Phys. Rev. Lett. 89, 200404 (2002) 9. Salasnich L., Parola A. and Reatto L.: Phys. Rev. A 66, 043603 (2002) 10. Salasnich L., Parola A. and Reatto L.: Phys. Rev. Lett. 91 080405 (2003) 11. Dalfovo F., Giorgini S., Pitaevskii L. P. et al: Rev. Mod. Phys. 71, 463 (1999) 12. Adhikari S. K.: Phys. Rev. E 65, 016703 (2002) 13. Adhikari S. K. and Muruganandam P.: J. Phys. B 35, 2831 (2002); Muru- ganandam P. and Adhikari S. K.: J. Phys. B 36, 2501 (2003) 14. Adhikari S. K.: New J. Phys 5, 137 (2003) 15. Gammal A., Tomio L. and Frederico T.: Phys. Rev. A 66, 043619 (2002) 16. Salasnich L.: e-print cond-mat/0309071 (2003) http://arxiv.org/abs/cond-mat/0309071 Introduction Mean-field Model and Results Conclusion