Mei-Xia Chu, Bo Tian^,?, Yu-Qiang Yuan^?, Ze Zhang^?, and He-Yuan Tian^?State Key Laboratory of Information Photonics and Optical Communications, and School of Science, Beijing University of Posts and Telecommunications, Beijing 100876, China

Corresponding authors: ?E-mail:tian_bupt@163.com

Received:2019-08-15Online:2019-12-1

Fund supported:

*Supported by the National Natural Science Foundation of China under Grant .Nos. 11772017 Supported by the National Natural Science Foundation of China under Grant .Nos. 11805020 Supported by the National Natural Science Foundation of China under Grant .Nos. 11272023 Supported by the National Natural Science Foundation of China under Grant .Nos. 11471050 the Fund of State Key Laboratory of Information Photonics and Optical Communications Beijing University of Posts and Telecommunications, China IPOC.2017ZZ05 the Fundamental Research Funds for the Central Universities of China under Grant .No. 2011BUPTYB02

Abstract
Twin-core optical fibers are applied in such fields as the optical sensing and optical communication, and propagation of the pulses, Gauss beams and laser beams in the non-Kerr media is reported.Studied in this paper are the coupled cubic-quintic nonlinear Schr?dinger equations with variable coefficients, which describe the effects of quintic nonlinearity for the ultrashort optical pulse propagation in a twin-core optical fiber or non-Kerr medium. Based on the integrable conditions, bilinear forms are derived, and dark-dark soliton solutions can be constructed in terms of the Gramian via the Kadomtsev-Petviashvili hierarchy reduction. Propagation and interaction of the dark-dark solitons are presented and discussed through the graphic analysis. With different values of the delayed nonlinear response effect $b(z)$, where $z$ represents direction of the propagation, the linear- and parabolic-shaped one dark-dark soltions can be derived. Interactions between the parabolic- and periodic-shaped two dark-dark solitons are presented with $b(z)$ as the linear and periodic functions, respectively. Directions of velocities of the two dark-dark solitons vary with $z$ and the amplitudes of the solitons remain unchanged can be observed. Interactions between the two dark-dark solitons of different types are displayed, and we observe that the velocity of one soliton is zero and direction of the velocity of the other soliton vary with $z$. We find that those interactions are elastic.
Keywords： optical fiber;coupled cubic-quintic nonlinear Schr?;dinger equations;dark-dark solitons;Kadomtsev-Petviashvili hierarchy reduction


PDF (1473KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Mei-Xia Chu, Bo Tian, Yu-Qiang Yuan, Ze Zhang, and He-Yuan Tian. Bilinear Forms and Dark-Dark Solitons for the Coupled Cubic-Quintic Nonlinear Schr?dinger Equations with Variable Coefficients in a Twin-Core Optical Fiber or Non-Kerr Medium*. [J], 2019, 71(12): 1393-1398 doi:10.1088/0253-6102/71/12/1393

1 Introduction

Twin-core optical fibers have been applied in such fields as the optical sensing^[1] and optical communication,^[2] which have embodied in the optical fiber sensors,^[1] optical interferometers^[2] and optical switches.^[3] Propagation of the pulses,^[4] Gauss beams^[5] and laser beams^[6-8] in the non-Kerr media has attracted the researchers' attention.

Applications of the nonlinear SchrÕdinger (NLS)-type equations in nonlinear optics,^{[6, 9-10]} plasma physics,^[11-13] hydrodynamics^[14-15] and Bose-Einstein condensation^{[11, 16]} have attracted the researchers' attention. In nonlinear optical fibers, optical solitons, which are produced by the balance between the nonlinear effect and group velocity dispersion, have been described by the NLS-type equations.^{[6, 17-20]} With the consideration of the two solitons simultaneously propagate in the two-mode fiber or birefringent, the coupled NLS equations have been investigated.^[21-24]

When the intensity of the incident light field becomes stronger and non-Kerr nonlinearity plays a role, the NLS-type equations with higher-order nonlinearity terms have been used for describing the pulses in optical fibers.^[25-29] Combination of the cubic and quintic nonlinearity terms has allowed stable multidimensional structure exist of the cubic-quintic NLS equation.^[30-31] So as to describe the effects of quintic nonlinearity for the ultrashort optical pulse propagation in a twin-core optical fiber or non-Kerr medium, people have studied the following coupled cubic-quintic NLS equations with variable coefficients,^[32-35]

(1)$$ i{}q_{1,z}\!+\!k(z)q_{1,tt}\!+\!\varphi(z)(|q_1|^2\!+\!|q_2|^2)q_1\!+\!\varsigma(z)(|q_1|^2\!+\!|q_2|^2)^2q_1 \ \ -i{}l(z)[(|q_1|^2\!+\!|q_2|^2)q_1]_t\!+\! i h(z)(q_1^*q_{1,t}\!+\!q_2^*q_{2,t})q_1\!=\!0\,, i{}q_{2,z}\!+\!m(z)q_{2,tt}\!+\!\pi(z)(|q_1|^2\!+\!|q_2|^2)q_2\!+\!\omega(z)(|q_1|^2\!+\!|q_2|^2)^2q_2\! \ \ -i{}n(z)[(|q_1|^2\!+\!|q_2|^2)q_2]_t\!+\!i b(z)(q_1^*q_{1,t}\!+\!q_2^*q_{2,t})q_2\!=\!0\,, $$
where $q_1$ and $q_2$ represent the two components of the electromagnetic fields, the subscripts $z$ and $t$ denote the partial derivatives with respect to the scaled distance along the direction of the propagation and retarded time, respectively, $"*"$ represents the complex conjugate, $ i=\sqrt{-1}$, $k(z)$ and $m(z)$ depict the group velocity dispersions, $\varphi(z)$ and $\pi(z)$ denote the nonlinearity parameters, $\varsigma(z)$ and $\omega(z)$ are the saturations of the nonlinear refractive indexes, $l(z)$ and $n(z)$ represent the self-steepening effects, $h(z)$ and $b(z)$ are the delayed nonlinear response effects. To our knowledge, the conservation laws, Lax pair, Darboux transformation (DT),^[32] rogue wave solutions,^[33] mixed-type vector solitons,^[34] bright and dark solitons and B\"{a}cklund transformations^[35] for Eqs. (1) have been studied. Refs. ^{[32, 35]} have presented the integrable conditions for Eqs. (1), which can be rewritten as

(2) $$ n(z)=h(z)=l(z)=4Rk(z)=4Rm(z)=\frac{\varsigma(z)}{R}\\ =\frac{\omega(z)}{R}=\frac{4R\varphi(z)}{v}=\frac{4R\pi(z)}{v}=b(z)\,, $$
where $R$ and $v$ are the real constants.

Special cases of Eqs. (1) have been given as

(i) When $k(z)=m(z)=1$, $\varphi(z)=\pi(z)=2,$ $\varsigma(z)=\omega(z)=\rho_1^2,$ $l(z)=n(z)=h(z)=b(z)=2\rho_1,$ Eqs. (1) have been reduced to the coupled cubic-quintic NLS equations,^[36]

(3) $$ i{}q_{1,z}\!+\!q_{1,tt}\!+\!2(|q_1|^2+|q_2|^2)q_1\!+\!\rho_1^2(|q_1|^2\!+\!|q_2|^2)^2q_1\\ \ \ -2i\rho_1[(|q_1|^2\!+\!|q_2|^2)q_1]_t\!+\!2i\rho_1(q_1^*q_{1,t}\!+\!q_2^*q_{2,t})q_1\!=\!0\,,\\ i{}q_{2,z}\!+\!q_{2,tt}\!+\!2(|q_1|^2\!+\!|q_2|^2)q_2\!+\!\rho_1^2(|q_1|^2\!+\!|q_2|^2)^2q_2\\ \ \ -2i\rho_1[(|q_1|^2\!+\!|q_2|^2)q_2]_t\!+\!2i\rho_1(q_1^*q_{1,t}\!+\!q_2^*q_{2,t})q_2\!=\!0\,. $$
Soliton solutions and DT for this case have been investigateded.^[36]

(ii) When $q_1=q$ and $q_2=0$, Eqs. (1) have been reduced to the variable coefficients Kundu-Eckhaus equation,^[37]

$ i{}q_z\!+\!k(z)q_{tt}\!+\!\varphi(z)|q|^2q\!+\!\varsigma(z)|q|^4q\!-\!i{}b(z)(|q|^2)_tq\!=\!0\,,(4) $

which depicts the ultra-short femtosecond pulses propagate in optical fibers, where $q(z,t)$ are the electromagnetic wave, $b(z)$ is the nonlinear dispersion. Soliton interactions for this case have been studied via the Hirota method.^[37]

However, dark-dark soliton solutions for Eqs. (1) have not been studied. In Sec. 2 of this paper, under Integrable Conditions (2), bilinear forms will be constructed, and dark-dark soliton solutions for Eqs. (1) will be derived through the Kadomtsev-Petviashvili (KP) hierarchy reduction. In Sec. 3, we will graphically analyze the interaction and propagation of dark solitons. Conclusions will be written in Sec. 4.

2 Bilinear Forms and Dark Soliton Solutions for Eqs. (1)

In this section, we will consider the following variable transformations,

(5) $$ q_{1}=\mu_1\frac{g}{f}e^{i{}c_1t+i{}F_1(z)+2i{}R\theta}\,,\\q_{2}=\mu_2\frac{h}{f}e^{i{}c_2t+i{}F_2(z)+2i{}R\theta}\,, $$
where

(6) $ F_j(z)=\frac{1}{4R}[v(|\mu_1|^2\!+|\mu_2|^2)-c_j^2]\int b(z)\textrm{d}z\,,\ \ (j=1,2)\,,\\ \theta=\int \Big(|\mu_1|^2\frac{|g|^2}{f^2}+|\mu_2|^2\frac{|h|^2}{f^2}\Big)\textrm{d}t + S(z)\,, \\ \\ S(z)=\int\Big[-\int(|\mu_1|^2+|\mu_2|^2)_z d t+\frac{i b(z)}{4R} \\ \hphantom{S(z)=} \times (u_{1,t}u^*_1-u_{1}u^*_{1,t}+u_{2,t}u^*_2-u_{2}u^*_{2,t})\Big]d z\,, \\ u_j=\mu_j \frac{g}{f}e^{i c_j t}+i F_j(z)\,, $
$g$ and $h$ represent the complex functions of $z$ and $t$, $f$ denotes a real function, $c_j$'s and $\mu_j$'s are the real and complex constants, respectively. Substituting Variable Transformations (5) into Eqs. (1), we can derive the bilinear forms for Eqs. (1) with Integrable Conditions (2) as

$ \Big[i{}D_z+\frac{b(z)}{4R}D_t^2+i{}c_1\frac{b(z)}{2R}D_t\Big]g\cdot f=0\,,\\ \Big[i{}D_z+\frac{b(z)}{4R}D_t^2+i{}c_2\frac{b(z)}{2R}D_t\Big]h\cdot f=0\,,\\ [D_t^2\!+\!v(|\mu_1|^2\!+\!|\mu_2|^2)]f\cdot f\!=\!v(|\mu_1|^2|g|^2\!+\!|\mu_2|^2|h|^2), $

where the Hirota bilinear operators $D_z, D_t$ are defined as^[38]

$ D_{z}^{l_1}D_{t}^{l_2}(F\cdot H) \\ \;\;=\Big(\frac{\partial}{\partial z}-\frac{\partial}{\partial z^{\prime}}\Big)^{l_1}\Big(\frac{\partial}{\partial t}-\frac{\partial}{\partial t^{\prime}}\Big)^{l_2} F(z,t)H(z^{\prime},t^{\prime})| _{z^{\prime}=z, t^{\prime}=t}\,,$

$F(z, t)$ denotes an analytic function of $z$ and $t$, $H(z^{\prime}, t^{\prime})$ represents an analytic function of the formal variables $z^{\prime}$ and $t^{\prime}$, $l_1$ and $l_2$ are the non-negative integers.

Via Integrable Conditions (2), dark-dark soliton solutions in terms of the Gramian for Eqs. (1) can be constructed as

(7) $ q_1=\mu_1\frac{g}{f}{e}^{i c_1t+i F_1(z)+2i R\theta},\\ q_2=\mu_2\frac{h}{f}{e}^{i c_2t+i F_2(z)+2i R\theta}$
where

(8) $$ f=|I+M|,\quad g=|I+V|,\quad h=|I+W|,$$
$I$ represents the $N \times N$ identity matrix, $N$ is an non-negative integer, $M, V$ and $W$ are the $N\times N$ matrices, of which elements are $m_{kl}, v_{kl}$ and $w_{kl}$, respectively, defined as

(9) $$ m_{kl}=\frac{1}{p_k+p_l^*}e^{\xi_k+\xi_l^*},\quad v_{kl}=-\frac{p_k-i{}c_1}{p_l^*+i{}c_1}m_{kl},\\ w_{kl}=-\frac{p_k-i{}c_2}{p_l^*+i{}c_2}m_{kl},$$
$\xi_k=p_kt+({i}/{4R})p_k^2\int b(z)\textrm{d}z+\xi_{k0}$, $(k, l=1,2,\ldots,N)$, $p_k$'s and $\xi_{k0}$'s are the complex constants, where $p_k$ satisfies the constraint

(10) $$ \frac{|\mu_1|^2}{(p_k-i{}c_1)(p_k^*+i{}c_1)}+\frac{|\mu_2|^2}{(p_k-i{}c_2)(p_k^*+i{}c_2)}=-\frac{2}{v}. $$
It is worth noting that the proof of the above process is similar to that in Ref. ^[39]. It is observed that $v$ should be negative according to Constraint (10).

When $N=1$, one dark-dark soliton solutions via Solutions (7) can be derived as

(11) $ q_1=\frac{\mu_1}{2}e^{i{}c_1t+i{}F_1(z)+2i{}R\theta} \hphantom{q_1=} \times \Big[1+\varrho_1+(\varrho_1-1)\tanh\Big(\frac{\xi_1+\xi_1^*+ s_1}{2}\Big)\Big], \hphantom{q_1=} q_2=\frac{\mu_2}{2}e^{i{}c_2t+i{}F_2(z)+2i{}R\theta}\hphantom{q_2=}\times \Big[1+\iota_1+(\iota_1-1)\tanh\Big(\frac{\xi_1+\xi_1^*+ s_1}{2}\Big)\Big], $
where

$ \varrho_1\!=\!-\frac{p_1-i{}c_1}{p_1^*+i{}c_1},\ \ \iota_1=-\frac{p_1-i{}c_2}{p_1^*+i{}c_2},s_1=\ln \Big(\frac{p_1+p_1^*}{2}\Big), $

$ \xi_1=p_1t+\frac{i}{4R}p_1^2\int b(z)\textrm{d}z+\xi_{k0}$,

$v$, $c_1$, $c_2$, $\mu_1$, $\mu_2$, and $p_1$ satisfy

(12) $$ \frac{|\mu_1|^2}{(p_1-i{}c_1)(p_1^*+i{}c_1)}+\frac{|\mu_2|^2}{(p_1-i{}c_2)(p_1^*+i{}c_2)}=-\frac{2}{v}. $$
Based on Solutions (11), we require that $s_1>0$, i.e., $p_1+p_1^*>0$, to guarantee Solutions (11) against the singularity.

When $N=2$, two dark-dark solutions via Solutions (7) can be presented as

(13) $$ q_1=\mu_1\frac{g}{f}e^{i c_1t+i F_1(z)+2i R\theta}\,,\\\ \\ q_2=\mu_2\frac{h}{f}e^{i c_2t+i F_2(z)+2i R\theta}\,, $$
Where

(14)$$ \displaylines{f\!\!=\!\!1\!+\!e^{\xi_1\!+\!\xi_1^*\!+\!s_1}\!+\!e^{\xi_2\!+\!\xi_2^*\!+\!\ln s_2}\!+\!e^{\xi_1\!+\!\xi_1^*\!+\!\xi_2\!+\!\xi_2^*\!+\!\ln s_3},\hfill{}\cr g\!\!=\!\!1\!+\!\varrho_1e^{\xi_1\!+\!\xi_1^*\!+\!s_1}\!+\!\varrho_2e^{\xi_2\!+\!\xi_2^*\!+\!\ln s_2}\!+\!\varrho_1\varrho_2e^{\xi_1\!+\!\xi_1^*\!+\!\xi_2\!+\!\xi_2^*\!+\!\ln s_3}\,,\hfill{}\cr h\!\!=\!\!1\!+\!iota_1e^{\xi_1\!+\!\xi_1^*\!+\!s_1}\!+\!iota_2e^{\xi_2\!+\!\xi_2^*\!+\!\ln s_2}\!+\!iota_1iota_2e^{\xi_1\!+\!\xi_1^*\!+\!\xi_2\!+\!\xi_2^*\!+\!\ln s_3}\,,\hfill{}\cr \varrho_2=-\frac{p_2-i{}c_1}{p_2^*+i{}c_1},\quad iota_2=-\frac{p_2-i{}c_2}{p_2^*+i{}c_2},\hfill{}\cr s_2=\frac{1}{p_2+p_2^*},\quad s_3=\frac{|p_1-p_2|^2}{(p_1+p_1^*)(p_2+p_2^*)|p_1+p_2^*|^2},\hfill{}\cr \xi_k=p_kt+\frac{i}{4R}p_k^2int b(z)\textrm{d}z+\xi_{k0}\,, \quad (k=1, 2)\,,\hfill{}\cr}$$
$v$, $c_1$, $c_2$, $\mu_1$, $\mu_2$, $p_1$, and $p_2$ satisfy Constraint (10). Similarly, we require that $s_1>0$ and $s_2>0$, i.e., $p_1+p_1^*>0$ and $p_2+p_2^*>0$, to guarantee Solutions (11) against the singularity.

3 Discussions

Based on Solutions (11) and (13), dark-dark solitons for Eqs. (1) will be displayed graphically. Via Solutions (11), amplitude $A_k$ and velocity $\vec{v}$ of $q_1$ and $q_2$ can be given as

(15) $$ A_k=\frac{|\mu_1|}{2}\Big[2-\Big(2-\frac{p_1-i{}c_k}{p_1^*+i{}c_k}-\frac{p_1^*+i{}c_k}{p_1-i{}c_k}\Big)^{{1}/{2}}\Big],\\ v=\Big|-\frac{i}{4R}(p_1-p_1^*)b(z)\Big|\,, $$
which is the magnitude of $\vec{v}$.

Direction of $\vec{v}$: When $-({i}/{4R})(p_1-p_1^*)b(z) > 0$, the soliton propagates along the positive direction of $t$ axis; When $-({i}/{4R})(p_1-p_1^*)b(z) < 0$, the soliton propagates along the negative direction of $t$ axis.

We can see that the amplitude is irrelevant to $b(z)$, the velocity of $q_k$ is relevant to $b(z)$. When $({p_1-i{}c_k})/({p_1^*+i{}c_k})=1$, amplitude of $q_k$ reaches the maximum, and minimum intensity drops to zero at the same time. The $q_k$ component is black with $p_1-p_1^*=2i{}c_k$, or else it is grey.

Hereafter we mainly discuss the case of $c_1\neq c_2$ for Eqs. (1). Based on Solutions (11), $q_1=\chi q_2$ (where $\chi$ is a constant) when $c_1=c_2$, and the soliton is black in both the $q_1$ and $q_2$ components; when $c_1\neq c_2$, $q_1\neq\chi q_2$.

More on the solitonic issues can be seen, e.g., in Refs. ^[40].

In Fig. 1, we choose $c_1\neq c_2$, the $q_2$ component is black when $p_1-p_1^*=2i{}c_2$. When the variable coefficient $b(z)$ is a constant, velocity of the one dark-dark soliton remains unchanged when the soliton propagates, as shown in Fig. 1.

When $b(z)=2z$, Fig. 2 depicts the parabolic-shaped dark-dark soliton, of which the direction of the velocity changes with $z$.

New window|Download| PPT slide
Fig. 1(Color online) Linear-shaped dark-dark soliton via Solutions (11) with the parameters as $p_1={1}/{2}+{i}/{2},$ $c_1=0, c_2={1}/{2},$ $\mu_1=1,$ $\mu_2={3}/{2},$ $v=-{2}/{11}$, $R=1$, $\xi_1^{(0)}=0$ and $b(z)={1}/{7}$.

New window|Download| PPT slide
Fig. 2(Color online) Parabolic-shaped dark-dark soliton via Solutions (11) with the same parameters as those in Fig. 1 except that $b(z)=2z$.



When Im$(p_1)$Im$(p_2)<0$, Figs. 3--4 show the head-on interactions between the two solitons of the same type with Im$(p_k)$ being the imaginary part of $p_k$. The soliton $S_1$ is black and $S_2$ is grey in the $q_1$ component when $p_1-p_1^*=2i{}c_1$, and just the opposite, in the $q_2$ component. When $b(z)={8}/{5}+z$, Fig. 3 depicts interaction of the parabolic-shaped two dark-dark solitons, of which the directions of the velocities vary with $z$.

The amplitudes in both the $q_1$ and $q_2$ components remain unchanged during the propagation. When $b(z)={5}/{2}+({5}/{2})\sin[({11}/{10})z]$, interaction of the periodic-shaped two dark-dark solitons is displayed in Fig. 4. Amplitudes of the two solitons keep invariant during those interactions except that some phase shifts vary, implying that those interactions are elastic.

New window|Download| PPT slide
Fig. 3(Color online) Head-on interaction between the parabolic-shaped dark-dark solitons via Solutions (13), with $p_1=1+{i}/{2},$ $p_2=1-{i}/{2}$, $c_1={1}/{2},$ $c_2=-{1}/{2},$ $\mu_1=\mu_2=({2\sqrt{3}})/{3},$ $v=-1$, $R=2$, $\xi_1^{(0)}=\xi_2^{(0)}=0$ and $b(z)=3+2z$.

New window|Download| PPT slide
Fig. 4(Color online) Head-on interaction between the periodic-shaped dark-dark solitons via Solutions (13) with the same parameters as those in Fig. 3 except that $b(z)={5}/{2}+{5}/{2}\sin[({11}/{10})z]$.

New window|Download| PPT slide
Fig. 5(Color online) Overtaking interaction between the parabolic-shaped two dark-dark soltions via Solutions (13), with $p_1=1+{i}/{3}$, $p_2=1+i$, $c_1=1$, $c_2=-1$, $\mu_1=({\sqrt{13}}/{3})$, $\mu_2={5}/{3}$, $v=-1$, $R={3}/{10}$, $\xi_1^{(0)}=2$, $\xi_2^{(0)}=0$, and $b(z)={1}/{7}+({7}/{10})z$.



When the two dark-dark solitons have the same velocity, which means that Im$(p_1)$ = Im$(p_2)$, we derive the bound-state dark-dark solitons based on Solutions (13). Since the coefficients of the same nonlinear terms for Eqs. (1) are the same, similar to the analysis in Ref. ^[41], we cannot obtain the bound-state dark-dark solitons. When Im$(p_1)$Im$(p_2)>0$, Figs. 5--6 illustrate the overtaking interactions of the two solitons with the same type. In Fig. 5, parabolic-shaped two dark-dark solitons are displayed. Directions of the velocities of the two solitons change with $z$ simultaneously. Fig. 6 displays the periodic-shaped two dark-dark solitons, of which the directions of the velocities vary periodically. As we can see, the amplitudes of the above two dark-dark solitons are unaffected by those interactions, namely, those interactions are elastic.

New window|Download| PPT slide
Fig. 6(Color online) Overtaking interaction between the periodic-shaped two dark-dark soltions via Solutions (13) with the same parameters as those in Fig. 5 except that $b(z)=1+\sin({z}/{2})$.



Figures 7--8 show the interactions between the two solitons of different types when Im$(q_1)=0$. Velocities of the linear-shaped dark soliton are zero. Fig. 7 illustrates the interaction between the linear-shaped dark soliton $S_1$ and parabolic-shaped dark soliton $S_2$. Interaction between the linear-shaped dark soliton and periodic-shaped dark soliton is displayed in Fig. 8. Similar to the above, those interactions are elastic.

New window|Download| PPT slide
Fig. 7(Color online) Interaction between the two dark-dark soltions via Solutions (13), with $p_1=1,$ $p_2={3}/{4}+{1}/{2}i$, $c_1=0$, $c_2=-1$, $\mu_1=({13\sqrt{38}})/{76}$, $\mu_2=({3\sqrt{285}})/{38}$, $v=-1$, $R={1}/{2}$, $\xi_1^{(0)}=\xi_2^{(0)}=0$, and $b(z)=z$.

New window|Download| PPT slide
Fig. 8(Color online) The same as Fig. 7 except that $b(z)=-2+2\sin({3z}/{2})$.

4 Conclusions

Twin-core optical fibers have been applied in such fields as the optical sensing and optical communication, and propagation of the pulses, Gauss beams and laser beams in the non-Kerr media has been reported. The coupled cubic-quintic NLS equations with variable coefficients have been investigated, i.e., Eqs. (1), which describe the effects of quintic nonlinearity for the ultrashort optical pulse propagation in a twin-core optical fiber or non-Kerr medium. Through the KP hierarchy reduction, Dark-dark Soliton Solutions (7) have been constructed in terms of the Gramian and bilinear forms have been derived with Integrable Conditions (2). When the delayed nonlinear response effect $b(z)$ is the constant and linear functions, the linear- and parabolic-shaped disproportional one dark-dark soltions have been displayed in Figs. 1$-$2, respectively. With $b(z)$ as the linear and periodic functions and Im$(p_1)\neq0$, head-on interactions between the two dark-dark solitons of the parabolic- and periodic-shaped have been presented in Figs. 3--4, respectively, overtaking interactions between the two dark-dark solitons of the same type have been presented in Figs. 5--6, while the velocities of the two dark-dark solitons have altered with $z$ changing and the amplitudes of the two dark-dark solitons have been seen to remain unchanged. Besides, when Im$(p_1)=0$, interactions between the two dark-dark solitons of the different types have been displayed in Figs. 7$-$8, and we have seen that the velocities of the one soliton is zero and direction of the velocity of the other soliton vary with $z$. Graphically analyzing the interactions between the two solitons in Figs. 3$-$8 has indicated that those interactions are elastic.

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1] P. A. Kumar, S. A. Krishno, R. A. B Siddiqur , and K. Abdul , IEEE Sens.J. 18(2018) 1558.
[Cited within: 2]

[2] D. Kowal, W. Urbanczyk, P. Mergo , Sensors-Basel 18 ( 2018) .
[Cited within: 2]

[3] L. Curilla, I. Astrauskas, A. Pugzlys , et al., Opt. Fiber Technol. 42(2018) 39.
DOI:10.1016/j.yofte.2018.02.020URL [Cited within: 1]

[4] M. S. M. Rajan , Z. Naturforsch. A 71 ( 2016) 751.
[Cited within: 1]

[5] B. Tang, L. R. Bian, X. Zhou, K. Chen , Laser Phys. 28(2018) 015001.
DOI:10.1088/1555-6611/aa9628URL [Cited within: 1]

[6] M. Jerome and N. Alan, Nonlinear Optics , CRC, Boca Raton ( 2018).
[Cited within: 3]

[7] W. J. Liu, M. L. Liu, B. Liu , et al., Opt. Express 27
( 2019) 6689.


[8] M. L. Liu, Y. Y. Ouyang, H. R. Hou , et al., Chin. Opt.Lett. 17(2019) 020006.
[Cited within: 1]

[9] M. Li, J. J. Shui, T. Xu , Appl. Math. Lett. 83
( 2018) 110; X. Y. Gao, Appl. Math. Lett. 73(2017) 143.
DOI:10.1016/j.aml.2017.03.020URL [Cited within: 1]

[10] M. L. Liu, W. J. Liu, Z. Y.Wei , J. Lightwave Technol. 37( 2019) 3100;
[J]. J. Su and Y. T. Gao, Wave Random Comples 28 (2019) 708; X. Y. Xie and G. Q. Meng, Eur Phys. J. Plus. 134(2019) 359.
[Cited within: 1]

[11] W. Z. Bao The Nonlinear Schrodinger Equation and Applications in Bose-Einstein Condensation and Plasma Physics, World Scienti?c, Singapore (2007).
[Cited within: 2]

[12] Z. Z. Lan , Appl. Math. Lett. 86(2018) 243.
DOI:10.1016/j.aml.2018.05.014URL

[13] Z. Z. Lan , Opt. Laser Technol. 113(2019) 462.
DOI:10.1016/j.optlastec.2018.12.040URL [Cited within: 1]

[14] E. Heifetz and E. Cohen , Found. Phys. 45(2015) 1514; T. T. Jia, Y. T. Gao, Y. J. Feng, et al
., Nonlinear Dyn. 96(2019) 229.
DOI:10.1007/s11071-019-04786-0URL [Cited within: 1]

[15] W. X. Ma , J. Geom. Phys. 133( 2018) 10; T. T. Jia, X. T. Gao, G. F. Deng, et al., Nonlinear Dyn. 98 (2019) 269; D. Y. Liu and W. R. Sun, Appl. Math. Lett. 84(2018) 63.
[Cited within: 1]

[16] S. Nazarenko, M. Onorato, D. Proment , Phys. Rev. A 90 ( 2014) 144.
[Cited within: 1]

[17] V. E. Zakharov and A. B. Shabat , Sov. Phys. JETP 34 ( 1972) 62.
[Cited within: 1]

[18] D. W. Zuo, H. X. Jia, D. M. Shan , l Superlattice. Microst. 101(2017) 522.
DOI:10.1016/j.spmi.2016.11.051URL

[19] A. H. Kara, A. Biswas, Q. Zhou , et al., Optik. 174(2018) 195.
DOI:10.1016/j.ijleo.2018.08.067URL

[20] Y.S Kivsharand B. Luther-Davies, . Phys.Rep. 298(1998) 81.
[Cited within: 1]

[21] T. Kanna, M. Vijayajayanthi, M. Lakshmanan , J. Phys. A 43( 2010) 434018.
[Cited within: 1]

[22] S. V. Manakov , Sov. Phys. JETP 38 (1974) 248.


[23] V. I. Kruglov, A. C. Peacock, J. D. Harvey , Phys. Rev. E 71(2015) 056619.


[24] X. Y. Xie and G. Q. Meng , Chaos,Solitons and Fractals 107 ( 2018) 143.
[Cited within: 1]

[25] Z. Y. Yan, J. Phys.Soc Jpn. .73(2004) 2397.
[Cited within: 1]

[26] S. Gatz and J. Herrmann , J. Opt. Soc. Amer. B 11 ( 1991) 2296.


[27] J. Herrmann , Opt. Commun. 87(1992) 161.
DOI:10.1016/0030-4018(92)90006-DURL

[28] L. Min and X. Tao , Phys. Rev. E 91 ( 2015) 033202.
DOI:10.1103/PhysRevE.91.033202URLPMID:25871233
Via the Nth Darboux transformation, a chain of nonsingular localized-wave solutions is derived for a nonlocal nonlinear Schr?dinger equation with the self-induced parity-time (PT) -symmetric potential. It is found that the Nth iterated solution in general exhibits a variety of elastic interactions among 2N solitons on a continuous-wave background and each interacting soliton could be the dark or antidark type. The interactions with an arbitrary odd number of solitons can also be obtained under different degenerate conditions. With N=1 and 2, the two-soliton and four-soliton interactions and their various degenerate cases are discussed in the asymptotic analysis. Numerical simulations are performed to support the analytical results, and the stability analysis indicates that the PT-symmetry breaking can also destroy the stability of the soliton interactions.

[29] W. J. Liu, C. Y. Yang, M. L. Liu , et al., Phys. Rev. E 96 ( 2017) 042201.
DOI:10.1103/PhysRevE.96.042201URLPMID:29347524 [Cited within: 1]
The interactions of multiple solitons show different properties with two-soliton interactions. For the difficulty of deriving multiple soliton solutions, it is rare to study multiple soliton interactions analytically. In this paper, three-soliton interactions in inhomogeneous optical fibers, which are described by the variable coefficient Hirota equation, are investigated. Via the Hirota bilinear method and symbolic computation, analytic three-soliton solutions are obtained. According to the obtained solutions, properties and features of three-soliton interactions are discussed by changing the third-order dispersion (TOD) and other relevant coefficients, and some plentiful structure of three-soliton interactions are presented for the first time. The influences of TOD on the intensity and propagation distance of solitons are described, which can be used to realize the soliton control. Besides, the method that can achieve the phase reverse of solitons is suggested, and bound states of three solitons are observed, which have potential applications in the mode-locked fiber lasers. Furthermore, comparing to two-soliton interactions, a novel phenomenon of three-soliton interactions with a strong phase shift at x=0 is revealed, which is potentially useful for optical logic switches.

[30] N. Fabien, M. Alidou, C. K. Timoléeon , Opt. Commun. 275(2007) 421.
DOI:10.1016/j.optcom.2007.02.036URL [Cited within: 1]

[31] S. Arshed, A. Biswas, Q. Zhou , et al., Optik. 172(2018) 826.
DOI:10.1016/j.ijleo.2018.07.086URL [Cited within: 1]

[32] F. H. Qi, H. M. Ju, X. H. Meng, J. Li , Nonlinear Dyn. 77(2014) 1331.
DOI:10.1007/s11071-014-1382-5URL [Cited within: 3]

[33] F. H. Qi, X. G. Xu, P. Wang , Appl. Math. Lett. 54(2016) 60.
DOI:10.1016/j.aml.2015.11.001URL [Cited within: 1]

[34] J. Chai, B. Tian, Y. F. Wang , et al., Phys. A 434 ( 2015) 296.
[Cited within: 1]

[35] J. Chai, B. Tian, H. L. Zhen , et al., Phys. Scr. 90(2015) 045206.
DOI:10.1088/0031-8949/90/4/045206URL [Cited within: 3]

[36] F. H. Qi, B. Tian ,X. Lu, et al., Commun. Nonlinear Sci. Numer. Simul. 17(2012) 2372.
DOI:10.1016/j.cnsns.2011.10.017URL [Cited within: 2]

[37] X.Y Xie . and Z. H. Yan, Appl. Math. Lett. 80(2018) 48.
DOI:10.1016/j.aml.2018.01.003URL [Cited within: 2]

[38] R. Hirota , The Direct Method in Soliton Theory, Cambridge Univ. Press, Cambridge( 2004).
[Cited within: 1]

[39] L. Liu, B. Tian, X. Y. Wu, Y. Q. Yuan, Z. Naturforsch. A . 72(2017) 9.
[Cited within: 1]

[40] J. J. Su, Y. T. Gao and C. C. Ding , Appl. Math. Lett. 88 ( 2019) 201; X. Y. Gao, Appl. Math. Lett. 91 (2019) 165; C. C. Ding, Y. T. Gao, L. Hu, et al., Eur. Phys. J. Plus 133 (2019) 406; C. C. Ding, Y. T. Gao, L. Q. Li. Chaos, Solitons and Fractals 120 (2019) 259; G. F. Deng, Y. T. Gao, and X. Y. Gao, Wave. Random Complex 28 (2019) 468; G. F. Deng, Y. T. Gao, J. J. Su, et al., Appl. Math. Lett. 98 (2019) 177; Y. J. Feng, Y. T. Gao, X. Yu, Nonlinear Dyn. 91 (2018) 29; Y. J. Feng, Y. T. Gao, L. Q. Li, et al.,
Applicable Analysis (2019) DOI: 10.1080/00036811.2019.1652734.
[Cited within: 1]

[41] Y. Ohta, D. S. Wang, J. Yang , Stud. Appl. Math. 127(2014) 345.
DOI:10.1111/sapm.2011.127.issue-4URL [Cited within: 1]