The rogue waves of the nonlinear Schrödinger equation with time-dependent linear potential function are investigated by using the similarity transformation in this paper. The first-order and second-order rogue waves solutions are obtained and the nonlinear dynamic behaviors of these solutions are discussed in detail. In addition, the amplitudes of the rogue waves under the effect of the gravity field and external magnetic field changing with the time are analyzed by using numerical simulation. The results can be used to study the matter rogue waves in the Bose-Einstein condensates and other fields of nonlinear science.
1. Introduction
The rogue wave is giant single wave which was firstly found in the ocean [1]. The amplitude of this wave is two to three times higher than those of its surrounding waves. The key feature of the rogue wave is that it “appears from nowhere” and “disappears without a trace” [2]. The most terrible thing is that it is very dangerous for sailors because it can appear unexpectedly and form larger amplitude in one minute to shred a boat. Beyond oceanic expanses, the rogue wave has been also found in optical fibers [3], Bose-Einstein condensates (BECs) [4], superfluids [5], and so on. However, it is very difficult to explain the rogue waves by using the linear theories based on the superposition principles. The nonlinear theories of ocean waves [6–8] can be used to explain why the rogue waves can appear from nowhere. In recent years, it becomes an important issue for ones to study the rogue waves theoretically in the fields of the nonlinear science [9–13]. The Darboux transformation (DT) [14, 15], the similarity transformation, and the numerical simulation [13, 16–18] were used to analyze the occurrence of the rogue waves and the larger amplitudes. One of the important known models for the rogue waves is the nonlinear Schrödinger (NLS) equation which is a foundational model in describing numerous nonlinear physical phenomena and the first-order rational solution was derived by Peregrine [19] and the second-order one was obtained by using the modified Darboux transformation. Cheng et al. [20] show the controllable rogue waves in coupled NLS with varying potential functions and nonlinearities. Wu et al. [21] derive the evolution of optical solitary waves in a generalized NLS equation with variable coefficients.
Since three affirmative observations on Bose-Einstein condensates in ultracold rubidium [22], lithium [23], and sodium [24] atomic gases were reported, respectively, in 1995, researches on the behaviors of the macroscopical quantum and dynamic evolution of BECs are more active and competitive than before in both experiment and theory aspects. Ivancevic and Reid [25] show that the NLS equation is fundamental for quantum turbulence, while its closed form solutions include shock-waves, solitons, and rogue waves. Zhu [26] gives nonlinear tunneling for controllable rogue waves in two-dimensional graded-index waveguides.
A new area of physics, namely, nonlinear matter waves and nonlinear atomic optics, was originated. Generation and dynamics of solitary wave pulses in BECs are one of the most important related problems. In the approximation of the mean field, the evolution characteristics of the nonlinear excitations for BECs can be described by the Gross-Pitaevskii equation [27, 28] which is also called the nonlinear Schrödinger equation. Particularly, in the gravity field, the evolution characteristics of the nonlinear excitations for BECs, which change with the time under the external magnetic field, can be described by the continuous Gross-Pitaevskii equation [29]:(1)iħ∂Ψ∂T=-ħ22m∂2Ψ∂X2+Xmg+Bcosω1TΨ+4πħ2amΨ2Ψ,where Ψ(X,T) is the wave function of Bose-Einstein condensates, m is atomic mass, a is the scattering length between atoms, mg is the gravity field, Bcos(ω1T) is the time-dependent external magnetic field, B is the time-dependent field parameter, and ω1 is the frequency of external magnetic field.
To simplify (1), we introduce the following dimensionless transformation: (2)t=ħml02T,x=1L0X,φ=1n0Ψ.Then, (1) can be rewritten as(3)i∂φ∂t+12∂2φ∂x2+xb+lcosωtφ+μφ2φ=0,where (4)b=-m2l03għ2,l=-ml03Bħ2,ω=ml02ω1ħ,and x is measured in units of l0=1μm, a characteristic length unit in this type of experiment, t in units of ml02/ħ, and φ in units of the square root of n0a3, the maximum density in the initial distribution of the condensate and the interaction constant is defined as μ=-4πn0l02a, which is taken positive to the attractive potential function (a<0).
Equation (3) is NLS equation with time-dependent linear potential function. The single dark soliton and single bright soliton [30] are obtained by applying the classical inverse scattering method [31]. The soliton and periodic solutions are also obtained by applying the transformation from the corresponding solution of the NLS equation [32].
In this paper, we apply the similarity transformation [33] and direct hypothesis [16, 34] to obtain the analytical rational-like solution of (3) which can describe the possible formation mechanisms of optical rogue waves. In order to make the problem more general, we consider the general form of (3):(5)i∂φ∂t+12∂2φ∂x2+xftφ+μφ2φ=0.Without loss of generality, we set μ=1 for a>0 and μ=-1 for a<0.
2. Explicit Solutions through Similarity Transformation
Firstly, we construct the following transformation for the envelope field φ(x,t) in the gauge form [33]:(6)φx,t=φRx,t+iφIx,teiϕx,t,whose intensity can be written as(7)φx,t2=φRx,t2+φIx,t2,where φRx,t, φI(x,t) and phase ϕ(x,t) are real functions with respect to space x at time t.
Substituting (6) into (5) yields the following coupled partial differential equations with variable coefficients:(8a)-φI,t-φRϕt+12φR,xx-φI,xϕx-12φRϕx2-12φIϕxx+μφR2+φI2φR+xfφR=0,(8b)φR,t-φIϕt+12φI,xx+φR,xϕx-12φIϕx2+12φRϕxx+μφR2+φI2φI+xfφI=0.
For the real functions φR(x,t), φI(x,t), and ϕ(x,t), introducing the new variables η(x,t) and τ(t) and further utilizing the similarity transformations, we have(9a)φR=At+BtPηx,t,τt,(9b)φI=CtQηx,t,τt,(9c)ϕ=χx,t+λτt,where λ is a constant.
Substituting (9a), (9b), and (9c) into (8a) and (8b), we obtain the following equations:(10)-CtQ-CQηηt+Qττt-A+BPχt+λτt+B2Pηηηx2+Pηηxx-CQηηxχx-12A+BPχx2-12CQχxx+μA+BP2+C2Q2A+BP+xfA+BP=0,At+BtP-BPηηt+Pττt-CQχt+λτt+C2Pηηηx2+Pηηxx-BPηηxχx-12CQχx2-12A+BPχxx+μA+BP2+C2Q2CQ+xfCQ=0.
Simplifying the above equations, we obtain the similarity reduction(11a)ηxx=0,(11b)ηt+ηxχx=0,(11c)σt+σ2χxx=0,σ=A,B,C,(11d)χt+12χx2-xf=0,(11e)-CQττt-A+BPλτt+B2Pηηηx2+μA+BPA+BP2+C2Q2=0,(11f)-BPττt-CQλτt+C2Qηηηx2+μCQA+BP2+C2Q2=0,where η(x,t), χ(x,t), A(t), B(t), C(t), P(η,τ), and Q(η,τ) are the functions to be determined.
After some algebra computation, it follows from (11a)–(11d) that we have(12a)ηx,t=αtx+βt,(12b)χx,t=-αt2αx2-βtαx+χ0t,(12c)At=a0α,Bt=bA,Ct=dA,where a0, b, and d are constants, α(t) is the inverse of the wave width, and -β(t)/α(t) is the position of its center of mass. α(t), β(t), and χ0(t) are all free functions with respect to time t.
To further reduce (11e)-(11f) to the coupled partial differential equations, we require (13)τt=12∫0tα2sds,α2=2A2.
Therefore, (11a), (11b), (11c), (11d), (11e), and (11f) are reduced to (14a)-dQτ-1+bPλ+bPηη+μ1+bP1+bP2+d2Q2=0,(14b)-bPτ-dQλ+dQηη+μdQ1+bP2+d2Q2=0.
For simplicity, we set μ=1. According to the direct method developed in [14, 15], we obtain the first-order rational solution of (14a) and (14b) (15)Pη,τ=-4R1η,τb,Qη,τ=-8τR1η,τd,where R1(η,τ)=1+2η2+4τ2 and the second-order rational solution of (14a) and (14b) is given(16)Pη,τ=P1η,τR2η,τb,Qη,τ=Q1η,ττR2η,τd,where (17a)P1η,τ=38-9τ2-32η2-6η2τ2-10τ4-12η4,(17b)Q1η,τ=-154+2τ2-3η2+4η2τ2+4τ4+η4,(17c)R2η,τ=332+338τ2+916η2-32η2τ2+92τ4+18η4+23τ6+η2τ6+12η4τ2+112η6.
Thus, a direct reduction solution of (5) can be derived as (18)φ=A1+bP+idQeiχ+τ,where A(t), P(η,τ), Q(η,τ), χ(x,t), and τ(t) are expressed by (12a), (12b), (12c), (13), and (16), respectively.
3. Rogue Waves Solutions3.1. First-Order Solution
Now we focus on the rogue wave structures in the solutions of the nonlinear Schrödinger equation with time-dependent linear potential function. Substituting (15) into (18), we have the first-order rational-like solution of (5):(19)φ=a0α1-4+8iτ1+2η2+4τ2eiχ+τ,whose amplitude can be written as (20)φ2=a02αt2αtx+βt2+4τ2-32+64τ2t1+2αtx+βt2+4τ2t2.
It can be seen that the structure of solution (19) is different from one of the solutions obtained in [29, 31]. In what follows, we will choose the function αt=ft=b+lcos(ωt) to exhibit the nonlinear dynamic behaviors of the rogue waves which change with the gravity field b=δmg and time-dependent external magnetic field lcos(ωt), where δ is a constant.
We study the nonlinear dynamic behaviors of the rogue waves when there is only the gravity field; namely, l=0. Thus, the amplitude given by (20) is changed to(21)φ2=a02b22bx+βt2+4τ2-32+64τ2t1+2bx+βt2+4τ2t2,where τ(t)=1/2bt.
We plot figures of the amplitude using (21) with different value b and function β(t).
Case 1.
Supposing that free function β(t) is a constant and setting the gravity field b=0.5, the structure of the rogue wave described by (21) is plotted, as shown in Figure 1. For better understanding the structure of the rogue wave, b=0.5 and 1.2. The contour plots of the first-order rogue wave propagation are depicted, as shown in Figures 2(a)-2(b), respectively. It can be observed from Figures 2(a)-2(b) that the pattern of the rogue waves is the same when the gravity field b takes any value. However, it is found from Figures 2(a)-2(b) that the position x, at which the highest amplitudes of the rogue waves occur, is different. Suppose b=0.8 and b=1.5; the relations on the rogue waves amplitude via coordinate and time are, respectively, obtained, as shown in Figure 3. It is clearly seen from Figure 3(a) that, with increasing of the gravity field, the amplitude of the rogue waves has an increase in the peak density. The maximum amplitude of the rogue waves occurs at the points x=-2.4, -1.5, -1.2, and -1. It is also observed from Figure 3(b) that the time, on which the maximum amplitude of the rogue waves occurs, is the same, where t=0.
The first-order rogue wave propagations for the intensity φ2 with the gravity field b=0.5.
Contour plots of the first-order rogue wave propagation with the gravity fields (a) b=0.5; (b) b=1.2.
(a) Coordinate variation of the amplitude of the first-order rogue waves with the gravity fields b=0.5, b=0.8, b=1.2, and b=1.5. (b) Time variation of the amplitude of the first-order rogue waves with the gravity fields b=0.5, b=0.8, b=1.2, and b=1.5.
Case 2.
Suppose that the free function β(t) is the polynomial with respect to time t, namely, β(t)=0.2t2, and the compound function βt=0.5exp(sech(0.2t3)). Then, Figures 4(a) and 4(c) give the nonlinear dynamical behaviors of the rational-like solution (21). When t=0, the changing cases of φ(x,0)2 via x are depicted, as shown in Figures 4(b) and 4(d). In Figure 4(b), the position of maximum amplitude φ2 appears at x=0. However, in Figure 4(d), the position of maximum amplitude φ2 appears at x=-2.7. In two cases, it is found from Figures 4(b) and 4(d) that the maximum amplitudes equal 4.5. If different free functions β(t) are chosen, the solution determined by (21) is of different dynamic properties.
The first-order rogue waves propagation for the intensity φ2 with the gravity field b=0.5 and (a)-(b) βt=0.2t2, (c)-(d) βt=0.5exp(sech(0.2t3)).
In the following analysis, we study the nonlinear dynamic behaviors of the rogue waves when there exist the gravity field and external magnetic field; namely, l≠0. For simplicity, we suppose b=0.8+1.2cos(0.005t) and β=0.2t2. The nonlinear dynamical behavior of the rational-like solution (21) is plotted, as shown in Figure 5(a). When t=0, the changing case of φ(x,0)2 via x is depicted, as shown in Figure 5(b). In Figure 5(b), the position of maximum amplitude φ2 appears at x=0 and the maximum amplitude is equal to 18. Comparing the amplitudes of Figure 5(b) with Figure 4(b), we observe that the amplitude of the rogue wave has a large and quick increase when the external magnetic field exists, as shown in Figure 5(c). The amplitudes φ(x,0)2 reach the largest at the point x=0. Furthermore, we adjust the intensity l of the external magnetic field to analyze the changing cases of the amplitude for the rogue waves.
(a) The first-order rogue wave propagations for the intensity φ2 with the gravity field and external magnetic field b=0.8+1.2cos(0.005t) and β=0.2t2. (b) Rogue wave propagations for the intensity φx,02 with t=0. (c) The comparison of the amplitude between b=0.8 and b=0.8+1.2cos(0.005t).
3.2. Second-Order Solution
Substituting (16) into (18), we obtain the second-order rational-like solution of (5)(22)φ=a0α1+P1R2+iQ1τR2eiχ+τ,whose intensity is written as (23)φ2=a02αt1+P1R22+Q1τR22,where P1(η,τ), Q1(η,τ), and R2(η,τ) are expressed by (17a)–(17c), respectively.
The effect of the gravity field on the rogue wave is studied. Similar to the first-order solution, we suppose b=0.5 and β(t)=0.2. The nonlinear dynamic behaviors of the second-order rogue wave are depicted, as shown in Figures 6(a)-6(b). Compared with the first-order solution, it is found that there are four small peaks around one high peak in the second-order rogue waves and all energy of the rogue wave is focused on the high peak. The amplitude of the second-order rational-like solution is even larger than the one of the first-order solution. It is observed from Figure 6(a) that the maximum amplitude equals 12.5.
(a) The second-order rogue wave propagations for the intensity φ2. (b) Contour plot with the gravity field b=0.5, β(t)=0.2.
Suppose b=0.5 and βt=0.5exp(sech(0.2t3)); the nonlinear dynamic behaviors of the second-order rogue waves are also depicted, as shown in Figures 7(a)-7(b). When t=0, the changing procedure of φ(x,0)2 via x is depicted, as shown in Figure 7(c). In Figure 7(c), the position of maximum amplitude φ2 appears at x=-2.7. It is observed that the maximum amplitude equals 12.48.
(a) The second-order rogue wave propagations for the intensity φ2. (b) Contour plot with the gravity field b=0.5, βt=0.5expsech0.2t3. (c) Rogue wave propagations for the intensity φ(x,0)2 with t=0.
Then, we also study the nonlinear dynamic behaviors of the second-order rogue waves when there exist the gravity field and external magnetic field. Suppose αt=0.5+1.2cos(8t) and β(t)=0.2t; the nonlinear dynamic behaviors of the second-order rogue waves are depicted, as shown in Figure 8. It is observed from Figure 8 that the second-order rogue wave is different from the ones displayed in Figures 6(a) and 7(a). In this case, the wave width is thin and the rogue wave changes more quickly than the first-order one displayed in Figure 5(a). The maximum amplitude equals 45.
The second-order rogue wave propagations for the intensity φ2 with the gravity field and external magnetic field αt=0.5+1.2cos(8t).
4. Conclusion
In this paper, we present the rogue waves of the nonlinear Schrödinger equation with time-dependent linear potential function. Using the similarity transformation and direct hypothesis, we obtain the first-order and second-order solutions of the rogue waves and discuss the nonlinear dynamic behaviors of these solutions in detail. In addition, we also analyze effect of the gravity field and external magnetic field changing with the time on the amplitudes of the rogue waves. It is found that the gravity field and external magnetic field have significant effect on the amplitude of the rogue waves. The numerical results obtained in this paper are useful to study the matter rogue waves in the Bose-Einstein condensates and other fields of nonlinear science.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors thank sincerely Professor Wei Zhang for many useful discussions. This work is supported by the National Sciences Foundation of Shanxi Province (2015011009).
MüllerP.GarrettC.OsborneA.Rogue waves-the fourteenth 'Aha Huliko'a Hawaiian Winter Workshop2005183667510.5670/oceanog.2005.302-s2.0-84855264670AkhmedievN.AnkiewiczA.TakiM.Waves that appear from nowhere and disappear without a trace2009373667567810.1016/j.physleta.2008.12.0362-s2.0-58149490756SolliD. R.RopersC.KoonathP.JalaliB.Optical rogue waves200745071721054105710.1038/nature064022-s2.0-37249058537BludovY. V.KonotopV. V.AkhmedievN.Matter rogue waves200980303361010.1103/physreva.80.0336102-s2.0-70349332917GanshinA. N.EfimovV. B.KolmakovG. V.Mezhov-DeglinL. P.McClintockP. V. E.Observation of an inverse energy cascade in developed acoustic turbulence in superfluid helium2008101606530310.1103/physrevlett.101.0653032-s2.0-49749126376KharifC.PelinovskyE.Physical mechanisms of the rogue wave phenomenon200322660363410.1016/j.euromechflu.2003.09.002MR20238632-s2.0-3042781172JanssenP. A. E. M.Nonlinear four-wave interactions and freak waves200333486388410.1175/1520-0485(2003)33<863:nfiafw>2.0.co;2OnoratoM.OsborneA. R.SerioM.BertoneS.Freak waves in random oceanic sea states200186255831583410.1103/physrevlett.86.58312-s2.0-0035908014DaiC.-Q.WangY.-Y.TianQ.ZhangJ.-F.The management and containment of self-similar rogue waves in the inhomogeneous nonlinear Schrödinger equation2012327251252110.1016/j.aop.2011.11.0162-s2.0-84855558235WangX.-L.ZhangW.-G.ZhaiB.-G.ZhangH.-Q.Rogue waves of the higher-order dispersive nonlinear Schrödinger equation201258453153810.1088/0253-6102/58/4/15MR30977662-s2.0-84867696828WangY. Y.DaiC. Q.Spatiotemporal Rogue waves for the Variable-coefficient (3 + 1)-dimensional nonlinear Schrödinger equation201258225526010.1088/0253-6102/58/2/15WangX.-C.HeJ.-S.LiY.-S.Rogue wave with a controllable center of nonlinear Schrödinger equation201156463163710.1088/0253-6102/56/4/06ZBL1247.351552-s2.0-80054809764YanZ.-Y.Financial rogue waves201054594794910.1088/0253-6102/54/5/312-s2.0-78649761042PeregrineD. H.Water waves, nonlinear Schrödinger equations and their solutions1983251164310.1017/s0334270000003891MR702914AkhmedievN.AnkiewiczA.Soto-CrespoJ. M.Rogue waves and rational solutions of the nonlinear Schrödinger equation200980202660110.1103/physreve.80.0266012-s2.0-68949175669AkhmedievN.Soto-CrespoJ. M.AnkiewiczA.Extreme waves that appear from nowhere: on the nature of rogue waves2009373252137214510.1016/j.physleta.2009.04.023MR25271412-s2.0-67349165528YanZ. Y.Nonautonomous ‘rogons’ in the inhomogeneous nonlinear Schrödinger equation with variable coefficients2010374467267910.1016/j.physleta.2009.11.0302-s2.0-73349088704SongS. Y.WangJ.MengJ. M.WangJ. B.HuP. X.Nonlinear Schrödinger equation for internal waves in deep sea201059211231129MR2664400PeregrineD. H.Water waves, nonlinear Schrödinger equations and their solutions1983251164310.1017/s0334270000003891MR702914ChengX. P.WangJ. Y.LiJ. Y.Controllable rogue waves in coupled nonlinear Schrödinger equations with varying potentials and nonlinearities201477354555210.1007/s11071-014-1316-22-s2.0-84905081235WuX.-F.HuaG.-S.MaZ.-Y.Evolution of optical solitary waves in a generalized nonlinear Schrödinger equation with variable coefficients20127032259226710.1007/s11071-012-0616-7MR29922132-s2.0-84870923441AndersonM. H.EnsherJ. R.MatthewsM. R.WiemanC. E.CornellE. A.Observation of Bose-Einstein condensation in a dilute atomic vapor1995269522119820110.1126/science.269.5221.1982-s2.0-11944274056BradleyC. C.SackettC. A.TollettJ. J.HuletR. G.Evidence of Bose-Einstein condensation in an atomic gas with attractive interactions19957591687169010.1103/PhysRevLett.75.16872-s2.0-4243132347DavisK. B.MewesM.-O.AndrewsM. R.Van DrutenN. J.DurfeeD. S.KurnD. M.KetterleW.Bose-Einstein condensation in a gas of sodium atoms199575223969397310.1103/physrevlett.75.39692-s2.0-4244115335IvancevicV. G.ReidD. J.Turbulence and shock-waves in crowd dynamics2012681-228530410.1007/s11071-011-0227-8MR29041752-s2.0-84858705725ZhuH.-P.Nonlinear tunneling for controllable rogue waves in two dimensional graded-index waveguides201372487388210.1007/s11071-013-0759-1MR30645062-s2.0-84878588921DalfovoF.GiorginiS.PitaevskiiL. P.StringariS.Theory of Bose-Einstein condensation in trapped gases199971346351210.1103/revmodphys.71.4632-s2.0-0033246313LeggettA. J.Bose-Einstein condensation in the alkali gases: some fundamental concepts200173230735610.1103/revmodphys.73.3072-s2.0-0035535327FanW. B.2001Beijing, ChinaInstitute of Applied Physics and Computational MathematicsLiuX.-Q.JiangS.FanW.-B.LiuW.-M.Soliton solutions in linear magnetic field and time-dependent laser field20049336136510.1016/s1007-5704(02)00109-02-s2.0-0842327666AblowitzM. J.ClarksonP. A.1991Cambridge, UKCambridge University Press10.1017/cbo9780511623998MR1149378HeJ.-S.JiM.LiY.-S.Solutions of two kinds of non-isospectral generalized nonlinear schrödinger equation related to Bose-Einstein condensates2007248215710.1088/0256-307x/24/8/0012-s2.0-34547593669YanZ. Y.2007Beijing, ChinaScience PressAkhmedievN.Soto-CrespoJ. M.AnkiewiczA.How to excite a rogue wave200980404381810.1103/physreva.80.0438182-s2.0-70350736734