Appendix/Mathematics/ToroidalConfusion

From JETohlineWiki
Jump to navigation Jump to search


Confusion Regarding Whipple Formulae[edit]

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[edit]

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[edit]

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[edit]

Proper Interpretation of DLMF Expression[edit]

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[edit]

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[edit]

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[edit]

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[edit]

Tiled Menu

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