^{1}

^{*}

^{1}

In this paper, we solve chiral nonlinear Schrodinger equation (CNSE) numerically. Two numerical methods are derived using the explicit Runge-Kutta method of order four and the linear multistep method (Predictor-Corrector method of fourth order). The resulting schemes of fourth order accuracy in spatial and temporal directions. The CNSE is non-integrable and has two kinds of soliton solutions: bright and dark soliton. The exact solutions and the conserved quantities of CNSE are used to display the efficiency and robustness of the numerical methods we derived. Interaction of two bright solitons for different parameters is also displayed.

The chiral nonlinear Schrödinger equation (CNSE) [

i ∂ ψ ∂ t + 1 2 ∂ 2 ψ ∂ x 2 + i λ ( ψ ∗ ∂ ψ ∂ x − ψ ∂ ψ ∂ x ∗ ) ψ = 0 , − ∞ < x < ∞ (1)

where t is time, x is spatial coordinate, ψ ( x , t ) is a complex-valued function, and λ is a nonlinear coupling constant appears through derivative coupling. This kind of nonlinearity is also known as the current density, unlike the case of cubic nonlinearity which is also known as the Kerr nonlinearity. The Chiral nonlinear Schrödinger is a non-integrable equation by the classical method of Inverse Scattering method. The single bright soliton solution of Equation (1) is given by [

ψ ( x , t ) = A sech ( β ( x − v t ) ) exp ( i ( v x + ω t ) ) , (2)

where

β = A 2 λ v and ω = v ( 2 λ A 2 − v ) 2 , (3)

A is the amplitude of the soliton, β is the inverse width of the soliton, v is the soliton velocity, and ω is the wave number. The dark soliton solution for (1) exists for λ v < 0 , and has the form

ψ ( x , t ) = A tanh ( β ( x − v t ) ) exp ( i ( v x + ω t ) ) (4)

where

β = A − 2 λ v and ω = v ( 4 λ A 2 − v ) 2 . (5)

Thus Equation (1) has bright or dark solitons that are given by (2) and (4), respectively, depending on the sign of λ v . This phenomenon makes the solitons chiral. Equation (1) has at least four conserved quantities [

I 1 ( t ) = ∫ − ∞ ∞ | ψ ( x , t ) | 2 d x , (6)

I 2 ( t ) = i ∫ − ∞ ∞ ( ψ ¯ ( x , t ) ψ x ( x , t ) − ψ ( x , t ) ψ ¯ x ( x , t ) ) d x , (7)

I 3 ( t ) = ∫ − ∞ ∞ [ λ | ψ ( x , t ) | 4 + 2 i ( ψ ¯ ( x , t ) ψ x ( x , t ) − ψ ( x , t ) ψ ¯ x ( x , t ) ) ] d x , (8)

and

I 4 = ∫ − ∞ ∞ [ λ | ψ x ( x , t ) | 2 + 2 i ( ψ ¯ ( x , t ) ψ x ( x , t ) − ψ ( x , t ) ψ ¯ x ( x , t ) ) + λ | ψ ( x , t ) | 6 ] d x (9)

Due to the exponential decay of the bright soliton solution (2) when | x | → ∞ , the conserved quantities (6)-(9) are well defined. By using (2), the exact values of the conserved quantities are

I 1 = β λ v , I 2 = 2 β λ v , I 3 = 4 A 2 3 β ( λ A 2 + 6 v ) , I 4 = 4 A 2 3 b ( 15 v 2 + 5 β 2 + 40 λ v A 2 + 8 λ 2 A 4 ) . (10)

The conserved quantities (6)-(9) using dark soliton solution (4) are not well defined due to the nonzero boundary condition as | x | → ∞ .

By assuming ψ ( x , t ) = u ( x , t ) + i v ( x , t ) , where u ( x , t ) , v ( x , t ) are real functions, the CNSE (1) can be written as the nonlinear coupled system [

∂ u ∂ t + 1 2 ∂ 2 v ∂ x 2 + 2 λ ( u ∂ v ∂ x − v ∂ u ∂ x ) v = 0, (11)

∂ v ∂ t − 1 2 ∂ 2 u ∂ x 2 − 2 λ ( u ∂ v ∂ x − v ∂ u ∂ x ) u = 0. (12)

The resulting system obtained (11)-(12) can be displayed in a matrix vector form as

∂ w ∂ t + 1 2 A ∂ 2 w ∂ x 2 + G ( w ) A w = 0 , (13)

where

w = [ u v ] , A = ( 0 1 − 1 0 ) , G ( w ) = 2 λ ( u ∂ v ∂ x − v ∂ u ∂ x ) . (14)

There are many theoretical and numerical studies in the literature about the Nonlinear Schrödinger Equations (NLS). Most of these works are motivated to single NLS and the coupled NLS (see [

Step 1: We convert the partial differential equation into a system of first order ordinary differential equation by approximating the space derivatives in the CNSE using fourth order finite difference approximation order.

Step 2: We solve the ordinary differential system obtained in Step 1 using Runge-Kutta method of fourth order and Predictor Corrector method of fourth order.

Several Numerical examples single and interaction of two solitons will be investigated to confirm the efficiency of the derived schemes.

The paper is organized as follows: In Section 2, two methods of lines are derived to solve the chiral NLS (1). In Section 3, some numerical results single soliton and interaction of two solitons to check the efficiency of the proposed schemes. Finally, the conclusion is given in Section 4.

We will consider the numerical solution of the nonlinear system (10)-(11) in a finite interval [ x L , x R ] . We assume x m = x L + m h , where m = 1,2, ⋯ , M − 1 , and h is called the space grid size, also we assume t n = n k , k is the time step size. We denote the exact and numerical solutions at the grid point ( x m , t n ) by ( w m n , v m n ) and ( U m n , V m n ) , respectively. Also we approximate the space derivatives using the fourth order approximation [

∂ u ( x m ) ∂ x = 1 12 h [ U m − 2 − 8 U m − 1 + 8 U m + 1 − U m + 2 ] (15)

∂ 2 u ( x m ) ∂ x 2 = 1 12 h 2 [ − U m − 2 + 16 U m − 1 − 30 U m + 16 U m + 1 − U m + 2 ] . (16)

By using these approximations into (11)-(12), we will get the first order ordinary differential coupled system

U ˙ m + 1 24 h 2 [ − U m − 2 + 16 U m − 1 − 30 U m + 16 U m + 1 − U m + 1 ] + λ 6 h ( U m [ V m − 2 − 8 V m − 1 + 8 V m + 1 − V m + 2 ] − V m [ U m − 2 − 8 U m − 1 + 8 U m + 1 − U m + 2 ] ) V m = 0, (17)

U − 1 = U 0 = 0 , U M = U M + 1 = 0 , for m = 1 , 2 , ⋯ , M − 1

V ˙ m − 1 24 h 2 [ − U m − 2 + 16 U m − 1 − 30 U m + 16 U m + 1 − U m + 1 ] − λ 6 h ( U m [ V m − 2 − 8 V m − 1 + 8 V m + 1 − V m + 2 ] − V m [ U m − 2 − 8 U m − 1 + 8 U m + 1 − U m + 2 ] ) U m = 0, (18)

V − 1 = V 0 = 0 , V M = V M + 1 = 0 , for m = 1 , 2 , ⋯ , M − 1 ,

The numerical solution of the coupled system (17)-(18) can be obtained by using the Runge-Kutta method of fourth order or by using a linear multistep method in the following manner. First we write the previous system (17)-(18) in the coupled form

U ˙ m = F ( U m , V m ) , m = 1 , 2 , ⋯ , M − 1 (19)

V ˙ m = F ( U m , V m ) , m = 1 , 2 , ⋯ , M − 1 (20)

The basic idea of the method of lines is to replace the spatial derivatives in the partial differential equation with algebraic approximations. Once this is carried out, the spatial derivatives are no longer stated explicitly in terms of the spatial independent variables. Thus, in effect we have a coupled system of ordinary differential equations that approximate the original PDE. Two methods for solving the ODEs coupled system (19)-(20) are presented; Runge-Kutta method of fourth order and the fourth order Predictor-Corrector methods. The details of these methods will be given next.

In this method we use one-step method, to achieve high accuracy it is typically necessary to use multistage method, where the intermediate values of the solution and its derivatives are generated and used within a single time step. To apply this approach, we use Runge-Kutta method of fourth order to solve the ordinary differential coupled systems (ODEs) [

First we define

U n = [ U 1 n , U 2 n , ⋯ , U M − 1 n ] t and V n = [ V 1 n , V 2 n , ⋯ , V M − 1 n ] t ,

and then we define the Runge-Kutta method as follows

U n + 1 = U n + 1 6 [ K 11 + 2 K 12 + 2 K 13 + K 14 ] (21)

V m n + 1 = V n + 1 6 [ K 21 + 2 K 22 + 2 K 23 + K 24 ] (22)

where

K 11 = h F ( U n , V m n ) , K 21 = h G ( U n , V n ) , (23)

K 12 = h F ( U n + 1 2 K 11 , V n + 1 2 K 21 ) , (24)

K 22 = h G ( U n + 1 2 K 11 , V n + 1 2 K 21 ) , (25)

K 13 = h F ( U n + 1 2 K 12 , V n + 1 2 K 22 ) , (26)

K 23 = h G ( U n + 1 2 K 12 , V n + 1 2 K 22 ) , (27)

K 14 = h F ( U n + K 13 , V n + K 23 ) , (28)

K 24 = h G ( U n + K 13 , V n + K 23 ) (29)

Note to apply the Runge-Kutta method, we need to evaluate eight vector functions at each time step.

One approach to use the linear multistep method is to use a predictor-corrector method, in which an explicit Adams-Basforth method of fourth order is used to predict a value for numerical solution U ( p ) , and then the Adams-Moulton of method of fourth order is used to correct this value. This is done by using U ( p ) on the right hand side of the Adams -Moulton method inside the function evaluation, so that the Adams-Moulton is no longer implicit. Now to apply this method in our problem, we execute the following steps.

Step 1 Starting Values

Use the explicit Runge-Kutta method of fourth order

U m n + 1 = U m n + 1 6 [ K 11 + 2 K 12 + 2 K 13 + K 14 ] (30)

V m n + 1 = V m n + 1 6 [ K 21 + 2 K 22 + 2 K 23 + K 24 ] (31)

where

K 11 = h F ( U n , V m n ) , K 21 = h G ( U n , V n ) , (32)

K 12 = h F ( U n + 1 2 K 11 , V n + 1 2 K 21 ) , (33)

K 22 = h G ( U n + 1 2 K 11 , V n + 1 2 K 21 ) , (34)

K 13 = h F ( U n + 1 2 K 12 , V n + 1 2 K 22 ) , (35)

K 23 = h G ( U n + 1 2 K 12 , V n + 1 2 K 22 ) , (36)

K 14 = h F ( U n + K 13 , V n + K 23 ) , (37)

K 24 = h G ( U n + K 13 , V n + K 23 ) (38)

to calculate the missing values

U n , V m , for n = 1,2,3 and m = 1,2, ⋯ , M − 1 (39)

Step 2: Predictor (Adams Bashforth of fourth order)

U ( p ) = U n + k 24 [ 55 F ( U n , V n ) − 59 F ( U n − 1 , V n − 1 ) + 37 F ( U n − 2 , V n − 2 ) − 9 F ( U n − 3 , V n − 3 ) ] (40)

V ( p ) = V n + k 24 [ 55 G ( U n , V n ) − 59 G ( U n − 1 , V n − 1 ) + 37 G ( U n − 2 , V n − 2 ) − 9 G ( U n − 3 , V n − 3 ) ] (41)

for n = 3,4, ⋯ , N .

Step 3: Corrector (Adams Moulton of fourth order)

U n + 1 = U n + k 24 [ 9 F ( U p , V p ) + 19 F ( U n , V n ) − 5 F ( U n − 1 , V n − 1 ) + F ( U n − 2 , V n − 2 ) ] (42)

V n + 1 = V n + k 24 [ 9 G ( U p , V p ) + 19 G ( U n , V n ) − 5 G ( U n − 1 , V n − 1 ) + G ( U n − 2 , V n − 2 ) ] (43)

for n = 3,4, ⋯ , N .

Note that step 1 is used to find the starting values using Runge-Kutta method of fourth order. We conclude that, at each time step, the numerical solution obtained by executing Step 2 and Step 3. Note that we need only four function evaluation at each time step, which is almost half of the Runge-Kutta method (21)-(22).

In this section, we will test the efficiency of the numerical schemes we derived, by considering different numerical tests. Trapezoidal rule is used to calculate the conserved quantities. Define the error norms

L ∞ ( u ) = max 1 ≤ m ≤ M − 1 | u m n − U m n | (44)

L 2 ( u ) = [ h ∑ m = 1 M − 1 | u m n − U m n | 2 ] . (45)

To study the behavior of the single bright soliton solution, we choose the initial condition

ψ ( x , 0 ) = A sech ( β x ) exp ( i v x ) (46)

together with the homogenous Dirichlet boundary conditions ψ ( x , 0 ) = 0 at x = x L , x R . The following set of parameters are used

x L = − 50.0 , x R = 50.0.0 , h = 0.1 , k = 0.001 , A = 0.5 , λ = 0.5 , v = 0.5 , t = 0 , 1 , 2 , ⋯ , 20. (47)

T | I_{1} | I_{2} | I_{3} | I_{4} | L_{∞} | L_{2} |
---|---|---|---|---|---|---|

0 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | - | - |

5 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.00000 | 0.000001 |

10 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000001 |

15 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000002 |

20 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000003 |

T | I_{1} | I_{2} | I_{3} | I_{4} | L_{∞} | L_{2} |
---|---|---|---|---|---|---|

0 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | - | - |

5 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.00000 | 0.000001 |

10 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000001 |

15 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000002 |

20 | 1.414214 | 1.413330 | 2.944511 | 3.044373 | 0.000001 | 0.000003 |

In non-integrable equations, stable solitons may exist, but their collisions are generally inelastic. To study the interaction behavior of the two bright solitons in the chiral nonlinear Schrodinger, we choose the initial condition [

ψ ( x , 0 ) = ψ 1 ( x , 0 ) + e δ ψ 2 ( x , 0 ) , (48)

where

ψ 1 ( x , 0 ) = A 1 sech ( β 1 ( x − x 1 ) ) exp ( i ( v 1 ( x − x 1 ) ) ) , (49)

ψ 2 ( x , 0 ) = A 2 sech ( β 1 ( x − x 2 ) ) exp ( i ( v 2 ( x − x 2 ) ) ) , (50)

and δ ∈ [ − π , π ] is the phase difference between the solitons. The initial condition represents the sum of two well separated single bright solitons. To study the interaction scenario the following parameters are selected

h = 0.1 , k = 0.001 , x l = − 50 , x 1 = − x 2 = 20 , A 1 = 1.0 , v 1 = 1.0 , A 2 = 0.5 , v 2 = 0.1 , λ = 0.5 , t = 0 , 1 , 2 , ⋯ , 80. (51)

The conserved quantities of the interaction scenario are given in _{1}. In

In Figures 4-9, we display the interaction scenario of two solitons for different parameters. Different scenario obtained, the interaction are similar to those produced by conventional nonlinear Schrodinger equation dynamics. Also we have noticed in

t | I_{1} | I_{2} | I_{3} | I_{4} |
---|---|---|---|---|

0 | 5.164771 | 4.618162 | 10.169401 | 10.881909 |

20 | 5.164771 | 4.575337 | 10.105110 | 10.811868 |

40 | 5.164771 | 4.551480 | 10.070448 | 10.725963 |

60 | 5.164771 | 4.503832 | 9.998313 | 10.703303 |

80 | 5.164769 | 4.479468 | 9.499313 | 10.661822 |

To study the behavior of the single dark soliton solution, we choose the initial condition

ψ ( x , 0 ) = A tanh ( β x ) exp ( i v x ) (52)

and use the Dirichlet boundary conditions at x = x L , x R which are extracted from the exact solution (Boundary conditions are nonzero). The following set of parameters are used

x L = − 50.0 , x R = 50 , h = 0.1 , k = 0.001 , A = 0.5 , λ = − 0.5 , v = 0.5 , t = 0 , 0.4 , 0.8 , ⋯ , 20. (53)

In _{∞} and L_{2} error norms and the conserved quantity I_{1}, we found that the proposed methods are highly accurate. The simulation of the single dark soliton displayed in

The method of lines is used to solve the Chiral nonlinear Schrodinger equation, using two methods of fourth order accuracy in both directions space and time. The numerical results using single soliton, and the interaction of two solitons of the CNSE confirmed the efficiency, reliability, and accuracy of the present methods. Concerning the cpu time required for execution, the merit goes for the Predictor-Corrector method. We can conclude that, the methods presented in this paper can be applied to similar type of equations, like nonlinear Schrodinger and coupled nonlinear Schrodinger equations.

T | I_{1} | L_{∞} | L_{2} |
---|---|---|---|

0 | 23586786 | 0.000000 | 0.000000 |

5 | 23.585682 | 0.000062 | 0.00000 |

10 | 23.585674 | 0.000065 | 0.000001 |

15 | 23.585675 | 0.000067 | 0.000001 |

20 | 23585677 | 0.000066 | 0.000001 |

The authors declare no conflicts of interest regarding the publication of this paper.

AL-Basyouni, K.S. and Ismail, M.S. (2020) Method of Lines for the Chiral Nonlinear Schrödinger Equation. Applied Mathematics, 11, 447-459. https://doi.org/10.4236/am.2020.116032