12 Parabolic Cylinder FunctionsProperties

§12.10 Uniform Asymptotic Expansions for Large Parameter

Contents
  1. §12.10(i) Introduction
  2. §12.10(ii) Negative a, 2⁢−a<x<∞
  3. §12.10(iii) Negative a, −∞<x<−2⁢−a
  4. §12.10(iv) Negative a, −2⁢−a<x<2⁢−a
  5. §12.10(v) Positive a, −∞<x<∞
  6. §12.10(vi) Modifications of Expansions in Elementary Functions
  7. §12.10(vii) Negative a, −2⁢−a<x<∞. Expansions in Terms of Airy Functions
  8. §12.10(viii) Negative a, −∞<x<2⁢−a. Expansions in Terms of Airy Functions

§12.10(i) Introduction

In this section we give asymptotic expansions of PCFs for large values of the parameter a that are uniform with respect to the variable z, when both a and z (=x) are real. These expansions follow from Olver (1959), where detailed information is also given for complex variables.

With the transformations

12.10.1 a =±12⁢μ2,
x =μ⁢t⁢2,

(12.2.2) becomes

12.10.2 d2wdt2=μ4⁢(t2±1)⁢w.

With the upper sign in (12.10.2), expansions can be constructed for large μ in terms of elementary functions that are uniform for t∈(−∞,∞) (§2.8(ii)). With the lower sign there are turning points at t=±1, which need to be excluded from the regions of validity. These cases are treated in §§12.10(ii)–12.10(vi).

The turning points can be included if expansions in terms of Airy functions are used instead of elementary functions (§2.8(iii)). These cases are treated in §§12.10(vii)–12.10(viii).

Throughout this section the symbol δ again denotes an arbitrary small positive constant.

§12.10(ii) Negative a, 2⁢−a<x<∞

As a→−∞

12.10.3 U⁡(−12⁢μ2,μ⁢t⁢2)∼g⁡(μ)⁢e−μ2⁢ξ(t2−1)14⁢∑s=0∞𝒜s⁡(t)μ2⁢s,
12.10.4 U′⁡(−12⁢μ2,μ⁢t⁢2)∼−μ2⁢g⁡(μ)⁢(t2−1)14⁢e−μ2⁢ξ⁢∑s=0∞ℬs⁡(t)μ2⁢s,
12.10.5 V⁡(−12⁢μ2,μ⁢t⁢2)∼2⁢g⁡(μ)Γ⁡(12+12⁢μ2)⁢eμ2⁢ξ(t2−1)14⁢∑s=0∞(−1)s⁢𝒜s⁡(t)μ2⁢s,
12.10.6 V′⁡(−12⁢μ2,μ⁢t⁢2)∼2⁢μ⁢g⁡(μ)Γ⁡(12+12⁢μ2)⁢(t2−1)14⁢eμ2⁢ξ⁢∑s=0∞(−1)s⁢ℬs⁡(t)μ2⁢s,

uniformly for t∈[1+δ,∞), where

12.10.7 ξ=12⁢t⁢t2−1−12⁢ln⁡(t+t2−1).

The coefficients are given by

12.10.8 𝒜s⁡(t)=us⁡(t)(t2−1)32⁢s,ℬs⁡(t)=vs⁡(t)(t2−1)32⁢s,

where us⁡(t) and vs⁡(t) are polynomials in t of degree 3⁢s, (s odd), 3⁢s−2 (s even, s≥2). For s=0,1,2,

12.10.9 u0⁡(t) =1,
u1⁡(t) =t⁢(t2−6)24,
u2⁡(t) =−9⁢t4+249⁢t2+1451152,
12.10.10 v0⁡(t) =1,
v1⁡(t) =t⁢(t2+6)24,
v2⁡(t) =15⁢t4−327⁢t2−1431152.

Higher polynomials us⁡(t) can be calculated from the recurrence relation

12.10.11 (t2−1)⁢us′⁡(t)−3⁢s⁢t⁢us⁡(t)=rs−1⁡(t),

where

12.10.12 8⁢rs⁡(t)=(3⁢t2+2)⁢us⁡(t)−12⁢(s+1)⁢t⁢rs−1⁡(t)+4⁢(t2−1)⁢rs−1′⁡(t),

and the vs⁡(t) then follow from

12.10.13 vs⁡(t)=us⁡(t)+12⁢t⁢us−1⁡(t)−rs−2⁡(t).

Lastly, the function g⁡(μ) in (12.10.3) and (12.10.4) has the asymptotic expansion:

12.10.14 g⁡(μ)∼h⁡(μ)⁢(1+12⁢∑s=1∞γs(12⁢μ2)s),

where

12.10.15 h⁡(μ)=2−14⁢μ2−14⁢e−14⁢μ2⁢μ12⁢μ2−12,

and the coefficients γs are defined by

12.10.16 Γ⁡(12+z)∼2⁢π⁢e−z⁢zz⁢∑s=0∞γszs;

compare (5.11.8). For s≤4

12.10.17 γ0 =1,
γ1 =−124,
γ2 =11152,
γ3 =10034 14720,
γ4 =−4027398 13120.

§12.10(iii) Negative a, −∞<x<−2⁢−a

When μ→∞, asymptotic expansions for the functions U⁡(−12⁢μ2,−μ⁢t⁢2) and V⁡(−12⁢μ2,−μ⁢t⁢2) that are uniform for t∈[1+δ,∞) are obtainable by substitution into (12.2.15) and (12.2.16) by means of (12.10.3) and (12.10.5). Similarly for U′⁡(−12⁢μ2,−μ⁢t⁢2) and V′⁡(−12⁢μ2,−μ⁢t⁢2).

§12.10(iv) Negative a, −2⁢−a<x<2⁢−a

As a→−∞

12.10.18 U⁡(−12⁢μ2,μ⁢t⁢2) ∼2⁢g⁡(μ)(1−t2)14⁢(cos⁡κ⁢∑s=0∞(−1)s⁢𝒜~2⁢s⁡(t)μ4⁢s−sin⁡κ⁢∑s=0∞(−1)s⁢𝒜~2⁢s+1⁡(t)μ4⁢s+2),
12.10.19 U′⁡(−12⁢μ2,μ⁢t⁢2) ∼μ⁢2⁢g⁡(μ)⁢(1−t2)14⁢(sin⁡κ⁢∑s=0∞(−1)s⁢ℬ~2⁢s⁡(t)μ4⁢s+cos⁡κ⁢∑s=0∞(−1)s⁢ℬ~2⁢s+1⁡(t)μ4⁢s+2),
12.10.20 V⁡(−12⁢μ2,μ⁢t⁢2) ∼2⁢g⁡(μ)Γ⁡(12+12⁢μ2)⁢(1−t2)14⁢(cos⁡χ⁢∑s=0∞(−1)s⁢𝒜~2⁢s⁡(t)μ4⁢s−sin⁡χ⁢∑s=0∞(−1)s⁢𝒜~2⁢s+1⁡(t)μ4⁢s+2),
12.10.21 V′⁡(−12⁢μ2,μ⁢t⁢2) ∼μ⁢2⁢g⁡(μ)⁢(1−t2)14Γ⁡(12+12⁢μ2)⁢(sin⁡χ⁢∑s=0∞(−1)s⁢ℬ~2⁢s⁡(t)μ4⁢s+cos⁡χ⁢∑s=0∞(−1)s⁢ℬ~2⁢s+1⁡(t)μ4⁢s+2),

uniformly for t∈[−1+δ,1−δ]. The quantities κ and χ are defined by

12.10.22 κ =μ2⁢η−14⁢π,
χ =μ2⁢η+14⁢π,

where

12.10.23 η=12⁢arccos⁡t−12⁢t⁢1−t2,

and the coefficients 𝒜~s⁡(t) and ℬ~s⁡(t) are given by

12.10.24 𝒜~s⁡(t) =us⁡(t)(1−t2)32⁢s,
ℬ~s⁡(t) =vs⁡(t)(1−t2)32⁢s;

compare (12.10.8).

§12.10(v) Positive a, −∞<x<∞

As a→∞

12.10.25 U⁡(12⁢μ2,μ⁢t⁢2)∼g¯⁡(μ)⁢e−μ2⁢ξ¯(t2+1)14⁢∑s=0∞u¯s⁡(t)(t2+1)32⁢s⁢1μ2⁢s,

uniformly for t∈ℝ. Here bars do not denote complex conjugates; instead

12.10.26 ξ¯=12⁢t⁢t2+1+12⁢ln⁡(t+t2+1),
12.10.27 u¯s⁡(t)=is⁢us⁡(−i⁢t),

and the function g¯⁡(μ) has the asymptotic expansion

12.10.28 g¯⁡(μ)∼1μ⁢2⁢h⁡(μ)⁢(1+12⁢∑s=1∞(−1)s⁢γs(12⁢μ2)s),

where h⁡(μ) and γs are as in §12.10(ii).

With the same conditions

12.10.29 U′⁡(12⁢μ2,μ⁢t⁢2)∼−μ2⁢g¯⁡(μ)⁢(t2+1)14⁢e−μ2⁢ξ¯⁢∑s=0∞v¯s⁡(t)(t2+1)32⁢s⁢1μ2⁢s,

where

12.10.30 v¯s⁡(t)=is⁢vs⁡(−i⁢t).

§12.10(vi) Modifications of Expansions in Elementary Functions

In Temme (2000) modifications are given of Olver’s expansions. An example is the following modification of (12.10.3)

12.10.31 U⁡(−12⁢μ2,μ⁢t⁢2)∼h⁡(μ)⁢e−μ2⁢ξ(t2−1)14⁢∑s=0∞𝖠s⁡(τ)μ2⁢s,

where ξ and h⁡(μ) are as in (12.10.7) and (12.10.15) ,

12.10.32 τ=12⁢(tt2−1−1),

and the coefficients 𝖠s⁡(τ) are the product of τs and a polynomial in τ of degree 2⁢s. They satisfy the recursion

12.10.33 𝖠s+1⁡(τ)=−4⁢τ2⁢(τ+1)2⁢ddτ⁡𝖠s⁡(τ)−14⁢∫0τ(20⁢u2+20⁢u+3)⁢𝖠s⁡(u)⁢du,
s=0,1,2,…,

starting with 𝖠0⁡(τ)=1. Explicitly,

12.10.34 𝖠1⁡(τ) =−112⁢τ⁢(20⁢τ2+30⁢τ+9),
𝖠2⁡(τ) =1288⁢τ2⁢(6160⁢τ4+18480⁢τ3+19404⁢τ2+8028⁢τ+945).

The modified expansion (12.10.31) shares the property of (12.10.3) that it applies when μ→∞ uniformly with respect to t∈[1+δ,∞). In addition, it enjoys a double asymptotic property: it holds if either or both μ and t tend to infinity. Observe that if t→∞, then 𝖠s⁡(τ)=O⁡(t−2⁢s), whereas 𝒜s⁡(t)=O⁡(1) or O⁡(t−2) according as s is even or odd. The proof of the double asymptotic property then follows with the aid of error bounds; compare §10.41(iv).

For additional information see Temme (2000). See also Olver (1997b, pp. 206–208) and Jones (2006).

§12.10(vii) Negative a, −2⁢−a<x<∞. Expansions in Terms of Airy Functions

The following expansions hold for large positive real values of μ, uniformly for t∈[−1+δ,∞). (For complex values of μ and t see Olver (1959).)

12.10.35 U⁡(−12⁢μ2,μ⁢t⁢2) ∼2⁢π12⁢μ13⁢g⁡(μ)⁢ϕ⁡(ζ)⁢(Ai⁡(μ43⁢ζ)⁢∑s=0∞As⁡(ζ)μ4⁢s+Ai′⁡(μ43⁢ζ)μ83⁢∑s=0∞Bs⁡(ζ)μ4⁢s),
12.10.36 U′⁡(−12⁢μ2,μ⁢t⁢2) ∼(2⁢π)12⁢μ23⁢g⁡(μ)ϕ⁡(ζ)⁢(Ai⁡(μ43⁢ζ)μ43⁢∑s=0∞Cs⁡(ζ)μ4⁢s+Ai′⁡(μ43⁢ζ)⁢∑s=0∞Ds⁡(ζ)μ4⁢s),
12.10.37 V⁡(−12⁢μ2,μ⁢t⁢2) ∼2⁢π12⁢μ13⁢g⁡(μ)⁢ϕ⁡(ζ)Γ⁡(12+12⁢μ2)⁢(Bi⁡(μ43⁢ζ)⁢∑s=0∞As⁡(ζ)μ4⁢s+Bi′⁡(μ43⁢ζ)μ83⁢∑s=0∞Bs⁡(ζ)μ4⁢s),
12.10.38 V′⁡(−12⁢μ2,μ⁢t⁢2) ∼(2⁢π)12⁢μ23⁢g⁡(μ)ϕ⁡(ζ)⁢Γ⁡(12+12⁢μ2)⁢(Bi⁡(μ43⁢ζ)μ43⁢∑s=0∞Cs⁡(ζ)μ4⁢s+Bi′⁡(μ43⁢ζ)⁢∑s=0∞Ds⁡(ζ)μ4⁢s).

The variable ζ is defined by

12.10.39 23⁢ζ32 =ξ,1≤t,(ζ≥0);
23⁢(−ζ)32 =η,−1<t≤1,(ζ≤0),

where ξ,η are given by (12.10.7), (12.10.23), respectively, and

12.10.40 ϕ⁡(ζ)=(ζt2−1)14.

The function ζ=ζ⁡(t) is real for t>−1 and analytic at t=1. Inversely, with w=2−13⁢ζ,

12.10.41 t=1+w−110⁢w2+11350⁢w3−82363000⁢w4+1 50653242 55000⁢w5+⋯,
|ζ|<(34⁢π)23.

For g⁡(μ) see (12.10.14). The coefficients As⁡(ζ) and Bs⁡(ζ) are given by

12.10.42 As⁡(ζ) =ζ−3⁢s⁢∑m=02⁢sβm⁢(ϕ⁡(ζ))6⁢(2⁢s−m)⁢u2⁢s−m⁡(t),
ζ2⁢Bs⁡(ζ) =−ζ−3⁢s⁢∑m=02⁢s+1αm⁢(ϕ⁡(ζ))6⁢(2⁢s−m+1)⁢u2⁢s−m+1⁡(t),

where ϕ⁡(ζ) is as in (12.10.40), uk⁡(t) is as in §12.10(ii), α0=1, and

12.10.43 αm =(2⁢m+1)⁢(2⁢m+3)⁢⋯⁢(6⁢m−1)m!⁢(144)m,
βm =−6⁢m+16⁢m−1⁢αm.

The coefficients Cs⁡(ζ) and Ds⁡(ζ) in (12.10.36) and (12.10.38) are given by

12.10.44 Cs⁡(ζ) =χ⁡(ζ)⁢As⁡(ζ)+As′⁡(ζ)+ζ⁢Bs⁡(ζ),
Ds⁡(ζ) =As⁡(ζ)+χ⁡(ζ)⁢Bs−1⁡(ζ)+Bs−1′⁡(ζ),

where

12.10.45 χ⁡(ζ)=ϕ′⁡(ζ)ϕ⁡(ζ)=1−2⁢t⁢(ϕ⁡(ζ))64⁢ζ.

Explicitly,

12.10.46 ζ⁢Cs⁡(ζ) =−ζ−3⁢s⁢∑m=02⁢s+1βm⁢(ϕ⁡(ζ))6⁢(2⁢s−m+1)⁢v2⁢s−m+1⁡(t),
Ds⁡(ζ) =ζ−3⁢s⁢∑m=02⁢sαm⁢(ϕ⁡(ζ))6⁢(2⁢s−m)⁢v2⁢s−m⁡(t),

where vk⁡(t) is as in §12.10(ii).

Modified Expansions

The expansions (12.10.35)–(12.10.38) can be modified, again see Temme (2000), and the new expansions hold if either or both μ and t tend to infinity. This is provable by the methods used in §10.41(v).

§12.10(viii) Negative a, −∞<x<2⁢−a. Expansions in Terms of Airy Functions

When μ→∞, asymptotic expansions for U⁡(−12⁢μ2,−μ⁢t⁢2) and V⁡(−12⁢μ2,−μ⁢t⁢2) that are uniform for t∈[−1+δ,∞) are obtained by substitution into (12.2.15) and (12.2.16) by means of (12.10.35) and (12.10.37). Similarly for U′⁡(−12⁢μ2,−μ⁢t⁢2) and V′⁡(−12⁢μ2,−μ⁢t⁢2).