SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research Article

Hierarchies of Difference Boundary Value Problems

Sonja Currie* and AnneD Love

Author Affiliations

School of Mathematics, University of the Witwatersrand, Private Bag 3, Wits 2050, South Africa

For all author emails, please log on.

Boundary Value Problems 2011, 2011:743135  doi:10.1155/2011/743135


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2011/1/743135


Received:25 November 2010
Accepted:11 January 2011
Published:18 January 2011

© 2011 Sonja Currie and Anne D. Love.

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.

This paper generalises the work done in Currie and Love (2010), where we studied the effect of applying two Crum-type transformations to a weighted second-order difference equation with various combinations of Dirichlet, non-Dirichlet, and affine -dependent boundary conditions at the end points, where is the eigenparameter. We now consider general -dependent boundary conditions. In particular we show, using one of the Crum-type transformations, that it is possible to go up and down a hierarchy of boundary value problems keeping the form of the second-order difference equation constant but possibly increasing or decreasing the dependence on of the boundary conditions at each step. In addition, we show that the transformed boundary value problem either gains or loses an eigenvalue, or the number of eigenvalues remains the same as we step up or down the hierarchy.

1. Introduction

Our interest in this topic arose from the work done on transformations and factorisations of continuous Sturm-Liouville boundary value problems by Binding et al. [1] and Browne and Nillsen [2], notably. We make use of analogous ideas to those discussed in [35] to study difference equations in order to contribute to the development of the theory of discrete spectral problems.

Numerous efforts to develop hierarchies exist in the literature, however, they are not specifically aimed at difference equations per se and generally not for three-term recurrence relations. Ding et al., [6], derived a hierarchy of nonlinear differential-difference equations by starting with a two-parameter discrete spectral problem, as did Luo and Fan [7], whose hierarchy possessed bi-Hamiltonian structures. Clarkson et al.'s, [8], interest in hierarchies lay in the derivation of infinite sequences of systems of difference equations by using the Bcklund transformation for the equations in the second Painlev equation hierarchy. Wu and Geng, [9], showed early on that the hierarchy of differential-difference equations possesses Hamiltonian structures while a Darboux transformation for the discrete spectral problem is shown to exist.

In this paper, we consider a weighted second-order difference equation of the form

(11)

where represents a weight function and a potential function.

Our aim is to extend the results obtained in [10, 11] by establishing a hierarchy of difference boundary value problems. A key tool in our analysis will be the Crum-type transformation (2.1). In [10], it was shown that (2.1) leaves the form of the difference equation (1.1) unchanged. For us, the effect of (2.1) on the boundary conditions will be crucial. We consider (eigenparameter)-dependent boundary conditions at the end points. In particular, the eigenparameter dependence at the initial end point will be given by a positive Nevanlinna function, say, and at the terminal end point by a negative Nevanlinna function, say. The case of was covered in [10] and the the case of constant was studied in [11]. Applying transformation (2.1) to the boundary conditions results in a so-called transformed boundary value problem, where either the new boundary conditions have more -dependence, less -dependence, or the same amount of -dependence as the original boundary conditions. Consequently the transformed boundary value problem has either one more eigenvalue, one less eigenvalue, or the same number of eigenvalues as the original boundary value problem. Thus, it is possible to construct a chain, or hierarchy, of difference boundary value problems where the successive links in the chain are obtained by applying the variations of (2.1) given in this paper. For instance, it is possible to go from a boundary value problem with -dependent boundary conditions to a boundary value problem with -independent boundary conditions or vice versa simply by applying the correct variation of (2.1) an appropriate number of times. Moreover, at each step, we can precisely track the eigenvalues that have been lost or gained. Hence, this paper provides a significant development in the theory of three-term difference boundary value problems in regard to singularities and asymptotics in the hierarchy structure. For similar results in the continuous case, see [12].

There is an obvious connection between the three-term difference equation and orthogonal polynomials. In fact, the three-term recurrence relation satisfied by orthogonal polynomials is perhaps the most important information for the constructive and computational use of orthogonal polynomials [13].

Difference equations and operators and results concerning their existence and construction of their solutions have been discussed in [14, 15]. Difference equations arise in numerous settings and have applications in diverse areas such as quantum field theory, combinatorics, mathematical physics and biology, dynamical systems, economics, statistics, electrical circuit analysis, computer visualization, and many other fields. They are especially useful where recursive computations are required. In particular see [16] [9, Introduction] for three physical applications of the difference equation (1.1), namely, the vibrating string, electrical network theory and Markov processes, in birth and death processes and random walks.

It should be noted that G. Teschl's work, [17, Chapter  11], on spectral and inverse spectral theory of Jacobi operators, provides an alternative factorisation, to that of [10], of a second-order difference equation, where the factors are adjoints of one another.

This paper is structured as follows.

In Section 2, all the necsessary results from [10] are recalled, in particular how (1.1) transforms under (2.1). In addition, we also recap some important properties of Nevanlinna functions.

The focus of Section 3 is to show exactly the effect that (2.1) has on boundary conditions of the form

(12)

We give explicitly the new boundary conditions which are obeyed, from which it can be seen whether the -dependence has increased, decreased, or remained the same.

Lastly, in Section 4, we compare the spectrum of the original boundary value problem with that of the transformed boundary value problem and show under which conditions the transformed boundary value problem has one more eigenvalue, one less eigenvalue, or the same number of eigenvalues as the original boundary value problem.

2. Preliminaries

In [10], we considered (1.1) for , where the values of and are given by boundary conditions, that is, is defined for .

Let the mapping be defined by

(21)

where, throughout this paper, is a solution to (1.1) for such that for all . Whether or not obeys the various given boundary conditions (to be specified later) is of vital importance in obtaining the results that follow.

From [10], we have the following theorem.

Theorem 2.1.

Under the mapping (2.1), (1.1) transforms to

(22)

where for

(23)

We now recall some properties of Nevanlinna functions.

(I) The inverse of a positive Nevanlinna function is a negative Nevanlinna function, that is

(24)

where are positive Nevanlinna functions. This follows directly from the fact that if and only if .

(II) If

(25)

then

(26)

This follows by (I) together with the fact that since has zeros has poles. Also as so as . Thus, if is a positive Nevanlinna function of the form (2.5), then for , is a negative Nevanlinna function of the same form.

(III) If

(27)

then

(28)

since has zeros so has poles and as so as .

For the remainder of the paper, will denote a Nevanlinna function where

is the number of terms in the sum;

indicates the value of at which the boundary condition is imposed and

(29)

3. General -Dependent Boundary Conditions

In this section, we show how obeying general -dependent boundary conditions transforms, under (2.1), to obeying various types of -dependent boundary conditions. The exact form of these boundary conditions is obtained by considering the number of zeros and poles (singularities) of the various Nevanlinna functions under discussion and these correlations are illustrated in the different graphs depicted in this section.

Lemma 3.1.

If obeys the boundary condition

(31)

then the domain of may be extended from to by forcing the condition

(32)

where

(33)

with .

Proof.

The transformed equation (2.2), for , together with (3.2) gives

(34)

Also the mapping (2.1), together with (3.1), yields

(35)

Substituting (3.5) into (3.4), we obtain

(36)

Now (2.1), with , gives

(37)

which when substituted into (3.6) and dividing through by results in

(38)

This may be rewritten as

(39)

Using (1.1), with , together with (3.1), gives

(310)

Subtracting (3.10) from (3.9) results in

(311)

Rearranging the above equation and dividing through by yields

(312)

and hence

(313)

Thus obeys the equation on the extended domain.

The remainder of this section illustrates why it is so important to distinguish between the two cases of obeying or not obeying the boundary conditions.

Theorem 3.2.

Consider obeying the boundary condition (3.1) where is a positive Nevanlinna function, that is, for . Under the mapping (2.1), obeying (3.1) transforms to obeying (3.2) as follows.

If does not obey (3.1) then obeys

(314)

(315)

If does obey (3.1) for then obeys

(316)

(317)

where , that is, are positive Nevanlinna functions.

In (A) and (B), is not possible.

Proof.

The fact that is by construction, see Lemma 3.1. We now examine the form of in Lemma 3.1. Let , , and then

(318)

But

(319)

thus

(320)

Now has the expansion

(321)

where and the 's correspond to where , that is, the singularities of (3.20).

Since is a positive Nevanlinna function it has a graph of the form shown in Figure 1.

Clearly, the gradient of at is positive for all , that is,

(322)

If does not obey (3.1), then the zeros of

(323)

are the poles of , that is, the 's and where for . It is evident, from Figure 1, that the number of 's is equal to the number of 's, thus in (3.21), .

We now examine the form of in (3.21). As it follows that . Thus

(324)

Therefore

(325)

Hence, substituting into (3.20) gives

(326)

Let

(327)

Then since , and we have that and clearly if then giving (3.14), that is,

(328)

If then we want so that we have a positive Nevanlinna function, that is

(329)

which means that either,

(330)

giving that, since ,

(331)

which is as shown in Figure 1, or,

(332)

giving that

(333)

but this means that which is not possible.

Thus, for , that is, given , the ratio must be chosen suitably to ensure that is a positive Nevanlinna function as required. Hence we obtain (3.15), that is

(334)

If obeys (3.1), for , then . Thus in Figure 1, one of the 's is equal to and since is less than the least eigenvalue of the boundary value problem (1.1), (3.1) together with a boundary condition at (specified later) it follows that , as for all .

Now

(335)

and as

(336)

Thus is a removable singularity. Alternatively,

(337)

which illustrates that the singularity at is removable.

We now have that the number of nonremovable singularities, , in (3.20) is one less than the number of 's , see Figure 1. Thus (3.21) becomes

(338)

which may be rewritten as

(339)

where , for .

We now examine the form of in (3.39). As , we have that, as before, . Thus

(340)

Hence, from (3.20),

(341)

Let

(342)

Then since , and we have that and clearly if then giving (3.16), that is,

(343)

If then we need so that we have a positive Nevanlinna function, that is

(344)

which means that either

(345)

giving that, since ,

(346)

which is as shown in Figure 1, or,

(347)

giving that

(348)

but this means that which is not possible.

Thus, for , that is, given , the ratio must be chosen suitably to ensure that is a positive Nevanlinna function as required. Hence, we obtain (3.17), that is,

(349)

In the theorem below, we increase the dependence by introducing a nonzero term in the original boundary condition. As in Theorem 3.2, the dependence of the transformed boundary condition depends on whether or not obeys the given boundary condition. In addition, to ensure that the dependence of the transformed boundary condition is given by a positive Nevanlinna function it is necessary that the transformed boundary condition is imposed at 0 and 1 as opposed to −1 and 0. Thus the interval under consideration shrinks by one unit at the initial end point. By routine calculation it can be shown that the form of the dependence of the transformed boundary condition, if imposed at −1 and 0, is neither a positive Nevalinna function nor a negative Nevanlinna function.

Theorem 3.3.

Consider obeying the boundary condition

(350)

where is a positive Nevanlinna function, that is, and for . Under the mapping (2.1), obeying (3.50) transforms to obeying the following.

(1) If does not obey (3.50) then obeys

(351)

(2) If does obey (3.50), for , then obeys

(352)

where .

Proof.

Since and are defined we do not need to extend the domain in order to impose the boundary conditions (3.51) or (3.52).

The mapping (2.1), at , together with (3.50) gives

(353)

Also (2.1), at , is

(354)

Substituting in for from (1.1), with , and using (3.50), we obtain that

(355)

From (3.53) and (3.55), it now follows that

(356)

As in Theorem 3.2, let and . Then (3.56) becomes

(357)

From Theorem 3.2, we have that so

(358)

Also, as in Theorem 3.2,

(359)

has the expansion

(360)

where corresponds to , that is, the singularities of (3.59). Now is a positive Nevanlinna function with graph given in Figure 2.

Clearly, the gradient of at is positive for all , that is,

(361)

If does not obey (3.50) then the zeros of

(362)

are the poles of , that is, the 's and where for . It is evident, from Figure 2, that the number of 's is one more than the number of 's, thus in (3.60), .

We now examine the form of in (3.60). As it follows that , thus

(363)

Hence, .

Using (3.58) we now obtain

(364)

Note that . Let

(365)

then

(366)

Now since if then , that is, but and so this is not possible. Therefore by Section 2, Nevanlinna result (II), we have that

(367)

that is, (3.51) holds.

If does obey (3.50) for then . Thus, in Figure 2, one of the 's, is equal to and since is less than the least eigenvalue of the boundary value problem (1.1), (3.50) together with a boundary condition at (specified later) it follows that , as for all .

Now (3.59) can be written as

(368)

and as

(369)

Thus is a removable singularity. Alternatively, we could substitute in for and to illustrate that the singularity at is removable, see Theorem 3.2. Hence the number of nonremovable 's is the same as the number of 's, see Figure 2. So (3.60) becomes

(370)

which may be rewritten as

(371)

where and for .

We now examine the form of in (3.70). As , we have that , thus

(372)

Hence, . So, from (3.58) with , we obtain

(373)

where, as before,

(374)

Thus, by Section 2, Nevanlinna result (II), we have that

(375)

that is, (3.52) holds.

In Theorem 3.4, we impose a boundary condition at the terminal end point and show how it is transformed according to whether or not obeys the given boundary condition.

Theorem 3.4.

Consider obeying the boundary condition at given by

(376)

where is a negative Nevanlinna function, that is, and for . Under the mapping (2.1), obeying (3.76) transforms to obeying the following.

(I) If does not obey (3.76) then obeys

(377)

(II) If does obey (3.76) then obeys

(378)

where .

Proof.

Since and are defined we do not need to extend the domain of in order to impose the boundary conditions (3.77) or (3.78).

The mapping (2.1), at , gives

(379)

From (1.1), with , we can substitute in for in the above equation to get

(380)

Using (3.76), we obtain

(381)

But obeys (1.1) at , for , so that (3.81) becomes

(382)

Also, for , (2.1) together with (3.76) yields

(383)

Therefore,

(384)

Let , then (3.84) may be rewritten as

(385)

By Section 2, Nevanlinna result (I), since is a negative Nevanlinna function it follows that is a positive Nevanlinna function, which has the form

(386)

by Section 2, Nevanlinna result (III).

As before has expansion

(387)

where , , corresponds to the singularities of (3.85), that is, where . The graph of is as shown in Figure 3.

As before, the gradient of at is positive for all , that is

(388)

If does not obey (3.76) then the zeros of

(389)

are the poles of , that is, the 's and where for . Clearly, from Figure 3, the number of 's is the same as the the number of 's, thus in (3.87), .

Next, we examine the form of in (3.87). As it follows that . Thus

(390)

Therefore, . Hence,

(391)

where , , and for , which is precisely (3.77).

If does obey (3.76) for then . Thus in Figure 3, one of the 's, is equal to and since is less than the least eigenvalue of the boundary value problem (1.1), (3.76) together with a boundary condition at −1 (as given in Theorems 3.2 or 3.3) it follows that , as for all .

Now

(392)

and as

(393)

Thus is a removable singularity. Again, alternatively, we could have substituted in for and to illustrate that the singularity at is removable, see Theorem 3.2. Hence the number of nonremovable 's is one less than the number of 's, see Figure 3.

So (3.87) becomes

(394)

which may be rewritten as

(395)

where and for .

Now as ,

(396)

So, we obtain

(397)

where , , , and for all , that is, we obtain (3.78).

4. Comparison of the Spectra

In this section, we investigate how the spectrum of the original boundary value problem compares to the spectrum of the transformed boundary value problem. This is done by considering the degree of the eigenparameter polynomial for the various eigenconditions.

Lemma 4.1.

Consider the boundary value problem given by (1.1) for together with boundary conditions

(41)

(42)

Then the boundary value problem (1.1), (4.1), (4.2) has eigenvalues. (Note that the number of unit intervals considered is .)

Proof.

From (1.1), with , we obtain

(43)

Substituting in for from (4.1) yields

(44)

which may be rewritten as

(45)

where , are real constants.

Now (1.1), for , together with (4.5) results in

(46)

where , are real constants.

Thus, by induction,

(47)

for real constants , . Similarly

(48)

for real constants , .

Since , using boundary condition (4.2) we obtain the following eigencondition:

(49)

where , , are real constants.

Thus, the numerator is a polynomial, in , of order . Note that, none of the roots of this polynomial are given by , or , since, from Figures 1 to 3, it is easy to see that none of the eigenvalues of the boundary value problem are equal to the poles of the boundary conditions. Also is not a problem as the curve of the Nevanlinna function never intersects with the horizontal or oblique asymptote. This means that there are no common factors to cancel out. Hence the eigencondition has roots giving that the boundary value problem has eigenvalues.

As a direct consequence of Theorems 2.1, 3.2, 3.3, 3.4, and Lemma 4.1 we have the following theorem.

Theorem 4.2.

For the original boundary value problem we consider twelve cases, (see Table 1 in the Appendix), each of which has s+l+m+1 eigenvalues. The corresponding transformed boundary value problem for each of the twelve cases, together with the number of eigenvalues for that transformed boundary value problem, is given in Table 1 (see the appendix).

Remark 4.3.

To summarise we have the following.

(a) If obeys the boundary conditions at both ends the transformed boundary value problem will have one less eigenvalue than the original boundary value problem, namely, .

(b) If obeys the boundary condition at one end only the transformed boundary value problem will have the same eigenvalues as the original boundary value problem.

(c) If does not obey any of the boundary conditions the transformed boundary value problem will have one more eigenvalue than the original boundary value problem, namely, .

Corollary 4.4.

If are the eigenvalues of any one of the original boundary value problems (1)–(9), in Theorem 4.2, with corresponding eigenfunctions then

(i) are the eigenvalues of the corresponding transformed boundary value problems (1)–(3), in Theorem 4.2, with corresponding eigenfunctions ;

(ii) are the eigenvalues of the corresponding transformed boundary value problems (4)–(9), in Theorem 4.2, with corresponding eigenfunctions .

Also, if are the eigenvalues of any one of the original boundary value problems (10)–(12), in Theorem 4.2, with corresponding eigenfunctions then are the eigenvalues of the corresponding transformed boundary value problems (10)–(12), in Theorem 4.2, with corresponding eigenfunctions .

Proof.

By Theorems 2.1, 3.2, 3.3, and 3.4, we have that (2.1) transforms eigenfunctions of the original boundary value problems (1)–(9) to eigenfunctions of the corresponding transformed boundary value problems. In particular, if are the eigenvalues of one of the original boundary value problems, (1)–(9), with eigenfunctions then

(i) are the eigenfunctions of the corresponding transformed boundary value problem, (1)–(3), with eigenvalues . Since the transformed boundary value problems, (1)–(3), have eigenvalues it follows that constitute all the eigenvalues of the transformed boundary value problem;

(ii) are the eigenfunctions of the corresponding transformed boundary value problem, (4)–(9), with eigenvalues . Since the transformed boundary value problems, (4)–(9), have eigenvalues it follows that constitute all the eigenvalues of the transformed boundary value problem.

Also, again by Theorems 2.1, 3.2, 3.3, and 3.4, we have that (2.1) transforms eigenfunctions of the original boundary value problems (10)–(12) to eigenfunctions of the corresponding transformed boundary value problems. In particular, if are the eigenvalues of one of the original boundary value problems, (10)–(12), with eigenfunctions then are the eigenfunctions of the corresponding transformed boundary value problem, (10)–(12), with eigenvalues . Since the transformed boundary value problems, (10)–(12), have eigenvalues it follows that constitute all the eigenvalues of the transformed boundary value problem.

Acknowledgments

The authors would like to thank Professor Bruce A. Watson for his useful input. S. Currie is supported by NRF Grant nos. TTK2007040500005 and FA2007041200006.

References

  1. Binding, PA, Browne, PJ, Watson, BA: Spectral isomorphisms between generalized Sturm-Liouville problems. Linear Operators and Matrices, Operator Theory: Advances and Applications, pp. 135–152. Birkhäuser, Basel, Switzerland (2002)

  2. Browne, PJ, Nillsen, RV: On difference operators and their factorization. Canadian Journal of Mathematics. 35(5), 873–897 (1983). Publisher Full Text OpenURL

  3. Binding, PA, Browne, PJ, Watson, BA: Sturm-Liouville problems with boundary conditions rationally dependent on the eigenparameter. II. Journal of Computational and Applied Mathematics. 148(1), 147–168 (2002). Publisher Full Text OpenURL

  4. Binding, PA, Browne, PJ, Watson, BA: Sturm-Liouville problems with reducible boundary conditions. Proceedings of the Edinburgh Mathematical Society. Series II. 49(3), 593–608 (2006). Publisher Full Text OpenURL

  5. Binding, PA, Browne, PJ, Watson, BA: Transformations between Sturm-Liouville problems with eigenvalue dependent and independent boundary conditions. Bulletin of the London Mathematical Society. 33(6), 749–757 (2001). Publisher Full Text OpenURL

  6. Ding, H-Y, Sun, Y-P, Xue, F-C: A hierarchy of differential-difference equations, conservation laws and new integrable coupling system. Communications in Nonlinear Science and Numerical Simulation. 15(8), 2037–2043 (2010). Publisher Full Text OpenURL

  7. Luo, L, Fan, E-G: A hierarchy of differential-difference equations and their integrable couplings. Chinese Physics Letters. 24(6), 1444–1447 (2007). Publisher Full Text OpenURL

  8. Clarkson, PA, Hone, ANW, Joshi, N: Hierarchies of difference equations and Bäcklund transformations. Journal of Nonlinear Mathematical Physics. 10(supplement 2), 13–26 (2003). Publisher Full Text OpenURL

  9. Wu, Y, Geng, X: A new hierarchy of integrable differential-difference equations and Darboux transformation. Journal of Physics A. 31(38), L677–L684 (1998). Publisher Full Text OpenURL

  10. Currie, S, Love, AD: Transformations of difference equations I. Advances in Difference Equations. 2010, (2010)

  11. Currie, S, Love, AD: Transformations of difference equations II. Advances in Difference Equations. 2010, (2010)

  12. Binding, PA, Browne, PJ, Watson, BA: Sturm-Liouville problems with boundary conditions rationally dependent on the eigenparameter. I. Proceedings of the Edinburgh Mathematical Society. Series II. 45(3), 631–645 (2002)

  13. Gautschi, W: Orthogonal Polynomials: Computation and Approximation, Numerical Mathematics and Scientific Computation,p. x+301. Oxford University Press, New York, NY, USA (2004)

  14. Miller, KS: Linear Difference Equations,p. x+105. W. A. Benjamin, New York, NY, USA (1968)

  15. Miller, KS: An Introduction to the Calculus of Finite Differences and Difference Equations,p. viii+167. Dover, New York, NY, USA (1966)

  16. Atkinson, FV: Discrete and Continuous Boundary Problems, Mathematics in Science and Engineering,p. xiv+570. Academic Press, New York, NY, USA (1964)

  17. Teschl, G: Jacobi Operators and Completely Integrable Nonlinear Lattices, Mathematical Surveys and Monographs,p. xvii+351. American Mathematical Society, Providence, RI, USA (2000)

Appendix

Twelve Cases for Theorem 4.2

See Table 1.