arxiv: v1 [math.sp] 15 Dec 2014

Transcription

1 An Application of Value Distribution Theory to Schrödinger Operators with Absolutely Continuous Spectrum arxiv: v [math.sp] 5 Dec 204 CHARLES FULTON Department of Mathematical Sciences Florida Insitute of Technology Melbourne, FL DAVID PEARSON Department of Mathematics University of Hull Cottingham Road Hull, HU6 7RX England United Kingsom STEVEN PRUESS 33 N. Desert Deer Pass Green Valley, Arizona Keywords: Value distribution, Schrödinger Operator, absolutely continuous spectrum, uniformly asymptotically distributed, Bessel function Abstract We use the Value Distribution theory developed by Pearson and Breimesser to obtain a sequence of functions in the eigenvalue parameter for some Sturm-Liouville problems which have the property of being uniformly asymptotically distributed. Introduction In this paper we apply a general result of Breimesser and Pearson[] to obtain a very general asymptotic result for Schrödinger equations with one regular endpoint and one Limit Point endpoint which generates some absolutely continuous spectrum. The software package SLEDGE [2, 4], developed by Pruess and Fulton, has as one of its options the computation of the spectral function for Sturm-Liouville problems having continuous spectrum. For cases where the spectrum is known to be absolutely continuous, more recent research codes described in [5, 6, 7] have proven to be much better than the SLEDGE package. This paper will, in the opinion of the authors, provide a backgound for development of even further formulas and numerical algorithms for the spectral density function in cases having some a.c. spectrum. The general theory of value distribution which is providing the new approach described here has been developed by Pearson and Breimesser [8, 9, 0, ]. For simplicity we restrict our attention to the problem for x [a, ), a RI, with a Dirichlet boundary condition at x = a:

2 y +qx)y = y, a x <,.) ya) = 0..2) We pose the assumptions: i) qx) is continuous in [a, ), ii) the LP/O-N or LP/O case occurs at x =, iii) there exist one or more intervals where a.c. spectrum is known to exist. See Fulton and Pruess [3, p. 4] for these general endpoint classifications: here O means equation.) is oscillatory at for all real and O-N means there exists RI such that equation.) is oscillatory at for all, ) and nonoscillatory at for all,). The LP/O-N case includes potentials which satisfy lim qx) = L for some finite limit = L in this case) and periodic potentials qx), as well as all quasi-periodic and almost periodic potentials. The LP/O case includes potentials which are LP at and satisfy lim qx) = In all of the above mentioned cases, it is possible for a.c. spectrum to occur. We now state a very general theorem on Herglotz functions of Breimesser and Pearson[, Thm,p.40] which will be applied to the Titchmarsh-Weyl m-function associated with the above self-adjoint Sturm- Liouville problem.)-.2) in the next section: Theorem Breimesser and Pearson, 2003) Let fz) be an arbitrary Herglotz function. Let S be an arbitrary Borel subset of RI, and let be any set of finite measure. Define the angle θ subtended by S at z for all z CI CI = complex numbers) by ) θs,z) = Im dt..3) t z Letting the boundary value of the Herglotz function on the real -axis be S f + ) := lim ǫ 0 [f+iǫ)].4) we define, and Then we have ω,s,f + ) := π θs,f+ ),.5) ω,s,f+iǫ)) := θs,f+iǫ))..6) π π ω,s,f + )d ω,s,f+iǫ)d.7) θ c,+iǫ)d 0, as ǫ 0. where θ C,+iǫ) is the angle subtended at +iǫ by C = RI \. 2

3 2 A Uniformly Asymptotically Distributed Sequence of Functions We define a fundamental system {u,v} of solutions of the Schrödinger equation.) at x = a by the initial conditions [ ] [ ] ua,) va,) 0 u a,) v =, CI. 2.) a,) 0 The Titchmarsh-Weyl m function for the problem.)-.2) is defined for Im z 0) by ux,z)+mz)vx,z) L 2 a, ). 2.2) The boundary value of the m-function on the real -axis is then defined by m + ) := lim ǫ 0 [m+iǫ)] = A)+iB), 2.3) wherewe haveintroducedthe real-valuedfunctions A), B) asthe realandimaginarypartsofm + ). The standard method of approximating spectral quantities for the half line problem.)-.2) is to make use of spectral quantities associated with the truncated regular problem on [a,b] for b a, ): y +qx)y = y, a x b 2.4) ya) = 0 2.5) yb) = 0 2.6) We let m b z) be the m-function associated with the Sturm-Liouville problem 2.4)-2.6), that is, for all z not eigenvalues, ub,z)+m b z)vb,z) = 0, or m b z) := ub,z) vb,z). 2.7) The boundary value of this meromorphic m b function on the real -axis is then well-defined for all not eigenvalues of 2.4)-2.6) as m + b := lim [m b+iǫ)] = ub,) ǫ 0 vb,). 2.8) Here, m + b ) is necessarily real-valued. As in Theorem, we let S be any Borel subset of RI and be any finite interval in RI. For x [a, ) and all RI, not eigenvalues with b = x) of 2.4)-2.6), we put and define Fx,) := ux,) vx,) = m+ x), 2.9) F x S) = { RI Fx,) S}. 2.0) For any real number X it follows from.3) that the angle subtended by S at X is θx,s) = { π if X S 0 if X / S. 2.) Hence for X = Fx,) = m + x ) RI we have ω,s,m + x) := θs,fx,)) = χf x S)) 2.2) π 3

4 where χfx S)) is the characteristic function of the set Fx S) on the real -axis. Hence, ω,s,m + x)d = χfx S))d = µ Fx S) ), 2.3) where µ Fx S)) denotes the Lebesgue measure of F x S). We then have the following theorem: Theorem 2. i) For any Borel set S and any finite interval we have lim µ Fx S) ) = ω,s,m + )d 2.4) where ω,s,m + ) is defined as in.5) for the boundary value 2.3) of the Titchmarsh-Weyl m-function for the half line problem.)-.2). ii) Let S = α),β)), where α) and β) are bounded, Lebesgue measurable functions on the interval. Then we have lim µ Fx S ) ) = ω,s,m + )d. 2.5) Proof of i): Replacing the left-hand side of 2.4) by the left-hand side of 2.3) we observe that µ Fx S)) d ω,s,m + )d = ω,s,m + x)d ω,s,m + )d ω,s,m + x ) ω,s,m x+iǫ))d + ω,s,m x +iǫ) ω,s,m+iǫ))d + ω,s,m x +iǫ) ω,s,m + ))d. Since m x z) and mz) are Herglotz functions by virtue of being Titchmarsh-Weyl m-functions for the problems.)-.2) and 2.4)-2.6), Theorem applies to give the limit of the first and third terms to be zero as ǫ 0. For the second term we have for the integrand for fixed ǫ > 0), as x, that: ω,s,m x +iǫ)) ω,s,m+iǫ)) [ ) )] = Im Im dt S t m x +iǫ) t m+iǫ) ) m+iǫ)+m x +iǫ) = Im dt S t m x +iǫ))t m+iǫ)) r x +iǫ) S t m x +iǫ))t m+iǫ)) dt ) r x +iǫ) 2 S t m+iǫ)) dt + ) K ǫ max r x+iǫ), where K ǫ is independent of, and r x +iǫ) is the radius of the circle C x,+iǫ in the m-plane which contracts to the limit point m+iǫ). Hence, for fixed ǫ > 0, lim ω,s,m x+iǫ)) = ω,s,m+iǫ)) 4

5 uniformly for. It follows that the above three terms can be made arbitrarily small by fixing ǫ > 0 sufficiently small and then choosing x sufficiently large. Proof of ii). First we assume that S = α),β)), α) < β) for all, where α) and β) are step functions on, each taking only finitely many values. Letting {α j },j =,...,N, and {β k },k =,...,M, denote the values which α) and β) may have for, we define C j,k = { α) = α j and β) = β k }. Then the C j,k are disjoint sets and C j,k = α j < β k. Putting S j,k := α j,β k ) whenever α j < β k, we define for x a, ), Fx S j,k ) = { RI Fx,) S j,k }. Then the statement 2.4) from parti) applies with S = S j,k, that is, lim µ F x S j,k) ) = ω,s j,k,m + )d. Since the C j,k are disjoint sets on the -axis, it follows that µ Fx S ) ) = µ C j,k Fx S j,k) ) j,k ω,s j,k,m + )d j,k = j,k = = βk α j β) α) [ Im [ Im t m + ) )] dt d )] t m + dt d ) ω,s,m + )d, as x, which proves the statement 2.5) when S = α),β)) and α) and β) are simple functions on. The case where α) and β) are general bounded Lebesgue measurable functions can then be handled by approximating α) and β) above and below by simple functions. Similarly to the analysis in Fulton, Pearson and Pruess [7] and Al-Naggar and Pearson [], we now introduce the space of quadratic forms of solutions of equation.): For RI, let R = { Rx,) = au 2 +buv +cv 2 a,b,c RI } 2.6) be the vector space of all quadratic forms of solutions of equation.). If y = c u + d v and y 2 = c 2 u+d 2 v, it is readily verified that y 2, y y 2, y2 2 also belong to R. Hence R is a 3-dimensional space with {u 2,uv,v 2 } as a basis. Also associated with R is the solution space of the first order Appell system see [7, 2]) of equations du dx = d P Q = 0 q 0 P q) Q, 2.7) dx R 0 0 R which was used extensively in [5, 6, 7] and turns out to be a useful companion system to the Schrödinger equation.) in the study of cases involving some a.c. spectrum. One choice of basis for the solution space of 2.7) is [U,U 2,U 3 ] = u ) 2 u v v ) 2 2uu [u v +uv ] 2vv. 2.8) u 2 uv v 2 5

6 Also on the solution space of 2.7) we can define an indefinite inner product by U,Ũ := 2P R+ PR) Q Q = const, independent of x [a, ), 2.9) and the above solutions {U,U 2,U 3 } are orthogonal with respect to this inner product. If then this inner product can also be expressed as Ux,) = a)u x,)+b)u 2 x,)+c)u 3 x,), Ũx,) = ã)u x,)+ b)u 2 x,)+ c)u 3 x,), U,Ũ = 2a c+cã) b b. 2.20) A solution U, or equivalently, its third component R is said to be normalized with respect to this indefinite inner product compare [7, Equa 4.0] and [, Equa 7]) if U,U := 4PR Q 2 = 2RR ) 2 4q )R 2 = 4ac b 2 = ) We call 4ac b 2 the discriminant of the quadratic form R = au 2 +buv+cv 2. Also, it is readily shown from 2.7) that the third component of R of any solution U of 2.7) must satisfy the third order linear equation R +4 q)r 2q R = ) The solution space of 2.22) is the quadratic form space R and so there is a one-to-one correspondence between the solution space of 2.22) and 2.7). Indeed, for every element, Rx, ), of the quadratic form space R we can generate from the Appell equations 2.7) unique functions Px,) and Qx,) such that P,Q,R) T is a solution of 2.7). We now cite the following lemma see Al-Naggar and Pearson [3, Lemma and Theorem 2]) which gives a sufficient condition for existence of a.c. spectrum in a given interval I, ): Definition. The Sturm-Liouville equation.) satisfies Condition A for a given real value of if and only if there exists a complex-valued solution yx, ) of.) for which lim N N 0 N 0 yx,) 2 dx yx,) 2 dx = ) Lemma Al-Naggar and Pearson). Let I, ) be an interval on which Condition A holds for the general equation.), and let the fundamental system {u,),v,)} be defined by the initial conditions 2.) at x = a. Then i) There exists a complex valued function M) on I which is uniquely defined by the properties: N ux,)+m)vx,)) 2 dx. x a) Im[M)] > 0 and b) lim 0 N N = ) ux,)+m)vx,) 2 dx x 0 ii) For I the function M in i) is the boundary value of the Titchmarsh-Weyl m-fuction defined by 2.3), that is, M) = lim[m+iǫ)] = A)+iB). 2.25) ǫ 0 iii)i σ ac, whereσ ac is the absolutelycontinuousspectrum ofthe Sturm-Liouville problem.)-.2), and on I the spectral density function is given by f) := π lim Im[m+iǫ)]) = Im[m+ )] = B) ǫ 0 π π. 2.26) 6

7 In accordance with the above lemma, we can formulate the assumption iii) of a.c. spectrum as follows: iii) Let I, ) be an interval on which B) = Im[m + )] > 0. Under this assumption, we now introduce a polar coordinaterepresentationof u+m + v = u+av)+ibv, namely u+av = r 0 x,)cosθ 0 x,)) 2.27) Bv = r 0 x,)sinθ 0 x,)) so that θ 0 x,) may be defined for I and x [a, ) by ) ux,)+a)vx,) cotθ 0 x,)) =. 2.28) B)vx, ) Taking the inverse of the cotangent, we have θ 0 x,) = cot ux,)+a)vx,) B)vx, ) ) 2.29) from which we deduce that θ 0 x,) = B u+av) 2 +B 2 v 2 > ) for all I and all x [a, ). This suggests that we look at the quadratic form defined by R 0 x,) := u+m+ )v 2 Im m + ) = u+a+ib)v 2 B = u+av)2 +B 2 v 2 B = B u2 + 2A B uv + A2 +B 2 v 2. B 2.3) From 2.30) and 2.3) we see that the monotone increasing function θ 0 x,), x, is related to R 0 x,) by x x θ 0 x,) = R 0 t,) dt = B a u+av) 2 dt. 2.32) +B 2 v2 a To keep θ 0 x,) lying in [0,π] for all x [a, ), we restrict it by setting it back to zero whenever x passes through a zero of vx,), that is, we set θ 0 x,) := θ 0 x,) modπ := θ 0 x,) nπ for x x n,x n+ ), 2.33) where x 0 = a, x n = x n ) = nth zero of v. We can now prove the following theorem concerning the asymptotic behavior of the sequence of functions of, { θ 0 x,)} as x. Theorem 3. Let I be an interval on which B) > 0. Let I be a bounded interval and assume C) and D) are measurable functions on satisfying for all 0 C) < D) π. Then the set of functions of, { θ 0 x, ) : I [0,π] x a, )} has the asymptotic property as x lim µ { C) < θ 0 x,) < D)} = {D) C)}d. 2.34) π 7

8 Proof. For y = cotβ with β 0,π), we define βy) = cot y) mapping, ) onto 0,π). Since cotβ is monotonically descreasing from β = 0 to β = π, it follows that for x x n,x n+ ), or or Hence, cotc)) > cot θ 0 x,)) = u+av Bv > cotd)), B)cotC)) A) > ux,) vx,) > B)cotD)) A), A) B)cotC)) < Fx,) = ux,) vx,) < A) B)cotD)). { µ RI C) < θ }) 0 x,) < D) = µ { Fx,) S }) where S := A) B)cotC)),A) B)cotD))), so it follows from Theorem 2ii) that { lim µ RI C) < θ }) 0 x,) < D) = ω,s,m + )d. But, making use of the definitions.3) and.5) we observe that ω,s,m + ) = π θs,m + ) = ) Im dt π S t A+iB) = B π t A) 2 +B 2 dt, and making the change of variable t = A Bcotθ, we obtain ω,s,m + ) = π S D) C) = π D) C). B 2 csc 2 θ B 2 cot 2 θ +) dt Substitution of this into the integral of ω,s,m + ) thus yields the result 2.34). Theorem 3 gives rise to the concept of a sequence of uniformly asymptotically distributed functions modulo π. Namely, when g n ) is a sequence of functions defined on I σ a c and mapping into [0,π] such as θ 0 x,)) which satisfy, for any two Lebesgue measurable functions C),D) with 0 C) < D) π), and any interval I, the property 2.34), that is, lim µ { : C) < g n) < D)} = {D) C)}d, n π we will say that g n ) is uniformly asymptotically distributed modulo π. Further investigations of such functions and how they give rise to formulas for spectral density functions defined as in 2.26) are currently in progress. 3 An example: Bessel s equation In this section we give an example for which the spectral density function is known, and show that it can be obtained from the formula in [7, Thm 4.2]. Consider the Bessel equation of order ν [0, ) on [a, ), a > 0, 8

9 ν 2 ) /4 y + y = y, a x < 3.) x 2 ya) = 0. Since this problem has a.c. spectrum in 0, ) we may take I = 0, ) in Theorem 3. The solutions of this equation defined by the initial conditions 2.) at x = a are compare [4, p. 86]), where ux,) = D xjν x )+D 2 xyν x ) 3.2) vx,) = D 3 xjν x )+D 4 xyν x ), 3.3) D = π 2 a 2 Yν a ), D 2 = π 2 a 2 Jν a ), D 3 = π 4 a 2 Yν a ) + π 2 a 2 Y ν a ), D 4 = π 4 a 2 Jν a ) π 2 a 2 J ν a ). Using the fundamental system of solutions of the Appell system 2.8), we can find the unique solution, for 0, ), of the initial value problem at infinity see [7, Equa.6)-.7)] or [5, Thm ]), lim [ãu + bu 2 + cu 3 ] = 0, 3.4) by employing asymptotics for the Bessel functions J ν and Y ν, as x, and this yields the formulas ã) = πa 2 π [ b) = Jν 2 2 c) = π 8a + π 2 a ) +Yν 2 a )), a ) a )] +πa [J ν a ) J 2 ν [ Jν 2 a ) +Y 2 ν +Y 2 ν [J ν a ) J ν a ) a )] + πa 2 +Y ν a )Y ν J ν a ) [ J νa ) 2 + Y νa ) ] 2 ) a )]. From [7, Thm 4.2, Equa 4.29] it follows that in agreement with [4, p. 86]), f) = πã) = 2 π 2 a [ Jν 2 a ) +Y 2 ν +Y ν a )Y ν a )]. a )] 3.5) From [7, Equas 4.22 & 4.26] we also have A) + ib) = b 2ã + i ã. Hence putting these formulas, together with the solutions {u, v} from 3.2)-3.3) into 2.29), we have produced, in the definition 2.33) of θ 0 x, ), a concrete example of a uniformly asymptotically distributed sequence of functions modulo π. Acknowledgments The authors thank Jeremy Mandelkern for carrying out the calculations yielding the above formulas for the connection coefficients ã, b, c. 9

10 References [] S.V. Breimesser and D.B. Pearson,2003) Geometrical Aspects of spectral theory and Value Distribution for Herglotz functions. Math Physics, Anal and Geom, 6, [2] S. Pruess and C. Fulton, 993) Mathematical software for Sturm-Liouville problems. ACM Trans. Math. Software, 9, [3] C. Fulton and S. Pruess, The computation of spectral density functions for singular Sturm-Liouville problems involving simple continuous spectra. 998) ACM Trans. Math. Software, 34, [4] J.D. Pryce, 993) Numerical Solution of Sturm-Liouville Problems, Oxford, Clarendon Press. [5] C. Fulton, D. Pearson, and S. Pruess, 2008) New characterizations of spectral density functions for singular Sturm-Liouville problems. J. Comp. Appl. Math, 22, [6] C. Fulton, D. Pearson, and S. Pruess, 2008) Efficient calculation of spectral density functions for specific classes of singular Sturm-Liouville problems. J. Comp. Appl. Math, 22, [7] C. Fulton, D. Pearson, and S. Pruess, 204) Estimating spectral density functions for Sturm- Liouville problems with two singular endpoints. IMA Jour. of Numer. Anal ), [8] D.B. Pearson, 993) Value Distribution and spectral analysis of differential operators. J. Phys. A: Math Gen, 26, [9] D.B. Pearson,994) Value distribution and spectral theory. Proc. London Math. Soc, 68 3), [0] S.V. Breimesser and D.B. Pearson,2000) Asymptotic value distribution for solutions of the Schrödinger Equation. Math Physics, Anal and Geom, 3, [] I. Al-Naggar and D.B. Pearson, 996) Quadratic forms and solutions of the Schrödinger equation. J. Phys A: Math Gen, 29, [2] M. Appell, 880) Sur la transformation des équations différentielles linéaires. Comptes rendus hebdomadaires des séances de l Académie des sciences, 9 4), [3] I. Al-Naggar and D.B. Pearson, 994) A new asymptotic condition for absolutely continuous spectrum of the Sturm-Liouville operator on the half-line. Helvetica Physica Acta 67, [4] Titchmarsh, E.C., 962) Eigenfunction Expansions associated with Second-order Differential Equations, 2nd edn., Oxford, Clarendon Press. 0