Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Wronskian-type determinant solutions of the nonlocal derivative nonlinear Schrödinger equation

  • A nonlocal derivative nonlinear Schrödinger (DNLS) equation is analytically studied in this paper. By constructing Darboux transformations (DTs) of arbitrary order, new determinant solutions of the nonlocal DNLS equation in the form of Wronskian-type are derived from both zero and nonzero seed solutions. Periodic solitons are obtained with different parameter choices. When one eigenvalue tends to another one, generalized DTs are constructed, leading to rogue waves. Due to complex parametric constraints, the derived solutions may have singularities. Despite this, the work presented in this paper can still provide a valuable reference for the study of nonlocal integrable systems.

    Citation: Dexin Meng. Wronskian-type determinant solutions of the nonlocal derivative nonlinear Schrödinger equation[J]. AIMS Mathematics, 2025, 10(2): 2652-2667. doi: 10.3934/math.2025124

    Related Papers:

    [1] Li Cheng, Yi Zhang, Ying-Wu Hu . Linear superposition and interaction of Wronskian solutions to an extended (2+1)-dimensional KdV equation. AIMS Mathematics, 2023, 8(7): 16906-16925. doi: 10.3934/math.2023864
    [2] Yaqiong Liu, Yunting Li, Qiuping Liao, Yunhui Yi . Classification of nonnegative solutions to fractional Schrödinger-Hatree-Maxwell type system. AIMS Mathematics, 2021, 6(12): 13665-13688. doi: 10.3934/math.2021794
    [3] Yong-Chao Zhang . Least energy solutions to a class of nonlocal Schrödinger equations. AIMS Mathematics, 2024, 9(8): 20763-20772. doi: 10.3934/math.20241009
    [4] Amjad Hussain, Muhammad Khubaib Zia, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar, Ilyas Khan . Lie analysis, conserved vectors, nonlinear self-adjoint classification and exact solutions of generalized (N+1)-dimensional nonlinear Boussinesq equation. AIMS Mathematics, 2022, 7(7): 13139-13168. doi: 10.3934/math.2022725
    [5] Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, M. Elsaid Ramadan, Islam Samir, Soliman Alkhatib . Influence of the β-fractional derivative on optical soliton solutions of the pure-quartic nonlinear Schrödinger equation with weak nonlocality. AIMS Mathematics, 2025, 10(3): 7489-7508. doi: 10.3934/math.2025344
    [6] Zh. Myrzakulova, Z. Zakariyeva, K. Suleimenov, U. Uralbekova, K. Yesmakhanova . Explicit solutions of nonlocal reverse-time Hirota-Maxwell-Bloch system. AIMS Mathematics, 2024, 9(12): 35004-35015. doi: 10.3934/math.20241666
    [7] Wen-Xin Zhang, Yaqing Liu . Solitary wave solutions and integrability for generalized nonlocal complex modified Korteweg-de Vries (cmKdV) equations. AIMS Mathematics, 2021, 6(10): 11046-11075. doi: 10.3934/math.2021641
    [8] Weerawat Sudsutad, Wicharn Lewkeeratiyutkul, Chatthai Thaiprayoon, Jutarat Kongson . Existence and stability results for impulsive (k,ψ)-Hilfer fractional double integro-differential equation with mixed nonlocal conditions. AIMS Mathematics, 2023, 8(9): 20437-20476. doi: 10.3934/math.20231042
    [9] Xinyi Zhang, Jian Zhang . On Schrödinger-Poisson equations with a critical nonlocal term. AIMS Mathematics, 2024, 9(5): 11122-11138. doi: 10.3934/math.2024545
    [10] Zeliha Korpinar, Mustafa Inc, Ali S. Alshomrani, Dumitru Baleanu . The deterministic and stochastic solutions of the Schrodinger equation with time conformable derivative in birefrigent fibers. AIMS Mathematics, 2020, 5(3): 2326-2345. doi: 10.3934/math.2020154
  • A nonlocal derivative nonlinear Schrödinger (DNLS) equation is analytically studied in this paper. By constructing Darboux transformations (DTs) of arbitrary order, new determinant solutions of the nonlocal DNLS equation in the form of Wronskian-type are derived from both zero and nonzero seed solutions. Periodic solitons are obtained with different parameter choices. When one eigenvalue tends to another one, generalized DTs are constructed, leading to rogue waves. Due to complex parametric constraints, the derived solutions may have singularities. Despite this, the work presented in this paper can still provide a valuable reference for the study of nonlocal integrable systems.



    Integrable evolution equations are very common and play a key role in nonlinear science, as shown in studies like [1,2,3]. Specifically, the Korteweg-de Vries (KdV) equation, which comes from fluid dynamics, explains how waves with slight dispersion and small amplitude evolve in media with quadratic nonlinearity [2,4]. Additionally, the nonlinear Schrödinger (NLS) equation describes the behavior of waves that are slightly nonlinear and quasi-monochromatic in media with cubic nonlinearity [3,5,6]. Moreover, the derivative nonlinear Schrödinger (DNLS) equation is used to describe Alfvén waves in plasmas with low-β values, as well as large-amplitude waves that run parallel to the magnetic field [7,8,9]. These equations are mathematically integrable, meaning they can be solved exactly.

    Recently, Ablowitz and Musslimani introduced the following parity time (PT)-symmetric nonlocal NLS equation and derived its explicit solutions[10]:

    iut(x,t)=uxx(x,t)+2εu(x,t)2u(x,t)(ε=±1), (1.1)

    where ε=±1 denotes the focusing (+) and defocusing () nonlinearity, and the asterisk denotes the complex conjugate. The PT symmetry implies that the self-induced potential V(x,t)=u(x,t)u(x,t) fulfills the relation V(x,t)=V(x,t), while the nonlocality indicates that the value of the potential V(x,t) at x requires the information on u(x,t) at x, as well as at x[11,12,13]. Since Eq (1.1) exhibits properties akin to those of the standard NLS equation, such as integrability in the Lax and inverse scattering senses, a significant amount of theoretical research has been conducted[14,15,16], and the study on the equations has extended to higher-order and coupled ones[17,18,19]. Furthermore, there has been experimental evidence observed in the fields of optics and materials science[20,21,22].

    Subsequently, some other nonlocal integrable equations have been proposed[23,24,25]. One of these nonlocal equations is the nonlocal DNLS equation[8,9],

    iut(x,t)=uxx(x,t)+ε[u(x,t)2u(x,t)]x(ε=±1), (1.2)

    which is obtained from the standard DNLS equation by replacing |u|2u with u(x,t)2u(x,t), and the conserved density V(x,t)=u(x,t)u(x,t) satisfies the relation V(x,t)=V(x,t). In [8], hierarchies of nonlocal DNLS-type equations have been derived utilizing Lie algebra splitting, with the specific nonlocal DNLS equation being termed as nonlocal DNLSI. Zhou[9] formulates the 2n-fold Darboux transformation (DT) and derives global bounded solutions from zero seed solutions for Eq (1.2). In this paper, we concentrate on how to construct arbitrary order DT of Eq (1.2), what happens if the eigenvalues degenerate, and how to construct solutions from different seeds.

    This paper is organized as follows: In Section 2, we give a detailed derivation of the Wronskian-type determinant representation of arbitrary order DT. Moreover, different choice of the eigenvalues may lead different Darboux matrix, even in the same order. If there are duplicate eigenvalues, a generalized form of the DT may prove useful. In Section 3, we obtain solutions of Eq (1.2) from zero and non-zero seed solutions. Periodic solutions and rogue waves are constructed despite of the existence of singularity. Finally, conclusions and discussions are given in Section 4.

    In this section, we will construct the DT of the coupled nonlocal DNLS equations

    iutuxx+(u2v)x=0, (2.1a)
    ivt+vxx+(uv2)x=0, (2.1b)

    which is the compatibility condition of the following Lax pair [9]:

    Φx=UΦ,Φt=VΦ, (2.2)

    where Φ=(f(x,t),g(x,t))T,

    U=λ2σ3+λP=(λ2λuλvλ2), (2.3)
    V=i(λ2uv2λ4)σ32iλ3P+iλ(P3σ3Px)=(iλ2uv2iλ4iλ(u2vux)2iλ3uiλ(uv2+vx)2iλ3viλ2uv+2iλ4), (2.4)

    with

    P=(0uv0), (2.5)

    and Pauli matrix

    σ3=(1001). (2.6)

    Under the reduction condition v(x,t)=εˉu(x,t)=εu(x,t), Eq (2.1) will be reduced to the nonlocal DNLS Eq (1.2). Here, ˉf(x,t) means f(x,t). Unlike the local DNLS equation, there are four linear independence solutions of (2.2) associate with one complex eigenvalue. Specifically, if (fg) are solutions of (2.2) with λ, then (fg), (εˉgˉf), and (εˉgˉf) are solutions of (2.2) with λ, λ, and λ, respectively.

    We can suppose Φ[1]=TΦ as the DT, where

    T=(λa1(x,t)+a0(x,t)λb1(x,t)+b0(x,t)λc1(x,t)+c0(x,t)λd1(x,t)+d0(x,t)) (2.7)

    is the first-order DT matrix, and there exist U[1] and V[1] possessing the same form as U and V with u and v replaced by the new potentials u[1] and v[1], which satisfy Eq (2.2). Therefore, T[λ] must satisfy conditions

    Tx+TUU[1]T=0, (2.8a)
    Tt+TVV[1]T=0. (2.8b)

    We set Φ1=(f1g1) are solutions of (2.2) with λ=λ1. Comparing the coefficients of λj,(j=3,2,1,0) in Eq (2.8), we obtain b1(x,t)=c1(x,t)=a0(x,t)=d0(x,t)=0, b0(x,t)=b0,c0(x,t)=c0 and a1(x,t)=g1f1,d1(x,t)=f1g1 by setting b0=c0=λ1, so we have

    u[1]=ρ12u+2λ1ρ1, (2.9a)
    v[1]=1ρ12v2λ11ρ1, (2.9b)

    where ρ1=g1f1.

    Theorem 1. Suppose λ10, which is real when ε=1 or purely imaginary when ε=1. Define Φ[1] and ρ1 as above such that ρ1¯ρ1=1. Then Eq (2.9) are new solutions of (1.2), and the corresponding Darboux matrix is

    T=(λρ1λ1λ1λρ1)=λ(ρ11ρ1)(λ1λ1). (2.10)

    Proof. We only need to prove v[1]=εˉu[1], which can be easily verified by using (2.9) and ρ1¯ρ10.

    As λ1 will be real or purely imaginary, one-fold DT only allows two eigenvalues λ1 and λ1, so one-fold DT is degenerated in a sense. If we choose λ1, λ1, λ1, and λ1 as four different eigenvalues, the real and imaginary part of λ1 are both not equal to zero and we need to construct multi-fold DT.

    In this subsection, we will construct the two-fold DT of (1.2). The multi-fold DT can be iterated by the one-fold DT, so we set second-order Darboux matrix as

    T2=λ2(a21a2)+λ(b1c1)+(a0a0). (2.11)

    We set Φk=(fkgk) are solutions of (2.2) with λ=λk, ρk=gkfk(k=1,2), and then T2 maps (u,v,Φ) to (u[2],v[2],Φ[2]). By following the previous steps, we obtain

    u[2]=a22u2a2b1,v[2]=va22+2c1a2 (2.12)

    with

    a0=λ1λ2,a2=|λ1g1f1λ2g2f2||λ1f1g1λ2f2g2|,b1=|f1λ12f1f2λ22f2||λ1f1g1λ2f2g2|,c1=|g1λ12g1g2λ22g2||λ1g1f1λ2g2f2|. (2.13)

    If λ2=λ1 and λ1±λ1, we can obtain Theorem 1 in [9]. Moreover, if λ1,λ2 are both real or purely imaginary, we can also obtain new solutions of (1.2).

    From above results, we know the DT will be iterated two times if λ is a complex eigenvalue and one times if λ is real or purely imaginary. Therefore, if we choose λk(k=0,1,2,,n1) as complex eigenvalues and λj(j=0,1,,n2) as real (ε=1) or purely imaginary (ε=1) eigenvalues, then we can construct the (2n1+n2)-fold DT of (1.2). For example, we can construct third-order DT by choosing (n1,n2)=(1,1) or (0,3) and fourth-order DT by choosing (n1,n2)=(2,0), (1,2), or (0,4). Even though the 2n-fold DT have been derived in [9] by choosing n1=n,n2=0, we construct the Wronskian-type determinant representation of the Darboux matrix.

    2n-order DT

    We set 2n-order Darboux matrix as

    T2n(λ)=λ2n(a2nd2n)+λ2n1(b2n1c2n1)++(a0d0), (2.14)

    and Φk=(fkgk) are solutions of (2.2) with λ=λk, ρk=gkfk(k=1,2,,2n), and then T2n maps (u,v,Φ) to (u[2n],v[2n],Φ[2n]), where a2k,d2k(k=0,1,,n), b2k1,c2k1(k=1,2,,n) are 4n+2 unknown functions. Solving 4n equations T2n(λk)Φ[2n]=0(k=1,2,,2n) by setting d2n=1a2n and d0=a0, we can uniquely determine the 4n+2 unknown functions, and obtain

    u[2n]=a2n2u2a2nb2n1,v[2n]=va2n2+2c2n1a2n (2.15)

    with

    a2n=(1)na0D3D1,d2n=(1)na0D4D2,b2n1=a0D5D1,c2n1=a0D6D2, (2.16)

    and

    D1=|F(2,4,,2n)(1,2,,2n)G(1,3,,2n1)(1,2,,2n)|, (2.17a)
    D2=|F(1,3,,2n1)(1,2,,2n)G(2,4,,2n)(1,2,,2n)|, (2.17b)
    D3=|F(0,2,,2n2)(1,2,,2n)G(1,3,,2n1)(1,2,,2n)|, (2.17c)
    D4=|F(1,3,,2n1)(1,2,,2n)G(0,2,,2n2)(1,2,,2n)|, (2.17d)
    D5=|F(0,2,,2n)(1,2,,2n)G(1,3,,2n3)(1,2,,2n)|, (2.17e)
    D6=|F(1,3,,2n3)(1,2,,2n)G(0,2,,2n)(1,2,,2n)|, (2.17f)

    where

    F(2,4,,2n)(1,2,,2n)=[f1λ12f1λ14f1λ12nf2λ22f2λ24f2λ22nf2nλ2n2f2nλ2n4f2nλ2n2n], (2.18)
    G(1,3,,2n1)(1,2,,2n)=[g1λ1g1λ13g1λ12n1g2λ2g2λ23g2λ22n1g2nλ2ng2nλ2n3g2nλ2n2n1], (2.19)

    and other matrix blocks with the following blocks have been similarly defined.

    As D1=2nj=1λjD4, D2=2nj=1λjD3, and a0=2nj=1λj, then d2n=1a2n will be satisfied automatically.

    If we choose λn+j=λj and Φn+j=(ε¯gj¯fj)(j=1,2,,n), then a0=nj=1|λj|2, D1=|F(2,4,,2n)(1,2,,n)G(1,3,,2n1)(1,2,,n)εˉG(2,4,,2n)(1,2,,n)ˉF(1,3,,2n1)(1,2,,n)| and so are D2,,D6, where ˉG(2,4,,2n)(1,2,,n) means doing complex conjugate and space-reverse operation to G(2,4,,2n)(1,2,,n).

    Theorem 2. Suppose λj(j=1,2,,n) are complex, which are all not real or purely imaginary. Then Eq (2.15) are new solutions of (1.2). Especially, when n=1, (2.15) become (2.12).

    Proof. We only need to prove v[2n]=εˉu[2n], which can be easily verified by using ¯D3=D4, and ¯D5=εD6.

    Remark 1. It is a little different if we choose λk+j=λj, Φk+j=(ε¯gj¯fj)(j=1,2,,k) and λs(s=2k+1,,2n) are all real (e.g.,ε=1), here a0=kj=1|λj|22ns=2k+1λs, D1=|F(2,4,,2n)(1,2,,k)G(1,3,,2n1)(1,2,,k)εˉG(2,4,,2n)(1,2,,k)ˉF(1,3,,2n1)(1,2,,k)F(2,4,,2n)(2k+1,,2n)G(1,3,,2n1)(2k+1,,2n)| and so are D2,,D6. Equation (2.15) is new solutions of (1.2) with constraints ρs¯ρs=1(s=2k+1,,2n).

    (2n1)-order DT

    We set (2n1)-order Darboux matrix as

    T2n1(λ)=(nk=1a2k1λ2k1n1k=0b2kλ2kn1k=0c2kλ2knk=1d2k1λ2k1), (2.20)

    and Φk=(fkgk) are solutions of (2.2) with λ=λk, ρk=gkfk(k=1,2,,2n1). Then T2n1 maps (u,v,Φ) to (u[2n1],v[2n1],Φ[2n1]), where a2k1,d2k1(k=1,2,,n), b2k,c2k(k=0,1,,n1) are 4n unknown functions. Solving 4n2 equations T2n1(λk)Φ[2n1]=0(k=1,2,,2n1) by setting d2n1=1a2n1 and c0=b0, we can uniquely determine the 4n unknown functions, and obtain

    u[2n1]=a2n12u2a2n1b2n2,v[2n1]=va2n12+2c2n2a2n1 (2.21)

    with

    a2n1=b0D3D1,d2n1=b0D4D2,b2n2=(1)n1b0D5D1,c2n2=(1)n1b0D6D2, (2.22)

    where

    D1=|F(1,3,,2n1)(1,2,,2n1)G(2,4,2n2)(1,2,,2n1)|, (2.23a)
    D2=|F(2,4,,2n2)(1,2,,2n1)G(1,3,,2n1)(1,2,,2n1)|, (2.23b)
    D3=|F(1,3,,2n3)(1,2,,2n1)G(0,2,,2n2)(1,2,,2n1)|, (2.23c)
    D4=|F(0,2,,2n2)(1,2,,2n1)G(1,3,,2n3)(1,2,,2n1)|, (2.23d)
    D5=|F(1,3,,2n1)(1,2,,2n1)G(0,2,,2n4)(1,2,,2n1)|, (2.23e)
    D6=|F(0,2,,2n4)(1,2,,2n1)G(1,3,,2n1)(1,2,,2n1)|. (2.23f)

    Noticing D1=2n1j=1λjD4, D2=2n1j=1λjD3 and b0=2n1j=1λj, then d2n=1a2n will be satisfied automatically.

    If we choose λk+j=λj and Φk+j=(ε¯gj¯fj)(j=1,2,,k), then b0=kj=1|λj|22n1s=2k+1λs, D1=|F(1,3,,2n1)(1,2,,k)G(2,4,,2n2)(1,2,,k)εˉG(1,3,,2n1)(1,2,,k)ˉF(2,4,,2n2)(1,2,,k)F(1,3,,2n1)(2k+1,2k+2,,2n1)G(2,4,,2n2)(2k+1,2k+2,,2n1)| and so are D2,,D6.

    Theorem 3. Suppose λj(j=1,2,,k) are complex, and λs(s=2k+1,2k+2,,2n1) are real when ε=1 or purely imaginary when ε=1, then Eq (2.21) is new solutions of (1.2) with constraints ρs¯ρs=1(s=2k+1,,2n1). Especially, when n=1, (2.21) become (2.9).

    Proof. Similarly, we only need to prove v[2n1]=εˉu[2n1], which can be easily verified by using ¯D3=(1)(n21)εnD4, and ¯D5=(1)(n21)εn1D6.

    Until now, we have constructed the arbitrary DT of (1.2) with Wronskian-type determinant representation. Compared with [9], we have not only constructed the DT of arbitrary order, but also adopted a simpler and more feasible way.

    In this section, we will construct several explicit solutions of Eq (1.2).

    • If we set seed solution u=v=0, then solution of (2.2) with λ=λk is Φk=(fkgk)=(αkeθkβkeθk), where θk=λ2kx2iλ4kt(k=1,2,). As λ1 is real when ε=1 or purely imaginary when ε=1 in Theorem 1, background of the new solution of (1.2) will go to infinity, which appears meaningless.

    • If we set seed solution u=ρeερ2x+iϕ, the background will also go to infinity. So, we take u=ρeiϕ as the seed solution eclectically to avoid infinity of the background. Solution of (2.2) with λ=λk is Φk=(fkgk), where

    fk=αkeθk+βkeθk,gk=λkeiϕρfk+eiϕλkρfkx (3.1)

    with θk=λ4kλ2kρ2ε[xit(2λ2k+ρ2ε)], ρk=gkfk, αk and βk(k=1,2,) are none-zero complex constants.

    Using Theorem 1, we obtain

    u[1]=eiϕ[β21ρ2ε+2α1β1(2λ21ρ2ε)e2θ1+α21ρ2εe4θ1]ρ[β1+α1e2θ1]2 (3.2)

    under the constraint

    |α1|2|β1|2=λ21ρ2ε+λ21(ρ2ελ21)λ21+ρ2ε+λ21(ρ2ελ21), (3.3)

    which means |λ1|>|ρ|. This solution has been shown in Figure 1, which displays the periodic one-soliton with ε=1,ρ=4,λ1=5,α1=2,β1=1, ϕ=π, and ε=1,ρ=4,λ1=5i,α1=1/2,β1=1, ϕ=0. However, there are periodic singularities in any case for ε. The reason for the existence of the singularity may be that θk in solution (3.1) is too complex. When ε takes the values of -1 and 1, respectively, eigenvalue λ will differ due to the constraints, which subsequently causes variations in other parameters. Additionally, we found that the positions of the singular points have changed.

    Figure 1.  The periodic one-soliton from (3.2) with (a) ε=1,ρ=4,λ1=5,α1=2,β1=1, ϕ=π, and (b) ε=1,ρ=4,λ1=5i,α1=1/2,β1=1, ϕ=0.

    When the eigenvalue λ1=±ρ(ε=1) or λ1=±iρ(ε=1), the new solution will be complex constant and trivial, so we should resolve (2.2). Now, the solution of (2.2) (ε=1) is Φ1=(f1g1), where

    f1=α1+α2(x3iρ2t),g1=eiϕf1eiϕρ2f1x (3.4)

    with αk(k=1,2) are none-zero complex constants. We can obtain a new fundamental rational solution,

    u[1]=eiϕ[α21ρ4+α22(9ρ8t26iρ6tx+ρ4x21)+2α1α2ρ4(x3iρ2t)]ρ3[α1+α2(x3iρ2t)]2. (3.5)

    The graph of this solution is displayed in Figure 2(a) with ρ=1,α1=1,α2=2 and ϕ=π. It is seen that this is a rogue wave, similar to Eq (40) in Yang's article[26], which blows up to infinity at x=1/2 and finite time t=0. Furthermore, rising from a constant background, rogue wave (3.5) can be considered as resulting from the degradation of periodic soliton (3.2). If we set ε=1, then new rogue wave will blow up to infinity at x=1/2 and t=0, which is shown in Figure 2(b). Similarly, variations in the value of parameter ε only affects the position of solitons or rogue waves, without altering the type or mechanical properties of the solutions. Meanwhile, parameter ϕ do not affect the modulus of the solutions.

    Figure 2.  The one rogue wave from (3.5) with (a) ε=1,ρ=1,λ1=1,α1=1,α2=2, ϕ=π, and (b) ε=1,ρ=1,λ1=i,α1=2,α2=4, ϕ=0.

    • When we set seed solution u=v=0, the background of new solutions will not go to infinity if and only if λ2 is pure imaginary, which means 4argλ/π is an odd integer [9]. Hence, we can obtain a new soliton solution with λ1=a(1i),

    u[2]=(4+4i)εaα1β1e4ia2(4a2t+x)(α1α1+iεβ1β1e8ia2x)(εβ1β1e8ia2x+iα1α1)2. (3.6)

    This solution does not have singularity and global if |α1||β1|, which is periodic in both x and t. The graph is displayed in Figure 3 with ε=1,a=1,α1=1 and β1=2, from which we can see that the norm of the solution is periodic in x but independent of t. However, if λ1 and λ2 are both real or purely imaginary, the background of new solution will go to infinity.

    Figure 3.  Left panel is the periodic two solitons from (3.6) with ε=1,λ1=1i,α1=1 and β1=2; Right panel is the corresponding density plot.

    • If we set seed solution u=ρeερ2x+iϕ, the background will go to infinity. So, we take u=ρeiϕ as the seed solution, and the solution of (2.2) with λ=λk is (3.1). Using Eq (2.12), we can obtain new solutions of Eq (1.2) and display it in Figure 4. Figure 4 show the propagation of two solitons, the interaction of which seems elastic. However, the solution blows up to infinity in the interaction region at one fixed point. The reason for the existence of singularity may be due to nonlocal effects. Figure 5(a) shows the interaction of two periodic soliton with periodic singularities, as long as λ1,λ2 are both real numbers. This solution seems to be a superposition and interaction of two solutions (3.2). Furthermore, if λ1=±ρ, then we can see that the interaction of periodic soliton and rogue wave, which has been shown in Figure 5(b), will degenerate into the seed solution if λ2=±ρ.

    Figure 4.  Left panel is the interaction of two solitons from (2.12) with ε=1,ρ=1,λ1=(1+i3)/2,α1=1,β1=i,ϕ=0; Right panel is the corresponding density plot.
    Figure 5.  (a) Inteaction of two periodic solitons from (2.12) with ε=1,ρ=1,λ1=2, λ2=3,α1=23+3233,α2=62+8628,β1=β2=1,ϕ=0; (b) Interaction of one periodic soliton and rogue wave from (2.12) with ε=1,ρ=1,λ1=1,λ2=2,α1=1+i,α2=23+3233,β1=2,β2=1,ϕ=0.

    • Considering the third-order DT, three eigenvalues can be one complex and one real (or pure imaginary), or three real. New solutions derived from seeds u=ρeiϕ would be the superposition of previous solutions (2.9) and (2.12). For example, if we choose three eigenvalues as one complex and one real, new solutions is displayed in Figure 6(a). We can see Figure 6(a) is a superposition and interaction of Figure 1(a) and Figure 4. There is one periodic soliton and two elastic-interacted solitons. If the real eigenvalue equals ±ρ, then the solution will be different from Figure 6(a), which is shown in Figure 6(b). The phases of the solitons undergo changes, and two of them remain parallel both before and after the collision.

    Figure 6.  Interaction of three solitons from (2.21) with (a) ε=1,ρ=1,λ1=(1+i3)/2,λ2=2,α1=1,α2=23+3233,β1=β2=1,ϕ=0; (b) ε=1,ρ=λ1=1,λ2=(1+i3)/2,α1=α2=β1=β2=1,ϕ=0; (c) ε=1,ρ=1,λ1=2,λ2=3, λ3=4, α1=23+3233,α2=62+8628,α3=415+1541515,β1=β2=β3=1,ϕ=0; (d) ε=1,ρ=1, λ1=1,λ2=2,λ3=2,α1=1,α2=23+3233,α3=62+8628,β1=2,β2=β3=1,ϕ=0.

    If three eigenvalues are all selected real, then the new solution is displayed in Figure 6(c), which are three times superposition and interaction of Figure 1(a). Further more, if one of three real eigenvalues equals ±ρ, then there will be a rogue wave, which has been shown in Figure 6(d).

    • Considering the fourth-order DT, we choose zero seed solutions and two complex eigenvalues λ1, λ2, where λ2k(k=1,2) is pure imaginary to avoid infinity of the background. If λ2λ1, we can obtain periodic solitons, which have been shown in Figure 7. If λ2=λ1, then solution of (2.2) with λ=λ2, says Φ2, will linearly dependent on Φ1. Following the idea of generalized DT in [7], we set Φ[1]1=Φ1(λ1+δ),δ1, and obtain

    Φ2=limδ0Φ1(λ1+δ)Φ1(λ1)δ, (3.7)

    which is another solution of (2.2) with λ=λ1 but linearly independent on Φ1. Using (2.15), we obtain new solutions of (1.2),

    (4+4i)(64te8ix64itie8ix8e8ixx8ix1)(4096it2e8ix+128te8ix+64ie8ixx2ie8ix+e16ix1)e16it4ix(4096it2e8ix+128te8ix64ie8ixx2+ie8ix+e16ix1)2 (3.8)

    with ε=1,α1=β1=1. This solution is a mixed-type of exponential and rational functions, whose norm is displayed in Figure 8. From it we can see there are three solitons, one x-periodic soliton and two elastic-interacted solitons, and no rogue waves are presented obviously. However, the interaction blows up to infinity at two fixed points. As it has been seen that higher order solution can be the superposition and interaction of the lower order solution, we do not show the solutions derived from four real eigenvalues or one complex and two real eigenvalues. At the end, we display new periodic solitons in Figure 9 obtained by using six-order DT from zero seeds.

    Figure 7.  Left panel is the periodic four solitons from (2.15) with ε=1,λ1=1+i, λ2=2(1i), α1=α2=1/5 and β1=β2=1; Right panel is the corresponding density plot.
    Figure 8.  Left panel is interaction of three solitons from (3.8) with ε=1,λ1=λ2=1+i,α1=β1=1; Right panel is the corresponding density plot.
    Figure 9.  Left panel is interaction of periodic multi-solitons from (2.15) (n = 3) with ε=1, λ1=1+i,λ2=2(1i),λ3=3(1+i),α1=0.1,α2=0.2,α3=0.3,β1=β2=β3=1; Right panel is the corresponding density plot.

    In this paper, the nonlocal DNLS Eq (1.2) has been analytically studied. Firstly, using the Lax pair of the coupled Eq (2.1), from which the reduction conditions lead to Eq (1.2), we have constructed arbitrary order DTs in a more streamlined manner. From zero and non-zero seed solutions, new Wronskian-type determinant solutions are obtained. Specifically, we have obtained new solutions by applying the first four orders DT.

    Starting from zero seeds, due to the nonlocality of the equation, odd-order DTs will cause the background of the solutions to diverge towards infinity. For even-order DTs, the background of the solutions will not go to infinity unless the square of eigenvalues λ2 are pure imaginary, in order to maintain the reduction conditions, which could increase the difficulty in selecting parameter values. Soliton solutions without singularities have been obtained for specific parameter values, as demonstrated in Figures 3, 7, and 9. These new solutions are all periodic and bounded, similar to the results reported in [9]. Besides, when one eigenvalue tends to another one, generalized DT can be constructed. Taking the solitons in Figure 7 as an example, we can obtain new solutions which are mixed-type of exponential and rational functions with λ2=λ1, which can be seen Figure 8. Although the solution in Figure 7 is bounded with λ2λ1, there are two singularities, causing the solitons to blow up to infinity.

    Starting from non-zero seeds, we have obtained various solutions including the periodic solutions and rogue waves. If the eigenvalues are real or pure imaginary numbers, there will be constraints on the parameters. However, these solutions inherently possess singularities. Nevertheless, we have constructed soliton solutions for the first few orders. Figure 1 displays the periodic one-soliton, which degrades into a rogue wave (see Figure 2) when the magnitude of the eigenvalue is identical to the amplitude of the seed solution. Similarly, periodic two- and three-soliton solutions can also degrade or partially degrade into rogue waves, which can be observed in Figures 5, 6(c), and 6(d). When the eigenvalues are taken as complex numbers, the interaction of two solitons has been shown in Figure 4. Figure 6(a) shows a superposition and interaction of Figures 1(a) and 4. When the real eigenvalue equals to ±ρ, the periodic soliton will disappear, as shown in Figure 6(b).

    In addition, the value of ε does not affect much except the location of the singularities, which is different from the nonlocal NLS Eq (1.1). We hope that the work presented in this paper can provide a valuable reference and assistance for the study of nonlocal integrable systems.

    The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

    Thanks to the valuable suggestions of the editor and reviewers. Work of this paper was supported by State Key Laboratory of Heavy Oil Processing (SKLHOP2024115808).

    The author declares no conflict of interest.



    [1] M. Ablowitz, H. Segur, Solitons and inverse scattering transform, Philadelphia: SIAM, 1981. https://doi.org/10.1137/1.9781611970883
    [2] Z. Li, Z. Zhao, N-periodic wave solutions of the N = 2 supersymmetric KdV equation, Appl. Math. Lett., 160 (2025), 109313. https://doi.org/10.1016/j.aml.2024.109313 doi: 10.1016/j.aml.2024.109313
    [3] P. Xin, Z. Zhao, Y. Wang, Quasi-periodic breathers and their dynamics to the Fokas system in nonlinear optics, Wave Motion, 133 (2025), 103449. https://doi.org/10.1016/j.wavemoti.2024.103449 doi: 10.1016/j.wavemoti.2024.103449
    [4] Y. Liu, Z. Zhao, Periodic line wave, rogue waves and the interaction solutions of the (2+1)-dimensional integrable Kadomtsev-Petviashvili-based system, Chaos Soliton. Fract., 183 (2024), 114883. https://doi.org/10.1016/j.chaos.2024.114883 doi: 10.1016/j.chaos.2024.114883
    [5] M. Li, T. Xu, Dark and antidark soliton interactions in the nonlocal nonlinear Schrödinger equation with the self-induced parity-time-symmetric potential, Phys. Rev. E, 91 (2015), 033202. https://doi.org/10.1103/PhysRevE.91.033202 doi: 10.1103/PhysRevE.91.033202
    [6] M. Li, T. Xu, D. Meng, Rational solitons in the parity-time-symmetric nonlocal nonlinear Schrödinger model, J. Phys. Soc. Jpn., 85 (2016), 124001. https://doi.org/10.7566/JPSJ.85.124001 doi: 10.7566/JPSJ.85.124001
    [7] B. Guo, L. Ling, Q. Liu, High-order solutions and generalized Darboux transformations of derivative nonlinear Schrödinger equations, Stud. Appl. Math., 130 (2013), 317–344. https://doi.org/10.1111/j.1467-9590.2012.00568.x doi: 10.1111/j.1467-9590.2012.00568.x
    [8] Z. Wu, J. He, New hierarchies of derivative nonlinear Schrödinger-type equation, Rom. Rep. Phys., 68 (2016), 1425–1446. https://doi.org/10.48550/arXiv.1704.02070 doi: 10.48550/arXiv.1704.02070
    [9] Z. Zhou, Darboux transformations and global solutions for a nonlocal derivative nonlinear Schrödinger equation, Commun. Nonlinear Sci., 62 (2018), 480–488. https://doi.org/10.1016/j.cnsns.2018.01.008 doi: 10.1016/j.cnsns.2018.01.008
    [10] M. Ablowitz, Z. Musslimani, Integrable nonlocal nonlinear Schrödinger equation, Phys. Rev. Lett., 110 (2013), 064105. https://doi.org/10.1103/PhysRevLett.110.064105 doi: 10.1103/PhysRevLett.110.064105
    [11] C. Bender, S. Boettcher, Real spectra in non-Hermitian Hamiltonians having PT symmetry, Phys. Rev. Lett., 80 (1998), 5243. https://doi.org/10.1103/PhysRevLett.80.5243 doi: 10.1103/PhysRevLett.80.5243
    [12] T. Xu, H. Li, H. Zhang, M. Li, S. Lan, Darboux transformation and analytic solutions of the discrete PT-symmetric nonlocal nonlinear Schrödinger equation, Appl. Math. Lett., 63 (2017), 88–94. https://doi.org/10.1016/j.aml.2016.07.024 doi: 10.1016/j.aml.2016.07.024
    [13] D. Meng, K. Li, Darboux transformation of the second-type nonlocal derivative nonlinear Schrödinger equation, Mod. Phys. Lett. B, 33 (2019), 1950123. https://doi.org/10.1142/S0217984919501239 doi: 10.1142/S0217984919501239
    [14] A. Fokas, Integrable multidimensional versions of the nonlocal nonlinear Schrödinger equation, Nonlinearity, 29 (2016), 319. https://doi.org/10.1088/0951-7715/29/2/319 doi: 10.1088/0951-7715/29/2/319
    [15] X. Huang, L. Ling, Soliton solutions for the nonlocal nonlinear Schrödinger equation, Eur. Phys. J. Plus, 131 (2016), 148. https://doi.org/10.1140/epjp/i2016-16148-9 doi: 10.1140/epjp/i2016-16148-9
    [16] M. Ablowitz, Z. Musslimani, Inverse scattering transform for the integrable nonlocal nonlinear Schrödinger equation, Nonlinearity, 29 (2016), 915. https://doi.org/10.1088/0951-7715/29/3/915 doi: 10.1088/0951-7715/29/3/915
    [17] A. Khare, A. Saxena, New solutions of coupled nonlocal NLS and coupled nonlocal mKdV equations, Ann. Phys., 463 (2024), 169627. https://doi.org/10.1016/j.aop.2024.169627 doi: 10.1016/j.aop.2024.169627
    [18] H. Chen, H. Zhu, Higher-dimensional vector two-component solitons of a nonautonomous partially nonlocal coupled NLS model in a linear and harmonic potential, Nonlinear Dyn., 111 (2023), 581–590. https://doi.org/10.1007/s11071-022-07629-7 doi: 10.1007/s11071-022-07629-7
    [19] J. Chen, Y. Lu, B. Wang, Global Cauchy problems for the nonlocal (derivative) NLS in Esσ, J. Differ. Equations, 344 (2023), 767–806. https://doi.org/10.1016/j.jde.2022.11.002 doi: 10.1016/j.jde.2022.11.002
    [20] A. Regensburger, M. Miri, C. Bersch, J. Näger, G. Onishchukov, D. Christodoulides, et al., Observation of defect states in PT-symmetric optical lattices, Phys. Rev. Lett., 110 (2013), 223902. https://doi.org/10.1103/PhysRevLett.110.223902 doi: 10.1103/PhysRevLett.110.223902
    [21] B. Peng, Ş. Özdemir, F. Lei, F. Monifi, M. Gianfreda, G. Long, et al., Nonreciprocal light transmission in parity-time-symmetric whispering-gallery microcavities, Nature Phys., 10 (2014), 394–398. https://doi.org/10.1038/nphys2927 doi: 10.1038/nphys2927
    [22] L. Feng, Y. Xu, W. Fegadolli, M. Lu, J. Oliveira, V. Almeida, et al., Experimental demonstration of a unidirectional reflectionless parity-time metamaterial at optical frequencies, Nature Mater., 12 (2013), 108–113. https://doi.org/10.1038/nmat3495 doi: 10.1038/nmat3495
    [23] K. Chen, X. Deng, S. Lou, D. Zhang, Solutions of nonlocal equations reduced from the AKNS hierarchy, Stud. Appl. Math., 141 (2018), 113–141. https://doi.org/10.1111/sapm.12215 doi: 10.1111/sapm.12215
    [24] W. Liu, D. Qiu, Z. Wu, J. He, Dynamical behavior of solution in integrable nonlocal Lakshmanan-Porsezian-Daniel equation, Commun. Theor. Phys., 65 (2016), 671. https://doi.org/10.1088/0253-6102/65/6/671 doi: 10.1088/0253-6102/65/6/671
    [25] J. Rao, Y. Cheng, J. He, Rational and semirational solutions of the nonlocal Davey-Stewartson equations, Stud. Appl. Math., 139 (2017), 568–598. https://doi.org/10.1111/sapm.12178 doi: 10.1111/sapm.12178
    [26] B. Yang, J. Yang, Transformations between nonlocal and local integrable equations, Stud. Appl. Math., 140 (2018), 178–201. https://doi.org/10.1111/sapm.12195 doi: 10.1111/sapm.12195
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(703) PDF downloads(64) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog