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*

Mei-Xia Chu,Bo Tian,Yu-Qiang Yuan,Ze Zhang,and He-Yuan Tian

Communications in Theoretical Physics ›› 2019, Vol. 71 ›› Issue (12) : 1393-1398.

PDF(1473 KB)
Welcome to visit Communications in Theoretical Physics, June 5, 2025
PDF(1473 KB)
Communications in Theoretical Physics ›› 2019, Vol. 71 ›› Issue (12) : 1393-1398. DOI: 10.1088/0253-6102/71/12/1393
Mathematical Physics

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*

Author information +
History +

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.

Key words

optical fiber / coupled cubic-quintic nonlinear SchrÕ / dinger equations / dark-dark solitons / Kadomtsev-Petviashvili hierarchy reduction

Cite this article

Download Citations
Chu Mei-Xia, Tian Bo, Yuan Yu-Qiang, et al. 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]. Communications in Theoretical Physics, 2019, 71(12): 1393-1398 https://doi.org/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]
iq1,z+k(z)q1,tt+φ(z)(|q1|2+|q2|2)q1+ς(z)(|q1|2+|q2|2)2q1  il(z)[(|q1|2+|q2|2)q1]t+ih(z)(q1q1,t+q2q2,t)q1=0,iq2,z+m(z)q2,tt+π(z)(|q1|2+|q2|2)q2+ω(z)(|q1|2+|q2|2)2q2  in(z)[(|q1|2+|q2|2)q2]t+ib(z)(q1q1,t+q2q2,t)q2=0,
(1)
where q1 and q2 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=1, k(z) and m(z) depict the group velocity dispersions, φ(z) and π(z) denote the nonlinearity parameters, ς(z) and ω(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
n(z)=h(z)=l(z)=4Rk(z)=4Rm(z)=ς(z)R=ω(z)R=4Rφ(z)v=4Rπ(z)v=b(z),
(2)
where R and v are the real constants.
Special cases of Eqs. (1) have been given as
(i) When k(z)=m(z)=1, φ(z)=π(z)=2, ς(z)=ω(z)=ρ12, l(z)=n(z)=h(z)=b(z)=2ρ1, Eqs. (1) have been reduced to the coupled cubic-quintic NLS equations,[36]
iq1,z+q1,tt+2(|q1|2+|q2|2)q1+ρ12(|q1|2+|q2|2)2q1  2iρ1[(|q1|2+|q2|2)q1]t+2iρ1(q1q1,t+q2q2,t)q1=0,iq2,z+q2,tt+2(|q1|2+|q2|2)q2+ρ12(|q1|2+|q2|2)2q2  2iρ1[(|q1|2+|q2|2)q2]t+2iρ1(q1q1,t+q2q2,t)q2=0.
(3)
Soliton solutions and DT for this case have been investigateded.[36]
(ii) When q1=q and q2=0, Eqs. (1) have been reduced to the variable coefficients Kundu-Eckhaus equation,[37]
iqz+k(z)qtt+φ(z)|q|2q+ς(z)|q|4qib(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,
q1=μ1gfeic1t+iF1(z)+2iRθ,q2=μ2hfeic2t+iF2(z)+2iRθ,
(5)
where
Fj(z)=14R[v(|μ1|2+|μ2|2)cj2]b(z)dz,  (j=1,2),θ=(|μ1|2|g|2f2+|μ2|2|h|2f2)dt+S(z),S(z)=[(|μ1|2+|μ2|2)zdt+ib(z)4RS(z)=×(u1,tu1u1u1,t+u2,tu2u2u2,t)]dz,uj=μjgfeicjt+iFj(z),
(6)
g and h represent the complex functions of z and t, f denotes a real function, cj's and μ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
Bracket argument to \\ must be a dimension
where the Hirota bilinear operators Dz,Dt are defined as[38]
Dzl1Dtl2(FH)=(zz)l1(tt)l2F(z,t)H(z,t)|z=z,t=t,
F(z,t) denotes an analytic function of z and t, H(z,t) represents an analytic function of the formal variables z and t, l1 and l2 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
q1=μ1gfeic1t+iF1(z)+2iRθ,q2=μ2hfeic2t+iF2(z)+2iRθ
(7)
where
f=|I+M|,g=|I+V|,h=|I+W|,
(8)
I represents the N×N identity matrix, N is an non-negative integer, M,V and W are the N×N matrices, of which elements are mkl,vkl and wkl, respectively, defined as
mkl=1pk+pleξk+ξl,vkl=pkic1pl+ic1mkl,wkl=pkic2pl+ic2mkl,
(9)
ξk=pkt+(i/4R)pk2b(z)dz+ξk0, (k,l=1,2,,N), pk's and ξk0's are the complex constants, where pk satisfies the constraint
|μ1|2(pkic1)(pk+ic1)+|μ2|2(pkic2)(pk+ic2)=2v.
(10)
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
q1=μ12eic1t+iF1(z)+2iRθq1=×[1+ϱ1+(ϱ11)tanh(ξ1+ξ1+s12)],q1=q2=μ22eic2t+iF2(z)+2iRθq2=×[1+ι1+(ι11)tanh(ξ1+ξ1+s12)],
(11)
where
ϱ1=p1ic1p1+ic1,  ι1=p1ic2p1+ic2,s1=ln(p1+p12),
ξ1=p1t+i4Rp12b(z)dz+ξk0,
v, c1, c2, μ1, μ2, and p1 satisfy
|μ1|2(p1ic1)(p1+ic1)+|μ2|2(p1ic2)(p1+ic2)=2v.
(12)
Based on Solutions (11), we require that s1>0, i.e., p1+p1>0, to guarantee Solutions (11) against the singularity.
When N=2, two dark-dark solutions via Solutions (7) can be presented as
q1=μ1gfeic1t+iF1(z)+2iRθ, q2=μ2hfeic2t+iF2(z)+2iRθ,
(13)
Where
f=1+eξ1+ξ1+s1+eξ2+ξ2+lns2+eξ1+ξ1+ξ2+ξ2+lns3,g=1+ϱ1eξ1+ξ1+s1+ϱ2eξ2+ξ2+lns2+ϱ1ϱ2eξ1+ξ1+ξ2+ξ2+lns3,h=1+iota1eξ1+ξ1+s1+iota2eξ2+ξ2+lns2+iota1iota2eξ1+ξ1+ξ2+ξ2+lns3,ϱ2=p2ic1p2+ic1,iota2=p2ic2p2+ic2,s2=1p2+p2,s3=|p1p2|2(p1+p1)(p2+p2)|p1+p2|2,ξk=pkt+i4Rpk2intb(z)dz+ξk0,(k=1,2),
(14)
v, c1, c2, μ1, μ2, p1, and p2 satisfy Constraint (10). Similarly, we require that s1>0 and s2>0, i.e., p1+p1>0 and p2+p2>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 Ak and velocity v of q1 and q2 can be given as
Ak=|μ1|2[2(2p1ickp1+ickp1+ickp1ick)1/2],v=|i4R(p1p1)b(z)|,
(15)
which is the magnitude of v.
Direction of v: When (i/4R)(p1p1)b(z)>0, the soliton propagates along the positive direction of t axis; When (i/4R)(p1p1)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 qk is relevant to b(z). When (p1ick)/(p1+ick)=1, amplitude of qk reaches the maximum, and minimum intensity drops to zero at the same time. The qk component is black with p1p1=2ick, or else it is grey.
Hereafter we mainly discuss the case of c1c2 for Eqs. (1). Based on Solutions (11), q1=χq2 (where χ is a constant) when c1=c2, and the soliton is black in both the q1 and q2 components; when c1c2, q1χq2.
More on the solitonic issues can be seen, e.g., in Refs. [40].
In Fig. 1, we choose c1c2, the q2 component is black when p1p1=2ic2. 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.
Fig. 1 (Color online) Linear-shaped dark-dark soliton via Solutions (11) with the parameters as p1=1/2+i/2, c1=0,c2=1/2, μ1=1, μ2=3/2, v=2/11, R=1, ξ1(0)=0 and b(z)=1/7.

Full size|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.

Full size|PPT slide

When Im(p1)Im(p2)<0, Figs. 3--4 show the head-on interactions between the two solitons of the same type with Im(pk) being the imaginary part of pk. The soliton S1 is black and S2 is grey in the q1 component when p1p1=2ic1, and just the opposite, in the q2 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 q1 and q2 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.
Fig. 3 (Color online) Head-on interaction between the parabolic-shaped dark-dark solitons via Solutions (13), with p1=1+i/2, p2=1i/2, c1=1/2, c2=1/2, μ1=μ2=(23)/3, v=1, R=2, ξ1(0)=ξ2(0)=0 and b(z)=3+2z.

Full size|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/2sin[(11/10)z].

Full size|PPT slide

Fig. 5 (Color online) Overtaking interaction between the parabolic-shaped two dark-dark soltions via Solutions (13), with p1=1+i/3, p2=1+i, c1=1, c2=1, μ1=(13/3), μ2=5/3, v=1, R=3/10, ξ1(0)=2, ξ2(0)=0, and b(z)=1/7+(7/10)z.

Full size|PPT slide

When the two dark-dark solitons have the same velocity, which means that Im(p1) = Im(p2), 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(p1)Im(p2)>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.
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).

Full size|PPT slide

Figures 7--8 show the interactions between the two solitons of different types when Im(q1)=0. Velocities of the linear-shaped dark soliton are zero. Fig. 7 illustrates the interaction between the linear-shaped dark soliton S1 and parabolic-shaped dark soliton S2. 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.
Fig. 7 (Color online) Interaction between the two dark-dark soltions via Solutions (13), with p1=1, p2=3/4+1/2i, c1=0, c2=1, μ1=(1338)/76, μ2=(3285)/38, v=1, R=1/2, ξ1(0)=ξ2(0)=0, and b(z)=z.

Full size|PPT slide

Fig. 8 (Color online) The same as Fig. 7 except that b(z)=2+2sin(3z/2).

Full size|PPT slide

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. 12, respectively. With b(z) as the linear and periodic functions and Im(p1)0, 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(p1)=0, interactions between the two dark-dark solitons of the different types have been displayed in Figs. 78, 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. 38 has indicated that those interactions are elastic.

References

[1]
P. A. Kumar, S. A. Krishno, R. A. B Siddiqur , and K. Abdul , IEEE Sens.J. 18(2018) 1558.
[2]
D. Kowal, W. Urbanczyk, P. Mergo , Sensors-Basel 18 ( 2018) .
[3]
L. Curilla, I. Astrauskas, A. Pugzlys , et al., Opt. Fiber Technol. 42(2018) 39.
[4]
M. S. M. Rajan , Z. Naturforsch. A 71 ( 2016) 751.
[5]
B. Tang, L. R. Bian, X. Zhou, K. Chen , Laser Phys. 28(2018) 015001.
[6]
M. Jerome and N. Alan, Nonlinear Optics , CRC, Boca Raton ( 2018).
[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.
[9]
M. Li, J. J. Shui, T. Xu , Appl. Math. Lett. 83 ( 2018) 110; X. Y. Gao, Appl. Math. Lett. 73(2017) 143.
[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.
[11]
W. Z. Bao The Nonlinear Schrodinger Equation and Applications in Bose-Einstein Condensation and Plasma Physics, World ScientiÕc, Singapore (2007).
[12]
Z. Z. Lan , Appl. Math. Lett. 86(2018) 243.
[13]
Z. Z. Lan , Opt. Laser Technol. 113(2019) 462.
[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.
[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.
[16]
S. Nazarenko, M. Onorato, D. Proment , Phys. Rev. A 90 ( 2014) 144.
[17]
V. E. Zakharov and A. B. Shabat , Sov. Phys. JETP 34 ( 1972) 62.
[18]
D. W. Zuo, H. X. Jia, D. M. Shan , l Superlattice. Microst. 101(2017) 522.
[19]
A. H. Kara, A. Biswas, Q. Zhou , et al., Optik. 174(2018) 195.
[20]
Y.S Kivsharand B. Luther-Davies, . Phys.Rep. 298(1998) 81.
[21]
T. Kanna, M. Vijayajayanthi, M. Lakshmanan , J. Phys. A 43( 2010) 434018.
[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.
[25]
Z. Y. Yan, J. Phys.Soc Jpn. .73(2004) 2397.
[26]
S. Gatz and J. Herrmann , J. Opt. Soc. Amer. B 11 ( 1991) 2296.
[27]
J. Herrmann , Opt. Commun. 87(1992) 161.
[28]
L. Min and X. Tao , Phys. Rev. E 91 ( 2015) 033202.
Abstract
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.
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.
[31]
S. Arshed, A. Biswas, Q. Zhou , et al., Optik. 172(2018) 826.
[32]
F. H. Qi, H. M. Ju, X. H. Meng, J. Li , Nonlinear Dyn. 77(2014) 1331.
[33]
F. H. Qi, X. G. Xu, P. Wang , Appl. Math. Lett. 54(2016) 60.
[34]
J. Chai, B. Tian, Y. F. Wang , et al., Phys. A 434 ( 2015) 296.
[35]
J. Chai, B. Tian, H. L. Zhen , et al., Phys. Scr. 90(2015) 045206.
[36]
F. H. Qi, B. Tian ,X. Lu, et al., Commun. Nonlinear Sci. Numer. Simul. 17(2012) 2372.
[37]
X.Y Xie . and Z. H. Yan, Appl. Math. Lett. 80(2018) 48.
[38]
R. Hirota , The Direct Method in Soliton Theory, Cambridge Univ. Press, Cambridge( 2004).
[39]
L. Liu, B. Tian, X. Y. Wu, Y. Q. Yuan, Z. Naturforsch. A . 72(2017) 9.
[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.
[41]
Y. Ohta, D. S. Wang, J. Yang , Stud. Appl. Math. 127(2014) 345.

Funding

*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)

RIGHTS & PERMISSIONS

Copyright reserved © 2019.
PDF(1473 KB)

415

Accesses

0

Citation

Detail

Sections
Recommended

/