In this paper, we investigate non-traveling wave solutions of the (3+1)-dimensional variable coefficients Date-Jimbo-Kashiwara-Miwa (VC-DJKM) equation, which describes the real physical phenomena owing to the inhomogeneities of media. By combining the extended homoclinic test approach with variable separation method, we obtain abundant new exact non-traveling wave solutions of the (3+1)-dimensional VC-DJKM equation. These results with a parabolic tail or linear tail reveal the complex structure of the solutions for (3+1)-dimensional VC-DJKM equation. Moreover, the tail in these solutions maybe give a prediction of physical phenomenon. When arbitrary functions contained in these non-traveling wave solutions are taken as some special functions, we can get the kink-type solitons, singular solitary wave solutions, and periodic solitary wave solutions, and so on. As the special cases of our work, the corresponding results of (3+1)-dimensional DJKM equation, (2+1)-dimensional DJKM equation, (2+1)-dimensional VC-DJKM equation are also given.
Citation: Yuanqing Xu, Xiaoxiao Zheng, Jie Xin. New non-traveling wave solutions for (3+1)-dimensional variable coefficients Date-Jimbo-Kashiwara-Miwa equation[J]. AIMS Mathematics, 2021, 6(3): 2996-3008. doi: 10.3934/math.2021182
[1] | Ricardo Almeida . Variational problems of variable fractional order involving arbitrary kernels. AIMS Mathematics, 2022, 7(10): 18690-18707. doi: 10.3934/math.20221028 |
[2] | Xiaojing Du, Xiaotong Liang, Yonghong Xie . Integral expressions of solutions to higher order λ-weighted Dirac equations valued in the parameter dependent Clifford algebra. AIMS Mathematics, 2025, 10(1): 1043-1060. doi: 10.3934/math.2025050 |
[3] | Yongjian Hu, Huifeng Hao, Xuzhou Zhan . On the solvability of the indefinite Hamburger moment problem. AIMS Mathematics, 2023, 8(12): 30023-30037. doi: 10.3934/math.20231535 |
[4] | Hong Li, Keyu Zhang, Hongyan Xu . Solutions for systems of complex Fermat type partial differential-difference equations with two complex variables. AIMS Mathematics, 2021, 6(11): 11796-11814. doi: 10.3934/math.2021685 |
[5] | Kun Li, Peng Wang . Properties for fourth order discontinuous differential operators with eigenparameter dependent boundary conditions. AIMS Mathematics, 2022, 7(6): 11487-11508. doi: 10.3934/math.2022640 |
[6] | Tuba Gulsen, Emrah Yilmaz, Ayse Çiğdem Yar . Proportional fractional Dirac dynamic system. AIMS Mathematics, 2024, 9(4): 9951-9968. doi: 10.3934/math.2024487 |
[7] | Valérie Gauthier-Umaña, Henryk Gzyl, Enrique ter Horst . Decoding as a linear ill-posed problem: The entropy minimization approach. AIMS Mathematics, 2025, 10(2): 4139-4152. doi: 10.3934/math.2025192 |
[8] | Yong Liu, Chaofeng Gao, Shuai Jiang . On meromorphic solutions of certain differential-difference equations. AIMS Mathematics, 2021, 6(9): 10343-10354. doi: 10.3934/math.2021599 |
[9] | Clara Burgos, Juan Carlos Cortés, Elena López-Navarro, Rafael Jacinto Villanueva . Probabilistic analysis of linear-quadratic logistic-type models with hybrid uncertainties via probability density functions. AIMS Mathematics, 2021, 6(5): 4938-4957. doi: 10.3934/math.2021290 |
[10] | Noureddine Bahri, Abderrahmane Beniani, Abdelkader Braik, Svetlin G. Georgiev, Zayd Hajjej, Khaled Zennir . Global existence and energy decay for a transmission problem under a boundary fractional derivative type. AIMS Mathematics, 2023, 8(11): 27605-27625. doi: 10.3934/math.20231412 |
In this paper, we investigate non-traveling wave solutions of the (3+1)-dimensional variable coefficients Date-Jimbo-Kashiwara-Miwa (VC-DJKM) equation, which describes the real physical phenomena owing to the inhomogeneities of media. By combining the extended homoclinic test approach with variable separation method, we obtain abundant new exact non-traveling wave solutions of the (3+1)-dimensional VC-DJKM equation. These results with a parabolic tail or linear tail reveal the complex structure of the solutions for (3+1)-dimensional VC-DJKM equation. Moreover, the tail in these solutions maybe give a prediction of physical phenomenon. When arbitrary functions contained in these non-traveling wave solutions are taken as some special functions, we can get the kink-type solitons, singular solitary wave solutions, and periodic solitary wave solutions, and so on. As the special cases of our work, the corresponding results of (3+1)-dimensional DJKM equation, (2+1)-dimensional DJKM equation, (2+1)-dimensional VC-DJKM equation are also given.
We consider the system of Dirac equations
ℓy(x):=By′(x)+Q(x)y(x)=λy(x), x∈[a,b], | (1) |
where B=(01−10), Q(x)=(p(x)q(x)q(x)−p(x)), y(x)=(y1(x)y2(x)), p(x), q(x) are real valued functions in L2(a,b) and λ is a spectral parameter, with boundary conditions
U(y):=y2(a)+f1(λ)y1(a)=0 | (2) |
V(y):=y2(b)+f2(λ)y1(b)=0 | (3) |
and with transmission conditions
{y1(wi+0)=αiy1(wi−0)y2(wi+0)=α−1iy2(wi−0)+hi(λ)y1(wi−0)(i=1,2) | (4) |
where fi(λ), hi(λ)(i=1,2) are rational functions of Herglotz-Nevanlinna type such that
fi(λ)=aiλ+bi−Ni∑k=1fikλ−gik | (5) |
hi(λ)=miλ+ni−Pi∑k=1uikλ−tik (i=1,2) | (6) |
ai, bi, fik, gik,mi,ni,uik and tik are real numbers, a1<0, f1k<0, a2>0, f2k>0,mi>0, uik>0 and gi1<gi2<...<giNi, ti1<ti2<...<tiPi, αi>0 and a<w1<w2<b. In special case, when fi(λ)=∞, conditions (2) and (3) turn to Dirichlet conditions y1(a)=y1(b)=0 respectively. Moreover, when hi(λ)=∞, conditions (4) turn to y1(w2+0)=α2y1(w2−0), y2(w2+0)=α−12y2(w2−0)+h2(λ)y1(w2−0) and y1(w1+0)=α1y1(w1−0), y2(w1+0)=α−11y2(w1−0)+h1(λ)y1(w1−0) according to order i=1,2.
Inverse problems of spectral analysis compose of recovering operators from their spectral data. Such problems arise in mathematics, physics, geophysics, mechanics, electronics, meteorology and other branches of natural sciences. Inverse problems also play important role in solving many equations in mathematical physics.
R1(λ)y1(a)+R2(λ)y2(a)=0 is a boundary condition depending spectral parameter where R1(λ) and R2(λ) are polynomials. When degR1(λ)=degR2(λ)=1, this equality depends on spectral parameter as linearly. On the other hand, it is more difficult to search for higher orders of polynomials R1(λ) and R2(λ). When R1(λ)R2(λ) is rational function of Herglotz-Nevanlinna type such that f(λ)=aλ+b−N∑k=1fkλ−gk in boundary conditions, direct and inverse problems for Sturm-Liouville operator have been studied [1,2,3,4,5,6,7,8,9,10,11]. In this paper, direct and inverse spectral problem is studied for the system of Dirac equations with rational function of Herglotz-Nevanlinna in boundary and transmission conditions.
On the other hand, inverse problem firstly was studied by Ambarzumian in 1929 [12]. After that, G. Borg was proved the most important uniqueness theorem in 1946 [13]. In the light of these studies, we note that for the classical Sturm-Liouville operator and Dirac operator, the inverse problem has been studied fairly (see [14,15,16,17,18,19,20], where further references and links to applications can be found). Then, results in these studies have been extended to other inverse problems with boundary conditions depending spectral parameter and with transmission conditions. Therefore, spectral problems for differential operator with transmission conditions inside an interval and with eigenvalue dependent boundary and transmission conditions as linearly and non-linearly have been studied in so many problems of mathematics as well as in applications (see [21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43] and other works, and see [44,45,46,47,48,49,50,51,52,53,54] and other works cited therein respectively).
The aim of this article is to get some uniqueness theorems for mentioned above Dirac problem with eigenvalue dependent as rational function of Herglotz-Nevanlinna type in both of the boundary conditions and also transmission conditions at two different points. We take into account inverse problem for reconstruction of considered boundary value problem by Weyl function and by spectral data {λn,ρn}n∈Z and {λn,μn}n∈Z. Although the boundary and transmission conditions of the problem are not linearly dependent on the spectral parameter, this allows the eigenvalues to be real and to define normalizing numbers.
Consider the space H:=L2(a,b)⊕L2(a,b)⊕CN1+1⊕CN2+1⊕CP1+1⊕ CP2+1 and element Y in H is in the form of Y=(y1(x),y2(x),τ,η,β,γ), such that τ=(Y1,Y2,…,YN1,YN1+1), η=(L1,L2,…,LN2,LN2+1), β=(R1,R2,…,RP1,RP1+1), γ=(V1,V2,…,VP2,VP2+1). H is a Hilbert space with the inner product defined by
<Y,Z>:=b∫a(y1(x)¯z1(x)+y2(x)¯z2(x))dx−YN1+1¯Y′N1+1a1+LN2+1¯L′N2+1a2+α1m1RP1+1¯R′P1+1+α2m2VP2+1¯V′P2+1+N1∑k=1Yk¯Y′k(−1f1k)+N2∑k=1Lk¯L′kf2k+P1∑k=1α1Rr¯R′ru1k+P2∑k=1α2Vr¯V′ru2k | (7) |
for Y=(y1(x),y2(x),τ,η,β,γ) ve Z=(z1(x),z2(x),τ′,η′,β′,γ′) in H. Define the operator T on the domain
D(T)={Y∈H:y1(x),y2(x)∈AC(a,b),
ly∈L2(a,b), y1(w+i)=αiy1(w−i),i=1,2
YN1+1:=−a1y1(a), LN2+1:=−a2y1(b),
RP1+1:=−m1y1(w−1),VP2+1:=−m2y1(w−2)}
such that
TY:=(ly,Tτ,Tη,Tβ,Tγ) | (8) |
where
Tτ=TYi={g1iYi−f1iy1(a), i=¯1,N1y2(a)+b1y1(a)+N1∑k=1Yk, i=N1+1 | (9) |
Tη=TLi={g2iLi−f2iy1(b), i=¯1,N2y2(b)+b2y1(b)+N2∑k=1Lk, i=N2+1 | (10) |
Tβ=TRi={t1iRi−u1iy1(w−1), i=¯1,P1−y2(w+1)+α−11y2(w−1)+n1y1(w−1)+P1∑k=1Rk, i=P1+1 | (11) |
Tγ=TVi={t2iVi−u2iy1(w−2), i=¯1,P2−y2(w+2)+α−12y2(w−2)+n2y1(w−2)+∑k=1P2Vk, i=P2+1. | (12) |
Accordingly, equality TY=λY corresponds to problem (1)-(4) under the domain D(T)⊂H.
Theorem 1. The eigenvalues of the operator T and the problem (1)-(4) coincide.
Proof. Assume that λ is an eigenvalue of T and Y(x)=(y1(x),y2(x),τ,η,β,γ)∈H is the eigenvector corresponding to λ. Since Y∈D(T), it is obvious that the condition y1(wi+0)−αiy1(wi−0)=0 and Eq (1) hold. On the other hand, boundary conditions (2)-(3) and the second condition of (4) are satisfied by the following equalities;
Tτ=TYi=g1i−Yi−f1iy1(a)=λYi, i=¯1,N1
TYN1+1=y2(a)+b1y1(a)+N1∑k=1Yk=−a1y1(a)λ
Tη=TLi=g2iLi−f2iy1(b)=λLi, i=¯1,N2
TLN2+1=y2(b)+b2y1(b)+N2∑k=1Lk=−a2y1(b)λ
Tβ=TRi=t1iRi−u1iy1(w−1), i=¯1,P1
TRP1+1=−y2(w+1)+α−11y2(w−1)+n1y1(w−1)+P1∑k=1Rk=−m1y1(w−1)λ
Tγ=TVi=t2iVi−u2iy1(w−2), i=¯1,P2
TVP2+1=−y2(w+2)+α−12y2(w−2)+n2y1(w−2)+∑k=1P2Vk=−m2y1(w−2)λ.
If λ=gik(i=1,2 and k={1,2,…Ni}) are eigenvalues of operator T, then, from above equalities and the domain of T, equalities (1), y1(a,g1k)=0, y1(b,g2k)=0 and (4) are satisfied.
Moreover, If λ=tik(i=1,2 and k={1,2,…Pi}) are eigenvalues of operator T, from above equalities and the domain of T, Eqs (1)-(3) and y1(w−i,tik)=0=y1(w+i,tik) are valid. In that case, λ is also an eigenvalue of L.
Conversely, let λ be an eigenvalue of L and (y1(x)y2(x)) be an eigenfunction corresponding to λ. If λ≠gik(i=1,2, k={1,2,…Ni}) and λ≠tik(i=1,2, k={1,2,…Pi}) then, it is clear that λ is an eigenvalue of T and the vector
Y=(y1(x),y2(x),f11g11−λy1(a),f12g12−λy1(a),…,f1N1g1N1−λy1(a),−a1y1(a),
f21g21−λy1(b),f22g22−λy1(b),…,f2N2g2N2−λy1(b),−a2y1(b),
u11t11−λy1(w−1),u12t12−λy1(w−1),…,u1P1t1P1−λy1(w−1),−m1y1(w−1),
u21t21−λy1(w−2),u22t22−λy1(w−2),…,u2P2t2P2−λy1(w−2),−m2y1(w−2)) is the eigenvector corresponding to λ.
If λ=g1k(k={1,2,…N1}), then,
Y=(y1(x),y2(x),Y1,Y2,…,YN1,0,L1,L2,…,LN2,LN2+1,R1,R2,…,RP1,RP1+1, V1,V2,…,VP2,VP2+1),
Yi={0, i≠k−y2(a),i=k,i=1,2,…,N1 is the eigenvector of T corresponding to g1k.
If λ=g2k(k={1,2,…N2}), then,
Y=(y1(x),y2(x),Y1,Y2,…,YN1,YN1+1,L1,L2,…,LN2,0,R1,R2,…,RP1,RP1+1,V1,V2,…,VP2,VP2+1), Li={0, i≠k−y2(b),i=k,i=1,2,…,N2 is the eigenvector of T corresponding to g2k.
Furthermore, if λ=t1k(k={1,2,…P1}), then,
Y=(y1(x),y2(x),Y1,Y2,…,YN1,YN1,L1,L2,…,LN2,LN2+1,R1,R2,…,RP1,0,V1,V2,…,VP2,VP2+1), Ri={0, i≠ky2(w+1)−α−11y2(w−1),i=k,i=1,2,…,P1 is the eigenvector corresponding to t1k.
If λ=t2k(k={1,2,…P2}), then, Y=(y1(x),y2(x),Y1,Y2,…,YN1,YN1,L1,L2,…,LN2,LN2+1,R1,R2,…,RP1,RP1+1, V1,V2,…,VP2,0), Vi={0, i≠ky2(w+2)−α−12y2(w−2),i=k,i=1,2,…,P2 is the eigenvector corresponding to t2k.
It is possible to write fi(λ) as follows:
fi(λ)=ai(λ)bi(λ), i=1,2 where
ai(λ)=(aiλ+bi)Ni∏k=1(λ−gik)−Ni∑k=1Ni∏j=1(j≠k)fik(λ−gij)
bi(λ)=Ni∏k=1(λ−gik).
Assume that a2(λ) and b2(λ) do not have common zeros.
Let functions φ(x,λ) and ψ(x,λ) be the solutions of (1) under the initial conditions
φ(a,λ)=(−b1(λ)a1(λ)),ψ(b,λ)=(−b2(λ)a2(λ)) | (13) |
as well as the transmission conditions (4) respectively such that
φ(x,λ)={φ1(x,λ), x<w1φ2(x,λ), w1<x<w2φ3(x,λ), w2<x<b and ψ(x,λ)={ψ3(x,λ), x<w1ψ2(x,λ), w1<x<w2ψ1(x,λ), w2<x<b.
Then it can be easily proven that φi(x,λ) and ψi(x,λ), i=¯1,3 are the solutions of the following integral equations;
φi+1,1(x,λ)=αiφi1(wi,λ)cosλ(x−wi)
−[α−1iφi2(wi,λ)+hi(λ)φi1(wi,λ)]sinλ(x−wi)
+x∫wi[p(t)sinλ(x−t)+q(t)cosλ(x−t)]φi+1,1(t,λ)dt
+x∫wi[q(t)sinλ(x−t)−p(t)cosλ(x−t)]φi+1,2(t,λ)dt,
φi+1,2(x,λ)=αiφi1(wi,λ)sinλ(x−wi)+[α−1iφi2(wi,λ)+hi(λ)φi1(wi,λ)]cosλ(x−wi)+x∫wi[−p(t)cosλ(x−t)+q(t)sinλ(x−t)]φi+1,1(t,λ)dt+x∫wi[−q(t)cosλ(x−t)−p(t)sinλ(x−t)]φi+1,2(t,λ)dt,for i=1,2 |
and
ψi1(x,λ)=α−1iψi+1,1(wi,λ)cosλ(x−wi)+(−αiψi+1,2(wi,λ)+hi(λ)ψi+1,1(wi,λ))sinλ(x−wi)−wi∫x[p(t)sinλ(x−t)+q(t)cosλ(x−t)]ψi1(t,λ)dt+wi∫x[−q(t)sinλ(x−t)+p(t)cosλ(x−t)]ψi2(t,λ)dtψi2(x,λ)=α−1iψi+1,1(wi,λ)sinλ(x−wi)+(αiψi+1,2(wi,λ)−hi(λ)ψi+1,1(wi,λ))cosλ(x−wi)+wi∫x[p(t)cosλ(x−t)−q(t)sinλ(x−t)]ψi1(t,λ)dt+w2∫x[q(t)cosλ(x−t)+p(t)sinλ(x−t)]ψi2(t,λ)dt,for i=2,1 |
Lemma 1. For the solutions φi(x,λ) and ψi(x,λ), i=¯1,3 as |λ|→∞, the following asymptotic estimates hold;
φ11(x,λ)={a1λN1+1sinλ(x−a)+o(|λ|N1+1exp|Imλ|[(x−a)]),
φ12(x,λ)={a1λN1+1cosλ(x−a)+o(|λ|N1+1exp|Imλ|[(x−a)]),
φ21(x,λ)={a1m1λL1+N1+2sinλ(w1−a)sinλ(x−w1)+o(|λ|L1+N1+2exp|Imλ|[(w1−a)+(x−w1)])
φ22(x,λ)={a1m1λL1+N1+2sinλ(w1−a)cosλ(x−w1)+o(|λ|L1+N1+2exp|Imλ|[(w1−a)+(x−w1)])
φ31(x,λ)={−m2m1a1λL1+L2+N1+3sinλ(w1−a)sinλ(w2−w1)sinλ(x−w2)+o(|λ|L1+L2+N1+3exp|Imλ|[(w1−a)+(w2−w1)+(x−w2)])
φ32(x,λ)={m2m1a1λL1+L2+N1+3sinλ(w1−a)sinλ(w2−w1)cosλ(x−w2)+o(|λ|L1+L2+N1+3exp|Imλ|[(w1−a)+(w2−w1)+(x−w2)])
ψ11(x,λ)={−a2λN2+1sinλ(x−b)+o(|λ|N2+1exp|Imλ|[(x−b)])
ψ12(x,λ)={a2λN2+1cosλ(x−b)+o(|λ|N2+1exp|Imλ|[(x−b)])
ψ21(x,λ)={−m2a2λN2+L2+2sinλ(w2−b)sinλ(x−w2)+o(|λ|N2+L2+2exp|Imλ|[(w2−b)+(x−w2)])
ψ22(x,λ)={m2a2λN2+L2+2sinλ(w2−b)cosλ(x−w2)+o(|λ|N2+L2+2exp|Imλ|[(w2−b)+(x−w2)])
ψ31(x,λ)={−m1m2a2λN2+L1+L2+3sinλ(w2−b)sinλ(w1−w2)sinλ(x−w1)+o(|λ|N2+L1+L2+3exp|Imλ|[(w2−b)+(w1−w2)+(x−w2)])
ψ32(x,λ)={m1m2a2λN2+L1+L2+3sinλ(w2−b)sinλ(w1−w2)cosλ(x−w1)+o(|λ|N2+L1+L2+3exp|Imλ|[(w2−b)+(w1−w2)+(x−w1)])
Theorem 2. The eigenvalues {λn}n∈Z of problem L are real numbers.
Proof. It is enough to prove that eigenvalues of operator T are real. By using inner product (7), for Y in D(T), we compute that
⟨TY,Y⟩=b∫alyˉydx−1a1TYN1+1¯YN1+1+1a2TLN2+1¯LN2+1+α1m1TRP1+1¯RP1+1+α2m2TVP2+1¯VP2+1N1−∑k=1TYk¯Yk(1f1k)+N2∑k=1TLk¯Lk(1f2k)+P1∑k=1α1TRk¯Rk(1u1k)+P2∑k=1α2TVk¯Vk(1u2k). |
If necessary arrangements are made, we get
⟨TY,Y⟩=b∫ap(x)(|y1|2−|y2|2)dx+b∫aq(x)2Re(y2¯y1)dx+b1|y1(a)|+N1∑k=12Re(Yk¯y1(a))−b2|y1(b)|2−N2∑k=12Re(Lk¯y1(b))−a1n1|y1(w−1)|2−P1∑k=1a12Re(Rky1(w−1))−a2n2|y1(w−2)|2−P2∑k=1a22Re(Vky1(w−2))−N1∑k=1g1k|Yk|21f1k+N2∑k=1g2kf2k|Lk|2+P1∑k=1a1t1ku1k|Rk|2+P2∑k=1a2t2ku2k|Vk|2−b∫a2Re(y2¯y1′)dx. |
Accordingly, since ⟨TY,Y⟩ is real for each Y in D(T), λ∈R is obtained.
Lemma 2. The equality ‖Yn‖2=ρn is valid such that Yn is eigenvector corresponding to eigenvalue λn of T.
Proof. Let λn≠gik. When λn=gik, following proof is done with minor changes. By using the structure of D(T) and the Eqs (8)-(12), we get
‖Yn‖2=b∫a(φ21(x,λn)+φ22(x,λn))dx−|YN1+1|2a1+|LN2+1|2a2+α1m1|RP1+1|2+α2m2|VP2+1|2−N1∑k=1|Yk|2f1k+N2∑k=1|Lk|2f2k+P1∑k=1α1u1k|Rk|2+P2∑k=1α2u2k|Vk|2 | (14) |
=b∫a(φ21(x,λn)+φ22(x,λn))dx−a1φ21(a,λn)+a2φ21(b,λn)+m1α1φ21(w1−0,λn) +m2α2φ21(w2−0,λn)−N1∑k=1f1kφ21(a,λn)(λn−g1k)2+N2∑k=1f2kφ21(b,λn)(λn−g2k)2 +P1∑k=1α1u1kφ21(w1−0,λn)(λn−t1k)2+P2∑k=1α2u2kφ21(w2−0,λn)(λn−t2k)2 =b∫a(φ21(x,λn)+φ22(x,λn))dx−φ21(a,λn)(a1+N1∑k=1f1k(λn−g1k)2) +φ21(b,λn)(a2+N2∑k=1f2k(λn−g2k)2)+α1φ21(w1−0,λn)(m1+P1∑k=1u1k(λn−t1k)2) +α2φ21(w2−0,λn)(m2+P2∑k=1u2k(λn−t2k)2) =b∫a(φ21(x,λn)+φ22(x,λn))dx−φ21(a,λn)f′1(λn)+φ21(b,λn)f′2(λn) +α1φ21(w−1,λn)h′1(λn)+α2φ21(w−2,λn)h′2(λn)=ρn. |
On the other hand, the expression
W(φ,ψ)=φ1(x,λ)ψ2(x,λ)−φ2(x,λ)ψ1(x,λ) |
is called characteristic function of problem (1)-(4). Moreover, since solutions φ(x,λ) and ψ(x,λ) satisfy the problem L,
for ∀x∈[a,b]
∂∂xW(φ,ψ) =φ′1(x,λ)ψ2(x,λ)+ψ′2(x,λ)φ1(x,λ)−φ′2(x,λ)ψ1(x,λ)−ψ′1(x,λ)φ2(x,λ) =[q(x)φ1(x,λ)−p(x)φ2(x,λ)−λφ2(x,λ)]ψ2(x,λ) +[−p(x)ψ1(x,λ)−q(x)ψ2(x,λ)+λψ1(x,λ)]φ1(x,λ) −[−p(x)φ1(x,λ)−q(x)φ2(x,λ)+λφ1(x,λ)]ψ1(x,λ) −[q(x)ψ1(x,λ)−p(x)ψ2(x,λ)−λψ2(x,λ)]φ2(x,λ)=0 |
is obtained. Furthermore, since solutions φ(x,λ) and ψ(x,λ) also satisfy transmission conditions (4), we get
W(wi+0)=φ1(wi+0,λ)ψ2(wi+0,λ)−φ2(wi+0,λ)ψ1(wi+0,λ)=αiφ1(wi−0,λ)[α−1iψ2(wi−0,λ)+hi(λ)ψ1(wi−0,λ)]−[α−1iφ2(wi−0,λ)+hi(λ)φ1(wi−0,λ)]αiψ1(wi−0,λ)=φ1(wi−0,λ)ψ2(wi−0,λ)−φ2(wi−0,λ)ψ1(wi−0,λ)=W(wi−0). |
Therefore, since characteristic function W(φ,ψ) is independent from x,
W{φ,ψ}:=Δ(λ)=φ1(x,λ)ψ2(x,λ)−φ2(x,λ)ψ1(x,λ)=a2(λ)φ1(b,λ)+b2(λ)φ2(b,λ)=−b1(λ)ψ2(a,λ)−a1(λ)ψ1(a,λ) |
can be written.
It is clear that Δ(λ) is an entire function and its zeros namely {λn}n∈Zcoincide with the eigenvalues of the problem L.
Accordingly, for each eigenvalue λn equality ψ(x,λn)=snφ(x,λn) is valid where sn=ψ1(a,λn)−b1(λn)=ψ2(a,λn)a1(λn).
On the other hand, since ai(gik)≠0 ve bi(gik)=0 for ∀i∈{1,2} and k={1,2,…,Ni}, gik is an eigenvalue if and only if φ1(b,g2k)=0, φ1(a,g1k)=0 i.e., Δ(gik)=0.
At the same time, tik is an eigenvalue if and only if φ1(w−i,tik)=0=φ1(w+i,tik) i.e., Δ(tik)=0 such that i=1,2 and k={1,2,…Pi}.
Theorem 3. Eigenvalues of problem L are simple.
Proof. Let λn≠gik and φ(x,λn) be eigenfunction corresponds to the eigenvalue λn. In that case, the Eq (1) can be written for ψ(x,λ) and φ(x,λn) as follows;
Bψ′(x,λ)+Q(x)ψ(x,λ)=λψ(x,λ)Bφ′(x,λn)+Q(x)φ(x,λn)=λnφ(x,λn). |
If we multiply these equations by φ(x,λn) and ψ(x,λ) respectively and add side by side, we get the following equality;
(ψ2(x,λ)φ1(x,λn)−ψ1(x,λ)φ2(x,λn))′=(λ−λn)(ψ1(x,λ)φ1(x,λn)+ψ2(x,λ)φ2(x,λn)). |
Then if last equality is integrated over the interval [a,b] and the initial conditions (13) and transmission conditions (4) are used to get
b∫a(ψ1(x,λ)φ1(x,λn)+ψ2(x,λ)φ2(x,λn))dx+α2φ1(w−2,λn)ψ1(w−2,λ)h2(λ)−h2(λn)λ−λn+α1ψ1(w−1,λ)φ1(w−1,λn)h1(λ)−h1(λn)λ−λn+ψ1(b,λn)φ1(b,λn)f2(λ)−f2(λn)λ−λn−φ1(a,λn)ψ1(a,λ)f1(λ)−f1(λn)λ−λn=−(Δ(λ)−Δ(λn)(λ−λn)). |
Then, considering that ψ(x,λn)=snφ(x,λn)
if the limit is passed when λ→λn, snρn=−˙Δ(λn) is obtained.
If g1k and g2k are non-simple eigenvalues then φ1(a,g1k)=0, φ1(b,g2k)=0 and so b∫a(φ21(x,λn)+φ22(x,λn))dx=−[α1φ21(w−1,λn)h′1(λn)+α2φ21(w−2,λn)h′2(λn)] is obtained. Since α1, α2 and for all λn, h′1(λn), h′2(λn) are positive, we have a contradiction. Therefore, eigenvalues gik are also simple.
Using expressions a2(λ), b2(λ) and asymptotic behaviour of solution φ(x,λ), we obtain the following asymptotic of characteristic function Δ(λ) as |λ|→∞; Δ(λ)=−a1a2m1m2λN1+N2+L1+L2+4sinλ(w1−a)sinλ(w2−w1)sinλ(b−w2)+o(|λ|N1+N2+L1+L2+4e|Imλ|(b−a)).
Let Φ(x,λ):=(Φ1(x,λ)Φ2(x,λ)) be the solution of Eq (1) under the conditions U(Φ)=1, V(Φ)=0 as well as the transmission conditions (4).
Since V(Φ)=0=V(ψ), it can be supposed that Φ(x,λ)=kψ(x,λ) (k≠0) where k is a constant.
W(φ,Φ)=φ1(x,λ)Φ2(x,λ)−φ2(x,λ)Φ1(x,λ)|x=a=−b1(λ)Φ2(a,λ)−a1(λ)Φ1(a,λ)=−U(Φ)=−1. |
By the relation U(Φ)=1, we get k[b1(λ)ψ2(a,λ)+a1(λ)ψ1(a,λ)]=1. Since U(ψ)=−Δ(λ), we obtain Φ(x,λ)=kψ(x,λ)=−ψ(x,λ)Δ(λ) for λ≠λn.
Let S(x,λ)=(S1(x,λ)S2(x,λ)) and C(x,λ)=(C1(x,λ)C2(x,λ)) be solutions of (1) satify the conditions S(a,λ)=(01), C(a,λ)=(10) and transmission conditions (4).
Accordingly, the following equalities are obtained:
φ1(x,λ)=−b1(λ)C(x,λ)+a1(λ)S(x,λ) | (15) |
Φ(x,λ)=1b1(λ)(S(x,λ)−Φ1(a,λ)φ(x,λ)). | (16) |
The function Φ(x,λ) is called Weyl solution and the function M(λ)=−Φ1(a,λ) is called Weyl function of problem L. Therefore, since Φ(x,λ)=−ψ(x,λ)Δ(λ), we set M(λ):=ψ1(a,λ)Δ(λ).
Consider the boundary value problem ˜L in the same form with L but different coefficients. Here, the expressions related to the L problem are shown with s and the ones related to ˜L are shown with ˜s. According to this statement, we set the problem ˜L as follows:
˜ℓ[y(x)]:=By′(x)+˜Q(x)y(x)=λy(x), x∈[a,b] |
˜U(y):=y2(a)+˜f1(λ)y1(a)=0˜V(y):=y2(b)+˜f2(λ)y1(b)=0y1(wi+0)=˜αiy1(wi−0)y2(wi+0)=−˜α−1iy2(wi−0)+˜hi(λ)y1(wi−0) |
where ˜Q(x)=(˜p(x)q(x)q(x)−˜p(x)).
Theorem 4. If M(λ)=˜M(λ), f1(λ)=˜f1(λ), then Q(x)=˜Q(x) almost everywhere in (a,b), f2(λ)=˜f2(λ), hi(λ)=˜hi(λ),and αi(λ)=˜αi(λ) (i=1,2).
Proof. Introduce a matrix P(x,λ)=[Pij(x,λ)]i,j=1,2 by the equality as follows;
(P11P12P21P22)(˜φ1˜Φ1˜φ2˜Φ2)=(φ1Φ1φ2Φ2).
According to this, we get
P11(x,λ)=−φ1(x,λ)˜Φ2(x,λ)+Φ1(x,λ)˜φ2(x,λ)P12(x,λ)=−˜φ1(x,λ)Φ1(x,λ)+φ1(x,λ)˜Φ1(x,λ)P21(x,λ)=−φ2(x,λ)˜Φ2(x,λ)+Φ2(x,λ)˜φ2(x,λ)P22(x,λ)=−˜φ1(x,λ)Φ2(x,λ)+φ2(x,λ)˜Φ1(x,λ) | (17) |
or by using the relation Φ(x,λ)=−ψ(x,λ)Δ(λ),
we obtain
P11(x,λ)=φ1(x,λ)˜ψ2(x,λ)˜Δ(λ)−˜ψ2(x,λ)ψ1(x,λ)Δ(λ)P12(x,λ)=−φ1(x,λ)˜ψ1(x,λ)˜Δ(λ)+˜φ1(x,λ)ψ1(x,λ)Δ(λ)P21(x,λ)=φ2(x,λ)˜ψ2(x,λ)˜Δ(λ)−˜φ2(x,λ)ψ2(x,λ)Δ(λ)P22(x,λ)=˜φ1(x,λ)ψ2(x,λ)Δ(λ)−φ2(x,λ)˜ψ1(x,λ)˜Δ(λ). | (18) |
Taking into account the Eqs (15) and (16) and M(λ)=˜M(λ), we can easily get
P11(x,λ)=C1(x,λ)˜S2(x,λ)−S1(x,λ)˜C2(x,λ)P12(x,λ)=˜C1(x,λ)S1(x,λ)−C1(x,λ)˜S1(x,λ)P21(x,λ)=C2(x,λ)˜S2(x,λ)−S2(x,λ)˜C2(x,λ)P22(x,λ)=˜C1(x,λ)S2(x,λ)−C2(x,λ)˜S1(x,λ). |
Hence, the functions Pij(x,λ) are entire in λ. Denote
Gδ:={λ:|λ−λn|≥δ, n=0,±1,±2,…},δ>0 and
˜Gδ:={λ:|λ−˜λn|≥δ, n=0,±1,±2,…} where δ>0 is sufficiently small and fixed.
Clearly, for λ∈Gδ∩˜Gδ, |sinλx|≥Cδe|Imλ|x, |λ|→∞.
Therefore, |Δ(λ)|≥CδλN1+N2+L1+L2+4e|Imλ|(b−a), \ λ∈Gδ∩˜Gδ, |λ|≥λ∗ for sufficiently large λ∗=λ∗(δ) and from (18) we see that Pij(x,λ) are bounded with respect to λ where λ∈Gδ∩˜Gδ and |λ| sufficiently large. From Liouville's theorem, it is obtained that these functions do not depend on λ.
On the other hand, from (18)
P11(x,λ)−1=φ1(x,λ)(˜ψ2(x,λ)˜Δ(λ)−ψ2(x,λ)Δ(λ))−ψ1(x,λ)Δ(λ)(˜φ2(x,λ)−φ2(x,λ))
P12(x,λ)=˜φ1(x,λ)(ψ1(x,λ)Δ(λ)−˜ψ1(x,λ)˜Δ(λ))−˜ψ1(x,λ)˜Δ(λ)(φ1(x,λ)−˜φ1(x,λ))
P21(x,λ)=φ2(x,λ)(˜ψ2(x,λ)˜Δ(λ)−ψ2(x,λ)Δ(λ))−ψ2(x,λ)Δ(λ)(˜φ2(x,λ)−φ2(x,λ))
P22(x,λ)−1=ψ2(x,λ)Δ(λ)(˜φ1(x,λ)−φ1(x,λ))−φ2(x,λ)(˜ψ1(x,λ)˜Δ(λ)−ψ1(x,λ)Δ(λ)).
If it is considered that Pij(x,λ) do not depend on λ and asymptotic formulas of solutions φ(x,λ) and ψ(x,λ), we obtain
limλ→−∞φ1(x,λ)(˜ψ2(x,λ)˜Δ(λ)−ψ2(x,λ)Δ(λ))=0,
limλ→−∞ψ1(x,λ)Δ(λ)(˜φ2(x,λ)−φ2(x,λ))=0
for all x in [a,b]. Hence, limλ→−∞[P11(x,λ)−1]=0.
Thus, P11(x,λ)=1 and similarly, P22(x,λ)=1 and P12(x,λ)=P21(x,λ)=0.
Substitute these relations in (17), to obtain
φ1(x,λ)=˜φ1(x,λ), ψ1(x,λ)Δ(λ)=˜ψ1(x,λ)˜Δ(λ)
φ2(x,λ)=˜φ2(x,λ), ψ2(x,λ)Δ(λ)=˜ψ2(x,λ)˜Δ(λ) for all x and λ.
Taking into account these results and Eq (1), we have
(Q(x)−˜Q(x))φ(x,λ)=0. |
Therefore, Q(x)=˜Q(x) i.e., p(x)=˜p(x). Moreover, it is considered that
ψ1(x,λ)Δ(λ)=˜ψ1(x,λ)˜Δ(λ), ψ2(x,λ)Δ(λ)=˜ψ2(x,λ)˜Δ(λ)
and
b2(λ)ψ2(x,λ)+a2(λ)ψ1(x,λ)=0
˜b2(λ)˜ψ2(x,λ)+˜a2(λ)˜ψ1(x,λ)=0
we get a2(λ)˜b2(λ)−b2(λ)˜a2(λ)=0. As we have said above, a2(λ), b2(λ) as well as ˜a2(λ), ˜b2(λ) do not have common zeros. Hence, a2(λ)=˜a2(λ), b2(λ)=˜b2(λ), i.e., f2(λ)=˜f2(λ).
On the other hand, substituting φ1 and φ2 into transmission conditions (4), we get
φ1(w+i,λ)=αiφ1(w−i,λ), ˜φ1(w+i,λ)=˜ai˜φ1(w−i,λ)
φ2(w+i,λ)=α−1iφ2(w−1,λ)+hi(λ)φ1(w−i,λ),
˜φ2(w+i,λ)=˜α−1i˜φ2(w−i,λ)+˜hi(λ)˜φ1(w−i,λ), i=1,2.
Therefore, since φ1(x,λ)=˜φ1(x,λ), φ2(x,λ)=˜φ2(x,λ), these yield that α1=˜α1, α2=˜α2
and h1(λ)=˜h1(λ), h2(λ)=˜h2(λ).
Theorem 5. If {λn,ρn}n∈Z={˜λn,˜ρn}n∈Z, f1(λ)=˜f1(λ) then Q(x)=˜Q(x) almost everywhere in (a,b), f2(λ)=˜f2(λ), hi(λ)=˜hi(λ),and αi(λ)=˜αi(λ) (i=1,2).
Proof. Since λn=˜λn, Δ(λ)=c˜Δ(λ). On the other hand, also since snρn=−˙Δ(λn) and ρn=˜ρn, we get that sn=c˜sn. Therefore, ψ1(a,λn)=c˜ψ1(a,λn) is obtained.
Denote H(λ):=ψ1(a,λ)−c˜ψ1(a,λ)Δ(λ) which is an entire function in λ. Since lim|λ|→∞H(λ)=0, H(λ)≡0 and so ψ1(a,λ)=c˜ψ1(a,λ). Hence, M(λ)=˜M(λ). As a result, the proof of theorem is finished by Theorem 4.
We examine the boundary value problem L1 with the condition y1(a)=0 instead of (2) in problem L. Let {μn}n∈Z be eigenvalues of the problem L1. It is clear that {μn}n∈Z are zeros of Δ1(μ):=−ψ1(a,μ).
Theorem 6. If {λn,μn}n∈Z={˜λn,˜μn}n∈Z, f1(λ)=˜f1(λ) and K=˜K such that K=a2m1m2, ˜K=˜a2˜m1˜m2 then Q(x)=˜Q(x) almost everywhere in (a,b), f2(λ)=˜f2(λ), hi(λ)=˜hi(λ),and αi(λ)=˜αi(λ) (i=1,2).
Proof. Since for all n∈Z, λn=˜λn and μn=˜μn, Δ(λ)˜Δ(λ) and Δ1(μ)˜Δ1(μ) are entire functions in λ and in μ respectively. On the other hand, taking into account the asymptotic behaviours of Δ(λ), Δ1(μ) and K=˜K, we obtain limλ→−∞Δ(λ)˜Δ(λ)=1 and limμ→−∞Δ1(μ)˜Δ1(μ)=1. Therefore, since λn=˜λn and μn=˜μn, we get Δ(λ)=˜Δ(λ) and Δ1(μ)=˜Δ1(μ). If we consider the case Δ1(μ)=˜Δ1(μ), then ψ1(a,μ)=˜ψ1(a,μ) is obtained. Furthermore, since M(λ)=ψ1(a,λ)Δ(λ), M(λ)=˜M(λ). Hence, the proof is completed by Theorem 4.
The purpose of this paper is to state and prove some uniqueness theorems for Dirac equations with boundary and transmission conditions depending rational function of Herglotz-Nevanlinna. Accordingly, it has been proved that while f1(λ) in condition (2) is known, the coefficients of the boundary value problem (1)-(4) can be determined uniquely by each of the following;
i) The Weyl function M(λ)
ii) Spectral data {λn,ρn} forming eigenvalues and normalizing constants respectively
iii) Two given spectra {λn,μn}
These results are the application of the classical uniqueness theorems of Marchenko, Gelfand, Levitan and Borg to such Dirac equations. Considering this study, similar studies can be made for classical Sturm-Liouville operators, the system of Dirac equations and diffusion operators with finite number of transmission conditions depending spectral parameter as Herglotz-Nevanlinna function.
There is no conflict of interest.
[1] | Z. Z. Lan, Periodic, breather and rogue wave solutions for a generalized (3+1)-dimensional variable-coefficient B-type Kadomtsev-Petviashvili equation in fluid dynamics, Appl. Math. Lett., 94 (2019), 126–132. |
[2] | M. H. Huang, M. A. S. Murad, O. A. Ilhan, One-, two- and three-soliton, periodic and cross-kink solutions to the (2+1)-D variable-coefficient KP equation, Modern Phys. Lett. B, 34 (2020), 2050045. |
[3] |
X. Y. Gao, Y. J. Guo, W. R. Shan, Magneto-optical/ferromagnetic-material computation: Bäcklund transformations, bilinear forms and N solitons for a generalized (3+1)-dimensional variable-coefficient modified Kadomtsev-Petviashvili system, Appl. Math. Lett., 111 (2021), 106627. doi: 10.1016/j.aml.2020.106627
![]() |
[4] |
M. Arshad, A. R. Seadawy, D. Lu, Modulation stability and optical soliton solutions of nonlinear Schrödinger equation with higher order dispersion and nonlinear terms and its applications, Superlattices Microst., 112 (2017), 422–434. doi: 10.1016/j.spmi.2017.09.054
![]() |
[5] |
A. R. Seadawy, K. El-Rashidy, Dispersive solitary wave solutions of Kadomtsev-Petviashvili and modified Kadomtsev-Petviashvili dynamical equations in unmagnetized dust plasma, Results Phys., 8 (2018), 1216–1222. doi: 10.1016/j.rinp.2018.01.053
![]() |
[6] |
Y. S. Özkan, E. Yaşar, On the exact solutions of nonlinear evolution equations by the improved tan(φ/2)-expansion method, Pramana J. Phys., 94 (2020), 37. doi: 10.1007/s12043-019-1883-3
![]() |
[7] |
M. Iqbal, A. R. Seadawy, O. H. Khalil, D. Lu, Propagation of long internal waves in density stratified ocean for the (2+1)-dimensional nonlinear Nizhnik-Novikov-Vesselov dynamical equation, Results Phys., 16 (2020), 102838. doi: 10.1016/j.rinp.2019.102838
![]() |
[8] |
X. M. Peng, Y. D. Shang, X. X. Zheng, New non-travelling wave solutions of Calogero equation, Adv. Appl. Math. Mech., 8 (2016), 1036–1049. doi: 10.4208/aamm.2015.m1121
![]() |
[9] |
L. B. Lv, Y. D. Shang, Abundant new non-travelling wave solutions for the (3+1)-dimensional potential-YTSF equation, Appl. Math. Lett., 107 (2020), 106456. doi: 10.1016/j.aml.2020.106456
![]() |
[10] | J. Wang, H. L. An, B. Li, Non-traveling lump solutions and mixed lump-kink solutions to (2+1)-dimensional variable-coefficient Caudrey-Dodd-Gibbon-Kotera-Sawada equation, Modern Phys. Lett. B, 33 (2019), 1950262. |
[11] |
S. M. Guo, L. Q. Mei, Y. B. Zhou, The compound (G′/G)-expansion method and double non-traveling wave solutions of (2+1)-dimensional nonlinear partial differential equations, Comput. Math. Appl., 69 (2015), 804–816. doi: 10.1016/j.camwa.2015.02.016
![]() |
[12] |
J. G. Liu, Y. Tian, J. G. Hu, New non-traveling wave solutions for the (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equation, Appl. Math. Lett., 79 (2018), 162–168. doi: 10.1016/j.aml.2017.12.011
![]() |
[13] |
A. M. Wazwaz, New (3+1)-dimensional Date-Jimbo-Kashiwara-Miwa equations with constant and time-dependent coefficients: Painlevé integrability, Phys. Lett. A, 384 (2020), 126787. doi: 10.1016/j.physleta.2020.126787
![]() |
[14] |
Y. Q. Yuan, B. Tian, W. R. Sun, J. Chai, L. Liu, Wronskian and Grammian solutions for a (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation, Comput. Math. Appl., 74 (2017), 873–879. doi: 10.1016/j.camwa.2017.06.008
![]() |
[15] |
L. Cheng, Y. Zhang, M. J. Lin, Lax pair and lump solutions for the (2+1)-dimensional DJKM equation associated with bilinear Bäcklund transformations, Anal. Math. Phys., 9 (2019), 1741–1752. doi: 10.1007/s13324-018-0271-3
![]() |
[16] |
Y. H. Wang, H. Wang, C. L. Temuer, Lax pair, conservation laws, and multi-shock wave solutions of the DJKM equation with Bell polynomials and symbolic computation, Nonlinear Dynam., 78 (2014), 1101–1107. doi: 10.1007/s11071-014-1499-6
![]() |
[17] |
A. R. Adem, Y. Yildirim, E. Yaşar, Complexiton solutions and soliton solutions: (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation, Pramana-J. Phys., 92 (2019), 36. doi: 10.1007/s12043-018-1707-x
![]() |
[18] |
Z. Z. Kang, T. C. Xia, Construction of abundant solutions of the (2+1)-dimensional time-dependent Date-Jimbo-Kashiwara-Miwa equation, Appl. Math. Lett., 103 (2020), 106163. doi: 10.1016/j.aml.2019.106163
![]() |
[19] |
A. M. Wazwaz, A (2+1)-dimensional time-dependent Date-Jimbo-Kashiwara-Miwa equation: Painlevé integrability and multiple soliton solutions, Comput. Math. Appl., 79 (2020), 1145–1149. doi: 10.1016/j.camwa.2019.08.025
![]() |
1. | Mehmet Kayalar, A uniqueness theorem for singular Sturm-liouville operator, 2023, 34, 1012-9405, 10.1007/s13370-023-01097-x |