Difference between revisions of "User:Tohline/2DStructure/ToroidalGreenFunction"

From VistrailsWiki
Jump to navigation Jump to search
 
(44 intermediate revisions by the same user not shown)
Line 158: Line 158:
and <math>~\theta</math> has the same sign as <math>~(z-Z_0)</math>.
and <math>~\theta</math> has the same sign as <math>~(z-Z_0)</math>.


<span id="DiffVolumeElement">&nbsp;</span>
<table border="1" cellpadding="10" align="center" width="85%">
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 16 August 2017:  In equation (2.17) of his &sect;IIB &#8212; when Wong (1973) introduces the differential volume element &#8212; the variable used to represent the azimuthal coordinate angle switches from &psi; to &#x003A6;We will stick with the &psi; notation, here.]]According to p. 1301, eq. (10.3.75) of [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>] &#8212; or, for example, as found in [https://en.wikipedia.org/wiki/Toroidal_coordinates#Scale_factors Wikipedia] &#8212; in toroidal coordinates the differential volume element is,
<tr>
  <td align="center">[[File:Hicks1881TitlePage.png|500px|Title Page of Hicks (1881)]]</td>
</tr>
<tr><td align="left">
[http://rstl.royalsocietypublishing.org/content/172/609.full.pdf+html W. M. Hicks (1881)] presents one of the first &#8212; if not ''the'' first &#8212; discussions of toroidal coordinates and associated toroidal functionsEquation (4) on p. 614 of his article provides the following definition of the pair of meridional-plane coordinates <math>~(u,v)</math>, written in terms of the traditional cylindrical-coordinate pair <math>~(\rho,z)</math> and the specified ''anchor ring'' radius, <math>~a</math>:
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~d^3x</math>
<math>~u</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 170: Line 174:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~h_\eta h_\theta h_\psi d\eta d\theta d\psi</math>
<math>~
\frac{1}{2} \ln \biggl[ \frac{z^2 + (\rho+a)^2}{z^2 + (\rho - a)^2} \biggr] \, ,
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~v</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 176: Line 188:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{a^3 \sinh\eta}{(\cosh\eta - \cos\theta)^3} \biggr] d\eta~ d\theta~ d\psi \, .</math>
<math>~
\tan^{-1} \biggl[ \frac{2az}{\rho^2 + z^2 - a^2} \biggr] \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
 
His equation (5) presents the reverse mapping, namely,
 
==Green's Function Expression==
 
===As presented by Wong (1973)===
Referencing [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] states that, in toroidal coordinates, the Green's function is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
<math>~\rho</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 198: Line 206:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
\frac{a\sinh u}{\cosh u - \cos v} \, ,
\sum\limits^\infty_{m,n=0} (-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
</math>
</math>
   </td>
   </td>
Line 206: Line 213:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~z</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~\begin{cases}P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^') ~~~\eta^' > \eta \\P^m_{n-1 / 2}(\cosh\eta^') ~Q^m_{n-1 / 2}(\cosh\eta)~~~\eta^' < \eta
\frac{a\sin v}{\cosh u - \cos v} \, .
\end{cases}\, ,
</math>
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.53)<br />
[see also: [http://adsabs.harvard.edu/abs/1997JMP....38.3679B J. W. Bates (1997)], p. 3685, Eq. (31)]<br />
[see also: [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl, Tohline, Rau, &amp; Srivastava (2000)], &sect;6.2, Eq. (48)]
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
where, <math>~P^m_{n-1 / 2}, Q^m_{n-1 / 2}</math> are ''Associated Legendre Functions'' of the first and second kind with order <math>~n - \tfrac{1}{2}</math> and degree <math>~m</math> (toroidal harmonics), and <math>~\epsilon_m</math> is the Neumann factor, that is, <math>~\epsilon_0 = 1</math> and <math>~\epsilon_m = 2</math> for all <math>~m \ge 1</math>.  This Green's function expression can indeed be found as eq. (10.3.81) on p. 1304 of [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], but it should be noted that the MF53 expression differs from Wong's in two respects (see footnote 2 on p. 370 of [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl et al. (2000)] for a proposed explanation):  First, the factor, <math>~(-1)^m</math>, appears as <math>~(-i)^m</math> in MF53; and, second, in the term that is composed of a ratio of gamma functions, the denominator appears in MF53 as <math>~\Gamma(n - m + \tfrac{1}{2})</math>, whereas it should be <math>~\Gamma(n + m + \tfrac{1}{2})</math>, as presented here.


===Rearranging Terms and Incorporating Special-Function Relations===
</td></tr></table>


Let's focus on the situation when <math>~\eta^' > \eta</math>, and begin rearranging or substituting terms.
<span id="DiffVolumeElement">&nbsp;</span>
<div align="center">
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 16 August 2017:  In equation (2.17) of his &sect;IIB &#8212; when Wong (1973) introduces the differential volume element &#8212; the variable used to represent the azimuthal coordinate angle switches from &psi; to &#x003A6;.  We will stick with the &psi; notation, here.]]According to p. 1301, eq. (10.3.75) of [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>] &#8212; or, for example, as found in [https://en.wikipedia.org/wiki/Toroidal_coordinates#Scale_factors Wikipedia] &#8212; in toroidal coordinates the differential volume element is,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
<math>~d^3x</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~h_\eta h_\theta h_\psi d\eta d\theta d\psi</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 243: Line 246:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{a^3 \sinh\eta}{(\cosh\eta - \cos\theta)^3} \biggr] d\eta~ d\theta~ d\psi \, .</math>
\frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] 
</math>
   </td>
   </td>
</tr>
</tr>
</table>
==Green's Function Expression==
===As presented by Wong (1973)===
Referencing [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] states that, in toroidal coordinates, the Green's function is,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')]  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
~P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^')
\sum\limits^\infty_{m,n=0} (-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
</math>
</math>
   </td>
   </td>
Line 270: Line 278:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~\begin{cases}P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^') ~~~\eta^' > \eta \\P^m_{n-1 / 2}(\cosh\eta^') ~Q^m_{n-1 / 2}(\cosh\eta)~~~\eta^' < \eta
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] 
\end{cases}\, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.53)<br />
[see also: [http://adsabs.harvard.edu/abs/1997JMP....38.3679B J. W. Bates (1997)], p. 3685, Eq. (31)]<br />
[see also: [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl, Tohline, Rau, &amp; Srivastava (2000)], &sect;6.2, Eq. (48)]
  </td>
</tr>
</table>
</div>
where, <math>~P^m_{n-1 / 2}, Q^m_{n-1 / 2}</math> are ''Associated Legendre Functions'' of the first and second kind with degree <math>~n - \tfrac{1}{2}</math> and order <math>~m</math> (toroidal harmonics), and <math>~\epsilon_m</math> is the Neumann factor, that is, <math>~\epsilon_0 = 1</math> and <math>~\epsilon_m = 2</math> for all <math>~m \ge 1</math>.  This Green's function expression can indeed be found as eq. (10.3.81) on p. 1304 of [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], but it should be noted that the MF53 expression differs from Wong's in two respects (see footnote 2 on p. 370 of [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl et al. (2000)] for a proposed explanation):  First, the factor, <math>~(-1)^m</math>, appears as <math>~(-i)^m</math> in MF53; and, second, in the term that is composed of a ratio of gamma functions, the denominator appears in MF53 as <math>~\Gamma(n - m + \tfrac{1}{2})</math>, whereas it should be <math>~\Gamma(n + m + \tfrac{1}{2})</math>, as presented here.
<!--
For later reference, note that after drawing from, for example, [https://en.wikipedia.org/wiki/Gamma_function#General Wikipedia's account of the general properties of gamma functions],  the collection of factors immediately inside the double summation may be more explicitly written as,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~L_n^m \equiv (-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~
(-1)^m \epsilon_m \epsilon_n \biggl[ \frac{\sqrt{\pi}[2(n-m)]!}{4^{n-m}(n-m)!} \biggr]
\sqrt{ \frac{\pi}{2} }~\Gamma(n-m+\tfrac{1}{2}) \sqrt{\sinh\eta}~P^m_{n-1 / 2}(\cosh\eta)  
\biggl[ \frac{4^{n+m}(n+m)!}{\sqrt{\pi}[2(n+m)]!} \biggr]
\biggl\}\biggr\{
~\sqrt{ \frac{2}{\pi} }~\frac{\sqrt{\sinh\eta^'}}{\Gamma(n + m + \tfrac{1}{2})} Q^m_{n-1 / 2}(\cosh\eta^')
\biggr\}
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
The term contained within the first set of curly braces on the right-hand side of this expression can be replaced by the derived expression labeled <font color="green" size="+1">&#x2460;</font> in the [[#A.1|Appendix, below]], and simultaneously the term contained within the second set of curly braces can be replaced by the derived expression labeled <font color="green" size="+1">&#x2461;</font> in the [[#A.2|same Appendix]].  After making these substitutions, we have,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 312: Line 326:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
(-1)^m \epsilon_m \epsilon_n \biggl[ \frac{[2(n-m)]!}{(n-m)!} \biggr]
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] 
\biggl[ \frac{(n+m)!}{[2(n+m)]!} \biggr]2^{4m}
</math>
</math>
   </td>
   </td>
Line 323: Line 337:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~
(-1)^m \epsilon_m \epsilon_n \biggl[ \frac{\cdots [2(n-m)-4]\cdot [2(n-m)-3]\cdot[2(n-m)-2]\cdot[2(n-m)-1]\cdot[2(n-m)] }{\cdots [(n-m)-4]\cdot [(n-m)-3]\cdot [(n-m)-2]\cdot [(n-m)-1]\cdot (n-m) } \biggr]
(-1)^{-n}Q^n_{m-1 / 2}(\coth\eta)  
\biggl\}\biggr\{
~(-1)^{-m} P^{-n}_{m-1 / 2}(\coth\eta^')
\biggr\}
</math>
</math>
   </td>
   </td>
Line 341: Line 351:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\times
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\biggl[ \frac{\cdots [(n+m)-4]\cdot [(n+m)-3]\cdot [(n+m)-2]\cdot [(n+m)-1]\cdot (n+m) }{\cdots [2(n+m)-4]\cdot [2(n+m)-3]\cdot [2(n+m)-2]\cdot [2(n+m)-1]\cdot [2(n+m)] } \biggr]2^{4m}
\sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')]
</math>
</math>
   </td>
   </td>
Line 356: Line 365:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \sum\limits^\infty_{n=0} (-1)^{n} \epsilon_n \cos[n(\theta - \theta^')]  
(-1)^m \epsilon_m \epsilon_n \biggl[ \frac{\cdots 2[(n-m)-2]\cdot [2(n-m)-3]\cdot 2[(n-m)-1]\cdot [2(n-m)-1] \cdot 2(n-m) }{\cdots [(n-m)-4]\cdot [(n-m)-3]\cdot [(n-m)-2]\cdot [(n-m)-1]\cdot (n-m) } \biggr]
Q^n_{m-1 / 2}(\coth\eta) P^{-n}_{m-1 / 2}(\coth\eta^') \, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
where, in writing this last expression we have acknowledged that, since <math>~n</math> is either zero or a positive integer, <math>~(-1)^{-n} = (-1)^n</math>.  Next we draw upon the "Key Equation" relation,
{{ User:Tohline/Math/EQ_Toroidal01 }}
which, after making the substitutions, <math>~\nu \rightarrow (m - \tfrac{1}{2})</math> and <math>~\psi \rightarrow (\theta - \theta^')</math>, and incorporating the Neumann factor, <math>~\epsilon_n</math>, becomes,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
&nbsp;
Q_{m - \frac{1}{2} }\ [t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^') ]
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\times
\sum_{n=0}^\infty (-1)^n \epsilon_n Q^n_{m - \frac{1}{2} }(t) P^{-n}_{m - \frac{1}{2} }(t^') \cos[(n(\theta- \theta^')] \, .
\biggl[ \frac{\cdots [(n+m)-4]\cdot [(n+m)-3]\cdot [(n+m)-2]\cdot [(n+m)-1]\cdot (n+m) }{\cdots 2[(n+m)-2]\cdot [2(n+m)-3]\cdot 2[(n+m)-1]\cdot [2(n+m)-1]\cdot 2[(n+m)] } \biggr]2^{4m}
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
Finally, after making the associations, <math>~t \rightarrow \coth\eta</math> and <math>~t^' \rightarrow \coth\eta^'</math>, this last expression allows us to rewrite Wong's (1973) Green's function in a significantly more compact form, namely,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 403: Line 397:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
(-1)^m \epsilon_m \epsilon_n \biggl[ \frac{\cdots 2^{n-m}  [2(n-m)-3]\cdot  [2(n-m)-1] }{1} \biggr]
\sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')]  Q_{m - \frac{1}{2}}(\Chi) \, ,
\biggl[ \frac{1}{\cdots 2^{n+m} [2(n+m)-3]\cdot [2(n+m)-1]] } \biggr]2^{4m}
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
where the argument, <math>~\Chi</math>, of the toroidal function, <math>~Q_{m - \frac{1}{2}}</math>, is,
-->
===As Presented by Hicks (1881)===
Notably, more than 135 years ago, [http://rstl.royalsocietypublishing.org/content/172/609.full.pdf+html W. M. Hicks (1881)] already had constructed a Green's function expression for the reciprocal distance between two points in toroidal coordinates.  The following boxed-in table contains a snapshot image of  equation (35) from Hicks (1881), along with its associated text; notice that he refers to <math>~\phi</math> as the "potential."  Although notations are different &#8212; for example, <math>~C</math> is shorthand for <math>~\cosh u</math> and <math>~c</math> is shorthand for <math>~\cos v</math> &#8212; one can easily see factor-by-factor agreement when compared with the Green's function [[#As_presented_by_Wong_.281973.29|presented by Wong (1973)]].
<table border="1" cellpadding="10" align="center">
<tr>
  <td align="center">[[File:Hicks1881TitlePage.png|500px|Title Page of Hicks (1881)]]</td>
</tr>
<tr><td align="center">
[[File:Hicks1881GreenFunction.png|750px|Green's Function from Hicks (1881)]]
</td></tr></table>
 
===Rearranging Terms and Incorporating Special-Function Relations===


Let's focus on the situation when <math>~\eta^' > \eta</math>, and begin rearranging or substituting terms.
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Chi</math>
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^')  
\frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')
</math>
</math>
   </td>
   </td>
Line 432: Line 440:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\coth\eta \coth\eta^' - (\coth^2\eta-1)^{1 / 2} (\coth^2\eta'- 1)^{1 / 2} \cos(\theta- \theta^')  
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')]  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
~P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^')
</math>
</math>
   </td>
   </td>
Line 450: Line 459:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{\cosh\eta \cosh\eta^'}{\sinh\eta \sinh\eta^'} - \biggl[ \frac{1}{\sinh^2\eta} \biggr]^{1 / 2} \biggl[ \frac{1}{\sinh^2\eta'}\biggr]^{1 / 2} \cos(\theta- \theta^')  
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')
</math>
</math>
   </td>
   </td>
Line 460: Line 470:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{\cosh\eta \cosh\eta^' -  \cos(\theta- \theta^') }{\sinh\eta \sinh\eta^'} \, .
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~
\sqrt{ \frac{\pi}{2} }~\Gamma(n-m+\tfrac{1}{2}) \sqrt{\sinh\eta}~P^m_{n-1 / 2}(\cosh\eta)
\biggl\}\biggr\{
~\sqrt{ \frac{2}{\pi} }~\frac{\sqrt{\sinh\eta^'}}{\Gamma(n + m + \tfrac{1}{2})} Q^m_{n-1 / 2}(\cosh\eta^')
\biggr\}
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>


===As Presented in Cohl &amp; Tohline (1999)===
The term contained within the first set of curly braces on the right-hand side of this expression can be replaced by the derived expression labeled <font color="green" size="+1">&#x2460;</font> in the [[#A.1|Appendix, below]], and simultaneously the term contained within the second set of curly braces can be replaced by the derived expression labeled <font color="green" size="+1">&#x2461;</font> in the [[#A.2|same Appendix]].  After making these substitutions, we have,
<table border="0" cellpadding="5" align="center">


This last, compact Green's function expression &#8212; which we have derived, here, from [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] published Green's function by drawing strategically upon a variety of special function relations &#8212; precisely matches the "compact cylindrical Green's function expression" that has been derived independently by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)] via a less tortuous route, namely,
<table border="0" align="center" cellpadding="5">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~ \frac{1}{|\vec{x} - \vec{x}^{~'}|}</math>
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 482: Line 496:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~  
<math>~
\frac{1}{\pi \sqrt{\varpi \varpi^'}} \sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi)
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')]
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="2">&nbsp;</td>
   <td align="right">
   <td align="left" colspan="1">
&nbsp;
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (17)
  </td>
  <td align="center">
&nbsp;
  </td>
   <td align="left">
<math>~
\times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~
(-1)^{-n}Q^n_{m-1 / 2}(\coth\eta)
\biggl\}\biggr\{
~(-1)^{-m} P^{-n}_{m-1 / 2}(\coth\eta^')
\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
Line 501: Line 529:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~  
<math>~
\frac{1}{a\pi}
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\biggl[  \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}
\sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')]
\sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi) \, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Chi</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{(\varpi^')^2 + \varpi^2 + (z^' - z)^2}{2\varpi^' \varpi}
\times \sum\limits^\infty_{n=0} (-1)^{n} \epsilon_n \cos[n(\theta - \theta^')]
</math>
Q^n_{m-1 / 2}(\coth\eta) P^{-n}_{m-1 / 2}(\coth\eta^') \, ,
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
</table>
   <td align="center" colspan="2">
 
&nbsp; &nbsp;
where, in writing this last expression we have acknowledged that, since <math>~n</math> is either zero or a positive integer, <math>~(-1)^{-n} = (-1)^n</math>.  Next we draw upon the "Key Equation" relation,
 
{{ User:Tohline/Math/EQ_Toroidal01 }}
 
which, after making the substitutions, <math>~\nu \rightarrow (m - \tfrac{1}{2})</math> and <math>~\psi \rightarrow (\theta - \theta^')</math>, and incorporating the Neumann factor, <math>~\epsilon_n</math>, becomes,
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
<math>~
Q_{m - \frac{1}{2} }\ [t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^') ]
</math>
  </td>
  <td align="center">
<math>~=</math>
   </td>
   </td>
   <td align="left" colspan="3">
   <td align="left">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (16)
<math>~
\sum_{n=0}^\infty (-1)^n \epsilon_n Q^n_{m - \frac{1}{2} }(t) P^{-n}_{m - \frac{1}{2} }(t^') \cos[(n(\theta- \theta^')] \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>


Finally, after making the associations, <math>~t \rightarrow \coth\eta</math> and <math>~t^' \rightarrow \coth\eta^'</math>, this last expression allows us to rewrite Wong's (1973) Green's function in a significantly more compact form, namely,
<table border="0" cellpadding="5" align="center">


'''Note from J. E. Tohline (June, 2018)''':  &nbsp;This is the first time that I have been able to formally demonstrate to myself that these two separately derived Green's function expressions are identical. See, however, the earlier identification of ''new'' addition theorems in association with equations (49) and (50) of [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl et al. (2000)].
<tr>
 
  <td align="right">
==Gravitational Potential==
<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>
 
  </td>
Quite generally, then, the gravitational potential can be obtained at any coordinate location, <math>~(\eta,\theta,\psi)</math> &#8212; both inside and outside of a specified mass distribution &#8212; by carrying out three nested spatial integrals over the product of:&nbsp; <math>~\rho(\vec{x}^{~'})</math>, the [[#DiffVolumeElement|differential volume element]], and the Green's function as specified ''either'' by [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] or by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)]. 
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} }
\sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')] Q_{m - \frac{1}{2}}(\Chi) \, ,
</math>
  </td>
</tr>
</table>
where the argument, <math>~\Chi</math>, of the toroidal function, <math>~Q_{m - \frac{1}{2}}</math>, is,


In what follows we will make an effort to elucidate the pros and cons of adopting one Green's function expression over the other.  In each case we begin by writing the expression for the potential in such a way that variations in the azimuthal coordinate, <math>~\psi</math>, are described by ''Fourier components,'' <math>~\Phi_m^{(1)}(\eta,\theta)</math> and <math>~\Phi_m^{(2)}(\eta,\theta)</math>, of the potential, such that,
<table border="0" cellpadding="5" align="center">
<div align="center">
<math>~\Phi(\vec{x}) = \tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \Phi_m^{(1)}(\eta,\theta)  + \sum_{m=1}^\infty \sin (m\psi) \Phi_m^{(2)}(\eta,\theta) \, .</math>
</div>


Each Fourier component of the potential depends explicitly on the corresponding ''Fourier component'' of the density distribution, defined such that, 
<div align="center">
<math>~\rho(\vec{x}) = \tfrac{1}{2}\rho_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \rho_m^{(1)}(\eta,\theta)  + \sum_{m=1}^\infty \sin (m\psi) \rho_m^{(2)}(\eta,\theta) \, .</math>
</div>
<table border="1" align="center" cellpadding="8" width="70%">
<tr>
  <th align="center" bgcolor="yellow">
LaTeX mathematical expressions cut-and-pasted directly from
<br />
NIST's ''Digital Library of Mathematical Functions''
  </th>
</tr>
<tr>
<tr>
  <td align="right">
<math>~\Chi</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
   <td align="left">
   <td align="left">
As an additional primary point of reference, note that according to [https://dlmf.nist.gov/1.8 &sect;1.8(i) of NIST's ''Digital Library of Mathematical Functions''], a Fourier Series is defined as follows:
<math>~
<table border="0" cellpadding="5" align="center">
t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^')  
</math>
  </td>
</tr>


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~f(x)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 581: Line 620:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\tfrac{1}{2}a_{0}+\sum^{\infty}_{n=1}\biggl[ a_{n}\cos\bigl(nx\bigr)+b_{n}\sin\bigl(nx\bigr) \biggr],</math>
<math>~
\coth\eta \coth\eta^' - (\coth^2\eta-1)^{1 / 2} (\coth^2\eta'- 1)^{1 / 2} \cos(\theta- \theta^')  
</math>
   </td>
   </td>
</tr>
</tr>
Line 587: Line 628:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~a_n</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 593: Line 634:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\cos\bigl(nx\bigr)\mathrm{d}x,</math>
<math>~
\frac{\cosh\eta \cosh\eta^'}{\sinh\eta \sinh\eta^'} - \biggl[ \frac{1}{\sinh^2\eta} \biggr]^{1 / 2} \biggl[ \frac{1}{\sinh^2\eta'}\biggr]^{1 / 2} \cos(\theta- \theta^')  
</math>
   </td>
   </td>
</tr>
</tr>
Line 599: Line 642:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~b_n</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 605: Line 648:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\sin\bigl(nx\bigr)\mathrm{d}x.</math>
<math>~
\frac{\cosh\eta \cosh\eta^' -  \cos(\theta- \theta^') }{\sinh\eta \sinh\eta^'} \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>


  </td>
===As Presented in Cohl &amp; Tohline (1999)===
</tr>
</table>
 
Notice, therefore, that,
<table border="0" cellpadding="5" align="center">


This last, compact Green's function expression &#8212; which we have derived, here, from [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] published Green's function by drawing strategically upon a variety of special function relations &#8212; precisely matches the "compact cylindrical Green's function expression" that has been derived independently by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)] via a less tortuous route, namely,
<table border="0" align="center" cellpadding="5">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\rho_m^{(1)}(\eta,\theta)</math>
<math>~ \frac{1}{|\vec{x} - \vec{x}^{~'}|}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 625: Line 667:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\cos\bigl(m\psi\bigr)\mathrm{d}\psi,</math>
<math>~  
\frac{1}{\pi \sqrt{\varpi \varpi^'}} \sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi)
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="2">&nbsp;</td>
  <td align="left" colspan="1">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (17)
   </td>
   </td>
  <td align="center">&nbsp; &nbsp; and, &nbsp; &nbsp;</td>
</tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\rho_m^{(2)}(\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 635: Line 686:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\sin\bigl(m\psi\bigr)\mathrm{d}\psi \, .</math>
<math>~  
\frac{1}{a\pi}  
\biggl[  \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}
\sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi) \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
where,
<table border="0" cellpadding="5" align="center">


===The CT99 Expression for the Potential===
====In Three-Dimensional Generality====
Employing the Green's function expression derived by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], the gravitational potential for any three-dimensional matter distribution is,
<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~ \Phi(\eta,\theta,\psi)</math>
<math>~\Chi</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ -G
<math>~
\iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\}  \biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^'
\frac{(\varpi^')^2 + \varpi^2 + (z^' - z)^2}{2\varpi^' \varpi}
</math>
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 668: Line 713:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ -\frac{Ga^2}{\pi} \int d\eta^' \int d\theta^' \int d\psi^'
<math>~
\iiint \rho(\eta^',\theta^',\psi^') \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr]
\frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="center" colspan="2">
&nbsp;
&nbsp; &nbsp;
   </td>
   </td>
   <td align="center">
   <td align="left" colspan="3">
&nbsp;
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (16)
  </td>
  <td align="left">
<math>~ \times
\biggl[ \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}
\sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi)  
</math>
   </td>
   </td>
</tr>
</tr>
</table>
'''Note from J. E. Tohline (June, 2018)''':  &nbsp;This is the first time that I have been able to formally demonstrate to myself that these two separately derived Green's function expressions are identical.  See, however, the earlier identification of ''new'' addition theorems in association with equations (49) and (50) of [http://adsabs.harvard.edu/abs/2000AN....321..363C Cohl et al. (2000)].
==Gravitational Potential==


<tr>
Quite generally, then, the gravitational potential can be obtained at any coordinate location, <math>~(\eta,\theta,\psi)</math> &#8212; both inside and outside of a specified mass distribution &#8212; by carrying out three nested spatial integrals over the product of:&nbsp; <math>~\rho(\vec{x}^{~'})</math>, the [[#DiffVolumeElement|differential volume element]], and the Green's function as specified ''either'' by [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] or by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)].  
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \sum_{m=0}^{\infty} \epsilon_m
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
\int d\psi^' \rho(\eta^',\theta^',\psi^')   
</math>
  </td>
</tr>


<tr>
In what follows we will make an effort to elucidate the pros and cons of adopting one Green's function expression over the other.  In each case we begin by writing the expression for the potential in such a way that variations in the azimuthal coordinate, <math>~\psi</math>, are described by ''Fourier components,'' <math>~\Phi_m^{(1)}(\eta,\theta)</math> and <math>~\Phi_m^{(2)}(\eta,\theta)</math>, of the potential, such that,
<div align="center">
<math>~\Phi(\vec{x}) = \tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \Phi_m^{(1)}(\eta,\theta)  + \sum_{m=1}^\infty \sin (m\psi) \Phi_m^{(2)}(\eta,\theta) \, .</math>
</div>
 
Each Fourier component of the potential depends explicitly on the corresponding ''Fourier component'' of the density distribution, defined such that, 
<div align="center">
<math>~\rho(\vec{x}) = \tfrac{1}{2}\rho_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \rho_m^{(1)}(\eta,\theta)  + \sum_{m=1}^\infty \sin (m\psi) \rho_m^{(2)}(\eta,\theta) \, .</math>
</div>
 
<table border="1" align="center" cellpadding="8" width="70%">
<tr>
  <th align="center" bgcolor="yellow">
LaTeX mathematical expressions cut-and-pasted directly from
<br />
NIST's ''Digital Library of Mathematical Functions''
  </th>
</tr>
<tr>
  <td align="left">
As an additional primary point of reference, note that according to [https://dlmf.nist.gov/1.8 &sect;1.8(i) of NIST's ''Digital Library of Mathematical Functions''], a Fourier Series is defined as follows:
<table border="0" cellpadding="5" align="center">
 
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~f(x)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ \times
<math>~\tfrac{1}{2}a_{0}+\sum^{\infty}_{n=1}\biggl[ a_{n}\cos\bigl(nx\bigr)+b_{n}\sin\bigl(nx\bigr) \biggr],</math>
[\cos(m\psi)\cos(m\psi^') + \sin(m\psi)\sin(m\psi^') ]
   </td>
</math>
   </td>
</tr>
</tr>


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~a_n</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 726: Line 778:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\cos\bigl(nx\bigr)\mathrm{d}x,</math>
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi)
\int_{-\pi}^{\pi} d\psi^' \rho(\eta^',\theta^',\psi^') 
</math>
   </td>
   </td>
</tr>
</tr>
Line 735: Line 784:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~b_n</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi)
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\sin\bigl(nx\bigr)\mathrm{d}x.</math>
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
  </td>
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \cos(m\psi^')
</tr>
</math>
</table>
 
   </td>
   </td>
</tr>
</tr>
</table>
Notice, therefore, that,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\rho_m^{(1)}(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi)
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\cos\bigl(m\psi\bigr)\mathrm{d}\psi,</math>
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^')
\biggr\}
</math>
   </td>
   </td>
</tr>
  <td align="center">&nbsp; &nbsp; and, &nbsp; &nbsp;</td>
 
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\rho_m^{(2)}(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 772: Line 820:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{
<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\sin\bigl(m\psi\bigr)\mathrm{d}\psi \, .</math>
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi)  \rho_0^{(1)}(\eta^',\theta^')
</math>
   </td>
   </td>
</tr>
</tr>
</table>
===The CT99 Expression for the Potential===
====In Three-Dimensional Generality====
Employing the Green's function expression derived by [http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], the gravitational potential for any three-dimensional matter distribution is,


<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~ \Phi(\eta,\theta,\psi)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi)
<math>~ -G
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1)}(\eta^',\theta^')
\iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\\biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^'
</math>
</math>
   </td>
   </td>
Line 797: Line 850:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi)
<math>~ -\frac{Ga^2}{\pi} \int d\eta^' \int d\theta^' \int d\psi^'
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(2)}(\eta^',\theta^')
\iiint \rho(\eta^',\theta^',\psi^')  \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr]
\biggr\} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,


<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_m^{(1),(2)}(\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}  
<math>~ \times
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1),(2)}(\eta^',\theta^') \, ,
\biggl[ \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}
\sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi)  
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (20)
  </td>
</tr>
</table>
where, [[#As_Presented_in_Cohl_.26_Tohline_.281999.29|as above]],
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Chi</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \sum_{m=0}^{\infty} \epsilon_m
\frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, .
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
\int d\psi^' \rho(\eta^',\theta^',\psi^') 
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
====For Axisymmetric Systems====
For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>.  Hence, for all <math>~m > 0</math>, the ''Fourier components'' of the density, <math>~\rho_m^{(1),(2)}</math>, are zero.  The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>.  For axisymmetric systems, then, the gravitational potential is,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi(\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>
<math>~ \times
[\cos(m\psi)\cos(m\psi^') + \sin(m\psi)\sin(m\psi^') ]
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
Line 872: Line 911:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}  
<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- 1 / 2}(\Chi) \rho(\eta^',\theta^') \, .
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi)
\int_{-\pi}^{\pi} d\psi^' \rho(\eta^',\theta^',\psi^')
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="3">
   <td align="right">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eqs. (31) &amp; (32a)
&nbsp;
  </td>
</tr>
</table>
 
===Wong's Expression for the Potential===
 
====Fully Three-Dimensional Case====
 
Employing [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Green's function expression, the gravitational potential for any three-dimensional matter distribution is,
<table border="0" align="center">
<tr>
  <td align="right">
<math>~ \Phi(\vec{x})</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ -G
<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi)
\iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\}  \biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^'
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^')  \cos(m\psi^')
</math>
</math>
   </td>
   </td>
Line 909: Line 938:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - \frac{a^2G}{\pi} \int d\eta^' \int d\theta^' \int d\psi^'
<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi)
\biggl[ \frac{\rho(\eta^',\theta^',\psi^') ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr]  
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi)
\biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^')
\sum\limits^\infty_{n=0} \sum\limits^\infty_{m=0}(-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\biggr\}
</math>
</math>
   </td>
   </td>
Line 925: Line 954:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ - Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{
\times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)  
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi)  \rho_0^{(1)}(\eta^',\theta^')
</math>
</math>
   </td>
   </td>
Line 939: Line 968:
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2}
<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi)
\sum\limits^\infty_{n=0} \epsilon_n \sum\limits^\infty_{m=0}(-1)^m \epsilon_m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1)}(\eta^',\theta^')
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
</math>
</math>
   </td>
   </td>
Line 957: Line 985:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi)
\times~
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(2)}(\eta^',\theta^')
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\biggr\} \, .
\int d\psi^' \rho(\eta^',\theta^',\psi^') [\cos(m\psi)\cos(m\psi^') + \sin(m\psi) \sin(m\psi^')]
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,


<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Phi_m^{(1),(2)}(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 973: Line 1,004:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{
<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}  
\int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1),(2)}(\eta^',\theta^') \, ,
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^')
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eq. (20)
  </td>
</tr>
</table>
where, [[#As_Presented_in_Cohl_.26_Tohline_.281999.29|as above]],
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Chi</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
+~
\frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, .
\sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \cos(m\psi^')
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
====For Axisymmetric Systems====
For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>.  Hence, for all <math>~m > 0</math>, the ''Fourier components'' of the density, <math>~\rho_m^{(1),(2)}</math>, are zero.  The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>.  For axisymmetric systems, then, the gravitational potential is,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Phi(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>
+~
   </td>
\sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int_{-\pi}^\pi  d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^')
\biggr\}
</math>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="right">
Line 1,026: Line 1,057:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - a^2G (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{
<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}  
\int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- 1 / 2}(\Chi) \rho(\eta^',\theta^') \, .
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\rho_0^{(1)}(\eta^',\theta^')  
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1999ApJ...527...86C Cohl &amp; Tohline (1999)], p. 88, Eqs. (31) &amp; (32a)
  </td>
</tr>
</table>
===Wong's Expression for the Potential===


====Fully Three-Dimensional Case====
Employing [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Green's function expression, the gravitational potential for any three-dimensional matter distribution is,
<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~ \Phi(\vec{x})</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ -G
+~
\iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\} \biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^'
\sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\rho_m^{(1)}(\eta^',\theta^')
</math>
</math>
   </td>
   </td>
Line 1,057: Line 1,094:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ - \frac{a^2G}{\pi} \int d\eta^' \int d\theta^' \int d\psi^'
+~
\biggl[ \frac{\rho(\eta^',\theta^',\psi^') ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr]  
\sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 }
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\sum\limits^\infty_{n=0} \sum\limits^\infty_{m=0}(-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\rho_m^{(2)} (\eta^',\theta^')
\biggr\} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,


<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_m^{(1),(2)} (\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
- 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)  
\sum\limits^\infty_{n=0} \epsilon_n (-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
</math>
</math>
   </td>
   </td>
Line 1,095: Line 1,124:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\times
<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2}  
\int d\eta^' ~\sinh\eta^' ~ P^m_{n-1 / 2}(\cosh\eta_<) ~ Q^m_{n-1 / 2}(\cosh\eta_>)
\sum\limits^\infty_{n=0} \epsilon_n \sum\limits^\infty_{m=0}(-1)^m \epsilon_m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\theta^' \biggl\{ \frac{ \cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\rho_m^{(1),(2)}(\eta^',\theta^') \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
====Axisymmetric Systems====
For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>.  Hence, for all <math>~m > 0</math>, the ''Fourier components'' of the density, <math>~\rho_m^{(1),(2)}</math>, are zero.  The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>.  For axisymmetric systems, then, the gravitational potential is,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi(\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>
<math>~
\times~
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int d\psi^' \rho(\eta^',\theta^',\psi^') [\cos(m\psi)\cos(m\psi^') + \sin(m\psi) \sin(m\psi^')]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,132: Line 1,158:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{
- 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\sum\limits^\infty_{n=0} \epsilon_n  
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int d\eta^' ~\sinh\eta^'~P^0_{n-1 / 2}(\cosh\eta_<) ~Q^0_{n-1 / 2}(\cosh\eta_>)
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^')  
\int d\theta^' \biggl\{ \frac{ ~\cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\}
\rho(\eta^',\theta^') \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.55)
  </td>
</tr>
</table>
==Uniform-Density Torus==
From an [[User:Tohline/Apps/DysonWongTori#Introducing_Toroidal_Coordinates|accompanying review]], we see that for a torus of major radius <math>~R</math> and minor radius <math>~d</math>, the most accommodating (cylindrical-coordinate-based) radial  coordinate location of the toroidal coordinate system's ''anchor ring'', <math>~a</math>, is defined such that,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~a^2</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~R^2 - d^2 \, .</math>
<math>~
+~
\sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \cos(m\psi^')
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="3">
  <td align="right">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eq. (2.8)
&nbsp;
  </td>
   <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~
\sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\int_{-\pi}^\pi  d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^')
\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</div>
Then the corresponding (toroidal-coordinate-based) radial coordinate location <math>~\eta_0</math> of the surface of the torus is,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\eta_0</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,184: Line 1,211:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\cosh^{-1}\biggl(\frac{R}{d}\biggr) \, .</math>
<math>~ - a^2G (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{
\int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\rho_0^{(1)}(\eta^',\theta^')
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eq. (2.9)
  </td>
</tr>
</table>
</div>
Alternatively, given <math>~\eta_0</math> and the value of the parameter <math>~a</math>, we have,
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~R</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~a \coth\eta_0 \, ,</math>
<math>~
+~
\sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
\rho_m^{(1)}(\eta^',\theta^')
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,212: Line 1,239:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~d</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{a}{\sinh\eta_0} \, .</math>
<math>~
  </td>
+~
</tr>
\sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
<tr>
\int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>)
  <td align="center" colspan="3">
\int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')]
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eqs. (2.10) &amp; (2.11)
\rho_m^{(2)} (\eta^',\theta^')
\biggr\} \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>


If such an axisymmetric torus has a uniform density, <math>~\rho_0</math>, then according to the Cohl &amp; Tohline (1999) expression, the gravitational potential is,
We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,
 
<table border="0" cellpadding="5" align="center">


<table border="0" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{CT}(\eta,\theta)</math>
<math>~\Phi_m^{(1),(2)} (\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,241: Line 1,268:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ - 2G\rho_0 a^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2}  
<math>~
\int_{\eta_0}^\infty d\eta^' \int_{-\pi}^\pi d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- 1 / 2}(\Chi) \, .
- 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\sum\limits^\infty_{n=0} \epsilon_n (-1)^m  ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})}
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
In contrast, according to Wong's (1973) expression, the gravitational potential is,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
&nbsp;
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\times
- 2G \rho_0 a^2 (\cosh\eta - \cos\theta)^{1 / 2}
\int d\eta^' ~\sinh\eta^' ~ P^m_{n-1 / 2}(\cosh\eta_<) ~ Q^m_{n-1 / 2}(\cosh\eta_>)
\sum\limits^\infty_{n=0} \epsilon_n
\int d\theta^' \biggl\{ \frac{ \cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\}  
\int_{\eta_0}^\infty d\eta^' ~\sinh\eta^'~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\rho_m^{(1),(2)}(\eta^',\theta^') \, .
\int_{-\pi}^\pi \frac{ \cos[n(\theta - \theta^')] d\theta^' }{(\cosh\eta^' - \cos\theta^')^{5/2}}  
\, .
</math>
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.55)
   </td>
   </td>
</tr>
</tr>
</table>
</table>


===Which Is Easier to Evaluate?===
====Axisymmetric Systems====
Upon evaluation, the expressions for <math>~\Phi_\mathrm{CT}</math> and for <math>~\Phi_\mathrm{W}</math> should give identical values for the gravitational potential at all coordinate locations <math>~(\eta,\theta)</math> across the entire meridional planeBut which of the two expressions is easier to evaluate?
For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>.  Hence, for all <math>~m > 0</math>, the ''Fourier components'' of the density, <math>~\rho_m^{(1),(2)}</math>, are zero.  The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>.  For axisymmetric systems, then, the gravitational potential is,


====CT====
If the desire is to perform the evaluation numerically, then the <math>~\Phi_\mathrm{CT}</math> expression is almost certainly easier to contend with.  It ''only'' requires a double integration; and as has been detailed in an [[User:Tohline/2DStructure/ToroidalCoordinates#Expression_for_the_Axisymmetric_Potential|accompanying discussion]], the only relatively unfamiliar special function that appears in the integrand, <math>~Q_{-\frac{1}{2}}(\Chi)</math>, can be re-expressed in terms of (the more familiar) complete elliptic integral of the first kind, namely,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{-\frac{1}{2}}(\Chi)</math>
<math>~\Phi(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,291: Line 1,305:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>
\mu K(\mu) \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\mu</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\equiv</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{2}{\Chi + 1} \biggr]^{1 / 2}</math>
<math>~
- 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\sum\limits^\infty_{n=0} \epsilon_n
\int d\eta^' ~\sinh\eta^'~P^0_{n-1 / 2}(\cosh\eta_<) ~Q^0_{n-1 / 2}(\cosh\eta_>)
\int d\theta^' \biggl\{ \frac{ ~\cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\}
\rho(\eta^',\theta^') \, .
</math>
   </td>
   </td>
  <td align="center">
</tr>
<math>~=</math>
<tr>
  </td>
   <td align="center" colspan="3">
   <td align="left">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.55)
<math>~\biggl[ \frac{2 \sinh\eta \cdot \sinh\eta^' }{ \cosh(\eta^' + \eta) - \cos(\theta^' - \theta) } \biggr]^{1 / 2} \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>


The task is daunting, however, if the desire is to evaluate the expression for <math>~\Phi_\mathrm{CT}</math> ''analytically.'' To our knowledge, no one has yet been able to start from the expression as presented here for <math>~\Phi_\mathrm{CT}</math> and complete the integral over the angular coordinate, <math>~\theta^'</math>, analytically, let alone analytically evaluate the second, ''radial'' integral.  This is principally because the integrand of the first (angular) integral contains a nontrivial function multiplied by a special function whose argument is, itself, a nontrivial function of both <math>~\theta^'</math> and <math>~\eta^'</math>.
==Uniform-Density Torus==
 
In an [[User:Tohline/Apps/Wong1973Potential#Wong.27s_.281973.29_Analytic_Potential|accompanying discussion]], we build upon the above technical foundation and detail how [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] was able to complete both integrals to derive an analytic expression for the potential (inside as well as outside) of axisymmetric, uniform-density tori having an arbitrarily specified ratio of the major to minor (cross-sectional) radii, <math>~R/d</math>.  This is an outstanding accomplishment that has received little attention in the astrophysics literature and, therefore, has heretofore been underappreciated.
 
=Appendix A: Selected Toroidal Function Relationships=


====W====
Here, we draw from the set of toroidal function relationships that have been identified as "Key Equations" in our [[User:Tohline/Appendix/Equation_templates#Torus|accompanying ''Equations'' appendix]].  Two related appendices &#8212; numbered '''A.3''' and '''A.4''' &#8212; may be found in an [[User:Tohline/Apps/Wong1973Potential#Appendix:_Selected_Toroidal_Function_Relationships|accompanying discussion]].


At first glance, the expression for <math>~\Phi_\mathrm{W}</math> is much more complex than the expression for <math>~\Phi_\mathrm{CT}</math>. It not only requires a double integration, but also an infinite-series summation over the integer index, <math>~n</math>; and the radial integrand contains a product of two (relatively unfamiliar) associated Legendre functions of various half-integer orders (toroidal functions).  Developing a numerical algorithm to evaluate <math>~\Phi_\mathrm{W}</math> to a certain accuracy &#8212; presumably set by the number of terms included in the summation over <math>~n</math> &#8212; would certainly be a more challenging task than developing a numerical algorithm to evaluate <math>~\Phi_\mathrm{CT}</math>.
==A.1==
 
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Beginning with the identified "Key Equation",
{{ User:Tohline/Math/EQ_Toroidal02 }}
 
we'll identify <math>~x</math> with <math>~\cosh\eta</math> &#8212; in which case we have <math>~\lambda = \coth\eta</math> &#8212; and switch the index notations, <math>~n \leftrightarrow m</math>. This gives,


=====Analytic Integration Over the Angular Coordinate=====
However, focusing only on the integration over the angular coordinate, we see that the integrand in the expression for <math>~\Phi_\mathrm{W}</math> is significantly less imposing than the one that appears in the expression for <math>~\Phi_\mathrm{CT}</math>.  [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)] was able to evaluate this definite integral in closed form, analytically.  While Wong does not record the detailed steps that he used to evaluate this definite integral, he does indicate that he received guidance from Volume I of [https://authors.library.caltech.edu/43491/1/Volume%201.pdf A. Erd&eacute;lyi's (1953)] ''Higher Transcendental Functions''.  We therefore presume that he adopted the line of reasoning that we have [[#A.3|detailed in the Appendix, below, in deriving the expression labeled]] <font color="green" size="+1">&#x2462;</font>.  Wong recognized, what we have explicitly demonstrated, namely,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\int_{-\pi}^\pi d\theta^' (\cosh\eta^' - \cos\theta^')^{- 5 / 2} \cos[n(\theta - \theta^')]</math>
<math>~Q_{m-1 / 2}^n (\coth\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~=</math>
   </td>
   </td>
   <td align="right">
   <td align="left">
<math>~2\cos(n\theta)
<math>~(-1)^m
\int_0^\pi \frac{ \cos(n\theta^')~d\theta^' }{ (\cosh\eta^' - \cos \theta^')^{\frac{5}{2}} }
\frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(m-n+\frac{1}{2})} (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, .
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{8\sqrt{2}}{3} \biggl[ \frac{\cos(n\theta)}{\sinh^2\eta^'} \biggr] Q^2_{n- \frac{1}{2}} (\cosh\eta^') \,.
</math>
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="5">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 293, Eq. (2.56)
   </td>
   </td>
</tr>
</tr>
</table>
</table>


<table border="1" cellpadding="10" align="center" width="90%"><tr><td align="left">
Drawing upon the ''Euler reflection formula for gamma functions'', namely,
<font color="red">'''CAUTION:'''</font>&nbsp; It is important to appreciate that, in this expression as well as in the expressions to follow, the term, <math>~Q^2_{n-\frac{1}{2}}(z)</math>, is ''not'' the square of the zero-order toroidal function, <math>~Q^0_{n - \frac{1}{2}}(z)</math>, but is instead the toroidal function of order two.  In an [[User:Tohline/Appendix/Mathematics/ToroidalSynopsis01#Evaluating_Q2.CE.BD|accompanying discussion]] we present an analytic expression for <math>~Q^2_{-\frac{1}{2}}(z)</math> &#8212; and a separate analytic expression for <math>~Q^1_{-\frac{1}{2}}(z)</math> &#8212; in terms of complete elliptic integrals, as well as a recurrence relation that can be used to generate analytic expressions for all other order-two (and all other order-one) toroidal functions that have higher half-integer degrees, <math>~n-\tfrac{1}{2}</math> for <math>~n \ge 1</math>.
{{ User:Tohline/Math/EQ_Gamma01 }}
</td></tr></table>
 
 
Hence, Wong was able to simplify the expression for <math>~\Phi_\mathrm{W}</math> to one that &#8212; albeit, in addition to an infinite summation over the index, <math>~n</math> &#8212; only requires integration over the radial coordinate, <math>~\eta^'</math>.  Specifically, he obtained,


where it is understood that <math>~m</math> and <math>~n</math> are each either zero or a positive integer, this toroidal-function relation becomes,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~Q_{m-1 / 2}^n (\coth\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,375: Line 1,381:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(-1)^m~
- 2G \rho_0 a^2 (\cosh\eta - \cos\theta)^{1 / 2}
\frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{ \Gamma(n - m +\frac{1}{2}) }{ \pi (-1)^{m-n} } \biggr] (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta)  
\sum\limits^\infty_{n=0} \epsilon_n
\int_{\eta_0}^\infty d\eta^' ~\sinh\eta^'~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)\biggl\{
\frac{8\sqrt{2}}{3} \biggl[ \frac{\cos(n\theta)}{\sinh^2\eta^'} \biggr] Q^2_{n- \frac{1}{2}} (\cosh\eta^')
\biggr\}
</math>
</math>
  </td>
  <td align="right" width="15%">
&nbsp;
   </td>
   </td>
</tr>
</tr>
Line 1,394: Line 1,399:
   <td align="left">
   <td align="left">
<math>~
<math>~
- \biggl( \frac{16\sqrt{2}}{3}  \biggr) G \rho_0 a^2 (\cosh\eta - \cos\theta)^{1 / 2}  
(-1)^n \sqrt{ \frac{\pi}{2} } ~\Gamma(n - m + \tfrac{1}{2} )(\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, .
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)
\int_{\eta_0}^\infty d\eta^' \Biggl[ \frac{Q^2_{n- \frac{1}{2}} (\cosh\eta^')}{\sinh\eta^'} \Biggr] ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>)
\, .
</math>
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2460;</font>
   </td>
   </td>
</tr>
</tr>
</table>


<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], p. 294, Eq. (2.57)
   </td>
   </td>
</tr>
</tr>
</table>
</table>


=====Analytic Integration Over the Radial Coordinate=====


Now, in considering how to handle integration over the radial coordinate, <math>~\eta^'</math>, let's examine, first, the case where <math>~\eta^' \ge \eta_0 > \eta</math>, that is, the potential is being evaluated at a location that is entirely outside of the torusIn this case,
==A.2==
 
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Again, beginning with the identified "Key Equation",
{{ User:Tohline/Math/EQ_Toroidal02 }}
 
this time, without switching index notations, we'll identify  <math>~x</math> with <math>~\coth\eta</math> &#8212; in which case we have <math>~\lambda = \cosh\eta</math>.  This gives,


<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 1,417: Line 1,427:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,423: Line 1,433:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(-1)^n
- \biggl( \frac{16\sqrt{2}}{3}  \biggr) G \rho_0 a^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(n-m+\frac{1}{2})} \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta) \, .
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta) ~P_{n-1 / 2}(\cosh\eta)  
\int_{\eta_0}^\infty d\eta^' \Biggl[ \frac{Q^2_{n- \frac{1}{2}} (\cosh\eta^')}{\sinh\eta^'} \Biggr] ~Q_{n-1 / 2}(\cosh\eta^')
\, .
</math>
</math>
   </td>
   </td>
Line 1,433: Line 1,440:
</table>
</table>


So we are interested in carrying out the "radial" integral,
Drawing upon the same ''Euler reflection formula for gamma functions'', as quoted above, this toroidal function relation can be rewritten as,


<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 1,439: Line 1,446:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
\int_{\eta_0}^\infty d\eta^' \sinh\eta^' \Biggl[ \frac{Q^2_{n- \frac{1}{2}} (\cosh\eta^')}{\sinh^2\eta^'} \Biggr] ~Q_{n-1 / 2}(\cosh\eta^')
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,447: Line 1,452:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(-1)^n
\int_{\cosh(\eta_0)}^{\cosh(\infty)} dt \Biggl[ \frac{Q^2_{n- \frac{1}{2}} (t)~Q_{n-\frac{1}{2}}(t)}{t^2 - 1} \Biggr]  \, ,
\frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\pi (-1)^{m-n}} \biggr] \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta)
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
where we have made the association, <math>~\cosh\eta^' \rightarrow t</math>, in which case, <math>~dt = \sinh\eta^' d\eta^'</math>.  Following the line of reasoning that we have [[#A.4|detailed in the Appendix, below, in deriving the expression labeled]] <font color="green" size="+1">&#x2463;</font>, this integral can be evaluated in closed form to give,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
&nbsp;
\int_a^b\biggl[  \frac{Q_{n - \frac{1}{2}}^2(t) ~Q_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,467: Line 1,466:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(-1)^{-m}~
\frac{1}{4} \biggl[
\sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr]  P_{m-1 / 2}^n(\coth\eta) \, .
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ Q_{n - \frac{1}{2}}(t)
- (n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)  
\biggr]_a^b \, .
</math>
</math>
  </td>
</tr>
<tr>
  <td align="center" colspan="3">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eq. (2.58)
   </td>
   </td>
</tr>
</tr>
</table>
</table>
Hence, we have,
Finally, calling upon the "Key Equation" relation,
 
{{ User:Tohline/Math/EQ_Toroidal00 }}


making the index notation substitution, <math>~\nu \rightarrow (m-\tfrac{1}{2})</math>, and associating <math>~z</math> with <math>~ \coth\eta</math> gives,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~P^n_{m-\frac{1}{2}}(\coth\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,493: Line 1,487:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\Gamma(m-n+\frac{1}{2})} \biggr]P^{-n}_{m-\frac{1}{2}}(\coth\eta) \, .</math>
- \biggl( \frac{4\sqrt{2}}{3}  \biggr) G \rho_0 a^2 (\cosh\eta - \cos\theta)^{1 / 2}  
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
</math>
   </td>
   </td>
</tr>
</tr>
</table>
As a result, we can write,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(-1)^{-m}~
\times P_{n-1 / 2}(\cosh\eta)
\sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr] P_{m-1 / 2}^{-n}(\coth\eta)
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ Q_{n - \frac{1}{2}}(t)
-  (n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)
\biggr]_{t = \cosh \eta_0}^{t = \cosh\infty} \, .
</math>
</math>
  </td>
  <td align="right" width="15%">
&nbsp;
   </td>
   </td>
</tr>
</tr>
</table>
Or, given that, in this uniform-density configuration, the [[User:Tohline/Apps/DysonWongTori#DensityFormula|density is exactly the mass divided by the torus volume]], that is,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\rho_0 = \frac{M}{V}</math>
<math>~\Rightarrow ~~~
P_{m-1 / 2}^{-n}(\coth\eta)
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,529: Line 1,522:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{M \sinh^3\eta_0}{2\pi^2 a^3 \cosh{\eta_0}} \, ,</math>
<math>~(-1)^m~\sqrt{\frac{2}{\pi}} \biggl[ \frac{\sqrt{\sinh\eta}}{\Gamma(m+n+\frac{1}{2})} \biggr] 
Q_{n-1 / 2}^m (\cosh\eta) \, .
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2461;</font>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
we have,
  </td>
</tr>
</table>
 
 
==A.3==


<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Here, our objective is to evaluate the definite integral,
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~\int_{-\pi}^{\pi} \frac{\cos[n(\theta - \theta^')] d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} </math>
   </td>
   </td>
   <td align="center">
   <td align="right">
<math>~=</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{2\sqrt{2}~ GM}{3\pi^2 a} \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}  
\cos(n\theta^') \int_{-\pi}^{\pi} \frac{ \cos(n\theta) ~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}}  
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
+ \sin(n\theta^') \cancelto{0}{ \int_{-\pi}^{\pi} \frac{ \sin(n\theta) d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} }
</math>
</math>
   </td>
   </td>
Line 1,556: Line 1,565:
&nbsp;
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="right">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times P_{n-1 / 2}(\cosh\eta)  
2 \cos(n\theta^') \int_{0}^{\pi} \frac{ \cos(n\theta)~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} \, .
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ Q_{n - \frac{1}{2}}(t)
(n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)
\biggr]_{t = \cosh \eta_0}^{t = \cosh\infty} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
Notice that, since the limits of the integration are <math>~-\pi</math> to <math>~+\pi</math>: &nbsp; The second integral on the right-hand-side goes to zero because the numerator of its integrand &#8212; ''i.e.,'' <math>~\sin(n\theta)</math> &#8212; is an [https://en.wikipedia.org/wiki/Even_and_odd_functions#Odd_functions odd function]; and, with regard to the first integral on the right-hand-side, the lower integration limit can be set to zero and the result doubled because the numerator of its integrand &#8212; ''i.e.,'' <math>~\cos(n\theta)</math> &#8212; is an [https://en.wikipedia.org/wiki/Even_and_odd_functions#Even_functions even function].
Drawing from Volume I of [https://authors.library.caltech.edu/43491/1/Volume%201.pdf A. Erd&eacute;lyi's (1953)] ''Higher Transcendental Functions,'' we find the following "Key Equation":


In a [[User:Tohline/Appendix/Mathematics/ToroidalFunctions#Asymptotic_Behavior|separate chapter]], we have examined the asymptotic behavior of the toroidal functions, <math>~Q_{n-\tfrac{1}{2}}(t)</math>, which is,
{{ User:Tohline/Math/EQ_Toroidal03 }}


Next we adopt the associations, <math>~z \rightarrow \cosh\eta</math>, <math>~t \rightarrow \theta</math>, <math>~\mu \rightarrow 2</math>, and, <math>~\nu \rightarrow n - \tfrac{1}{2}</math>, where <math>~n</math> is zero or a positive integer.  In this case we have,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\lim_{\chi\rightarrow \infty} Q_{n-\frac{1}{2}}(t)</math>
<math>~Q_{n - \frac{1}{2}}^2 (\cosh\eta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~\propto</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{1}{t^{n+\frac{1}{2}} } \, .
(2\pi)^{-\frac{1}{2}} (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~\biggl\{
\int_0^\pi (\cosh\eta - \cos \theta)^{-\frac{5}{2}} \cos(n\theta) ~d\theta
- \cancelto{0}{\cos[(n-\tfrac{1}{2})\pi] }~~\int_0^\infty (\cosh\eta + \cosh \theta)^{- \frac{5}{2}} e^{-n\theta} ~d\theta
\biggr\} \, ,
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
This means that, at the upper integration limit, these toroidal functions go to zero.  As a result, we can write,
where the prefactor of the second term &#8212; that is, <math>~\cos[(n-\tfrac{1}{2})\pi] </math> &#8212; goes to zero for all allowable values of the integer, <math>~n</math>.  Hence, we conclude that,
 
<table border="1" align="center" cellpadding="8">
<tr><td align="center">Exterior Solution: &nbsp;<math>~\eta_0 \ge \eta</math></td></tr>
<tr><td align="left">
 
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~2\cos(n\theta^')
\int_0^\pi \frac{ \cos(n\theta)~d\theta }{ (\cosh\eta - \cos \theta)^{\frac{5}{2}} }
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left" colspan="2">
<math>~
<math>~\frac{ 2(2\pi)^{\frac{1}{2}} Q_{n - \frac{1}{2}}^2 (\cosh\eta) \cos(n\theta^')}{ (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~ }
\frac{2\sqrt{2}~ GM}{3\pi^2 a}  \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
</math>
</math>
   </td>
   </td>
Line 1,618: Line 1,625:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{ 2^3 \sqrt{2} }{ 3 } \biggr] \frac{  Q_{n - \frac{1}{2}}^2(\cosh\eta) \cos(n\theta^')}{ \sinh^2\eta  }
\, .
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2462;</font>
  </td>
</tr>
</table>
where we have set,
<div align="center">
<math>~
<math>~
\times P_{n-1 / 2}(\cosh\eta)  
\Gamma(\tfrac{5}{2}) = \Gamma(\tfrac{1}{2} + 2) = \frac{ \sqrt{\pi} \cdot 4! }{4^2 \cdot 2!}
\biggl[
= \frac{\sqrt{\pi} \cdot 2^3\cdot 3}{ 2^5 } = \frac{3 \sqrt{\pi}}{2^2} \, .
(n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(\cosh \eta_0) Q_{n - \frac{1}{2}}^2(\cosh \eta_0)
- (n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(\cosh \eta_0) ~ Q_{n - \frac{1}{2}}(\cosh \eta_0)
\biggr] \, .
</math>
</math>
</div>
   </td>
   </td>
</tr>
</tr>
</table>
</table>


</td></tr>
==A.4==
</table>


This expression matches the potential for the ''exterior region'' that is obtained by combining [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Eqs. (2.59), (2.61) and (2.63).


Now, let's examine the potential evaluated at a radial location, <math>~\eta</math>, that is positioned inside the surface of the torus.  That is, <math>~\eta_0 < \eta \le \infty</math>.  We need to add expressions that have two different sets of integration limits as follows.  First, over the subregion, <math>~\eta \le \eta^' \le \infty</math>, we use the same expression for <math>~\Phi_\mathrm{W}</math> but employ the new limits, as indicated; and second, over the subregion <math>~\eta_0 \le \eta^' \le \eta</math>, we swap the roles of <math>~P_{n-\frac{1}{2}}</math> and <math>~Q_{n-\frac{1}{2}}</math>.
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Beginning with the identified "Key Equation",


{{ User:Tohline/Math/EQ_Toroidal05 }}
we will adopt the associations: &nbsp; <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math>, <math>~\nu \rightarrow (n - \tfrac{1}{2})</math>, <math>~\rho \rightarrow 0</math>, and <math>~\sigma \rightarrow ( n - \tfrac{1}{2})</math>.  As a result, [https://authors.library.caltech.edu/43491/1/Volume%201.pdf Erd&eacute;lyi's (1953)] expression becomes,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~
\int_a^b\biggl[ -4(1 - t^2)^{-1} \biggr] Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) ~dt
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,651: Line 1,673:
   <td align="left">
   <td align="left">
<math>~
<math>~
\Phi_\mathrm{W}|_\mathrm{subregion1} + \Phi_\mathrm{W}|_\mathrm{subregion2}
\biggl[
(n - \tfrac{1}{2} ) Q_{n - \frac{1}{2}}^2(t) ~ X_{n - \frac{3}{2}}(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t)
\biggr]_a^b \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
Drawing upon the ''recurrence'' "Key Equation,"
{{ User:Tohline/Math/EQ_Toroidal04 }}
which means, after making the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 0</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, that,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~(n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,664: Line 1,696:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~2n t X_{n-\frac{1}{2}}(t) - (n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t)</math>
- \frac{2\sqrt{2}~ GM}{3\pi^2 a} \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}  
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,673: Line 1,702:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow ~~~(n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t) </math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>
\times P_{n-1 / 2}(\cosh\eta)
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ Q_{n - \frac{1}{2}}(t)
- (n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)
\biggr]_{t = \cosh \eta}^{t = \cosh\infty}
</math>
   </td>
   </td>
</tr>
</tr>
</table>
the integral can be rewritten as,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
\int_a^b\biggl[  \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{2\sqrt{2}~ GM}{3\pi^2 a}  \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}  
\frac{1}{4} \biggl\{
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
\biggl[ 2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) \biggr] Q_{n - \frac{1}{2}}^2(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t)
\biggr\}_a^b
</math>
</math>
   </td>
   </td>
Line 1,709: Line 1,739:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times Q_{n-1 / 2}(\cosh\eta)  
\frac{1}{4} \biggl\{
\biggl[
\biggl[ 2n t  Q_{n - \frac{1}{2}}^2(t)  
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ P_{n - \frac{1}{2}}(t)
- (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) \biggr] X_{n - \frac{1}{2}}(t)
-  (n+\tfrac{1}{2})P_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)  
-  (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)  
\biggr]_{t = \cosh \eta_0}^{t = \cosh\eta}
\biggr\}_a^b \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
Returning to the same ''recurrence'' "Key Equation," but this time adopting the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, we can write,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,731: Line 1,764:
   <td align="left">
   <td align="left">
<math>~
<math>~
- \frac{2\sqrt{2}~ GM}{3\pi^2 a} \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}
2n t Q^2_{n - \frac{1}{2}}(t) - (n + \tfrac{3}{2}) Q^2_{n - \frac{3}{2}} (t) \, ,
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
in which case the integral becomes,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
\int_a^b\biggl[  \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times \biggl\{ P_{n-1 / 2}(\cosh\eta)
\frac{1}{4} \biggl\{
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ X_{n - \frac{1}{2}}(t)
(n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(\cosh\eta) Q_{n - \frac{1}{2}}^2(\cosh\eta)  
- (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)  
- (n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(\cosh\eta) ~ Q_{n - \frac{1}{2}}(\cosh\eta)
\biggr\}_a^b \, .
\biggr]
</math>
</math>
   </td>
   </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2463;</font>
  </td>
</tr>
</table>
</td>
</tr>
</tr>
</table>
=Appendix B: Elliptic Integrals=
==Summary Table==
If the default argument of a [[User:Tohline/Appendix/Equation_templates#Complete_Elliptic_Integrals|complete elliptic integral function]] is the parameter, <math>~k</math>, then we should also appreciate that the following additional arguments may be useful:


<table border="1" cellpadding="8" align="center">
<tr>
  <td align="center" rowspan="2" width="10%">parameter</td>
  <td align="center" colspan="4">In terms of &hellip;</td>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="center" width="20%"><math>~k</math><br />
&nbsp;
<font size="-1">([http://dlmf.nist.gov/22.7.i Descending Landen Transformation])</font>
   </td>
   </td>
   <td align="center">
   <td align="center" width="20%"><math>~k_1</math><br />
&nbsp;
<font size="-1">([http://dlmf.nist.gov/22.7.ii Ascending Landen Transformation])</font>
  </td>
  <td align="left">
<math>~
+ Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(\cosh\eta) ~ P_{n - \frac{1}{2}}(\cosh\eta)
-  (n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh\eta) Q_{n - \frac{1}{2}}^2(\cosh\eta)
\biggr]
</math>
   </td>
   </td>
  <td align="center" width="20%><math>~\frac{R}{R_1} = e^{-\eta}</math></td>
  <td align="center"><math>~\eta \equiv \ln\biggl(\frac{R_1}{R}\biggr)</math></td>
</tr>
<tr>
  <td align="center"><math>~k</math></td>
  <td align="center"><math>~k</math></td>
  <td align="center"><math>~\frac{2\sqrt{k_1}}{1+k_1}</math></td>
  <td align="center"><math>~\biggl[1 - \biggl( \frac{R}{R_1}\biggr)^2 \biggr]^{1 / 2}</math></td>
  <td align="center"><math>~\biggl[ \frac{2}{\coth\eta + 1} \biggr]^{1 / 2}= \biggl[ 1 - e^{-2\eta} \biggr]^{1 / 2}</math></td>
</tr>
</tr>
<tr>
<tr>
   <td align="right">
   <td align="center"><math>~k^{\prime}</math></td>
&nbsp;
   <td align="center"><math>~\sqrt{1 - k^2}</math></td>
   </td>
   <td align="center"><math>~\frac{1-k_1}{1+k_1}</math></td>
   <td align="center">
   <td align="center"><math>~\frac{R}{R_1}</math></td>
&nbsp;
   <td align="center"><math>~\frac{1}{\cosh\eta + \sinh\eta} = e^{-\eta}</math></td>
   </td>
</tr>
   <td align="left">
<tr>
<math>~
  <td align="center"><math>~k_1 = \frac{1-k^{\prime}}{1+k^{\prime}}</math></td>
- Q_{n-1 / 2}(\cosh\eta)
  <td align="center"><math>~\frac{1-\sqrt{1 - k^2}}{1+\sqrt{1 - k^2}}</math></td>
\biggl[
  <td align="center"><math>~k_1</math></td>
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(\cosh{\eta_0}) ~ P_{n - \frac{1}{2}}(\cosh{\eta_0})
  <td align="center"><math>~\frac{R_1 - R}{R_1 +R}</math></td>
-  (n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh{\eta_0}) Q_{n - \frac{1}{2}}^2(\cosh{\eta_0})
  <td align="center"><math>~\tanh\biggl(\frac{\eta}{2}\biggr) = \frac{1 - e^{-\eta}}{1 + e^{-\eta}}</math></td>
\biggr]
\biggr\}
</math>
  </td>
</tr>
</tr>
</table>
<span id="GandR65">As published, for example,</span> in &sect;8.126 of [http://www.mathtable.com/gr/ Gradshteyn &amp; Ryzhik (1965)], the following mappings are valid:
<table border="1" cellpadding="8" width="85%" align="center">
<tr>
<th>
Explicitly Published in &sect;8.126 of [http://www.mathtable.com/gr/ G &amp; R (1965)]
</th>
<th>
Additional Implications
</th>
</tr>
<tr>
  <td>
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
K\biggl( \frac{1-k^\prime}{1+k^\prime} \biggr)
=
K(k_1)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,800: Line 1,869:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{2} (1+k^{\prime}) K(k)</math>
\frac{2\sqrt{2}~ GM}{3\pi^2 a} \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)  
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,809: Line 1,875:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
E\biggl( \frac{1-k^\prime}{1+k^\prime} \biggr)
=
E(k_1)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{1+k^\prime} \biggl[E(k) + k^\prime K(k) \biggr]</math>
\times \biggl\{ P_{n-1 / 2}(\cosh\eta)
\biggl[
(n+\tfrac{1}{2})Q_{n+\frac{1}{2}}(\cosh\eta) Q_{n - \frac{1}{2}}^2(\cosh\eta)
\biggr]
- Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh\eta) Q_{n - \frac{1}{2}}^2(\cosh\eta)  
\biggr]
</math>
   </td>
   </td>
</tr>
</tr>
 
</table>
  </td>
  <td>
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
\frac{4K(k_1)}{R_1 + R}
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{2K(k)}{R_1}</math>
- Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(\cosh{\eta_0}) ~ P_{n - \frac{1}{2}}(\cosh{\eta_0})
-  (n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh{\eta_0}) Q_{n - \frac{1}{2}}^2(\cosh{\eta_0})
\biggr]
\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,849: Line 1,908:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~E(k)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,855: Line 1,914:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(1+k^\prime)E(k_1) - k^\prime K(k)</math>
-  \frac{2\sqrt{2}~ GM}{3\pi^2 a}  \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}
  </td>
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)
</tr>
</math>
 
<tr>
  <td align="center" colspan="3">[http://dlmf.nist.gov/19.8.E12 DLMS &sect;19.8.12]</td>
</tr>
</table>
   </td>
   </td>
</tr>
</tr>
<tr>
  <td>
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
K\biggl( \frac{2\sqrt{k}}{1+k}  \biggr)
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(1+k)K(k)</math>
\times \biggl\{ 
Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh{\eta_0}) Q_{n - \frac{1}{2}}^2(\cosh{\eta_0})
- (n - \tfrac{3}{2}) ~ P_{n - \frac{1}{2}}(\cosh{\eta_0})~Q^2_{n + \frac{1}{2}}(\cosh{\eta_0})
\biggr]
</math>
   </td>
   </td>
</tr>
</tr>
Line 1,883: Line 1,945:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~
E\biggl( \frac{2\sqrt{k}}{1+k} \biggr)
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{1+k} \biggl[2E(k) - (k^\prime)^2 K(k) \biggr]</math>
- (n+\tfrac{1}{2})Q_{n - \frac{1}{2}}^2(\cosh\eta)  \biggl[
Q_{n-1 / 2}(\cosh\eta) P_{n+\frac{1}{2}}(\cosh\eta)
- Q_{n+\frac{1}{2}}(\cosh\eta) P_{n-1 / 2}(\cosh\eta)
\biggr]
\biggr\}
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
 
  </td>
Drawing on the identified "Key Equation" expression,
  <td>
 
{{ User:Tohline/Math/EQ_Toroidal06 }}
 
and adopting the associations, <math>~z = \xi = \cosh\eta</math> and <math>~\ell \rightarrow (n-\tfrac{1}{2})</math>, we recognize that,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~(n+\tfrac{1}{2})[P_{n+\frac{1}{2}}(\cosh\eta) Q_{n-\frac{1}{2}} (\cosh\eta) - P_{n-\frac{1}{2}}(\cosh\eta)Q_{n+\frac{1}{2}}(\cosh\eta)] </math>
<math>~
K( k )
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,915: Line 1,971:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~(1+k_1)K(k_1)</math>
1 \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>


Hence, we can write,
<tr>
 
  <td align="center" colspan="3">[http://dlmf.nist.gov/19.8.E12 DLMS &sect;19.8.12]</td>
<table border="1" align="center" cellpadding="8">
</tr>
<tr><td align="center">Interior Solution: &nbsp;<math>~\eta \ge \eta_0 </math></td></tr>
<tr><td align="left">
 
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)</math>
<math>~
E\biggl( \frac{2\sqrt{k_1}}{1+k_1} \biggr) = E(k)
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,938: Line 1,989:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{1}{1+k_1} \biggl[2E(k_1) - (1-k_1^2) K(k_1) \biggr]</math>
\frac{2\sqrt{2}~ GM}{3\pi^2 a} \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] (\cosh\eta - \cos\theta)^{1 / 2}
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta)
</math>
   </td>
   </td>
</tr>
</tr>
 
</table>
  </td>
</tr>
</table>
 
==One Set of Details==
Following advice that we have received from Howard Cohl (private communication), we can demonstrate this most efficiently by employing the [https://dlmf.nist.gov/19.8#ii ''Descending Landen Transformation'' for the complete elliptic integral of the first kind], that is,
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~K(k)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(1 + k_1)K(k_1) \, ,
</math>
  </td>
<td align="center">&nbsp; &nbsp; &nbsp; where, &nbsp; &nbsp; &nbsp;</td>
  <td align="right">
<math>~k_1</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~  
\times \biggl\{ 
\frac{1-\sqrt{1-k^2}}{1+\sqrt{1-k^2}} \, .
Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh{\eta_0}) Q_{n - \frac{1}{2}}^2(\cosh{\eta_0})
- (n - \tfrac{3}{2}) ~ P_{n - \frac{1}{2}}(\cosh{\eta_0})~Q^2_{n + \frac{1}{2}}(\cosh{\eta_0})
\biggr] - Q_{n - \frac{1}{2}}^2(\cosh\eta) \biggr\} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
 
Given our just-stated definition of the parameter, <math>~k</math>, we have,
</td></tr>
</table>
 
 
This expression matches the potential for the ''interior region'' that is obtained by combining [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Eqs. (2.59), (2.60) and (2.62).  Finally, drawing on yet another identified "Key Equation," namely,
 
{{ User:Tohline/Math/EQ_Toroidal08 }}
 
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 2 July 2018:  Note that this matches [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Eq. (2.64), except that the numerator on the right-hand side of his equation contains  the hyperbolic sine function raised to the first power whereas in our expression it is squared.  We are confident that it should be squared, in part, because the term ''is'' squared when it reappears in Wong's Eq. (2.65).]]and adopting the associations, <math>~\mu \rightarrow 2</math>, <math>~\nu \rightarrow \theta</math>, and <math>~z \rightarrow \cosh\eta</math>, we see that,
 
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~
<math>~\sqrt{1-k^2}</math>
\sum_{n=0}^\infty \epsilon_n  \cos(n\theta) ~Q^2_{n-\frac{1}{2}}(\cosh\eta)
</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 1,988: Line 2,039:
   <td align="left">
   <td align="left">
<math>~
<math>~
(-1)^2~\biggl(\frac{\pi}{2} \biggr)^{1 / 2} \Gamma(2 + \tfrac{1}{2})\biggl[ \frac{ \sinh^2\eta }{(\cosh\eta - \cos\theta)^{\frac{5}{2}}} \biggr]
\biggl[ 1 - \frac{2}{\coth\eta + 1} \biggr]^{1 / 2} = \biggl[ \frac{\coth\eta - 1 }{\coth\eta + 1} \biggr]^{1 / 2}
</math>
</math>
   </td>
   </td>
Line 1,995: Line 2,046:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow ~~~ 1 \pm \sqrt{1-k^2}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,002: Line 2,053:
   <td align="left">
   <td align="left">
<math>~
<math>~
\biggl(\frac{\pi}{2} \biggr)^{1 / 2} \biggl[ \frac{4! \sqrt{\pi}}{ 4^2 \cdot 2! } \biggr] \biggl[ \frac{ \sinh^2\eta }{(\cosh\eta - \cos\theta)^{\frac{5}{2}} } \biggr]
1 \pm  \biggl[ \frac{\coth\eta - 1 }{\coth\eta + 1} \biggr]^{1 / 2} = \frac{ [\coth\eta + 1]^{1 / 2} \pm  [\coth\eta - 1]^{1 / 2}{ [\coth\eta + 1]^{1 / 2} }
</math>
</math>
   </td>
   </td>
Line 2,009: Line 2,060:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~\Rightarrow ~~~k_1 \equiv \frac{1 - \sqrt{1-k^2}}{1 + \sqrt{1-k^2}}</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,016: Line 2,067:
   <td align="left">
   <td align="left">
<math>~
<math>~
\biggl(\frac{1}{2} \biggr)^{1 / 2} \biggl[ \frac{ 3\pi}{ 2^2} \biggr] \biggl[ \frac{ \sinh^2\eta }{(\cosh\eta - \cos\theta)^{\frac{5}{2}} } \biggr] \, .
\frac{ [\coth\eta + 1]^{1 / 2} -  [\coth\eta - 1]^{1 / 2} }{ [\coth\eta + 1]^{1 / 2} +  [\coth\eta - 1]^{1 / 2} }
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="3">
  <td align="right">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eq. (2.64)<br />except, see accompanying (pink, scroll-over balloon) COMMENT
<math>~\Rightarrow ~~~k_1^2 </math>
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ [\coth\eta + 1] -2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} +  [\coth\eta - 1]  }{ [\coth\eta + 1] + 2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} +  [\coth\eta - 1] }
</math>
   </td>
   </td>
</tr>
</tr>
</table>
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 2 July 2018:  Note that this matches [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong's (1973)] Eq. (2.65), except that in his expression the first term inside the curly braces on the right-hand side contains an extra factor of &pi;.  We are confident that our derived expression is correct because it is consistent with the immediately preceding expression &#8212; labeled Eq. (2.64) in [http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)].]]Hence, the ''interior solution'' can be rewritten as,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~\Phi_\mathrm{W}(\eta,\theta)\bigr|_\mathrm{interior}</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,039: Line 2,094:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~ \frac{ GM}{2\pi a}  \biggl[ \frac{\sinh^3\eta_0}{ \cosh{\eta_0}} \biggr] \biggl\{ \biggl[ \frac{ \sinh\eta }{(\cosh\eta - \cos\theta) } \biggr]^2
<math>~
~ -  ~ \frac{2^2\sqrt{2}~ }{3\pi}  (\cosh\eta - \cos\theta)^{1 / 2}  
\frac{ 2\coth\eta -2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} }{ 2\coth\eta + 2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2}  }
\sum\limits^\infty_{n=0} \epsilon_n \cos(n\theta) \, .
</math>
</math>
   </td>
   </td>
Line 2,051: Line 2,105:
   </td>
   </td>
   <td align="center">
   <td align="center">
&nbsp;
<math>~=</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~
\times
\frac{ \coth\eta - [\coth^2\eta - 1]^{1 / 2} }{ \coth\eta + [\coth^2\eta - 1]^{1 / 2} }
Q_{n-1 / 2}(\cosh\eta)
\biggl[
(n+\tfrac{1}{2})P_{n+\frac{1}{2}}(\cosh{\eta_0}) Q_{n - \frac{1}{2}}^2(\cosh{\eta_0})
- (n - \tfrac{3}{2}) ~ P_{n - \frac{1}{2}}(\cosh{\eta_0})~Q^2_{n + \frac{1}{2}}(\cosh{\eta_0})
\biggr] \biggr\} \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center" colspan="3">
  <td align="right">
[http://adsabs.harvard.edu/abs/1973AnPhy..77..279W Wong (1973)], Eq. (2.65)<br />except, see accompanying (pink, scroll-over balloon) COMMENT
&nbsp;
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{ \coth\eta/[\coth^2\eta - 1]^{1 / 2}  - 1 }{ \coth\eta/[\coth^2\eta - 1]^{1 / 2}  +  1 }
</math>
   </td>
   </td>
</tr>
</tr>
</table>


=Appendix: Selected Toroidal Function Relationships=
Here, we draw from the set of toroidal function relationships that have been identified as "Key Equations" in our [[User:Tohline/Appendix/Equation_templates#Torus|accompanying ''Equations'' appendix]].
==A.1==
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
   <td align="left">
   <td align="left">
Beginning with the identified "Key Equation",
<math>~\frac{\cosh\eta - 1}{\cosh\eta + 1}</math>
{{ User:Tohline/Math/EQ_Toroidal02 }}
  </td>
 
</tr>
we'll identify  <math>~x</math> with <math>~\cosh\eta</math> &#8212; in which case we have <math>~\lambda = \coth\eta</math> &#8212; and switch the index notations, <math>~n \leftrightarrow m</math>.  This gives,
 
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{m-1 / 2}^n (\coth\eta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,095: Line 2,148:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~(-1)^m
<math>~\tanh^2\biggl(\frac{\eta}{2}\biggr) \, .</math>
\frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(m-n+\frac{1}{2})} (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, .
</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>


Drawing upon the ''Euler reflection formula for gamma functions'', namely,
Note as well that,
{{ User:Tohline/Math/EQ_Gamma01 }}
<table border="0" cellpadding="5" align="center">


where it is understood that <math>~m</math> and <math>~n</math> are each either zero or a positive integer, this toroidal-function relation becomes,
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{m-1 / 2}^n (\coth\eta)</math>
<math>~1 + k_1</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,115: Line 2,164:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~(-1)^m~
<math>~1+ \frac{1 - \sqrt{1-k^2}}{1 + \sqrt{1-k^2}} = \frac{2}{1 + \sqrt{1-k^2}} \, .</math>
\frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{ \Gamma(n - m +\frac{1}{2}) }{ \pi (-1)^{m-n} } \biggr] (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta)  
  </td>
</tr>
</table>
 
Hence, the "second attempt" expression for <math>~P_{1 / 2}(z)</math> becomes,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~P_{-1 / 2}(\cosh\eta)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\sqrt{2}}{\pi}~ (\sinh\eta)^{-1 / 2} \biggl[ \frac{2k}{1 + \sqrt{1-k^2}} \biggr] K\biggl(\tanh \frac{\eta}{2} \biggr)  
</math>
</math>
  </td>
  <td align="right" width="15%">
&nbsp;
   </td>
   </td>
</tr>
</tr>
Line 2,133: Line 2,195:
   <td align="left">
   <td align="left">
<math>~
<math>~
(-1)^n \sqrt{ \frac{\pi}{2} } ~\Gamma(n - m + \tfrac{1}{2} )(\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, .
\frac{2^{2}}{\pi}~ (\sinh\eta)^{-1 / 2} \biggl[ \frac{1}{\coth\eta + 1}\biggr]^{1 / 2} \biggl\{ \frac{ [\coth\eta + 1]^{1 / 2} }{ [\coth\eta + 1]^{1 / 2} +  [\coth\eta - 1]^{1 / 2} } \biggr\}K\biggl(\tanh \frac{\eta}{2} \biggr)  
</math>
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2460;</font>
   </td>
   </td>
</tr>
</tr>
</table>


<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
   </td>
   </td>
</tr>
</table>
==A.2==
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
   <td align="left">
   <td align="left">
Again, beginning with the identified "Key Equation",
<math>~
{{ User:Tohline/Math/EQ_Toroidal02 }}
\frac{2^{2}}{\pi}~ \biggl[ ( \cosh\eta + \sinh\eta )^{1 / 2} +  (\cosh\eta - \sinh\eta)^{1 / 2}  \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr)
 
</math>
this time, without switching index notations, we'll identify  <math>~x</math> with <math>~\coth\eta</math> &#8212; in which case we have <math>~\lambda = \cosh\eta</math>.  This gives,
  </td>
 
</tr>
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,167: Line 2,222:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~(-1)^n
<math>~
\frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(n-m+\frac{1}{2})} \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta) \, .
\frac{2^{2}}{\pi}~ \biggl[ e^{\eta/2} +  e^{-\eta/2} \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr)  
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
Drawing upon the same ''Euler reflection formula for gamma functions'', as quoted above, this toroidal function relation can be rewritten as,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
&nbsp;
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,186: Line 2,236:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~(-1)^n
<math>~
\frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\pi (-1)^{m-n}} \biggr] \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta)
\biggl[ \frac{\pi}{2} \cdot \cosh\frac{\eta}{2} \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr) \, .
</math>
</math>
   </td>
   </td>
</tr>
</tr>
</table>
This is, indeed, identical to the "Attempt #1" expression for  <math>~P_{-1 / 2}(z)</math> that was used by Wong.
'''Q.E.D.'''
==Series Expansions==
From our [[User:Tohline/Appendix/Equation_templates#Complete_Elliptic_Integrals|accompanying set of Key mathematical relations]] we find the following series expansions for complete elliptic integrals of, respectively, the first and second kind:


<tr>
<div align="center">
  <td align="right">
{{ User:Tohline/Math/EQ_EllipticIntegral01 }}<br />
&nbsp;
{{ User:Tohline/Math/EQ_EllipticIntegral02 }}
  </td>
</div>
  <td align="center">
If you don't know how to handle a ''double factorial'' &#8212; also referred to as a ''semifactorial'' &#8212; see, for example, [https://en.wikipedia.org/wiki/Double_factorial this wikipedia page].
<math>~=</math>
  </td>
  <td align="left">
<math>~(-1)^{-m}~
\sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr]  P_{m-1 / 2}^n(\coth\eta) \, .
</math>
  </td>
</tr>
</table>
Finally, calling upon the "Key Equation" relation,


{{ User:Tohline/Math/EQ_Toroidal00 }}
=See Also=


making the index notation substitution, <math>~\nu \rightarrow (m-\tfrac{1}{2})</math>, and associating <math>~z</math> with <math>~ \coth\eta</math> gives,
<ul>
<li>[http://user.xmission.com/~rimrock/Documents/The%20Charged%20Bowl%20in%20Toroidal%20Coordinates.pdf P. Lucht (2016)], ''The Charged Bowl in Toroidal Coordinates'' &#8212; Note from Tohline:  On 3 July 2018 I stumbled on this article by Phil Lucht, who lists his affiliation as Rimrock Digital Technology (rimrock at emission.com), Salt Lake City, Utah. As I have done over the past approximately half a year, Lucht appears to have spent quite a bit of time investigating how toroidal coordinates might be used to solve a few specific potential problems.  He draws on many of the same, rich technical resource publications as I have done and, at one point, explicitly expresses amazement at the extensive amount of useful material that can be found in [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>].  Anyone digging into this topic would be well-advised to read Lucht's '''Overview and Summary''' &#8212; especially including the subsection titled, ''Legendre Functions'' &#8212; and to recognize that his document contains a large number of technically detailed appendices.
<br />&nbsp;<br />
To provide one point of comparison with [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], Lucht derived an analytic expression for the potential outside of a torus (with toroidal-coordinate radius, <math>~\eta_0</math>) whose surface is held at uniform potential, <math>~\Phi_0</math>.  This is in contrast to Wong's (1973) derivation &#8212; which we have verified, above &#8212; of the potential outside (as well as inside) of a torus having uniform charge/mass.  For his particular example problem, Lucht obtains,  
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~P^n_{m-\frac{1}{2}}(\coth\eta)</math>
<math>~\Phi(\eta,\theta)</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 2,221: Line 2,272:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\Gamma(m-n+\frac{1}{2})} \biggr]P^{-n}_{m-\frac{1}{2}}(\coth\eta)  \, .</math>
<math>~
\Phi_0 \frac{\sqrt{2}}{\pi} \sqrt{\cosh\eta - \cos\theta} \sum_{n=0}^\infty \epsilon_n P_{n-\frac{1}{2}}(\cosh\eta)
\Biggl[ \frac{Q_{n-\frac{1}{2}}(\cosh\eta_0)}{P_{n-\frac{1}{2}}(\cosh\eta_0)} \Biggr] \cos(n\theta) \, ,
</math>
   </td>
   </td>
</tr>
</tr>
</table>
As a result, we can write,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(-1)^{-m}~
\sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr]  P_{m-1 / 2}^{-n}(\coth\eta)
</math>
  </td>
  <td align="right" width="15%">
&nbsp;
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~
P_{m-1 / 2}^{-n}(\coth\eta)
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(-1)^m~\sqrt{\frac{2}{\pi}} \biggl[ \frac{\sqrt{\sinh\eta}}{\Gamma(m+n+\frac{1}{2})} \biggr] 
Q_{n-1 / 2}^m (\cosh\eta) \, .
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2461;</font>
  </td>
</tr>
</table>
  </td>
</tr>
</table>
==A.3==
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Here, our objective is to evaluate the definite integral,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\int_{-\pi}^{\pi} \frac{\cos[n(\theta - \theta^')] d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} </math>
  </td>
  <td align="right">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\cos(n\theta^') \int_{-\pi}^{\pi} \frac{ \cos(n\theta) ~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}}
+ \sin(n\theta^') \cancelto{0}{ \int_{-\pi}^{\pi} \frac{ \sin(n\theta)  d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} }
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2 \cos(n\theta^') \int_{0}^{\pi} \frac{ \cos(n\theta)~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} \, .
</math>
  </td>
</tr>
</table>
</div>
Notice that, since the limits of the integration are <math>~-\pi</math> to <math>~+\pi</math>: &nbsp; The second integral on the right-hand-side goes to zero because the numerator of its integrand &#8212; ''i.e.,'' <math>~\sin(n\theta)</math> &#8212; is an [https://en.wikipedia.org/wiki/Even_and_odd_functions#Odd_functions odd function]; and, with regard to the first integral on the right-hand-side, the lower integration limit can be set to zero and the result doubled because the numerator of its integrand &#8212; ''i.e.,'' <math>~\cos(n\theta)</math> &#8212; is an [https://en.wikipedia.org/wiki/Even_and_odd_functions#Even_functions even function].
Drawing from Volume I of [https://authors.library.caltech.edu/43491/1/Volume%201.pdf A. Erd&eacute;lyi's (1953)] ''Higher Transcendental Functions,'' we find the following "Key Equation":
{{ User:Tohline/Math/EQ_Toroidal03 }}
Next we adopt the associations, <math>~z \rightarrow \cosh\eta</math>, <math>~t \rightarrow \theta</math>, <math>~\mu \rightarrow 2</math>, and, <math>~\nu \rightarrow n - \tfrac{1}{2}</math>, where <math>~n</math> is zero or a positive integer.  In this case we have,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~Q_{n - \frac{1}{2}}^2 (\cosh\eta)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
(2\pi)^{-\frac{1}{2}} (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~\biggl\{
\int_0^\pi (\cosh\eta - \cos \theta)^{-\frac{5}{2}} \cos(n\theta) ~d\theta
- \cancelto{0}{\cos[(n-\tfrac{1}{2})\pi] }~~\int_0^\infty (\cosh\eta + \cosh \theta)^{- \frac{5}{2}} e^{-n\theta} ~d\theta
\biggr\} \, ,
</math>
  </td>
</tr>
</table>
where the prefactor of the second term &#8212; that is, <math>~\cos[(n-\tfrac{1}{2})\pi] </math> &#8212; goes to zero for all allowable values of the integer, <math>~n</math>.  Hence, we conclude that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~2\cos(n\theta^')
\int_0^\pi \frac{ \cos(n\theta)~d\theta }{ (\cosh\eta - \cos \theta)^{\frac{5}{2}} }
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left" colspan="2">
<math>~\frac{ 2(2\pi)^{\frac{1}{2}} Q_{n - \frac{1}{2}}^2 (\cosh\eta) \cos(n\theta^')}{ (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~ }
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \frac{ 2^3 \sqrt{2} }{ 3 } \biggr] \frac{  Q_{n - \frac{1}{2}}^2(\cosh\eta) \cos(n\theta^')}{ \sinh^2\eta  }
\, .
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2462;</font>
  </td>
</tr>
</table>
where we have set,
<div align="center">
<math>~
\Gamma(\tfrac{5}{2}) = \Gamma(\tfrac{1}{2} + 2) = \frac{ \sqrt{\pi} \cdot 4! }{4^2 \cdot 2!}
= \frac{\sqrt{\pi} \cdot 2^3\cdot 3}{  2^5 } = \frac{3 \sqrt{\pi}}{2^2} \, .
</math>
</div>
  </td>
</tr>
</table>
==A.4==
<table border="0" align="center" cellpadding="10" width="90%">
<tr>
  <td align="left">
Beginning with the identified "Key Equation",
{{ User:Tohline/Math/EQ_Toroidal05 }}
we will adopt the associations: &nbsp; <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math>, <math>~\nu \rightarrow (n - \tfrac{1}{2})</math>, <math>~\rho \rightarrow 0</math>, and <math>~\sigma \rightarrow ( n - \tfrac{1}{2})</math>.  As a result, [https://authors.library.caltech.edu/43491/1/Volume%201.pdf Erd&eacute;lyi's (1953)] expression becomes,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\int_a^b\biggl[ -4(1 - t^2)^{-1} \biggr] Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) ~dt
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[
(n - \tfrac{1}{2} ) Q_{n - \frac{1}{2}}^2(t) ~ X_{n - \frac{3}{2}}(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t)
\biggr]_a^b \, .
</math>
  </td>
</tr>
</table>
Drawing upon the ''recurrence'' "Key Equation,"
{{ User:Tohline/Math/EQ_Toroidal04 }}
which means, after making the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 0</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, that,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~(n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2n t X_{n-\frac{1}{2}}(t) - (n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t)</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~(n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t) </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>
  </td>
</tr>
</table>
the integral can be rewritten as,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
\int_a^b\biggl[  \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{4} \biggl\{
\biggl[ 2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) \biggr]  Q_{n - \frac{1}{2}}^2(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t)
\biggr\}_a^b
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{4} \biggl\{
\biggl[ 2n t  Q_{n - \frac{1}{2}}^2(t)
- (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) \biggr] X_{n - \frac{1}{2}}(t)
-  (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)
\biggr\}_a^b \, .
</math>
  </td>
</tr>
</table>
Returning to the same ''recurrence'' "Key Equation," but this time adopting the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, we can write,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2n t Q^2_{n - \frac{1}{2}}(t) - (n + \tfrac{3}{2}) Q^2_{n - \frac{3}{2}} (t) \, ,
</math>
  </td>
</tr>
</table>
in which case the integral becomes,
<table border="0" cellpadding="5" align="center">


<tr>
<tr>
   <td align="right">
   <td align="center" colspan="3">
<math>~
[http://user.xmission.com/~rimrock/Documents/The%20Charged%20Bowl%20in%20Toroidal%20Coordinates.pdf Lucht (2016)], p. 68, Eq. (10.1.11)
\int_a^b\biggl[ \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{4} \biggl\{
(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ X_{n - \frac{1}{2}}(t)
-  (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t)
\biggr\}_a^b \, .
</math>
  </td>
  <td align="right" width="15%">
<font color="green" size="+3">&#x2463;</font>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
then he points out that the associated expression for the potential, which appears as equation (10.3.80) on p. 1304 of [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>], is missing the Neumann factor, <math>~\epsilon_n</math>.
</li>
</ul>


</td>
</tr>
</table>
=See Also=
* [http://user.xmission.com/~rimrock/Documents/The%20Charged%20Bowl%20in%20Toroidal%20Coordinates.pdf P. Lucht (2016)], ''The Charged Bowl in Toroidal Coordinates'' &#8212; Note from Tohline:  On 3 July 2018 I stumbled on this article by Phil Lucht, who lists his affiliation as Rimrock Digital Technology (rimrock at emission.com), Salt Lake City, Utah. As I have done over the past approximately half a year, Lucht appears to have spent quite a bit of time investigating how toroidal coordinates might be used to solve a "potential" problem.  He draws on many of the same, rich technical resource publications as I have done and, at one point, explicitly expresses amazement at the extensive amount of useful material that can be found in [<b>[[User:Tohline/Appendix/References#MF53|<font color="red">MF53</font>]]</b>].  Anyone digging into this topic would be well-advised to read Lucht's '''Overview and Summary''' and to recognize that his document contains a large number of technically detailed appendices.




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 22:48, 15 July 2020

Using Toroidal Coordinates to Determine the Gravitational Potential

NOTE:   An earlier version of this chapter has been shifted to our "Ramblings" Appendix.

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |


Here we build upon our accompanying review of the types of numerical techniques that various astrophysics research groups have developed to solve for the Newtonian gravitational potential, <math>~\Phi(\vec{x})</math>, given a specified, three-dimensional mass distribution, <math>~\rho(\vec{x})</math>. Our focus is on the use of toroidal coordinates to solve the integral formulation of the Poisson equation, namely,

<math>~ \Phi(\vec{x})</math>

<math>~=</math>

<math>~ -G \iiint \frac{\rho(\vec{x}^{~'})}{|\vec{x}^{~'} - \vec{x}|} d^3x^' \, .</math>

For the most part, we will adopt the notation used by C.-Y. Wong (1973, Annals of Physics, 77, 279); in an accompanying discussion, we review additional results from this insightful 1973 paper, as well as a paper of his that was published the following year in The Astrophysical Journal, namely, Wong (1974).

In order to accomplish this task, we first present the expressions that define how toroidal coordinates, <math>~(\eta,\theta,\psi)</math>, map to and from Cartesian coordinates <math>~(x, y, z)</math>, and present the toroidal-coordinate expression for the differential volume element, <math>~d^3 x</math>.

Basic Elements of a Toroidal Coordinate System

Given the meridional-plane coordinate location of a toroidal-coordinate system's axisymmetric anchor ring, <math>~(\varpi,z) = (a,Z_0)</math>, the relationship between toroidal coordinates <math>~(\eta,\theta,\psi) </math>and Cartesian coordinates <math>~(x, y, z)</math> is,

<math>~x</math>

<math>~=</math>

<math>~\frac{a \sinh\eta \cos\psi}{(\cosh\eta - \cos\theta)} \, ,</math>

<math>~y</math>

<math>~=</math>

<math>~\frac{a \sinh\eta \sin\psi}{(\cosh\eta - \cos\theta)} \, ,</math>

<math>~z - Z_0</math>

<math>~=</math>

<math>~\frac{a \sin\theta}{(\cosh\eta - \cos\theta)} \, .</math>

This set of coordinate relations appears as equations 2.1 - 2.3 in Wong (1973). This set of relations may also be found, for example, on p. 1301 within eq. (10.3.75) of [MF53]; in §14.19 of NIST's Digital Library of Mathematical Functions; or even within Wikipedia. (In most cases the implicit assumption is that <math>~Z_0 = 0</math>.) It is clear, of course, that the cylindrical radial coordinate is,

<math>~\varpi = (x^2 + y^2)^{1 / 2}</math>

<math>~=</math>

<math>~\frac{a \sinh\eta}{(\cosh\eta - \cos\theta)} \, .</math>


Mapping the other direction [see equations 2.13 - 2.15 of Wong (1973) ], we have,

<math>~\eta</math>

<math>~=</math>

<math>~\ln\biggl(\frac{r_1}{r_2} \biggr) \, ,</math>

<math>~\cos\theta</math>

<math>~=</math>

<math>~\frac{(r_1^2 + r_2^2 - 4a^2)}{2r_1 r_2} \, ,</math>

<math>~\tan\psi</math>

<math>~=</math>

<math>~\frac{y}{x} \, ,</math>

where,

<math>~r_1^2 </math>

<math>~\equiv</math>

<math>~[(x^2 + y^2)^{1 / 2} + a]^2 + (z-Z_0)^2 \, ,</math>

<math>~r_2^2 </math>

<math>~\equiv</math>

<math>~[(x^2 + y^2)^{1 / 2} - a]^2 + (z-Z_0)^2 \, ,</math>

and <math>~\theta</math> has the same sign as <math>~(z-Z_0)</math>.

Title Page of Hicks (1881)

W. M. Hicks (1881) presents one of the first — if not the first — discussions of toroidal coordinates and associated toroidal functions. Equation (4) on p. 614 of his article provides the following definition of the pair of meridional-plane coordinates <math>~(u,v)</math>, written in terms of the traditional cylindrical-coordinate pair <math>~(\rho,z)</math> and the specified anchor ring radius, <math>~a</math>:

<math>~u</math>

<math>~=</math>

<math>~ \frac{1}{2} \ln \biggl[ \frac{z^2 + (\rho+a)^2}{z^2 + (\rho - a)^2} \biggr] \, , </math>

<math>~v</math>

<math>~=</math>

<math>~ \tan^{-1} \biggl[ \frac{2az}{\rho^2 + z^2 - a^2} \biggr] \, . </math>

His equation (5) presents the reverse mapping, namely,

<math>~\rho</math>

<math>~=</math>

<math>~ \frac{a\sinh u}{\cosh u - \cos v} \, , </math>

<math>~z</math>

<math>~=</math>

<math>~ \frac{a\sin v}{\cosh u - \cos v} \, . </math>

 

Comment by J. E. Tohline on 16 August 2017: In equation (2.17) of his §IIB — when Wong (1973) introduces the differential volume element — the variable used to represent the azimuthal coordinate angle switches from ψ to Φ. We will stick with the ψ notation, here.

According to p. 1301, eq. (10.3.75) of [MF53] — or, for example, as found in Wikipedia — in toroidal coordinates the differential volume element is,

<math>~d^3x</math>

<math>~=</math>

<math>~h_\eta h_\theta h_\psi d\eta d\theta d\psi</math>

<math>~=</math>

<math>~\biggl[ \frac{a^3 \sinh\eta}{(\cosh\eta - \cos\theta)^3} \biggr] d\eta~ d\theta~ d\psi \, .</math>

Green's Function Expression

As presented by Wong (1973)

Referencing [MF53], Wong (1973) states that, in toroidal coordinates, the Green's function is,

<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>

<math>~=</math>

<math>~ \frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 } \sum\limits^\infty_{m,n=0} (-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} </math>

 

 

<math>~ \times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~\begin{cases}P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^') ~~~\eta^' > \eta \\P^m_{n-1 / 2}(\cosh\eta^') ~Q^m_{n-1 / 2}(\cosh\eta)~~~\eta^' < \eta \end{cases}\, , </math>

Wong (1973), p. 293, Eq. (2.53)
[see also: J. W. Bates (1997), p. 3685, Eq. (31)]
[see also: Cohl, Tohline, Rau, & Srivastava (2000), §6.2, Eq. (48)]

where, <math>~P^m_{n-1 / 2}, Q^m_{n-1 / 2}</math> are Associated Legendre Functions of the first and second kind with degree <math>~n - \tfrac{1}{2}</math> and order <math>~m</math> (toroidal harmonics), and <math>~\epsilon_m</math> is the Neumann factor, that is, <math>~\epsilon_0 = 1</math> and <math>~\epsilon_m = 2</math> for all <math>~m \ge 1</math>. This Green's function expression can indeed be found as eq. (10.3.81) on p. 1304 of [MF53], but it should be noted that the MF53 expression differs from Wong's in two respects (see footnote 2 on p. 370 of Cohl et al. (2000) for a proposed explanation): First, the factor, <math>~(-1)^m</math>, appears as <math>~(-i)^m</math> in MF53; and, second, in the term that is composed of a ratio of gamma functions, the denominator appears in MF53 as <math>~\Gamma(n - m + \tfrac{1}{2})</math>, whereas it should be <math>~\Gamma(n + m + \tfrac{1}{2})</math>, as presented here.

As Presented by Hicks (1881)

Notably, more than 135 years ago, W. M. Hicks (1881) already had constructed a Green's function expression for the reciprocal distance between two points in toroidal coordinates. The following boxed-in table contains a snapshot image of equation (35) from Hicks (1881), along with its associated text; notice that he refers to <math>~\phi</math> as the "potential." Although notations are different — for example, <math>~C</math> is shorthand for <math>~\cosh u</math> and <math>~c</math> is shorthand for <math>~\cos v</math> — one can easily see factor-by-factor agreement when compared with the Green's function presented by Wong (1973).

Title Page of Hicks (1881)

Green's Function from Hicks (1881)

Rearranging Terms and Incorporating Special-Function Relations

Let's focus on the situation when <math>~\eta^' > \eta</math>, and begin rearranging or substituting terms.

<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>

<math>~=</math>

<math>~ \frac{1}{\pi a} \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 } \sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] </math>

 

 

<math>~ \times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} ~P^m_{n-1 / 2}(\cosh\eta) ~Q^m_{n-1 / 2}(\cosh\eta^') </math>

 

<math>~=</math>

<math>~ \frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} } \sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] </math>

 

 

<math>~ \times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~ \sqrt{ \frac{\pi}{2} }~\Gamma(n-m+\tfrac{1}{2}) \sqrt{\sinh\eta}~P^m_{n-1 / 2}(\cosh\eta) \biggl\}\biggr\{ ~\sqrt{ \frac{2}{\pi} }~\frac{\sqrt{\sinh\eta^'}}{\Gamma(n + m + \tfrac{1}{2})} Q^m_{n-1 / 2}(\cosh\eta^') \biggr\} </math>

The term contained within the first set of curly braces on the right-hand side of this expression can be replaced by the derived expression labeled in the Appendix, below, and simultaneously the term contained within the second set of curly braces can be replaced by the derived expression labeled in the same Appendix. After making these substitutions, we have,

<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>

<math>~=</math>

<math>~ \frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} } \sum\limits^\infty_{m=0} (-1)^m \epsilon_m \cos[m(\psi - \psi^')] </math>

 

 

<math>~ \times \sum\limits^\infty_{n=0} \epsilon_n \cos[n(\theta - \theta^')] \biggl\{ ~ (-1)^{-n}Q^n_{m-1 / 2}(\coth\eta) \biggl\}\biggr\{ ~(-1)^{-m} P^{-n}_{m-1 / 2}(\coth\eta^') \biggr\} </math>

 

<math>~=</math>

<math>~ \frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} } \sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')] </math>

 

 

<math>~ \times \sum\limits^\infty_{n=0} (-1)^{n} \epsilon_n \cos[n(\theta - \theta^')] Q^n_{m-1 / 2}(\coth\eta) P^{-n}_{m-1 / 2}(\coth\eta^') \, , </math>

where, in writing this last expression we have acknowledged that, since <math>~n</math> is either zero or a positive integer, <math>~(-1)^{-n} = (-1)^n</math>. Next we draw upon the "Key Equation" relation,

LSU Key.png

<math>~ Q_\nu[t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos\psi] </math>

<math>~=</math>

<math>~ Q_\nu(t) P_\nu(t^') + 2\sum_{n=1}^\infty (-1)^n Q^n_\nu(t) P^{-n}_\nu(t^') \cos(n\psi) </math>

A. Erdélyi (1953):  Volume I, §3.11, p. 169, eq. (4)

Valid for:    

<math>~t, t^'</math>  real

       

<math>~1 < t^' < t</math>

       

<math>~\nu \ne -1, -2, -3, </math> …

       

<math>~\psi</math>   real

which, after making the substitutions, <math>~\nu \rightarrow (m - \tfrac{1}{2})</math> and <math>~\psi \rightarrow (\theta - \theta^')</math>, and incorporating the Neumann factor, <math>~\epsilon_n</math>, becomes,

<math>~ Q_{m - \frac{1}{2} }\ [t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^') ] </math>

<math>~=</math>

<math>~ \sum_{n=0}^\infty (-1)^n \epsilon_n Q^n_{m - \frac{1}{2} }(t) P^{-n}_{m - \frac{1}{2} }(t^') \cos[(n(\theta- \theta^')] \, . </math>

Finally, after making the associations, <math>~t \rightarrow \coth\eta</math> and <math>~t^' \rightarrow \coth\eta^'</math>, this last expression allows us to rewrite Wong's (1973) Green's function in a significantly more compact form, namely,

<math>~\frac{1}{|~\vec{x} - {\vec{x}}^{~'} ~|} </math>

<math>~=</math>

<math>~ \frac{ [ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^')]^{1 /2 } }{\pi a \sqrt{\sinh\eta^'} \sqrt{\sinh\eta} } \sum\limits^\infty_{m=0} \epsilon_m \cos[m(\psi - \psi^')] Q_{m - \frac{1}{2}}(\Chi) \, , </math>

where the argument, <math>~\Chi</math>, of the toroidal function, <math>~Q_{m - \frac{1}{2}}</math>, is,

<math>~\Chi</math>

<math>~\equiv</math>

<math>~ t t^' - (t^2-1)^{1 / 2} (t^{'2} - 1)^{1 / 2} \cos(\theta- \theta^') </math>

 

<math>~=</math>

<math>~ \coth\eta \coth\eta^' - (\coth^2\eta-1)^{1 / 2} (\coth^2\eta'- 1)^{1 / 2} \cos(\theta- \theta^') </math>

 

<math>~=</math>

<math>~ \frac{\cosh\eta \cosh\eta^'}{\sinh\eta \sinh\eta^'} - \biggl[ \frac{1}{\sinh^2\eta} \biggr]^{1 / 2} \biggl[ \frac{1}{\sinh^2\eta'}\biggr]^{1 / 2} \cos(\theta- \theta^') </math>

 

<math>~=</math>

<math>~ \frac{\cosh\eta \cosh\eta^' - \cos(\theta- \theta^') }{\sinh\eta \sinh\eta^'} \, . </math>

As Presented in Cohl & Tohline (1999)

This last, compact Green's function expression — which we have derived, here, from Wong's (1973) published Green's function by drawing strategically upon a variety of special function relations — precisely matches the "compact cylindrical Green's function expression" that has been derived independently by Cohl & Tohline (1999) via a less tortuous route, namely,

<math>~ \frac{1}{|\vec{x} - \vec{x}^{~'}|}</math>

<math>~=</math>

<math>~ \frac{1}{\pi \sqrt{\varpi \varpi^'}} \sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi) </math>

 

Cohl & Tohline (1999), p. 88, Eq. (17)

 

<math>~=</math>

<math>~ \frac{1}{a\pi} \biggl[ \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi) \, , </math>

where,

<math>~\Chi</math>

<math>~\equiv</math>

<math>~ \frac{(\varpi^')^2 + \varpi^2 + (z^' - z)^2}{2\varpi^' \varpi} </math>

<math>~=</math>

<math>~ \frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, . </math>

   

Cohl & Tohline (1999), p. 88, Eq. (16)


Note from J. E. Tohline (June, 2018):  This is the first time that I have been able to formally demonstrate to myself that these two separately derived Green's function expressions are identical. See, however, the earlier identification of new addition theorems in association with equations (49) and (50) of Cohl et al. (2000).

Gravitational Potential

Quite generally, then, the gravitational potential can be obtained at any coordinate location, <math>~(\eta,\theta,\psi)</math> — both inside and outside of a specified mass distribution — by carrying out three nested spatial integrals over the product of:  <math>~\rho(\vec{x}^{~'})</math>, the differential volume element, and the Green's function as specified either by Wong (1973) or by Cohl & Tohline (1999).

In what follows we will make an effort to elucidate the pros and cons of adopting one Green's function expression over the other. In each case we begin by writing the expression for the potential in such a way that variations in the azimuthal coordinate, <math>~\psi</math>, are described by Fourier components, <math>~\Phi_m^{(1)}(\eta,\theta)</math> and <math>~\Phi_m^{(2)}(\eta,\theta)</math>, of the potential, such that,

<math>~\Phi(\vec{x}) = \tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \Phi_m^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \sin (m\psi) \Phi_m^{(2)}(\eta,\theta) \, .</math>

Each Fourier component of the potential depends explicitly on the corresponding Fourier component of the density distribution, defined such that,

<math>~\rho(\vec{x}) = \tfrac{1}{2}\rho_0^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \cos (m\psi) \rho_m^{(1)}(\eta,\theta) + \sum_{m=1}^\infty \sin (m\psi) \rho_m^{(2)}(\eta,\theta) \, .</math>

LaTeX mathematical expressions cut-and-pasted directly from
NIST's Digital Library of Mathematical Functions

As an additional primary point of reference, note that according to §1.8(i) of NIST's Digital Library of Mathematical Functions, a Fourier Series is defined as follows:

<math>~f(x)</math>

<math>~=</math>

<math>~\tfrac{1}{2}a_{0}+\sum^{\infty}_{n=1}\biggl[ a_{n}\cos\bigl(nx\bigr)+b_{n}\sin\bigl(nx\bigr) \biggr],</math>

<math>~a_n</math>

<math>~=</math>

<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\cos\bigl(nx\bigr)\mathrm{d}x,</math>

<math>~b_n</math>

<math>~=</math>

<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}f(x)\sin\bigl(nx\bigr)\mathrm{d}x.</math>

Notice, therefore, that,

<math>~\rho_m^{(1)}(\eta,\theta)</math>

<math>~=</math>

<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\cos\bigl(m\psi\bigr)\mathrm{d}\psi,</math>

    and,    

<math>~\rho_m^{(2)}(\eta,\theta)</math>

<math>~=</math>

<math>~\frac{1}{\pi}\int^{\pi}_{-\pi}\rho(\eta,\theta,\psi)\sin\bigl(m\psi\bigr)\mathrm{d}\psi \, .</math>


The CT99 Expression for the Potential

In Three-Dimensional Generality

Employing the Green's function expression derived by Cohl & Tohline (1999), the gravitational potential for any three-dimensional matter distribution is,

<math>~ \Phi(\eta,\theta,\psi)</math>

<math>~=</math>

<math>~ -G \iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\} \biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^' </math>

 

<math>~=</math>

<math>~ -\frac{Ga^2}{\pi} \int d\eta^' \int d\theta^' \int d\psi^' \iiint \rho(\eta^',\theta^',\psi^') \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] </math>

 

 

<math>~ \times \biggl[ \frac{(\cosh\eta^' - \cos\theta^')}{\sinh\eta^' } \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \sum_{m=0}^{\infty} \epsilon_m \cos[m(\psi - \psi^')] Q_{m- 1 / 2}(\Chi) </math>

 

<math>~=</math>

<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \sum_{m=0}^{\infty} \epsilon_m \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \int d\psi^' \rho(\eta^',\theta^',\psi^') </math>

 

 

<math>~ \times [\cos(m\psi)\cos(m\psi^') + \sin(m\psi)\sin(m\psi^') ] </math>

 

<math>~=</math>

<math>~ -\frac{Ga^2}{\pi} \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{ \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi) \int_{-\pi}^{\pi} d\psi^' \rho(\eta^',\theta^',\psi^') </math>

 

 

<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi) \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \cos(m\psi^') </math>

 

 

<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi) \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^') \biggr\} </math>

 

<math>~=</math>

<math>~ - Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \biggl\{ \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- \frac{1}{2}}(\Chi) \rho_0^{(1)}(\eta^',\theta^') </math>

 

 

<math>~ + \sum_{m=1}^{\infty} 2\cos(m\psi) \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1)}(\eta^',\theta^') </math>

 

 

<math>~ + \sum_{m=1}^{\infty} 2\sin(m\psi) \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(2)}(\eta^',\theta^') \biggr\} \, . </math>

We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,

<math>~\Phi_m^{(1),(2)}(\eta,\theta)</math>

<math>~=</math>

<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{m- 1 / 2}(\Chi) \rho_m^{(1),(2)}(\eta^',\theta^') \, , </math>

Cohl & Tohline (1999), p. 88, Eq. (20)

where, as above,

<math>~\Chi</math>

<math>~\equiv</math>

<math>~ \frac{\cosh\eta \cdot \cosh\eta^' - \cos(\theta^' - \theta) }{ \sinh\eta \cdot \sinh\eta^'} \, . </math>

For Axisymmetric Systems

For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>. Hence, for all <math>~m > 0</math>, the Fourier components of the density, <math>~\rho_m^{(1),(2)}</math>, are zero. The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>. For axisymmetric systems, then, the gravitational potential is,

<math>~\Phi(\eta,\theta)</math>

<math>~=</math>

<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>

 

<math>~=</math>

<math>~ - 2Ga^2 \biggl[ \frac{(\cosh\eta - \cos\theta)}{\sinh\eta } \biggr]^{1 / 2} \int d\eta^' \int d\theta^' \biggl[ \frac{\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^5} \biggr]^{1 / 2} Q_{- 1 / 2}(\Chi) \rho(\eta^',\theta^') \, . </math>

Cohl & Tohline (1999), p. 88, Eqs. (31) & (32a)

Wong's Expression for the Potential

Fully Three-Dimensional Case

Employing Wong's (1973) Green's function expression, the gravitational potential for any three-dimensional matter distribution is,

<math>~ \Phi(\vec{x})</math>

<math>~=</math>

<math>~ -G \iiint \rho(\eta^',\theta^',\psi^') \biggl\{ \frac{1}{|\vec{x}^{~'} - \vec{x}|} \biggr\} \biggl[ \frac{a^3 \sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] d\eta^'~ d\theta^'~ d\psi^' </math>

 

<math>~=</math>

<math>~ - \frac{a^2G}{\pi} \int d\eta^' \int d\theta^' \int d\psi^' \biggl[ \frac{\rho(\eta^',\theta^',\psi^') ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^3} \biggr] \biggl[ (\cosh\eta - \cos\theta)(\cosh \eta^' - \cos\theta^') \biggr]^{1 /2 } \sum\limits^\infty_{n=0} \sum\limits^\infty_{m=0}(-1)^m \epsilon_m \epsilon_n ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} </math>

 

 

<math>~ \times \cos[m(\psi - \psi^')]\cos[n(\theta - \theta^')] ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) </math>

 

<math>~=</math>

<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \sum\limits^\infty_{m=0}(-1)^m \epsilon_m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} \int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) </math>

 

 

<math>~ \times~ \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \int d\psi^' \rho(\eta^',\theta^',\psi^') [\cos(m\psi)\cos(m\psi^') + \sin(m\psi) \sin(m\psi^')] </math>

 

<math>~=</math>

<math>~ - \frac{a^2G}{\pi} (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{ \int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') </math>

 

 

<math>~ +~ \sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} \int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \cos(m\psi^') </math>

 

 

<math>~ +~ \sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} \int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \int_{-\pi}^\pi d\psi^' \rho(\eta^',\theta^',\psi^') \sin(m\psi^') \biggr\} </math>

 

<math>~=</math>

<math>~ - a^2G (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \biggl\{ \int d\eta^' ~P_{n-1 / 2}(\cosh\eta_<) ~Q_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \rho_0^{(1)}(\eta^',\theta^') </math>

 

 

<math>~ +~ \sum\limits^\infty_{m=1} 2\cos(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} \int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \rho_m^{(1)}(\eta^',\theta^') </math>

 

 

<math>~ +~ \sum\limits^\infty_{m=1} 2 \sin(m\psi)(-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} \int d\eta^' ~P^m_{n-1 / 2}(\cosh\eta_<) ~Q^m_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl[ \frac{ ~\sinh\eta^'}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr]\cos[n(\theta - \theta^')] \rho_m^{(2)} (\eta^',\theta^') \biggr\} \, . </math>

We conclude, therefore, that each one of the Fourier components of the gravitational potential is given by the expression,

<math>~\Phi_m^{(1),(2)} (\eta,\theta)</math>

<math>~=</math>

<math>~ - 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n (-1)^m ~\frac{\Gamma(n-m+\tfrac{1}{2})}{\Gamma(n + m + \tfrac{1}{2})} </math>

 

 

<math>~\times \int d\eta^' ~\sinh\eta^' ~ P^m_{n-1 / 2}(\cosh\eta_<) ~ Q^m_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl\{ \frac{ \cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\} \rho_m^{(1),(2)}(\eta^',\theta^') \, . </math>

Axisymmetric Systems

For axisymmetric systems, the density distribution has no dependence on the azimuthal coordinate, <math>~\psi</math>. Hence, for all <math>~m > 0</math>, the Fourier components of the density, <math>~\rho_m^{(1),(2)}</math>, are zero. The only nonzero component is, <math>~\rho_0^{(1)}(\eta,\theta) = 2\rho(\eta,\theta)</math>. For axisymmetric systems, then, the gravitational potential is,

<math>~\Phi(\eta,\theta)</math>

<math>~=</math>

<math>~\tfrac{1}{2}\Phi_0^{(1)}(\eta,\theta)</math>

 

<math>~=</math>

<math>~ - 2Ga^2 (\cosh\eta - \cos\theta)^{1 / 2} \sum\limits^\infty_{n=0} \epsilon_n \int d\eta^' ~\sinh\eta^'~P^0_{n-1 / 2}(\cosh\eta_<) ~Q^0_{n-1 / 2}(\cosh\eta_>) \int d\theta^' \biggl\{ \frac{ ~\cos[n(\theta - \theta^')]}{(\cosh\eta^' - \cos\theta^')^{5/2}} \biggr\} \rho(\eta^',\theta^') \, . </math>

Wong (1973), p. 293, Eq. (2.55)

Uniform-Density Torus

In an accompanying discussion, we build upon the above technical foundation and detail how Wong (1973) was able to complete both integrals to derive an analytic expression for the potential (inside as well as outside) of axisymmetric, uniform-density tori having an arbitrarily specified ratio of the major to minor (cross-sectional) radii, <math>~R/d</math>. This is an outstanding accomplishment that has received little attention in the astrophysics literature and, therefore, has heretofore been underappreciated.

Appendix A: Selected Toroidal Function Relationships

Here, we draw from the set of toroidal function relationships that have been identified as "Key Equations" in our accompanying Equations appendix. Two related appendices — numbered A.3 and A.4 — may be found in an accompanying discussion.

A.1

Beginning with the identified "Key Equation",

LSU Key.png

<math>~Q_{n-1 / 2}^m (\lambda)</math>

<math>~=</math>

<math>~(-1)^n \frac{\pi^{3/2}}{\sqrt{2}~ \Gamma(n-m+1 / 2)} (x^2-1)^{1 / 4} P_{m-1 / 2}^n(x) \, , </math>

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

where:    

<math>~\lambda \equiv x/\sqrt{x^2-1}</math>

we'll identify <math>~x</math> with <math>~\cosh\eta</math> — in which case we have <math>~\lambda = \coth\eta</math> — and switch the index notations, <math>~n \leftrightarrow m</math>. This gives,

<math>~Q_{m-1 / 2}^n (\coth\eta)</math>

<math>~=</math>

<math>~(-1)^m \frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(m-n+\frac{1}{2})} (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, . </math>

Drawing upon the Euler reflection formula for gamma functions, namely,

LSU Key.png

<math>~ \Gamma(z) ~\Gamma(1-z) </math>

<math>~=</math>

<math>~ \frac{\pi}{\sin(\pi z)} </math>

<math>~\biggl|</math>

for example, if
<math>~z \rightarrow (m-n + \tfrac{1}{2})</math>

<math>~\Rightarrow ~~~\Gamma(m-n+\tfrac{1}{2})~\Gamma(n-m+\tfrac{1}{2})</math>

<math>~=</math>

<math>~\pi \biggl\{\sin\biggl[ \frac{\pi}{2} + \pi(m-n) \biggr] \biggr\}^{-1}</math>

 

<math>~=</math>

<math>~\pi (-1)^{m-n} </math>

DLMF §5.5(ii)

<math>~\biggl|</math>
Valid for:

   <math>~z \ne0, \pm 1, \pm 2, </math> …

<math>~\biggl|</math>

where it is understood that <math>~m</math> and <math>~n</math> are each either zero or a positive integer, this toroidal-function relation becomes,

<math>~Q_{m-1 / 2}^n (\coth\eta)</math>

<math>~=</math>

<math>~(-1)^m~ \frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{ \Gamma(n - m +\frac{1}{2}) }{ \pi (-1)^{m-n} } \biggr] (\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) </math>

 

 

<math>~=</math>

<math>~ (-1)^n \sqrt{ \frac{\pi}{2} } ~\Gamma(n - m + \tfrac{1}{2} )(\sinh\eta)^{1 / 2} P_{n-1 / 2}^m(\cosh\eta) \, . </math>


A.2

Again, beginning with the identified "Key Equation",

LSU Key.png

<math>~Q_{n-1 / 2}^m (\lambda)</math>

<math>~=</math>

<math>~(-1)^n \frac{\pi^{3/2}}{\sqrt{2}~ \Gamma(n-m+1 / 2)} (x^2-1)^{1 / 4} P_{m-1 / 2}^n(x) \, , </math>

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

where:    

<math>~\lambda \equiv x/\sqrt{x^2-1}</math>

this time, without switching index notations, we'll identify <math>~x</math> with <math>~\coth\eta</math> — in which case we have <math>~\lambda = \cosh\eta</math>. This gives,

<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>

<math>~=</math>

<math>~(-1)^n \frac{\pi^{3/2}}{\sqrt{2} ~\Gamma(n-m+\frac{1}{2})} \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta) \, . </math>

Drawing upon the same Euler reflection formula for gamma functions, as quoted above, this toroidal function relation can be rewritten as,

<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>

<math>~=</math>

<math>~(-1)^n \frac{\pi^{3/2}}{\sqrt{2} } \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\pi (-1)^{m-n}} \biggr] \biggl( \frac{1}{\sinh\eta} \biggr)^{1 / 2} P_{m-1 / 2}^n(\coth\eta) </math>

 

<math>~=</math>

<math>~(-1)^{-m}~ \sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m-n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr] P_{m-1 / 2}^n(\coth\eta) \, . </math>

Finally, calling upon the "Key Equation" relation,

LSU Key.png

<math>~ P_\nu^n(z) </math>

<math>~=</math>

<math>~ \frac{\Gamma(\nu + n + 1)}{\Gamma(\nu - n + 1)} P_\nu^{-n}(z) </math>

A. Erdélyi (1953):  Volume I, §3.3.1, p. 140, eq. (7)

making the index notation substitution, <math>~\nu \rightarrow (m-\tfrac{1}{2})</math>, and associating <math>~z</math> with <math>~ \coth\eta</math> gives,

<math>~P^n_{m-\frac{1}{2}}(\coth\eta)</math>

<math>~=</math>

<math>~\biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\Gamma(m-n+\frac{1}{2})} \biggr]P^{-n}_{m-\frac{1}{2}}(\coth\eta) \, .</math>

As a result, we can write,

<math>~Q_{n-1 / 2}^m (\cosh\eta)</math>

<math>~=</math>

<math>~(-1)^{-m}~ \sqrt{\frac{\pi}{2}} \biggl[ \frac{\Gamma(m+n+\frac{1}{2})}{\sqrt{\sinh\eta}} \biggr] P_{m-1 / 2}^{-n}(\coth\eta) </math>

 

<math>~\Rightarrow ~~~ P_{m-1 / 2}^{-n}(\coth\eta) </math>

<math>~=</math>

<math>~(-1)^m~\sqrt{\frac{2}{\pi}} \biggl[ \frac{\sqrt{\sinh\eta}}{\Gamma(m+n+\frac{1}{2})} \biggr] Q_{n-1 / 2}^m (\cosh\eta) \, . </math>


A.3

Here, our objective is to evaluate the definite integral,

<math>~\int_{-\pi}^{\pi} \frac{\cos[n(\theta - \theta^')] d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} </math>

<math>~=</math>

<math>~ \cos(n\theta^') \int_{-\pi}^{\pi} \frac{ \cos(n\theta) ~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} + \sin(n\theta^') \cancelto{0}{ \int_{-\pi}^{\pi} \frac{ \sin(n\theta) d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} } </math>

 

<math>~=</math>

<math>~ 2 \cos(n\theta^') \int_{0}^{\pi} \frac{ \cos(n\theta)~ d\theta}{(\cosh\eta - \cos\theta)^{5 / 2}} \, . </math>

Notice that, since the limits of the integration are <math>~-\pi</math> to <math>~+\pi</math>:   The second integral on the right-hand-side goes to zero because the numerator of its integrand — i.e., <math>~\sin(n\theta)</math> — is an odd function; and, with regard to the first integral on the right-hand-side, the lower integration limit can be set to zero and the result doubled because the numerator of its integrand — i.e., <math>~\cos(n\theta)</math> — is an even function.

Drawing from Volume I of A. Erdélyi's (1953) Higher Transcendental Functions, we find the following "Key Equation":

LSU Key.png

<math>~Q_\nu^\mu(z)</math>

<math>~=</math>

<math>~ e^{i \mu \pi} ~ (2\pi)^{-\frac{1}{2}} (z^2-1)^{\mu/2} ~\Gamma(\mu + \tfrac{1}{2})~\biggl\{ \int_0^\pi (z - \cos t)^{-\mu - \frac{1}{2}} \cos[(\nu + \tfrac{1}{2})t] ~dt -\cos(\nu\pi) \int_0^\infty (z + \cosh t)^{-\mu - \frac{1}{2}} e^{-(\nu + \frac{1}{2})t} ~dt \biggr\} </math>

A. Erdélyi (1953):  Volume I, §3.7, p. 156, eq. (10)

Valid for:    

<math>~\mathrm{Re} ~\nu > -\tfrac{1}{2}</math> 

    and    

<math>~\mathrm{Re} (\nu + \mu + 1) > 0 \, .</math>

Next we adopt the associations, <math>~z \rightarrow \cosh\eta</math>, <math>~t \rightarrow \theta</math>, <math>~\mu \rightarrow 2</math>, and, <math>~\nu \rightarrow n - \tfrac{1}{2}</math>, where <math>~n</math> is zero or a positive integer. In this case we have,

<math>~Q_{n - \frac{1}{2}}^2 (\cosh\eta)</math>

<math>~=</math>

<math>~ (2\pi)^{-\frac{1}{2}} (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~\biggl\{ \int_0^\pi (\cosh\eta - \cos \theta)^{-\frac{5}{2}} \cos(n\theta) ~d\theta - \cancelto{0}{\cos[(n-\tfrac{1}{2})\pi] }~~\int_0^\infty (\cosh\eta + \cosh \theta)^{- \frac{5}{2}} e^{-n\theta} ~d\theta \biggr\} \, , </math>

where the prefactor of the second term — that is, <math>~\cos[(n-\tfrac{1}{2})\pi] </math> — goes to zero for all allowable values of the integer, <math>~n</math>. Hence, we conclude that,

<math>~2\cos(n\theta^') \int_0^\pi \frac{ \cos(n\theta)~d\theta }{ (\cosh\eta - \cos \theta)^{\frac{5}{2}} } </math>

<math>~=</math>

<math>~\frac{ 2(2\pi)^{\frac{1}{2}} Q_{n - \frac{1}{2}}^2 (\cosh\eta) \cos(n\theta^')}{ (\cosh^2\eta-1) ~\Gamma(\tfrac{5}{2})~ } </math>

 

<math>~=</math>

<math>~\biggl[ \frac{ 2^3 \sqrt{2} }{ 3 } \biggr] \frac{ Q_{n - \frac{1}{2}}^2(\cosh\eta) \cos(n\theta^')}{ \sinh^2\eta } \, . </math>

where we have set,

<math>~ \Gamma(\tfrac{5}{2}) = \Gamma(\tfrac{1}{2} + 2) = \frac{ \sqrt{\pi} \cdot 4! }{4^2 \cdot 2!} = \frac{\sqrt{\pi} \cdot 2^3\cdot 3}{ 2^5 } = \frac{3 \sqrt{\pi}}{2^2} \, . </math>

A.4

Beginning with the identified "Key Equation",

LSU Key.png

<math>~ \int_a^b\biggl[(\nu - \sigma)(\nu + \sigma + 1) + (\rho^2 - \mu^2)(1 - z^2)^{-1} \biggr] w_\nu^\mu ~w_\sigma^\rho ~dz </math>

<math>~=</math>

<math>~ \biggl[ z(\nu-\sigma) w_\nu^\mu ~w_\sigma^\rho + (\sigma+\rho) w_\nu^\mu ~ w_{\sigma-1}^\rho - (\nu + \mu) w_{\nu - 1}^\mu ~w_\sigma^\rho \biggr]_a^b </math>

A. Erdélyi (1953):  Volume I, §3.12, p. 169, eq. (1)

where, <math>~w_\nu^\mu(z)</math> and <math>~w_\sigma^\rho(z)</math> denote any solutions of Legendre's differential equation

we will adopt the associations:   <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math>, <math>~\nu \rightarrow (n - \tfrac{1}{2})</math>, <math>~\rho \rightarrow 0</math>, and <math>~\sigma \rightarrow ( n - \tfrac{1}{2})</math>. As a result, Erdélyi's (1953) expression becomes,

<math>~ \int_a^b\biggl[ -4(1 - t^2)^{-1} \biggr] Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) ~dt </math>

<math>~=</math>

<math>~ \biggl[ (n - \tfrac{1}{2} ) Q_{n - \frac{1}{2}}^2(t) ~ X_{n - \frac{3}{2}}(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) \biggr]_a^b \, . </math>

Drawing upon the recurrence "Key Equation,"

LSU Key.png

<math>~(\nu - \mu + 1)P^\mu_{\nu + 1} (z)</math>

<math>~=</math>

<math>~ (2\nu + 1)z P_\nu^\mu(z) - (\nu + \mu)P^\mu_{\nu-1}(z) </math>

Abramowitz & Stegun (1995), p. 334, eq. (8.5.3)

NOTE: <math>~Q_\nu^\mu</math>, as well as <math>~P_\nu^\mu</math>, satisfies this same recurrence relation.

which means, after making the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 0</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, that,

<math>~(n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>

<math>~=</math>

<math>~2n t X_{n-\frac{1}{2}}(t) - (n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t)</math>

<math>~\Rightarrow ~~~(n-\tfrac{1}{2})X_{n - \frac{3}{2}}(t) </math>

<math>~=</math>

<math>~2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t)</math>

the integral can be rewritten as,

<math>~ \int_a^b\biggl[ \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt </math>

<math>~=</math>

<math>~ \frac{1}{4} \biggl\{ \biggl[ 2n t X_{n-\frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) \biggr] Q_{n - \frac{1}{2}}^2(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) \biggr\}_a^b </math>

 

<math>~=</math>

<math>~ \frac{1}{4} \biggl\{ \biggl[ 2n t Q_{n - \frac{1}{2}}^2(t) - (n + \tfrac{3}{2}) Q_{n - \frac{3}{2}}^2(t) \biggr] X_{n - \frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t) \biggr\}_a^b \, . </math>

Returning to the same recurrence "Key Equation," but this time adopting the associations, <math>~z \rightarrow t</math>, <math>~\mu \rightarrow 2</math> and <math>~\nu \rightarrow (n-\tfrac{1}{2})</math>, we can write,

<math>~(n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t)</math>

<math>~=</math>

<math>~ 2n t Q^2_{n - \frac{1}{2}}(t) - (n + \tfrac{3}{2}) Q^2_{n - \frac{3}{2}} (t) \, , </math>

in which case the integral becomes,

<math>~ \int_a^b\biggl[ \frac{Q_{n - \frac{1}{2}}^2(t) ~X_{n - \frac{1}{2}}(t) }{(t^2-1)}\biggr]~dt </math>

<math>~=</math>

<math>~ \frac{1}{4} \biggl\{ (n - \tfrac{3}{2})Q^2_{n + \frac{1}{2}}(t) ~ X_{n - \frac{1}{2}}(t) - (n+\tfrac{1}{2})X_{n+\frac{1}{2}}(t) Q_{n - \frac{1}{2}}^2(t) \biggr\}_a^b \, . </math>

Appendix B: Elliptic Integrals

Summary Table

If the default argument of a complete elliptic integral function is the parameter, <math>~k</math>, then we should also appreciate that the following additional arguments may be useful:

parameter In terms of …
<math>~k</math>

(Descending Landen Transformation)

<math>~k_1</math>

(Ascending Landen Transformation)

<math>~\frac{R}{R_1} = e^{-\eta}</math> <math>~\eta \equiv \ln\biggl(\frac{R_1}{R}\biggr)</math>
<math>~k</math> <math>~k</math> <math>~\frac{2\sqrt{k_1}}{1+k_1}</math> <math>~\biggl[1 - \biggl( \frac{R}{R_1}\biggr)^2 \biggr]^{1 / 2}</math> <math>~\biggl[ \frac{2}{\coth\eta + 1} \biggr]^{1 / 2}= \biggl[ 1 - e^{-2\eta} \biggr]^{1 / 2}</math>
<math>~k^{\prime}</math> <math>~\sqrt{1 - k^2}</math> <math>~\frac{1-k_1}{1+k_1}</math> <math>~\frac{R}{R_1}</math> <math>~\frac{1}{\cosh\eta + \sinh\eta} = e^{-\eta}</math>
<math>~k_1 = \frac{1-k^{\prime}}{1+k^{\prime}}</math> <math>~\frac{1-\sqrt{1 - k^2}}{1+\sqrt{1 - k^2}}</math> <math>~k_1</math> <math>~\frac{R_1 - R}{R_1 +R}</math> <math>~\tanh\biggl(\frac{\eta}{2}\biggr) = \frac{1 - e^{-\eta}}{1 + e^{-\eta}}</math>


As published, for example, in §8.126 of Gradshteyn & Ryzhik (1965), the following mappings are valid:

Explicitly Published in §8.126 of G & R (1965)

Additional Implications

<math>~ K\biggl( \frac{1-k^\prime}{1+k^\prime} \biggr) = K(k_1)</math>

<math>~=</math>

<math>~\frac{1}{2} (1+k^{\prime}) K(k)</math>

<math>~ E\biggl( \frac{1-k^\prime}{1+k^\prime} \biggr) = E(k_1)</math>

<math>~=</math>

<math>~\frac{1}{1+k^\prime} \biggl[E(k) + k^\prime K(k) \biggr]</math>

<math>~ \frac{4K(k_1)}{R_1 + R} </math>

<math>~=</math>

<math>~\frac{2K(k)}{R_1}</math>

<math>~E(k)</math>

<math>~=</math>

<math>~(1+k^\prime)E(k_1) - k^\prime K(k)</math>

DLMS §19.8.12

<math>~ K\biggl( \frac{2\sqrt{k}}{1+k} \biggr) </math>

<math>~=</math>

<math>~(1+k)K(k)</math>

<math>~ E\biggl( \frac{2\sqrt{k}}{1+k} \biggr) </math>

<math>~=</math>

<math>~\frac{1}{1+k} \biggl[2E(k) - (k^\prime)^2 K(k) \biggr]</math>

<math>~ K( k ) </math>

<math>~=</math>

<math>~(1+k_1)K(k_1)</math>

DLMS §19.8.12

<math>~ E\biggl( \frac{2\sqrt{k_1}}{1+k_1} \biggr) = E(k) </math>

<math>~=</math>

<math>~\frac{1}{1+k_1} \biggl[2E(k_1) - (1-k_1^2) K(k_1) \biggr]</math>

One Set of Details

Following advice that we have received from Howard Cohl (private communication), we can demonstrate this most efficiently by employing the Descending Landen Transformation for the complete elliptic integral of the first kind, that is,

<math>~K(k)</math>

<math>~=</math>

<math>~ (1 + k_1)K(k_1) \, , </math>

      where,      

<math>~k_1</math>

<math>~\equiv</math>

<math>~ \frac{1-\sqrt{1-k^2}}{1+\sqrt{1-k^2}} \, . </math>

Given our just-stated definition of the parameter, <math>~k</math>, we have,

<math>~\sqrt{1-k^2}</math>

<math>~=</math>

<math>~ \biggl[ 1 - \frac{2}{\coth\eta + 1} \biggr]^{1 / 2} = \biggl[ \frac{\coth\eta - 1 }{\coth\eta + 1} \biggr]^{1 / 2} </math>

<math>~\Rightarrow ~~~ 1 \pm \sqrt{1-k^2}</math>

<math>~=</math>

<math>~ 1 \pm \biggl[ \frac{\coth\eta - 1 }{\coth\eta + 1} \biggr]^{1 / 2} = \frac{ [\coth\eta + 1]^{1 / 2} \pm [\coth\eta - 1]^{1 / 2} }{ [\coth\eta + 1]^{1 / 2} } </math>

<math>~\Rightarrow ~~~k_1 \equiv \frac{1 - \sqrt{1-k^2}}{1 + \sqrt{1-k^2}}</math>

<math>~=</math>

<math>~ \frac{ [\coth\eta + 1]^{1 / 2} - [\coth\eta - 1]^{1 / 2} }{ [\coth\eta + 1]^{1 / 2} + [\coth\eta - 1]^{1 / 2} } </math>

<math>~\Rightarrow ~~~k_1^2 </math>

<math>~=</math>

<math>~ \frac{ [\coth\eta + 1] -2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} + [\coth\eta - 1] }{ [\coth\eta + 1] + 2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} + [\coth\eta - 1] } </math>

 

<math>~=</math>

<math>~ \frac{ 2\coth\eta -2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} }{ 2\coth\eta + 2[\coth\eta + 1]^{1 / 2}[\coth\eta - 1]^{1 / 2} } </math>

 

<math>~=</math>

<math>~ \frac{ \coth\eta - [\coth^2\eta - 1]^{1 / 2} }{ \coth\eta + [\coth^2\eta - 1]^{1 / 2} } </math>

 

<math>~=</math>

<math>~ \frac{ \coth\eta/[\coth^2\eta - 1]^{1 / 2} - 1 }{ \coth\eta/[\coth^2\eta - 1]^{1 / 2} + 1 } </math>

 

<math>~=</math>

<math>~\frac{\cosh\eta - 1}{\cosh\eta + 1}</math>

 

<math>~=</math>

<math>~\tanh^2\biggl(\frac{\eta}{2}\biggr) \, .</math>

Note as well that,

<math>~1 + k_1</math>

<math>~=</math>

<math>~1+ \frac{1 - \sqrt{1-k^2}}{1 + \sqrt{1-k^2}} = \frac{2}{1 + \sqrt{1-k^2}} \, .</math>

Hence, the "second attempt" expression for <math>~P_{1 / 2}(z)</math> becomes,

<math>~P_{-1 / 2}(\cosh\eta)</math>

<math>~=</math>

<math>~ \frac{\sqrt{2}}{\pi}~ (\sinh\eta)^{-1 / 2} \biggl[ \frac{2k}{1 + \sqrt{1-k^2}} \biggr] K\biggl(\tanh \frac{\eta}{2} \biggr) </math>

 

<math>~=</math>

<math>~ \frac{2^{2}}{\pi}~ (\sinh\eta)^{-1 / 2} \biggl[ \frac{1}{\coth\eta + 1}\biggr]^{1 / 2} \biggl\{ \frac{ [\coth\eta + 1]^{1 / 2} }{ [\coth\eta + 1]^{1 / 2} + [\coth\eta - 1]^{1 / 2} } \biggr\}K\biggl(\tanh \frac{\eta}{2} \biggr) </math>

 

<math>~=</math>

<math>~ \frac{2^{2}}{\pi}~ \biggl[ ( \cosh\eta + \sinh\eta )^{1 / 2} + (\cosh\eta - \sinh\eta)^{1 / 2} \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr) </math>

 

<math>~=</math>

<math>~ \frac{2^{2}}{\pi}~ \biggl[ e^{\eta/2} + e^{-\eta/2} \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr) </math>

 

<math>~=</math>

<math>~ \biggl[ \frac{\pi}{2} \cdot \cosh\frac{\eta}{2} \biggr]^{-1} K\biggl(\tanh \frac{\eta}{2} \biggr) \, . </math>

This is, indeed, identical to the "Attempt #1" expression for <math>~P_{-1 / 2}(z)</math> that was used by Wong.

Q.E.D.

Series Expansions

From our accompanying set of Key mathematical relations we find the following series expansions for complete elliptic integrals of, respectively, the first and second kind:

LSU Key.png

<math>~\frac{2K(k)}{\pi}</math>

<math>~=</math>

<math>~ 1 + \biggl( \frac{1}{2} \biggr)^2k^2 + \biggl( \frac{1\cdot 3}{2\cdot 4}\biggr)^2 k^4 + \biggl( \frac{1\cdot 3\cdot 5}{2^4\cdot 3}\biggr)^2 k^6 + \biggl( \frac{1\cdot 3\cdot 5 \cdot 7}{2^7 \cdot 3}\biggr)^2 k^8 + \cdots + \biggl[ \frac{(2n-1)!!}{2^n n!} \biggr]^2 k^{2n} + \cdots </math>

Gradshteyn & Ryzhik (1965), §8.113.1


LSU Key.png

<math>~\frac{2E(k)}{\pi}</math>

<math>~=</math>

<math>~ 1 - \frac{1}{2^2} ~k^2 - \frac{1^2\cdot 3}{2^2\cdot 4^2}~ k^4 - \biggl(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}\biggr)^2~\frac{ k^6 }{5} - \biggl( \frac{1\cdot 3\cdot 5 \cdot 7}{2^7 \cdot 3}\biggr)^2 \frac{k^8}{7} ~-~ \cdots -

\biggl[ \frac{(2n-1)!!}{2^n n!} \biggr]^2 \frac{k^{2n}}{2n-1}

~-~ \cdots </math>

Gradshteyn & Ryzhik (1965), §8.114.1

If you don't know how to handle a double factorial — also referred to as a semifactorial — see, for example, this wikipedia page.

See Also

  • P. Lucht (2016), The Charged Bowl in Toroidal Coordinates — Note from Tohline: On 3 July 2018 I stumbled on this article by Phil Lucht, who lists his affiliation as Rimrock Digital Technology (rimrock at emission.com), Salt Lake City, Utah. As I have done over the past approximately half a year, Lucht appears to have spent quite a bit of time investigating how toroidal coordinates might be used to solve a few specific potential problems. He draws on many of the same, rich technical resource publications as I have done and, at one point, explicitly expresses amazement at the extensive amount of useful material that can be found in [MF53]. Anyone digging into this topic would be well-advised to read Lucht's Overview and Summary — especially including the subsection titled, Legendre Functions — and to recognize that his document contains a large number of technically detailed appendices.
     
    To provide one point of comparison with [MF53], Lucht derived an analytic expression for the potential outside of a torus (with toroidal-coordinate radius, <math>~\eta_0</math>) whose surface is held at uniform potential, <math>~\Phi_0</math>. This is in contrast to Wong's (1973) derivation — which we have verified, above — of the potential outside (as well as inside) of a torus having uniform charge/mass. For his particular example problem, Lucht obtains,

    <math>~\Phi(\eta,\theta)</math>

    <math>~=</math>

    <math>~ \Phi_0 \frac{\sqrt{2}}{\pi} \sqrt{\cosh\eta - \cos\theta} \sum_{n=0}^\infty \epsilon_n P_{n-\frac{1}{2}}(\cosh\eta) \Biggl[ \frac{Q_{n-\frac{1}{2}}(\cosh\eta_0)}{P_{n-\frac{1}{2}}(\cosh\eta_0)} \Biggr] \cos(n\theta) \, , </math>

    Lucht (2016), p. 68, Eq. (10.1.11)

    then he points out that the associated expression for the potential, which appears as equation (10.3.80) on p. 1304 of [MF53], is missing the Neumann factor, <math>~\epsilon_n</math>.


Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation