arxiv: v2 [math.ca] 30 May 2011

Transcription

1 Orthogonal polynomials associated to a certain fourth order differential equation arxiv:97.61v [math.ca] 3 May 11 Joachim Hilgert, Toshiyuki Kobayashi, Gen Mano, Jan Möllers Abstract We introduce orthogonal polynomials (x as eigenfunctions of a certain self-adoint fourth order differential operator depending on two parameters µ C and l N. These polynomials arise as K-finitevectors inthel -model oftheminimal unitary representations of indefinite orthogonal groups, and reduce to the classical Laguerre polynomials L µ (x for l. We establish various recurrence relations and integral representations for our polynomials, as well as a closed formula for the L -norm. Further we show that they are uniquely determined as polynomial eigenfunctions. Mathematics Subect Classification: Primary 33C45; Secondary E46, 34A5, 4C15. Key words and phrases: orthogonal polynomials, generating functions, Bessel functions, Laguerre polynomials, fourth order, recurrence formulas, minimal representation, Meier s G-function. Contents 1 Introduction 1 Statement of the main theorems 3 3 Proofs of the main theorems 6 4 Comparison with Laguerre polynomials 1 References 14 1 Introduction Many classical sequences of special polynomials (P N such as Hermite polynomials, Laguerre polynomials or Jacobi polynomials have the following standard properties: P (x is a polynomial of degree (,1,,..., P (x is an eigenfunction of a second order differential operator. Partially supported by Grant-in-Aid for Scientific Research (B (183437, 346, Japan Society for the Promotion of Science, and the Alexander Humboldt Foundation. Partially supported by the International Research Training Group 1133 Geometry and Analysis of Symmetries, and the GCOE program of the University of Tokyo. 1

2 In this article we introduce a family of polynomials ( N depending on two parameters µ C and l N. A distinguishing feature in our setting is that these polynomials have the following properties: (x is a polynomial of degree +l (,1,,..., (x is an eigenfunction of a fourth order differential operator (Theorem.3. It turns out that the polynomials ( N enoy a number of good properties one typically finds for classical orthogonal polynomials: They are unique as polynomial solutions of differential equations (Theorem.3, form orthonormal bases for L (R +,x µ l e x dx (Theorem.4, admit integral representations (Theorem.7, satisfy recurrence relations (Theorem.5. Ouranalysisofthesespecialpolynomials is motivated by recent progress on the minimal representation of the non-compact semisimple Lie group G O(p,q (the indefinite orthogonal group on the Hilbert space L (C where C is an isotropic cone in R p+q (see [6, 7]. In contrast to traditional analysis on homogeneous spaces, the group G in our setting is too large to act geometrically on C. This very feature lets one expect many functional equations in the representation space arising from the action of the large group G on the small representation space. In fact, some of the formulas here are predicted by unitary representation theory with (p,q (µ+3,l+4 where µ is an odd integer. The polynomials arise as K-finite vectors in the representation, and indeed generate all the K-types parametrized by non-negative integers. In the bottom case l our polynomials M µ, reduce to Laguerre polynomials L µ and the fourth order differential operator P µ,l is of the form P µ, Q µ +const with the Laguerre operator Q µ having L µ as eigenfunctions. As its remarkable consequence, Laguerre polynomials (multiplied by elementary functions give a basis of K-finite vectors for the minimal representation of O(m, 4. Similar results were previously known only for the conformal group O(p, (in physics terms the minimal representation of O(4, appears as the bound states of the Hydrogen atom, and incidentally as the quantum Kepler problem. Our proofs of the main results rely largely on purely analytic methods, in particular the analysis of the corresponding fourth order differential equation, which we started in [5]. The anonymous referee pointed out that in the 194s H. L. Krall [8] introduced families of orthogonal polynomials associated to certain fourth order differential operators. Later A. M. Krall and L. L. Littleohn took up the study (see e.g. [9] for a survey. While it is not immediate that Krall s Laguerre type polynomials are in any way related to our new polynomials, this issue seems to deserve closer study. Notation: N {,1,,...}, R + {x R : x > }.

3 Statement of the main theorems We introduce the family of polynomials { (x : N } for µ C and l N via meromorphic generating functions G µ,l (t,x of two variables t and x defined by G µ,l (t,x : 1 (1 t l+µ+3 ( x ( ( l+1e x tx Ĩ K µ (1 t l+ 1 x (1 t, (.1 where Ĩα(z ( z α I α (z and K α (z ( z α K α (z denote the normalized I- and K-Bessel functions. More precisely, we set Γ( +µ+1 (x :! µ Γ( + µ+1 t G µ,l (t,x. (. t The coefficient in front is chosen to produce suitable normalizations of the top term. Theorem.1. Suppose µ 1,, 3,... and l N. (1 (x is a polynomial of degree +l (,1,,... ( (Top term (3 (Constant term (x ( 1 x +l +lower order terms.! ( l µ Γ(l+ 1!Γ( µ+ ( Γ( +µ+1 µ+1 l Γ( + µ+1, where (a n a(a+1 (a+n 1 is the Pochhammer symbol. A combinatorial formula for intermediate terms of (x will be given in Proposition 3.1. Here are some further special values of the polynomials (x. Example.. (1 (l The polynomials (x for l reduce to the Laguerre polynomials M µ, (x L µ (x. ( ( The bottom of the series with amounts to (x l k (l k! k!(l k! xk. These polynomials appear in the explicit formula for the K-Bessel functions with half-integer parameter. In fact, for l N and any µ we have (see e.g. [1, III.71 (1] K l+ 1 (z πz (l+1 e z (z. 3

4 To state the differential equationfor the polynomials (x we set θ : d dx and introduce the fourth order differential operator P µ,l : 1 ( ( x θ+µ l 1 x ( θ +µ x on R +. ( x ( ( θ l 1 x ( θ x ( x Theorem.3 (Differential equation. Let µ 1,, 3,... and l N. For every N the polynomial (x is a solution of the fourth order differential equation P µ,l u ( +µ+1u. (.3 is the unique polyno- Moreover, if µ l+1, then, up to scalar multiple, mial solution of this equation. Theorem.4 (Orthonormal basis. For µ > l 1 we have (x k (xxµ l e x dx Γ( +µ+1γ( +l+ µ+3 µ+1 Γ( l+!( +µ+1γ( + µ+1 ( k, ( k. If further µ l+1 is an odd integer, then the sequence ( N forms an orthogonal basis of L (R +,x µ l e x dx. Theorem.5 (Recurrence relations. Suppose µ 1,, 3,... and l N. Then the polynomials are subect to the following recurrence relations: (1 The three-term recurrence relation for (x d dx xmµ,l (x: (θ x (x ( (x (µ l+1mµ,l (x ( +µ+l+1( +µ l 1 ( +µ 1 ( +µ+1( +µ 1 (x. ( The five-term recurrence relation for x (x: x (x k a µ,l,k Mµ,l +k (x 4

5 with coefficients a µ,l, ( +1( +, a µ,l,1 ( +1( +µ+, a µ,l, (6 +6(µ+1 +(µ+1(µ+ a µ,l, 1 4( +(µ+1 +(µ 1(µ+ l(l+1, ( +µ 1( +µ+3 +µ( +µ( +µ+l+1( +µ l 1 (, ( +µ 1( +µ+1 a µ,l ( +µ 1( +µ( +µ+l 1( +µ l 3, ( +µ 3( +µ 1 (3 The recurrence relation in µ: µ( +µ 1 (x µ( +µ 1 (x (4 The recurrence relation in l (l 1: ( +µ+l+1( +µ l 1. ( +µ+1 ( +µ 1( +µm µ,l (x x M µ+,l (x. ν( +µ 1 (x (l+1( +µ 1 (x 1 ( +µ 1Mµ,l+1 (x ( +µ 1 (5 The recurrence relation in µ and l: ( x M µ,l 1 (x. ( θ l 1 x ( ( +µ 1 (x 4( +µ 1 (x x M µ+,l (x ( +µ 1+1 (x. Remark.6. The recurrence relation in Theorem.5(5 is useful for the actual computation of the polynomials +1 for fixed l from the polynomials. In fact, since in the bottom case l the polynomials are simply the Laguerre polynomials, Theorem.5(5 gives us an inductive method to calculate the series of special orthogonal polynomials. For example, for l and µ 1,, 3,... we find M µ,1 (x x +µ 1 Lµ+ (x ( θ 1 x ( L µ 4( +µ (x +µ 1 Lµ 1 (x Next we turn to integral representations of the polynomials in terms of Laguerre polynomials. Note that for these integrals it is a priory not clear that they are polynomial in x.. 5

6 Theorem.7 (Integral representation. For Re µ > 1 and x > the integral x l+1 π e x(coshϕ 1 L l+µ+1 (x(cosθ +coshϕsin µ θsinh l+1 ϕdϕdθ. becomes a polynomial in x of degree +l. Further, it is equal to µ l!γ( µ+1 µ+1 Γ( + Γ( +µ+1 (x. If µ l+1 is an odd integer, the polynomials satisfy a remarkable reproducing property with respect to Meier s G-transform which is built from Meier s G-function G 4 (z b 1,b,b 3,b 4 (see [5, Section 7]. Theorem.8 (Meier s G-transform. If µ l + 1 is an odd integer and x >, then ( G 4 (xy, l 1, µ µ+1, l e y y µ+1 (4y dy 3 Proofs of the main theorems ( 1 e x x (l+1 (4x. The proof of Theorem.1 will be obtained from an explicit calculation which leadsusto acombinatorialexpressionofthe functions in terms of Laguerre polynomials. Then the strategy is to relate the polynomials to a series of functions which was already studied thoroughly in [5]. This will help to prove Theorem.3 and reduce the remaining theorems to results from [5]. We start with a proof of the fact that the functions are indeed polynomials. Proposition 3.1. Suppose µ 1,, 3,... and l N. Then (x (,1,,... is given by Γ( +µ+1 (x Γ( + µ+1 k l k k i ( 1 k Γ( k + µ+1 (l i! k!γ( k +µ+1(l i k!i! Lµ k (xxi +l β µ,l,k xk, (3.1 where L α n denotes the Laguerre polynomials and with β µ,l,k Γ( +µ+1 Γ( + µ+1 (m,n S µ,l,k µ+1 ( 1 m+nγ( m+ Γ(n+µ+1 S µ,l,k { (m,n N : (l+n k! m!n!(k n!( m n!(l+n k m! n m k n l m }. (3. 6

7 Proof. Let us first assume Re µ > 1. The K-Bessel functions with halfinteger parameter can be written explicitly as (see e.g. [1, III.71 (1] K l+ 1 (z πz (l+1 e z l i (l i! (l i! i! (zi. (3.3 Using the following integral representation for the I-Bessel function (cf. [1, III.71 (9] we obtain Ĩ α (z G µ,l (t,x 1 πγ(α+ 1 1 (1 t l+µ+3 π 1 Γ( µ+1 π ( x e xcosϕ sin α ϕdϕ, Re α > 1, ( ( l+1e x tx x Ĩ K µ (1 t l+ 1 (1 t 1 tx (1 t µ 1 +1e (1 t (cosϕ+1 sin µ ϕdϕ l i (l i! (l i! i! xi (1 t l i. Next, we compute the derivatives of the first factor with respect to t at t. Using the formula of the generating function for the Laguerre polynomials (see e.g. [1, (6..4] L α 1 n (xtn tx 1 t (3.4 (1 t α+1e we find that! t π t L µ 1 n [ π and substituting y 1 (cosϕ+1 yields! µ 1 1 tx (1 t µ 1 (1 t (cosϕ+1 sin µ ϕdϕ +1e ( x (cosϕ+1 sin µ ϕdϕ (1 y µ 1 y µ 1!µ Γ( + µ+1 Γ(µ+1 Γ( +µ+1 µ 1 L L µ (x, (xydy where the last equality is the integral formula[4, 16.6(5]. Now we can compute ] 7

8 the Taylor coefficients of G µ,l (t,x at t explicitly as follows t G µ,l (t,x t 1 ( [ k ] π 1 tx Γ( µ+1 k t k k t (1 t µ 1 (1 t (cosϕ+1 sin µ ϕdϕ +1e [ k l ] (l i! t k (l i! i! xi (1 t l i k! µ Γ( k + µ+1 k!γ( k +µ+1 Lµ k (x l k k i l i t i (l i! (l i! i! xi ( 1 k (l i (l i k +1 ( 1 k! µ Γ( k + µ+1 (l i! k!γ( k +µ+1(l i k!i! Lµ k (xxi. This gives the first expression for (x. Inserting the explicit formula (cf. [1, (6..] L α n(x Γ(n+α+1 n! n ( n ( 1 k x k k Γ(k +α+1 k one obtains the expressions for the coefficients β µ,l,k given in the proposition. Since these clearly have meromorphic continuation for µ C with poles at most at µ 1,, 3,..., the claim follows. Proof of Theorem.1. It remains to compute top and bottom term of the polynomials (x. For k +l the set S µ,l,k defined in (3. only contains the tuple (, and we obtain the top term β µ,l,+l ( 1.! Tocalculate the bottom term ( simply observe that (see e.g. [1, Chapter III.7] so that Ĩ µ ( 1 Γ( µ+ and ( x l+1 Kl+ 1 G µ,l (t, l+1 Γ(l+ 1 Γ( µ+ (1 t l µ+1 Together with (. this proves the claim. (x x Γ(l+ 1 l Γ(l+ 1 (µ+1 l!γ( µ+ t. In order to be able to apply results from [5] in our context we need two observations relating the polynomials (x and the differential operators P µ,l to corresponding obects introduced in [5]., 8

9 Lemma 3.. In the notation of [5] we have, (x µ Γ( + µ+1 Γ( +µ+1 x (l+1 e x (x. Λ µ,l+1 Proof. Directly from [5, equations (3. and (4.]. Next, we relate the operator P µ,l to the fourth order differential operator D µ,ν 1 x ( (θ +ν(θ +µ+ν x ( θ(θ +µ x introduced in [5]. Lemma 3.3. For f C (R + we have (µ ν(µ+ν + (3.5 D µ,l+1 (x (l+1 e x f(x (( x (l+1 e x 4P µ,l + (µ l 1(µ+l+3 f (x. Proof. Since D µ,ν is symmetric in µ and ν by [5, Proposition.1 (1], we can rewrite the expression (3.5 for D µ,ν as D µ,ν 1 x ( (θ +µ(θ +µ+ν x ( θ(θ +ν x + (µ ν(µ+ν +. Then for ν l+1 the claim follows by using the commutator relation [ θ,x ν e x] (ν +xx ν e x. Now, in view of Lemmas 3. and 3.3, it follows from [5, Theorem 4.6] that (x solves the differential equation (.3. The crucial observation for the proof of uniqueness in Theorem.3 is that the subspace of solutions of (.3 bounded near x can be wholly described in terms of the polynomials (x. Note, that the operator P µ,l is invariant under the transformation f(x f(x e x f( x, i.e. P µ,l f P µ,l f. Then clearly N µ,l µ,l (x : M (x e x ( x is another non-trivial solution of (.3. Since N µ,l (x grows exponentially as x and (x ust polynomially, the two solutions are linearly independent. We even have the following lemma: Lemma 3.4. Suppose µ l + 1. Then the subspace of solutions of (.3 which are bounded near x is two-dimensional and spanned by the functions (x and N µ,l (x. Proof. By the previous considerations it remains to show that the subspace of solutions of (.3 which are bounded near x is two-dimensional. We note first that the differential operator P µ,l has a regular singularity at x with characteristic exponents {, µ,l+1,l+1 µ}. In fact, an easy computation shows that x P µ,l (θ +µ l 1(θ +µ(θ l 1θ (mod x C[x,θ], 9

10 where C[x,θ] denotes the left C[x]-module generated by 1,θ,θ,... in the Weyl d algebra C[x, dx ]. Therefore, the differential equation P µ,lu λu is of regular singularity at x, and its characteristic equation is given by (s+µ l 1(s+µ(s l 1s. Since l + 1 > l+1 µ > µ for µ l+1, the theory of regular singularities (see e.g. [, Chapter 4] assures that the subspace of solutions of (.3 bounded near x is two-dimensional. In fact, for µ > l+1 there are exactly two non-negative exponents, namely l + 1 and. Thus, the subspace of solutions which are bounded near x is spanned by functions f 1 (x and f (x with f 1 (x x l+1 and f (x x 1 for x. If µ l+1, then there are three non-negative exponents, but the exponent has multiplicity two. Hence, the subspace of solutions with asymptotic behavior at x given by these exponents is three-dimensional and spanned by functions f 1 (x, f (x and f 3 (x with asymptotic behaviorat x given by f 1 (x x l+1, f (x x 1 and f 3 (x log(x. But f 3 (x is clearly not bounded near x, so again the subspace of solutions bounded near x is two-dimensional. Proof of Theorem.3. Only the uniqueness remains to be shown. To do that, observe that by Lemma 3.4 every polynomial solution of (.3 has to be a linear combination of (x and N µ,l (x. But N µ,l grows exponentially as x and hence cannot be a polynomial. This leaves as the only polynomial solution of (.3 (up to scalar multiples. With these preparations the remaining results of Section are immediate consequences of results from [5]: Proof of Theorem.4. This follows from [5, Corollaries 4.8 and 6.]. Proof of Theorem.5. In view of Lemma 3. the recurrence relations follow from [5, Propositions 6.1, 6.4 and 6.6]. Proof of Theorem.7. This is essentially a reformulation of [5, Theorem 5.1 (1]. Proof of Theorem.8. The stated formula follows immediately from [5, Theorem 7.]. 4 Comparison with Laguerre polynomials We saw in Example. that the polynomials are Laguerre polynomials if l. Thus our results on the specialize to results on Laguerre polynomials. In this section we examine how these results are related to standard results on Laguerre polynomials. We start by collecting a number of formulas for Laguerre polynomials (cf. [1, Chapter 6.] and [3, II, Chapter 1.1]. The polynomial L α n (x is defined by n L α n (x ( 1 k ( n+α x k. (4.1 k! n k k 1

11 We start with the three independent formulas xl α n (n+1l α n+1 +(n+α+1l α n (n+αl α n 1, (4. d dx Lα n d dx Lα n 1 L α n 1, (4.3 d dx Lα n L α+1 n 1. (4.4 These imply the four additional identities θl α n nlα n (n+αlα n 1 (x, (4.5 d dx +(α+1 x d dx +n L α n, (4.6 xl α+1 n (n+α+1l α n (n+1l α n+1, (4.7 L α n 1 L α n L α 1 n. (4.8 We will also need the summation formula n d n dx Lα n (x ( 1 L α+ n (x Lα 1 n (x. (4.9 k which follows from (4.4 and (4.1 by a simple calculation. Now we can examine how the various types of results from Section specialize to Laguerre polynomials. Generating function. For l the expression (.1 for the generating function G µ,l (t,x of the polynomials (x can be simplified using (see [1, 3.71 (13] π K 1 (x x e x. Then for l combining the equations (.1 and (. yields Γ( + µ+1 Γ( +µ+1 t L µ (x π µ (1 t µ+1 e tx (1 t Ĩ µ ( tx (1 t (4.1 To see that (4.1 is in agreement with the standard generating function for the Laguerre polynomials given by (3.4 one can use the following two formulas ([3, I, 6.1 (5] and [3, II, 7.. (1] n (c n (α+1 n t n L α n(x (1 t c 1F 1 Ĩ α (x e x Γ(α+1 1 F 1 where 1 F 1 (a;b;x denotes the hypergeometric function. ( c;α+1; tx 1 t, (α+ 1 ;α+1;x, Differential equation. For l the fourth order differential operator P µ,l essentially degenerates to the square of a second order operator. In fact, P µ, Q µ ( µ+1 11

12 with Q µ 1 ( θ +(µ xθ µ+1 x x d x dx +(µ+1 x d dx µ+1. But the Laguerre differential equation (4.6 shows that the polynomials M µ, L µ are eigenfunctions of Q µ for the eigenvalue ( + µ+1. The fourth order differential equation (.3 follows from this by applying Q µ twice. Orthonormal basis. In the special case l, we have (x L µ (x (see Example. (1, and therefore Theorem.4 reduces to the well-known fact that the Laguerre polynomials (L µ N form a complete orthogonal system of L (R +,x µ e x dx with norms L µ L (R +,x µ e x dx Γ( +µ+1. (4.11! Ourresultinthis specialcaseisslightlyweakerin the sensethat weassumedµ Z+1 for the completeness of the sequence ( N in L (R +,x µ l e x dx whereas this remains true for the Laguerre polynomials for arbitrary µ > 1 (see e.g. [1, Chapter 6]. Recurrence relations. We examine three of the different recurrence relations given in Theorem.5. (1 For l the three-term recurrence relation of Theorem.5 (1 simplifies to (θ xl µ ( +1Lµ +1 (µ+1lµ ( +µlµ 1. One can use the identities (4. and (4.5 to give an independent proof of this equation. ( In the case that l the five-term recurrence relation of Theorem.5 ( arises from the three-term recurrence relation (4.. The existence of this three term recurrence relation is predicted by the general theory of orthogonal polynomials. In fact, every sequence of orthogonal polynomials startingin degreeand increasingdegree by 1 in everystep satisfies threeterm recurrence relations for the multiplication by x (see e.g. [1, Chapter 5.]. The polynomial, however, is of degree + l and hence the degrees of the polynomials in this sequence start with only in the case where l. Already for l 1 it is possible to prove that there are no three-term recurrence relations for the multiplication by x. Indeed, if a three-term recurrence relation existed, there would be constants a, b C such that xm µ,1 (x am µ,1 1 (x+bm µ,1 (x. (4.1 Since in this case the polynomials have the specific form M µ,1 (x x+, M µ,1 1 (x x +(µ 1x+(µ 1, the equation (4.1 is equivalent to x(x+ a( x +(µ 1x+(µ 1+b(x+ which cannot hold for all x. 1

13 (3 The formula in Theorem.5 (3 for l reduces to x L µ+ ( +µ 1( +µlµ µ( +µ 1L µ +µ( +µlµ 1. This identity also follows by applying (4.7 twice and (4.8 three times (in this order: [ ] x L µ+ x ( +µl µ+1 ( 1Lµ+1 1 ( +µ 1( +µl µ ( 1( +µlµ 1 +( 1Lµ ( +µ 1( +µl µ 1 1 +( +µ+1( +µlµ 1 +( 1Lµ ( +µ 1( +µl µ ( +µ 1( +µl µ 1 +( +µ+1( +µl µ 1 +( 1Lµ ( +µ 1( +µl µ µ( +µ 1L µ +µ( +µlµ 1. Integral representation. In the case where l the integral representation of Theorem.7 amounts to L µ (x Γ( +µ+1 µ Γ( µ+1 µ+1 Γ( + π ( e x(coshϕ 1 L µ+1 (x(cosθ +coshϕxsinhϕdϕ sin µ θdθ. Iterated integration by parts for the inner integral, using yields L µ (x Γ( +µ+1 µ Γ( µ+1 µ+1 Γ( + [ ϕ e x(coshϕ 1 xe x(coshϕ 1 sinhϕ, π L µ d dx µ+1 L ] (x(cosθ+1sin µ θdθ. Using the summation formula (4.9 and the substitution y : cosθ+1 we finally obtain the integral formula [4, 16.6 (5] L µ (x Γ( +µ+1 µ Γ( µ+1 µ+1 Γ( + Γ( +µ+1 µ Γ( µ+1 µ+1 Γ( + π 1 L µ 1 (x(cosθ +1sin µ θdθ L µ 1 (xyy µ 1 µ 1 (1 y dy. Meier s G-transform. For l the G-function appearing in Theorem.8 reduces to a J-Bessel function, namely (cf. [3, I, 5.6 (11] ( G 4 t, 1, µ, µ+1 t µ 4 Jµ (4t 1 4, t >. 13

14 Then (after a suitable substitution the reproducing property of Theorem.8 can be rewritten as J µ (xy(xy 1 y µ+ 1 e 1 y L µ (y dy ( 1 x µ+1 e 1 x L µ (x, which is the well known formula for the Hankel transform of Laguerre polynomials (see e.g. [4, 8.9 (3]. References [1] G. E. Andrews, R. Askey, and R. Roy, Special functions, Encyclopedia of Mathematics and its Applications, vol. 71, Cambridge University Press, Cambridge, [] E. A. Coddington and N. Levinson, Theory of ordinary differential equations, McGraw-Hill Book Company, Inc., New York, [3] A. Erdélyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Higher transcendental functions. Vols. I, II, McGraw-Hill Book Company, Inc., New York, [4], Tables of integral transforms. Vol. II, McGraw-Hill Book Company, Inc., New York, [5] J. Hilgert, T. Kobayashi, G. Mano, and J. Möllers, Special functions associated to a certain fourth order differential equation, to appear in Ramanuan J. Math. (arxiv:97.68 [6] T. Kobayashi and G. Mano, Integral formula of the unitary inversion operator for the minimal representation of O(p, q, Proc. Japan Acad. Ser. A (7, [7] The Schrödinger model for the minimal representation of the indefinite orthogonal group O(p,q, the Mem. Amer. Math. Soc. (11vol. 1, no. 1 DOI: 1.19/S (available at arxiv: [8] H. L. Krall, On orthogonal polynomials satisfying a certain fourth order differential equation, Pennsylvania State College Studies (194, no. 6. [9] L. L. Littleohn and A. M. Krall, Orthogonal polynomials and higher order singular Sturm-Liouville systems, Acta Appl. Math. 17 (1989, no., [1] G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, Cambridge, England, Joachim Hilgert Institut für Mathematik, Universität Paderborn, Warburger Str. 1, 3398 Paderborn, Germany. 14

15 address: Toshiyuki Kobayashi Home address: Graduate School of Mathematical Sciences, IPMU, the University of Tokyo, Komaba, Meguro, Tokyo, , Japan. Current address: Max-Planck-Institut für Mathematik, Vivatsgasse 7, Bonn, Germany. address: Gen Mano Graduate School of Mathematical Sciences, the University of Tokyo, Komaba, Meguro, Tokyo, , Japan. Jan Möllers Institut für Mathematik, Universität Paderborn, Warburger Str. 1, 3398 Paderborn, Germany. address: 15