Travelling wave solutions and modulation instability analysis of the nonlinear Manakov-system

Abstract

In this paper, the accurate closed-form solutions of the Manakov-system are extracted via the extended $(\frac{G^{\prime }}{G^2})$-expansion method, the $\exp (-\varphi (\mathrm{\Omega }))$-expansion method and generalized Kudryashov method. The solutions are obtained in the form of trigonometric, hyperbolic and rational functions. The dynamical behaviour of obtained solutions is discussed by plotting 3D graphs, 2D contour graphs and density graphs. The symbolic softwares such as Maple and Mathematica are used to plot the graphs of retrieved solutions. The modulation instability (MI) of Manakov-system is also discussed in this paper.

Keywords:

• Manakov-system
• extended $(\frac{G^{\prime }}{G^2})$-expansion method
• $\exp (-\varphi (\mathrm{\Omega }))$-expansion method
• generalized Kudryashov method
• modulation instability
• exact solutions

1. Introduction

The Lorentz force law and a set of coupled partial differential equations known as Maxwell's equations form the basis of classical electromagnetism, classical optics and electric circuits. The coupled nonlinear Schrödinger equations are obtained by converting Maxwell's equations into cylindrical coordinates and accounting for an optical fibre's boundary conditions. When the inverse scattering transform is applied to the resultant equations, the Manakov-system is obtained.

Manakov-system [1,2] has been studied many researchers in recent years and novel results have been achieved. It can be written as (1) $\begin{array}{rl}& x_1{\vartheta }_{\mathrm{xx}}+\iota {\vartheta }_t+y_1(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\vartheta =0,\end{array}$(1) (2) $\begin{array}{rl}& x_2{\mathrm{\Lambda }}_{\mathrm{xx}}+\iota {\mathrm{\Lambda }}_t+y_2(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\mathrm{\Lambda }=0,\end{array}$(2) where $\vartheta (x,t)$ and $\mathrm{\Lambda }(x,t)$ are complex valued functions and $x_i$'s and $y_i$'s are constants.

The Manakov-system describes the wave propagation in two mode optical fibres and photo refractive materials [3]. Multicolour solitons are vector solitons that maintain their form throughout propagation and have two or more polarization components coupled to each other. Only coupled nonlinear Schrödinger equations can explain the behaviour of both vector and multicolour solitons. Manakov initially presented the vector soliton in 1974, when he deduced the one-soliton solution and performed an asymptotic analysis for the two-soliton solution using the inverse scattering transform (IST) approach. The auxiliary equation method is used to obtain the V-shaped, periodic singular, periodic bright, etc., soliton solutions of Manakov-system [2]. The physical interpretation of Manakov-system is confirmed in 1996 in the experimental observation of Manakov-solitons [4]. The Manakov-system is a special case of coupled nonlinear Schrödinger equation [5–7].

Many physical phenomena and dynamical processes are modelled by nonlinear evolution equations. The graphical demonstration and stability of NLEEs urged many researchers to find exact solutions using different exact methods. The multiple lump and rogue waves solutions of some nonlinear coupled equations are obtained by different exact methods. For example, the solution of nematic liquid crystals model is obtained by Kerr law [8], the Hirota bilinear approach is used to find the solution of time-fractional resonant nonlinear Schrödinger equation [9], the solution of Zabolotskaya–Khokholov model is obtained by new modified extended direct algebraic [10], the solution of $(3+1)$-dimensional generalized Kortewegde–VriesZakharovKuznetsov equation is obtained by new modified extended direct algebraic method [11], to obtain some chirped periodic waves with self-steeping phenomena we investigate the resonant nonlinear Schrödinger equation [12], the soliton solution of nonlinear telegraph equation is obtained by Hirota bilinear method [13], the solution of chlorite iodide malonic acid chemical equation is obtained by residual power series method [14], the chirped periodic wave type solutions of cubic-quintic nonlinear Schrödinger equation is obtained by Jacobi elliptic function [15], the solution of Ablowitz–Kaup–Newell–Segur water waves equation is obtained by painlevé analysis [16], the solution of nonlinear Biswas–Milovic coupled system is obtained in [17] and the solution of time-fractional ion sound and Langmuir waves system is obtained in [18].

The graphical demonstration and stability of Manakov-system urged many researchers to find its exact solutions using different exact methods. The exact solutions of Manakov-system are extracted using bilinear method [19], Darboux transformation method [20,21], extended simplest equation method [22], Hirota method [23,24], modified physics-informed neural network method [25], Hirota's bilinearization method [26], finite difference and finite element method [27], complex Toda chain model [28,29], inverse scattering transform [30], deep learning approach [31], Hamiltonian boundary value method [32], dressing technique [33], Riemann–Hilbert method [34], double-hump soliton solutions with the help of Hirota method [35], trial equation method [1], perturbed complex Toda chain [36], integrable decomposition of Manakov-system [37], polarization modulation instability in a Manakov fibre system [38], extended auxiliary equation method [39], generalization of Manakov equations [40].

The Manakov-system is investigated in the article by employing three most favourite, reliable and authentic integrating strategies. The three integrating techniques are the $\exp (-\varphi (\mathrm{\Omega }))$-expansion method, $(\frac{G^{\prime }}{G^2})$-expansion method and the generalized Kudryashov method. In order to utilize the proposed techniques, the nonlinear PDEs are first converted into nonlinear ODEs by applying travelling wave transformation. The proposed techniques efficiently extract trigonometric function solutions, hyperbolic function solutions and rational solutions. Trigonometric function solutions can be characterized as periodic solutions. Hyperbolic function solutions can be characterized as dark soliton, kink soliton, bright soliton, singular soliton and complexitons.

Integration schemes applied in this article are most efficient and reliable, but the three have their own limitations. In fact, it is important to state here that every analytical method for nonlinear evolution equations has to satisfy certain constraint conditions. The existence and validity of the obtained solutions depend on these constraint conditions. For example, the classic inverse scattering transform cannot be extended to the cases of power law linearity, dual-power law or even log-law nonlinearity. The method of semi-inverse variational principle is restricted to the retrieval of bright soliton solutions only. Otherwise, for dark or singular solitons, the corresponding stationary integral would be rendered divergent and that prevents universal applicability.

The proposed methods have been successfully utilized in various studies to extract the soliton solutions of some important nonlinear partial differential equations (NLPDEs) such as KDV equation [41], time-fractional parabolic equation [42], strain wave equation [43], Gerdjikov–Ivanov equation [44], Lakshmanan–Porsezian–Daniel model [45,46], Triki–Biswas equation [47], the Zakharov Kuznetsov–Benjamin Bona Mahony equation and the ill-posed Boussinesq equation [48], the nonlinear Schrödinger equation [49], the $(1+1)$-dimensional nonlinear dispersive modified Benjamin–Bona–Mahony and the seventh-order Sawada–Kotera–Ito equations [50], $(2+1)$-dimensional DLW equations and Maccari's equations [51], the $(1+1)$-dimensional generalized Broer–Kaup equation [52], the Korteweg–de Vries equation [53], PHI-four equation and the Fisher equation [54], Westervelt equation [55] and $(2+1)$-dimensional Biswas–Milovic equation [56].

This paper has seven sections. Section 2 presents demarcation of Manakov-system. Suggested are elaborated in Section 3. Suggested methods are applied to extract the solutions of Manakov-system in Section 4. The modulation instability of Manakov-system is explained in Section 5. The work is summarized in Section 6.

2. Demarcation of Manakov-system

Considering the following travelling wave transformations: (3) $\begin{array}{r}\begin{array}{l}\vartheta (x,t)=q_1(r_1)\exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ \mathrm{\Lambda }(x,t)=q_2(r_2)\exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}\}\end{array}$(3) where (4) $\begin{array}{r}{r}_i=x-{\rho }_{i}t\end{array}$(4) and (5) $\begin{array}{r}{\mathrm{\Delta }}_i(x,t)=-{\mathrm{\Psi }}_{i}x+{\mathrm{\Phi }}_{i}t+{ϝ}_i,\end{array}$(5) where $q_i(r_i)$ and ${\mathrm{\Delta }}_i(x,t)$ for i = 1, 2 are the amplitude and phase component, respectively. Substituting Equation (3(3) $\begin{array}{r}\begin{array}{l}\vartheta (x,t)=q_1(r_1)\exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ \mathrm{\Lambda }(x,t)=q_2(r_2)\exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}\}\end{array}$(3) ) into Equations (1(1) $\begin{array}{rl}& x_1{\vartheta }_{\mathrm{xx}}+\iota {\vartheta }_t+y_1(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\vartheta =0,\end{array}$(1) ) and (2(2) $\begin{array}{rl}& x_2{\mathrm{\Lambda }}_{\mathrm{xx}}+\iota {\mathrm{\Lambda }}_t+y_2(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\mathrm{\Lambda }=0,\end{array}$(2) ), following real and imaginary parts are achieved as

For imaginary part (6) $\begin{array}{r}{\rho }_i=-2x_i{\mathrm{\Psi }}_i.\end{array}$(6) Real part takes the following form: (7) $\begin{array}{rl}& y_i{q}_i^3+x_i{q}_i^{″}-(x_i{\mathrm{\Psi }}_i^2\\ & +{\mathrm{\Phi }}_i)q_i+y_i{q}_i{q}_{\tilde{i}}^2=0,i=1,2,\tilde{i}=-i+3.\end{array}$(7) By applying principle of balance, then (8) $\begin{array}{r}{q}_i=q_{\tilde{i}}.\end{array}$(8) Equation (7(7) $\begin{array}{rl}& y_i{q}_i^3+x_i{q}_i^{″}-(x_i{\mathrm{\Psi }}_i^2\\ & +{\mathrm{\Phi }}_i)q_i+y_i{q}_i{q}_{\tilde{i}}^2=0,i=1,2,\tilde{i}=-i+3.\end{array}$(7) ) is transformed as (9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9)

3. Description of methods

In order to utilize the proposed techniques, the nonlinear PDEs are first converted into nonlinear ODEs by applying travelling wave transformation.

3.1. The extended $(\frac{G^{\prime }}{G^2})$-expansion method

The exact solution of transformed ODE is considered as (10) $\begin{array}{r}q(\mathrm{\Omega })=d_0+\sum _{i=1}^N[d_i{\left(\frac{G^{\prime }}{G^2}\right)}^i+e_i{\left(\frac{G^{\prime }}{G^2}\right)}^{-i}],\end{array}$(10) where $G=G(\mathrm{\Omega })$ satisfies (11) $\begin{array}{r}{\left(\frac{G^{\prime }}{G^2}\right)}^{\prime }=\varpi +\chi {\left(\frac{G^{\prime }}{G^2}\right)}^2.\end{array}$(11) In Equations (10(10) $\begin{array}{r}q(\mathrm{\Omega })=d_0+\sum _{i=1}^N[d_i{\left(\frac{G^{\prime }}{G^2}\right)}^i+e_i{\left(\frac{G^{\prime }}{G^2}\right)}^{-i}],\end{array}$(10) ) and (11(11) $\begin{array}{r}{\left(\frac{G^{\prime }}{G^2}\right)}^{\prime }=\varpi +\chi {\left(\frac{G^{\prime }}{G^2}\right)}^2.\end{array}$(11) ), $d_0$, $d_i$'s, $e_i$'s are unknown constants and $\varpi \ne 1$ and $\chi \ne 0$ are integers. The formal solutions of Equation (11(11) $\begin{array}{r}{\left(\frac{G^{\prime }}{G^2}\right)}^{\prime }=\varpi +\chi {\left(\frac{G^{\prime }}{G^2}\right)}^2.\end{array}$(11) ) are given in [57].

3.2. The $\exp (-\varphi (\mathrm{\Omega }))$-expansion method

According to $\exp (-\varphi (\mathrm{\Omega }))$-expansion method, solution of ODE is considered as (12) $\begin{array}{r}q(\mathrm{\Omega })=\sum _{i=0}^N\left[v_i\exp (-\varphi (\mathrm{\Omega }))\right],\end{array}$(12) where $v_i$'s are constants to be determined. In Equation (12(12) $\begin{array}{r}q(\mathrm{\Omega })=\sum _{i=0}^N\left[v_i\exp (-\varphi (\mathrm{\Omega }))\right],\end{array}$(12) ), $\varphi (\mathrm{\Omega })$ satisfies the following ODE: (13) $\begin{array}{r}{\varphi }^{\prime }(\mathrm{\Omega })=\exp (-\varphi (\mathrm{\Omega }))+\varpi \exp (\varphi (\mathrm{\Omega }))+\chi ,\end{array}$(13) where ϖ and χ are arbitrary constants. The general solutions of Equation (13(13) $\begin{array}{r}{\varphi }^{\prime }(\mathrm{\Omega })=\exp (-\varphi (\mathrm{\Omega }))+\varpi \exp (\varphi (\mathrm{\Omega }))+\chi ,\end{array}$(13) ) are given in [58].

3.3. The generalized Kudryashov method

The following general solution of Kudryashov method is used to obtained the exact solution of transformed ODE (14) $\begin{array}{r}q(\mathrm{\Omega })=\frac{\sum _{i=0}^N[g_i{u}^i(\mathrm{\Omega })]}{\sum _{j=0}^M[h_j{u}^j(\mathrm{\Omega })]},\end{array}$(14) where $g_i(i=0,1,2,\dots ,N)$ and $h_j(j=0,1,2,\dots ,M)$ are constant such that $g_N\ne 0$ and $h_M\ne 0$. Consider the ODE (15) $\begin{array}{r}\frac{\mathrm{du}(\mathrm{\Omega })}{d\mathrm{\Omega }}=u^2(\mathrm{\Omega })-u(\mathrm{\Omega }),\end{array}$(15) where $u=u(\mathrm{\Omega })$ satisfied. The solution of Equation (15(15) $\begin{array}{r}\frac{\mathrm{du}(\mathrm{\Omega })}{d\mathrm{\Omega }}=u^2(\mathrm{\Omega })-u(\mathrm{\Omega }),\end{array}$(15) ) is considered as (16) $\begin{array}{r}u(\mathrm{\Omega })=\frac{1}{1+D{e}^{\mathrm{\Omega }}},\end{array}$(16) where D is an integrating function.

4. Determination of travelling wave solutions

This section provides novel solutions of the proposed model via suggested integrating techniques.

4.1. The extended $(\frac{G^{\prime }}{G^2})$-expansion method

The exact solution of Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ) is investigated in Equation (17(17) $\begin{array}{r}{q}_i(v_i)=d_0+[d_1{\left(\frac{G^{\prime }}{G^2}\right)}^1+e_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}],i=1,2,\end{array}$(17) ). We obtained N = 1 by using homogenous balance principle. The solution is considered in Equation (17(17) $\begin{array}{r}{q}_i(v_i)=d_0+[d_1{\left(\frac{G^{\prime }}{G^2}\right)}^1+e_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}],i=1,2,\end{array}$(17) ) as follows: (17) $\begin{array}{r}{q}_i(v_i)=d_0+[d_1{\left(\frac{G^{\prime }}{G^2}\right)}^1+e_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}],i=1,2,\end{array}$(17) where in Equation (17(17) $\begin{array}{r}{q}_i(v_i)=d_0+[d_1{\left(\frac{G^{\prime }}{G^2}\right)}^1+e_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}],i=1,2,\end{array}$(17) ) the unknown constants are $d_0$, $d_1$ and $e_1$. Equations (17(17) $\begin{array}{r}{q}_i(v_i)=d_0+[d_1{\left(\frac{G^{\prime }}{G^2}\right)}^1+e_1{\left(\frac{G^{\prime }}{G^2}\right)}^{-1}],i=1,2,\end{array}$(17) ) and (11(11) $\begin{array}{r}{\left(\frac{G^{\prime }}{G^2}\right)}^{\prime }=\varpi +\chi {\left(\frac{G^{\prime }}{G^2}\right)}^2.\end{array}$(11) ) are substituted in Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ). The following system is obtained by collecting all the coefficients of $(\frac{G^{\prime }}{G^2}{)}^i$ along with i = $-3,-2,-1,0,1,2,3$ and putting them equal to zero (18) $\begin{array}{r}\begin{array}{l}2{\varpi }^2{x}_i{e}_1+2y_i{e}_1^2=0,\\ 6d_0{y}_i{e}_1^2=0,\\ 2\mathrm{\chi \varpi }(x_i{e}_1)+6y_i{d}_0^2{e}_1+6y_i{e}_1^2{d}_1\\ -x_i{\mathrm{\Psi }}_i^2{e}_1-e_1{\mathrm{\Phi }}_i=0,\\ 2d_0^3{y}_i+12d_0{d}_1{e}_1{y}_i-x_i{d}_0{\mathrm{\Psi }}_i^2-d_0{\mathrm{\Phi }}_i=0,\\ 2\mathrm{\chi \varpi }(x_i{d}_1)+6y_i{d}_0^2{d}_1+6y_i{d}_1^2{e}_1\\ -x_i{\mathrm{\Psi }}_i^2{d}_1-d_1{\mathrm{\Phi }}_i=0,\\ 6d_0{d}_1^2{y}_i=0,\\ 2{\chi }^2{x}_i{d}_1+2y_i{d}_1^2=0.\end{array}\}\end{array}$(18) By solving these algebraic system of equations simultaneously using Mathematica, the following sets of solutions are obtained: $\begin{array}{rl}& \mathbf{S}\mathbf{1}\mathbf{:}\{d_0=0,d_1=0,e_1=\pm \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}},\\ & {\mathrm{\Phi }}_i=x_i(2\mathrm{\chi \varpi }-{\mathrm{\Psi }}_i^2)\phantom{\frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}}\},\\ & \mathbf{S}\mathbf{2}\mathbf{:}\{d_0=0,d_1=\pm \frac{\mathrm{\iota \chi }\sqrt{x_i}}{\sqrt{y_i}},e_1=0,\\ & {\mathrm{\Phi }}_i=x_i(2\mathrm{\varpi \chi }-{\mathrm{\Psi }}_i^2)\phantom{\frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}}\},\\ & \mathbf{S}\mathbf{3}\mathbf{:}\{d_0=0,d_1=\pm \frac{\mathrm{\iota \chi }\sqrt{x_i}}{\sqrt{y_i}},e_1=\mp \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}},\\ & {\mathrm{\Phi }}_i=x_i(8\mathrm{\chi \varpi }-{\mathrm{\Psi }}_i^2)\phantom{\frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}}\},\\ & \mathbf{S}\mathbf{4}\mathbf{:}\{d_0=0,d_1=\pm \frac{\mathrm{\iota \chi }\sqrt{x_i}}{\sqrt{y_i}},e_1=\pm \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}},\\ & {\mathrm{\Phi }}_i=-4x_i\mathrm{\chi \varpi }-x_i{\mathrm{\Psi }}_i^2\phantom{\frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}}\}.\end{array}$

Family 1

According to S1 the trigonometric, hyperbolic and rational functions solution are obtained.

If ($\mathrm{\varpi \chi }>0$), then following periodic solutions are obtained: $\begin{array}{rl}& {\vartheta }_a(x,t)\\ & =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}\\ & \times {\left[\sqrt{\frac{\varpi }{\chi }}\frac{n_1\cos (\sqrt{\mathrm{\varpi \chi }}{r}_1)+n_2\sin (\sqrt{\mathrm{\varpi \chi }}{r}_1)}{n_2\cos (\sqrt{\mathrm{\varpi \chi }}{r}_1)-n_1\sin (\sqrt{\mathrm{\varpi \chi }}{r}_1)}\right]}^{-1}\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_a(x,t)\\ & =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}\\ & \times {\left[\sqrt{\frac{\varpi }{\chi }}\frac{n_1\cos (\sqrt{\mathrm{\varpi \chi }}{r}_2)+n_2\sin (\sqrt{\mathrm{\varpi \chi }}{r}_2)}{n_2\cos (\sqrt{\mathrm{\varpi \chi }}{r}_2)-n_1\sin (\sqrt{\mathrm{\varpi \chi }}{r}_2)}\right]}^{-1}\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$The graphical simulation of $|{\vartheta }_a(x,t)|$ for choosing arbitrary values of parameters as $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$.

If ($\mathrm{\varpi \chi }<0$), then the following dark soliton solutions are obtained: $\begin{array}{rl}& {\vartheta }_b(x,t)\\ & =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}\\ & \times {[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]}^{-1}\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_b(x,t)\\ & =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}\\ & \times {[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]}^{-1}\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$The graphical simulation of $|{\vartheta }_b(x,t)|$ for $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$.

If $\varpi =0$ and $\chi \ne 0$, then $\begin{array}{rl}{\vartheta }_c(x,t)& =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}\\ & \times {(-\frac{n_1}{\chi (n_1{r}_1+n_2)})}^{-1}\exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_c(x,t)& =\pm \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}\\ & \times {(-\frac{n_1}{\chi (n_1{r}_2+n_2)})}^{-1}\exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}$where $\begin{array}{rl}& {\mathrm{\Delta }}_1(x,t)=-{\mathrm{\Psi }}_{1}x+{ϝ}_1+\mathrm{tx}(-{\mathrm{\Psi }}_1^2+2\mathrm{\varpi \chi }),\\ & {\mathrm{\Delta }}_2(x,t)=-{\mathrm{\Psi }}_{2}x+{ϝ}_2+\mathrm{tx}(-{\mathrm{\Psi }}_2^2+2\mathrm{\varpi \chi }).\end{array}$

Family 2

According to S2 the trigonometric, hyperbolic and rational functions solution are obtained.

If ($\mathrm{\varpi \chi }>0$), then the following periodic solutions are obtained: $\begin{array}{rl}& {\vartheta }_d(x,t)=\pm \frac{\mathrm{\iota \sigma }\sqrt{x_1}}{\sqrt{y_1}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_1)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_1)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}\right]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_d(x,t)=\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_2)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_2)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}\right]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$The graphical simulation of $|{\vartheta }_d(x,t)|$ for $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$.

If ($\mathrm{\varpi \chi }<0$), then following singular soliton solutions are obtained: $\begin{array}{rl}& {\vartheta }_e(x,t)\\ & =\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_e(x,t)\\ & =\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If $\varpi =0$ and $\chi \ne 0$, then following rational solutions are obtained: $\begin{array}{rl}& {\vartheta }_f(x,t)=\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}(-\frac{n_1}{\chi (n_1{r}_1+n_2)})\exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_f(x,t)=\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}(-\frac{n_1}{\chi (n_1{r}_2+n_2)})\exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}$where $\begin{array}{rl}& {\mathrm{\Delta }}_1(x,t)=-{\mathrm{\Psi }}_{1}x+{ϝ}_1+\mathrm{tx}(-{\mathrm{\Psi }}_1^2+2\mathrm{\varpi \chi }),\\ & {\mathrm{\Delta }}_2(x,t)=-{\mathrm{\Psi }}_{2}x+{ϝ}_2+\mathrm{tx}(-{\mathrm{\Psi }}_2^2+2\mathrm{\varpi \chi }).\end{array}$

Family 3

According to S3 the trigonometric, hyperbolic and rational functions solution are obtained. If($\mathrm{\varpi \chi }>0$), then following periodic solutions are obtained: $\begin{array}{rl}{\vartheta }_g(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_1)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_1)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}\right]\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}{\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_1)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_1)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}\right]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_g(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_2)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_2)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}\right]\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}{\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_2)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_2)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}\right]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)).\end{array}$If ($\mathrm{\varpi \chi }<0$), then following dark-singular combo soliton solutions are obtained: $\begin{array}{rl}& {\vartheta }_h(x,t)\\ & =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}{[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_h(x,t)\\ & =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}{[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If $\varpi =0$ and $\chi \ne 0$, then following rational solutions are obtained: $\begin{array}{rl}{\vartheta }_i(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}(-\frac{n_1}{\chi (n_1{r}_1+n_2)})\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}{(-\frac{n_1}{\chi (n_1{r}_1+n_2)})}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_i(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}(-\frac{n_1}{\chi (n_1{r}_2+n_2)})\\ & \mp \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}{(-\frac{n_1}{\chi (n_1{r}_2+n_2)})}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}$where $\begin{array}{rl}& {\mathrm{\Delta }}_1(x,t)=-{\mathrm{\Psi }}_{1}x+{ϝ}_1+\mathrm{tx}(-{\mathrm{\Psi }}_1^2+8\mathrm{\varpi \chi }),\\ & {\mathrm{\Delta }}_2(x,t)=-{\mathrm{\Psi }}_{2}x+{ϝ}_2+\mathrm{tx}(-{\mathrm{\Psi }}_2^2+8\mathrm{\varpi \chi }).\end{array}$

Family 4

According to S4 the trigonometric, hyperbolic and rational functions solution are obtained.

If($\mathrm{\varpi \chi }>0$), then following periodic solutions are obtained: $\begin{array}{rl}{\vartheta }_j(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_1)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_1)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}\right]\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}{\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_1)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_1)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_1)\end{array}}\right]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_j(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_2)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_2)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}\right]\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}{\left[\sqrt{\frac{\varpi }{\chi }}\frac{\begin{array}{c}{n}_1\cos (\sqrt{\varpi \chi }{r}_2)\\ +n_2\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}{\begin{array}{c}{n}_2\cos (\sqrt{\varpi \chi }{r}_2)\\ -n_1\sin (\sqrt{\varpi \chi }{r}_2)\end{array}}\right]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If ($\mathrm{\varpi \chi }<0$), then following dark-singular combo soliton solutions are obtained: $\begin{array}{rl}& {\vartheta }_k(x,t)\\ & =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_1}}{\sqrt{y_1}}{[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_1)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_1)-n_2\end{array}}]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_k(x,t)\\ & =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh \\ (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_2}}{\sqrt{y_2}}{[-\frac{\sqrt{|\varpi \chi |}}{\chi }\frac{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)+n_2\end{array}}{\begin{array}{c}{n}_1\cosh (2\sqrt{|\varpi \chi |}{r}_2)\\ +n_1\sinh \\ (2\sqrt{|\varpi \chi |}{r}_2)-n_2\end{array}}]}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If $\varpi =0$ and $\chi \ne 0$, then following rational solutions are obtained: $\begin{array}{rl}{\vartheta }_l(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_1}}{\sqrt{y_1}}(-\frac{n_1}{\chi (n_1{r}_1+n_2)})\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}{(-\frac{n_1}{\chi (n_1{r}_1+n_2)})}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_l(x,t)& =[\pm \frac{\mathrm{\iota \chi }\sqrt{x_2}}{\sqrt{y_2}}(-\frac{n_1}{\chi (n_1{r}_2+n_2)})\\ & \pm \frac{\mathrm{\iota \varpi }\sqrt{x_i}}{\sqrt{y_i}}{(-\frac{n_1}{\chi (n_1{r}_2+n_2)})}^{-1}]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}$where $\begin{array}{rl}& {\mathrm{\Delta }}_1(x,t)=-{\mathrm{\Psi }}_{1}x+{ϝ}_1-\mathrm{tx}({\mathrm{\Psi }}_1^2+4\mathrm{\varpi \chi }),\\ & {\mathrm{\Delta }}_2(x,t)=-{\mathrm{\Psi }}_{2}x+{ϝ}_2-\mathrm{tx}({\mathrm{\Psi }}_2^2+4\mathrm{\varpi \chi }).\end{array}$

4.2. The $\exp (-\varphi (\mathrm{\Omega }))$-expansion method

According to $\exp (-\varphi (\mathrm{\Omega }))$-expansion method, solution of Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ) is considered as (19) $\begin{array}{r}{q}_i({\mathrm{\Omega }}_i)=v_0+v_1\exp (-\mathrm{\Psi }({\mathrm{\Omega }}_i)),i=1,2,\end{array}$(19) where the unknown constants are $v_0$ and $v_1$. Equations (19(19) $\begin{array}{r}{q}_i({\mathrm{\Omega }}_i)=v_0+v_1\exp (-\mathrm{\Psi }({\mathrm{\Omega }}_i)),i=1,2,\end{array}$(19) ) and (13(13) $\begin{array}{r}{\varphi }^{\prime }(\mathrm{\Omega })=\exp (-\varphi (\mathrm{\Omega }))+\varpi \exp (\varphi (\mathrm{\Omega }))+\chi ,\end{array}$(13) ) are substituted in Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ). The following system of algebraic equations is obtained by collecting all the coefficients of $\exp (-\varphi (\mathrm{\Omega }){)}^i$ with i = −3, −2, −1, 0 and putting them equal to zero (20) $\begin{array}{r}\begin{array}{l}2x_i{v}_1+2y_i{v}_1^3=0,\\ 3x_i{v}_1\chi +6y_i{v}_0{v}_1^2=0,\\ {\chi }^2{x}_i{v}_1+2x_i{v}_1\varpi -x_i{v}_1{\mathrm{\Psi }}_i^2+6y_i{v}_0^2{r}_1-t_1{\mathrm{\Phi }}_i=0,\\ -v_0{\mathrm{\Phi }}_i+2y_i{v}_0^3-x_i{v}_0{\mathrm{\Psi }}_i^2+v_1{x}_i\mathrm{\varpi \chi }=0.\end{array}\}\end{array}$(20) The following set of solution is obtained by solving these algebraic equations using mathematica: $\begin{array}{rl}\mathbf{S}\mathbf{5}\mathbf{:}& \{v_1=\frac{2v_0}{\chi },x_i=-\frac{4y_i{v}_0^2}{{\chi }^2},\\ & {\mathrm{\Phi }}_i=\frac{2y_i({\chi }^2-4\varpi +2{\mathrm{\Psi }}_i^2)v_0^2}{{\chi }^2}\}.\end{array}$

Family 5

The following solutions of Manakov-system are obtained corresponding to S5.

If ${\chi }^2-4\varpi >0$ and $\varpi \ne 0$, then $\begin{array}{rl}& {\vartheta }_m(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{2\varpi }{\begin{array}{c}-\sqrt{{\chi }^2-4\varpi }\\ \tanh (\frac{\sqrt{{\chi }^2-4\varpi }({\mathrm{\Omega }}_1+c)}{2})-\chi \end{array}}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_m(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{2\varpi }{\begin{array}{c}-\sqrt{{\chi }^2-4\varpi }\\ \tanh (\frac{\sqrt{{\chi }^2-4\varpi }({\mathrm{\Omega }}_2+c)}{2})-\chi \end{array}}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$The graphical simulation of $|{\vartheta }_m(x,t)|$ for $v_0=1$, c = 0, $\varpi =1$, $\chi =3$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$ is represented in Figure 4.

If ${\chi }^2-4\varpi >0$ and $\varpi =0$, then $\begin{array}{rl}{\vartheta }_n(x,t)& =[v_0+\frac{2v_0}{\chi }(\frac{\chi }{\exp (\chi ({\mathrm{\Omega }}_1+c))-1}))]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_n(x,t)& =[v_0+\frac{2v_0}{\chi }(\frac{\chi }{\exp (\chi ({\mathrm{\Omega }}_2+c))-1}))]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If ${\chi }^2-4\varpi =0$ and $\chi \ne 0$ and $\varpi \ne 0$, then $\begin{array}{rl}{\vartheta }_o(x,t)& =[v_0-\frac{2v_0}{\chi }\left(\frac{{\chi }^2({\mathrm{\Omega }}_1+c)}{2(\chi ({\mathrm{\Omega }}_1+c)+2)}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ {\mathrm{\Lambda }}_o(x,t)& =[v_0-\frac{2v_0}{\chi }\left(\frac{{\chi }^2({\mathrm{\Omega }}_1+c)}{2(\chi ({\mathrm{\Omega }}_1+c)+2)}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$The graphical simulation of $|{\vartheta }_o(x,t)|$ for $v_0=1$, c = 1, $\varpi =1$, $\chi =2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$.

If ${\chi }^2-4\varpi =0$ and $\chi =0$ and $\varpi =0$, then $\begin{array}{rl}& {\vartheta }_p(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{1}{{\mathrm{\Omega }}_1+c}\right)]\exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_p(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{1}{{\mathrm{\Omega }}_2+c}\right)]\exp (\iota {\mathrm{\Delta }}_2(x,t)).\end{array}$If ${\chi }^2-4\varpi <0$, then $\begin{array}{rl}& {\vartheta }_s(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{2\varpi }{\begin{array}{c}\sqrt{4\varpi -{\chi }^2}\\ \tan (\frac{\sqrt{4\varpi -{\chi }^2}({\mathrm{\Omega }}_1+c)}{2})-\chi \end{array}}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_1(x,t)),\\ & {\mathrm{\Lambda }}_s(x,t)=[v_0+\frac{2v_0}{\chi }\left(\frac{2\varpi }{\begin{array}{c}\sqrt{4\varpi -{\chi }^2}\\ \tan (\frac{\sqrt{4\varpi -{\chi }^2}({\mathrm{\Omega }}_2+c)}{2})-\chi \end{array}}\right)]\\ & \times \exp (\iota {\mathrm{\Delta }}_2(x,t)),\end{array}$where $\begin{array}{rl}& {\mathrm{\Delta }}_1(x,t)=-{\mathrm{\Psi }}_{1}x+{ϝ}_1+\frac{2t{y}_1(2{\mathrm{\Psi }}_1^2-4\varpi +{\chi }^2)v_0^2}{{\chi }^2},\\ & {\mathrm{\Delta }}_2(x,t)=-{\mathrm{\Psi }}_{2}x+{ϝ}_2+\frac{2t{y}_2(2{\mathrm{\Psi }}_2^2-4\varpi +{\chi }^2)v_0^2}{{\chi }^2}.\end{array}$

4.3. The generalized Kudryashov method

In this section, Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ) is investigated by generalized Kudryashov method. Using homogenous balance principle we obtained N = M + 1. Choosing M = 1 we obtain N = 1. The general solution is assumed as (21) $\begin{array}{r}q(\mathrm{\Omega })=\frac{g_0+g_{1}u+g_2{u}^2}{h_0+h_{1}u},\end{array}$(21) where $g_0$, $g_1$, $g_2$, $h_0$ and $h_1$ are constant. Equation (21(21) $\begin{array}{r}q(\mathrm{\Omega })=\frac{g_0+g_{1}u+g_2{u}^2}{h_0+h_{1}u},\end{array}$(21) ) is substituted in Equation (9(9) $\begin{array}{r}{x}_i{q}_i^{″}+2y_i{q}_i^3-(x_i{\mathrm{\Psi }}_i^2+{\mathrm{\Phi }}_i)q_i=0.\end{array}$(9) ), we attain a polynomial in $u(\mathrm{\Omega })$. By collecting all the coefficient of $u(\mathrm{\Omega })$ of same power and putting them equal to zero, the following system is obtained: (22) $\begin{array}{r}\begin{array}{l}2g_2^3{y}_i+2g_2{h}_1^2{x}_i=0,\\ 6g_1{g}_2^2{y}_i+6g_2{h}_0{h}_1{x}_i-3g_2{h}_1^2{x}_i=0,\\ -{\mathrm{\Psi }}_i^2{g}_2{h}_1^2{x}_i-{\mathrm{\Phi }}_i{g}_2{h}_1^2+6g_0{g}_2^2{y}_i+6g_1^2{g}_2{y}_i\\ +6g_2{h}_0^2{x}_i-9g_2{h}_0{h}_1{x}_i+g_2{h}_1^2{x}_i=0,\\ -{\mathrm{\Psi }}_i^2{g}_1{h}_1^2{x}_i-2{\mathrm{\Psi }}_i^2{g}_2{h}_0{h}_1{x}_i-{\mathrm{\Phi }}_i{g}_1{h}_1^2\\ -2{\mathrm{\Phi }}_i{g}_2{h}_0{h}_1+12g_0{g}_1{g}_2{y}_i-2g_0{h}_0{h}_1{x}_i\\ -g_0{h}_1^2{x}_i+2g_1^3{y}_i+2g_1{h}_0^2{x}_i+g_1{h}_0{h}_1{x}_i\\ -10g_2{h}_0^2{x}_i+3g_2{h}_0{h}_1{x}_i=0,\\ -{\mathrm{\Psi }}_i^2{g}_0{h}_1^2{x}_i-2{\mathrm{\Psi }}_i^2{g}_1{h}_0{h}_1{x}_i-{\mathrm{\Psi }}_i^2{g}_2{h}_0^2{x}_i\\ -{\mathrm{\Phi }}_i{g}_0{h}_1^2-2{\mathrm{\Phi }}_i{g}_1{h}_0{h}_1-{\mathrm{\Phi }}_i{g}_2{h}_0^2\\ +6g_0^2{g}_2{y}_i+6g_0{g}_1^2{y}_i+3g_0{h}_0{h}_1{x}_i\\ +g_0{h}_1^2{x}_i-3g_1{h}_0^2{x}_i-g_1{h}_0{h}_1{x}_i\\ +4g_2{h}_0^2{x}_i=0,\\ -2{\mathrm{\Psi }}_i^2{g}_0{h}_0{h}_1{x}_i-{\mathrm{\Psi }}_i^2{g}_1{h}_0^2{x}_i\\ -2{\mathrm{\Phi }}_i{g}_0{h}_0{h}_1-{\mathrm{\Phi }}_i{g}_1{h}_0^2+6g_0^2{g}_1{y}_i\\ -g_0{h}_0{h}_1{x}_i+g_1{h}_0^2{x}_i=0,\\ -{\mathrm{\Psi }}_i^2{g}_0{h}_0^2{x}_i-{\mathrm{\Phi }}_i{g}_0{h}_0^2+2y_i{g}_0^3=0.\end{array}\}\end{array}$(22) The following set of solutions are obtained by solving above algebraic system of equations: $\begin{array}{rl}& \mathbf{S}\mathbf{6}\mathbf{:}\{h_0=h_0,h_1=h_1,\\ & g_0=\frac{1}{12}\frac{\begin{array}{c}(2{\mathrm{\Phi }}_i+x_i)x_i{h}_1^2-6h_0{h}_1{x}_i^2\\ -2x_i{h}_1^2{\mathrm{\Phi }}_i-h_1^2{x}_i^2\end{array}}{h_1{y}_i^2(-\frac{x_i}{y_i}{)}^{\frac{3}{2}}},\\ & g_1=-\frac{1}{2}\frac{x_i(2h_0-h_1)}{\sqrt{-\frac{x_i}{y_i}}{y}_i},g_2=\sqrt{-\frac{x_i}{y_i}}{h}_1,\\ & {\mathrm{\Psi }}_i=\sqrt{-\frac{1}{2}\frac{2{\mathrm{\Phi }}_i+x_i}{x_i}}\},\\ & \mathbf{S}\mathbf{7}\mathbf{:}\{h_0=-\frac{1}{2}{h}_1,h_1=h_1,g_0=0,\\ & g_1=\frac{h_1{x}_i}{\sqrt{-\frac{x_i}{y_i}}{y}_i},g_2=\sqrt{-\frac{x_i}{y_i}}{h}_1,{\mathrm{\Psi }}_i=\sqrt{-\frac{{\mathrm{\Phi }}_i-x_i}{x_i}}\},\\ & \mathbf{S}\mathbf{8}\mathbf{:}\{h_0=-\frac{1}{2}{h}_1,h_1=h_1,\\ & g_0=\frac{1}{6}\frac{h_1(({\mathrm{\Phi }}_i+2x_i)x_i-x_i{\mathrm{\Phi }}_i+x_i^2)}{y_i^2(-\frac{x_i}{y_i}{)}^{\frac{3}{2}}},\\ & g_1=\frac{x_i{h}_1}{\sqrt{-\frac{x_i}{y_i}}{y}_i},g_2=\sqrt{-\frac{x_i}{y_i}}{h}_1,{\mathrm{\Psi }}_i=\sqrt{-\frac{{\mathrm{\Phi }}_i+2x_i}{x_i}}\}.\end{array}$

Family 6

The hyperbolic function solution of Manakov-system is obtained corresponding to S6. (23) $\begin{array}{rl}{\vartheta }_1(x,t)& =\frac{1}{2}{\exp }^{\frac{1}{2}\iota (2t{\mathrm{\Phi }}_1-x\sqrt{-2-\frac{4{\mathrm{\Phi }}_1}{x_1}}+2{ϝ}_1)}\\ & \times \sqrt{-\frac{x_1}{y_1}}\tanh \left(\frac{t-x}{2}\right),\end{array}$(23) (24) $\begin{array}{rl}{\mathrm{\Lambda }}_1(x,t)& =\frac{1}{2}{\exp }^{\frac{1}{2}\iota (2t{\mathrm{\Phi }}_2-x\sqrt{-2-\frac{4{\mathrm{\Phi }}_2}{x_2}}+2{ϝ}_2)}\\ & \times \sqrt{-\frac{x_2}{y_2}}\tanh \left(\frac{t-x}{2}\right).\end{array}$(24)

Family 7

The hyperbolic function solution of Manakov-system is obtained corresponding to S7. (25) $\begin{array}{rl}& {\vartheta }_2(x,t)=-{\exp }^{\iota (t{\mathrm{\Phi }}_1-x\sqrt{1-\frac{{\mathrm{\Phi }}_1}{x_1}}+{ϝ}_1)}\mathrm{csch}(t-x)\sqrt{-\frac{x_1}{y_1}},\end{array}$(25) (26) $\begin{array}{rl}& {\mathrm{\Lambda }}_2(x,t)=-{\exp }^{\iota (t{\mathrm{\Phi }}_2-x\sqrt{1-\frac{{\mathrm{\Phi }}_2}{x_2}}+{ϝ}_2)}\mathrm{csch}(t-x)\sqrt{-\frac{x_2}{y_2}}.\end{array}$(26)

Family 8

The hyperbolic function solution of Manakov-system is obtained corresponding to S8. (27) $\begin{array}{rl}& {\vartheta }_3(x,t)={\exp }^{\iota (t{\mathrm{\Phi }}_1-x\sqrt{-2-\frac{{\mathrm{\Phi }}_1}{x_1}}+{ϝ}_1)}\coth (t-x)\sqrt{-\frac{x_1}{y_1}},\end{array}$(27) (28) $\begin{array}{rl}& {\mathrm{\Lambda }}_3(x,t)={\exp }^{\iota (t{\mathrm{\Phi }}_2-x\sqrt{-2-\frac{{\mathrm{\Phi }}_2}{x_2}}+{ϝ}_2)}\coth (t-x)\sqrt{-\frac{x_2}{y_2}}.\end{array}$(28)

5. Graphical representation

The presented graphs show the diverse wave behaviour corresponding to the Manakov-system. Figure 1 represents the periodic solution ${\vartheta }_a(x,t)$. Periodic solutions are the solutions expressed in terms of trigonometric functions. The graphical simulations in Figure 2 represent a dark soliton. Physically, dark soliton is a localized surface soliton that causes a temporary decrease in an associated wave amplitude. The periodic solution ${\vartheta }_d(x,t)$ is represented in Figure 3. A bright soliton is a localized surface soliton that causes a temporary increase in an associated wave amplitude (Figure 4). Thus, Figure 5 shows a dark-bright soliton. The graphical illustrations presented in Figure 6 exhibit a kink soliton. Figures 7 and 8 represent singular solitons.

Figure 1. The graphical simulation of $|{\vartheta }_a(x,t)|$ for $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 2. The graphical simulation of $|{\vartheta }_b(x,t)|$ for $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 3. The graphical simulation of $|{\vartheta }_d(x,t)|$ for $n_1=1$, $n_2=1$, $\varpi =1$, $\chi =1$, $x_1=2$, $y_1=-2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 4. The graphical simulation of $|{\vartheta }_m(x,t)|$ for $v_0=1$, c = 0, $\varpi =1$, $\chi =3$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 5. The graphical simulation of $|{\vartheta }_o(x,t)|$ for $v_0=1$, c = 1, $\varpi =1$, $\chi =2$, ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 6. The graphical simulation of ${\vartheta }_1(x,t)$ for $h_0=1$, $h_1=1$, c = 1, k = 1, $x_1=-4$, $y_1=2$ ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 7. The graphical simulation of ${\vartheta }_2(x,t)$ for $h_0=1$, $h_1=1$, c = 1, k = 1, $x_1=-4$, $y_1=2$ ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

Figure 8. The graphical simulation of ${\vartheta }_3(x,t)$ for $h_0=1$, $h_1=1$, c = 1, k = 1, $x_1=-4$, $y_1=2$ ${\mathrm{\Psi }}_1=1$, ${ϝ}_1=1$. (a) Surface graph. (b) Contour graph. (c) Density graph.

6. Modulation instability

The modulation instability (MI) of Equations (1(1) $\begin{array}{rl}& x_1{\vartheta }_{\mathrm{xx}}+\iota {\vartheta }_t+y_1(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\vartheta =0,\end{array}$(1) ) and (2(2) $\begin{array}{rl}& x_2{\mathrm{\Lambda }}_{\mathrm{xx}}+\iota {\mathrm{\Lambda }}_t+y_2(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\mathrm{\Lambda }=0,\end{array}$(2) ) is finalized in this section on the basis of linear stability technique. The steady-state solution of Equations (1(1) $\begin{array}{rl}& x_1{\vartheta }_{\mathrm{xx}}+\iota {\vartheta }_t+y_1(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\vartheta =0,\end{array}$(1) ) and (2(2) $\begin{array}{rl}& x_2{\mathrm{\Lambda }}_{\mathrm{xx}}+\iota {\mathrm{\Lambda }}_t+y_2(|\vartheta {|}^2+|\mathrm{\Lambda }{|}^2)\mathrm{\Lambda }=0,\end{array}$(2) ) is given as (29) $\begin{array}{rl}& \vartheta (x,t)=(p_1+M(x,t))e^{\iota (p_{1}x)},\end{array}$(29) (30) $\begin{array}{rl}& \mathrm{\Lambda }(x,t)=(p_2+N(x,t))e^{\iota (p_{2}x)},\end{array}$(30) where $M(x,t)$ and $N(x,t)$ represents the perturbed terms and the ratio of power incident is $f=\frac{p_1}{p_2}$. Inserting Equations (29(29) $\begin{array}{rl}& \vartheta (x,t)=(p_1+M(x,t))e^{\iota (p_{1}x)},\end{array}$(29) ) and (30(30) $\begin{array}{rl}& \mathrm{\Lambda }(x,t)=(p_2+N(x,t))e^{\iota (p_{2}x)},\end{array}$(30) ) to the Manakov-system and linearizing the result, the following equation are obtained: (31) $\begin{array}{rl}& x_1{M}_{\mathrm{xx}}+y_1{p}_1{p}_2(N+N^{\ast })+\iota M_t+y_1{p}_1^2(2M+M^{\ast })=0,\end{array}$(31) (32) $\begin{array}{rl}& x_2{N}_{\mathrm{xx}}+y_2{p}_1{p}_2(M+M^{\ast })+\iota N_t+y_2{p}_2^2(2N+N^{\ast })=0,\end{array}$(32) where $M(x,t{)}^{\ast }$ and $N(x,t{)}^{\ast }$ are complex conjugate of $M(x,t)$ and $N(x,t)$, respectively. The wave transformations are applied as (33) $\begin{array}{rl}& M(x,t)=H_1{e}^{-\iota (\mathrm{kx}-\mathrm{wt})}+G_1{e}^{\iota (\mathrm{kx}-\mathrm{wt})},\end{array}$(33) (34) $\begin{array}{rl}& M(x,t)=H_2{e}^{-\iota (\mathrm{kx}-\mathrm{wt})}+G_2{e}^{\iota (\mathrm{kx}-\mathrm{wt})}.\end{array}$(34) By substituting Equations (33(33) $\begin{array}{rl}& M(x,t)=H_1{e}^{-\iota (\mathrm{kx}-\mathrm{wt})}+G_1{e}^{\iota (\mathrm{kx}-\mathrm{wt})},\end{array}$(33) ) and (34(34) $\begin{array}{rl}& M(x,t)=H_2{e}^{-\iota (\mathrm{kx}-\mathrm{wt})}+G_2{e}^{\iota (\mathrm{kx}-\mathrm{wt})}.\end{array}$(34) ) into Equations (31(31) $\begin{array}{rl}& x_1{M}_{\mathrm{xx}}+y_1{p}_1{p}_2(N+N^{\ast })+\iota M_t+y_1{p}_1^2(2M+M^{\ast })=0,\end{array}$(31) ) and (32(32) $\begin{array}{rl}& x_2{N}_{\mathrm{xx}}+y_2{p}_1{p}_2(M+M^{\ast })+\iota N_t+y_2{p}_2^2(2N+N^{\ast })=0,\end{array}$(32) ), the following system of equations is obtained: (35) $\begin{array}{rl}& w{G}_1-x_1{k}^2{G}_1+2y_1{f}^2{p}^2{G}_2+y_{1}f{p}^2{G}_2\\ & +y_1{f}^2{p}^2{H}_1+y_{1}f{p}^2{H}_2=0,\end{array}$(35) (36) $\begin{array}{rl}& y_1{f}^2{p}^2{G}_1+y_{1}f{p}^2{G}_2-w{H}_1-x_1{k}^2{H}_1\\ & +2y_1{f}^2{p}^2{H}_1+y_{1}f{p}^2{H}_2=0,\end{array}$(36) (37) $\begin{array}{rl}& y_{2}f{p}^2{G}_1+w{G}_2-x_2{k}^2{G}_2+2y_2{p}^2{G}_2\\ & +y_{2}f{p}^2{H}_1+y_2{p}^2{H}_2=0,\end{array}$(37) (38) $\begin{array}{rl}& y_{2}f{p}^2{G}_1+y_2{p}^2{G}_2+y_{2}f{p}^2{H}_1-w{H}_2\\ & -x_2{k}^2{H}_2+2y_2{p}^2{H}_2=0.\end{array}$(38) The $4\times 4$ matrix is obtained by arranging the system of equations. This matrix is solved for unknowns $H_1$, $H_2$, $G_1$, $G_2$. The result of the determinant of matrix is given as follows: (39) $\begin{array}{rl}& -w^4+k^4{w}^2{x}_1^2+k^4{w}^2{x}_2^2-k^8{x}_1^2{x}_2^2-4f^2{k}^2{p}^2{w}^2{x}_1{y}_1\\ & +4f^2{k}^6{p}^2{x}_1{x}_2^2{y}_1-f^2{w}^2{p}^4{y}_1^2\\ & +4f^4{w}^2{p}^4{y}_1^2+k^4{f}^2{p}^4{x}_2^2{y}_1^2\\ & -4k^4{f}^4{p}^4{x}_2^2{y}_1^2-4k^2{p}^2{w}^2{x}_2{y}_2+4k^6{p}^2{x}_1^2{x}_2{y}_2\\ & +4f{k}^4{p}^4{x}_1{x}_2{y}_1{y}_2-16f^2{k}^4{p}^4{x}_1{x}_2{y}_1{y}_2\\ & +2f{k}^2{p}^6{x}_2{y}_1^2{y}_2-6f^3{k}^2{p}^6{x}_2{y}_1^2{y}_2-4f^2{k}^2{p}^6{x}_2{y}_1^2{y}_2\\ & +16f^4{k}^2{p}^2{x}_1{y}_1{y}_2^2+12f^2{k}^2{p}^6{x}_1{y}_1{y}_2^2-2f{p}^8{y}_1^2{y}_2^2\\ & +6f^3{p}^8{y}_1^2{y}_2^2+3f^2{p}^8{y}_1^2{y}_2^2\\ & -12f^4{p}^8{y}_1^2{k}^4{w}^2{y}_2^2=0.\end{array}$(39) the instability occurred and grows exponentially when the wave number is imaginary. Consequently, we must have the following condition for modulation instability to exist: (40) $\begin{array}{rl}& k^4(x_1^2+x_2^2)-4k^2{p}^2(f^2{x}_1{y}_1+x_2{y}_2)\\ & +p^4(f^2(-1+4f^2)y_1^2+3y_2^2)\\ & -((k^4(x_1^2-x_2^2)-4f^2{k}^2{p}^2{x}_1{y}_1\\ & +f^2(-1+4f^2)p^4{y}_1^2{)}^2+8k^2{x}_2(k^4{p}^2\\ & (x_1^2-x_2^2)+2k^2(f{p}^4-2f^2{p}^4)x_1{y}_1\\ & +p^2((1-3f^2)f{p}^4+f^2(-1+4f^2)\\ & p^4)y_1^2)y_2+2p^2(k^4{p}^2(-3x_1^2+11x_2^2)\\ & +4k^2(-2f{p}^4+3f^2{p}^4)x_1{y}_1+p^2\\ & (4(-1+3f^2)f{p}^4+p^8{y}_2^4+3f^2(1-4f^2{p}^4)y_1^2)y_2^2\\ & -{24k^2{p}^6{x}_2{y}_2^3)}^{\frac{1}{2}})<0.\end{array}$(40) The MI gain spectrum is retrieved as (41) $\begin{array}{rl}g(k)& =2\mathrm{Im}(w)\\ & =2\sqrt{\begin{array}{c}(p^4(f^2(-1+4f^2)y_1^2+3y_2^2)+k^4\\ (x_1^2+x_2^2)-4k^2{p}^2(f^2{x}_1{y}_1+x_2{y}_2)-\sqrt{S})\end{array}},\end{array}$(41) where (42) $\begin{array}{rl}S& =((k^4(x_1^2-x_2^2)-4f^2{k}^2{p}^2{x}_1{y}_1\\ & +f^2(-1+4f^2)p^4{y}_1^2{)}^2\\ & +2k^2(f{p}^4-2f^2{p}^4)\\ & \times x_1{y}_1+8k^2{x}_2(k^4{p}^2(x_1^2-x_2^2)+p^2((1-3f^2)f{p}^4\\ & +f^2(-1+4f^2)p^4)y_1^2)\\ & \times y_2+2p^2(k^4{p}^2(-3x_1^2+11x_2^2)\\ & +4k^2(-2f{p}^4+3f^2{p}^4)x_1{y}_1+p^2(4(-1+3f^2)f{p}^4\\ & +3f^2(1-4f^2{p}^4)y_1^2)y_2^2\\ & -24k^2{p}^6{x}_2{y}_2^3+9p^8{y}_2^4))).\end{array}$(42) It can be noticed that the modulation instability gain spectrum is dependent on group velocity dispersion and incident power (Figures 9 and 10) The MI gain is plotted in Figures 9–16 with different parameter values. The results are increasing curve by decreasing the values of p and f and vice versa. The results are decreasing curve if the value of p increases and the value of f decreases. Similarly, the results are increasing curve if the value of p decreases and the value of f increases.

Figure 9. Increasing value of f gives decreasing curve.

Figure 10. Decreasing values of f results in increasing curve.

Figure 11. Increasing in p results in decreasing curve.

Figure 12. Decrease in p results in increasing curve.

Figure 13. Increasing in both p and f results in decreasing curve.

Figure 14. Decreasing in both p and f both results in increasing curve.

Figure 15. Increase in p and decrease in f result in decreasing curve.

Figure 16. Decrease in p and increase in f result in increasing curve.

7. Conclusion

The three proposed methods namely, the extended $(\frac{G^{\prime }}{G^2})$-expansion method, the $\exp (-\varphi (\mathrm{\Omega }))$-expansion method and the generalized Kudryashov method are used to extract the travelling wave solutions of Manakov-system. Symbolic software, specifically Maple and Mathematica, is used to plot the 3D graphs, density graphs and 2D contour graphs of certain derived solutions. The obtained exact solutions are in the form of dark soliton, dark-singular combo soliton, singular soliton, Kink and anti-Kink soliton, periodic solutions and rational solutions. It is important to mention here that the solutions obtained in this paper by the application of three state of the art integrating techniques are novel and reliable. The dark-singular combo soliton solutions, Kink soliton and rational solutions have been retrieved for the first time in this article as compared with [1,2]. Moreover, the modulation instability analysis has been carried out for the very first time in this paper as compared with [1,2]. On the basis of linear stability analysis, the MI analysis of Manakov-system is investigated. To identify the region for the occurrence of instability, a dispersion relation is derived. The dispersion relation reflects stability of steady-state nature. The unstable solution arises when the wave number is complex, because the perturbation grows exponentially. Moreover, the stable solution occurs when the wave number is real. The MI gain is significantly affected by the incident power. Due to novelty of results, this article becomes a good contribution in the theory of soliton in future.

Disclosure statement

No potential conflict of interest was reported by the authors.