Diverse chirped optical solitons and new complex traveling waves in nonlinear optical fibers
本站小编 Free考研考试/2022-01-02
Savaissou Nestor,1,6, Souleymanou Abbagari,2, Alphonse Houwe,1, Mustafa Inc,3,4, Gambo Betchewe,1, Serge Y Doka,51Department of Physics, Faculty of Science, the University of Maroua, P O Box 814, Maroua, Cameroon 2Department of Basic Science, Faculty of Mines and Petroleum Industries, University of Maroua, Maroua, Cameroon 3Firat University, Science Faculty, Department of Mathematics, 23119 Elazig, Turkey 4Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan, China 5Department of Physics, Faculty of Science, the University of Ngaoundere, P O Box 454, Cameroon
First author contact:6Author to whom any correspondence should be addressed. Received:2019-12-4Revised:2020-02-29Accepted:2020-03-3Online:2020-05-18
Abstract This paper studies chirped optical solitons in nonlinear optical fibers. However, we obtain diverse soliton solutions and new chirped bright and dark solitons, trigonometric function solutions and rational solutions by adopting two formal integration methods. The obtained results take into account the different conditions set on the parameters of the nonlinear ordinary differential equation of the new extended direct algebraic equation method. These results are more general compared to Hadi et al (2018 Optik172 545–53) and Yakada et al (2019 Optik 197 163108). Keywords:diverse chirped optical solitons;nonlinear optical fibers;new extended direct algebraic method
PDF (932KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite Cite this article Savaissou Nestor, Souleymanou Abbagari, Alphonse Houwe, Mustafa Inc, Gambo Betchewe, Serge Y Doka. Diverse chirped optical solitons and new complex traveling waves in nonlinear optical fibers. Communications in Theoretical Physics, 2020, 72(6): 065501- doi:10.1088/1572-9494/ab7ecd
1. Introduction
Today, research in the field of optical solitons has reached cruising speed. However, a lot of the attention has been concentrated on the investigatation of solitons of different suitable modeling equations describing sub-picosecond pulse propagation in optical fibers, with a plethora of mathematical integration to manipulate them. Moreover, the interesting aspect of optical solitons lies in optic fiber communication, which has revolutionized long-distance communication (i.e trans-continental and trans-oceanic). Hence, the transport of data has become fluid, and even reliable.
Over the last five years, a handful of categories of solitons like traveling-wave solutions and array breathers have been reported in the literature, such as variable sinh-Gaussian solitons, dark-combo solitons, soliton pairs, bright solitons, complexiton solitons, Gaussian-shaped solitons and singular solitons [1–12]. In addition, one of the environments that have shown themselves favorable to the propagation of these optical solitons is metamaterials.
Metamaterials are artificial materials which display unusual electromagnetic (EM) properties [13]. Because of these attractive properties, several studies on localized wave research have been successfully conducted. Furthermore, these localized solutions in nonlinear metamaterials usually take the forms of shift solitons, spatiotemporal solitons, temporal solitons and spatial solitons [13]. More recently, many works have been conducted in nonlinear metamaterials. For example, chirped soliton and periodic solutions of the generalized nonlinear Schrödinger equation, bright and dark soliton solutions for the (3+1)-dimensional coupled nonlinear Schrödinger equation with electric and magnetic fields, and the gray solitary wave of an extended nonlinear Schrödinger equation with third-order dispersion [3–5, 8, 13, 14] were investigated analytically.
The hard task today is to adopt some relevant methods to build exact traveling-wave solutions of these nonlinear Schrödinger equations which describe a short pulse propagation in nonlinear metamaterials. In this way, some mathematical tools were successfully used to investigate analytical solutions for partial differential equations (PDEs), such as the sine-Gordon expansion method [15], the auxiliary equation method [16], the first integral method [17], the new extended direct algebraic method [18, 19] and the generalized tanh method [18]. Thus, exact optical solitons in metamaterials with dispersive permeability have been reported by [13–15]. Recently, some authors have investigated optical solitons for higher-order nonlinear Schrödinger’s equations, describing a sub-picosecond pulse propagation in nonlinear metamaterials with cubic-quintic nonlinearity and fourth-order dispersion [12]. This work aims to study chirped optical solitons for the cubic nonlinear Schrödinger’s equation (CNLSE), describing a sub-picosecond pulse propagation in a nonlinear metamaterial with Kerr dispersion, low group velocity dispersion (GVD) and Kerr nonlinearity. The following governing equation will be adopted$\begin{eqnarray}{\rm{i}}{q}_{x}+\displaystyle \frac{1}{2}{\beta }_{2}{q}_{{tt}}+{\gamma }_{1}| q{| }^{2}q-{\rm{i}}{\alpha }_{1}{\left(| {q}^{2}| q\right)}_{t}-{\rm{i}}{\alpha }_{2}q{\left(| q{| }^{2}\right)}_{t}=0.\end{eqnarray}$Here, q(x, t) stands for a complex envelop amplitude, t represents the time (in the group velocity frame), x is the distance along the direction of propagation, while β2, α1 and α2 are, respectively, the GVD and the self-phase modulation (SPM) parameters. However, γ1 is the cubic nonlinearity term which compensates the dispersion and SPM terms to obtain the soliton.
Therefore, the paper is organized as follows: section 2 is devoted to the description of the methods. In section 3 the application of the methods will be carried out, and some graphical illustrations will follow. The last section will summarize the work.
2. The sight of methods
2.1. The new extended direct algebraic methods
The fundamental of the new extended direct algebraic method is given by the following steps [18–20].
Step 1: adopting the PDE in the following form$\begin{eqnarray}H(u,{u}_{t},{u}_{x},{u}_{{xx}},{u}_{{xt}},\ldots )=0,\end{eqnarray}$where u(x, t) is an unknown function, and H is a polynomial of u. Surmise the traveling-wave hypothesis as follows and then by adopting u(x, t)=F(ξ)$\begin{eqnarray}\xi =x-{vt},\end{eqnarray}$therefore, the PDE can turn into an ordinary differential equation$\begin{eqnarray}G(F,{F}^{{\prime} },{F}^{{\prime\prime} },{F}^{\prime\prime\prime },...)=0,\end{eqnarray}$and prime denotes the derivative with respect to ξ.
Step 2: considering that (4) has the solution in the following expression$\begin{eqnarray}F(\xi )=\sum _{j=0}^{N}{h}_{i}{Q}^{{\rm{i}}}(\xi ),\quad {g}_{n}\ne 0.\end{eqnarray}$Where gj(0≤j≤N) are constants to be determined later, and Q(ξ) satisfies the following ordinary differential equation (ODE)$\begin{eqnarray}{Q}^{{\prime} }(\xi )=\mathrm{ln}(A)(\lambda +\mu Q(\xi )+\sigma {Q}^{2}(\xi )),\end{eqnarray}$and $A\ne 0,1$. The solutions of the ODE (6) are:
Step 4: substituting (5) and (6) into (4), then collecting all the terms of Qj(ξ) to set to zero yields a system of algebraic equations.
Step 5: with the aid of software programme, MAPLE, the results of the system of algebraic equations can be obtained, and then we use the results of (6) to construct the exact solutions of (4).
However, the details of the generalized hyperbolic and trigonometric functions are given by [18, 19].
2.2. The extended auxiliary equation method
The pipe of the extended auxiliary equation method is described by the following step.
Step 1: the extended auxiliary equation method is based on the auxiliary ordinary equation method [12].
Therefore, an auxiliary ordinary equation (4) can be written as follows$\begin{eqnarray}{\left(\displaystyle \frac{\partial F}{\partial \xi }\right)}^{2}={h}_{0}{F}^{2}+{h}_{1}{F}^{3}+{h}_{2}{F}^{4},\end{eqnarray}$where h0, h1 and h2 are constants to be determined. Then, the solutions of (4) are given as follows:
To obtain the chirped optical solitons to the model, the following ansatz will be adopted like traveling-wave transformation$\begin{eqnarray}q(x,t)={e}^{{\rm{i}}[f(\xi )]}\phi (\xi )\quad \xi =x-{vt},\end{eqnarray}$φ(ξ) is the amplitude component of the soliton and f(ξ) is the chirp parameter, and $\delta \omega =-\tfrac{\partial }{\partial \xi }[f(\xi )]$ is the corresponding chirped optical solitons. Substituting (49) into (1) leads to$\begin{eqnarray}2{\phi }^{{\prime} }+2{\beta }_{2}{v}^{2}{\phi }^{{\prime} }{f}^{{\prime} }+{\beta }_{2}{v}^{2}\phi {f}^{{\prime\prime} }+6{\alpha }_{1}v{\phi }^{2}{\phi }^{{\prime} }+4{\alpha }_{2}v{\phi }^{2}{\phi }^{{\prime} }=0,\end{eqnarray}$and$\begin{eqnarray}2\phi {f}^{{\prime} }+2{\alpha }_{1}v{\phi }^{3}{f}^{{\prime} }+{\beta }_{2}{v}^{2}\phi {{f}^{{\prime} }}^{2}-2{\gamma }_{1}{\phi }^{3}-{\beta }_{2}{v}^{2}{\phi }^{{\prime\prime} }=0.\end{eqnarray}$From (51) the chirped component is recovered$\begin{eqnarray}\delta \omega =\displaystyle \frac{1}{{\beta }_{2}{v}^{2}}+\displaystyle \frac{1}{2{\beta }_{2}v}(3{\alpha }_{1}+2{\alpha }_{2}){\phi }^{2}.\end{eqnarray}$Consequently, (50) becomes$\begin{eqnarray}\begin{array}{l}-\displaystyle \frac{1}{{\beta }_{2}{v}^{2}}\phi -\left(\displaystyle \frac{2{\alpha }_{1}}{{\beta }_{2}v}+2{\gamma }_{1}\right){\phi }^{3}\\ \quad +\,\displaystyle \frac{1}{{\beta }_{2}}\left(\displaystyle \frac{{\alpha }_{1}}{4}(4{\alpha }_{2}-3{\alpha }_{1})+{\alpha }_{2}^{2}\right){\phi }^{5}-{\beta }_{2}{v}^{2}{\phi }^{{\prime\prime} }=0.\end{array}\end{eqnarray}$Next, multiplying (53) by ${\phi }^{{\prime} }$, then integrating once with the zero constant of integration, we obtain$\begin{eqnarray}\begin{array}{l}-\displaystyle \frac{1}{2{\beta }_{2}{v}^{2}}{\phi }^{2}-\displaystyle \frac{1}{4}\left(\displaystyle \frac{2{\alpha }_{1}}{{\beta }_{2}v}+2{\gamma }_{1}\right){\phi }^{4}\\ +\,\displaystyle \frac{1}{6{\beta }_{2}}\left(\displaystyle \frac{{\alpha }_{1}}{4}(4{\alpha }_{2}-3{\alpha }_{1})+{\alpha }_{2}^{2}\right){\phi }^{6}-\displaystyle \frac{1}{2}{\beta }_{2}{v}^{2}{{\phi }^{{\prime} }}^{2}=0.\end{array}\end{eqnarray}$To obtain a suitable and malleable form of (54), we surmise $\phi (\xi )=\sqrt{F(\xi )}$. Thence,$\begin{eqnarray}\begin{array}{l}-\displaystyle \frac{1}{2{\beta }_{2}{v}^{2}}{F}^{2}(\xi )-\displaystyle \frac{1}{4}\left(\displaystyle \frac{2{\alpha }_{1}}{{\beta }_{2}v}+2{\gamma }_{1}\right){F}^{3}(\xi )\\ +\,\displaystyle \frac{1}{6{\beta }_{2}}\left(\displaystyle \frac{{\alpha }_{1}}{4}(4{\alpha }_{2}-3{\alpha }_{1})+{\alpha }_{2}^{2}\right){F}^{4}(\xi )-\displaystyle \frac{1}{8}{\beta }_{2}{v}^{2}{{F}^{{\prime} }}^{2}(\xi )=0.\end{array}\end{eqnarray}$
4. Application of the new extended direct algebraic method to (55)
This section will apply the two methods described above to construct diverse exact traveling-wave solutions to the CNLSE describing a sub-picosecond pulse propagation in a nonlinear optic fiber.
Using the balance principle between F4 and ${{F}^{{\prime} }}^{2}(\xi )$, N=1 is gained. Thereafter, (5) can be expressed as follows$\begin{eqnarray}F(\xi )={h}_{0}+{h}_{1}Q(\xi ).\end{eqnarray}$Next, substituting (56) and (6) into (55) a system of algebraic equations in terms of ${(Q(\xi ))}^{j}$are recovered. Then, collecting all the coefficients of the obtained algebraic equations and setting to zero gives the set of results.
For μ=k and $\sigma ={mk}(m\ne 0)$ and λ=0, we obtain ${h}_{0}={\left(12{\theta }_{2}^{2}{\beta }_{2}^{2}\mathrm{ln}{\left(A\right)}^{2}\right)}^{\tfrac{1}{4}}$, ${h}_{1}={{mh}}_{0},v=\tfrac{\sqrt{-\tfrac{2{\theta }_{2}}{\sqrt{3}}(3{\alpha }_{1}+{\alpha }_{2})({\alpha }_{1}-2{\alpha }_{2})}}{{\beta }_{2}\mathrm{ln}(A)}$, α1=α1, α2=α2, β2=β2, ${\gamma }_{1}=\tfrac{3\left[2{\beta }_{2}\mathrm{ln}(A)+v\mathrm{ln}(A){\left(-12{\alpha }_{1}{\theta }_{2}^{2}\mathrm{ln}{\left(A\right)}^{2}{\sigma }^{2}\right)}^{\tfrac{1}{4}}\right]}{-12{\theta }_{2}^{2}{\sigma }^{2}{\beta }_{2}^{2}\mathrm{ln}{\left(A\right)}^{2}(3{\alpha }_{1}+2{\alpha }_{2})({\alpha }_{1}-2{\alpha }_{2})}$,$\begin{eqnarray}{q}_{35}(\xi )=\sqrt{{h}_{0}\left(-\displaystyle \frac{{{mpA}}^{k\xi }}{q-{{mpA}}^{k\xi }}\right)}\times {{\rm{e}}}^{{\rm{i}}[f(\xi )]},\end{eqnarray}$and lastly the chirped optical soliton to (1) is given$\begin{eqnarray}{\delta }_{35}=\displaystyle \frac{1}{{\beta }_{2}{v}^{2}}-\displaystyle \frac{1}{2{\beta }_{2}v}(3{\alpha }_{1}+2{\alpha }_{2})\displaystyle \frac{{{pA}}^{k\xi }}{q-{{mpA}}^{k\xi }}.\end{eqnarray}$Here, in this paper, we omit the generalized hyperbolic and triangular function which are given by [18, 19].
5. Application of the extended auxiliary equation method to (55)
To obtain the optical soliton solution, (55) can to be transformed as follows [12]:$\begin{eqnarray}{\left(\displaystyle \frac{\partial F}{\partial \xi }\right)}^{2}={h}_{0}{F}^{2}(\xi )+{h}_{1}{F}^{3}(\xi )+{h}_{2}{F}^{4}(\xi ).\end{eqnarray}$Where ${h}_{0}=-\tfrac{4}{{v}^{4}{\beta }_{2}^{2}}$, ${h}_{1}=\tfrac{4\left({\alpha }_{1}+{\gamma }_{1}{\beta }_{2}v\right)}{{\beta }_{2}^{2}{v}^{3}}$, ${h}_{2}=\tfrac{{\alpha }_{1}\left(4{\alpha }_{2}-3{\alpha }_{1}\right)+4{\alpha }_{1}^{2}}{3{v}^{2}{\beta }_{2}^{2}}$.
Figures 1–2 plot the spatiotemporal evolution of the dipole bright optical solitons and the corresponding chirped of $| {q}_{35}(x,t){| }^{2}$. Moreover, by using the extended algebraic method, for certain values of the parameters β2, α1 & α2, the dipole solitary waves are obtained (see figures 1, 2, 4). However, figure 3 depicts the bright optical soliton evolution, which corresponds to the spatiotemporal shift. Figure 5 depicts the dipole optical solitons which corresponds to the spatiotemporal shift. Three dimensional plots of these chirped optical solitons may be observed in figures 1–5. The results obtained are very important in optical fibre context. Therefore, the dark optical soliton can help to carry out the communication transmission shipped to where the renewal can be collected by the dark–bright soliton conversion.
Figure 1.
New window|Download| PPT slide Figure 1.Spatiotemporal plots of the combined bright optical solitons $| {q}_{35}(x,t){| }^{2}$ for β2=1, α1=−0.074, α2=0.0084, at (a) v=25.255 and (b) v=15.25 respectively.
Figure 2.
New window|Download| PPT slide Figure 2.Spatiotemporal plots of the chirped bright soliton $| \delta {\omega }_{35}{| }^{2}$ for β2=1, α1=−0.74, α2=0.25, γ1=−0.025 at (a) v=25.25, (b) v=15.25, (c) v=10.25, and (d) v=10.25, respectively.
Figure 3.
New window|Download| PPT slide Figure 3.Spatiotemporal plots of the combined bright optical solitons $| {q}_{36}(x,t){| }^{2}$ for β2=1, α1=−0.74, α2=0.25, γ1=−0.025 and (a) t = 0, (b) t = 5, (c) t = 10, (d) t = 15, (e) t = 20 at v = 5.25.
Figure 4.
New window|Download| PPT slide Figure 4.Spatiotemporal plots of dipole bright optical solitons $| {q}_{37}(x,t){| }^{2}$ for β2=1, α1=−0.0084, α2=0.0084, γ1=−0.025 at (a) v=10.25, (b) v=15.25, (c) v=20.25 and (d) v=25.25, respectively.
Figure 5.
New window|Download| PPT slide Figure 5.Spatiotemporal plots (2D) of dipole optical solitons $| {\delta }_{37}{| }^{2}$ for β2=1, α1=−0.0084, α2=0.0084, γ1=−0.025 at (a) v=5.25, (b) v=15.25, (c) v=25.25 and (d) v=30.25, respectively.
6. Conclusion
In this present work, through the governing equation, we describe a propagation of sub-picosecond pulses with cubic nonlinearity and SPM in optical fibers. We investigate the behavior of the new optical soliton solutions. To target the main objective of this paper, the traveling-wave transformation was adopted to obtain the chirped component and the resulting ordinary differential equation. To solve the obtained ODE, the new extended direct algebraic scheme and the extended algebraic equation have been used. In our results, chirped dark–bright soliton solutions, trigonometric function solutions, new complex traveling-wave solutions, some rational function solutions and optical dipole solitons also emerge. Also, cases (8), (9) and (10) of the new extended direct algebraic method did not satisfy the conditions set on nonlinear ordinary differential equation (NODE) parameters, therefore, no solution was obtained. We think that the obtained one-soliton solutions should be useful in the field of solitary-wave theory. The reader should look at [21, 22] for the nonlocal Schrödinger equation.
ZhouQZhuQLiuYBiswasABhrawyA HKhanK RMahmoodM FBelicM2014 Solitons in optical metamaterials with parabolic law nonlinearity and spatio-temporal dispersion 16 12211225
XuYSuarezPMilovicDKhanK RMahmoodM FBiswasABelicM2016 Raman solitons in nanoscale optical waveguides, with metamaterials, having polynomial law non-linearity 63 3237 DOI:10.1080/09500340.2016.1193240 [Cited within: 1]
LiminSZhenjunYShuminZXingliangL2019 Spiraling anomalous vortex beam arrays in strongly nonlocal nonlinear media 99063817 DOI:10.1103/PhysRevA.99.063817
TrikiHBiswasAMoshokoaS PBelicM2017 Dipole solitons in optical metamaterials with Kerr law nonlinearity 128 7176 DOI:10.1016/j.ijleo.2016.10.017
SongL MYangZ JPangZ GLiX LZhangS M2018 Interaction theory of mirror-symmetry soliton pairs in nonlocal nonlinear Schrödinger equation 90 4248 DOI:10.1016/j.aml.2018.10.008 [Cited within: 1]
SongL MYangZ JLiX LZhangS M2018 Controllable Gaussian-shaped soliton clusters in strongly nonlocal media 26 1918219198 DOI:10.1364/OE.26.019182
MaW X1993 Travelling wave solutions to a seventh order generalized KdV equation 180 221224
HadiRMirhosseini-AlizaminiS MEslamiMRezazadehMMirzazadehMAbbagariS2018 New optical solitons of nonlinear conformable fractional Schrodinger-Hirota equation 172 545553 DOI:10.1016/j.ijleo.2018.06.111 [Cited within: 5]
YakadaSDepelairBGamboBSergeY D2019 Miscellaneous new traveling waves in metamaterials by means of the new extended direct algebraic method 197163108 DOI:10.1016/j.ijleo.2019.163108 [Cited within: 3]
AbbagariSAlperKHadiRSerge PaulinT MAhmetB2019 Soliton solutions in different classes for the Kaup-Newell model equation 342050038 DOI:10.1142/S0217984920500384 [Cited within: 1]
SongL MYangZ JPangZ GLiX LZhangS M2019 Interaction theory of mirror-symmetry soliton pairs in nonlocal Schrodinger equation 90 4248 DOI:10.1016/j.aml.2018.10.008 [Cited within: 1]