Optical Solitons with Beta and M-Truncated Derivatives in Nonlinear Negative-Index Materials with Bohm Potential

Abstract

In this article, we explore solitary wave structures in nonlinear negative-index materials with beta and M-truncated fractional derivatives with the existence of a Bohm potential. The consideration of Bohm potential produced quantum phase behavior in electromagnetic waves. The applied technique is the New extended algebraic method. By use of this approach, acquired solutions convey various types of new families containing dark, dark-singular, dark-bright, and singular solutions of Type 1 and 2. Moreover, the constraint conditions for the presence of the obtained solutions are a side-effect of this technique. Finally, graphical structures are depicted.

1. Introduction

The subject of fractional calculus (calculus of integrals and derivatives of any arbitrary real or complex order) has attained great importance over the previous three decades or so, due to it having major applications in science and engineering. In reality, it offers many possibly valuable techniques for resolving differential and integral equations and numerous other concerns that include particular equations, mathematical physics features, as their augmentations and speculations in a single variable, and n is the limit from there. A portion of the regions of current utilization of fractional calculus incorporates rheology, electrical networks, probability and statistics, electrochemistry of corrosion, chemical physics, control theory of dynamical systems, optics and signal processing, etc. As of late, numerous attempts have been dedicated to this subject; a couple of them are accounted for in [1,2,3,4]. The examination of evaluating fractional derivative operators is consistently a hotly debated issue of research. Many attempts have been given lately, and numerous revelations have been made in this course; some of them are recorded in [5,6,7,8,9].

The model which represents the motion of electromagnetic waves is termed the perturbed nonlinear Schrodinger equation (NLSE). The new type of soliton solutions of time-fractional perturbed NLSE with conformable derivative in nonlinear negative-index materials with Bohm potential are discussed in [10]. The Perturbed NLSE with a conformable space–time fractional model is studied in [11]. Now in the present work, we employ the New extended algebraic method to find solitary wave solutions of fractional Perturbed NLSE with beta and M-truncated derivatives. The results are new and not seen in the literature.

The equation of consideration here depicts the elements of soliton propagation through optical meta-materials with self-steepening (SS), interemodal dispersion (IMD), Bohm potential and nonlinear dispersion (ND) and is of the type [12,13,14]

$$\begin{array}{cc}\hfill i{D}_{\tau }^{\mu }v+a{v}_{\xi \xi }+{(b|v|}^2)v=& i({\theta }_1\frac{v}{\left|v\right|}{\left|v\right|}_{\xi \xi }+{\theta }_2{v}_{\xi \xi }+{\theta }_3{(|v|}^2{v)}_{\xi }+{\theta }_4{(|v|}^2{)}_{\xi }v)+c_1{(|v|}^2{v)}_{\xi \xi }\hfill \\ & c_2{\left|v\right|}^2{(v)}_{\xi \xi }+c_3{v}^2{(v^{\star })}_{\xi \xi }.\hfill \end{array}$$

(1)

The complex function $v(\xi ,\tau )$ is the dependent variable which shows the wave profile, $\mu$ is the fractional parameter with a value between zero and one, $\xi$ indicates the non-dimensional distance across the fiber and the temporal component is $\tau$. Furthermore, coefficient parameter a is the group velocity dispersion while b describes the cubic nonlinearity known as self-phase modulation.

The paper divided as follows. In Section 1.1 and Section 1.2 main definitions for fractional derivatives are reported. In Section 2, we present governing equations and mathematical analysis of the given equation. In Section 3, we present the fundamentals of the new extended direct algebraic method and use this algorithm to find the DE. In Section 4, we describe necessary and sufficient conditions and physical features of some obtained solutions. The conclusion is stated at the end.

1.1. Beta Derivative

Definition 1.

The beta derivative is defined by [8]:

$${}_0^E{D}_{\xi }^{\mu }(\mathbb{G}(\xi ))=\underset{\epsilon \to 0}{lim}\frac{\mathbb{G}\left(\xi +\epsilon \left(\xi +\frac{1}{\Gamma (\alpha )}\right)\right)-\mathbb{G}(\xi )}{\epsilon },$$

(2)

with a few properties as labeled below:

Theorem 1.

Let $0<\mu \le 1$, $\delta ,\lambda \in \mathbb{R}$. f and g are ${\mu }^{th}$ order differentiable functions at $\tau >0$. Then, we have:

1:${}_0^E{D}_{\xi }^{\mu }\left(\delta f(\xi )+\lambda g(\xi )\right)=\delta {}_0^E{D}_{\xi }^{\mu }\left(f(\xi )\right)+\lambda {}_0^E{D}_{\xi }^{\mu }\left(g(\xi )\right)$,

2:$D_{\xi }^{\mu }(c)=0$, here c is constant.

3:${}_0^E{D}_{\xi }^{\mu }(f(\xi )\ast g(\xi ))=g(\xi ){}_0^E{D}_{\xi }^{\mu }\left(f(\xi )\right)+f(\xi ){}_0^E{D}_{\xi }^{\mu }\left(g(\xi )\right)$,

4:${}_0^E{D}_{\xi }^{\mu }\left(\frac{f(\xi )}{g(\xi )}\right)=\frac{g(\xi ){}_0^E{D}_{\xi }^{\mu }f(\xi )-f(\xi ){}_0^E{D}_{\xi }^{\mu }g(\xi )}{g^2(\xi )}$,

5:For $\epsilon ={\left(\xi +\frac{1}{\Gamma (\mu )}\right)}^{1-\mu }h,h\to 0$ when $\epsilon \to 0$; therefore, we have

$${}_0^E{D}_{\xi }^{\mu }f(\xi )={\left(\xi +\frac{1}{\Gamma (\mu )}\right)}^{1-\mu }\frac{df(\xi )}{d\xi }$$

with $\xi =\frac{\nu }{\mu }{\left(\xi +\frac{1}{\Gamma (\mu )}\right)}^{\mu }$, where ν is a constant.

6:${}_0^E{D}_{\xi }^{\mu }\left(\frac{f(\tau )}{g(\xi )}\right)=\tau \frac{df(\tau )}{d\tau }$.

1.2. Truncated M-Fractional Derivative

Definition 2.

The truncated Mittag–Leffler function with one parameter is defined as:

$${}_i{T}_{\Theta }(y)=\sum _{k=0}^i\frac{y^k}{\Gamma (\Theta k+1)},$$

(3)

where $\Theta >0$ & $y\in \mathbb{C}.$ It is characterized by a non-fuzzy idea as illustrated below.

Definition 3.

Assume $G:[0,\infty )\to \mathbb{R}$, and $\mu \in (0,1)$ the TMD of G of order μ is given by:

$${}_i{\mathbb{D}}_M^{\mu ,\Theta }G(\tau )=\underset{\epsilon \to 0}{lim}\frac{G(\tau +{}_i{T}_{\Theta }(\epsilon {\tau }^{-\mu }))-G(\tau )}{\epsilon },$$

(4)

for $\tau >0$ and ${}_i{T}_{\Theta }(.)$, $\Theta >0$ is defined in the above definition.

Theorem 2.

Let $\mu \in (0,1]$, $\Theta >0$. G and H are ${\mu }^{th}$ order differentiable functions at $\tau >0$. Then, we have:

1-${}_i{\mathbb{D}}_M^{\mu ,\Theta }(p_{1}G+p_{2}H)=p_1{}_i{\mathbb{D}}_M^{\mu ,\Theta }(G)+p_2{}_i{\mathbb{D}}_M^{\mu ,\Theta }(H),$ where $p_1,p_2\in \mathbb{R}$

2-${}_i{\mathbb{D}}_M^{\mu ,\Theta }({\tau }^{\sigma })=\sigma {\tau }^{\sigma -\mu },\sigma \in \mathbb{R}$,

3-${}_i{\mathbb{D}}_M^{\mu ,\Theta }(GH)=G{}_i{\mathbb{D}}_M^{\mu ,\Theta }(H)+H{}_i{\mathbb{D}}_M^{\mu ,\Theta }(G)$,

4-${}_i{\mathbb{D}}_M^{\mu ,\Theta }(\frac{G}{H})=\frac{G{}_i{\mathbb{D}}_M^{\mu ,\Theta }(H)-H{}_i{\mathbb{D}}_M^{\mu ,\Theta }(G)}{H^2},$

5-${}_i{\mathbb{D}}_M^{\mu ,\Theta }(G)(\tau )=\frac{{\tau }^{1-\mu }}{\Gamma (\Theta +1)}\frac{dG}{d\tau }$,

6-${}_i{\mathbb{D}}_M^{\mu ,\Theta }(G\circ H)(\tau )=f^{\prime }(H(\tau )){}_i{\mathbb{D}}_M^{\mu ,\Theta }H(\tau )$.

2. Governing Equations

By considering a beta derivative, the above equation can be composed as:

$$\begin{gathered} i{}_0^E{D}_{\tau }^{\mu }(v)+a{}_0^E{D}_{\xi }^{2\mu }(v)+{(b|u|}^{2}v)=i\left({\theta }_1\frac{v}{\left|v\right|}{}_0^E{D}_{\xi }^{2\mu }(|v|)+{\theta }_2{}_0^E{D}_{\xi }^{\mu }(v)+{\theta }_3{}_0^E{D}_{\xi }^{\mu }{(|v|}^{2}v)+{\theta }_4{}_0^E{D}_{\xi }^{\mu }{(|v|}^2)v\right) \\ +c_1{}_0^E{D}_{\xi }^{2\mu }{(|v|}^{2}v)+c_2{\left|v\right|}^2{}_0^E{D}_{\xi }^{2\mu }(v)+c_3{v}^2{}_0^E{D}_{\xi }^{2\mu }(v^{\star }), \end{gathered}$$

(5)

Meanwhile by taking the M-truncated derivative into account, the model under investigation takes the structure as below:

$$\begin{array}{cc}\hfill & i{}_0^E{D}_{M,\tau }^{\mu ,\Theta }(v)+a{}_0^E{D}_{M,\xi }^{2\mu ,\Theta }(v)+{(b|u|}^{2}v)=i({\theta }_1\frac{v}{\left|v\right|}{}_0^E{D}_{M,\xi }^{2\mu ,\Theta }(|v|)+{\theta }_2{}_0^E{D}_{M,\xi }^{\mu ,\Theta }(v)+{\theta }_3{}_0^E{D}_{M,\xi }^{\mu ,\Theta }{(|v|}^{2}v)\hfill \\ & +{\theta }_4{}_0^E{D}_{M,\xi }^{\mu ,\Theta }{(|v|}^2)v)+c_1{}_0^E{D}_{M,\xi }^{2\mu ,\Theta }{(|v|}^{2}v)+c_2{\left|v\right|}^2{}_0^E{D}_{M,\xi }^{2\mu ,\Theta }(v)+c_3{v}^2{}_0^E{D}_{M,\xi }^{2\mu ,\Theta }(v^{\star }),\hfill \end{array}$$

(6)

in above ${}_0^E{D}_{M,\tau }^{\mu ,\Theta }$ and ${}_0^E{D}_{M,\xi }^{\mu ,\Theta }$ are M-truncated derivatives with $\tau$ and $\xi$, respectively.

Mathematical Survey

For the solution of Equation (1), Equations (5) and (6) the first step is follows:

$$v(\xi ,\tau )=u(\eta )e^{i\psi (\xi ,\tau )},$$

(7)

where the $v(\xi ,\tau )$ represents pulse shape of soliton. In sense of the beta derivative, we have

$$\eta =\frac{1}{\mu }{\left(\xi +\frac{1}{\Gamma (\mu )}\right)}^{\mu }-\frac{\nu }{\mu }{\left(\tau +\frac{1}{\Gamma (\mu )}\right)}^{\mu }$$

(8)

and

$$\psi (\xi ,\tau )=-\frac{k}{\mu }{\left(\xi +\frac{1}{\Gamma (\mu )}\right)}^{\mu }+\frac{\omega }{\mu }{\left(\tau +\frac{1}{\Gamma (\mu )}\right)}^{\mu }+{\theta }_0(\eta ).$$

(9)

By the virtue of M-truncated derivative, we have:

$$\eta =\frac{\Gamma (\Theta +1)}{\mu }({\xi }^{\mu }-\nu {\tau }^{\mu }),$$

(10)

$$\psi (\xi ,\tau )=-\frac{\Gamma (\Theta +1)}{\mu }(\kappa {\xi }^{\mu }-\nu \omega {\tau }^{\mu })+{\theta }_0(\eta ),$$

(11)

where $\psi (\xi ,\tau )$, $\kappa$, $\omega$, $\nu$ and ${\theta }_0(\eta )$ are the phase component, frequency, wave number, speed and phase function of soliton, respectively. Putting these values in Equation (1), Equations (5) and (6), then the imaginary part results in:

$$\nu ={\theta }_1-2a\kappa$$

(12)

and

$$3{\theta }_2+2{\theta }_4-2\kappa (3c_1+c_2-c_3)=0,$$

(13)

while the real part implies

$$a{u}^{″}-(\omega +a{\kappa }^2+\kappa {\theta }_2)u+(b-\kappa {\theta }_3+{\kappa }^2{c}_1+{\kappa }^2{c}_2+{\kappa }^2{c}_3)u^3-(3c_1+c_2+c_3)u^2{u}^{″}-6c_1{u}^{{\prime }^2}=0.$$

(14)

We apply the accompanying changes $c_1=0$, and $c_2=-c_3$ in Equation (14) for its solution as

$$a{u}^{″}-(\omega +a{\kappa }^2+{\theta }_2\kappa )u+(b-\kappa {\theta }_3)u^3=0,$$

(15)

where

$$3{\theta }_2+2{\theta }_4+4\kappa c_3=0.$$

(16)

3. Applications

This section is devoted to the application of the method with two different definitions of the derivatives.

3.1. The New Extended Direct Algebraic Method

In this section, the general procedure of this method [15,16,17] is reported in two steps as given below:

Step: 1 

Suppose that the given nonlinear PDE is of the form:

$$G(v,v_{\xi },v_{\tau },v_{\xi \xi },v_{\tau \tau },...)=0,$$

(17)

where v represents the dependent variable and $\xi$, $\tau$ are the independent variables.

By using the wave transformation:

$$v(\xi ,\tau )=u(\eta ),\eta =\xi -c\tau .$$

Equation (17) can be transformed into nonlinear ODE:

$$J(u,u^{\prime },u^{″},u^{‴}...)=0.$$

(18)

Step: 2 

We assume the solution of the ODE (18) of the type:

$$u(\eta )=\sum _{i=0}^m{b}_i{Z}^i(\eta ),$$

(19)

where $b_i$$(0<i\le n)$ are the coefficients and $Z(\eta )$ satisfies the ODE of the type:

$$Z^{\prime }(\eta )=ln(B)\left(\alpha +\beta Z(\eta )+\gamma Z^2(\eta )\right),B\ne 0,1,$$

(20)

where $\alpha$, $\beta$ and $\gamma$ are the constants. The Equation (20) has the solution written as:

$\mathbf{1}$: When ${\beta }^2-4\alpha \gamma <0$ and $\gamma \ne 0$,

$Z_1(\eta )=-\frac{\beta }{2\gamma }+\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2\gamma }{tan}_B\left(\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2}\eta \right)$,

$Z_2(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2\gamma }{cot}_B\left(\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2}\eta \right)$,

$Z_3(\eta )=-\frac{\beta }{2\gamma }+\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2\gamma }\left({tan}_B\left(\sqrt{-({\beta }^2-4\alpha \gamma )}\eta \right)\pm \sqrt{mn}{sec}_B\left(\sqrt{-({\beta }^2-4\alpha \gamma )}\eta \right)\right)$,

$Z_4(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{2\gamma }\left({cot}_B\left(\sqrt{-({\beta }^2-4\alpha \gamma )}\eta \right)\pm \sqrt{mn}{csc}_B\left(\sqrt{-({\beta }^2-4\alpha \gamma )}\eta \right)\right)$,

$Z_5(\eta )=-\frac{\beta }{2\gamma }+\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{4\gamma }\left({tan}_B\left(\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{4}\eta \right)-{cot}_B\left(\frac{\sqrt{-({\beta }^2-4\alpha \gamma )}}{4}\eta \right)\right)$.

$\mathbf{2}$: When ${\beta }^2-4\alpha \gamma >0$ and $\gamma \ne 0$,

$Z_6(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2\gamma }{tan}_B\left(\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2}\eta \right)$,

$Z_7(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2\gamma }{cot}_B\left(\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2}\eta \right)$,

$Z_8(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2\gamma }\left({tanh}_B\left(\sqrt{{\beta }^2-4\alpha \gamma }\eta \right)\pm \iota \sqrt{mn}\mathrm{sech}{}_B\left(\sqrt{{\beta }^2-4\alpha \gamma }\eta \right)\right)$,

$Z_9(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{2\gamma }\left({coth}_B\left(\sqrt{{\beta }^2-4\alpha \gamma }\eta \right)\pm \sqrt{mn}\mathrm{csch}{}_B\left(\sqrt{{\beta }^2-4\alpha \gamma }\eta \right)\right)$,

$Z_{10}(\eta )=-\frac{\beta }{2\gamma }-\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{4\gamma }\left({tanh}_B\left(\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{4}\eta \right)+{coth}_B\left(\frac{\sqrt{{\beta }^2-4\alpha \gamma }}{4}\eta \right)\right)$.

$\mathbf{3}$: When $\alpha \gamma >0$ and $\beta =0$,

$Z_{11}(\eta )=\sqrt{\frac{\alpha }{\gamma }}{tan}_B(\sqrt{\alpha \gamma }\eta )$,

$Z_{12}(\eta )=-\sqrt{\frac{\alpha }{\gamma }}{cot}_B(\sqrt{\alpha \gamma }\eta )$,

$Z_{13}(\eta )=\sqrt{\frac{\alpha }{\gamma }}({tan}_B(2\sqrt{\alpha \gamma }\eta )\pm \sqrt{mn}{sec}_B(2\sqrt{\alpha \gamma }\eta ))$,

$Z_{14}(\eta )=-\left(\sqrt{\frac{\alpha }{\gamma }}({cot}_B(2\sqrt{\alpha \gamma }\eta )\pm \sqrt{mn}{csc}_B(2\sqrt{\alpha \gamma }\eta )\right)$,

$Z_{15}(\eta )=\frac{1}{2}\sqrt{\frac{\alpha }{\gamma }}\left({tan}_B\left(\frac{\sqrt{\alpha \gamma }}{2}\eta \right)-{cot}_B\left(\frac{\sqrt{\alpha \gamma }}{2}\eta \right)\right)$.

$\mathbf{4}$: When $\alpha \gamma <0$ and $\beta =0$,

$Z_{16}(\eta )=-\sqrt{-\frac{\alpha }{\gamma }}{tanh}_B(\sqrt{-\alpha \gamma }\eta )$,

$Z_{17}(\eta )=-\sqrt{-\frac{\alpha }{\gamma }}{coth}_B(-\sqrt{\alpha \gamma }\eta )$,

$Z_{18}(\eta )=-\sqrt{\frac{-\alpha }{\gamma }}({tanh}_B(2\sqrt{-\alpha \gamma }\eta )\pm \iota \sqrt{mn}\mathrm{sech}{}_B(2\sqrt{-\alpha \gamma }\eta ))$,

$Z_{19}(\eta )=-\sqrt{-\frac{\alpha }{\gamma }}\left({coth}_B(2\sqrt{-\alpha \gamma }\eta )\pm \sqrt{mn}\mathrm{csch}{}_B(2\sqrt{-\alpha \gamma }\eta )\right)$,

$Z_{20}(\eta )=-\frac{1}{2}\sqrt{\frac{-\alpha }{\gamma }}\left({tanh}_B\left(\frac{\sqrt{-\alpha \gamma }}{2}\eta \right)+{coth}_B\left(\frac{\sqrt{-\alpha \gamma }}{2}\eta \right)\right)$.

$\mathbf{5}$: When $\beta =0$ and $\gamma =\alpha$,

$Z_{21}(\eta )={tan}_B(\alpha \eta )$,

$Z_{22}(\eta )=-{cot}_B(\alpha \eta )$,

$Z_{23}(\eta )={tan}_B(2\alpha \eta )\pm \sqrt{mn}{sec}_B(2\alpha \eta )$,

$Z_{24}(\eta )=-{cot}_B(2\alpha \eta )\pm \sqrt{mn}{csc}_B(2\alpha \eta )$,

$Z_{25}(\eta )=\frac{1}{2}\left({tan}_B\left(\frac{\alpha }{2}\eta \right)-{cot}_B\left(\frac{\alpha }{2}\eta \right)\right)$.

$\mathbf{6}$: When $\beta =0$ and $\gamma =-\alpha$,

$Z_{26}(\eta )=-{tanh}_B(\alpha \eta )$,

$Z_{27}(\eta )=-{coth}_B(\alpha \eta )$,

$Z_{28}(\eta )=-{tanh}_B(2\alpha \eta )\pm \iota \sqrt{mn}\mathrm{sech}{}_B(2\alpha \eta )$,

$Z_{29}(\eta )=-{coth}_B(2\alpha \eta )\pm \sqrt{mn}\mathrm{csch}{}_B(2\alpha \eta )$,

$Z_{30}(\eta )=-\frac{1}{2}\left({tanh}_B\left(\frac{\alpha }{2}\eta \right)+{coth}_B\left(\frac{\alpha }{2}\eta \right)\right)$.

$\mathbf{7}$: When ${\beta }^2=4\alpha \gamma$,

$Z_{31}(\eta )=\frac{-2\alpha (\beta \eta ln(B)+2)}{{\beta }^2\eta ln(B)}$.

$\mathbf{8}$: When $\beta =\rho$, $\alpha =q\rho (q\ne 0)$ and $\gamma =0$,

$Z_{32}(\eta )=B^{\rho \eta }-q$.

$\mathbf{9}$: When $\beta =\gamma =0$,

$Z_{33}(\eta )=\alpha \eta ln(B)$.

$\mathbf{10}$: When $\beta =\alpha =0$,

$Z_{34}(\eta )=\frac{-1}{\gamma \eta ln(B)}$.

$\mathbf{11}$: When $\alpha =0$ and $\beta \ne 0$,

$Z_{35}(\eta )=-\frac{s\beta }{\gamma ({cosh}_B(\beta \eta )-{sinh}_B(\beta \eta )+s)}$,

$Z_{36}(\eta )=-\frac{\beta ({sinh}_B(\beta \eta )+{cosh}_B(\beta \eta ))}{\gamma ({sinh}_B(\beta \eta )+{cosh}_B(\beta \eta )+r)}$.

$\mathbf{12}$: When $\beta =\rho$, $\gamma =m\rho (m\ne 0)$ and $\alpha =0$,

$Z_{37}(\eta )=\frac{s{B}^{\rho \eta }}{r-ms{B}^{\rho \eta }}$.

Now, the hyperbolic and trigonometric functions are given as follows:

${sinh}_B(\eta )=\frac{q{B}^{\eta }-s{B}^{-\eta }}{2}$, ${cosh}_B(\eta )=\frac{r{B}^{\eta }+s{B}^{-\eta }}{2}$,

${tanh}_B(\eta )=\frac{r{B}^{\eta }-s{B}^{-\eta }}{r{B}^{\eta }+s{B}^{-\eta }}$, ${coth}_B(\eta )=\frac{r{B}^{\eta }+s{B}^{-\eta }}{r{B}^{\eta }-s{B}^{-\eta }}$,

$\mathrm{csch}{}_B(\eta )=\frac{2}{r{B}^{\eta }-s{B}^{-\eta }}$, $\mathrm{sech}{}_B(\eta )=\frac{2}{r{B}^{\eta }+s{B}^{-\eta }}$,

${sin}_B(\eta )=\frac{r{B}^{\iota \eta }-s{B}^{-\iota \eta }}{2\iota }$, ${cos}_B(\eta )=\frac{r{B}^{\iota \eta }+s{B}^{-\iota \eta }}{2}$,

${tan}_B(\eta )=-i\frac{r{B}^{\iota \eta }-s{B}^{-\iota \eta }}{r{B}^{\iota \eta }+s{B}^{-\iota \eta }}$, ${cot}_B(\eta )=i\frac{r{B}^{\iota \eta }+s{B}^{-\iota \eta }}{r{B}^{\iota \eta }-s{B}^{-\iota \eta }}$,

${csc}_B(\eta )=\frac{2\iota }{r{B}^{\iota \eta }-s{B}^{-\iota \eta }}$, $\mathrm{sech}{}_B(\eta )=\frac{2}{r{B}^{\iota \eta }+s{B}^{-\iota \eta }}$,

where r and s are constants which are known as the deformation parameters.

3.2. Application to the NLS Equation

Let us take the transformation of the form:

$$v(\xi ,\tau )=u(\eta )e^{i\psi (\xi ,\tau )},$$

(21)

by using the balancing scheme on these terms $u^3$ and $u^{″}$ of Equation (15), assigns $m=1$ to Equation (19). We obtain the following transformation:

$$u(\eta )=b_0+b_{1}Z(\eta ),$$

(22)

putting Equation (22) into (15), then collecting the coefficients of different powers of $Z(\eta )$, we obtain a system of the following algebraic equations:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill Z^3(\eta )& :2a{b}_1{\gamma }^{2}ln{(B)}^2-b_1^3\kappa {\theta }_3+b_1^{3}b=0,\hfill \\ \hfill Z^2(\eta )& :3a{b}_1\beta ln{(B)}^2\gamma -3b_0{b}_1^2\kappa {\theta }_3+3b_0{b}_1^{2}b=0,\hfill \\ \hfill Z^1(\eta )& :a{b}_1{\beta }^{2}ln{(B)}^2+2a{b}_1\gamma ln{(B)}^2\alpha -b_{1}a{\kappa }^2-b_1{\theta }_2\kappa -b_1\omega -3b_0^2{b}_1\kappa {\theta }_3+3b_0^2{b}_{1}b=0,\hfill \\ \hfill Z^0(\eta )& :a{b}_1\beta ln{(B)}^2\alpha -b_{0}a{\kappa }^2-b_0{\theta }_2\kappa -b_0\omega -b_0^3\kappa {\theta }_3+b_0^{3}b=0.\hfill \end{array}\end{array}$$

(23)

The solution of the system of Equation (23) by use of Maple for $b_0$, $b_1$ and $\omega$, we obtain the following values:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill b_0& =\pm \frac{\gamma a\beta ln(B)}{\sqrt{-\frac{2{\gamma }^{2}a}{b-\kappa {\theta }_3}}(b-\kappa {\theta }_3)},b_1=\pm \sqrt{-\frac{2{\gamma }^{2}a}{b-\kappa {\theta }_3}}ln(B),\hfill \\ \hfill \omega & =2\gamma ln{(B)}^{2}a\alpha -\frac{1}{2}ln{(B)}^{2}a{\beta }^2-a{\kappa }^2-\kappa {\theta }_2.\hfill \end{array}\end{array}$$

(24)

Let us consider $\Delta ={\beta }^2-4\alpha \gamma$ and $\Pi =\gamma \sqrt{-\frac{2a}{b-\kappa {\theta }_3}}$.

Then $b_0=\pm \frac{a\beta ln(B)}{\Pi (b-\kappa {\theta }_3)}$ and $b_1=\pm \Pi ln(B)$.

From Equations (21), (22) and (24) and the different cases of solutions of Equation (20), we obtain the solutions which come after:

$\mathbf{Case}\mathbf{1}$. When $\Delta <0$ and $\gamma \ne 0$,

$v_1(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(\sqrt{-\Delta }{tan}_B(\frac{\sqrt{-\Delta }}{2}\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_2(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(\sqrt{-\Delta }{cot}_B(\frac{\sqrt{-\Delta }}{2}\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_3(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{-\Delta }\left({tan}_B(\sqrt{-\Delta }\eta )\pm \sqrt{mn}{sec}_B(\sqrt{-\Delta }\eta )\right)\right)e^{i\psi (\xi ,\tau )}$,

$v_4(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{-\Delta }\left(-{cot}_B(\sqrt{-\Delta }\eta )\pm \sqrt{mn}{csc}_B(\sqrt{-\Delta }\eta )\right)\right)e^{i\psi (\xi ,\tau )}$,

$v_5(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{-\Delta }({tan}_B\left(\frac{\sqrt{-\Delta }}{4}\eta \right)\pm {cot}_B\left(\frac{\sqrt{-\Delta }}{4}\eta \right)\right))e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{2}$. When $\Delta >0$ and $\gamma \ne 0$,

$v_6(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(\sqrt{\Delta }{tanh}_B(\frac{\sqrt{\Delta }}{2}\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_7(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(\sqrt{\Delta }{coth}_B(\frac{\sqrt{\Delta }}{2}\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_8(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{\Delta }\left(-{tanh}_B(\sqrt{\Delta }\eta )\pm i\sqrt{mn}\mathrm{sech}{}_B(\sqrt{\Delta }\eta )\right)\right)e^{i\psi (\xi ,\tau )}$,

$v_9(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{\Delta }\left(-{coth}_B(\sqrt{\Delta }\eta )\pm \sqrt{mn}\mathrm{csch}{}_B(\sqrt{\Delta }\eta )\right)\right)e^{i\psi (\xi ,\tau )}$,

$v_{10}(\xi ,\tau )=\pm \frac{aln(B)}{\Pi (b-\kappa {\theta }_3)}\left(2\beta -\sqrt{\Delta }(-{tanh}_B\left(\frac{\sqrt{\Delta }}{4}\eta \right)\pm {coth}_B\left(\frac{\sqrt{\Delta }}{4}\eta \right)\right))e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{3}$. When $\gamma \alpha >0$ and $\beta =0$,

$v_{11}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{\alpha \gamma }{tan}_B(\sqrt{\alpha \gamma }\eta )e^{i\psi (\xi ,\tau )}$,

$v_{12}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{\alpha \gamma }{cot}_B(\sqrt{\alpha \gamma }\eta )e^{i\psi (\xi ,\tau )}$,

$v_{13}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{\alpha \gamma }\left({tan}_B(2\sqrt{\alpha \gamma }\eta )\pm \sqrt{mn}{sec}_B(2\sqrt{\alpha \gamma }\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{14}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{\alpha \gamma }\left(-{cot}_B(2\sqrt{\alpha \gamma }\eta )\pm \sqrt{mn}{csc}_B(2\sqrt{\alpha \gamma }\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{15}(\xi ,\tau )=\pm \frac{1}{2}\Pi ln(B)\sqrt{\alpha \gamma }\left({tan}_B\left(\frac{\sqrt{\alpha \gamma }}{2}\eta \right)\pm {cot}_B\left(\frac{\sqrt{\alpha \gamma }}{2}\eta \right)\right)e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{4}$. When $\gamma \alpha <0$ and $\beta =0$,

$v_{16}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{-\alpha \gamma }{tanh}_B(\sqrt{-\alpha \gamma }\eta )e^{i\psi (\xi ,\tau )}$,

$v_{17}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{-\alpha \gamma }{coth}_B(\sqrt{-\alpha \gamma }\eta )e^{i\psi (\xi ,\tau )}$,

$v_{18}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{-\alpha \gamma }\left(-{tanh}_B(2\sqrt{-\alpha \gamma }\eta )\pm i\sqrt{mn}\mathrm{sech}{}_B(2\sqrt{-\alpha \gamma }\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{19}(\xi ,\tau )=\pm \Pi ln(B)\sqrt{-\alpha \gamma }\left(-{coth}_B(2\sqrt{-\alpha \gamma }\eta )\pm \sqrt{mn}\mathrm{csch}{}_B(2\sqrt{-\alpha \gamma }\eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{20}(\xi ,\tau )=\pm \frac{1}{2}\Pi ln(B)\sqrt{-\alpha \gamma }\left(-{tanh}_B\left(\frac{\sqrt{-\alpha \gamma }}{2}\eta \right)\pm {coth}_B\left(\frac{\sqrt{-\alpha \gamma }}{2}\eta \right)\right)e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{5}$. When $\beta =0$ and $\gamma =\alpha$,

$v_{21}(\xi ,\tau )=\pm \Pi ln(B)\alpha {tan}_B(\alpha \eta )e^{i\psi (\xi ,\tau )}$,

$v_{22}(\xi \tau )=\pm \Pi ln(B)\alpha {cot}_B(\alpha \eta )e^{i\psi (\xi ,\tau )}$,

$v_{23}(\xi ,\tau )=\pm \Pi ln(B)\alpha ({tan}_B(2\alpha \eta )\pm \sqrt{mn}{sec}_B(2\alpha \eta ))e^{i\psi (\xi ,\tau )}$,

$v_{24}(\xi ,\tau )=\pm \Pi ln(B)\alpha \left(-{cot}_B(2\alpha \eta )\pm \sqrt{mn}{csc}_B(2\alpha \eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{25}(\xi ,\tau )=\pm \frac{1}{2}\Pi ln(B)\alpha \left({tan}_B(\frac{\alpha }{2}\eta )\pm {cot}_B(\frac{\alpha }{2}\eta )\right)e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{6}$. When $\beta =0$ and $\gamma =-\alpha$,

$v_{26}(\xi ,\tau )=\pm \Pi ln(B)\alpha {tanh}_B(\alpha \eta )e^{i\psi (\xi ,\tau )}$,

$v_{27}(\xi \tau )=\pm \Pi ln(B)\alpha {coth}_B(\alpha \eta )e^{i\psi (\xi ,\tau )}$,

$v_{28}(\xi ,\tau )=\pm \Pi ln(B)\alpha (-{tanh}_B(2\alpha \eta )\pm i\sqrt{mn}\mathrm{sech}{}_B(2\alpha \eta )e^{i\psi (\xi ,\tau )}$,

$v_{29}(\xi ,\tau )=\pm \Pi ln(B)\alpha \left(-{coth}_B(2\alpha \eta )\pm \sqrt{mn}\mathrm{csch}{}_B(2\alpha \eta )\right)e^{i\psi (\xi ,\tau )}$,

$v_{30}(\xi ,\tau )=\pm \frac{1}{2}\Pi ln(B)\alpha \left(-{tanh}_B(\frac{\alpha }{2}\eta )\pm {coth}_B(\frac{\alpha }{2}\eta )\right)e^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{7}$. When ${\beta }^2=4\alpha \gamma$,

$v_{31}(\xi ,\tau )=\pm \left(\frac{\alpha \sqrt{4\alpha \gamma }ln(B)}{\Pi (b-\kappa {\theta }_3)}+\frac{\Pi }{2\alpha \gamma \eta }(1+\sqrt{\alpha \gamma }ln(B)\eta )\right)e^{i\psi (\xi ,\tau )}$,

$\mathbf{Case}\mathbf{8}$. When $\beta =\rho$, $\alpha =q\rho (q\ne 0)$ and $\gamma =0$,

$v_{32}(\xi ,\tau )=0$.

$\mathbf{Case}\mathbf{9}$. When $\beta =\gamma =0$,

$v_{33}(\xi ,\tau )=0$.

$\mathbf{Case}\mathbf{10}$. When $\beta =\alpha =0$,

$v_{34}(\xi ,\tau )=\pm \frac{\Pi }{\gamma \eta }$.

$\mathbf{Case}\mathbf{11}$. When $\alpha =0$ and $\beta \ne 0$,

$v_{35}(\xi ,\tau )=\pm \Pi ln(B)\frac{s\beta }{\gamma \left({cosh}_B(\beta \eta )-{sinh}_B(\beta \eta )+s\right)}{e}^{i\psi (\xi ,\tau )}$,

$v_{36}(\xi ,\tau )=\pm \Pi ln(B)\frac{\beta \left({sinh}_B(\beta \eta )+{cosh}_B(\beta \eta )\right)}{\gamma \left({sinh}_B(\beta \eta )+{cosh}_B(\beta \eta )+r\right)}{e}^{i\psi (\xi ,\tau )}$.

$\mathbf{Case}\mathbf{12}$. When $\beta =\rho$, $\gamma =m\rho (m\ne 0)$ and $\alpha =0$,

$v_{37}(\xi ,\tau )=\pm \Pi ln(B)\frac{s{B}^{\rho \eta }}{r-ms{B}^{\rho \eta }}{e}^{i\psi (\xi ,\tau )}$,

where $\eta$ and $\psi$ are given by Equations (8) and (9) for the beta derivative and Equations (10) and (11) for M-truncated derivative.

4. Comparison and Discussion

Remark: The constraint conditions for the presence of the acquired solutions are given in type as below.

Proposition: If $v(\xi ,\tau )$ is the solitary wave solution of the considered Equation (1), then the conditions for the presence of its solutions are $a>0$ and $b-\kappa {\theta }_3<0$.

Remark: It is essential to mention that the acquired structures of the model under study (Equation (1)) speak to the various sorts of soliton solutions. As $v_6$, $v_{16}$, and $v_{26}$ represent dark soliton solutions, $v_8$, $v_{18}$, and $v_{28}$ are the dark–bright soliton solutions, $v_{10}$, $v_{20}$, and $v_{30}$ show the dark–singular solutions, $v_9$, $v_{19}$, and $v_{29}$ represent a family of singular solutions of type 1 and 2, while $v_7$, $v_{17}$, and $v_{27}$ are reported as a singular solution of type 2.

Comparison: Now, different solutions are taken into account in the sense of beta and M-truncated derivatives and are depicted in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 with various $\mu$’s values.

Figure 1: Figure 1a shows 3D-graph with beta derivative definition, the second indicates a 2D-structure of $v_3(\xi ,\tau )$ when $\tau =1$ by applying two various definitions. It is noted that when $\tau =1$ both definitions show various graphs and overlapping exits in the definite range of the values of the independent variable $\xi$.

Figure 2: In this figure, the first graph represents its 3D graph for the M-truncated derivative, while Figure 2b shows its structures with the same definition of the derivative but for various values of the fractional parameter $\mu$ and $\tau =1$. Here, it is intriguing to take note that the curves have the same structure, and overlapping exists in the definite range of the values of the independent variable $\xi$ similar to the case forgiven in Figure 1b.

Figure 3: In Figure 3a there is 3D graph for the beta derivative for $v_1(\xi ,\tau )$ and Figure 3b shows its structure with the beta definition at $\tau =1$; a 2D structure is seen by utilizing two definitions. It is intriguing to take note that now curves have the same shape and overlapping also exists. In Figure 4, both 3D and 2D structures show the same behavior as in Figure 3. In Figure 5, Figure 5a indicates 3D and Figure 5b shows 2D structures for the beta derivative. What is mentioning here is that the soliton’s amplitude increases for an increase in the value of $\mu$. In Figure 6, represents 3D and 2D graphs, respectively, for the M-truncated derivative. What is also worth mentioning here is that the soliton’s amplitude decreases for an increase in the value of $\mu$.

5. Conclusions

In the present work, Equation (1), which describes the propagation of waves in negative-index metamaterials with the Bohm potential term, is taken in the sense of the beta and M-truncated derivatives. The Bohm potential term is accounted for the quantum phase behavior in the NLSE. This exploration elaborates on new families of solitons in negative optical metamaterials. To obtain the various type of solutions, the new extended algebraic method is considered. The considered technique also yielded new families including dark–bright, dark, dark–singular, and singular solutions of Type 1 and 2 of the governing equations. The consequences of this paper are of incredible interest in the fiber optic industry and designed by and large.