# Rogue Waves of Nonlinear Schrödinger Equation with Time-Dependent Linear Potential Function

Authors: Ni Song; Yakui Xue
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

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.


---

## Body

## 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
=
-
ħ
2
2
m
∂
2
Ψ
∂
X
2
+
X
m
g
+
B
cos
⁡
ω
1
T
Ψ
+
4
π
ħ
2
a
m
Ψ
2
Ψ
,where Ψ
(
X
,
T
) is the wave function of Bose-Einstein condensates, m is atomic mass, a is the scattering length between atoms, m
g is the gravity field, B
cos
⁡
(
ω
1
T
) 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
=
ħ
m
l
0
2
T
,
x
=
1
L
0
X
,
φ
=
1
n
0
Ψ
.Then, (1) can be rewritten as(3)
i
∂
φ
∂
t
+
1
2
∂
2
φ
∂
x
2
+
x
b
+
l
cos
⁡
ω
t
φ
+
μ
φ
2
φ
=
0
,where (4)
b
=
-
m
2
l
0
3
g
ħ
2
,
l
=
-
m
l
0
3
B
ħ
2
,
ω
=
m
l
0
2
ω
1
ħ
,and x is measured in units of l
0
=
1
μ
m, a characteristic length unit in this type of experiment, t in units of m
l
0
2
/
ħ, and φ in units of the square root of n
0
a
3, the maximum density in the initial distribution of the condensate and the interaction constant is defined as μ
=
-
4
π
n
0
l
0
2
a, 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
+
1
2
∂
2
φ
∂
x
2
+
x
f
t
φ
+
μ
φ
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
=
φ
R
x
,
t
+
i
φ
I
x
,
t
e
i
ϕ
x
,
t
,whose intensity can be written as(7)
φ
x
,
t
2
=
φ
R
x
,
t
2
+
φ
I
x
,
t
2
,where φ
R
x
,
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
+
1
2
φ
R
,
x
x
-
φ
I
,
x
ϕ
x
-
1
2
φ
R
ϕ
x
2
-
1
2
φ
I
ϕ
x
x
+
μ
φ
R
2
+
φ
I
2
φ
R
+
x
f
φ
R
=
0
,
(8b)
φ
R
,
t
-
φ
I
ϕ
t
+
1
2
φ
I
,
x
x
+
φ
R
,
x
ϕ
x
-
1
2
φ
I
ϕ
x
2
+
1
2
φ
R
ϕ
x
x
+
μ
φ
R
2
+
φ
I
2
φ
I
+
x
f
φ
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
=
A
t
+
B
t
P
η
x
,
t
,
τ
t
,
(9b)
φ
I
=
C
t
Q
η
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)
-
C
t
Q
-
C
Q
η
η
t
+
Q
τ
τ
t
-
A
+
B
P
χ
t
+
λ
τ
t
+
B
2
P
η
η
η
x
2
+
P
η
η
x
x
-
C
Q
η
η
x
χ
x
-
1
2
A
+
B
P
χ
x
2
-
1
2
C
Q
χ
x
x
+
μ
A
+
B
P
2
+
C
2
Q
2
A
+
B
P
+
x
f
A
+
B
P
=
0
,
A
t
+
B
t
P
-
B
P
η
η
t
+
P
τ
τ
t
-
C
Q
χ
t
+
λ
τ
t
+
C
2
P
η
η
η
x
2
+
P
η
η
x
x
-
B
P
η
η
x
χ
x
-
1
2
C
Q
χ
x
2
-
1
2
A
+
B
P
χ
x
x
+
μ
A
+
B
P
2
+
C
2
Q
2
C
Q
+
x
f
C
Q
=
0
.Simplifying the above equations, we obtain the similarity reduction(11a)
η
x
x
=
0
,
(11b)
η
t
+
η
x
χ
x
=
0
,
(11c)
σ
t
+
σ
2
χ
x
x
=
0
,
σ
=
A
,
B
,
C
,
(11d)
χ
t
+
1
2
χ
x
2
-
x
f
=
0
,
(11e)
-
C
Q
τ
τ
t
-
A
+
B
P
λ
τ
t
+
B
2
P
η
η
η
x
2
+
μ
A
+
B
P
A
+
B
P
2
+
C
2
Q
2
=
0
,
(11f)
-
B
P
τ
τ
t
-
C
Q
λ
τ
t
+
C
2
Q
η
η
η
x
2
+
μ
C
Q
A
+
B
P
2
+
C
2
Q
2
=
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
=
α
t
x
+
β
t
,
(12b)
χ
x
,
t
=
-
α
t
2
α
x
2
-
β
t
α
x
+
χ
0
t
,
(12c)
A
t
=
a
0
α
,
B
t
=
b
A
,
C
t
=
d
A
,where a
0, 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
=
1
2
∫
0
t
α
2
s
d
s
,
α
2
=
2
A
2
.Therefore, (11a), (11b), (11c), (11d), (11e), and (11f) are reduced to (14a)
-
d
Q
τ
-
1
+
b
P
λ
+
b
P
η
η
+
μ
1
+
b
P
1
+
b
P
2
+
d
2
Q
2
=
0
,
(14b)
-
b
P
τ
-
d
Q
λ
+
d
Q
η
η
+
μ
d
Q
1
+
b
P
2
+
d
2
Q
2
=
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
η
,
τ
=
-
4
R
1
η
,
τ
b
,
Q
η
,
τ
=
-
8
τ
R
1
η
,
τ
d
,where R
1
(
η
,
τ
)
=
1
+
2
η
2
+
4
τ
2 and the second-order rational solution of (14a) and (14b) is given(16)
P
η
,
τ
=
P
1
η
,
τ
R
2
η
,
τ
b
,
Q
η
,
τ
=
Q
1
η
,
τ
τ
R
2
η
,
τ
d
,where (17a)
P
1
η
,
τ
=
3
8
-
9
τ
2
-
3
2
η
2
-
6
η
2
τ
2
-
10
τ
4
-
1
2
η
4
,
(17b)
Q
1
η
,
τ
=
-
15
4
+
2
τ
2
-
3
η
2
+
4
η
2
τ
2
+
4
τ
4
+
η
4
,
(17c)
R
2
η
,
τ
=
3
32
+
33
8
τ
2
+
9
16
η
2
-
3
2
η
2
τ
2
+
9
2
τ
4
+
1
8
η
4
+
2
3
τ
6
+
η
2
τ
6
+
1
2
η
4
τ
2
+
1
12
η
6
.Thus, a direct reduction solution of (5) can be derived as (18)
φ
=
A
1
+
b
P
+
i
d
Q
e
i
χ
+
τ
,where A
(
t
), P
(
η
,
τ
), Q
(
η
,
τ
), χ
(
x
,
t
), and τ
(
t
) are expressed by (12a), (12b), (12c), (13), and (16), respectively.


## 3. Rogue Waves Solutions

### 3.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)
φ
=
a
0
α
1
-
4
+
8
i
τ
1
+
2
η
2
+
4
τ
2
e
i
χ
+
τ
,whose amplitude can be written as (20)
φ
2
=
a
0
2
α
t
2
α
t
x
+
β
t
2
+
4
τ
2
-
3
2
+
64
τ
2
t
1
+
2
α
t
x
+
β
t
2
+
4
τ
2
t
2
.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
=
f
t
=
b
+
l
cos
⁡
(
ω
t
) to exhibit the nonlinear dynamic behaviors of the rogue waves which change with the gravity field b
=
δ
m
g and time-dependent external magnetic field l
cos
⁡
(
ω
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
=
a
0
2
b
2
2
b
x
+
β
t
2
+
4
τ
2
-
3
2
+
64
τ
2
t
1
+
2
b
x
+
β
t
2
+
4
τ
2
t
2
,where τ
(
t
)
=
1
/
2
b
t.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.Figure 1
The first-order rogue wave propagations for the intensityφ
2 with the gravity field b
=
0.5.Figure 2
Contour plots of the first-order rogue wave propagation with the gravity fields (a)b
=
0.5; (b) b
=
1.2.
(a)
(b)Figure 3
(a) Coordinate variation of the amplitude of the first-order rogue waves with the gravity fieldsb
=
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.
(a)
(b)Case 2.
Suppose that the free functionβ
(
t
) is the polynomial with respect to time t, namely, β
(
t
)
=
0.2
t
2, and the compound function β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
). 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.Figure 4
The first-order rogue waves propagation for the intensityφ
2 with the gravity field b
=
0.5 and (a)-(b) β
t
=
0.2
t
2, (c)-(d) β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
).
(a)
(b)
(c)
(d)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.2
cos
⁡
(
0.005
t
) and β
=
0.2
t
2. 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.Figure 5
(a) The first-order rogue wave propagations for the intensityφ
2 with the gravity field and external magnetic field b
=
0.8
+
1.2
cos
⁡
(
0.005
t
) and β
=
0.2
t
2. (b) Rogue wave propagations for the intensity φ
x
,
0
2 with t
=
0. (c) The comparison of the amplitude between b
=
0.8 and b
=
0.8
+
1.2
cos
⁡
(
0.005
t
).
(a)
(b)
(c)


### 3.2. Second-Order Solution

Substituting (16) into (18), we obtain the second-order rational-like solution of (5)(22)
φ
=
a
0
α
1
+
P
1
R
2
+
i
Q
1
τ
R
2
e
i
χ
+
τ
,whose intensity is written as (23)
φ
2
=
a
0
2
α
t
1
+
P
1
R
2
2
+
Q
1
τ
R
2
2
,where P
1
(
η
,
τ
), Q
1
(
η
,
τ
), and R
2
(
η
,
τ
) 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 supposeb
=
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.Figure 6
(a) The second-order rogue wave propagations for the intensityφ
2. (b) Contour plot with the gravity field b
=
0.5, β
(
t
)
=
0.2.
(a)
(b)Supposeb
=
0.5 and β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
); 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.Figure 7
(a) The second-order rogue wave propagations for the intensityφ
2. (b) Contour plot with the gravity field b
=
0.5, β
t
=
0.5
exp
⁡
sech
⁡
0.2
t
3. (c) Rogue wave propagations for the intensity φ
(
x
,
0
)
2 with t
=
0.
(a)
(b)
(c)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.2
cos
⁡
(
8
t
) and β
(
t
)
=
0.2
t; 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.Figure 8
The second-order rogue wave propagations for the intensityφ
2 with the gravity field and external magnetic field α
t
=
0.5
+
1.2
cos
⁡
(
8
t
).


## 3.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)
φ
=
a
0
α
1
-
4
+
8
i
τ
1
+
2
η
2
+
4
τ
2
e
i
χ
+
τ
,whose amplitude can be written as (20)
φ
2
=
a
0
2
α
t
2
α
t
x
+
β
t
2
+
4
τ
2
-
3
2
+
64
τ
2
t
1
+
2
α
t
x
+
β
t
2
+
4
τ
2
t
2
.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
=
f
t
=
b
+
l
cos
⁡
(
ω
t
) to exhibit the nonlinear dynamic behaviors of the rogue waves which change with the gravity field b
=
δ
m
g and time-dependent external magnetic field l
cos
⁡
(
ω
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
=
a
0
2
b
2
2
b
x
+
β
t
2
+
4
τ
2
-
3
2
+
64
τ
2
t
1
+
2
b
x
+
β
t
2
+
4
τ
2
t
2
,where τ
(
t
)
=
1
/
2
b
t.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.Figure 1
The first-order rogue wave propagations for the intensityφ
2 with the gravity field b
=
0.5.Figure 2
Contour plots of the first-order rogue wave propagation with the gravity fields (a)b
=
0.5; (b) b
=
1.2.
(a)
(b)Figure 3
(a) Coordinate variation of the amplitude of the first-order rogue waves with the gravity fieldsb
=
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.
(a)
(b)Case 2.
Suppose that the free functionβ
(
t
) is the polynomial with respect to time t, namely, β
(
t
)
=
0.2
t
2, and the compound function β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
). 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.Figure 4
The first-order rogue waves propagation for the intensityφ
2 with the gravity field b
=
0.5 and (a)-(b) β
t
=
0.2
t
2, (c)-(d) β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
).
(a)
(b)
(c)
(d)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.2
cos
⁡
(
0.005
t
) and β
=
0.2
t
2. 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.Figure 5
(a) The first-order rogue wave propagations for the intensityφ
2 with the gravity field and external magnetic field b
=
0.8
+
1.2
cos
⁡
(
0.005
t
) and β
=
0.2
t
2. (b) Rogue wave propagations for the intensity φ
x
,
0
2 with t
=
0. (c) The comparison of the amplitude between b
=
0.8 and b
=
0.8
+
1.2
cos
⁡
(
0.005
t
).
(a)
(b)
(c)


## 3.2. Second-Order Solution

Substituting (16) into (18), we obtain the second-order rational-like solution of (5)(22)
φ
=
a
0
α
1
+
P
1
R
2
+
i
Q
1
τ
R
2
e
i
χ
+
τ
,whose intensity is written as (23)
φ
2
=
a
0
2
α
t
1
+
P
1
R
2
2
+
Q
1
τ
R
2
2
,where P
1
(
η
,
τ
), Q
1
(
η
,
τ
), and R
2
(
η
,
τ
) 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 supposeb
=
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.Figure 6
(a) The second-order rogue wave propagations for the intensityφ
2. (b) Contour plot with the gravity field b
=
0.5, β
(
t
)
=
0.2.
(a)
(b)Supposeb
=
0.5 and β
t
=
0.5
exp
⁡
(
sech
⁡
(
0.2
t
3
)
); 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.Figure 7
(a) The second-order rogue wave propagations for the intensityφ
2. (b) Contour plot with the gravity field b
=
0.5, β
t
=
0.5
exp
⁡
sech
⁡
0.2
t
3. (c) Rogue wave propagations for the intensity φ
(
x
,
0
)
2 with t
=
0.
(a)
(b)
(c)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.2
cos
⁡
(
8
t
) and β
(
t
)
=
0.2
t; 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.Figure 8
The second-order rogue wave propagations for the intensityφ
2 with the gravity field and external magnetic field α
t
=
0.5
+
1.2
cos
⁡
(
8
t
).


## 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.



---

*Source: 7879517-2016-02-17.xml*