The chiral nonlinear Schrödinger equation (CNSE)   we are going to study is given by
where t is time, x is spatial coordinate, 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  
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 , and has the form
Thus Equation (1) has bright or dark solitons that are given by (2) and (4), respectively, depending on the sign of . This phenomenon makes the solitons chiral. Equation (1) has at least four conserved quantities  , namely
Due to the exponential decay of the bright soliton solution (2) when , the conserved quantities (6)-(9) are well defined. By using (2), the exact values of the conserved quantities are
The conserved quantities (6)-(9) using dark soliton solution (4) are not well defined due to the nonzero boundary condition as .
By assuming , where are real functions, the CNSE (1) can be written as the nonlinear coupled system    as
The resulting system obtained (11)-(12) can be displayed in a matrix vector form as
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    and reference therein). However, up to the authors knowledge, only a few numerical studies for the chiral nonlinear Schrodinger NLS. In , a conservative finite difference schemes for the chiral nonlinear schrodinger equation using finite difference method, three different schemes of second order accuracy are derived, a nonlinear implicit finite difference where we have to solve a nonlinear block tridigonal systems at each time step which is expensive to some extent, a linearization approach is used to overcome this difficulty and to obtain, two linearly implicit finite difference schemes. In this work, we are going to present numerical schemes of fourth order accuracy in both direction using method of lines. The method of lines has been used by several authors. A Method of lines is used to find the numerical solution of the Kortweg-de Vries Equation . The method of lines solution of the regularized long-wave equation using Runge-Kutta time discretization method presented in . The method of lines for solution of the one-dimensional wave equation subject to an integral conservation condition presented in . Application of the method of lines for solving the KdV-Burger equation, two methods are used to solve this equation MOL and the Adomian decomposition method the results reveal that, both methods are comparable . A generalization of the method of lines for the numerical solution of the coupled forced vibration of beams presented in . In this work, we are going to solve the chiral nonlinear Schrodinger equation using method of lines, which can be described in two major steps.
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.
2. Numerical Methods
We will consider the numerical solution of the nonlinear system (10)-(11) in a finite interval . We assume , where , and h is called the space grid size, also we assume , k is the time step size. We denote the exact and numerical solutions at the grid point by and , respectively. Also we approximate the space derivatives using the fourth order approximation .
By using these approximations into (11)-(12), we will get the first order ordinary differential coupled system
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
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.
2.1. Method of Lines Using Fourth Order Runge-Kutta Method
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
and then we define the Runge-Kutta method as follows
Note to apply the Runge-Kutta method, we need to evaluate eight vector functions at each time step.
2.2. Method of Lines Using Linear Multistep Method of Fourth Order
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 , and then the Adams-Moulton of method of fourth order is used to correct this value. This is done by using 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
to calculate the missing values
Step 2: Predictor (Adams Bashforth of fourth order)
Step 3: Corrector (Adams Moulton of fourth order)
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).
3. Numerical Results
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
3.1. Bright Soliton Solution
To study the behavior of the single bright soliton solution, we choose the initial condition
together with the homogenous Dirichlet boundary conditions at . The following set of parameters are used
Table 1 and Table 2 display the errors and the conserved quantities for our proposed schemes. The numerical results are highly accurate and conserved the conserved quantities exactly. The two methods produced almost the same results (both fourth order methods in space and time directions). In Figure 1, we display the motion of the single bright solitons for . We have noticed that the cpu execution time required for producing the results in Table 1, Table 2 are 4.14 and 2.65 seconds respectively…
Table 1. Single bright soliton using Runge-Kutta method (cpu = 4.14).
Table 2. Single bright soliton using predictor-corrector method (cpu = 2.656).
Figure 1. Single Bright soliton: .
3.2. Interaction of Two Bright Solitons
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 
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
The conserved quantities of the interaction scenario are given in Table 3. We have noticed that the numerical results are highly accurate and the exact conservation of I1. In Figure 2 and Figure 3, we display the interaction scenario, and we have noticed how the two solitons approach each other and the faster soliton emerges to the slower soliton. Eventually, the two original solitons leave the interaction region with their original shape.
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 Figure 6 and Figure 7, fussion of solitons is observed by choosing .
Figure 2. Interaction of two solitons: .
Figure 3. Interaction of two solitons: .
Figure 4. Interaction of two solitons: .
Figure 5. Interaction of two solitons: .
Figure 6. Interaction of two solitons: .
Figure 7. Interaction of two solitons: .
Figure 8. Interaction of two solitons: .
Figure 9. Interaction of two solitons: .
Table 3. Interaction of two bright solitons using Runge-Kutta method.
3.3. Single Dark Soliton
To study the behavior of the single dark soliton solution, we choose the initial condition
and use the Dirichlet boundary conditions at which are extracted from the exact solution (Boundary conditions are nonzero). The following set of parameters are used
In Table 4, we calculate the L∞ and L2 error norms and the conserved quantity I1, we found that the proposed methods are highly accurate. The simulation of the single dark soliton displayed in Figure 10 for . The cpu execution time is 4.14 section, which is the same as the bright soliton.
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.
Figure 10. Dark soliton: with .
Table 4. Single Dark soliton using Runge-Kutta method (cpu = 4.14).
 Ismail, M.S. and Taha, T.R. (2001) Numerical Simulation of Coupled Nonlinear SchrÖdinger Equation. Mathematics and Computers in Simulation, 56, 547-562.
 Ismail, M.S. and Alamri, S.Z. (2004) Highly Accurate Finite Difference Method for Coupled Nonlinear SchrÖdinger Equation. International Journal of Computer Mathematics, 81, 333-351.
 Ismail, M.S. and Taha, T.R. (2007) A Linearly Implicit Conservative Scheme for the Coupled Nonlinear SchrÖdinger Equation. Mathematics and Computers in Simulation, 74, 302-311.
 Ismail, M.S., Al-Basyouni, K.S. and Aydin, A. (2015) Conservative Finite Difference Schemes for the Chiral Nonlinear Schrodinger Equation. Boundary Value Problems, 2015, Article No. 89.
 Bakodah, H.O. and Banaja, M.A. (2013) The Method of Lines Solution of the Regularized Long-Wave Equation Using Runge-Kutta Time Discretization Method. Mathematical Problems in Engineering, 2013, Article ID: 804317
 Shakeri, F. and Dehghan, M. (2008) The Method of Lines for Solution of the One-Dimensional Wave Equation Subject to an Integral Conservation Condition. Computers and Mathematics with Applications, 56, 2175-2188.
 Sarker, P. and Chakravarty, U.K. (2020) A Generalization of the Method of Lines for the Numerical Solution of Coupled, Forced Vibration of Beams. Mathematics and Computers in Simulation, 170, 115-142.
 LeVeque, R.J. (2007) Finite Difference Methods for Ordinary and Partial Differential Equations, Steady-State and Time Dependent Problems. Siam, Philadelphia.
 Dingwell, R.D., Edmonds, M.J., Helm, J.L., Malomed, B.A. and Ohberg, P. (2018) Non-Integrable Dynamics of Matter-Wave Solitons in a Density-Dependent Gauag Theory. New Journal of Physics, 20, Article ID: 043004.