Open Access

On Basis Property of Root Functions For a Class Second Order Differential Operator


Cite

Introduction

We consider the boundary value problem for the differential equation on the interval [−1, 1] 𝓁(u):u''+q(x)u=λw(x)u,\ell (u): \equiv - {u^{''}} + q(x)u = \lambda w(x)u, with the boundary conditions L1(u):=u(1)+hu'(1)=0,{L_1}(u): = u( - 1) + h{u^'}( - 1) = 0,L2(u):=(β1u(1)β2u'(1))+λ(β˜1u(1)β˜2u'(1))=0,{L_2}(u): = ({\beta _1}u(1) - {\beta _2}{u^'}(1)) + \lambda \left( {{{\widetilde \beta }_1}u(1) - {{\widetilde \beta }_2}{u^'}(1)} \right) = 0, and the transmission conditions L3(y):=γ1u(0)δ1u(+0)=0,{L_3}(y): = {\gamma _1}u( - 0) - {\delta _1}u( + 0) = 0,L4(y):=γ2u'(0)δ2u'(+0)=0,{L_4}(y): = {\gamma _2}{u^'}( - 0) - {\delta _2}{u^'}( + 0) = 0, where w(x)={w12,1x<0,w22,0<x1,w(x) = \left\{ {\begin{array}{*{20}{c}} {w_1^2,}&{ - 1 \leqslant x < 0,} \\ {w_2^2,}&{0 < x \leqslant 1,} \end{array}} \right.w1w2, the real valued function q(x) ∈ C ([−1, 0) ∪ (0, 1] ) and has finite limits q(±0)=limx±0q(x),λq( \pm 0) = \mathop {\lim }\limits_{x \to \pm 0} q(x),\lambda , λ is complex parameter and we assume that h, wi, βi, β˜i{\widetilde \beta _i}, γi, δi (i = 1, 2 ) are real numbers, γi, δi are positive coefficients, |β1| + |β2| ≠ 0, |β˜1|+|β˜2|0\left| {{{\widetilde \beta }_1}} \right| + \left| {{{\widetilde \beta }_2}} \right| \ne 0 and ρ:=β1β˜2β˜1β2\rho : = {\beta _1}{\widetilde \beta _2} - {\widetilde \beta _1}{\beta _2}.

Spectral problems for the discontinuous Sturm-Liouville equations with eigenparameter dependent boundary conditions have been growing interest with physical applications and examined in [1, 2, 3, 4]. For example, boundary value problems with the form (1)(7) is encountered in vibrating string problems investigated in [19]. This happens whenever one applies the method of seperation of variables to solve the corresponding partial differential equation which contain a directional derivative in boundary conditions. On eigenvalues problems for second order equation with spectral parameter in the boundary conditions are considered in [5, 6, 7, 8, 9, 10, 11, 12, 13, 20, 21, 22]. Some self adjoint problems on eigenvalues for second order equation with spectral parameter in the boundary conditions are considered in [5, 6, 7, 8, 9, 10, 11]. The corresponding problems led to the eigenvalue problem for a linear operator acting on the space L2 ⊕ ℂN, where ℂN is N− dimensional Euclidean space of complex numbers. In [4] for distinct cases, it is shown that the eigenfunctions of the spectral problem formed a defect basis in L2 (0, 1). In [14] Rayleigh-Ritz formula is developed for eigenvalues.

The goal of this work is to investigate the problem of completeness, minimality and basis property of the eigenfunctions of the boundary value problem (1)(5). In this study, we introduce a special inner product in a special Hilbert space and construct a linear operator A in it so that the problem (1)(5) can be interpreted as the eigenvalue problem for A. Let us give the main definitions and theorem which will be used in our main results:

Definition 1

[23] A sequence {fj}j≥1 of vectors of a Hilbert space H is called a basis of this space if every vector fH can be expanded in a unique way in a series f=j=1cjfjf = \mathop {\mathop {\sum }\limits^\infty }\limits_{j = 1} {c_j}{f_j}, which converges in the norm of the space H.

Definition 2

[23] A basis {fj}j≥1 of H is called a Riesz basis if it is obtained from an orthonormal basis by means of a bounded linear invertible operator.

Theorem 1

[24] Let A be selfadjoint compact operator or selfadjoint operator which has discrete spectrums. Then the eigenfunctions of the operator A form orthonormal basis in the Hilbert space H.

For example; u''(x)+λu(x)=0,x(0,1)u(0)=0,u'(1)=dλu(1),d>0,\begin{gathered} {u^{''}}(x) + \lambda u(x) = 0,\;\;\;x \in \left( {0,1} \right) \hfill \\ u(0) = 0,{u^'}(1) = d\lambda u(1),\;d > 0, \hfill \\ \end{gathered} this boundary value problem has only the eigenfunctions un(x)=2sinλnx{u_n}(x) = \sqrt 2 \sin \sqrt {{\lambda _n}} x, n = 0, 1, 2.., with positive eigenvalues found from the equation cotλ=dλ\cot \sqrt \lambda = d\sqrt \lambda . By deleting an arbitrary eigenfunction, one obtains a basis in the space Lp(0, 1) p > 1, even a Riesz bases in the case p = 2. [4]

The Operator Formulation of the Problem and Main Results

It is convenient to represent the spectral problem (1)(5) as an eigenvalue problem for a linear problem in a Hilbert space. We denote by H = L2[−1, 1] ⊕ ℂ the special Hilbert space of all elements u˜=(u(x)u1),v˜=(v(x)v1)H,\widetilde u = \left( {\begin{array}{*{20}{c}} {u(x)} \\ {{u_1}} \end{array}} \right),\widetilde v = \left( {\begin{array}{*{20}{c}} {v(x)} \\ {{v_1}} \end{array}} \right) \in H, with the inner product (u˜,v˜)=w12γ1γ210u(x)v(x)¯dx+w22δ1δ201u(x)v(x)¯dx+δ1δ2ρu1v1¯.\left( {\widetilde u,\widetilde v} \right) = w_1^2{\gamma _1}{\gamma _2}\int\limits_{ - 1}^0 u(x)\overline {v(x)} dx + w_2^2{\delta _1}{\delta _2}\int\limits_0^1 u(x)\overline {v(x)} dx + \frac{{{\delta _1}{\delta _2}}}{\rho }{u_1}\overline {{v_1}} . In the space we define the operator Au˜=(1w(x)[u''+q(x)u]β1u(1)β2u'(1)),A\widetilde u = \left( {\begin{array}{*{20}{c}} {\frac{1}{{w(x)}}[ - {u^{''}} + q(x)u]} \\ {{\beta _1}u(1) - {\beta _2}{u^'}(1)} \end{array}} \right), on the domain D(A)={u˜|u˜=(u(x),u1)H,u(x),u'(x)AC([1,0)(0,1]),u'(±0)=limx±0u'(x),𝓁(u)L2[1,1],L1u=L3u=L4u=0,u1=β˜1u(1)β˜2u'(1)),}D(A) = \left\{ {\begin{array}{*{20}{c}} {\widetilde u|\widetilde u = \left( {u(x),{u_1}} \right) \in H,u(x),{u^'}(x) \in AC\left( {[ - 1,0) \cup \left( {0,1} \right]} \right),} \\ {{u^'}( \pm 0) = \mathop {\lim }\limits_{x \to \pm 0} {u^'}(x),\ell (u) \in {L_2}\left[ { - 1,1} \right],} \\ {{L_1}u = {L_3}u = {L_4}u = 0,{u_1} = {{\widetilde \beta }_1}u(1) - {{\widetilde \beta }_2}{u^'}(1)),} \end{array}} \right\} where AC ([−1, 1]) is the space of all absolutely continuous functions on the interval [−1, 1]. Obviously, the operator A is well defined in H. It is clear that the spectral problem (1)(5) is equivalent to operator equation Au˜=λu˜,A\widetilde u = \lambda \widetilde u, and the eigenvalues of A coincide with the eigenvalues of the problem (1)(5) (see Lemma 1.4. in [5]). Also, there exists a correspondence between eigenfunctions u˜k(x)(uk(x)β˜1u(1)β˜2u'(1))){\widetilde u_k}(x) \leftrightarrow \left( {\begin{array}{*{20}{c}} {{u_k}(x)} \\ {{{\widetilde \beta }_1}u(1) - {{\widetilde \beta }_2}{u^'}(1))} \end{array}} \right)

Lemma 2

The domain D(A) of the operator A is dense in the space H.

Proof

The proof is similar using the same method in [15]. Suppose that f˜H\widetilde f \in H is orthogonal to all g˜D(A)\widetilde g \in D(A) with respect to the inner product (6), where f˜=(f(x),f1)\widetilde f = (f(x),{f_1}), g˜=(g(x),g1)\widetilde g = (g(x),{g_1}). Let C˜0\widetilde C_0^\infty denote the set of functions Φ(x)={φ1(x),x[1,0),φ2(x),x(0,1],\;\;\;\Phi (x) = \left\{ {\begin{array}{*{20}{c}} {{\varphi _1}(x),x \in \left[ { - 1,0} \right),} \\ {{\varphi _2}(x),x \in \left( {0,1} \right],} \end{array}} \right. where φ1(x)C0[1,0){\varphi _1}(x) \in C_0^\infty \left[ { - 1,0} \right) and φ2(x)C0(0,1]{\varphi _2}(x) \in C_0^\infty \left( {0,1} \right]. Since C˜00D(A)\widetilde C_0^\infty \oplus 0 \subset D(A) (0 ∈ ℂ), any u˜=(u(x),0)C˜00\widetilde u = \left( {u(x),0} \right) \in \widetilde C_0^\infty \oplus 0 is orthogonal to f˜\widetilde f , namely, (f˜,u˜)=w12γ1γ210f(x)u(x)¯dx+w22δ1δ201f(x)u(x)¯dx=(f,u)1(\widetilde f,\widetilde u) = w_1^2{\gamma _1}{\gamma _2}\int\limits_{ - 1}^0 f(x)\overline {u(x)} dx + w_2^2{\delta _1}{\delta _2}\int\limits_0^1 f(x)\overline {u(x)} dx = (f,u{)_1} where (,)1 denotes inner product in L2[−1, 1]. This implies that f(x) is orthogonal to C˜0\widetilde C_0^\infty and (f, u)1 = 0. Hence, (f˜,g˜)=δ1δ2ρf1g1¯=0.(\widetilde f,\widetilde g) = \frac{{{\delta _1}{\delta _2}}}{\rho }{f_1}\overline {{g_1}} = 0. Thus f1 = 0 since g1=β˜1g(1)β˜2g'(1){g_1} = {\widetilde \beta _1}g(1) - {\widetilde \beta _2}{g^'}(1) can be chosen arbitrary. So f˜=(0,0).\widetilde f = (0,0).. Therefore, D(A) is dense in H.

Lemma 3

The operator A is selfadjoint.

Proof

In this case integrating by parts, we obtain that (Af˜,g˜)(A\widetilde f,\widetilde g) is real. Taking into account Lemma 2, we find that the operator A is symmetric in the space H. Boundary value problem (1)(5) is solvable for every non-eigenvalue λ and has discrete spectrum. Therefore the operator A is symmetric and has discrete spectrum. Hence, the operator A is selfadjoint in H.

Theorem 4

The eigenfunctions of the operator A form an orthonormal basis in the space H = L2[−1, 1] ⊕ ℂ.

Proof

The operator A has countable many eigenvalues {λn}n=1\left\{ {{\lambda _n}} \right\}_{n = 1}^\infty which have the asymptotic form [16]: λn=1w1+w2π(n1)+O(1n),n.{\lambda _n} = \frac{1}{{{w_1} + {w_2}}}\pi (n - 1) + O(\frac{1}{n}),n \to \infty .

Then for any number λ which is not an eigenvalue and arbitary f˜H\widetilde f \in H it can be found an element ũD(A) satisfying the condition (AλI)u˜=f˜.\left( {A - \lambda I} \right)\widetilde u = \widetilde f.. Thus the operator (AλI) is invertible except for the isolated eigenvalues. Without loss of generality we assume that the point λ = 0 is not an eigenvalue. Then we obtain that the bounded inverse operator A−1 is defined in H. Thus, the selfadjoint operator A−1 has at most countable many eigenvalues and each one of them converges to zero at the infinity. So, the selfadjoint operator A−1 is compact. Applying the Hilbert-Schmidt theorem to this operator we obtain that the eigenfunctions of the operator A form an orthonormal basis in H. Theorem is proved.

Now we consider the case ρ < 0. We assume that the operator A is defined by formula (7) on the domain D(A). In the space H = L2 ⊕ ℂ for ũ, H the scalar product is defined by formula (u˜,v˜)=w12γ1γ210u(x)v(x)¯dx+w22δ1δ201u(x)v(x)¯dxδ1δ2ρu1v1¯.\left( {\widetilde u,\widetilde v} \right) = w_1^2{\gamma _1}{\gamma _2}\int\limits_{ - 1}^0 u(x)\overline {v(x)} dx + w_2^2{\delta _1}{\delta _2}\int\limits_0^1 u(x)\overline {v(x)} dx - \frac{{{\delta _1}{\delta _2}}}{\rho }{u_1}\overline {{v_1}} . In this case the operator A is not selfadjoint in the space H. Therefore we introduce the operator J is defined by J=(I00I)J = \left( {\begin{array}{*{20}{c}} I&0 \\ 0&{ - I} \end{array}} \right) where I is the identity operator in H. Operator J is selfadjoint and invertible.

In this case, the boundary value problem (1)(5) is equivalent to eigenvalue problem for the operator pencil (BλJ)u˜=0,(B - \lambda J)\widetilde u = 0, in the space H such that B = JA. We obtain that (8) is equivalent to (10).

Lemma 5

The operator A is J− selfadjoint in the Hilbert space H.

Proof

Analogously to Lemma 2, we can show that the domain D(A) is dense in space H. From (9) and (10) applying two times integration by parts, (, ũ) is real. Hence, the operator B is symmetric. Therefore, the operator A is J− symmetric in the space H. In this case it can be proved that the operator A has a discrete spectrum. Taking into consideration that the operator B is symmetric we have that the operator JA is selfadjoint.

Corollary 6

From the system{un}0\left\{ {{u_n}} \right\}_0^\infty one can eliminate one element so that the remaining elements will form a complete and minimal system in the space L2[−1, 1].

Proof

By Theorem 4, the system of eigenfunctions u˜n(x)={un(x)u1},{\widetilde u_n}(x) = \left\{ {\begin{array}{*{20}{c}} {{u_n}(x)} \\ {{u_1}} \end{array}} \right\}, (u1 ∈ ℂ ) of the operator A forms an orthonormal basis in H. Hence, the system of the eigenfunctions {u˜n(x)}1\left\{ {{{\widetilde u}_n}(x)} \right\}_1^\infty is complete and minimal in the space H. Thus, of course, codimP = 1, then by Lemma 2.1 in [5], the system {n(x)} = {un(x)} whose one element is omitted from forms a complete and minimal system in P(H) = L2[−1, 1]. Hence, the eigenfunctions {un(x)}0\left\{ {{u_n}(x)} \right\}_0^\infty (nn0, n0 is an arbitrary nonnegative integer) of the boundary problem (1)(5) are complete and minimal system in L2[−1, 1].

Theorem 7

The eigenfunctions of the operator A form a Riesz basis in the Hilbert space H.

Proof

From the operator J is a bounded operator, using the idea in Theorem 4, it can be shown that the operator B = JA is invertible since the selfadjoint operator B−1 is compact. Since the selfadjoint operator B−1 has at most countable many eigenvalues which converge to zero at infinity. Hence the operator B−1 is compact. Then applying Azizov-Iokhvidov theorem in section IV of [18] to operator B we obtain that the eigenfunctions of the J− selfadjoint operator A form a Riesz basis in the space H = L2 ⊕ ℂ.

Let us consider another boundary value problem for the Sturm-Liouville equation 𝓁u:=p(x)u''+q(x)u=λu,\ell u: = - p(x){u^{''}} + q(x)u = \lambda u, on [−1, 0) ∪ (0, 1] with boundary condition, L1(u):=α1u(1)+α2u'(1)=0{L_1}(u): = {\alpha _1}u( - 1) + {\alpha _2}{u^'}( - 1) = 0L2(λ)u:=λ(β˜1u(1)β˜2u'(1))+(β1u(1)β2u'(1))=0,{L_2}(\lambda )u: = \lambda ({\widetilde \beta _1}u(1) - {\widetilde \beta _2}{u^'}(1)) + ({\beta _1}u(1) - {\beta _2}{u^'}(1)) = 0, with transmission conditions at the point of discontinuity x = 0, L3(u):=γ1u(0)δ1u(+0)=0,{L_3}(u): = {\gamma _1}u( - 0) - {\delta _1}u( + 0) = 0,L4(u):=γ2u'(0)δ2u'(+0)=0,{L_4}(u): = {\gamma _2}{u^'}( - 0) - {\delta _2}{u^'}( + 0) = 0, here p(x)={p12,1x<0,p22,0<x1,p(x) = \left\{ {\begin{array}{*{20}{c}} {p_1^2,}&{ - 1 \leqslant x < 0,} \\ {p_2^2,}&{0 < x \leqslant 1,} \end{array}} \right.p1p2, the real valued function q(x) is continuous in [−1, 0) ∪ (0, 1] and has finite limits q(±0)=limx±0q(x)q( \pm 0) = \mathop {\lim }\limits_{x \to \pm 0} q(x), λ is complex parameter, we assume that h, pi, βi, β˜i{\widetilde \beta _i}, γi, δi (i = 1, 2 ) are real numbers, γi, δi are positive coefficients, |β1|+ |β2| ≠ 0, |β˜1|+|β˜2|0\left| {{{\widetilde \beta }_1}} \right| + \left| {{{\widetilde \beta }_2}} \right| \ne 0 and ρ:=β1β˜2β˜1β2>0\rho : = {\beta _1}{\widetilde \beta _2} - {\widetilde \beta _1}{\beta _2} > 0.

It is convenient to represent the spectral problem (11)(15) as an eigenvalue problem for a linear problem in a special Hilbert space. We denote by H = L2 (−1, 1) ⊕ ℂ the special Hilbert space of all elements f=(f(x)f1),g=(g(x)g1)H,f = \left( {\begin{array}{*{20}{c}} {f(x)} \\ {{f_1}} \end{array}} \right),g = \left( {\begin{array}{*{20}{c}} {g(x)} \\ {{g_1}} \end{array}} \right) \in H, with the inner product (f,g)=1p12γ1γ210f(x)g(x)¯dx+1p22δ1δ201f(x)g(x)¯dx+δ1δ2ρf1g1¯\left( {f,g} \right) = \frac{1}{{p_1^2}}{\gamma _1}{\gamma _2}\int\limits_{ - 1}^0 f(x)\overline {g(x)} dx + \frac{1}{{p_2^2}}{\delta _1}{\delta _2}\int\limits_0^1 f(x)\overline {g(x)} dx + \frac{{{\delta _1}{\delta _2}}}{\rho }{f_1}\overline {{g_1}} In the space we define the operator Lu=(p(x)u''+q(x)uβ1u(1)β2u'(1)),Lu = \left( {\begin{array}{*{20}{c}} { - p(x){u^{''}} + q(x)u} \\ {{\beta _1}u(1) - {\beta _2}{u^'}(1)} \end{array}} \right), on the domain D(L)={u˜|u˜=(u(x),u1)H,u(x),u'(x)AC([1,0)(0,1]),u'(±0)=limx±0u'(x),𝓁(u)L2[1,1],L1u=L3u=L4u=0,u1=β˜1u(1)β˜2u'(1)),}D(L) = \left\{ {\begin{array}{*{20}{c}} {\widetilde u|\widetilde u = \left( {u(x),{u_1}} \right) \in H,u(x),{u^'}(x) \in AC\left( {[ - 1,0) \cup \left( {0,1} \right]} \right),} \\ {{u^'}( \pm 0) = \mathop {\lim }\limits_{x \to \pm 0} {u^'}(x),\ell (u) \in {L_2}\left[ { - 1,1} \right],} \\ {{L_1}u = {L_3}u = {L_4}u = 0,{u_1} = {{\widetilde \beta }_1}u(1) - {{\widetilde \beta }_2}{u^'}(1)),} \end{array}} \right\}u(±0) has finite limits. Now we can rewrite the problem (11)(15) in the operator form as Lu=λu.Lu = \lambda u. The eigenvalues of the operator L coincide with the eigenvalues of the spectral problem (11)(15).

Theorem 8

The eigenfunctions of the operator L form an orthonormal basis in the space H = L2[−1, 1] ⊕ ℂ.

Proof

The operator L has countable many eigenvalues {λn}n=1\left\{ {{\lambda _n}} \right\}_{n = 1}^\infty which have the asymptotic form [17]: λn=α1α2α1+α2π(n1)+O(1n),n.{\lambda _n} = \frac{{{\alpha _1}{\alpha _2}}}{{{\alpha _1} + {\alpha _2}}}\pi (n - 1) + O(\frac{1}{n}),n \to \infty . Then using the idea of the proof of Theorem 4, it can be proved analogusly.

Now for the boundary value problem (11)(15) we consider the case ρ < 0. In the space H = L2 ∈ ℂ for f, gH the scalar product is defined by formula (f,g)=1p12γ1γ210f(x)g(x)¯dx+1p22δ1δ201f(x)g(x)¯dxδ1δ2ρf1g1¯.\left( {f,g} \right) = \frac{1}{{p_1^2}}{\gamma _1}{\gamma _2}\int\limits_{ - 1}^0 f(x)\overline {g(x)} dx + \frac{1}{{p_2^2}}{\delta _1}{\delta _2}\int\limits_0^1 f(x)\overline {g(x)} dx - \frac{{{\delta _1}{\delta _2}}}{\rho }{f_1}\overline {{g_1}} . Since ρ < 0, the operator A is not selfadjoint in the space H. Therefore we introduce the operator J is defined by J=(I00I)J = \left( {\begin{array}{*{20}{c}} I&0 \\ 0&{ - I} \end{array}} \right) where I is the identity operator in H. Operator J is selfadjoint and invertible in H.

In this case, the boundary value problem (11)(15) is equivalent to eigenvalue problem for the operator pencil (BλJ)u˜=0,(B - \lambda J)\widetilde u = 0, where B = JL is symmetric and L is J− symmetric in the space H.. Similary to Lemma 5, The operator L is J− selfadjoint in the Hilbert space H and we obtain the following results:

Theorem 9

The eigenfunctions of the operator L form a Riesz basis in the Hilbert space H.

Corollary 10

From the system{un}0\left\{ {{u_n}} \right\}_0^\infty one can eliminate one element so that the remaining elements will form a complete and minimal system in the space L2[−1, 1].

eISSN:
2444-8656
Language:
English
Publication timeframe:
Volume Open
Journal Subjects:
Life Sciences, other, Mathematics, Applied Mathematics, General Mathematics, Physics