Optical Solitons and Conservation Laws for the Concatenation Model: Undetermined Coefficients and Multipliers Approach

Abstract

This paper retrieves an optical 1–soliton solution to a model that is written as a concatenation of the Lakshmanan–Porsezian–Daniel model and Sasa–Satsuma equation. The method of undetermined coefficients obtains a full spectrum of 1–soliton solutions. The multiplier approach yields the conserved densities, which subsequently lead to the conserved quantities from the bright 1–soliton solution.

1. Introduction

The concept and applications of optical solitons has been strong and steady ever since its inception in 1973. For the last five decades, the concept has paved the way for various advanced concepts, and many performance enhancements have been achieved. This gave way to the concepts of non-Kerr law solitons, dispersion-managed solitons, quiescent solitons, cubic–quartic solitons, highly dispersive solitons, various new forms of self-phase modulations and many others. This emerged into an ocean of novel concepts that is consistently being studied and analyzed [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. It is worth mentioning that topological solitons are a very wide class of solutions to non-linear Klein–Gordon equations (KGE) [21,22,23,24,25]. KGE appears in the study of relativistic Quantum Mechanics, while the study model in the current paper appears in the study of Quantum Optics and Optoelectronics.

Another form of new concept that has lately been proposed is the concatenation model formulated of a combination of two familiar models. These are the Lakshmanan–Porsezian–Daniel (LPD) model [1] and the Sasa–Satsuma equation (SSE) [13], while the linear temporal evolution, chromatic dispersion and the Kerr law of nonlinear refractive index change is the same as the familiar, nonlinear Schrödinger’s equation (NLSE). This model was first proposed in 2014 [2,3]. Later, it was revisited by Triki et al. during 2022 [16]. A Painleve analysis of this model was carried out very recently by Kudryashov et al. [6]. The current paper again considers the model to retrieve a full spectrum of 1–soliton solutions using the technique of undetermined coefficients. This is followed by the location of the conservation laws using the multipliers method. The soliton solutions, together with the parameter constraints for their existence, as well as the conservation laws, are displayed after a brief and concise introduction to the model.

Governing Model

The concatenation model given by (1) is true to its name. For $c_1=0$, (1) reduces to the familiar SSE, while for $c_2=0$, Equation (1) reduces to LPD model. However, for $c_1=c_2=0$, (1) collapses to the familiar NLSE. The bright 1–soliton solution, obtained using the indeterminate coefficients approach, will reveal the conserved quantities. The soliton existence criteria will also be guaranteed by the several constraint conditions that will naturally emerge from the scheme.

$$\begin{array}{ccc}\hfill i{q}_t+a{q}_{xx}+b{\left|q\right|}^{2}q& +& c_1\left[{\sigma }_1{q}_{xxxx}+{\sigma }_2{\left(q_x\right)}^2{q}^{*}+{\sigma }_3{\left|q_x\right|}^{2}q+{\sigma }_4{\left|q\right|}^2{q}_{xx}+{\sigma }_5{q}^2{q}_{xx}^{*}+{\sigma }_6{\left|q\right|}^{4}q\right]\hfill \\ & +& i{c}_2\left[{\sigma }_7{q}_{xxx}+{\sigma }_8{\left|q\right|}^2{q}_x+{\sigma }_9{q}^2{q}_x^{*}\right]=0.\hfill \end{array}$$

(1)

This model is an extended version of the familiar NLSE, SSE and LPD models. Thus, the current model describes the propagation of solitons through an optical fiber in an accurate manner. The presence of CD and self–phase modulation provides the basic ingredients needed for the solitons to exist. In case CD runs low, this is compensated by the third-order dispersion and fourth-order dispersion effects that stem from the SSE and LPD components of the concatenation model. The effect of the soliton self-frequency shift is also provided by the SSE component. The effect of nonlinear dispersion as well as the two-photon absorption effect are emerging from the LPD component of the model. Thus, this model carries an “universal” structure that can more accurately describe the soliton propagation dynamics through optical fibers and other form of waveguides. The details of the derivation of the soliton solutions and its conservation laws are enumerated and displayed in the subsequent sections.

2. Undetermined Coefficients

To construct soliton solutions for the system (1), we allow for it to have solutions in the form

$$q(x,t)=P(x,t)e^{i\Psi (x,t)}=P\left(\xi \right)e^{i\Psi (x,t)}$$

(2)

for which the arguments of the amplitude factors and the phase factor, respectively, are:

$$\xi (x,t)=x-vt,\mathrm{and}\Psi (x,t)=-\kappa x+\omega t+{\theta }_0,$$

(3)

where $P(x,t)$ denotes the waveform according to each type of nonlinear wave, $\kappa$ represents the soliton frequency, while $\omega$ and ${\theta }_0$ represent the wave number and phase constant, respectively. By substituting (2) into (1) and separating it into real and imaginary components, one obtains, for the real component,

$$\begin{array}{ccc}& & \left(c_1{\sigma }_1{\kappa }^4-c_2{\sigma }_7{\kappa }^3-a{\kappa }^2-\omega \right)P+\left\{b-c_1\left({\sigma }_2-{\sigma }_3+{\sigma }_4+{\sigma }_5\right){\kappa }^2+c_2\left({\sigma }_8-{\sigma }_9\right)\kappa \right\}{P}^3\hfill \\ & +& c_1{\sigma }_6{P}^5+\left(a-6c_1{\sigma }_1{\kappa }^2+3c_2{\sigma }_7\kappa \right)\frac{{\partial }^{2}P}{\partial x^2}+c_1{\sigma }_1\frac{{\partial }^{4}P}{\partial x^4}\hfill \\ & +& c_1\left({\sigma }_4+{\sigma }_5\right)P^2\frac{{\partial }^{2}P}{\partial x^2}+c_1\left({\sigma }_2+{\sigma }_3\right)P{\left(\frac{\partial P}{\partial x}\right)}^2=0.\hfill \end{array}$$

(4)

while the imaginary counterpart was

$$\begin{array}{ccc}\hfill & & \frac{\partial P}{\partial t}-\left(2a\kappa +3c_2{\sigma }_7{\kappa }^2-4c_1{\sigma }_1{\kappa }^3\right)\frac{\partial P}{\partial x}\hfill \\ & +& \left\{c_2\left({\sigma }_8+{\sigma }_9\right)-2\kappa c_1\left({\sigma }_2+{\sigma }_4-{\sigma }_5\right)\right\}{P}^2\frac{\partial P}{\partial x}+\left(c_2{\sigma }_7-4\kappa c_1{\sigma }_1\right)\frac{{\partial }^{3}P}{\partial x^3}=0.\hfill \end{array}$$

(5)

From the imaginary part, the soliton speed is,

$$v=-2\kappa \left(a+4c_1{\sigma }_1{\kappa }^2\right)$$

(6)

whenever

$$c_2\left({\sigma }_8+{\sigma }_9\right)=2\kappa c_1\left({\sigma }_2+{\sigma }_4-{\sigma }_5\right),$$

(7)

and

$$c_2{\sigma }_7=4\kappa c_1{\sigma }_1$$

(8)

hold. In view of the above two conditions (7) and (8), the real part of the equation reduces to

$$\begin{array}{ccc}\hfill & -& \left(3c_1{\sigma }_1{\kappa }^4+a{\kappa }^2+\omega \right)P+\left\{b-c_1\left({\sigma }_2-{\sigma }_3+{\sigma }_4+{\sigma }_5\right){\kappa }^2+c_2\left({\sigma }_8-{\sigma }_9\right)\kappa \right\}{P}^3+c_1{\sigma }_6{P}^5\hfill \\ & +& \left(a+6c_1{\sigma }_1{\kappa }^2\right)\frac{{\partial }^{2}P}{\partial x^2}+c_1{\sigma }_1\frac{{\partial }^{4}P}{\partial x^4}+c_1\left({\sigma }_4+{\sigma }_5\right)P^2\frac{{\partial }^{2}P}{\partial x^2}+c_1\left({\sigma }_2+{\sigma }_3\right)P{\left(\frac{\partial P}{\partial x}\right)}^2=0.\hfill \end{array}$$

(9)

In the following four subsections, the last real component Equation (9) will be employed to generate four distinct soliton solutions for the concatenation model (1).

2.1. Bright Soliton

In this subsection, a bright soliton solution is constructed with the assistance of the balancing principle. First, we assume the following form for the wave profile:

$$P(x,t)=A{\mathrm{sech}}^p\tau ,\tau =B(x-vt),$$

(10)

where A is the soliton’s amplitude, p is a parameter to be determined, and B is the corresponding inverse width, with v representing the soliton speed. The substitution of (10) into (9) leads to

$$\begin{array}{ccc}\hfill & & \left[p^2{Z}_3{B}^2+c_1{\sigma }_1{p}^4{B}^4-\left(\omega +a{\kappa }^2+3c_1{\sigma }_1{\kappa }^4\right)\right]A{sech}^p\tau \hfill \\ & & +\left[Z_2+p^2(Z_4+Z_5)B^2\right]A^3{sech}^{3p}\tau +c_1{\sigma }_6{A}^5{sech}^{5p}\tau \hfill \\ & & -p(1+p)\left[Z_3+2\left\{2+p(2+p)\right\}{c}_1{\sigma }_1{B}^2\right]A{B}^2{sech}^{p+2}\tau \hfill \\ & & +p(1+p)(2+p)(3+p)c_1{\sigma }_{1}A{B}^4{sech}^{p+4}\tau \hfill \\ & & -p\left[(1+p)Z_4+p{Z}_5\right]A^3{B}^2{sech}^{3p+2}\tau =0\hfill \end{array}$$

(11)

where, for convenience, we adopted the following notation:

$$\begin{array}{ccc}& & Z_2=b-c_1{\kappa }^2\left({\sigma }_2-{\sigma }_3+{\sigma }_4+{\sigma }_5\right)+\kappa c_2\left({\sigma }_8-{\sigma }_9\right)\hfill \end{array}$$

(12)

$$\begin{array}{ccc}& & Z_3=a+6c_1{\sigma }_1{\kappa }^2\hfill \end{array}$$

(13)

$$\begin{array}{ccc}& & Z_4=c_1\left({\sigma }_4+{\sigma }_5\right)\hfill \end{array}$$

(14)

$$\begin{array}{ccc}& & Z_5=c_1\left({\sigma }_2+{\sigma }_3\right)\hfill \end{array}$$

(15)

The balance between nonlinearity and dispersion leads to equating $3p$ and $p+2$, from which the value of the unknown parameter p becomes

$$p=1.$$

(16)

Substituting (16) into (11), by setting the coefficients of the linearly independent set of functions ${sech}^j\left[\tau \right]$ for $j=1,3,5$ to zero, we obtain:

$$\omega =-a{\kappa }^2+Z_3{B}^2-c_1{\sigma }_1\left(3{\kappa }^4-B^4\right)$$

(17)

and

$$A=\sqrt{\frac{2\left(Z_3+10c_1{\sigma }_1{B}^2\right)B^2}{Z_2+\left(Z_4+Z_5\right)B^2}},$$

(18)

subject to

$$\begin{array}{ccc}\hfill & & \left\{a+2c_1{\sigma }_1\left(3{\kappa }^2+5B^2\right)\right\}\hfill \\ & \times & \left\{b-c_1{\kappa }^2\left({\sigma }_2-{\sigma }_3+{\sigma }_4+{\sigma }_5\right)+c_2\kappa \left({\sigma }_8-{\sigma }_9\right)+c_1\left({\sigma }_2+{\sigma }_3+{\sigma }_4+{\sigma }_5\right)B^2\right\}>0\hfill \end{array}$$

(19)

where the soliton width can be written as

$$\begin{array}{ccc}\hfill B& =& \frac{1}{2}\left[\frac{Z_3\left\{\left(Z_4+Z_5\right)\left(2Z_4+Z_5\right)-40c_1^2{\sigma }_1{\sigma }_6\right\}-2c_1{\sigma }_1\left(2Z_4+7Z_5\right)Z_2}{c_1{\sigma }_1\left(100c_1^2{\sigma }_1{\sigma }_6-4Z_4^2-3Z_4{Z}_5+Z_5^2\right)}\right.\hfill \\ & \pm & {\left.\frac{\sqrt{{\left\{Z_3\left(Z_4+Z_5\right)-10c_1{\sigma }_1{Z}_2\right\}}^2\left\{{\left(2Z_4+Z_5\right)}^2-96c_1^2{\sigma }_1{\sigma }_6\right\}}}{c_1{\sigma }_1\left(100c_1^2{\sigma }_1{\sigma }_6-4Z_4^2-3Z_4{Z}_5+Z_5^2\right)}\right]}^{\frac{1}{2}}\hfill \end{array}$$

(20)

whenever the conditions

$${\sigma }_1{c}_1^3\left[100{\sigma }_1{\sigma }_6+\left({\sigma }_2+{\sigma }_3+{\sigma }_4+{\sigma }_5\right)\left\{{\sigma }_2+{\sigma }_3-4\left({\sigma }_4+{\sigma }_5\right)\right\}\right]\ne 0$$

(21)

and

$${\left\{2c_1\left({\sigma }_2+{\sigma }_3+{\sigma }_4+{\sigma }_5\right)\right\}}^2>96c_1^2{\sigma }_1{\sigma }_6$$

(22)

hold. Therefore, a bright soliton solution for the system (1) is

$$\begin{array}{c}\hfill q(x,t)=A\mathrm{sech}\left[B(x-vt)\right]e^{i\left(\kappa x+\omega t+{\theta }_0\right)}\end{array}$$

(23)

where the soliton speed v was computed early on (6), the amplitude A and wave numbers $\omega$ are expressed on (18) and (17) in terms of the soliton width, while the width B is given in (20) with corresponding constraints (21) and (22). The other solvability conditions needed for the bright soliton to exist are (7), (8), and the inequality given in (19).

Figure 1 shows the profile of a bright soliton, where the parameters are chosen as: $a=0.3$, $b=0.8$, $c_1=1.2$, $c_2=1.1$, ${\sigma }_1=0.5$, ${\sigma }_2=1.0$, ${\sigma }_3=2.1$, ${\sigma }_4=2.0$, ${\sigma }_5=-1.0$, ${\sigma }_6=3.2$, ${\sigma }_7=0.1$, ${\sigma }_8=0.4$, ${\sigma }_9=0.2$.

2.2. Dark Soliton

Then, we created a dark soliton solution following a similar procedure to the one used to construct the bright solution on the previous part. First, we assume the following form for the wave profile:

$$P(x,t)=A{\tanh }^p\tau ,\tau =B(x-vt),$$

(24)

where, in this case, A and B are free parameters for the soliton, while p is a parameter whose value will be calculated. The rest of the parameters have the same meaning as for the bright soliton. Now, substituting (24) into (9), and after adopting the notation (12)–(15), the real part of the Equation (9) emerges as:

$$\begin{array}{ccc}\hfill & & -\left[3c_1{\sigma }_1{\kappa }^4+a{\kappa }^2+\omega +2p^2(Z_3-c_1{\sigma }_1(5+3p^2)B^2)B^2\right]A{tanh}^p\tau \hfill \\ & & +\left[Z_2-2p^2(Z_4+Z_5)B^2\right]A^3{tanh}^{3p}\tau +c_1{\sigma }_6{A}^5{tanh}^{5p}\tau \hfill \\ & & +p(1+p)\left[Z_3-4c_1{\sigma }_1\left\{2+p(2+p)\right\}{B}^2\right]A{B}^2{tanh}^{p+2}\tau \hfill \\ & & +p(1+p)(2+p)(3+p)c_1{\sigma }_{1}A{B}^4{tanh}^{p+4}\tau \hfill \\ & & +p\left\{(1+p)Z_4+p{Z}_5\right\}{A}^3{B}^2{tanh}^{3p+2}\tau \hfill \\ & & +p\left\{(p-1)Z_4+p{Z}_5\right\}{A}^3{B}^2{tanh}^{3p-2}\tau \hfill \\ & & +p(p-1)\left\{Z_3-4c_1{\sigma }_1(2+p(p-2))B^2\right\}A{B}^2{tanh}^{p-2}\tau \hfill \\ & & +p(p-1)(p-2)(p-3)c_1{\sigma }_{1}A{B}^4{tanh}^{p-4}\tau =0.\hfill \end{array}$$

(25)

The appropriate balance leads to $p=1$, as for the case of the bright soliton. Substituting the resulting value of p into the Equation (25), and taking the set of linearly independent functions ${tanh}^p\tau$ for $p=1,3,5$ and putting all their coefficients zero, we can obtain the following three equalities:

$$\begin{array}{ccc}& & \omega +a{\kappa }^2+\left(2Z_3-Z_5{A}^2\right)B^2+c_1{\sigma }_1\left(3{\kappa }^4-16B^4\right)=0\hfill \end{array}$$

(26)

$$\begin{array}{ccc}& & \left\{Z_2-2\left(Z_4+Z_5\right)B^2\right\}{A}^2+2\left(Z_3-20c_1{\sigma }_1{B}^2\right)B^2=0\hfill \end{array}$$

(27)

$$\begin{array}{ccc}& & c_1{\sigma }_6{A}^4+24c_1{\sigma }_1{B}^4+\left(2Z_4+Z_5\right)A^2{B}^2=0\hfill \end{array}$$

(28)

From (27), the free parameter A becomes

$$A=\sqrt{-\frac{2\left(Z_3+10c_1{\sigma }_1{B}^2\right)B^2}{Z_2+\left(Z_4+Z_5\right)B^2}}$$

(29)

while from (26), the wave number is

$$\omega =-\frac{2Z_5\left(Z_3-20c_1{\sigma }_1{B}^2\right)B^4+\left\{Z_2-2\left(Z_4+Z_5\right)B^2\right\}\left\{a{\kappa }^2+2Z_3{B}^2+c_1{\sigma }_1\left(3{\kappa }^4-16B^4\right)\right\}}{Z_2-2\left(Z_4+Z_5\right)B^2}$$

(30)

in view of (29), and as long as,

$${\sigma }_1{c}_1^3\left[100{\sigma }_1{\sigma }_6+\left({\sigma }_2+{\sigma }_3+{\sigma }_4+{\sigma }_5\right)\left\{{\sigma }_2+{\sigma }_3-4\left({\sigma }_4+{\sigma }_5\right)\right\}\right]<0$$

(31)

and the constraint from $\omega$:

$$\left\{b-c_1\left({\sigma }_2-{\sigma }_3+{\sigma }_4+{\sigma }_5\right){\kappa }^2+c_2\left({\sigma }_8-{\sigma }_9\right)\kappa \right\}\ne 2c_1\left({\sigma }_2+{\sigma }_3+{\sigma }_4+{\sigma }_5\right)B^2$$

(32)

are satisfied. From (27) and (28), the free parameter B is

$$\begin{array}{ccc}\hfill B& =& \pm \frac{1}{2}\left[\frac{2c_1{\sigma }_1\left(2Z_4+7Z_5\right)Z_2-Z_3\left\{\left(Z_4+Z_5\right)\left(2Z_4+Z_5\right)-40c_1^2{\sigma }_1{\sigma }_6\right\}}{2c_1{\sigma }_1\left(100c_1^2{\sigma }_1{\sigma }_6-4Z_4^2-3Z_4{Z}_5+Z_5^2\right)}\right.\hfill \\ & \pm & {\left.\frac{\sqrt{{\left\{Z_3\left(Z_4+Z_5\right)-10c_1{\sigma }_1{Z}_2\right\}}^2\left\{{\left(2Z_4+Z_5\right)}^2-96c_1^2{\sigma }_1{\sigma }_6\right\}}}{2c_1{\sigma }_1\left(100c_1^2{\sigma }_1{\sigma }_6-4Z_4^2-3Z_4{Z}_5+Z_5^2\right)}\right]}^{\frac{1}{2}}\hfill \end{array}$$

(33)

where both constraints (21) and (22) are also required for dark soliton. Therefore, the single dark soliton solution for the concatenation model (1) is

$$\begin{array}{c}\hfill q(x,t)=Atanh\left[B(x-vt)\right]e^{i\left(\kappa x+\omega t+{\theta }_0\right)}\end{array}$$

(34)

where the parameters and associated restrictions are addressed previously.

Figure 2 shows the profile of a dark soliton, where the parameters are chosen as: $a=0.5$, $b=0.1$, $c_1=0.5$, $c_2=3.2$, ${\sigma }_1=1.2$, ${\sigma }_2=0.2$, ${\sigma }_3=2.3$, ${\sigma }_4=1.1$, ${\sigma }_5=1.7$, ${\sigma }_6=2.0$, ${\sigma }_7=-1.1$, ${\sigma }_8=1.3$, ${\sigma }_9=0.8$.

2.3. Singular Soliton (Type-I)

We accept the following wave form for the first type of singular soliton:

$$P(x,t)=A{\mathrm{csch}}^p\tau ,\tau =B(x-vt),$$

(35)

where A and B, are considered free parameters for the soliton. Substituting the chosen type-I form into (9) reduces the last into

$$\begin{array}{ccc}\hfill & & \left[p^2{Z}_3{B}^2+c_1{\sigma }_1{p}^4{B}^4-\left(\omega +a{\kappa }^2+3c_1{\sigma }_1{\kappa }^4\right)\right]A{csch}^p\tau \hfill \\ & & +\left[Z_2+p^2\left(Z_4+Z_5\right)B^2\right]A^3{csch}^{3p}\tau +c_1{\sigma }_6{A}^5{csch}^{5p}\tau \hfill \\ & & +p(1+p)\left[Z_3+2\left\{2+p\left(2+p\right)\right\}{c}_1{\sigma }_1{B}^2\right]A{B}^2{csch}^{p+2}\tau \hfill \\ & & +p(1+p)(2+p)(3+p)c_1{\sigma }_{1}A{B}^4{csch}^{p+4}\tau \hfill \\ & & +p\left[\left(1+p\right)Z_4+p{Z}_5\right]A^3{B}^2{csch}^{3p+2}\tau =0\hfill \end{array}$$

(36)

where the $Z^{\prime }s$ were defined above on (12)–(15). Balancing yields, as expected $p=1.$ Substituting this into the Equation (36), and taking the set of linearly independent functions ${csch}^p\left[\tau \right]$ for $p=1,3,5$ and setting all their coefficients as zero, we can obtain the following three equalities:

$$\begin{array}{ccc}& & \omega +a{\kappa }^2+3c_1{\sigma }_1{\kappa }^4-\left(Z_3+c_1{\sigma }_1{B}^2\right)B^2=0\hfill \end{array}$$

(37)

$$\begin{array}{ccc}& & \left[Z_2+\left(Z_4+Z_5\right)B^2\right]A^2+2\left(Z_3+10c_1{\sigma }_1{B}^2\right)B^2=0\hfill \end{array}$$

(38)

$$\begin{array}{ccc}& & c_1{\sigma }_6{A}^4+24c_1{\sigma }_1{B}^4+\left(2Z_4+Z_5\right)A^2{B}^2=0.\hfill \end{array}$$

(39)

Isolating $\omega$ on (37) leads to the same expression for the wave number as in the bright soliton (17), while solving for A (38) yields the same result as in (29), as long as the constraint (31) is satisfied. In view of the result for the parameter A, the free parameter B is the same as for the bright soliton given in (20), along with corresponding constraints. Therefore, the type-I singular soliton solution supported by the system (1) is

$$\begin{array}{c}\hfill q(x,t)=A\mathrm{csch}\left[B(x-vt)\right]e^{i\left(\kappa x+\omega t+{\theta }_0\right)}\end{array}$$

(40)

where the free parameters A and B are expressed in (29) and (20), respectively, as long as the solvability conditions (21), (22) and (31) are satisfied. The soliton speed was obtained earlier (6), while the wave number $\omega$ was the same as for bright soliton (17).

2.4. Singular Soliton (Type-II)

The second type of singular soliton solution will be the last type of soliton discussed in this work. To retrieve type-I singular solitons from our concatenation system (1), a wave form with the following structure is assumed:

$$\begin{array}{c}\hfill P(x,t)=A{\coth }^p\tau ,\tau =B(x-vt),\end{array}$$

(41)

where A and B are considered free parameters, p is a parameter whose value is to be determined, and v represents the soliton speed. The substitution of the hypothesis (41) into the real part of Equation (9) brings:

$$\begin{array}{ccc}& & -\left[3c_1{\sigma }_1{\kappa }^4+a{\kappa }^2+\omega +2p^2\left\{Z_3-c_1{\sigma }_1\left(5+3p^2\right)B^2\right\}{B}^2\right]A{coth}^p\tau \hfill \\ & & +\left[Z_2-2p^2(Z_4+Z_5)B^2\right]A^3{coth}^{3p}\tau +c_1{\sigma }_6{A}^5{coth}^{5p}\tau \hfill \\ & & +p(1+p)\left[Z_3-4c_1{\sigma }_1\left\{2+p\left(2+p\right)\right\}{B}^2\right]A{B}^2{coth}^{p+2}\tau \hfill \\ & & +p(1+p)(2+p)(3+p)c_1{\sigma }_{1}A{B}^4{coth}^{p+4}\tau \hfill \\ & & +p\left[(1+p)Z_4+p{Z}_5\right]A^3{B}^2{coth}^{3p+2}\tau \hfill \\ & & +p\left[(p-1)Z_4+p{Z}_5\right]A^3{B}^2{coth}^{3p-2}\tau \hfill \\ & & +p(p-1)\left[Z_3-4c_1{\sigma }_1\left\{2+p\left(p-2\right)\right\}{B}^2\right]A{B}^2{coth}^{p-2}\tau \hfill \\ & & +p(p-1)(p-2)(p-3)c_1{\sigma }_{1}A{B}^4{coth}^{p-4}\tau =0.\hfill \end{array}$$

(42)

Proper balancing leads to the same value of p as in the previous three cases. Substituting the resulting value of p into (43) and setting the coefficients of the resulting linearly independent functions as equal to zero yields the same three identities (26)–(28). Thus, the parameter findings for this subsection will be identical to those for the dark soliton case (29)–(33). Finally, the type-II soliton solution supported by the concatenation model is

$$\begin{array}{c}\hfill q(x,t)=Acoth\left[B(x-vt)\right]e^{i\left(\kappa x+\omega t+{\theta }_0\right)}\end{array}$$

(43)

where all parameter values are the same as the case of dark soliton, along with the corresponding solvability constraints.

3. Conservation Laws

It is important and imperative to investigate the conservation laws for any model to obtain a better understanding of its physical features. This gives a deeper perspective to the model and will lead to any future avenues that could be explored. For the conserved flow that renders a closed form of the respective nonlinear evolution equations (NLEEs), we let $q=u+iv$ and split the NLEE into a real system with conserved vectors; $(T^t,T^x)$ satisfies the $D_t{T}^t+D_x{T}^x=0$ along the solutions of the NLEEs. We have

(a)

conserved power density for ${\sigma }_2=0$ and ${\sigma }_4={\sigma }_5$, viz.,

$$\begin{array}{c}\hfill {\Phi }_{Pr}^t=\frac{1}{2}{\left|q\right|}^2\end{array}$$

(44)

and, additionally, for ${\sigma }_3=2{\sigma }_5$ and ${\sigma }_9=0$, we get

(b)

linear momentum density given by

$$\begin{array}{c}\hfill {\Phi }_M^t=a\mathcal{I}\left(q^{*}{q}_x\right)\end{array}$$

(45)

and

(c)

Hamiltonian density given by,

$$\begin{array}{ccc}\hfill {\Phi }_H^t& =& \frac{1}{2}a\mathcal{R}\left(q{q}_{xx}^{*}\right)+\frac{1}{4}b{\left|q\right|}^4\hfill \\ & +& c_1\left[\frac{1}{2}{\sigma }_1\mathcal{R}\left(q{q}_{xxxx}^{*}\right)+\frac{1}{2}\left({\sigma }_4+{\sigma }_5\right)\left\{\mathcal{R}\left(q{q}_{xx}^{*}\right)+{\left|q\right|}^2{\left|q_x\right|}^2\right\}+\frac{1}{6}{\sigma }_6{\left|q\right|}^6\right]\hfill \\ & +& c_2\left[\frac{1}{2}{\sigma }_7\mathcal{I}\left(q^{*}{q}_{xxx}\right)+\frac{1}{4}{\sigma }_8\mathcal{I}\left(q^{*}{q}_x\right)\right].\hfill \end{array}$$

(46)

The conserved quantity are therefore, respectively, given as:

$$\begin{array}{c}\hfill Pr={\int }_{-\infty }^{\infty }{\Phi }_{Pr}^{t}dx={\int }_{-\infty }^{\infty }{\left|q\right|}^{2}dx=\frac{2A^2}{B},\end{array}$$

(47)

$$\begin{array}{c}\hfill M=a{\int }_{-\infty }^{\infty }{\Phi }_M^{t}dx=a{\int }_{-\infty }^{\infty }\mathcal{I}\left(q^{*}{q}_x\right)dx=\frac{a}{2i}{\int }_{-\infty }^{\infty }\left(q^{*}{q}_x-q{q}_x^{*}\right)dx=\frac{2a\kappa A^2}{B},\end{array}$$

(48)

and

$$\begin{array}{ccc}\hfill H& =& {\int }_{-\infty }^{\infty }{\Phi }_H^{t}dx\hfill \\ & =& \frac{a{A}^2}{3B}\left(B^2+3{\kappa }^2\right)+\frac{b{A}^4}{3B}\hfill \\ & +& c_1\left[\frac{{\sigma }_1{A}^2}{15B}\left(7B^4+30{\kappa }^2{A}^2{B}^2+15{\kappa }^4\right)+\frac{\left({\sigma }_4+{\sigma }_5\right)A^2}{15B}\left\{5\left(B^2+3{\kappa }^2\right)+2A^2\left(B^2+5{\kappa }^2\right)\right\}\right]\hfill \\ & +& \frac{8c_1{\sigma }_6{A}^6}{45B}+c_2\left[\frac{{\sigma }_7\kappa A^2}{B}\left(B^2+{\kappa }^2\right)-\frac{{\sigma }_8\kappa A^2}{2B}\right]\hfill \end{array}$$

(49)

4. Conclusions

Using the method of indeterminate coefficients, the current paper recovers the 1-soliton solution to the recently constructed concatenation model. Later, the bright 1-soliton solution is utilized to establish the conserved quantities derived for the model using the multipliers method. Power, linear momentum, and the Hamiltonian are the three conserved quantities that were derived. These results open up a floodgate of opportunities for further studies of this model. The perturbation terms will be incorporated and, thus, the extended version of the concatenation model will be addressed by various means, including the semi-inverse variational principle. The variational principle will be addressed to recover the soliton parameter dynamics. The established conservation laws will next be implemented to obtain the adiabatic parameter dynamics of the solitons. The stochasticity effect will also be addressed by studying the corresponding Langevin equations to give way to the mean free velocity of the soliton. The collective variables approach would also lead to the parameter dynamics of the solitons. From a numerical standpoint, the model will later be examined numerically with the variational iteration method, Adomian decomposition scheme and Laplace–Adomian decomposition scheme. Evidently, a lot of work lies ahead and the results should subsequently be revealed.