By using the bifurcation theory of dynamical systems, we present the exact representation and topological classification of coherent matter waves in Bose-Einstein 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 Bose-Einstein 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 two-body interactions, based on the mean-field theory, the BEC is described by the Gross-Pitaevskii equation (GPE), which is a one-dimensional nonlinear Schrödinger equation [10] where is the mass of gas particle, and are the dilute gas parameter, and is the two-body s-wave 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 [1517] 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, Bose-Einstein 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 two-dimensional 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 constant-amplitude 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.

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.

Also, for each fixed , near the equilibriums and , there are infinitely many periodic orbits of (7) as shown in Figure 4.

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.


This work is supported by the National Natural Science Foundation of China (Grants nos. 11301106, 11226130, 11361017, and 11261013).