The purpose of this paper is to investigate some spectral properties of Sturm-Liouville type problems with interior singularities. Some of the mathematical aspects necessary for developing our own technique are presented. By applying this technique we construct some special solutions of the homogeneous equation and present a formula and the existence conditions of Green's function. Furthermore, based on these results and introducing operator treatment in adequate Hilbert space, we derive the resolvent operator and prove self-adjointness of the considered problem.

1. Introduction

For inhomogeneous linear systems, the basic superposition principle says that the response to a combination of external forces is the self-same combination of responses to the individual forces. In a finite-dimensional system, any forcing function can be decomposed into a linear combination of unit impulse forces, each applied to a single component of the system, and so the full solution can be written as a linear combination of the solutions to the impulse problems. This simple idea will be adapted to boundary value problems governed by differential equations, where the response of the system to a concentrated impulse force is known as Green's function. With Green's function in hand, the solution to the inhomogeneous system with a general forcing function can be reconstructed by superimposing the effects of suitably scaled impulses. Green's function method provides a powerful tool to solve linear problems consisting of a differential equation (partial or ordinary, with, possibly, an inhomogeneous term) and enough initial and/or boundary conditions (also possibly inhomogeneous) so that this problem has a unique solution. The history of Green's function dates back to 1828, when Green [1] published work in which he sought solutions of Poisson’s equation for the electric potential defined inside a bounded volume with specified boundary conditions on the surface of the volume. He introduced a function now identified as what Riemann later coined Green's function. In 1877, Neumann [2] embraced the concept of Green's function in his study of Laplace's equation, particularly in the plane. He found that the two-dimensional equivalent of Green's function was not described by singularity of the form as in the three-dimensional case but by a singularity of the form . With the function’s success in solving Laplace's equation, other equations began to be solved using Green's function. The heat equation and Green's function have a long association with each other. After discussing heat conduction in free space, the classic solutions of the heat equation in rectangular, cylindrical, and spherical coordinates are offered. In the case of the heat equation, Hobson [3] derived the free-space Green's function for one, two and three dimensions, and the French mathematician Appell [4] recognized that there was a formula similar to Green's for the one-dimensional heat equation. Green's function is particularly well suited for wave problems with the detailed analysis of electromagnetic waves in surface wave guides and water waves. The leading figure in the development of Green's function for the wave equation was Kirchhoff [5], who used it during his study of the three-dimensional wave. Starting with Green's second formula, he was able to show that the three-dimensional Green's function is where .

The application of Green's function to ordinary differential equations involving boundary-value problems began with the work of Burkhardt [6]. Determination of Green's function is also possible using Sturm-Liouville theory. This leads to series representation of Green's function. Sturm-Liouville problems which contained spectral parameter in boundary conditions form an important part of the spectral theory of boundary value problems. This type of problems has a lot of applications in mechanics and physics (see [79] and references cited therein). In the recent years, there has been increasing interest in this kind of problems which also may have discontinuities in the solution or its derivative at interior points (see [1018]). In this study, we will investigate some spectral properties of the Sturm-Liouville differential equation on two intervals: on , with eigenparameter-dependent boundary conditions at the end points and . One has, and the transmission conditions at the singular interior point where the potential is real continuous function in each of the intervals and and has finite limits , is a complex spectral parameter, ,   and ,  and   and are real numbers.

Our problem differs from the usual regular Sturm-Liouville problem in the sense that the eigenvalue parameter is contained in both differential equation and boundary conditions, and two supplementary transmission conditions at one interior point are added to boundary conditions. Such problems are connected with discontinuous material properties, such as heat and mass transfer, vibrating string problems when the string loaded additionally with points masses, diffraction problems [8, 9], and varied assortment of physical transfer problems. We develop our own technique for the investigation of some spectral properties of this problem. In particular, we construct the Green’s function and adequate Hilbert space for self-adjoint realization of the considered problem.

2. Some Basic Solutions and Green's Function

Denote the determinant of the th and th columns of the matrix by . For self-adjoint realization in adequate Hilbert space, everywhere below we will assume that With a view to construct the Green's function we will define two special solutions of (2) by our own technique as follows. At first, consider the next initial-value problem on the left interval It is known that this problem has an unique solution which is an entire function of for each fixed (see, e.g., [19]). By applying the similar method of [13], we can prove that (2) on the right interval has an unique solution satisfying the equalities which is also an entire function of the parameter for each fixed . Consequently, the solution defined by satisfies (2) on whole , the first boundary condition of (3), and both transmission conditions (5).

By the same technique, we can define the solution by so that Consequently, satisfies (2) on whole , the second boundary condition (4), and both transmission condition (5). By using (9), (10), (14), and (15) and the well-known fact that the Wronskians and are independent of variable , it is easy to show that . We will introduce the characteristic function for the problems (2)–(5) as Similar to [13], we can prove that there are infinitely many eigenvalues of the BVTP (2)–(5) which coincide with the zeros of characteristic function .

Now, let us consider the nonhomogenous differential equation on together with the same boundary and transmission conditions (2)–(5), when . We will search the solution of this problem in the form (see, for example, [13]): where are arbitrary constants. Putting in (3)–(5) we have ,   Now, by substituting these equalities in (18), the following formula is obtained for the solution of (17) under boundary and transmission conditions (3)–(5): From this formula, we find that the Green’s function of the problem (2)–(5) has the form: and the solution (20) can be rewritten in the terms of this Green’s function as

3. Construction of the Resolvent Operator by means of Green's Function in the Adequate Hilbert Space

In this section, we define a linear operator in suitable Hilbert space in such a way that the considered problem can be interpreted as the eigenvalue problem of this operator. For this, we assume that and introduce a new inner product in the Hilbert space by for ,  .

Remark 1. Note that this modified inner product is equivalent to standard inner product of ; so is also Hilbert space.

For convenience, denote and define a linear operator with the domain consisting of all elements such that and are absolutely continuous in each interval and and has a finite limit and , , and .

Consequently the problems (2)–(5) can be written in the operator form as in the Hilbert space . It is easy to see that the operator is well defined in . Let be defined as above and let not be an eigenvalue of this operator. For construction of the resolvent operator , we will solve the operator equation for . This operator equation is equivalent to the nonhomogeneous differential equation on subject to nonhomogeneous boundary conditions and homogeneous transmission conditions Let . We already know that the general solution of (28) has the form (18). Putting this general solution in (29) yields Thus, the problems (28)-(29) have a unique solution Consequently, where and are the same with (21) and (22), respectively. From the equalities (13) and (21), it follows that By using (22), (31), and (33), we deduce that Consequently, the solution of the operator equation (27) has the form: From (34) and (35), it follows that where under Green's vector we mean Now, making use of (21), (34), (35), (36), and (37), we see that if not an eigenvalue of operator , then Hence, each nonreal is a regular point of an operator and Because of (38) and (40),

Theorem 2. The Resolvent operator is compact in the Hilbert space .

Proof. Let us define the operators , and by respectively. Then we can expressed the resolvent operator as . Since the linear operator is compact in the Hilbert space , the linear operator is compact in the Hilbert space . Compactness in is obvious. Therefore, the resolvent operator is also compact in .

4. Self-Adjoint Realization of the Problem

At first, we will prove the following lemmas.

Lemma 3. The domain is dense in .

Proof. Suppose that the element is orthogonal to . Denote by the set of infinitely differentiable functions on , each of which vanishes on some neighborhoods of the end-points , and . Since for , we have for all . Since is arbitrary, for all   and  , respectively. Taking into account that and are dense in and , respectively, we get that the function is equal to zero as element of . By choosing an element such that and putting in , we have . So . The proof is completed.

Lemma 4. The linear operator is symmetric in the Hilbert space .

Proof. Let ,  . By partial integration, we get From the definition of domain , we see easily that . The direct calculation gives Substituting these equalities in (45), we have so the operator is symmetric in . The proof is completed.

Remark 5. By Lemma 4, all eigenvalues of the problems (2)–(5) are real. Therefore, it is enough to investigate only real-valued eigenfunctions. Taking in view this fact, we can assume that the eigenfunctions are real-valued.

Corollary 6. If and are distinct eigenvalues of the problems (2)–(5), then the corresponding eigenfunctions and are orthogonal in the sense of the following equality:

Proof. The proof is immediate from the fact that the eigenelements and of the symmetric linear operator are orthogonal in the Hilbert space .

Theorem 7. The operator is self-adjoint in .

Proof. It is clear that the symmetry of a densely defined is equivalent to the condition for all . Notice that this implies that . If, in addition, we also have that , then is self-adjoint. Let . Then, by definition of , Let be any complex number for which . From this, it follows that Since any nonreal complex number is a regular point of , we can define the vector as Hence, Inserting this in (50) and recalling that is symmetric and , we have Consequently, Since is regular point of , we can choose . Inserting this in the last equality yields , and so , and therefore, . The proof is completed.

Remark 8. The main results of this study are derived in modified Hilbert space under simple condition (7). We can show that these conditions cannot be omitted. Indeed, let us consider the next special case of the problems (2)–(5): for which the condition (7) is not valid (). It is easy to verify that the operator corresponding to this problem is not symmetric in the classic Hilbert space under standard inner-product. Consider the following: Moreover, it is well known that the standard Sturm-liouville problems have infinitely many real eigenvalues. But it can be shown by direct calculation that the problem (55) has only one eigenvalue .