Skip to main content
Have a personal or library account? Click to login
The Resolvent of Impulsive Singular Hahn–Sturm–Liouville Operators Cover

The Resolvent of Impulsive Singular Hahn–Sturm–Liouville Operators

Open Access
|Jan 2024

Full Article

1. Introduction

Impulsive differential equations are one of the interesting topics in the theory of differential equations. These equations serve as basic models to study the dynamics of processes that are subject to sudden changes in their states. These types of problems are especially encountered in heat and mass transfer problems ([17]). There are many studies on this subject in the literature [2, 6, 7, 8, 9, 12, 13, 14, 19, 23].

W. Hahn introduced the concept of the Hahn derivative to the literature in 1949 [10]. With this definition he made, he gathered two important operators under a single structure. These are the q-difference and forward difference operators. In 2018, Annaby et al. [4] using this definition instead of the classical derivative, investigated the fundamental properties of the Sturm–Liouville problems. In [5], the authors studied singular q-Sturm-Liouville equations. In [18], the author studied a q-analog of the singular Dirac problem. Recently in [21], the author proved a spectral expansion theorem by constructing the spectral function of the Hahn–Sturm–Liouville equation in the singular case under impulsive conditions.

In this paper, our aim is to consider Hahn–Sturm–Liouville problems under impulsive boundary conditions. The integral representation of the resolvent operator corresponding to this type of problem will be obtained using Weyl’s method [16, 22, 24].

2. Preliminaries

Now, we provide a concise overview of the Hahn calculus [3, 4, 10, 11]. Let q ∈ (0, 1), ω0 := ω/ (1 − q) , ω > 0, and let Ψ : J ⊂ ℝ ω ℝ be a function such that ω0J.

Definition 2.1 ([10], [11]).

The Hahn derivative Dω,q Ψ is defined by

Dω,qΨζ=Ψω+qζΨζω+q1ζ,ζω0,Ψω0 ,ζ=ω0,
where the expression Ψ (ω0) shows the ordinary derivative of Ψ at ω0.

Definition 2.2 ([3]).

Let a, b, ω0J. The Hahn integral (ω, q-integral) is defined by

abΨζdω,qζ:=ω0bΨζdω,qζω0aΨζdω,qζ,
where
ω0ζΨ(t)dω,qt:=((1-q)ζ-ω)n=0qnΨ(ω1-qn1-q+ζqn),ζJ,
provided that the series converges at ζ = a and ζ = b.

Definition 2.3.

The ω, q-Wronskian of Ψ1 and Ψ2 is defined by

Wω,qΨ1, Ψ2:=Ψ1Dω,qΨ2Ψ2Dω,qΨ1.

3. Main results

Let us consider the following impulsive boundary-value problem (BVP)

(3.1)
1qDωq,1qDω,qyζ+vζyζ=λyζ,ζω0, dd, 1qn,
(3.2)
yω0, λ cos β+Dωq,1qyω0, λ sin β=0,
(3.3)
yd=ηyd+,
(3.4)
Dωq,1qyd=1ηDωq,1qyd+,
(3.5)
y1qn, λ cos γ+Dωq,1qy1qn, λ sin γ=0,
where q ∈ (0, 1), ω0 := ω/ (1 − q), ω > 0, γ, β ∈ ℝ, 1qn>d , n ∈ ℕ := {1, 2, 3, . . .}, η > 0, λ ∈ ℂ, y() := limζd± y (ζ) , v is a real-valued continuous function on [ω0, d) ∪ (d, ∞), and has finite limits v().

A similar problem has been studied by the authors without impulsive boundary conditions ([1]).

Hn=Lω,q2ω0, d+˙Lω,q2d, 1qn , 1qn>d , nH=Lω,q2ω0, d+˙Lω,q2d,  is a Hilbert space endowed with the following inner product

y, zn:=ω0dy1z1¯dω,qζ+d1qny2z2¯dω,qζ,(y, z:=ω0dy1z1¯dω,qζ+dy2z2¯dω,qζ)
where
yζ=y1ζ , ζω0, d ,y2ζ , ζd,  ,
and
zζ=z1ζ, ζω0, d,z2ζ, ζd, .

Let

ψζ, λ=ψ1ζ, λ, ζω0, d,ψ2ζ, λ, ζd, ,
and
θζ, λ=θ1ζ, λ, ζω0, d,θ2ζ, λ, ζd, ,
be solutions of Eq. (3.1) satisfying the following conditions
ψ1ω0, λ= cos β, Dωq,1qψ1ω0, λ= sin β,θ1ω0, λ= sin β, Dωq,1qθ1ω0, λ= cos β,
and
θd, λ=ηθd+, λ,Dωq,1qθd, λ=1ηDωq,1qθd+, λ,ψd, λ=ηψd+, λ,Dωq,1qψd, λ=1ηDωq,1qψd+, λ.

Then the solution of Eq. (3.1) be represented

ψζ, λ+𝓁λ, 1qnθζ, λ
which satisfies the boundary condition
Dωq,1qψ21qn, λ+𝓁λ, 1qnDωq,1qθ21qn, λ sin γ+ψ21qn, λ+𝓁λ, 1qnθ21qn, λ cos γ=0.
Hence
𝓁λ, 1qn=ψ21qn,λ cot γ+Dωq,1qψ21qn,λθ21qn,λ cot γ+Dωq,1qθψ21qn,λ.

Lemma 3.1.

Let

Z1qnζ, λ=ψζ, λ+𝓁λ, 1qnθζ, λ,
where Z1qnHn and 1qn>d , n ∈ ℕ. Then, for each nonreal λ, the following relations hold:
Z1qnζ, λZζ, λ , n,ω0dZ1qnζ, λ2dω,qζ+d1qnZ1qnζ, λ2dω,qζω0dZζ, λ2dω,qζ+dζ, λ2dω,qζ, n.

Proof

It is immediate that

Z1qnζ, λ=Zζ, λ+𝓁λ, 1qnmλθζ, λ
where Z(·, λ) ∈ H and m(λ) is the Titchmarsh–Weyl function. 𝓁λ, 1qn varies on a circle with a finite radius r1qn in the plane. In the limit-circle case, 𝓁λ, 1qnmλn ; therefore
Z1qnζ, λZζ, λ n.

Hence

ω0dZ11qnζ, λ2dω,qζ+d1qnZ21qnζ, λ2dω,qζω0dZ1ζ, λ2dω,qζ+dZ2ζ, λ2dω,qζ n,
due to Z(·, λ) ∈ H. In the limit-point case, we find
𝓁λ, 1qnmλr1qn=2 Im λω0dθ1ζ, λ2dω,qζ+d1qnθ2ζ, λ2dω,qζ1,
where Im λ ≠ 0. As r1qn0 , Z1qnζ, λZζ, λn . Moreover, we have
ω0d|{(λ, 1qn)-m(λ)}θ(1)(ζ, λ)|2dω,qζ+d1qn|{𝓁(λ, 1qn)-m(λ)}θ(2)(ζ, λ)|2dω,qζ=|𝓁(λ, 1qn)-m(λ)|2(ω0d|θ(1)(ζ, λ)|2dω,qζ+d1qn|θ(2)(ζ, λ)|2dω,qζ)(4(Im λ)2[ω0d|θ(1)(ζ, λ)|2dω,qζ+d1qn|θ(2)(ζ, λ)|2dω,qζ])-1,
which implies that
ω0d|Z(1)1qn(ζ, λ)|2dω,qζ+d1qn|Z(2)1qn(ζ, λ)|2dω,qζω0d|Z(1)(ζ, λ)|2dω,qζ+d|Z(2)(ζ, λ)|2dω,qζ. 

Let fHn1qn>d, n . Define

(3.6)
G1qnζ, ς, λ=Z1qnζ, λθζ, λ, ςζ,θζ, λZ1qnς, λ, ς>ζ,R1qnfζ, λ=ω0dG1qnζ, ς, λf1ςdω,qς+d1qnG1qnζ, ς, λf2ςdω,qς, λ.

Without loss of generality, we can assume that λ = 0 is not an eigenvalue of the BVP (3.1)–(3.5). Now let us prove that the resolvent operator is compact.

Theorem 3.2.

G1qnζ, ςλ=0(1qn>d, n) defined as (3.6) is a ω, q-Hilbert–Schmidt kernel, i.e.,

ω0dω0d|G1qnζ, ς|2dω,qζdω,qς<+,d1qnd1qn|G1qnζ, ς|2dω,qζdω,qς<+.

Proof

By (3.6), it is obvious that

ω0ddω,qζω0d|G1qnζ, ς|2dω,qς<+,d1qndω,qζd1qn|G1qnζ, ς|2dω,qς<+,
due to Z1qn , (·, λ), θ (·, λ) ∈ Hn (1qn>d, n) . Hence
(3.7)
ω0dω0d|G1qnζ, ς|2dω,qζdω,qς<+,d1qnd1qn|G1qnζ, ς|2dω,qζdω,qς<+.

Theorem 3.3 ([20]).

Let A {ti} = {xi}, i ∈ ℕ, where

(3.8)
xi=k=1ηiktk, i, k.

If

(3.9)
i,k=1|ηik|2<+,
then the operator A is compact in l2.

Theorem 3.4.

Let 𝒯 be the ω, q-integral operator 𝒯: Hn → Hn ( 1qn>d,n ),

Tfζ=ω0dG1qnζ, ςf1ςdω,qς, ζω0, d ,d1qnG1qnζ, ςf2ςdω,qς, ζd, 1qn,
where
fζ=f1ζ , ζω0, d ,f2ζ , ζd, 1qn.

Then 𝒯 is a compact self-adjoint operator in space Hn.

Proof

Let

φi:=φiζ=φi1ζ , ζω0, d ,φi2ζ ,ζd, 1qn,(i, n, 1qn>d)
be a complete, orthonormal basis of Hn. Let i, k, n ∈ ℕ, 1qn>d . Write
ti=f, φin=ω0df(1)(ζ)φi(1)(ζ)¯dω,qζ+d1qnf(2)(ζ)φi(2)(ζ)¯dω,qζ,xi=g,φin=ω0dg(1)(ζ)φi(1)(ζ)¯dω,qζ+d1qng(2)(ζ)φi(2)(ζ)¯dω,qζ,ηik=ω0dω0dG1qn(ζ, ς)φi(1)(ζ)φk(1)(ς)¯dω,qζdω,qς+d1qnd1qnG1qn(ζ, ς)φi(2)(ζ)φk(2)(ς)¯dω,qζdω,qς.
Hn is mapped isometrically on to l2. By this mapping, 𝒯 transforms into the operator A defined by (3.8) in l2 and (3.7) is translated into (3.9). By Theorems 3.2 and 3.3, we see that A and 𝒯 are compact operators.

Let h, gHn and 1qn>d , n ∈ ℕ. Then we have

𝒯h,gn=dω0(𝒯h(1))(ζ)g(1)(ζ¯)dω,qζ+1qnd(𝒯h(2))(ζ)g(2)(ζ¯)dω,qζ=dω0dω0G1qn(ζ, ς)h(1)(ς)dω,qςg(1)(ζ)¯dω,qζ+1qnd1qndG1qn(ζ, ς)h(2)(ς)dω,qςg(2)(ζ)¯dω,qζ=dω0h(1)(ς)(dω0G1qn(ς, ζ)g(1)(ζ)¯dω,qζ)dω,qς+1qndh(2)(ς)(1qndG1qn(ς, ζ)g(2)(ζ)¯dω,qζ)dω,qς=h, Tgn,
since G1qnζ, γ is a symmetric function.

From Theorem 3.4, we conclude that 𝒯 has a discrete spectrum. Let λm,1qn and

θm,1qnζ:=θm,1qn1ζ, λm,1qn , ζω0, d ,θm,1qn2ζ, λm,1qn , ζd, 1qn,(m, n, 1qn>d)
be the eigenvalues and eigenfunctions of the BVP (3.1)–(3.5) and
αm,1qn2=ω0dθm,1qn12ζdω,qς+d1qnθm,1qn22ζdω,qς.

By Theorem 3.4 and the Hilbert–Schmidt theorem, we infer that

(3.10)
ω0d|f1ζ|2dω,qζ+d1qn|f2ζ|2dω,qζ=m=11αm,1qn2ω0df1ζϕm,1qn1ζdω,qζ+d1qnf2ζϕm,1qn2ζdω,qζ2.

Define

ϱ1qnλ=λ<λm,1n<01αm,1qn2 ,for λ0,0λm,1qn<λ1αm,1qn2,for λ>0.
Then, (3.10) can be written as
(3.11)
ω0df1ζ2dω,qζ+d1qnf2ζ2dω,qζ=|Fλ|2dϱ1qnλ,
where
Fλ=ω0df1ζϕm,1qn1ζdω,qζ+d1qnf2ζϕm,1qn2ζdω,qζ.

Lemma 3.5.

For any positive S, there is a positive number B = B (S) not depending on n so that

-SSϱ1qnλ=Sλm,1qn<S1αm,1qn2=ϱ1qnSϱ1qnS<B.

Proof

Let sin β ≠ 0. Since θ(ζ, λ) is continuous in domain −SλS, ω0, d(d, 1qn] , and the condition θ(1)(ω0, λ) = sin β, there exists a positive number h such that for |λ| < S,

(3.12)
1h2ω0ω0+hθ1ζ, λdω,qζ2>12sin2β.

Let

fhζ=1h,ω0ζω0+h,0,ζ>ω0+h.

From (3.12), we find

ω0ω0+hfh2ζdω,qζ=1h=1hω0ω0+hθ1ζ, λdω,qζ2dϱ1qnλSS1hω0ω0+hθ1ζ, λdω,qζ2dϱαλ>12sin2βϱ1qnSϱ1qnS.

If sin β = 0, then we define fh(ζ) as

fhζ=1h2, ω0ζω0+h,0,ζ>ω0+h. 

This proves the lemma.

Now, we will give an expansion into a Fourier series of resolvent. By ω, q-integration by parts, we obtain

dω0[1qD-ωq,1qDω,qy(1)(ζ, λ)-v(ζ)y(1)(ζ, λ)]θm,1qn(1)(ζ)dω,qζ+1qnd[1qD-ωq,1qDω,qy(2)(ζ, λ)-v(ζ)y(2)(ζ, λ)]θm,1qn(2)(ζ)dω,qζ=dω0[1qD-ωq,1qDω,qϕm,1qn(1)(ζ)-v(ζ)θm,1qn(1)(ζ)]y(1)(ζ, λ)dω,qζ+1qnd[1qD-ωq,1qDω,qϕm,1qn(2)(ζ)-v(ζ)θm,1qn(2)(ζ)]y(2)(ζ, λ)dω,qζ=-λm,1qndω0y(1)(ζ, λ)θm,1qn(1)(ζ)dω,qζ-λm,1qn1qndy(2)(ζ, λ)θm,1qn(2)(ζ)dω,qζ=-λm,1qnφm(λ),
where m ∈ ℕ. Let
yζ, λ=m=1φmλψm,1qnζ,am=ω0dfζψm,1qn1ζdω,qζ+d1qnfζψm,1qn2ζdω,qζ,
where m ∈ ℕ. Since y(ζ, λ) satisfies the equation
1qDωq,1qDω,qyζ, λ+vζλyζ, λ=fζ,
we find
am=ω0d1qDωq,1qDω,qy1ζ, λ+vζλy1ζ, λθm,1qn1ζdω,qζ+d1qn1qDωq,1qDω,qy2ζ, λ+vζλy2ζ, λθm,1qn2ζdω,qζ=λm,1qnφmλλφmλ,   m,n,1qn>d.

Thus, we get

φmλ=amλm,1qnλ(m, n, 1qn>d),
and
yζ, λ=G1qnζ, , λ, f¯n=m=1amθm,1qnζλm,1qnλ.

Hence

(3.13)
R1qnfζ, z=m=1θm,1qnζαm,1qn2λm,1qnzf, θm,1qnn=θζ,λλzf, θm,1qnndϱ1qnλ.

Lemma 3.6.

For each nonreal z and fixed ζ, the following relation holds

(3.14)
θζ,λzλ2dϱ1qnλ<S.

Proof

Writing

fς=θm,1qnςαm,1qn
yields
(3.15)
1αm,1qnG1qnζ, , λ, θm,1qnn=θm,1qnζαm,1qnλm,1qnz,
due to the eigenfunctions θm,1qn(ζ) are orthogonal. Combining (3.15) and (3.10), we see that
ω0dG1qnζ, ς, z2dω,qς+d1qnG1qnζ, ς, z2dω,qς=m=1θm,1qnζ2αm,1qn2λm,1qnz2=θζ,λλz2dϱ1qnλ.

By Lemma 3.1, the integral on the left converges and the result is immediate.

It follows from Lemma 8 that the set ϱ1qnλ is bounded. Using Helly’s theorems ([15]), one can find a sequence {1/qnk} such that ϱ1qnkλ converges to a monotone function ϱ(λ) (as nk → ∞).

Lemma 3.7.

Let z be a nonreal number and ζ be a fixed number. Then we have

(3.16)
θζ,λzλ2dϱλS.

Proof

For arbitrary η > 0, it follows from (3.14) that

ηηϕς,λzλ2dϱ1qnλ<S.

Letting η → ∞ and n → ∞, we get the desired result.

Lemma 3.8.

For arbitrary η > 0, we have

ηdϱλ|zλ|2<, ηdϱλ|zλ|2<.

Proof

Let sin β ≠ 0. Writing ζ = 0 in (3.16), we obtain

dϱλ|zλ|2<.

Let sin β = 0. Then

1αm,1qnDq,ζG1qnζ, , z, θm,1qnn=Dq,ζθm,1qnζαm,1qnλm,1qnz.

By (3.11), we find

ω0dDq,ζG1qnζ, ς, z2dω,qζ+d1qnDq,ζG1qnζ, ς, z2dω,qζ=Dq,ζθζ,λzλ2dϱ1qnλ.

Lemma 3.9.

Let

Rfζ, z=ω0Gζ, ς, zfςdω,qς,
where fH, and
Gζ, ς, z=Zζ, zθς, z, ςζ, ζd, ςd,θζ, zZς, z , ς>ζ, ζd, ςd.

Then, we have

ω0d|Rfζ, z|2dω,qζ+d|Rfζ, z|2dω,qζ1v2ω0df1ζ2dω,qζ+df2ζ2dω,qζ,
where v = Im z.

Proof

Combining (3.13) and (3.10), for each 1qn>d , n ∈ ℕ, we obtain

ω0dR1qnfζ, z2dω,qζ+d1qnR1qnfζ, z2dω,qζ=m=1f,θm,1qn,zn2αm,1qn2|λm,1qnz|2=1v2ω0df1ς2dω,qς+1v2d1qnf2ς2dω,qς.

Letting n → ∞, we get the desired result.

Theorem 3.10 (Integral Representation of the Resolvent).

For every non-real z and for each fH, we obtain

Rfζ, z=θζ,λλzFλdϱλ,
where
Fλ=ω0df1ζθ1ζ, λdω,qζ+limσdσf2ζθ2ζ, λdω,qζ.

Proof

Suppose that f(ζ) = fσ(ζ) satisfies (3.2)–(3.4) and vanishes outside the set [ω0, d) ∪ (d, σ], where d<σ<1qn , n ∈ ℕ. Let

Fσλ=ω0dfσ1ζθ1ζ, λdω,qζ+dσfσ2ζθ2ζ, λdω,qζ.

By (3.13), we see that

(3.17)
R1qnfσζ, z=θζ,λλzFσλdϱ1qnλ=aθζ,λλzFσλdϱ1qnλ+aaθζ,λλzFσλdϱ1qnλ+aθζ,λλzFσλdϱ1qnλ=I1+I2+I3.

Firstly, we will estimate I1. From (3.13), we deduce that

I1=--aθ(ζ,λ)z-λFσ(λ)dϱ1qn(λ)=λk,1qn<-aθk,1qn(ζ)αk,1qn2(z-λk,1qn){ω0dfσ(1)(ζ)θk,qn(1)(ζ)dω,qζ+dσfσ(2)(ζ)θk,qn(2)(ζ)dω,qζ}(λk,1qn<-aθk,1qn2(ζ)αk,1qn2|z-λk,1qn|2)1/2×(λk,1qn<-a1αk,1qn2|ω0dfσ(1)(ζ)θk,qn(1)(ζ)dω,qζ+dσfσ(2)(ζ)θk,qn(2)(ζ)dω,qζ|2)1/2.

Integrating twice by parts, we find

ω0dfσ(1)(ζ)θk,qn(1)(ζ)dω,qζ+σdfσ(2)(ζ)θk,qn(2)(ζ)dω,qζ=-1λk,1qnω0dfσ(1)(ζ){1qD-ωq,1qDω,qθk,qn(1)(ζ)-v(ζ)θk,qn(1)(ζ)}dω,qζ-1λk,1qndσfσ(2)(ζ){1qD-ωq,1qDω,qθk,qn(2)(ζ)-v(ζ)θk,qn(2)(ζ)}dω,qζ=-1λk,1qnω0d{1qD-ωq,1qDω,qfσ(1)(ζ)-v(ζ)fσ(1)(ζ)}θk,qn(1)(ζ)dω,qζ-1λk,1qndσ{1qD-ωq,1qDω,qfσ(2)(ζ)-v(ζ)fσ(2)(ζ)}θk,qn(2)(ζ)dω,qζ.

By Lemma 3.6, we get

I1K1/2a×λk,1qn<a1αk,1qn2ω0d{1qDωq,1qDω,qfσ1ζvζfσ1ζ}θk,1qn(1)ζdω,qζ+dσ1qDωq,1qDω,qfσ2ζvζfσ2ζθk,1qn(2)ζdω,qζ21/2.
Using Bessel inequality, we see that
I1K1/2aω0σ1qDωq,1qDω,qfσ1ζvζfσ1ζ2dω,qζ+dσ1qDωq,1qDω,qfσ2ζvζfσ2ζ2dω,qζ1/2=Ca.

It is proved similarly that I3Ca . Then I1 and I3 tend to zero as a → ∞, uniformly in 1qn . It follows from the Helly selection theorem and (3.17) that

(3.18)
Rfσζ, z=θζ,λzλFσλdϱλ.
As is known, if f(·) ∈ H, then we find a sequence fσςσ=1 that satisfies the previous conditions and tends to f(ζ) as σ → ∞. From (3.10), the sequence of Fourier transform converges to the transform of f(ζ). Using Lemmas 3.7 and 3.9, we can pass to the limit σ → ∞ in (3.18). Thus, we get the desired result.

Remark 3.11.

Using Theorem 3.10, we infer that

ω0(Rf1)ς, zg1ςdω,qς+d(Rf2)ς, zg2ςdω,qς=FλGλzλdϱλ,
where
Fλ=ω0df1ζθ1ζ, λdω,qζ+limσdσf2ζθ2ζ, λdω,qζ,
and
Gλ=ω0dg1ζθ1ζ, λdω,qζ+limσω0σg2ζθ2ζ, λdω,qζ.

Notes

[1] Statements and Declarations.

This work does not have any conflict of interested result.

[2] Availability of data and materials.

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

DOI: https://doi.org/10.2478/amsil-2024-0001 | Journal eISSN: 2391-4238 (formerly 0860-2107) | Journal ISSN: 0860-2107
Language: English
Page range: 23 - 41
Submitted on: Oct 5, 2023
Accepted on: Jan 3, 2024
Published on: Jan 18, 2024
Published by: University of Silesia in Katowice, Institute of Mathematics
In partnership with: Paradigm Publishing Services

© 2024 Bilender P. Allahverdiev, Hüseyin Tuna, Hamlet A. Isayev, published by University of Silesia in Katowice, Institute of Mathematics
This work is licensed under the Creative Commons Attribution 4.0 License.