Research Article | Open Access

# Analysis of Sturm-Liouville Eigenproblem with Interior Singularities and a Perturbation Parameter

**Academic Editor:**Abdon Atangana

#### Abstract

We devote this work to the discussion underpinning the derivation of eigenvalues and eigenfunction solutions for Sturm-Liouville boundary value problems. The study reveals that the parameter dependent nonstandard Sturm-Liouville boundary value problem with interior singularities may have more than two turning points. The Titchmarsh-Weyl -function theory is applied here to obtain eigenfunction solutions valid for the whole interval in which pole singularities and two turning points are present. For the first time, with minimal constraints, the validity of the eigenfunction solutions are discussed when there are more than two turning points present. The eigenvalues are subsequently derived.

#### 1. Introduction

The nonstandard Sturm-Liouville (SL) boundary value problems with interior singularities and transition points may require a refinement of some perturbation method or a completely new approach to enable us obtain eigenfunction solutions valid within the entire interval. There have been recent interest in the study of both the standard and nonstandard Sturm-Liouville boundary value problems and other similarly perturbed systems (see [1–9]). Cases where the interior singularities are simple poles or double poles with supplementary one or two turning points (transition points) have been extensively studied (see [10–16]). One method that has been shown to solve some of the SL eigenproblems with interior singularities is the Titchmarsh-Weyl -function method [10, 11] (we refer the reader to whom the Titchmarsh-Weyl -function theory is unfamiliar to [17]). However, there has not been any presentation where these interior singularities are associated to more than two turning points which arise naturally from a nonrestricted parameter. The discussion underpinning the turning point analysis will be presented in Section 2.2. The motivation of this study is based on the fact that a typical example of a nonstandard Sturm-Liouville equation with a nonrestricted parameter occurs naturally in the high energy resonance contributing equations, which in this case happens to be the corresponding radial Schrödinger operator in . The fact that SL eigenproblems with interior singularities (which could either be a simple pole or a double pole singularity) could have a possibility of the presence of more than two turning points requires the introduction of a suitable constraint to enable the use of the Titchmarsh-Weyl -function technique. We will now give a brief background below of the nature of this problem and in subsequent sections use this example to demonstrate how the Titchmarsh-Weyl -function method can be applied to obtain eigenfunction solutions for a nonrestricted parameter dependent SL eigenproblem.

The resonance phenomenon appears in most material sciences and partly explains certain cases where structures such as bridges and buildings collapse (see [18]). Resonances are defined as poles of an analytic continuation of the quadratic form of the resolvent to the Riemann sheet through the branch cut along the continuous spectrum (see [19]). There have been studies ([20] and more recently [19]) on the investigation of resonances for the Hamiltonian operator. Asymptotic solutions to resonance equations of the Hamiltonian operator, acting on for , , , with appropriate conditions on a spherical symmetric support function , were presented in [19]. Eigenfunctions solutions for these operators are valid in the region . These solutions are singular at . On considering the operator acting on the desired eigenfunction solutions should satisfy boundary conditions for , where and/or could be infinite. In this case the eigenfunctions solutions presented in [19] would not be valid for the entire interval because of the pole singularities at and the behaviour of these functions (which are the Whittaker functions) as . In this work we use the following condition on the perturbation term : and we examine the case when From [19], the radial Schrödinger equation then takes the form where is the spectral parameter introduced in [19], such that and the Riemann sheet corresponds to the lower half plane of . With specification of we can maintain the definition of resonances as in Theorem of [19] as those points of the second Riemann sheet that can be represented in the form , where belongs to the lower half plane (that is ). It may be worth mentioning that even in a more simple case like the vibration of a circular membrane [18], the state of resonance is determined by the resulting radial equation. The objective of this work is not to solve the resonance equations as in [19] but to obtain eigenfunctions solutions and eigenvalues valid for the entire interval. We find that (4) has Whittaker functions as the eigenfunction solutions for as well as for , where is finite and different from zero. However, the case for , with boundary conditions such that where , presents difficulties in obtaining the eigenfunctions solutions as highlighted above. In Section 3 we present these solutions and in Section 4 we derive the eigenvalue relations. The eigenfunction solutions will of course satisfy (4) for and the restricted interval .

#### 2. Preliminaries and Turning Point Analysis

##### 2.1. Preliminaries with a Restricted Perturbation Term

Sequel to our discussion in the previous section we can rewrite (4) with the associated conditions as the following boundary value problem: where , and .

The differential equation (5) has the following turning points:

The value of the angular momentum is different from zero and is complex; the turning points and are distinct (i.e., do not coincide). Therefore, the eigenfunctions solutions of (5) will be in terms of some asymptotic approximation of Whittaker functions (see [21–24]). Using the transformation

equation (5) becomes

Consider whose eigenfunctions solution in terms of the and Whittaker functions is where and are constants that would be determined. On comparison of (9) and (10), we obtain

From whence we may then write the solution of (5) to be of the form where and are as in (12). It is obvious that for and ; and as presented in (13) give the two linearly independent solutions of the unperturbed equation as in [19]. However, considering an interval that includes the double pole singularity at and the possibility that and/or could be infinite results in the fact that for small, the internal boundary layer in which the asymptotic series for and are inaccurate includes the whole interval . This is as a result of the prominence in the effect of the double pole for small and the presence of the turning points that lie in the complex plane. Some of the behavioural properties of the and Whittaker functions are presented in [10, 25].

The method of solution we employ requires that we determine the quadrant in the complex plane on which the turning points lie. Considering the resonance condition , which implies that either and the fact that appears in both the discriminant and the denominator of (7) means that the turning points could lie in any of the four quadrants. Our method is applicable with the choice of branch cuts as in the theory of atmospheric and ocean waves. The locations of the turning points on the complex plane are as in the Figures 1, 2, 3, and 4. These figures also indicate the stokes lines and antistokes lines and how one of the zeros of the Whitaker function does lie on the positive real axis and the other one on the negative real axis which is essential for the boundary conditions to be satisfied.

Our method would require the use of Titchmarsh-Weyl -function [10, 17], so we shall only be interested in the turning points that are separated such that one lies in the left half of the complex plane and the other in the right half. The Titchmarsh-Weyl -function theory for (5) with two singular points , requires that we choose a regular point . We then let Denote the fundamental solution of (5), defined by where in this case . The Titchmarsh-Weyl -function at and is defined, respectively, by for . Both and are analytic functions in both and and obey the following relations:

Further, if we define by then it follows that

We use the above to construct our solutions in Section 3. Before doing so we present in the proceeding subsection a turning point analysis for the nonrestricted perturbation parameter.

##### 2.2. Turning Point Analysis for the Nonrestricted Perturbation Parameter

The turning point analysis for (notation used in the preceding and proceeding sections) would be similar to that presented in [10], so for the sake of brevity it would not be presented here. The corresponding Schrödinger equation for the spectral parameter and the perturbation parameter , is given by

We rewrite this in the form

For , we may consider the generalized interval , , with no restrictions on and . Furthermore, for the purpose of proper analysis of the turning points we will choose the boundary conditions as

The boundary value problems (22) and (23) have a double pole singularity, at , and three turning points which we now determine. The turning points are located at for which

We make the substitution to reduce (24) to the depressed cubic where and .

One of the roots of the depressed cubic is

We then obtain the remaining roots as

From (24) to (28) the three turning points are as follows:

The turning points coalesce into a single turning point when and there are two distinct turning points if . For any other values of there will be three distinct turning points located on the complex plane.

In order to apply the Titchmarsh-Weyl -function technique for reasons stated in Section 2.1, there should be two distinct turning points; therefore, we must have the restriction that

This provides the split intervals required for the application of the Titchmarsh-Weyl -function theory which will be similar to that when , ; provided of course that one turning point lies in the left half and the other in the right half of the complex plane.

#### 3. Eigenfunction Solutions

The analysis in Section 2 is essential for the application of the method of solutions we now provide here and in the section that follows. We will assume that the turning point lies in the left half plane and lies in the right half plane. We split the interval as follows: where , .

We obtain our solution for (5) as a linear combination of eigenfunctions for the split intervals and it turns out that the boundary conditions (6) are not necessary as they are automatically satisfied by our solution. The following theorem summarizes our first results.

Theorem 1. *Let and be appropriate combinations of the eigenfunctions solutions of (5) with (6) for small such that
**Then the solutions to (5) would be
**
where , , , are as defined in Section 2.1.*

*Proof. *For convenience we suppress the , dependence of and . We start with the regular solutions and of (5) on , where , , and .

Thus we let

We now proceed to determine , , , and .

Substituting (34) into (16) we obtain

From (36) we obtain the following:
where

Similarly, substitution of (35) into (16) yields
where

Using (17) we may then find and , respectively, at and as follows:
where and are as given in (34) and (35) with as obtained in (37) to (40). From (19) we then have

Therefore, the eigenfunctions for are
which automatically satisfy the boundary conditions at .

Now, for the interval , we let

Proceeding as above we obtain
where

Using the same process we find that
where

By using (44) to (48) we obtain
where

Therefore, the eigenfunctions for are
and these automatically satisfy the boundary conditions at . From (43) and (51) we find that

We then conclude that the eigenfunctions solutions to (5) with the accompanying boundary conditions would be

This completes the proof of the theorem.

#### 4. Computation of the Eigenvalues

As in [11], the eigenvalues for our problem are . To derive the eigenvalue relations for boundary value problems for which the eigenfunctions are obtained using a split interval technique requires a careful examination of the behaviour of the eigenfunctions. To enable us achieve our goal for the interval , , we will need the following conditions:(i)that the eigenfunctions solutions are valid for (ii)that (iii)and that

We summarise our results in this section in the following theorem.

Theorem 2. *The nonstandard SL eigenproblem (or the Schrödinger radial operator in ), (5) with (6) with separated turning points, which have solutions (53) and satisfy conditions (54)–(56) have its eigenvalues given by*

*Proof. *The asymptotic expansions of the eigenfunctions that meet the requirement (54) are
where as usual and .

We then express the asymptotic approximation of the function as follows:

Through a careful limiting process we find that

a result that comes in handy.

Using (60) and (59), we apply the boundary conditions at and to obtain

By employing the conditions in (55) and (56), we obtain from (61) that

Thus, the eigenvalues may be expressed as where is the mode number. The eigenvalues are complex, which is consistent with the fact that in the radial Schrödinger equation is complex. This result as expected is similar to the result obtained in [11] since the perturbation term which is an integral part of is independent of .

#### 5. Concluding Remarks

The boundary value problem with interior singularities and turning points for which resonance equations result may be represented as nonstandard Sturm-Liouville eigenproblems. With the asymptotic expansions of the eigenfunctions solutions for being the same for , we have in this paper presented solutions valid for an unrestricted interval. The three turning points analysis presented here when the perturbation term is a function of demonstrate the applicability of the method in these cases. The beauty of this split interval technique which employs the use of the Titchmarsh-Weyl -function theory is that eigenfunctions solutions could also be obtained for the nonconstant perturbation term which happens to have three turning points provided one of these turning points coalesce with another, thus permitting us to use the branch cuts as in the theory of atmospheric and oceanic waves.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### References

- D. Eric,
*Exact WKB Analysis near a Simple Turning Point. Algebraic Analysis of Differential Equations from Micro-Local Analysis to Exponential Asymtotics*, Springer, Tokyo, Japan, 2008. - A. Fruchard and R. Schäfke, “Composite asymptotic expansions and turning points of singularly perturbed ordinary differential equations,”
*Comptes Rendus Mathematique*, vol. 348, no. 23-24, pp. 1273–1277, 2010. View at: Publisher Site | Google Scholar | MathSciNet - H.-R. Sun, Y.-N. Li, J. J. Nieto, and Q. Tang, “Existence of solutions for Sturm-Liouville boundary value problem of impulsive differential equations,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 707163, 19 pages, 2012. View at: Publisher Site | Google Scholar | MathSciNet - L. Zhang, X. Huang, and C. Xing, “Solvability of Sturm-Liouville boundary value problems with impulses,”
*Boundary Value Problems*, vol. 2013, article 192, 2013. View at: Publisher Site | Google Scholar | MathSciNet - M. Jiaqi, “Singularly perturbed solution of boundary value problem for nonlinear equations of fourth order with two parameters,”
*Advances in Mathematics*, vol. 39, no. 6, pp. 736–740, 2010. View at: Google Scholar | MathSciNet - M. Ni and Z. Wang, “On step-like contrast structure of singularly perturbed systems,”
*Boundary Value Problems*, vol. 2009, Article ID 634324, 17 pages, 2009. View at: Publisher Site | Google Scholar | MathSciNet - Y. Tian and W. Ge, “Multiple positive solutions for a second order Sturm-Liouville boundary value problem with a $p$-Laplacian via variational methods,”
*The Rocky Mountain Journal of Mathematics*, vol. 39, no. 1, pp. 325–342, 2009. View at: Publisher Site | Google Scholar | MathSciNet - Y. Tian and W. Ge, “Multiple solutions of impulsive Sturm-Liouville boundary value problem via lower and upper solutions and variational methods,”
*Journal of Mathematical Analysis and Applications*, vol. 387, no. 2, pp. 475–489, 2012. View at: Publisher Site | Google Scholar | MathSciNet - Y. Tian and W. Ge, “Variational methods to Sturm-Liouville boundary value problem for impulsive differential equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 72, no. 1, pp. 277–287, 2010. View at: Publisher Site | Google Scholar | MathSciNet - T. M. Acho and D. P. Clemence, “Sturm-Liouville eigenproblems with an interior double pole,”
*Zeitschrift für Angewandte Mathematik und Physik*, vol. 46, no. 3, pp. 459–474, 1995. View at: Publisher Site | Google Scholar | MathSciNet - T. M. Acho, “Derivation of eigenrelations for the Sturm-Liouville boundary value problems with interior singularities,”
*Quarterly of Applied Mathematics*, vol. 65, no. 2, pp. 375–383, 2007. View at: Google Scholar | MathSciNet - J. P. Boyd, “Sturm-Liouville eigenproblems with an interior pole,”
*Journal of Mathematical Physics*, vol. 22, no. 8, pp. 1575–1590, 1981. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. V. Atkinson, N. N. Everitt, and A. Zettl, “Regularisation of a Sturm-Liouville problem with an interior singularity,”
*Differential Integral Equations*, vol. 1, no. 2, pp. 213–221, 1988. View at: Google Scholar - J. P. Boyd, “Complex coordinate methods for hydrodynamic instabilities and Sturm-Liouville eigenproblems with an interior singularity,”
*Journal of Computational Physics*, vol. 57, no. 3, pp. 454–471, 1985. View at: Publisher Site | Google Scholar | MathSciNet - T. M. Acho, “A matching technique and structure theorems for Sturm-Liouville boundary value problems with interior singularities,”
*Quaestiones Mathematicae*, vol. 25, no. 3, pp. 275–287, 2002. View at: Publisher Site | Google Scholar | MathSciNet - T. M. Acho and D. P. Clemence, “The parameter dependent Sturm-Liouville eigenproblem with an interior simple or double pole,”
*The ANZIAM Journal*, vol. 43, no. 4, pp. 479–491, 2002. View at: Google Scholar | MathSciNet - E. A. Coddington and N. Levinson,
*Theory of Ordinary Differential Equations*, R.E. Krieger, Malabar, Fla, USA, 1955. View at: MathSciNet - Y. Pinchover and J. Rubinstein,
*An Introduction to Partial Differential Equations*, Cambridge University Press, Cambridge, UK, 2007. View at: MathSciNet - M. A. Shubov, “High-energy asymptotics of resonances for three-dimensional Schrödinger operator with screened Coulomb potential,”
*Journal of Mathematical Physics*, vol. 35, no. 2, pp. 656–679, 1994. View at: Publisher Site | Google Scholar | MathSciNet - R. G. Newton,
*Scattering Theory of Waves and Particles*, Springer, New York, NY, USA, 2nd edition, 1982. View at: MathSciNet - L. J. Slater,
*Confluent Hypergeometric Functions*, Cambridge University Press, 1960. View at: MathSciNet - N. Bleistein and R. A. Handelsman,
*Asymptotic Expansion of Integrals*, Dover Publications, Ney York, NY, USA, 1986. View at: MathSciNet - G. F. Carrier, M. Krook, and C. E. Pearson,
*Functions of a Complex Variable*, vol. 49 of*Classics in Applied Mathematics*, SIAM, Philadelphia, Pa, USA, 2005. View at: Publisher Site | MathSciNet - A. H. Nayfeh,
*Perturbation Methods*, John Wiley & Sons, New York, NY, USA, 1973. View at: MathSciNet - M. Abramowitz and I. Stegun,
*Handbook of Mathematical Functions*, Dover, New York, NY, USA, 1965.

#### Copyright

Copyright © 2014 Thomas Mbah Acho. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.