arxiv: v3 [math.ca] 16 Nov 2010

Transcription

1 A note on parameter derivatives of classical orthogonal polynomials Rados law Szmytowsi arxiv: v3 [math.ca] 16 Nov 010 Atomic Physics Division, Department of Atomic Physics and Luminescence, Faculty of Applied Physics and Mathematics, Gdańs University of Technology, Narutowicza 11/1, PL Gdańs, Poland Coefficients in the expansions of the form P n(λ;z) October 9, 018 Abstract a n (λ)p (λ;z), 0 where P n(λ;z) is thenthclassical (thegeneralized Laguerre, Gegenbauer or Jacobi) orthogonal polynomial of variable z andλis aparameter, are evaluated. Amethodwe adoptinthepresent paper differs from that used by Fröhlich [Integral Transforms Spec. Funct. (1994) 53] for the Jacobi polynomials and by Koepf [Integral Transforms Spec. Funct. 5 (1997) 69] for the generalized Laguerre and the Gegenbauer polynomials. KEY WORDS: orthogonal polynomials; generalized Laguerre polynomials; Gegenbauer polynomials; Jacobi polynomials; parameter derivatives AMS subject classification: 4C05, 33C45, 4C10 1 Introduction The problem of evaluation of coefficients in the expansions of the form P n (λ;z) a n (λ)p (λ;z), (1.1) 0 where P n (λ;z) is the nth classical (the generalized Laguerre, Gegenbauer or Jacobi [1 3]) orthogonal polynomial of variable z, with λ being a parameter, is not new. Fröhlich [4] investigated the case of the Jacobi polynomials, P n (λ;z) P n, with either of the two parameters α or β playing the role of λ. He started with a relationship between the Jacobi polynomials and the Gauss hypergeometric function, found partial derivatives of the latter with respect to the second and the third parameters, and then ingeniously used some combinatorial identities to obtain closed-form expressions for the coefficients in two (one corresponding to the choice λ α, the other to the choice λ β) relevant expansions of the form (1.1). Using a notation more suitable for the present purposes than that used originally in Ref. [4], Fröhlich s results may be written as (cf Ref. [5, Eqs. 1

2 (1..3.1) and (1..3.)]) P n α [ψ(n+α+β +1) ψ(n+α+β +1)]P n + Γ(n+β +1) +α+β +1 Γ(n+α+β +1) (n )( +n+α+β +1) 0 Γ( +α+β +1) P (α,β) (z) (1.) Γ( +β +1) and P n β [ψ(n+α+β +1) ψ(n+α+β +1)]P n +( ) n Γ(n+α+1) ( ) +α+β +1 Γ(n+α+β +1) (n )( +n+α+β +1) 0 Γ( +α+β +1) P (α,β) (z), (1.3) Γ( +α+1) where ψ(ζ) 1 dγ(ζ) Γ(ζ) dζ is the digamma function. Somewhat later, Koepf [6] pointed out that one might use the Fröhlich s results(1.) and(1.3) and the following nown relationship between the Gegenbauer and the Jacobi polynomials: C n (λ) (z) Γ(λ+ 1 ) Γ(n+λ) (1.5) Γ(λ) Γ(n+λ+ 1 )P(λ 1/,λ 1/) to obtain an expansion of the form (1.1) for C (λ). With some simplifying notational changes, Koepf s result is (cf Ref. [5, Eq. (1...)]) C n (λ) (z) [ψ(n+λ) ψ(λ)]c n (λ) 1+( ) +n +λ (z)+4 (n )( +n+λ) C(λ) (z), (1.6) 0 which may be also rewritten as (cf Ref. [5, Eq. (1...1)]) (1.4) C n (λ) int(n/) (z) [ψ(n+λ) ψ(λ)]c n (λ) n +λ (z)+ (n +λ) C(λ) n (z). (1.7) 1 Moreover, Koepf [6] observed that combining the Fröhlich s result (1.3) with the nown formula ( L (α) lim β P(α,β) n 1 z ), (1.8) β relating the generalized Laguerre and the Jacobi polynomials, one obtains the remarably simple expansion L (λ) 1 n L(λ) (z), (1.9) 0 which finds applications in relativistic quantum mechanics, where it is used to determine corrections of the order α (α is the fine-structure constant) to wave functions of particles bound in the Coulomb potential. In Ref. [6], Koepf gave an alternative proof of this formula.

3 To mae this brief survey complete, we have to mention that in Ref. [7] Ronveaux et al. found recurrence relations for coefficients in an expansion more general than the one in Eq. (1.1), namely m P n (λ;z) m a n (m,λ)p (λ;z) (m N), (1.10) 0 and in some particular cases they were able to obtain the coefficients a n (m,λ) in closed forms. Our interest in the matter discussed in the present paper arose and evolved in the course of carrying out research on derivatives of the associated Legendre function of the first ind, P ν µ (z), with respect to its parameters µ and ν [8 10]. In particular, we have found [9] that the derivatives (z)/ ν] νn, with M,N N, may be expressed in terms of the parameter derivatives of the Jacobi polynomials, P n / β, with suitably chosen n, α and β. Being at the time of doing the aforementioned research unaware of the Fröhlich s paper [4], we derived the relationships (1.) and (1.3) independently. Then it appeared to us that our way of reasoning might be applied to derive the relationships (1.6) and (1.9) as well. Since our approach turns out to be entirely different from the methods adopted by Fröhlich [4] and Koepf [6] (as well as from the one used by Ronveaux et al. [7]), we believe it is worthy to be presented. Thus, in the present paper, we re-derive the expansions (1.), (1.3), (1.6) and (1.9). The path we shall follow for each particular polynomial family comprises the following steps. First, using the lemma proved in the Appendix and explicit representations of the polynomials in terms of powers of z z 0, with suitably chosen z 0, for each family we shall deduce the coefficient a nn (λ) standing in Eq. (1.1) at P n (λ;z). Then, from the homogeneous differential equation satisfied by P n (λ;z), in each particular case we shall derive an inhomogeneous differential equation obeyed by P n (λ;z)/. Next, we shall show how, in each case, after exploiting a relevant Christoffel Darboux identity, the inhomogeneity in that differential equation may be cast to a form of a linear combination of the polynomials P (λ;z) with 0. In the final step, we shall insert this combination and the expansion (1.1) into the aforementioned inhomogeneous differential equation for P n (λ;z)/, and then straightforwardly deduce the coefficients a n (λ) with 0. All properties of the classical orthogonal polynomials we have exploited in the present paper may be found in Ref. [3]. [ P ±M ν The generalized Laguerre polynomials The explicit formula defining the generalized Laguerre polynomials is L (λ) 0 ( )! ( ) n+λ z. (.1) n The coefficient at z n is n (λ) ( )n (.) n! and is seen to be independent of λ. Hence, using Eq. (A.1), we deduce that the coefficient a nn (λ) in the expansion vanishes: L (λ) 0 a n (λ)l (λ) (z) (.3) a nn (λ) 0. (.4) The coefficients a n (λ) with 0 may be found through the following reasoning. The differential equation obeyed by L (λ) is [z d +(λ+1 z) d +n ] L (λ) 0, (.5) 3

4 which implies that [z d +(λ+1 z) d +n ] L (λ) dl(λ). (.6) We shall express the derivative dl (λ) / as a linear combination of the polynomials L (λ) (z) with 0. To this end, at first we invoe the nown relation dl (λ) nl(λ) (n+λ)l (λ) z (z) Next, we recall that the pertinent Christoffel Darboux identity is 0! Γ( +λ+1) L(λ) (z)l(λ) If we set here z 0 and use Eq. (.8) becomes 0 Combining Eqs. (.7) and (.10) yields n! (z ) Γ(n+λ). (.7) L n (λ) (z)l (λ) (z ) L (λ) (z)l(λ) n (z ) z z. (.8) L (λ) Γ( +λ+1) (0)!Γ(λ+1), (.9) L (λ) (n+λ)l (z) nl(λ) (λ) z dl (λ) 0 (z). (.10) L (λ) (z). (.11) Next, we insert Eqs. (.3) and (.11) into the left- and right-hand sides of Eq. (.6), respectively, and then simplify the left-hand side of the resulting equation with the aid of Eq. (.5). In this way, we obtain (n )a n (λ)l (λ) (z) L (λ) (z). (.1) 0 Equating coefficients at L (λ) (z) on both sides of Eq. (.1), we arrive at the sought expression for a n (λ): a n (λ) 1 (0 ). (.13) n 0 From Eqs. (.3), (.4) and (.13) the relationship (1.9) follows. 3 The Gegenbauer polynomials The Gegenbauer polynomials may be defined through the formula C (λ) Γ(λ+ 1 ) Γ(λ) 0 To find the coefficient a (λ) nn in the expansion C n (λ) (z) Γ( +n+λ)!(n )!Γ( +λ+ 1 ) 0 ( ) z 1. (3.1) a n (λ)c (λ) (z), (3.) 4

5 we infer from Eq. (3.1) that the coefficient at (z 1) n is n (λ) Γ(λ+ 1 ) Γ(λ) Hence, application of the lemma from the Appendix gives Γ(n+λ) n n!γ(n+λ+ 1 ) n Γ(n+λ). (3.3) n!γ(λ) a nn (λ) ψ(n+λ) ψ(λ). (3.4) To proceed further, we recall that the Gegenbauer polynomials obey the differential identity [(1 z ) d (λ+1)z d ] +n(n+λ) C n (λ) (z) 0. (3.5) Differentiating with respect to λ, we obtain [(1 z ) d (λ+1)z d +n(n+λ) The next step is a bit tricy. It consists of writing [ z dc(λ) nc n (λ) (z) Using the nown relation Eq. (3.7) may be cast to the form z dc(λ) ] C (λ) (z +1) dc(λ) nc n ]+ (λ) (z) z dc(λ) nc n (λ) (z). (3.6) [ ] (z 1) dc(λ) nc n (λ) (z). (3.7) (z 1) dc(λ) nzc n (λ) (z) (n+λ 1)C (λ) (z), (3.8) nc n (λ) nc(λ) (n+λ 1)C (λ) (z) (z) nc(λ) z 1 Now, the Christoffel Darboux identity for the Gegenbauer polynomials is 0!( +λ) Γ( +λ) C(λ) (z)c(λ) (z ) If we set in this identity z ±1 and use we obtain (±) ( +λ)c (λ) (±)n (z) 0 Combining this result with Eq. (3.9) yields the expansion +(n+λ 1)C (λ) z +1 (z). (3.9) n! C n (λ) (z)c (λ) (z ) C (λ) (z)c(λ) n (z ) Γ(n+λ 1) z z. (3.10) C (λ) +λ) (±1) (±)Γ(!Γ(λ), (3.11) z dc(λ) nc n (λ) (z) 4 ±nc (λ) (n+λ 1)C (λ) 0 1+( ) +n z 1 (z). (3.1) ( +λ)c (λ) (z). (3.13) If we plug the expansions (3.) and (3.13) into Eq. (3.6), use Eq. (3.5) and equate coefficients at (z) on both sides of the resulting equation, this gives C (λ) a n (λ) 4 1+( )+n +λ (n )( +n+λ) (0 ). (3.14) Insertion of Eqs. (3.4) and (3.14) into Eq. (3.) results in the expansion (1.6). 5

6 4 The Jacobi polynomials Finally, we shall determine coefficients in the expansion where P n α P n Γ(n+α+1) Γ(n+α+β +1) 0 0 a (β) n (α)p(α,β) (z), (4.1) Γ( +n+α+β +1)!(n )!Γ( +α+1) ( ) z 1 (4.) is the Jacobi polynomial. Observe that once these coefficients are determined, coefficients in the counterpart expansion of P n / β are also nown, because in virtue of the relationship one has P n ( ) n P n (β,α) ( z) (4.3) P n β ( ) +n a (α) 0 n (β)p(α,β) The coefficient at (z 1) n on the right-hand side of Eq. (4.) is n (α,β) hence, the coefficient a (β) nn(α) in the expansion (4.1) is (z). (4.4) Γ(n+α+β +1) n n!γ(n+α+β +1), (4.5) a (β) nn (α) ψ(n+α+β +1) ψ(n+α+β +1). (4.6) To find the remaining coefficients a (β) n (α) with 0, we differentiate the identity {(1 z ) d +[β α (α+β +)z] d } +n(n+α+β +1) P n 0 (4.7) with respect to α, obtaining {(1 z ) d +[β α (α+β +)z] d } (α,β) P +n(n+α+β +1) α Using the nown relationship (n+α+β)(z 1) dp(α,β) the right-hand side of Eq. (4.8) may be rewritten as (z +1) dp(α,β) (z +1) dp(α,β) np n. (4.8) np n (n+β) n+α+β The Christoffel Darboux formula for the Jacobi polynomials is!γ( +α+β +1) ( +α+β +1) Γ( +α+1)γ( +β +1) P(α,β) 0 n!γ(n+α+β +1) (n+α+β)γ(n+α)γ(n+β) n[β α+(n+α+β)z]p n (n+α)(n+β)p (α,β) (z), (4.9) np n (n+α)p (α,β) z 1 (z)p (α,β) (z ) (z). (4.10) P n P (α,β) (z ) P (α,β) (z)p(α,β) n (z ) z z. (4.11) 6

7 We set in this formula z 1 and then use This yields ( +α+β +1) 0 P (α,β) Γ( +α+1) (1)!Γ(α+1). (4.1) Γ( +α+β +1) Γ( +β +1) Γ(n+α+β +1) (n+α+β)γ(n+β) P (α,β) (z) np n (n+α)p (α,β) z 1 (z). (4.13) Combining this result with Eq. (4.10) gives us the following expansion of the expression standing on the right-hand side of Eq. (4.8): (z +1) dp(α,β) np n Γ(n+β +1) Γ(n+α+β +1) ( +α+β +1) 0 Γ( +α+β +1) P (α,β) (z). (4.14) Γ( +β +1) Substitution of Eq. (4.1) into the left-hand side and of Eq. (4.14) into the right-hand side of Eq. (4.8), followed by the use of Eq. (4.7), after some obvious movements leads to a (β) n (α) +α+β +1 Γ(n+β +1)Γ( +α+β +1) (n )( +n+α+β +1) Γ( +β +1)Γ(n+α+β +1) (0 ). (4.15) Plugging Eqs. (4.6) and (4.15) into Eq. (4.1) results in the relation (1.), while Eqs. (4.4), (4.6) and (4.15) imply Eq. (1.3). A Appendix In this Appendix we shall prove the following Lemma 1. The coefficient a nn (λ) in the expansion (1.1) is given by a nn (λ) ln n(λ), (A.1) where n (λ) is the coefficient at (z z 0 ) n (with z 0 C arbitrary) in the expansion of the polynomial P n (λ;z) in powers of z z 0. Proof. We have P n (λ;z) n (λ)(z z 0 ) n +Q (λ,z 0 ;z) (A.) (obviously, n (λ) is independent of z 0 ), where Q (λ,z 0 ;z) is a polynomial in z z 0 of degree. Hence, it follows that P n (λ;z) On the other side, from Eqs. (1.1) and (A.) we have P n (λ;z) n(λ) (z z 0) n + Q (λ,z 0 ;z). (A.3) a nn (λ)p n (λ;z)+ a n (λ)p (λ;z) 0 a nn (λ) n (λ)(z z 0 ) n +R (λ,z 0 ;z), (A.4) where R (λ,z 0 ;z) is a polynomial in z z 0 of degree. Equating coefficients at (z z 0 ) n on the right-hand sides of Eqs. (A.3) and (A.4) yields the expression (A.1). 7

8 References [1] G. Szegö, Orthogonal Polynomials, American Mathematical Society, New Yor, 1939 [] W. Magnus, F. Oberhettinger, R. P. Soni, Formulas and Theorems for the Special Functions of Mathematical Physics, 3rd ed., Springer, Berlin, 1966, chapter V [3] U. W. Hochstrasser, Orthogonal Polynomials, in: Handboo of Mathematical Functions, 10th printing, edited by M. Abramowitz and I. A. Stegun, Dover, New Yor, 197, chapter [4] J. Fröhlich, Parameter derivatives of the Jacobi polynomials and the Gaussian hypergeometric function, Integral Transforms Spec. Funct. (1994) [5] Yu. A. Brychov, Handboo of Special Functions. Derivatives, Integrals, Series and Other Formulas, Chapman & Hall/CRC, Boca Raton, FL, 008. The reader should notice that the factor at Cn λ (z) in Eq. (1...) may be simplified considerably [cf Eq. (1...1)] [6] W. Koepf, Identities for families of orthogonal polynomials and special functions, Integral Transforms Spec. Funct. 5 (1997) [7] A. Ronveaux, A. Zarzo, I. Area, E. Godoy, Classical orthogonal polynomials: dependence on parameters, J. Comput. Appl. Math. 11 (000) [8] R. Szmytowsi, On the derivative of the associated Legendre function of the first ind of integer degree with respect to its order (with applications to the construction of the associated Legendre function of the second ind of integer degree and order), J. Math. Chem. 46 (009) [9] R. Szmytowsi, On the derivative of the associated Legendre function of the first ind of integer order with respect to its degree (with applications to the construction of the associated Legendre function of the second ind of integer degree and order), preprint arxiv: [10] R. Szmytowsi, On parameter derivatives of the associated Legendre function of the first ind (with applications to the construction of the associated Legendre function of the second ind of integer degree and order), preprint arxiv: