Research Article  Open Access
Dynamics of BoseEinstein Condensates: Exact Representation and Topological Classification of Coherent Matter Waves
Abstract
By using the bifurcation theory of dynamical systems, we present the exact representation and topological classification of coherent matter waves in BoseEinstein condensates (BECs), such as solitary waves and modulate amplitude waves (MAWs). The existence and multiplicity of such waves are determined by the parameter regions selected. The results show that the characteristic of coherent matter waves can be determined by the “angular momentum” in attractive BECs while for repulsive BECs; the waves of the coherent form are all MAWs. All exact explicit parametric representations of the above waves are exhibited and numerical simulations support the result.
1. Introduction
Particles in a dilute gas reside in the same quantum (ground) state at low temperatures, which form a BoseEinstein condensate (BEC). This phenomenon was first observed experimentally in 1995 with vapors of rubidium and sodium [1]. Since then the dynamics of BEC have received intensive interest both in mathematics and physics, including superfluid and dissipative dynamics [2], Bloch waves [3], the dynamical stability [4, 5], chaos [6, 7], and dark and bright soliton [8, 9].
Considering only twobody interactions, based on the meanfield theory, the BEC is described by the GrossPitaevskii equation (GPE), which is a onedimensional nonlinear Schrödinger equation [10] where is the mass of gas particle, and are the dilute gas parameter, and is the twobody swave scattering length, which is determined by the atomic species in the condensate. Interactions between atoms are repulsive when and attractive when .
Realizable potentials of particular theoretical and experimental interest usually include harmonic traps [11], optical lattices and superlattices [12], and superpositions of lattices or superlattices with harmonic traps. However, also in many cases, weak potentials or approximately constant potentials are typically used, such as with . Therefore, one can use the Lindstedts method and the multiple scale perturbation theory [13, 14] or the method of averaging [15–17] to study dynamics of BEC with the high order approximation. For , the solutions of system (1) with coherent structures may be equilibrium points, standing waves, solitary waves, or even unbounded. While the small perturbed parameter , for example, BoseEinstein condensed atoms being perturbed by a weak optical lattice potential, the dynamical characters of these solutions may be preserved [14, 15] as periodic or quasiperiodic modulated amplitude waves. However, under some conditions of the adjustable parameters such characters may be destroyed, and the spatial chaos in the BEC occurs [6]. Consequently, it is important to investigate the dynamics of the unperturbed system drastically and entirely. Some related papers we can refer to are [18, 19].
In this paper, we will give a full topological classification of coherent structural solutions for system (1) with various adjustable parameters, while restricting the external potentials . The plane of parameters is divided into six regions. Unbounded solutions, modulated amplitude waves, and solitary waves are determined in different regions. Moreover, all exact explicit representations of solutions are exhibited in this paper. For generality, along with our analysis and derivation, we restrict the complex solutions of system (1) with nonconstant phases, which leads to nonlinear equations with singularities. Therefore, the solutions will lose smoothness in phase space. This phenomenon has been studied in other physical models by some authors (see [20, 21]).
2. Coherent Structures
We consider uniformly propagating coherent structures with the ansatz where gives the amplitude dynamics of the wave function with , gives the phase dynamics, is the particle velocity, and is the BEC’s chemical potential. When the (temporally periodic) coherent structure (2) is also spatially periodic, it is called periodic modulated amplitude waves (MAWs) [22]; while the coherent structure (2) is also spatially quasiperiodic, it is called quasiperiodic modulated amplitude waves (QMAWs) [15, 16].
Substituting (2) into the GP equation (1), and equating the real and the imaginary components of the desired equation, we obtain which yields the following nonlinear Duffing equation with regularities: The parameter is given by the relation which indicates conservation of angular momentum [23]. Null angular momentum solutions satisfying are standing waves of system (1), which have been described in detail [14]. Here, we restrict that the parameter is nonzero.To the best of our knowledge, we have not seen any overall classification of solutions in this case. Consequently, (4) becomes where and . Consequently, (6) is equivalent to the twodimensional system with the first integral
3. Bifurcations of Phase Portraits
The spatial dynamics of amplitudes for the wave functions in BECs are determined by nonlinear system (7), which depends on the parameter variables , , . The parameter is related to the chemical potential and the external potential , and has the same sign of the scatter length . The parameter comes from the constant of integration, playing the role of “angular momentum”. Without loss of generality, we assume that .
For a fixed positive constant , the characteristics of the orbits in the phase plane are affected by two parameters . Consequently, in the following we shall give an overfull classification of the orbits by the regions selecting for . The bifurcations of phase portraits will be also investigated.
Let . Then, excepting the straight line , the system (7) has the same topological phase portraits as the following system: Now, the straight line is an integral invariant straight line of (9). Denote that , then we have . To find equilibrium points of (9) which are corresponding to the constantamplitude waves, we have to solve the following algebraic equation: Note that when , . Thus, we have , which implies the relations between and in the parameter plane determining the existence and multiplicity of equilibrium points. An example of the number of roots depending on the sign of  is shown in Figure 1.
In the parameter plane, both the curve and the straight line partition the plane into six regions. Thus, we have to take six cases to be considered.(i)If , , there are two equilibrium points of (9): (ii)If , , there are two equilibrium points of (9): (iii)If , , there are two equilibrium points of (9): where (iv)If , , there are two equilibrium points of (9): where (v)If , , there are four equilibrium points of (9): where satisfies (vi)If , , there are no equilibrium points of (9).
Let be the coefficient matrix of the linearized system of (9) at an equilibrium point . Then, we have By the theory of planar dynamical systems, we know that for an equilibrium point of a planar integrable system if , then the equilibrium point is a saddle point; if and , then it is a center point; if and , then it is a node; if and the index of the equilibrium point is , then it is a cusp; otherwise, it is a high order equilibrium point.
For the function defined by (8), we denote that
We next use the above statements to consider the bifurcations of the phase portraits of (7). The possible qualitative dynamics (for ) are illustrated in Figure 2 and summarized in Table 1.

(a) (I)
(b) (II)
(c) (III)
(d) (IV)
(e) (V)
(f) (VI)
4. Exact Explicit Representations of Matter Waves
In this section, we give all exact explicit parametric representations of solitary wave solutions and periodic wave solutions (MAWs or QMAWs). Denote the Jacobian elliptic function with the modulus and the Legendre incomplete elliptic integral of the third kind (see [24]). Two types of BECs will be considered, respectively.
4.1. Attractive BECs ()
(a) Suppose that (V); that is, , , . In this case, we have the phase portrait of (7) as shown in Figure 2(e).
(1) Notice that with , corresponding to the curve defined by , and there is a homoclinic orbit connecting the equilibrium point or . The arch curve has the algebraic equation
By using the first equation of (7), for (6), we obtain the parametric representation of a smooth solitary wave solution of valley type and a smooth solitary wave solution of peak type as follows: where .
Thus, by integrating (5) on both sides we have where
Therefore, system (1) has the following solution: with
(2) Recall that with the equilibriums , , where the function is defined by (8). For each , defines two families of periodic orbits of system (7) enclosing the center or , which lie on both sides of the straight line . Those orbits determine infinitely many (positive measure set) periodic wave solutions of system (1). Using factorization of polynomial, the arch curve has the algebraic equation where . From (7), we obtain the parametric representation of the periodic wave solutions as follows: where The least positive period of is given by
Similarly, we have where .
Consequently, system (1) has the following solution of the MAW: with
(b) Suppose that (II) or (VI). In this situation, all solutions of (6) are unbounded except two equilibriums , in case of the parameters belonging to the region (II).
4.2. Repulsive BECs
For repulsive BECs, the phase portraits of amplitude equation (6) have been shown in Figures 2(a), 2(c), and 2(d), respectively. The equivalent system (7) always has two center equilibriums . Let . For each , corresponding to , system (6) have two families of periodic orbits enclosing the two centers, respectively. All these periodic orbits, filling the whole plane except for the line , determine infinitely many periodic wave solutions of (1). Similar to Section 4.1, system (1) have the parametric representation of these arches as (28) and (31) by the coherent form (2). In this case, none of solitary waves of the coherent form can exist while the external potential is constant.
5. Numerical Simulations
In order to explain the previous results better, we perform numerical simulations for solutions of solitary waves and MAWs. We emphasize the best complex case; that is, (V), and other cases can be dealt with similarly. Now let us fix the parameters , which are corresponding to the attractive interactions in BECs.
While taking the “angular momentum” , we have four equilibriums for system (7). For example, when taking the constant of integration , we have : (−6.71409, 0), : (−4.19004, 0), : (4.19004, 0), and : (6.71409, 0). As shown in Figure 2(e), there are two homoclinic orbits lying in the intersection of the stable manifold and the unstable manifold of the equilibriums and , respectively. More precisely, there exist two solutions , such that
For each fixed “angular momentum” , near the homoclinic orbit is corresponding to a dark soliton, while near the homoclinic orbit is corresponding to a bright soliton. The height of solitons, , is varying as the value of and the phase is also monotonous increasing with respect to the spatial , which can be shown in Figure 3.
(a) The portrait of the amplitudes
(b) The portrait of the phases
Also, for each fixed , near the equilibriums and , there are infinitely many periodic orbits of (7) as shown in Figure 4.
(a) The portrait of the amplitudes
(b) The portrait of the phases
While the constant exceeds or arrives at the threshold , all those periodic orbits or solitons are destroyed and become unbounded. This indicates that the characteristic of coherent matter waves can be determined by the “angular momentum”.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by the National Natural Science Foundation of China (Grants nos. 11301106, 11226130, 11361017, and 11261013).
References
 K. B. Davis, M.O. Mewes, M. R. Andrews et al., “BoseEinstein condensation in a gas of sodium atoms,” Physical Review Letters, vol. 75, no. 22, pp. 3969–3973, 1995. View at: Publisher Site  Google Scholar
 S. Burger, F. S. Cataliotti, C. Fort et al., “Superfluid and dissipative dynamics of a BoseEinstein condensate in a periodic optical potential,” Physical Review Letters, vol. 86, no. 20, pp. 4447–4450, 2001. View at: Publisher Site  Google Scholar
 B. Wu, R. B. Diener, and Q. Niu, “Bloch waves and bloch bands of BoseEinstein condensates in optical lattices,” Physical Review A, vol. 65, no. 2, Article ID 025601, 4 pages, 2002. View at: Google Scholar
 A. Chen, S. Wen, and W. Huang, “Existence and orbital stability of periodic wave solutions for the nonlinear Schrödinger equation,” The Journal of Applied Analysis and Computation, vol. 2, no. 2, pp. 137–148, 2012. View at: Google Scholar  MathSciNet
 B. Deconinck, B. A. Frigyik, and J. N. Kutz, “Dynamics and stability of BoseEinstein condensates: the nonlinear Schrödinger equation with periodic potential,” Journal of Nonlinear Science, vol. 12, no. 3, pp. 169–205, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Chong, W. Hai, and Q. Xie, “Spatial chaos of trapped BoseEinstein condensate in onedimensional weak optical lattice potential,” Chaos, vol. 14, no. 2, pp. 217–223, 2004. View at: Publisher Site  Google Scholar
 G. Chong, W. Hai, and Q. Xie, “Controlling chaos in a weakly coupled array of BoseEinstein condensates,” Physical Review E, vol. 71, no. 1, Article ID 016202, 2005. View at: Publisher Site  Google Scholar
 P. Dyke, S. Lei, and R. Hulet, “Phasedependent Interactions of Bright MatterWave Solitons,” Bulletin of the American Physical Society, vol. 57, no. 5, Article ID U6. 00001, 2012. View at: Google Scholar
 Y. Xu, Y. Zhang, and B. Wu, “Bright solitons in spinorbitcoupled BoseEinstein condensates,” Physical Review A, vol. 87, no. 1, Article ID 013614, 2013. View at: Google Scholar
 B. Schlein and H. T. Yau, “Derivation of the GrossPitaevskii equation for the dynamics of BoseEinstein condensate,” Annals of Mathematics, vol. 172, no. 1, pp. 291–370, 2010. View at: Google Scholar
 W. Krauth, “Quantum Monte Carlo calculations for a large number of bosons in a harmonic trap,” Physical Review Letters, vol. 77, no. 18, pp. 3695–3699, 1996. View at: Google Scholar
 S. Nascimbène, Y. Chen, M. Atala et al., “Experimental realization of plaquette resonating valencebond states with ultracold atoms in optical superlattices,” Physical Review Letters, vol. 108, no. 20, Article ID 205301, 2012. View at: Google Scholar
 V. P. Chua and M. A. Porter, “Spatial resonance overlap in BoseEinstein condensates in superlattice potentials,” International Journal of Bifurcation and Chaos, vol. 16, no. 4, pp. 945–959, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Porter and P. G. Kevrekidis, “BoseEinstein condensates in superlattices,” SIAM Journal on Applied Dynamical Systems, vol. 4, no. 4, pp. 783–807, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Q. Liu and D. Qian, “Modulated amplitude waves with nonzero phases in BoseEinstein condensates,” Journal of Mathematical Physics, vol. 52, no. 8, Article ID 082702, 11 pages, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Q. Liu and D. Qian, “Construction of modulated amplitude waves via averaging in collisionally inhomogeneous BoseEinstein condensates,” Journal of Nonlinear Mathematical Physics, vol. 19, no. 2, Article ID 1250017, 14 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Porter, P. G. Kevrekidis, and B. A. Malomed, “Resonant and nonresonant modulated amplitude waves for binary BoseEinstein condensates in optical lattices,” Physica D, vol. 196, no. 12, pp. 106–123, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. Li, J. Han, and F. Wang, “The extended Riccati equation method for travelling wave solutions of ZK equation,” The Journal of Applied Analysis and Computation, vol. 2, no. 4, pp. 423–430, 2012. View at: Google Scholar  MathSciNet
 S. Chow, M. Jiang, and X. Lin, “Traveling wave solutions in coupled Chua's circuits, part I: periodic solutions,” The Journal of Applied Analysis and Computation, vol. 3, no. 3, pp. 213–237, 2013. View at: Google Scholar
 J. Li and Z. Liu, “Traveling wave solutions for a class of nonlinear dispersive equations,” Chinese Annals of Mathematics. Series B, vol. 23, no. 3, pp. 397–418, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 J. Li and G. Chen, “On a class of singular nonlinear traveling wave equations,” International Journal of Bifurcation and Chaos, vol. 17, no. 11, pp. 4049–4065, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Brusch, A. Torcini, M. van Hecke, M. G. Zimmermann, and M. Bär, “Modulated amplitude waves and defect formation in the onedimensional complex GinzburgLandau equation,” Physica D, vol. 160, no. 34, pp. 127–148, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. C. Bronski, L. D. Carr, B. Deconinck, and J. N. Kutz, “BoseEinstein condensates in standing waves: the cubic nonlinear Schrödinger equation with a periodic potential,” Physical Review Letters, vol. 86, no. 8, pp. 1402–1405, 2001. View at: Publisher Site  Google Scholar
 P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientists, vol. 67 of Die Grundlehren der mathematischen Wissenschaften, Springer, New York, NY, USA, 2nd edition, 1971. View at: Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 Leilei Jia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.