Appendix/Mathematics/ToroidalConfusion

From jetwiki
Jump to navigation Jump to search


Confusion Regarding Whipple Formulae

May, 2018 (J.E.Tohline): I am trying to figure out what the correct relationship is between half-integer degree, associated Legendre functions of the first and second kinds. In order to illustrate my current confusion, here I will restrict my presentation to expressions that give Qn1/2m(coshη) in terms of Pm1/2n(cothη).

Published Expressions

From equation (34) of H. S. Cohl, J. E. Tohline, A. R. P. Rau, & H. M. Srivastiva (2000, Astronomische Nachrichten, 321, no. 5, 363 - 372) I find:

Expression #1

Qn1/2m(coshη)

=

(1)nπΓ(nm+12)[π2sinhη]1/2Pm1/2n(cothη).


From Howard Cohl's online overview of toroidal functions, I find:

Expression #2

Qm1/2n(coshα)

=

(1)nΓ(nm+12)[π2sinhα]1/2Pn1/2m(cothα),

Copying the Whipple's formula from §14.19 of DLMF,

Expression #3

Qn12m(coshξ)

=

Γ(mn+12)Γ(m+n+12)(π2sinhξ)1/2Pm12n(cothξ).

So far, this gives me three similar but not identical formulae for the same function mapping! As per equation (8) in (yet another source!) A. Gil, J. Segura, & N. M. Temme (2000, JCP, 161, 204 - 217), the relationship is:

Qn1/2m(λ)

=

(1)nπ3/22Γ(nm+1/2)(x21)1/4Pm1/2n(x),

Gil, Segura, & Temme (2000):  eq. (8)

where:    

λx/x21

This expression from Gil et al. (2000) means, for example, that by identifying x with cothη, we have λ=coshη, and,

Qn1/2m(coshη)

=

(1)nπ3/22Γ(nm+1/2)(coth2η1)1/4Pm1/2n(cothη)

 

=

(1)nπΓ(nm+1/2)(π2)1/2[cosh2ηsinh2η1]1/4Pm1/2n(cothη)

 

=

(1)nπΓ(nm+1/2)(π2)1/2[1sinhη]1/2Pm1/2n(cothη).

That is, we have,

Expression #4

Qn1/2m(coshη)

=

(1)nπΓ(nm+1/2)[π2sinhη]1/2Pm1/2n(cothη),

which matches the above expression #1 drawn from Cohl et al. (2000), but which appears not to match either of the other two "published" (online) formulae, expressions #2 or #3.

Specific Application

I stumbled into this dilemma when I tried to explicitly demonstrate how Q1/2(coshη) can be derived from P1/2(z) where, from §8.13 of M. Abramowitz & I. A. Stegun (1995), we find,

Q1/2(coshη)

=

2eη/2K(eη),

Abramowitz & Stegun (1995), eq. (8.13.4)

and,

P1/2(z)

=

2π[2z+1]1/2K(z1z+1).

Abramowitz & Stegun (1995), eq. (8.13.1)

When I used the Whipple formula as defined in §14.19 of DLMF (expression #3 reprinted above), the function mapping gave me the wrong result; I was off by a factor of Γ(12)=π. But, as demonstrated below, the Whipple formula provided by Cohl et al. (2000) and by Gil et al. (2000) — that is, expressions #1 and #4, above — does give the correct result.

Demonstration that Q12 can be derived from P12

Copying equation (34) from Cohl et al. (2000), we begin with,

Qn1/2m(coshη)

=

(1)nπΓ(nm+12)[π2sinhη]1/2Pm1/2n(cothη);

then setting m=n=0, we have,

Q12(coshη)

=

πΓ(12)[π2sinhη]1/2P12(cothη)

 

=

π2[1sinhη]1/2P12(cothη).

Step #1:   Associate … zcoshη. Then,

Q12(coshη)

=

π2[1z21]1/2P12(zz21).

Step #2:   Now making the association … Λz/z21, and drawing on eq. (8.13.1) from Abramowitz & Stegun (1995), we can write,

P12(Λ)

=

2π[2Λ+1]1/2K(Λ1Λ+1)

 

=

2π[2z21z+z21]1/2K(zz21z+z21).

Step #3:   Again, making the association … zcoshη, means,

P12(Λ)

=

2π[2sinhηcoshη+sinhη]1/2K(coshηsinhηcoshη+sinhη)

Q12(coshη)

=

π2[1sinhη]1/22π[2sinhηcoshη+sinhη]1/2K(coshηsinhηcoshη+sinhη)

 

=

2[1coshη+sinhη]1/2K(cosh2ηsinh2η[coshη+sinhη]2)

 

=

2[1eη]1/2K(1e2η)

 

=

2eη/2K(eη).

This, indeed, matches eq. (8.13.4) from Abramowitz & Stegun (1995).

Cohl's Response to My (May 2018) Email Query

Proper Interpretation of DLMF Expression

Most of the confusion expressed above stems from the DLMF's use of bold fonts, such as the function on the left-hand side of expression #3, above — that is, the Whipple formula from §14.19 of DLMF,

Qn12m(coshξ)

=

Γ(mn+12)Γ(m+n+12)(π2sinhξ)1/2Pm12n(cothξ).

What has been missing in my discussion is an appreciation of the following relationship between bold and plain-text function names,

Qνμ(x)=eμπiQνμ(x)Γ(ν+μ+1).

After making the substitutions, μm and ν(n12), the Whipple formula displayed above as expression #3 becomes,

emπiQn12m(coshξ)Γ(n+m+12)

=

Γ(mn+12)Γ(m+n+12)(π2sinhξ)1/2Pm12n(cothξ)

Qn12m(coshξ)

=

emπiΓ(mn+12)(π2sinhξ)1/2Pm12n(cothξ)

 

=

(1)mΓ(mn+12)(π2sinhξ)1/2Pm12n(cothξ),

which matches expression #2, above. But it does not appear to match expressions #1 or #4.

The standard "Euler reflection formula for gamma functions" is usually presented in the form,

Γ(z)Γ(1z)

=

πsin(πz)

|

for example, if
z(mn+12)

Γ(mn+12)Γ(nm+12)

=

π{sin[π2+π(mn)]}1

 

=

π(1)mn

DLMF §5.5(ii)

|
Valid for:

   z0,±1,±2,

|

If we make the association,

z(mn+12),

with m and n both being either zero or a positive integer, then, this Euler reflection formula becomes,

Γ(mn+12)Γ(nm+12)

=

π{sin[π(mn+12)]}1

 

=

π(1)m+n.

However, in our situation the so-called "Euler reflection formula for gamma functions" gives the relation,

π(1)m+nΓ(nm+12)

=

Γ(mn+12).

Hence, we may also write,

Qn12m(coshξ)

=

(1)m[π(1)m+nΓ(nm+12)](π2sinhξ)1/2Pm12n(cothξ)

 

=

(1)nπΓ(nm+12)(π2sinhξ)1/2Pm12n(cothξ),

which matches expressions #1 and #4. So everything appears to be in agreement! Hooray!

Derivation From Scratch

Whenever he deals with these types of relations, Cohl usually begins with,

Expression #5

Qνμ(coshη)

=

π2Γ(ν+μ+1)eiμπ[1sinh2η]1/4Pμ12ν12(cothη)

Making the pair of substitutions,

ν

=

n12,

     

n

0={0,1,2,},

μ

=

m,

     

m

0={0,1,2,},

we also have,

ν+μ+1

=

n12+m+1

=

n+m+12,

μ12

=

m12,

 

 

ν12

=

(n12)12

=

n,

eiμπ

=

eimπ

=

(1)m,

in which case,

Qn12m(coshη)

=

π2Γ(n+m+12)(1)m[1sinhη]Pm12n(cothη).


Now, since,

Pνμ(z)

=

Pν1μ(z),

if we make the substitution,

(ν+1)

(m+12)

       

ν

m12,

we also know that,

Pm12μ(z)

=

Pm12μ(z).

Hence, we can write,

Qn12m(coshη)

=

π2Γ(n+m+12)(1)m[1sinhη]Pm12n(cothη).


Finally, another relation states that, for n0,

Pm12n(z)

=

[Γ(mn+12)Γ(m+n+12)]Pm12n(z).

So, we obtain,

Qn12m(coshη)

=

(1)mπ2Γ(n+m+12)[1sinhη][Γ(mn+12)Γ(m+n+12)]Pm12n(cothη).

 

=

(1)mπ2Γ(mn+12)[1sinhη]Pm12n(cothη).

This matches expressions #2 and #3, above.

Index Values of Zero

Setting n=m=0 gives the following sought-for relationship:

Q120(coshη)

=

π2Γ(12)[1sinhη]P120(cothη).

 

=

π2[1sinhη]P120(cothη).

Joel's Additional Manipulations

From §14.19.6 of DLMF, we find the following summation expression:

Q12μ(coshξ)+2n=1Γ(μ+n+12)Γ(μ+12)Qn12μ(coshξ)cos(nϕ)=(12π)1/2(sinhξ)μ(coshξcosϕ)μ+(1/2)

Then, if we again employ the DLMF relationship between bold and plain-text function names, namely,

Qn12μ(x)=eμπiQn12μ(x)Γ(μ+n+12),

where we have made the substitution, ν(n12), the Sums expression becomes,

eμπiQ12μ(coshξ)Γ(μ+12)

=

(12π)1/2(sinhξ)μ(coshξcosϕ)μ+(1/2)2n=1Γ(μ+n+12)Γ(μ+12)[eμπiQn12μ(coshξ)Γ(μ+n+12)]cos(nϕ)

Q12μ(coshξ)

=

eμπiΓ(μ+12)[(12π)1/2(sinhξ)μ(coshξcosϕ)μ+(1/2)]2n=1Qn12μ(coshξ)cos(nϕ).

When dealing with Dyson-Wong tori, we will set μ=0, in which case the Sums expression becomes,

Q12(coshξ)

=

[π/2(coshξcosϕ)12]2n=1Qn12(coshξ)cos(nϕ).

But this can be rewritten in the form,

n=0ϵnQn12(coshξ)cos(nϕ)

=

[π/2(coshξcosϕ)12]

See Also

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |