# Soliton Solutions and Collisions for the Multicomponent Gross–Pitaevskii Equation in Spinor Bose–Einstein Condensates

Authors: Ming Wang; Guo-Liang He
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

In this paper, we investigate a five-component Gross–Pitaevskii equation, which is demonstrated to describe the dynamics of anF=2 spinor Bose–Einstein condensate in one dimension. By employing the Hirota method with an auxiliary function, we obtain the explicit bright one- and two-soliton solutions for the equation via symbolic computation. With the choice of polarization parameter and spin density, the one-soliton solutions are divided into four types: one-peak solitons in the ferromagnetic and cyclic states and one- and two-peak solitons in the polar states. For the former two, solitons share the similar shape of one peak in all components. Solitons in the polar states have the one- or two-peak profiles, and the separated distance between two peaks is inversely proportional to the value of polarization parameter. Based on the asymptotic analysis, we analyze the collisions between two solitons in the same and different states.


---

## Body

## 1. Introduction

Bose–Einstein condensates (BECs) of the alkali-metal-atom gases have attracted certain attention in both experimental [1, 2] and theoretical studies [3–5]. BECs can have the internal degrees of freedom associated with the hyperfine spin and such condensates are usually called the spinor BECs [6]. The spinor BECs are classified according to the relative values of certain characteristic scattering lengths [7]. For the spinor BEC with the hyperfine spin F=1, the condensate can have either ferromagnetic state and polar state, such as 87Rb and 23Na, respectively, and the corresponding ground-state structure and dynamical properties of those two cases are distinct [8]. In the description of [7, 8], a spinor condensate with the spin F=2 bosons in an optical trap can be one of the three types: ferromagnetic, polar, and cyclic. Dynamics of the spinor BEC have been addressed in [9–14].The dynamics of nonlinear phenomena can be analyzed by means of the corresponding nonlinear evolution equations. There has been considerable work carried out on the control problem of nonlinear systems, such as those in chaotic and stochastic systems [15–19]. Solitons, resulted from the balance between the effects of nonlinearity and dispersion [20], have been studied in fields such as nonlinear optics, plasma physics, and condensed matter physics [3, 21–23]. In addition, breathers and rogue waves have been concerned in recent research [24, 25]. In the spinor BECs, matter-wave solitons are thought to be useful for their application in the atom laser, atom interferometry, and coherent atom transport, which can contribute to the realization of quantum information processors or computation [26]. The dynamics of magnetic soliton, dark soliton, and dark-bright vector soliton in spinor Bose–Einstein condensates have been investigated [27–30]. A frontier is the model extension from the single-component condensate to the multicomponent, which opens up new fields for the study of quantum matter waves and fluids [11]. The spinor BECs with spin F is described by a 2F+1-component mean-field wave function [31].In this paper, we will consider a five-component Gross–Pitaevskii (GP) equation for the dynamics of anF=2 spinor BEC in one dimension [32]:(1a)i∂tΦ±2+∂x2Φ±2+2nΦ±2−ΘΦ∓2∗=0,(1b)i∂tΦ±1+∂x2Φ±1+2nΦ±1+ΘΦ∓1∗=0,(1c)i∂tΦ0+∂x2Φ0+2nΦ0−ΘΦ0∗=0,where Φ±j’s j=0,1,2 are the wave functions of the five spin components, t and x, respectively, denote the time and spatial co-ordinate, and the asterisk represents the complex conjugation, while the number density n and singlet-pair amplitude Θ are defined by [10, 33](2)n=∑k=−22ΦkΦk∗,Θ=2Φ2Φ−2−2Φ1Φ−1+Φ02.Equation (1) is the integrable version of the multicomponent GP equation mentioned in [31] and simplified under the condition 2m=1, ℏ=1, c2=0, c4=1, and c0=−2<0, which corresponds to the attractive interaction. Due to the integrability, equation (1) has the infinite conservation laws, which restricts the dynamics of the system in an important way. The normalization is imposed as ∫dxn=NT, where NT is the total number of atoms. For equation (1), bright one- and two-soliton solutions have been given via the Zakharov–Shabat dressing [32]. However, the solitons are not classified according to the forms of the solutions and then the properties of soliton collisions are not clearly analyzed, that is the aim of this paper.In this paper, we will concentrate on the soliton types and collisions in the same and different states in anF=2 spinor BEC. For equation (1), we will present the Hirota form with an auxiliary function and explicit one- and two-soliton solutions. With symbolic computation [34, 35], this paper will be organized as follows. In Section 2, the Hirota bilinearization procedure for equation (1) will be presented. In Section 3, based on the bilinear form, we will obtain and analyze the one-soliton solutions, which reveal the one- or two-peak solitons in three different states. In Section 4, making the asymptotic analysis on the obtained two-soliton solutions, we will discuss the soliton collisions in the same and different states. Section 5 will be our conclusions.


## 2. Hirota Bilinearization Method for Equation (1)

In order to understand the dynamics of equation (1), it is essential to obtain the soliton solutions. The Hirota bilinear method is a tool to construct the soliton solutions for certain nonlinear evolution equations [36]. For equation (1), we will utilize the Hirota method with an auxiliary function to obtain the one- and two-soliton solutions. The procedure can be extended to obtain N-soliton solutions [37].To begin with, equation (1) can be expressed in the bilinear form:(3a)iDt+Dx2G±2⋅F=SG∓2∗,(3b)iDt+Dx2G±1⋅F=−SG∓1∗,(3c)iDt+Dx2G0⋅F=SG0∗,(3d)Dx2F⋅F=2∑k=−22GkGk∗,(3e)SF=2G2G−2−2G1G−1+G02,with an auxiliary function S and the following transformations:(4)Φ±j=G±jF,forj=0,1,2,where Gj’s are the complex functions of t and x and F is a real one. Dx and Dt are the bilinear differential operators [36] defined by(5)DxlDtmf⋅g=∂∂x−∂∂x′l∂∂t−∂∂t′mfx,tgx′,t′|x′=x,t′=t.where l and m are the positive integers and x′ and t′ are the formal variables. Multisoliton solutions for equation (1) can be generated by solving equation (3) with the following power series expansions as(6)G±j=εg±j1+ε3g±j3+ε5g±j5+⋯,F=1+ε2f2+ε4f4+ε6f6+⋯,S=ε2s2+ε4s4+ε6s6+⋯,where ε is the formal parameter, g±jl’s (l=1,3,5,⋯) are the complex functions of t and x, and fm’s (m=2,4,6,⋯) are the real ones which will be determined later.


## 3. Bright One-Soliton Solutions

In order to obtain one-soliton solutions for equation (1), the power series expansions for G±j, F, and S are terminated as(7)G±j=εg±j1+ε3g±j3,S=ε2S2,F=1+ε2f2+ε4f4.It is noted that equation (7) are usually taken to obtain the two-soliton solutions. Substituting equation (7) into equation (3) and collecting the terms with the same power of ε, we obtain the one-soliton solutions for equation (1):(8)Φ±j=b1±jeθ1+e2θ1+θ1∗+d1±j1+eθ1+θ1∗+l1+e2θ1+2θ1∗+c11,and the auxiliary function S given as(9)S=Λ1e2θ1,with(10)θ1=k1x+ik12t,ed1±j=δjb1∓j∗Λ12k1+k1∗2,δj=−1j=11j=0,2,el1=∑k=02b1±k2k1+k1∗2,ec11=Λ124k1+k1∗4,Λ1=2b1+1b1−1+b102−2b1+2b1−2,where b1±j’s and k1 are all complex constants. Here, we give the spin densities f=fx,fy,fz with [32](11)fx=2Φ−2Φ−1∗+62Φ0∗Φ+1+Φ−1+Φ2Φ1∗R,fy=2iΦ−2Φ−1∗+62Φ0∗Φ−1+62Φ0Φ1∗+Φ1Φ2∗I,fz=2Φ−22+Φ−12−Φ12−2Φ22,where the subscripts R and I, respectively, denote the real and imaginary parts and the total spin can be given with ∫dxf=FT. For equation (1), there exist different types of solitons according to the values of Λ1. When Λ1=0, equation (8) admit the one-peak solitons in the ferromagnetic and cyclic states corresponding to the spin densities f>0 and f=0, respectively. For the case of Λ1≠0, we can find the one- and two-peak solitons in the polar state, while only the two-peak soliton in the polar state for the F=1 spinor BECs [38].


### 3.1. One-Peak Solitons in the Ferromagnetic and Cyclic States

WhenΛ1=0, equation (9) can be expressed as(12)Φ±j=A±jeiθ1Isechθ1R+l12,with(13)A±j=12b1±je−l1/2,k1=k1R+ik1I,θ1R=k1Rx−2k1It,θ1I=k1Ix+k1R2−k1I2t.Here,A±j, 2kI, and l1/2kR, respectively, characterize the amplitude, velocity, and initial phase of the soliton. For f>0, it corresponds to the one-peak soliton in the ferromagnetic state [31]. When choosing suitable parameters to make f=0, we find that the soliton is in the cyclic state. Those two types of solitons posses the similar intensity profile of one peak in all components, as shown in Figures 1 and 2. For simplicity, we choose the parameters to make Φ+j=Φ−j in Figure 1. It is noted that the auxiliary function S vanishes in this case, as given in equation (12).Figure 1
Intensity plots of the one soliton in the ferromagnetic state for the componentsΦ±j via solution (8). Parameters are k1=1.2+i, b1±2 = 0.5, b10 = 0.5, and b1±1 = 0.78.Figure 2
Intensity plots of the one soliton in the cyclic state for the componentsΦ2 and Ψ−1 via solution (8). Parameters are k1 = 1.2+i, b12 = 0.5+0.7i, b1−2 = b10 = b11 = 0, and b1−1 = 20.5−0.7i.


### 3.2. One- and Two-Peak Solitons in the Polar State

For the caseΛ1≠0, equation (8) can be expressed as(14)Φ±j=A±jeiθ1Icoshθ1R+c11/4+P±j4cosh2θ1R+c11/4+L,where A±j=2ed1±j+lnb1±j−c11/2, L=el1−c11/2−2, and P±j=e2d1±j−2lnb1±j−c11/4. The form of equation (14) will be of use in the following asymptotic analysis of the two-soliton solutions. When f=0 under certain constraint, equation (14) depicts the one-peak soliton in the polar state, as the short dashed lines shown in Figure 3. When f≠0 but the spin F=∫fdx=0, we can observe the two-peak soliton in the polar state, as the solid lines shown in Figure 3. In addition, as the value of Λ1 decreases, the separated distance between two peaks increases, and it looks like that there exist two ferromagnetic solitons traveling parallel in the same velocity and amplitudes as Λ1 approaches to zero, as the long dashed lines shown in Figure 3.Figure 3
Intensity plots of the one soliton in the polar state for the componentsΦ±j via Solutions (8). Parameters are k1 = 1.2+i with (1) b12 = b1−2∗ = 0.5+0.8i, b10 = 1.01, and b11 = −b1−1∗ = 0.75+0.6i, short dashed lines; (2) b1±2 = 0.6, b10 = 1, and b1±1 = 0.75, solid lines; (3) b1±2 = 0.6, b10 = 1, and b1±1 = 0.85999, long dashed lines.


## 3.1. One-Peak Solitons in the Ferromagnetic and Cyclic States

WhenΛ1=0, equation (9) can be expressed as(12)Φ±j=A±jeiθ1Isechθ1R+l12,with(13)A±j=12b1±je−l1/2,k1=k1R+ik1I,θ1R=k1Rx−2k1It,θ1I=k1Ix+k1R2−k1I2t.Here,A±j, 2kI, and l1/2kR, respectively, characterize the amplitude, velocity, and initial phase of the soliton. For f>0, it corresponds to the one-peak soliton in the ferromagnetic state [31]. When choosing suitable parameters to make f=0, we find that the soliton is in the cyclic state. Those two types of solitons posses the similar intensity profile of one peak in all components, as shown in Figures 1 and 2. For simplicity, we choose the parameters to make Φ+j=Φ−j in Figure 1. It is noted that the auxiliary function S vanishes in this case, as given in equation (12).Figure 1
Intensity plots of the one soliton in the ferromagnetic state for the componentsΦ±j via solution (8). Parameters are k1=1.2+i, b1±2 = 0.5, b10 = 0.5, and b1±1 = 0.78.Figure 2
Intensity plots of the one soliton in the cyclic state for the componentsΦ2 and Ψ−1 via solution (8). Parameters are k1 = 1.2+i, b12 = 0.5+0.7i, b1−2 = b10 = b11 = 0, and b1−1 = 20.5−0.7i.


## 3.2. One- and Two-Peak Solitons in the Polar State

For the caseΛ1≠0, equation (8) can be expressed as(14)Φ±j=A±jeiθ1Icoshθ1R+c11/4+P±j4cosh2θ1R+c11/4+L,where A±j=2ed1±j+lnb1±j−c11/2, L=el1−c11/2−2, and P±j=e2d1±j−2lnb1±j−c11/4. The form of equation (14) will be of use in the following asymptotic analysis of the two-soliton solutions. When f=0 under certain constraint, equation (14) depicts the one-peak soliton in the polar state, as the short dashed lines shown in Figure 3. When f≠0 but the spin F=∫fdx=0, we can observe the two-peak soliton in the polar state, as the solid lines shown in Figure 3. In addition, as the value of Λ1 decreases, the separated distance between two peaks increases, and it looks like that there exist two ferromagnetic solitons traveling parallel in the same velocity and amplitudes as Λ1 approaches to zero, as the long dashed lines shown in Figure 3.Figure 3
Intensity plots of the one soliton in the polar state for the componentsΦ±j via Solutions (8). Parameters are k1 = 1.2+i with (1) b12 = b1−2∗ = 0.5+0.8i, b10 = 1.01, and b11 = −b1−1∗ = 0.75+0.6i, short dashed lines; (2) b1±2 = 0.6, b10 = 1, and b1±1 = 0.75, solid lines; (3) b1±2 = 0.6, b10 = 1, and b1±1 = 0.85999, long dashed lines.


## 4. Two-Soliton Solutions and Soliton Collisions

In this paper, we will concentrate on the two-soliton solutions and analyze the collisions between two solitons. Employing the following expansions,(15)G±j=εg±j1+ε3g±j3+ε5g±j5+ε7g±j7,S=ε2s2+ε4s4+ε6s6,F=1+ε2f2+ε4f4+ε6f6+ε8f8.We obtain the two-soliton solutions for equation (1) in the form of(16)Φ±j=G±jF,j=0,1,2,where the functions G±j and F are expressed as(17)G±j=b1±jeθ1+b2±jeθ2+e2θ1+θ1∗+d1±j+e2θ1+θ2∗+d12±j+e2θ2+θ1∗+d21±j+e2θ2+θ2∗+d2±j+eθ1+θ1∗+θ2+D12±j+eθ2+θ2∗+θ1+D21±j+e2θ1+2θ1∗+θ2+h12±j+e2θ2+2θ2∗+θ1+h21±j+e2θ2+2θ1∗+θ1+H12±j+e2θ1+2θ2∗+θ2+H21±j+e2θ1+θ1∗+θ2+θ2∗+w12±j+e2θ2+θ2∗+θ1+θ1∗+w21±j+e2θ1+2θ1∗+2θ2+θ2∗+W12±j+e2θ1+2θ2+2θ2∗+θ1∗+W21±j,(18)F=1+eθ1+θ1∗+l1+eθ1+θ2∗+l0+eθ2+θ1∗+l0∗+eθ2+θ2∗+l2+e2θ1+2θ1∗+c11+e2θ1+2θ2∗+c12+e2θ2+2θ1∗+c12∗+e2θ2+2θ2∗+c22+e2θ1+θ1∗+θ2∗+l12+e2θ1∗+θ1+θ2+l12∗+e2θ2+θ1∗+θ2∗+l21+e2θ2∗+θ1+θ2+l21∗+eθ1+θ1∗+θ2+θ2∗+l3+e2θ1+2θ1∗+θ2+θ2∗+L12+e2θ2+2θ2∗+θ1+θ1∗+L21+e2θ1+2θ2∗+θ2+θ1∗+l4+e2θ2+2θ1∗+θ1+θ2∗+l4∗+e2θ1+2θ1∗+2θ2+2θ2∗+l5,and the auxiliary function S takes the form(19)S=Λ1e2θ1+Λ2e2θ2+eθ1+θ2+r1+eθ1+θ1∗+2θ2+r21+eθ2+θ2∗+2θ1+r12+eθ1+θ2∗+2θ2+R21+eθ2+θ1∗+2θ1+R12+e2θ1+2θ1∗+2θ2+r3+e2θ1+2θ2+2θ2∗+r4+e2θ1+θ1∗+2θ2+θ2∗+r5,with(20)θm=kmx+ikm2t,Λm=2bm+1bm−1+bm02−2bm+2bm−2,m=1,2,edm±j=δjbm∓j∗Λm2km+km∗2,edlm±j=δjbl∓j∗Λmkm+kl∗2,l=1,2,m≠l,δj=−1,j=1,1,j=0,2,eDlm±j=1km+km∗2kl+km∗22km−klbm±jkm+km∗∑k=02bm±k∗bl±k−bl±jkl+km∗∑k=02δkbm±k2+δjbm∓j∗km+km∗kl+km∗∑k=02δkbm±kbl±k,ehlm±j=bl±jΛm2km∗−kl∗24km+km∗4kl+km∗2,eHlm±j=bm∓j∗ΛmΛl∗km∗−kl∗24km+km∗2km+kl∗4,ewlm±j=Λmkm−kl22km+km∗2km+kl∗4kl+kl∗2bl±jkl+kl∗kl+km∗∑k=02δkbm∓k∗bl±k∗+δjkm∗−kl∗bm∓j∗kl+km∗∑k=02bl±k2−bl∓j∗kl+kl∗∑k=02bm±k∗bl±k,eWlm±j=δjbl∓j∗Λm2Λlkm−kl4km∗−kl∗28km+km∗4kl+km∗4km+kl∗2kl+kl∗2,el0=∑k=02b1±kb2±k∗k1+k2∗2,elm=∑k=02bm±k2km+km∗2,eclm=ΛmΛl∗4km+kl∗4,el3=2k1+k1∗+k2+k2∗2k12+k22−2k1∗k2ec02−ec1ec2−k12+k22−2k1k2ec02+ec1ec2−2∑m,l=1,m≠l2∑k=02bm±k∗eDlm±kR,ellm=Λm∑k=02δkb1∓k∗b2±k∗2km+km∗2km+kl∗2,eLlm=Λm2∑k=02bl±k2k1−k244km+km∗4km+kl4kl+kl∗2,el4=Λ1Λ2∗∑k=02b2±kb1±k∗k1−k244k1+k1∗2k1+k24k1+k2∗2k2+k2∗2,em5=Λ12Λ22k1−k2816k1+k1∗4k1+k2∗8k2+k2∗4,er1=2∑k=02δkb1∓kb2±k,erlm=Λm∑k=02bl±k2km−kl2km+kl∗2kl+kl∗2,eRlm=Λm∑k=02bl±kbm±k∗k1−k22kl+km∗2km+km∗2,er3=Λ2Λ12k1−k244k1+k1∗4k2+k1∗4,er4=Λ1Λ22k1−k244k2+k2∗4k1+k2∗4,er5=Λ1Λ2∑k=02δkb1∓k∗b2±k∗k1−k242k1+k1∗2k1+k2∗4k2+k2∗2,where k1, k2, b1±j’s, and b2±j’s are all complex constants. To understand the collision property of solitons in the F=2 spinor BEC, we will make the asymptotic analysis on two-soliton solution (16) as t⟶∓∞, which are associated with the initial and final state, respectively. Without loss of generality, we set(21)k1R>0,k1I>0,k2R>0,k2I<0.Following the classification on the one-soliton solutions, we naturally consider the collisions between the same and different types of solitons, which are determined byΛm and f. Since the one-peak soliton in the cyclic state is obtained under certain strong restriction, its collisions with other types of soliton are omitted here.


### 4.1. Collisions between Two Solitons in the Ferromagnetic State

At first, we will consider the collision between two solitons in the ferromagnetic state. Under the conditionΛm = 0, the asymptotic forms for equation (16) are given as follows:(A)
Before the collisiont⟶−∞.(1)
SolitonS1−θ1R∼0,θ2R⟶−∞:(22)Φ±j⟶A±j1−eiθ1Isechθ1R+l12.(2)
SolitonS2−θ2R∼0,θ1R⟶∞:(23)Φ±j⟶A±j2−eiθ2Isechθ2R+l3−l12.(B)
After the collisiont⟶+∞.(1)
SolitonS1+θ1R∼0,θ2R⟶∞:(24)Φ±j⟶A±j1+eiθ1Isechθ1R+l3−l22.(2)
SolitonS2+θ2R∼0,θ1R⟶−∞:(25)Φ±j⟶A±j2+eiθ2Isechθ2R+l22,with(26)A±j1−=b1±j2e−l1/2,A±j2−=12eD12±j−l2−l3/2,A±j1+=12eD21±j−l2−l3/2,A±j2+=b2±j2e−l2/2.The superscriptsm±m=1,2 represent the solitons indicated as S1 and S2 at t⟶±∞. Setting b1±j/b2±j = constant, we find that the amplitudes A±jm± before and after the collision can be related through A±j1+ = A±j1−k1−k2k1∗+k2/k1∗−k2∗k1+k2∗ and A±j2+ = A±j2−k1∗−k2∗k1∗+k2/k1−k2k1+k2∗, namely, A±jm+ = A±jm−, which means that the collision is elastic in this case. Besides, the solitons S1 and S2, respectively, experience a phase shift of η1=l3−l1−l2/4k1R=lnk1−k2k1∗−k2∗/k1∗+k2k1+k2∗/k1R and η2=−η1k1R/k2R during the collision. For other choices of bm±j satisfying Λm = 0, we will take an example and give Figure 4. The intensity of soliton S2 in component Φ0 decreases while those of S2 in components Φ±1 increase after the collision, as shown in Figures 4(b) and 4(c). This process might be of value for the realization of soliton switching devices [39]. The similar changes of soliton intensities can be observed for soliton S1. For this case, the collision is similar to that with an intensity redistribution among all components of each soliton in the multicomponent NLS equation [40].Figure 4
Collision between two solitons in the ferromagnetic state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.6, b10 = 0.9, b1±1 = 0.9, b2±2 = 1.1, b20 = 0.08, and b2±1 = 1.25.


### 4.2. Collisions between Soliton in the Ferromagnetic State and Solitons in the Polar State

We will analyze the collisions between soliton in the ferromagnetic state and one- and two-peak solitons in the polar state. Here, we chooseS1 to be the soliton in the ferromagnetic state (Λ1= 0), and S2 to be the soliton in the polar state corresponding to the condition Λ2≠0. The asymptotic forms of solitons S1 and S2 are given as follows:(A)
Before the collision.(1)
SolitonS1−:(27)Φ±j⟶A±j1−eiθ1Isechθ1R+l12,whereA±j1−=1/2b1±je−l1/2.(2)
SolitonS2−:(28)Φ±j⟶A±j2−eiθ2Icoshθ2R−+P±j2−4cosh2θ2R−+L2−,whereA±j2−=2eW21±j+D12±j−L21−l1/2, L2−=el3−L21+l1/2−2, P±j2−=e2W21±j−2D12±j−L21+l1/4, and θ2R−=θ2R+L21−l1/4, with(29)k2=k2R+ik2I,θ2I=k2Ix+k2R2−k2I2t,θ2R=k2Rx−2k2Rk2It.(B)
After the collision.(1)
SolitonS1+:(30)Φ±j⟶A±j1+eiθ1Isechθ1R+L12−c222,whereA±j1+=1/2e2h12±j−L12−c22/2.(2)
SolitonS2+:(31)Φ±j⟶A±j2+eiθ2Icoshθ2R++P±j2+4cosh2θ2R++L2+,whereA±j2−=2ed2±j+lnb2±j−c22/2, L2−=el2−c22/2−2, P±j2−=e2d2±j−2lnb2±j−c22/4, and θ2R−=θ2R+c22/4. It is noted that equations (32) and (34) have the similar forms as equation (14). Choosing suitable parameters, we give Figure 5 illustrating the collision between one-peak solitons in the ferromagnetic and polar states. The solitons S1 in the ferromagnetic state retain their shapes invariant after the collision, and only has a phase shift at the moment of collision in all components, due to the total spin conservation. However, the amplitudes of one-peak soliton S2 in the polar state increase in components Φ±2 after the collision, accompanied by the decrease of amplitude for S2 in component Φ0. For the case of collision between one-peak soliton in the ferromagnetic state and two-peak soliton in the polar state, the collision occurs without change of intensities for the solitons in the ferromagnetic state, while with an intensity redistribution among the two-peak solitons in the polar state in all components, as shown in Figure 6. Therefore, we can conclude that there exists no mixing collision in the internal states for the solitons in the ferromagnetic state outside the collision region.Figure 5
Collision between one-peak solitons in the ferromagnetic and polar states for the componentsΦ±j via solution (16). Parameters are the same as those in Figure 4 except b22 = b2−2∗ = 0.8+0.5i, b20 = 0.3, and b21 = −b2−1∗ = 0.4−0.9i.Figure 6
Collision between one-peak soliton in the ferromagnetic state and two-peak soliton in the polar state for the componentsΦ±j via solutions (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.5, b10 = 1, b1±1 = 0.75, b2±2 = 0.8, b20 = 0.4, and b2±1 = 0.65.Figure 7
Collision between one-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b12 = b1−2∗ = 0.5+0.8i, b10 = 1.01, b11 = −b1−1∗ = 0.75−0.6i, b22 = b2−2∗ = 0.8+0.5i, b20 = 0.3, and b21 = −b2−1∗ = 0.4−0.9i.Figure 8
Collision between two-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.5, b10 = 1.1, b1±1 = 0.75, b2±2 = 0.8, b20 = 0.3, and b2±1 = 0.66.Figure 9
Collision between one- and two-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are the same as those in Figure 7 except b1±2 = 0.5, b10 = 1.01, and b1±1 = 0.75.


### 4.3. Collisions between Two Solitons in the Polar State

In this section, we will consider the collisions between the solitons in the polar state (Λm≠0). In this case, the asymptotic forms are given as(A)
Before the collision.(1)
SolitonS1−:(32)Φ±j⟶A±j1−eiθ1Icoshθ1R−+P±j1−4cosh2θ1R−+L1−,whereA±j1−=2ed1±j+lnb1±j−c11/2, L1−=el1−c11/2−2, P±j1−=e2d1±j−2lnb1±j−c11/4, and θ1R−=θ1R+c11/4.(2)
SolitonS2−:(33)Φ±j⟶2A±j2−eiθ2Icoshθ2R−+P±j2−4cosh2θ2R−+L2−,whereA±j2−=2eW12±j+h12±j−l5−c11/2, L2−=eL12−l5+c11/2−2, P±j2−=e2W12±j−2h12±j−l5+c11/4, and θ2R−=θ2R+l5−c11/4.(B)
After the collision.(1)
SolitonS1+:(34)Φ±j⟶2A±j1+eiθ1Icoshθ1R++P±j1+4cosh2θ1R++L1+,whereA±j1+=2eW21±j+h21±j−l5−c22/2=k1−k2k1∗+k2/k1∗−k2∗k1+k2∗A±j1−, L1+=eL21−l5+c22/2−2, θ1R+=θ1R+l5−c22/4, and P±j1+=e2W21±j−2h21±j−l5+c22/4.(2)
SolitonS2+:(35)Φ±j⟶2A±j2+eiθ2Icoshθ2R++P±j2+4cosh2θ2R++L2+,whereA±j2+=2ed2±j+lnb2±j−c22/2=k1∗−k2∗k1∗+k2/k1−k2k1+k2∗A±j2−,L2+=el2−c22/2−2, P±j2+=e2d2±j−2lnb2±j−c22/4, and θ2R+=θ2R+c22/4. Figures 7–9, respectively, describe the collisions between two one-peak solitons, two two-peak solitons, and one-peak and two-peak solitons in the polar state. From equations (32)–(35), one can find that the collisions are elastic in this case. As shown in Figures 7–9, solitons with one or two peaks collide with each other and the intensity and shape of each soliton keeps invariant after the collision.


## 4.1. Collisions between Two Solitons in the Ferromagnetic State

At first, we will consider the collision between two solitons in the ferromagnetic state. Under the conditionΛm = 0, the asymptotic forms for equation (16) are given as follows:(A)
Before the collisiont⟶−∞.(1)
SolitonS1−θ1R∼0,θ2R⟶−∞:(22)Φ±j⟶A±j1−eiθ1Isechθ1R+l12.(2)
SolitonS2−θ2R∼0,θ1R⟶∞:(23)Φ±j⟶A±j2−eiθ2Isechθ2R+l3−l12.(B)
After the collisiont⟶+∞.(1)
SolitonS1+θ1R∼0,θ2R⟶∞:(24)Φ±j⟶A±j1+eiθ1Isechθ1R+l3−l22.(2)
SolitonS2+θ2R∼0,θ1R⟶−∞:(25)Φ±j⟶A±j2+eiθ2Isechθ2R+l22,with(26)A±j1−=b1±j2e−l1/2,A±j2−=12eD12±j−l2−l3/2,A±j1+=12eD21±j−l2−l3/2,A±j2+=b2±j2e−l2/2.The superscriptsm±m=1,2 represent the solitons indicated as S1 and S2 at t⟶±∞. Setting b1±j/b2±j = constant, we find that the amplitudes A±jm± before and after the collision can be related through A±j1+ = A±j1−k1−k2k1∗+k2/k1∗−k2∗k1+k2∗ and A±j2+ = A±j2−k1∗−k2∗k1∗+k2/k1−k2k1+k2∗, namely, A±jm+ = A±jm−, which means that the collision is elastic in this case. Besides, the solitons S1 and S2, respectively, experience a phase shift of η1=l3−l1−l2/4k1R=lnk1−k2k1∗−k2∗/k1∗+k2k1+k2∗/k1R and η2=−η1k1R/k2R during the collision. For other choices of bm±j satisfying Λm = 0, we will take an example and give Figure 4. The intensity of soliton S2 in component Φ0 decreases while those of S2 in components Φ±1 increase after the collision, as shown in Figures 4(b) and 4(c). This process might be of value for the realization of soliton switching devices [39]. The similar changes of soliton intensities can be observed for soliton S1. For this case, the collision is similar to that with an intensity redistribution among all components of each soliton in the multicomponent NLS equation [40].Figure 4
Collision between two solitons in the ferromagnetic state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.6, b10 = 0.9, b1±1 = 0.9, b2±2 = 1.1, b20 = 0.08, and b2±1 = 1.25.


## 4.2. Collisions between Soliton in the Ferromagnetic State and Solitons in the Polar State

We will analyze the collisions between soliton in the ferromagnetic state and one- and two-peak solitons in the polar state. Here, we chooseS1 to be the soliton in the ferromagnetic state (Λ1= 0), and S2 to be the soliton in the polar state corresponding to the condition Λ2≠0. The asymptotic forms of solitons S1 and S2 are given as follows:(A)
Before the collision.(1)
SolitonS1−:(27)Φ±j⟶A±j1−eiθ1Isechθ1R+l12,whereA±j1−=1/2b1±je−l1/2.(2)
SolitonS2−:(28)Φ±j⟶A±j2−eiθ2Icoshθ2R−+P±j2−4cosh2θ2R−+L2−,whereA±j2−=2eW21±j+D12±j−L21−l1/2, L2−=el3−L21+l1/2−2, P±j2−=e2W21±j−2D12±j−L21+l1/4, and θ2R−=θ2R+L21−l1/4, with(29)k2=k2R+ik2I,θ2I=k2Ix+k2R2−k2I2t,θ2R=k2Rx−2k2Rk2It.(B)
After the collision.(1)
SolitonS1+:(30)Φ±j⟶A±j1+eiθ1Isechθ1R+L12−c222,whereA±j1+=1/2e2h12±j−L12−c22/2.(2)
SolitonS2+:(31)Φ±j⟶A±j2+eiθ2Icoshθ2R++P±j2+4cosh2θ2R++L2+,whereA±j2−=2ed2±j+lnb2±j−c22/2, L2−=el2−c22/2−2, P±j2−=e2d2±j−2lnb2±j−c22/4, and θ2R−=θ2R+c22/4. It is noted that equations (32) and (34) have the similar forms as equation (14). Choosing suitable parameters, we give Figure 5 illustrating the collision between one-peak solitons in the ferromagnetic and polar states. The solitons S1 in the ferromagnetic state retain their shapes invariant after the collision, and only has a phase shift at the moment of collision in all components, due to the total spin conservation. However, the amplitudes of one-peak soliton S2 in the polar state increase in components Φ±2 after the collision, accompanied by the decrease of amplitude for S2 in component Φ0. For the case of collision between one-peak soliton in the ferromagnetic state and two-peak soliton in the polar state, the collision occurs without change of intensities for the solitons in the ferromagnetic state, while with an intensity redistribution among the two-peak solitons in the polar state in all components, as shown in Figure 6. Therefore, we can conclude that there exists no mixing collision in the internal states for the solitons in the ferromagnetic state outside the collision region.Figure 5
Collision between one-peak solitons in the ferromagnetic and polar states for the componentsΦ±j via solution (16). Parameters are the same as those in Figure 4 except b22 = b2−2∗ = 0.8+0.5i, b20 = 0.3, and b21 = −b2−1∗ = 0.4−0.9i.Figure 6
Collision between one-peak soliton in the ferromagnetic state and two-peak soliton in the polar state for the componentsΦ±j via solutions (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.5, b10 = 1, b1±1 = 0.75, b2±2 = 0.8, b20 = 0.4, and b2±1 = 0.65.Figure 7
Collision between one-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b12 = b1−2∗ = 0.5+0.8i, b10 = 1.01, b11 = −b1−1∗ = 0.75−0.6i, b22 = b2−2∗ = 0.8+0.5i, b20 = 0.3, and b21 = −b2−1∗ = 0.4−0.9i.Figure 8
Collision between two-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are k1 = 1.2+i, k2 = 1.2−0.6i, b1±2 = 0.5, b10 = 1.1, b1±1 = 0.75, b2±2 = 0.8, b20 = 0.3, and b2±1 = 0.66.Figure 9
Collision between one- and two-peak solitons in the polar state for the componentsΦ±j via solution (16). Parameters are the same as those in Figure 7 except b1±2 = 0.5, b10 = 1.01, and b1±1 = 0.75.


## 4.3. Collisions between Two Solitons in the Polar State

In this section, we will consider the collisions between the solitons in the polar state (Λm≠0). In this case, the asymptotic forms are given as(A)
Before the collision.(1)
SolitonS1−:(32)Φ±j⟶A±j1−eiθ1Icoshθ1R−+P±j1−4cosh2θ1R−+L1−,whereA±j1−=2ed1±j+lnb1±j−c11/2, L1−=el1−c11/2−2, P±j1−=e2d1±j−2lnb1±j−c11/4, and θ1R−=θ1R+c11/4.(2)
SolitonS2−:(33)Φ±j⟶2A±j2−eiθ2Icoshθ2R−+P±j2−4cosh2θ2R−+L2−,whereA±j2−=2eW12±j+h12±j−l5−c11/2, L2−=eL12−l5+c11/2−2, P±j2−=e2W12±j−2h12±j−l5+c11/4, and θ2R−=θ2R+l5−c11/4.(B)
After the collision.(1)
SolitonS1+:(34)Φ±j⟶2A±j1+eiθ1Icoshθ1R++P±j1+4cosh2θ1R++L1+,whereA±j1+=2eW21±j+h21±j−l5−c22/2=k1−k2k1∗+k2/k1∗−k2∗k1+k2∗A±j1−, L1+=eL21−l5+c22/2−2, θ1R+=θ1R+l5−c22/4, and P±j1+=e2W21±j−2h21±j−l5+c22/4.(2)
SolitonS2+:(35)Φ±j⟶2A±j2+eiθ2Icoshθ2R++P±j2+4cosh2θ2R++L2+,whereA±j2+=2ed2±j+lnb2±j−c22/2=k1∗−k2∗k1∗+k2/k1−k2k1+k2∗A±j2−,L2+=el2−c22/2−2, P±j2+=e2d2±j−2lnb2±j−c22/4, and θ2R+=θ2R+c22/4. Figures 7–9, respectively, describe the collisions between two one-peak solitons, two two-peak solitons, and one-peak and two-peak solitons in the polar state. From equations (32)–(35), one can find that the collisions are elastic in this case. As shown in Figures 7–9, solitons with one or two peaks collide with each other and the intensity and shape of each soliton keeps invariant after the collision.


## 5. Conclusions

In this paper, we have investigated the solitons and their collisions for the five-component GP equations, i.e., equation (1), which can describe the dynamics of an F=2 spinor BECs in one dimension. Through the Hirota method and symbolic computation, we have derived one-soliton solutions (8) and two-soliton solutions (16) for equation (1). With different polarization parameters and spin densities, one-soliton solutions (8) have been classified as (a) the one-peak soliton in the ferromagnetic and cyclic states (as seen in Figures 1 and 2); (b) one- and two-peak solitons in the polar state (as seen in Figure 3). For the two-peak soliton in the polar state, the distance between the two peaks increases inversely with the polarization parameter Λ1 (as seen in Figure 3). Performing the asymptotic analysis on two-soliton solutions (16), we have analyzed the collisions between two solitons in the same and different states. For the collisions between two solitons in the ferromagnetic state, we have shown the constrained condition, under which the collision is elastic. In addition, we have given the example of the inelastic collision with intensity redistribution (as seen in Figure 4), and it is feasible to exploit the inelastic collision in the design of soliton switching architectures. For the collisions between soliton in the ferromagnetic state and solitons in the polar state, the solitons in the ferromagnetic state keep their characteristics invariant while the shapes and intensities change for those in the polar state after the collision (as seen in Figures 5 and 6). For the collisions between two solitons in the polar state, the procedure occurs only with a phase shift for each soliton at the moment of collision, and the collisions are inelastic (as seen in Figures 7–9). We hope that the results given in this paper will be of value for the applications in the atom laser and coherent atom transport, which can contribute to the realization of quantum information processors or computation.



---

*Source: 4632434-2020-07-31.xml*