Appendix/Ramblings/Bordeaux

Université de Bordeaux (Part 1)

Spheroid-Ring Systems

Through a research collaboration at the Université de Bordeaux, B. Basillais & J. -M. Huré (2019), MNRAS, 487, 4504-4509 have published a paper titled, Rigidly Rotating, Incompressible Spheroid-Ring Systems: New Bifurcations, Critical Rotations, and Degenerate States.

We discuss this topic in a separate, accompanying chapter.

Exterior Gravitational Potential of Toroids

J. -M. Huré, B. Basillais, V. Karas, A. Trova, & O. Semerák (2020), MNRAS, 494, 5825-5838 have published a paper titled, The Exterior Gravitational Potential of Toroids. Here we examine how their work relates to the published work by C.-Y. Wong (1973, Annals of Physics, 77, 279), which we have separately discussed in detail.

Our Presentation of Wong's (1973) Result

Summary:  First three terms in Wong's (1973) expression for the gravitational potential at any point, P(ϖ, z), outside of a uniform-density torus.

Wong diagram $\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$ $=$ $-\left(\frac{2^3}{3{\pi }^3}\right){\mathrm{{\rm Y}}}_{W0}({\eta }_0)\{\frac{a}{r_1}\cdot \mathit{K}(k)\},$ $\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$ $=$ $-\left(\frac{2^3}{3{\pi }^3}\right){\mathrm{{\rm Y}}}_{W1}({\eta }_0)\times \cos \theta \{\frac{a}{r_2}\cdot \mathit{E}(k)\},$ $\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{2}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$ $=$ $-\left(\frac{2^3}{3{\pi }^3}\right){\mathrm{{\rm Y}}}_{W2}({\eta }_0)\times \cos (2\theta )\left\{\right[\frac{r_1^2+r_2^2}{2r_1{r}_2}]\frac{a}{r_2}\cdot \mathit{E}(k)-\frac{a}{r_1}\cdot \mathit{K}(k)\},$ where, once the major ( R ) and minor ( d ) radii of the torus — as well as the vertical location of its equatorial plane (Z0) — have been specified, we have, $a^2$ $\equiv$ $R^2-d^2$       and,       $\cosh {\eta }_0\equiv \frac{R}{d}\Rightarrow \sinh {\eta }_0=\frac{a}{d},$ $r_1^2$ $\equiv$ $(\varpi +a{)}^2+(z-Z_0{)}^2,$ $r_2^2$ $\equiv$ $(\varpi -a{)}^2+(z-Z_0{)}^2,$ $\cos \theta$ $\equiv$ $\left[\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2}\right],$ $k$ $\equiv$ $[\frac{r_1^2-r_2^2}{r_1^2}{]}^{1/2}=[\frac{4a\varpi }{r_1^2}{]}^{1/2}=[\frac{4a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2}{]}^{1/2}.$   Leading Coefficient Expressions … … evaluated for:    $\frac{R}{d}=\cosh {\eta }_0=3$ ${\mathrm{{\rm Y}}}_{W0}({\eta }_0)$ $\equiv$ $\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[{\cosh }^{}{\eta }_0+1]-E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\},$ 7.134677 ${\mathrm{{\rm Y}}}_{W1}({\eta }_0)$ $\equiv$ $\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[(3{\cosh }^{}{\eta }_0-1)]-5E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\},$ 0.130324 ${\mathrm{{\rm Y}}}_{W2}({\eta }_0)$ $\equiv$ $\frac{2^{3/2}}{3^2}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\{K(k_0)\cdot K(k_0)\cdot \cosh {\eta }_0(1-\cosh {\eta }_0)[16{\cosh }^{}{\eta }_0-13]+2K(k_0)\cdot E(k_0)[16{\cosh }^{}{\eta }_0-13{\cosh }^{}{\eta }_0+3]$       $-E(k_0)\cdot E(k_0)\cdot \cosh {\eta }_0(1+\cosh {\eta }_0)[3+16{\cosh }^{}{\eta }_0]\},$ 0.003153 where, $k_0$ $\equiv$ $[\frac{2}{\cosh {\eta }_0+1}{]}^{1/2}.$ 0.707106781 NOTE: In evaluating these "leading coefficient expressions" for the case, $R/d=3$, we have used the complete elliptic integral evaluations, K(k0) = 1.854074677 and E(k0) = 1.350643881.

Setup

From our accompanying discussion of Wong's (1973) derivation, the exterior potential is given by the expression,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}}(\eta ,\theta )$

$=$

$-D_0(\cosh \eta -\cos \theta {)}^{1/2}\sum _{n=0}^{\mathrm{n}\mathrm{m}\mathrm{a}\mathrm{x}}{ϵ}_n\cos (n\theta )C_n(\cosh {\eta }_0)P_{n-\frac{1}{2}}(\cosh \eta ),$

Wong (1973), §II.D, p. 294, Eqs. (2.59) & (2.61)

where,

$D_0$	$\equiv$	$\frac{2^{3/2}}{3{\pi }^2}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]=\frac{2^{3/2}}{3{\pi }^2}\left[\frac{(R^2-d^2{)}^{3/2}}{d^{2}R}\right],$
$C_n(\cosh {\eta }_0)$	$\equiv$	$(n+\frac{1}{2})Q_{n+\frac{1}{2}}(\cosh {\eta }_0)Q_{n-\frac{1}{2}}^2(\cosh {\eta }_0)-(n-\frac{3}{2})Q_{n-\frac{1}{2}}(\cosh {\eta }_0)Q_{n+\frac{1}{2}}^2(\cosh {\eta }_0)$

Wong (1973), §II.D, p. 294, Eq. (2.63)

and where, in terms of the major ( R ) and minor ( d ) radii of the torus — or their ratio, ε ≡ d/R,

$\cosh {\eta }_0$	$=$	$\frac{R}{d}=\frac{1}{ϵ},$
$\sinh {\eta }_0$	$=$	$\frac{a}{d}=\frac{1}{d}[R^2-d^2{]}^{1/2}=\frac{1}{ϵ}[1-{ϵ}^2{]}^{1/2}.$

These expressions incorporate a number of basic elements of a toroidal coordinate system. In what follows, we will also make use of the following relations:

Once the primary scale factor, $a$, has been specified, the illustration shown at the bottom of this inset box — see also our accompanying set of similar figures used by other researchers — helps in explaining how transformations can be made between any two of the referenced coordinate pairs: $(\varpi ,z)$, $(\eta ,\theta )$, $(r_1,r_2)$. $\varpi$ $=$ $\frac{a\sinh \eta }{(\cosh \eta -\cos \theta )}$      $\Rightarrow$      $\cos \theta$ $=$ $\cosh \eta -\frac{a\sinh \eta }{\varpi }$ $z-Z_0$ $=$ $\frac{a\sin \theta }{(\cosh \eta -\cos \theta )}$      $\Rightarrow$      $\sin \theta$ $=$ $\frac{(z-Z_0)}{\varpi }\cdot \sinh \eta$ Given that (sin2θ + cos2θ) = 1, we have, $1$ $=$ $[\frac{(z-Z_0)}{\varpi }\cdot \sinh \eta {]}^2+[\cosh \eta -\frac{a\sinh \eta }{\varpi }{]}^2$ $\Rightarrow \coth \eta$ $=$ $\frac{1}{2a\varpi }[{\varpi }^2+a^2+(z-Z_0{)}^2].$ We deduce as well that, $\frac{2}{\coth \eta +1}$ $=$ $\frac{4a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2},$        and, $\sinh \eta +\cosh \eta$ $=$ $\frac{{\varpi }^2+a^2+(z-Z_0{)}^2}{(\varpi +a{)}^2+(z-Z_0{)}^2}.$ Given the definitions, $r_1^2$ $=$ $(\varpi +a{)}^2+(z-Z_0{)}^2,$ $r_2^2$ $=$ $(\varpi -a{)}^2+(z-Z_0{)}^2,$ we can use the transformations, $\varpi$ $=$ $\frac{(r_1^2-r_2^2)}{4a}$     and, $(z-Z_0{)}^2$ $=$ $r_2^2-\frac{1}{16a^2}[r_1^2-r_2^2-4a^2{]}^2,$     or, $(z-Z_0{)}^2$ $=$ $r_1^2-\frac{1}{16a^2}[r_1^2-r_2^2+4a^2{]}^2.$ Or we can use the transformations, $\eta$ $=$ $\ln \left(\frac{r_1}{r_2}\right),$ $\cos \theta$ $=$ $\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2}.$ Additional potentially useful relations can be found in an accompanying chapter wherein we present a variety of basic elements of a toroidal coordinate system. Wong diagram

Leading (n = 0) Term

Wong's Expression

Now, from our separate derivation we have,

$P_{-1/2}(\cosh \eta )$

$=$

$\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{-1/2}{Q}_{-1/2}(\coth \eta ).$

And if we make the function-argument substitution, $z\to \coth \eta$, in the "Key Equation,"

File:LSU Key.png	$Q_{-\frac{1}{2}}(z)$	$=$	$\sqrt{\frac{2}{z+1}}K\left(\sqrt{\frac{2}{z+1}}\right)$
Abramowitz & Stegun (1995), p. 337, eq. (8.13.3)

we can write,

$P_{-1/2}(\cosh \eta )$

$=$

$\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{-1/2}k\mathit{K}(k),$

where, from above, we recognize that,

$k\equiv [\frac{2}{\coth \eta +1}{]}^{1/2}=[\frac{4a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2}{]}^{1/2}.$

So, the leading (n = 0) term gives,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}(\eta ,\theta )$	$=$	$-D_0(\cosh \eta -\cos \theta {)}^{1/2}{C}_0(\cosh {\eta }_0)P_{-\frac{1}{2}}(\cosh \eta )$
	$=$	$-D_0{C}_0(\cosh {\eta }_0)[\frac{a\sinh \eta }{\varpi }{]}^{1/2}\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{-1/2}k\mathit{K}(k)$
	$=$	$-\frac{D_0{C}_0(\cosh {\eta }_0)}{\pi }[\frac{2a}{\varpi }{]}^{1/2}k\mathit{K}(k)$
	$=$	$-C_0(\cosh {\eta }_0)\cdot \frac{2^3}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]\frac{a}{[(\varpi +a{)}^2+(z-Z_0{)}^2{]}^{1/2}}\cdot \mathit{K}(k).$

Thin-Ring Evaluation of C₀

In an accompanying discussion of the thin-ring approximation, we showed that as $\cosh {\eta }_0\to \mathrm{\infty }$

$C_0(x){|}_{x\to \mathrm{\infty }}$

$=$

$\left(\frac{3{\pi }^2}{2^2}\right)\frac{1}{{\cosh }^2{\eta }_0}.$

Hence, in this limit we can write,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}(\eta ,\theta ){|}_{\mathrm{t}\mathrm{h}\mathrm{i}\mathrm{n}\mathrm{-}\mathrm{r}\mathrm{i}\mathrm{n}\mathrm{g}}$

$=$

$-\frac{2}{\pi }{\xcancel{[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}{]}^3}}^1\frac{a}{[(\varpi +a{)}^2+(z-Z_0{)}^2{]}^{1/2}}\cdot \mathit{K}(k).$

More General Evaluation of C₀

NOTE of CAUTION: In our above evaluation of the toroidal function, $Q_{-\frac{1}{2}}(z)$, we appropriately associated the function argument, $z$, with the hyperbolic-cotangent of $\eta$; that is, we made the substitution, $z\to \coth \eta$. Here, as we assess the behavior of, and evaluate, the leading coefficient, $C_0$, an alternate substitution is appropriate, namely, $z_0\to \cosh {\eta }_0$; we affix the subscript zero to this function argument in an effort to minimize possible confusion with the argument, $z$.

Drawing from our accompanying tabulation of Toroidal Function Evaluations, we have more generally,

$2C_0(\cosh {\eta }_0)$	$=$	$[Q_{+\frac{1}{2}}(\cosh {\eta }_0)][Q_{-\frac{1}{2}}^2(\cosh {\eta }_0)]+3[Q_{-\frac{1}{2}}(\cosh {\eta }_0)][Q_{+\frac{1}{2}}^2(\cosh {\eta }_0)]$
	$=$	$[\cosh {\eta }_0{k}_{0}K(k_0)-[2(\cosh {\eta }_0+1){]}^{1/2}E(k_0)]\times \{\frac{4\cosh {\eta }_{0}E(k_0)-(\cosh {\eta }_0-1)K(k_0)}{[2^3(\cosh {\eta }_0+1)(\cosh {\eta }_0-1{)}^2{]}^{1/2}}\}$
		$-\frac{3}{2^2}[k_{0}K(k_0)]\times \{\cosh {\eta }_0{k}_{0}K(k_0)-({\cosh }^{}{\eta }_0+3)[\frac{2}{(\cosh {\eta }_0-1)({\cosh }^{}{\eta }_0-1)}{]}^{1/2}E(k_0)\},$

where,

$k_0$

$\equiv$

$[\frac{2}{\cosh {\eta }_0+1}{]}^{1/2}\Rightarrow (\cosh {\eta }_0+1)=\frac{2}{k_0^2}.$

Looking back at our previous numerical evaluation of $C_0(\cosh {\eta }_0)$ when $z_0=\cosh {\eta }_0=3\Rightarrow k_0=2^{-1/2}$, we see that, Appendix Expression: $Q_{-\frac{1}{2}}(z_0)$ $=$ $k_{0}K(k_0)$ Hence MF53 value, $Q_{-\frac{1}{2}}(3)$ $=$ $1.311028777\Rightarrow K(k_0)=1.854074677$ Appendix Expression: $Q_{+\frac{1}{2}}(z_0)$ $=$ $z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)$ Hence MF53 value, $Q_{+\frac{1}{2}}(3)$ $=$ $0.1128885424\Rightarrow E(k_0)=1.350643881$ Appendix Expression: $Q_{-\frac{1}{2}}^2(z_0)$ $=$ $[2^3(z-1)(z^2-1){]}^{-1/2}[4zE(k_0)-(z-1)K(k_0)]$ Hence, $Q_{-\frac{1}{2}}^2(3)$ $=$ $1.104816977$, which matches MF53 value Additional derivation: $Q_{+\frac{1}{2}}^2(z_0)$ $=$ $-\frac{1}{2^2}\{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$ Hence, $Q_{+\frac{1}{2}}^2(3)$ $=$ $0.449302588$ $\Rightarrow C_0(3)$ $=$ $\frac{1}{2}{Q}_{+\frac{1}{2}}(3)\cdot Q_{-\frac{1}{2}}^2(3)+\frac{3}{2}{Q}_{-\frac{1}{2}}(3)\cdot Q_{+\frac{1}{2}}^2(3)=0.945933522.$

Attempting to simplify this expression, we have,

$2C_0(\cosh {\eta }_0)$	$=$	$\{\cosh {\eta }_0{k}_{0}K(k_0)-(\frac{2}{k_0})E(k_0)\}\times \left\{\frac{4\cosh {\eta }_{0}E(k_0)-(\cosh {\eta }_0-1)K(k_0)}{[2^2{k}_0^{-1}(\cosh {\eta }_0-1)]}\right\}$
		$-\frac{3}{2^2}[k_{0}K(k_0)]\times \{\cosh {\eta }_0{k}_{0}K(k_0)-({\cosh }^{}{\eta }_0+3)[\frac{k_0^2}{(\cosh {\eta }_0-1{)}^2}{]}^{1/2}E(k_0)\}$
$\Rightarrow 2^3(\cosh {\eta }_0-1)C_0(\cosh {\eta }_0)$	$=$	$\{\cosh {\eta }_0{k}_0^{2}K(k_0)-2E(k_0)\}\times \{4\cosh {\eta }_{0}E(k_0)-(\cosh {\eta }_0-1)K(k_0)\}$
		$-3k_{0}K(k_0)\times \{\cosh {\eta }_0(\cosh {\eta }_0-1)k_{0}K(k_0)-({\cosh }^{}{\eta }_0+3)k_{0}E(k_0)\}$
	$=$	$-K(k_0)\cdot K(k_0)[(\cosh {\eta }_0-1)\cdot \cosh {\eta }_0{k}_0^2+3\cosh {\eta }_0(\cosh {\eta }_0-1)k_0^2]$
		$+K(k_0)\cdot E(k_0)[2^2{\cosh }^2{\eta }_0{k}_0^2+2(\cosh {\eta }_0-1)+3k_0^2({\cosh }^{}{\eta }_0+3)]-E(k_0)\cdot E(k_0)[2^3\cosh {\eta }_0]$
$\Rightarrow \left[\frac{2^3(\cosh {\eta }_0-1)}{k_0^2}\right]C_0(\cosh {\eta }_0)$	$=$	$-K(k_0)\cdot K(k_0)[(\cosh {\eta }_0-1)\cdot \cosh {\eta }_0+3\cosh {\eta }_0(\cosh {\eta }_0-1)]$
		$+K(k_0)\cdot E(k_0)[2^2{\cosh }^2{\eta }_0+\frac{2}{k_0^2}(\cosh {\eta }_0-1)+3({\cosh }^{}{\eta }_0+3)]-E(k_0)\cdot E(k_0)\left[\frac{2^3\cosh {\eta }_0}{k_0^2}\right]$
$\Rightarrow ({\cosh }^2{\eta }_0-1)C_0(\cosh {\eta }_0)$	$=$	$K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[{\cosh }^{}{\eta }_0+1]-E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]$

This last, simplifed expression gives, as above, $C_0(3)=0.945933523$. TERRIFIC!

Finally then, for any choice of ${\eta }_0$,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}(\eta ,\theta ){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\frac{a}{[(\varpi +a{)}^2+(z-Z_0{)}^2{]}^{1/2}}\cdot \mathit{K}(k)$
		$\times \{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[{\cosh }^{}{\eta }_0+1]-E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\}.$

Second (n = 1) Term

The second (n = 1) term in Wong's (1973) expression for the exterior potential is,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\eta ,\theta )$

$=$

$-D_0(\cosh \eta -\cos \theta {)}^{1/2}\cdot 2\cos \theta \cdot C_1(\cosh {\eta }_0)P_{+\frac{1}{2}}(\cosh \eta ),$

where, $D_0$ is the same as above, and,

$C_1(\cosh {\eta }_0)$

$\equiv$

$\frac{3}{2}{Q}_{+\frac{3}{2}}(\cosh {\eta }_0)Q_{+\frac{1}{2}}^2(\cosh {\eta }_0)+\frac{1}{2}{Q}_{+\frac{1}{2}}(\cosh {\eta }_0)Q_{+\frac{3}{2}}^2(\cosh {\eta }_0).$

Now, from our accompanying table of "Toroidal Function Evaluations", it appears as though,

$P_{+\frac{1}{2}}(\cosh \eta )$

$=$

$\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{+1/2}{k}^{-1}E(k),$

where, as above,

$k$

$\equiv$

$[\frac{2}{\coth \eta +1}{]}^{1/2}.$

Hence, we have,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\eta ,\theta )$	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)[\cos \theta \cdot (\cosh \eta -\cos \theta {)}^{1/2}(\sinh \eta {)}^{+1/2}]k^{-1}E(k)$
	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\cdot \cos \theta \{\frac{a{\sinh }^2\eta }{\varpi }\cdot \frac{\coth \eta +1}{2}{\}}^{1/2}E(k)$
	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\cdot \cos \theta \left\{\right(\frac{a}{2\varpi }\left)\right[\frac{r_1^2-r_2^2}{2r_1{r}_2}{]}^2\cdot \left[\frac{2r_1^2}{r_1^2-r_2^2}\right]{\}}^{1/2}E(k)$
	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\cdot \cos \theta \left\{\right(\frac{a}{2}\left)\right[\frac{4a}{r_1^2-r_2^2}\left]\right[\frac{r_1^2-r_2^2}{2r_1{r}_2}{]}^2\cdot \left[\frac{2r_1^2}{r_1^2-r_2^2}\right]{\}}^{1/2}E(k)$
	$=$	$-\frac{2^{3}a}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\left[\frac{\cos \theta }{r_2}\right]E(k)=-\frac{2^{3}a}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\left[\frac{\cos \theta }{\sqrt{(\varpi -a{)}^2+(z-Z_0{)}^2}}\right]E(k)$
	$=$	$-\frac{2^{3}a}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\left[\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2^2}\right]\mathit{E}(k).$

 

From the above function tabulations & evaluations — for example, $K(k_0)=1.854074677$ and $E(k_0)=1.350643881$ — and a separate listing of Example Recurrence Relations, we have, Appendix Expression: $Q_{-\frac{1}{2}}(z_0)$ $=$ $k_{0}K(k_0)$ Appendix Expression: $Q_{+\frac{1}{2}}(z_0)$ $=$ $z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)$ $Q_{m-\frac{1}{2}}^0$ recurrence with m = 2: $Q_{+\frac{3}{2}}(z_0)$ $=$ $\frac{4}{3}z{Q}_{+\frac{1}{2}}(z_0)-\frac{1}{3}{Q}_{-\frac{1}{2}}(z_0)$   $=$ $\frac{4}{3}z\{z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)\}-\frac{1}{3}{k}_{0}K(k_0)$   $=$ $\frac{1}{3}[(4z^2-1)k_{0}K(k_0)-4z[2(z+1){]}^{1/2}E(k_0)]$ Hence, $Q_{+\frac{3}{2}}(3)$ $=$ $0.014544576.$ Appendix Expression: $Q_{-\frac{1}{2}}^2(z_0)$ $=$ $[2^3(z-1)(z^2-1){]}^{-1/2}[4zE(k_0)-(z-1)K(k_0)]$ Additional derivation: $Q_{+\frac{1}{2}}^2(z_0)$ $=$ $-\frac{1}{2^2}\{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$ Then, letting $\mu \to 2$ and, for all m ≥ 2, letting $\nu \to (m-\frac{1}{2})$ in the "Key Equation," $(\nu -\mu +1)P_{\nu +1}^{\mu }(z)$ $=$ $(2\nu +1)z{P}_{\nu }^{\mu }(z)-(\nu +\mu )P_{\nu -1}^{\mu }(z)$ Abramowitz & Stegun (1995), p. 334, eq. (8.5.3) NOTE: $Q_{\nu }^{\mu }$, as well as $P_{\nu }^{\mu }$, satisfies this same recurrence relation. we have, $(m-\frac{3}{2})Q_{m+\frac{1}{2}}^2(z)$ $=$ $(2m)z{Q}_{m-\frac{1}{2}}^2(z)-(m+\frac{3}{2})Q_{m-\frac{3}{2}}^2(z).$ Therefore, specifically for m = 1, we obtain the recurrence relation, $Q_{+\frac{3}{2}}^2(z_0)$ $=$ $5Q_{-\frac{1}{2}}^2(z_0)-4z{Q}_{+\frac{1}{2}}^2(z_0)$   $=$ $5\{[2^3(z-1)(z^2-1){]}^{-1/2}[4zE(k_0)-(z-1)K(k_0)]\}+z\{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$   $=$ $2^{1/2}[(z-1)(z^2-1){]}^{-1/2}\{[5z]-z(z^2+3)\}E(k_0)+\{z^2{k}_0-[(z-1)(z^2-1){]}^{-1/2}[2^{-3/2}\cdot 5(z-1)]\}K(k_0)$   $=$ $2^{-3/2}(z+1{)}^{-1/2}[4z^2-5]K(k_0)-2^{1/2}[(z-1)(z^2-1){]}^{-1/2}(z^2-2)zE(k_0)$ Hence, $Q_{+\frac{3}{2}}^2(3)$ $=$ $0.132453829.$ $\Rightarrow C_1(3)$ $=$ $\frac{3}{2}{Q}_{+\frac{3}{2}}(3)\cdot Q_{+\frac{1}{2}}^2(3)+\frac{1}{2}{Q}_{+\frac{1}{2}}(3)\cdot Q_{+\frac{3}{2}}^2(3)=0.017278633.$

While keeping in mind that,

$z_0$

$=$

$\cosh {\eta }_0,$

and,

$k_0^2$

$=$

$\frac{2}{\cosh {\eta }_0+1}=\frac{2}{z_0+1},$

let's attempt to express this leading coefficient, $C_1(\cosh {\eta }_0)$, entirely in terms of the pair of complete elliptic integral functions.

$2C_1(z_0)$	$=$	$3[Q_{+\frac{3}{2}}(z_0)]\times [Q_{+\frac{1}{2}}^2(z_0)]+[Q_{+\frac{1}{2}}(z_0)]\times [5Q_{-\frac{1}{2}}^2(z_0)-4z{Q}_{+\frac{1}{2}}^2(z_0)]$
	$=$	$[3Q_{+\frac{3}{2}}(z_0)-4z{Q}_{+\frac{1}{2}}(z_0)]\times [Q_{+\frac{1}{2}}^2(z_0)]+[5Q_{+\frac{1}{2}}(z_0)]\times [Q_{-\frac{1}{2}}^2(z_0)]$
	$=$	$-\frac{1}{2^2}\{(4z^2-1)k_{0}K(k_0)-4z[2(z+1){]}^{1/2}E(k_0)-4z[z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)\left]\right\}\times \{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$
		$+5[z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)]\times \{[2^3(z-1)(z^2-1){]}^{-1/2}[4zE(k_0)-(z-1)K(k_0)]\}$
	$=$	$\frac{1}{2^2}\cdot k_{0}K(k_0)\times \{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$
		$+5[2^3(z-1)(z^2-1){]}^{-1/2}[z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)]\times \{4zE(k_0)-(z-1)K(k_0)\}$
	$=$	$K(k_0)\cdot K(k_0)\{\frac{z{k}_0^2}{2^2}-5[2^3(z-1)(z^2-1){]}^{-1/2}\cdot z{k}_0(z-1)\}+E(k_0)\cdot E(k_0)\{-5[2^3(z-1)(z^2-1){]}^{-1/2}\cdot [2(z+1){]}^{1/2}\cdot 4z\}$
		$+K(k_0)\cdot E(k_0)\{-\frac{1}{2^2}\cdot k_0(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}+5[2^3(z-1)(z^2-1){]}^{-1/2}\cdot 4z^2{k}_0+5[2^3(z-1)(z^2-1){]}^{-1/2}\cdot [2(z+1){]}^{1/2}\cdot (z-1)\}.$

Hence,

$2[(z-1)(z^2-1){]}^{1/2}{C}_1(z_0)$	$=$	$z{k}_0\cdot K(k_0)\cdot K(k_0)\left\{\frac{k_0}{2^2}\right[(z-1)(z^2-1){]}^{1/2}-\frac{5(z-1)}{2^{3/2}}\}-10z(z+1{)}^{1/2}\cdot E(k_0)\cdot E(k_0)$
		$+2^{-3/2}K(k_0)\cdot E(k_0)\{k_0[19z^2-3]+5(z-1)[2(z+1){]}^{1/2}\}$
$\Rightarrow 2^{3/2}\left[\frac{(z-1)}{k_0}\right]C_1(z_0)$	$=$	$z{k}_0\cdot K(k_0)\cdot K(k_0)\left\{\frac{k_0}{2^2}\right[\frac{2^{1/2}(z-1)}{k_0}]-\frac{5(z-1)}{2^{3/2}}\}-\left[\frac{2^{3/2}\cdot 5z}{k_0}\right]E(k_0)\cdot E(k_0)$
		$+2^{-3/2}K(k_0)\cdot E(k_0)\{k_0[19z^2-3]+\frac{10(z-1)}{k_0}\}$
$\Rightarrow C_1(z_0)$	$=$	$\left[\frac{2(3z^2-1)}{(z^2-1)}\right]K(k_0)\cdot E(k_0)-\left[\frac{z}{(z+1)}\right]K(k_0)\cdot K(k_0)-\left[\frac{5z}{(z-1)}\right]E(k_0)\cdot E(k_0)$
$\Rightarrow (z_0^2-1)C_1(z_0)$	$=$	$2(3z^2-1)K(k_0)\cdot E(k_0)-z_0(z_0-1)K(k_0)\cdot K(k_0)-5z_0(z_0+1)E(k_0)\cdot E(k_0).$

Hence, we have,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\eta ,\theta )$

$=$

$-\frac{2^{3}a}{3{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_1(\cosh {\eta }_0)\left[\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2^2}\right]\mathit{E}(k).$

Third (n = 2) Term

Part A

The third (n = 2) term in Wong's (1973) expression for the exterior potential is,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{2}}(\eta ,\theta )$

$=$

$-D_0(\cosh \eta -\cos \theta {)}^{1/2}\cdot 2\cos (2\theta )\cdot C_2(\cosh {\eta }_0)P_{+\frac{3}{2}}(\cosh \eta ),$

where, $D_0$ is the same as above, and,

$C_2(\cosh {\eta }_0)$

$\equiv$

$\frac{5}{2}{Q}_{+\frac{5}{2}}(\cosh {\eta }_0)Q_{+\frac{3}{2}}^2(\cosh {\eta }_0)-\frac{1}{2}{Q}_{+\frac{3}{2}}(\cosh {\eta }_0)Q_{+\frac{5}{2}}^2(\cosh {\eta }_0).$

In order to evaluate $C_2(z)$, we will need the following pair of expressions in addition to the ones already used: $Q_{m-\frac{1}{2}}^0$ recurrence with m = 3, gives:     $Q_{+\frac{5}{2}}(z_0)$ $=$ $\frac{8}{5}z{Q}_{+\frac{3}{2}}(z_0)-\frac{3}{5}{Q}_{+\frac{1}{2}}(z_0)$ $\Rightarrow 15Q_{+\frac{5}{2}}(z_0)$ $=$ $8z[(4z^2-1)k_{0}K(k_0)-4z[2(z+1){]}^{1/2}E(k_0)]-9[z{k}_{0}K(k_0)-[2(z+1){]}^{1/2}E(k_0)]$   $=$ $z{k}_{0}K(k_0)[8(4z^2-1)-9]+[2(z+1){]}^{1/2}E(k_0)[-32z^2+9]$   $=$ $z{k}_{0}K(k_0)[32z^2-17]+[2(z+1){]}^{1/2}E(k_0)[9-32z^2].$ Hence,     $Q_{+\frac{5}{2}}(3)$ $=$ $0.002080867.$ And, setting m = 2 in the above recurrence relation for $Q_{m+\frac{1}{2}}^2(z)$ gives, $[(m-\frac{3}{2})Q_{m+\frac{1}{2}}^2(z){]}_{m=2}$ $=$ $[(2m)z{Q}_{m-\frac{1}{2}}^2(z)-(m+\frac{3}{2})Q_{m-\frac{3}{2}}^2(z){]}_{m=2}$ $\Rightarrow Q_{+\frac{5}{2}}^2(z)$ $=$ $8z{Q}_{+\frac{3}{2}}^2(z)-7Q_{+\frac{1}{2}}^2(z)$   $=$ $8z[5Q_{-\frac{1}{2}}^2(z_0)-4z{Q}_{+\frac{1}{2}}^2(z_0)]-7Q_{+\frac{1}{2}}^2(z)$   $=$ $40z{Q}_{-\frac{1}{2}}^2(z_0)-[32z^2+7]Q_{+\frac{1}{2}}^2(z_0)$   $=$ $40z\{[2^3(z-1)(z^2-1){]}^{-1/2}[4zE(k_0)-(z-1)K(k_0)]\}$     $+\frac{[32z^2+7]}{4}\{z{k}_{0}K(k_0)-(z^2+3)[\frac{2}{(z-1)(z^2-1)}{]}^{1/2}E(k_0)\}$ $\Rightarrow 4Q_{+\frac{5}{2}}^2(z)$ $=$ $2^5\cdot 5z\{2^{1/2}[(z-1)(z^2-1){]}^{-1/2}[zE(k_0)]-2^{-3/2}[(z-1)(z^2-1){]}^{-1/2}[(z-1)K(k_0)]\}$     $+[32z^2+7]\{z{k}_{0}K(k_0)-2^{1/2}(z^2+3)[(z-1)(z^2-1){]}^{-1/2}E(k_0)\}$   $=$ $\{2^{11/2}\cdot 5[z^2]-2^{1/2}[32z^2+7](z^2+3)\}[(z-1)(z^2-1){]}^{-1/2}E(k_0)$     $-2^{7/2}\cdot 5\{[(z-1)(z^2-1){]}^{-1/2}(z-1)\}zK(k_0)+[32z^2+7]\{2^{1/2}[z+1{]}^{-1/2}\}zK(k_0)$   $=$ $2^{1/2}\{32z^2-33\}z[z+1{]}^{-1/2}K(k_0)-2^{1/2}\{32z^4-57z^2+21\}[(z-1)(z^2-1){]}^{-1/2}E(k_0)$   $=$ $2^{1/2}[z+1{]}^{-1/2}[z-1{]}^{-1}[(32z^2-33)z(z-1)K(k_0)-(32z^4-57z^2+21)E(k_0)].$ Hence,     $Q_{+\frac{5}{2}}^2(3)$ $=$ $0.03377378.$

Part B

Let's evaluate $C_2(z)$ specifically for the case where $z=\cosh {\eta }_0=3$, using the already separately evaluated values of the four relevant toroidal functions. We find,

$2C_2(3)$	$=$	$5Q_{+\frac{5}{2}}(3)Q_{+\frac{3}{2}}^2(3)-Q_{+\frac{3}{2}}(3)Q_{+\frac{5}{2}}^2(3)$
	$=$	$5\cdot (0.002080867)\times (0.132453829)-(0.014544576)\times (0.03377378)$
	$=$	$8.868687\times 10^{-4}.$

Next, let's develop a consolidated expression for $C_2(z_0)$ that replaces all the toroidal functions with complete elliptic integrals of the first and second kind.

$2C_2(z_0)$	$=$	$5Q_{+\frac{5}{2}}(z_0)Q_{+\frac{3}{2}}^2(z_0)-Q_{+\frac{3}{2}}(z_0)Q_{+\frac{5}{2}}^2(z_0)$
	$=$	$\frac{1}{3}\{z{k}_{0}K(k_0)[32z^2-17]+[2(z+1){]}^{1/2}E(k_0)[9-32z^2]\}\times \{2^{-3/2}(z+1{)}^{-1/2}[4z^2-5]K(k_0)-2^{1/2}[(z-1)(z^2-1){]}^{-1/2}(z^2-2)zE(k_0)\}$
		$-\frac{1}{2^2\cdot 3}\{(4z^2-1)k_{0}K(k_0)-4z[2(z+1){]}^{1/2}E(k_0)\}\times \{2^{1/2}[z+1{]}^{-1/2}[z-1{]}^{-1}[(32z^2-33)z(z-1)K(k_0)-(32z^4-57z^2+21)E(k_0)]\}$

$\Rightarrow 2^2\cdot 3(z^2-1)C_2(z_0)$	$=$	$\{K(k_0)z[32z^2-17]+(z+1)E(k_0)[9-32z^2]\}\times \{(z-1)[4z^2-5]K(k_0)-4(z^2-2)zE(k_0)\}$
		$-\{(4z^2-1)K(k_0)-4z(z+1)E(k_0)\}\times \{(32z^2-33)z(z-1)K(k_0)-(32z^4-57z^2+21)E(k_0)\}$
	$=$	$\{(z-1)[32z^2-17][4z^2-5]zK(k_0)\cdot K(k_0)-4(z^2-2)z^2[32z^2-17]K(k_0)\cdot E(k_0)\}$
		$+\{(z-1)(z+1)[9-32z^2][4z^2-5]K(k_0)\cdot E(k_0)-4(z^2-2)z(z+1)[9-32z^2]E(k_0)\cdot E(k_0)\}$
		$+\{(32z^4-57z^2+21)(4z^2-1)K(k_0)\cdot E(k_0)-(32z^2-33)z(z-1)(4z^2-1)K(k_0)\cdot K(k_0)\}$
		$+\{4z(z+1)(32z^2-33)z(z-1)K(k_0)\cdot E(k_0)-4z(z+1)(32z^4-57z^2+21)E(k_0)\cdot E(k_0)\}$
	$=$	$(z-1)\left\{\right[(32z^2-17)(4z^2-5)z]-[(32z^2-33)z(4z^2-1)\left]\right\}K(k_0)\cdot K(k_0)$
		$+\left\{\right[(z-1)(z+1)(9-32z^2)(4z^2-5)]-[4(z^2-2)z^2(32z^2-17)]$
		$+[(32z^4-57z^2+21)(4z^2-1)]+[4z(z+1)(32z^2-33)z(z-1)]\}K(k_0)\cdot E(k_0)$
		$-2z(z+1)\left\{\right[2(32z^4-57z^2+21)]+2[(z^2-2)(9-32z^2)\left]\right\}E(k_0)\cdot E(k_0)$
	$=$	$z(z-1)\{[52-64z^2]\}K(k_0)\cdot K(k_0)-4z(z+1)\{[3+16z^2]\}E(k_0)\cdot E(k_0)$
		$+\left\{\right[5(32z^4-41z^2+9)]-[(32z^4-57z^2+21)]$
		$+4z^2[(32z^4-57z^2+21)+(32z^4-65z^2+33)+(-32z^4+41z^2-9)+(-32z^4+81z^2-34)]\}K(k_0)\cdot E(k_0)$
	$=$	$4z(z-1)\{13-16z^2\}K(k_0)\cdot K(k_0)-4z(z+1)\{[3+16z^2]\}E(k_0)\cdot E(k_0)+8\{16z^4-13z^2+3\}K(k_0)\cdot E(k_0).$

Finally, let's evaluate this consolidated expression for the specific case of $z_0=\cosh {\eta }_0=3$, remembering that in this specific case $k_0=2^{-1/2}$, $K(k_0)=1.854074677$, and $E(k_0)=1.350643881$. We find,

$2C_2(z_0)$	$=$	$[2\cdot 3(z^2-1){]}^{-1}\{4z(z-1)[13-16z^2]K(k_0)\cdot K(k_0)-4z(z+1)[3+16z^2]E(k_0)\cdot E(k_0)+8[16z^4-13z^2+3]K(k_0)\cdot E(k_0)\}$
	$=$	$[48{]}^{-1}\{-24[131]K(k_0)\cdot K(k_0)-48[147]E(k_0)\cdot E(k_0)+8[1182]K(k_0)\cdot E(k_0)\}$
	$=$	$8.8708\times 10^{-4}.$

This matches the numerically evaluated expression, from above (6/30/2020). There is a tremendous amount of cancellation between the three key terms in this expression, so the match is only to three significant digits.

Part C

Next …

Useful Relations from Above $\cosh \eta$ $=$ $\frac{r_1^2+r_2^2}{2r_1{r}_2};$ $\sinh \eta$ $=$ $\frac{r_1^2-r_2^2}{2r_1{r}_2};$ $\varpi$ $=$ $\frac{r_1^2-r_2^2}{2a};$ $\cosh \eta -\cos \theta$ $=$ $\frac{2a^2}{r_1{r}_2};$ $\cos \theta$ $=$ $\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2};$ $\frac{2}{\coth \eta +1}$ $=$ $\frac{4a\varpi }{r_1^2}.$

Now, from our tabulation of example recurrence relations, we see that,

$P_{+\frac{3}{2}}(\cosh \eta )$	$=$	$\frac{4}{3}\cdot \cosh \eta P_{+\frac{1}{2}}(\cosh \eta )-\frac{1}{3}{P}_{-\frac{1}{2}}(\cosh \eta )$
	$=$	$\frac{4}{3}\cdot \cosh \eta [\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{+1/2}{k}^{-1}\mathit{E}(k)]-\frac{1}{3}[\frac{\sqrt{2}}{\pi }(\sinh \eta {)}^{-1/2}k\mathit{K}(k)]$
	$=$	$\frac{2^{1/2}}{3\pi }[4\cosh \eta (\sinh \eta {)}^{+1/2}{k}^{-1}\mathit{E}(k)-(\sinh \eta {)}^{-1/2}k\mathit{K}(k)],$

where, as above,

$k\equiv [\frac{2}{\coth \eta +1}{]}^{1/2}=[\frac{4a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2}{]}^{1/2}=[\frac{4a\varpi }{r_1^2}{]}^{1/2}.$

So we have,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{2}}(\eta ,\theta )$	$=$	$-\frac{2^{5/2}}{3{\pi }^2}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_2(\cosh {\eta }_0)\cos (2\theta )\{(\cosh \eta -\cos \theta {)}^{1/2}{P}_{+\frac{3}{2}}(\cosh \eta )\}$
	$=$	$-\frac{2^3}{3^2{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_2(\cosh {\eta }_0)\cos (2\theta )\left(\frac{2a^2}{r_1{r}_2}{)}^{1/2}\right\{4\cosh \eta (\sinh \eta {)}^{+1/2}{k}^{-1}\mathit{E}(k)-(\sinh \eta {)}^{-1/2}k\mathit{K}(k)\}$
	$=$	$-\frac{2^3}{3^2{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_2(\cosh {\eta }_0)\cos (2\theta )(\frac{2a^2}{r_1{r}_2}{)}^{1/2}$
		$\times \left\{4\right[\frac{r_1^2+r_2^2}{2r_1{r}_2}\left]\right[\frac{r_1^2-r_2^2}{2r_1{r}_2}{]}^{+1/2}\left[\frac{4a}{r_1^2}\right(\frac{r_1^2-r_2^2}{2a}){]}^{-1/2}\mathit{E}(k)-[\frac{r_1^2-r_2^2}{2r_1{r}_2}{]}^{-1/2}\left[\frac{4a}{r_1^2}\right(\frac{r_1^2-r_2^2}{2a}){]}^{1/2}\mathit{K}(k)\}$
	$=$	$-\frac{2^{9/2}}{3^2{\pi }^3}\left[\frac{{\sinh }^3{\eta }_0}{\cosh {\eta }_0}\right]C_2(\cosh {\eta }_0)\cos (2\theta )\times \left\{\right[\frac{r_1^2+r_2^2}{2r_1{r}_2}]\frac{a}{r_2}\cdot \mathit{E}(k)-\frac{a}{r_1}\cdot \mathit{K}(k)\}.$

Finally, inserting the expression for ${\sinh }^2{\eta }_0{C}_2(\cosh {\eta }_0)=(z_0^2-1)C_2(z_0)$ that we have derived, above, gives,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{2}}(\eta ,\theta )$	$=$	$-\frac{2^{9/2}}{3^3{\pi }^3}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\times \cos (2\theta )\left\{\right[\frac{r_1^2+r_2^2}{2r_1{r}_2}]\frac{a}{r_2}\cdot \mathit{E}(k)-\frac{a}{r_1}\cdot \mathit{K}(k)\}$
		$\times \{z(z-1)[13-16z^2]K(k_0)\cdot K(k_0)-z(z+1)[3+16z^2]E(k_0)\cdot E(k_0)+2[16z^4-13z^2+3]K(k_0)\cdot E(k_0)\}.$

Summary

Once the major ( R ) and minor ( d ) radii of the torus have been specified, two key model parameters can be immediately determined, namely,

$a^2\equiv R^2-d^2,$       and,       $\cosh {\eta }_0\equiv \frac{R}{d},$

in which case also, $\sinh {\eta }_0=a/d.$ Once the mass-density ( ρ₀ ) of the torus has been specified, the torus mass is given by the expression,

$M=2{\pi }^2{\rho }_0{d}^{2}R.$

In addition to the principal pair of meridional-plane coordinates, $(\varpi ,z)$, it is useful to define the pair of distances,

$r_1^2$	$\equiv$	$(\varpi +a{)}^2+(z-Z_0{)}^2,$
$r_2^2$	$\equiv$	$(\varpi -a{)}^2+(z-Z_0{)}^2,$

where, the equatorial plane of the torus is located at $z=Z_0$. As we have shown above, the leading (n = 0) term in Wong's (1973) series-expression for the gravitational potential anywhere outside the torus is,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\frac{a}{r_1}\cdot \mathit{K}(k)$
		$\times \{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[{\cosh }^{}{\eta }_0+1]-E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\}.$

where, the two distinctly different arguments — one with, and one without a zero subscript — of the complete elliptic-integral functions are,

$k$	$\equiv$	$[\frac{2}{\coth \eta +1}{]}^{1/2}=[\frac{4a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2}{]}^{1/2}=[\frac{4a\varpi }{r_1^2}{]}^{1/2}=[\frac{r_1^2-r_2^2}{r_1^2}{]}^{1/2},$
$k_0$	$\equiv$	$[\frac{2}{\cosh {\eta }_0+1}{]}^{1/2}.$

As we also have shown above, the second (n = 1) term in Wong's (1973) series-expression for the exterior gravitational potential is,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\left[\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2}\right]\frac{a}{r_2}\cdot \mathit{E}(k)$
		$\times \{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[(3{\cosh }^{}{\eta }_0-1)]-5E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\}.$

Note that a transformation from the $(r_1,r_2)$ coordinate pair to the toroidal-coordinate pair $(\eta ,\theta )$ includes the expression,

$\cos \theta$

$=$

$\frac{r_1^2+r_2^2-4a^2}{2r_1{r}_2}.$

So this (n = 1) term's explicit dependence on "cos(nθ)" is clear. Finally, the third (n = 2) term in Wong's (1973) series-expression for the exterior gravitational potential is,

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{2}}(\eta ,\theta )$	$=$	$-\frac{2^3}{3{\pi }^3}\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\times \cos (2\theta )\left\{\right[\frac{r_1^2+r_2^2}{2r_1{r}_2}]\frac{a}{r_2}\cdot \mathit{E}(k)-\frac{a}{r_1}\cdot \mathit{K}(k)\}$
		$\times \frac{2^{3/2}}{3^2}\{K(k_0)\cdot K(k_0)\cdot \cosh {\eta }_0(1-\cosh {\eta }_0)[16{\cosh }^{}{\eta }_0-13]+2K(k_0)\cdot E(k_0)[16{\cosh }^{}{\eta }_0-13{\cosh }^{}{\eta }_0+3]$
		$-E(k_0)\cdot E(k_0)\cdot \cosh {\eta }_0(1+\cosh {\eta }_0)[3+16{\cosh }^{}{\eta }_0]\}.$

The Huré, et al (2020) Presentation

Notation

In Huré, et al. (2020), the major and minor radii of the torus surface ("shell") are labeled, respectively, R_c and b, and their ratio is denoted,

$e\equiv \frac{b}{R_c}.$

Huré, et al. (2020), §2, p. 5826, Eq. (1)

The authors work in cylindrical coordinates, $(R,Z)$, whereas we refer to this same coordinate-pair as, $({\varpi }_W,z_W)$. The quantity,

${\mathrm{\Delta }}^2$

$\equiv$

$[R+(R_c+b\cos {\theta }_H){]}^2+[Z-b\sin {\theta }_H{]}^2.$

Huré, et al. (2020), §2, p. 5826, Eqs. (5) & (7)

We have affixed the subscript "H" to their meridional-plane angle, θ, to clarify that it has a different coordinate-base definition from the meridional-plane angle, θ, that appears in our above discussion of Wong's (1973) work. In their paper, the subscript "0" is used in the case of an infinitesimally thin hoop $(b\to 0)$, that is to say,

${\mathrm{\Delta }}_0^2$

$=$

$[R+R_c{]}^2+Z^2.$

Huré, et al. (2020), §3, p. 5827, Eq. (13)

Generally, the argument (modulus) of the complete elliptic integral functions is,

$k_H$

$=$

$\frac{2}{\mathrm{\Delta }}[R(R_c+b\cos {\theta }_H){]}^{1/2},$

Huré, et al. (2020), §2, p. 5826, Eq. (4)

and, as stated in the first sentence of their §3, reference may also be made to the complementary modulus,

$k{\text{'}}_H$

$\equiv$

$[1-k_H^2{]}^{1/2}.$

(Again, we have affixed the subscript "H" in order to differentiate from the modulus employed by Wong (1973).) And in the case of an infinitesimally thin hoop $(b\to 0)$,

$[k_H^2{]}_0$

$=$

$\frac{4R{R}_c}{{\mathrm{\Delta }}_0^2}.$

Huré, et al. (2020), §3, p. 5827, Eq. (12)

Key Finding

On an initial reading, it appears as though the most relevant section of the Huré, et al. (2020) paper is §8 titled, The Solid Torus. They write the gravitational potential in terms of the series expansion,

${\mathrm{\Psi }}_{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{v}}(\overrightarrow{r})$

$\approx$

${\mathrm{\Psi }}_0+\sum _{n=1}^N{\mathrm{\Psi }}_n,$

Huré, et al. (2020), §7, p. 5831, Eq. (42)

where, after setting $M_{\mathrm{t}\mathrm{o}\mathrm{t}}=2{\pi }^2{\rho }_0{b}^2{R}_c$ and acknowledging that $V_{0,0}=1,$ we can write,

${\mathrm{\Psi }}_0$

$=$

$-\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{r}[\frac{r}{{\mathrm{\Delta }}_0}\cdot \frac{2}{\pi }\mathit{K}([k_H{]}_0)]$

Huré, et al. (2020), §8, p. 5832, Eqs. (52) & (53)

and,

$\frac{1}{e^2}[{\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2]$

$=$

$-\frac{G\pi {\rho }_0{R}_c{b}^2}{4(k{\text{'}}_H{)}^2{\mathrm{\Delta }}_0^3}\{[{\mathrm{\Delta }}_0^2-2R_c(R_c+R)]\mathit{E}(k_H)-(k{\text{'}}_H{)}^2{\mathrm{\Delta }}_0^2\mathit{K}(k_H)\}.$

Huré, et al. (2020), §8, p. 5832, Eq. (54)

Rewriting this last expression in a form that can more readily be compared with Wong's work, we obtain,

$\frac{2^3\pi }{e^2}\left[\frac{{\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2}{GM}\right]$	$=$	$-\frac{1}{(k{\text{'}}_H{)}^2{\mathrm{\Delta }}_0^3}\{[{\mathrm{\Delta }}_0^2-2R_c(R_c+R)]\mathit{E}(k_H)-(k{\text{'}}_H{)}^2{\mathrm{\Delta }}_0^2\mathit{K}(k_H)\}$
	$=$	$-\left[\frac{{\mathrm{\Delta }}_0^2-2R_c(R_c+R)}{(1-k_H^2){\mathrm{\Delta }}_0^2}\right]\frac{\mathit{E}(k_H)}{{\mathrm{\Delta }}_0}+\frac{\mathit{K}(k_H)}{{\mathrm{\Delta }}_0}.$

Hence, also,

$\frac{{\mathrm{\Psi }}_0}{GM}+\left[\frac{{\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2}{GM}\right]$	$=$	$-\frac{2}{\pi }\left\{\frac{\mathit{K}([k_H{]}_0)}{{\mathrm{\Delta }}_0}\right\}+\frac{e^2}{2^3\pi }\{\frac{\mathit{K}(k_H)}{{\mathrm{\Delta }}_0}-[\frac{{\mathrm{\Delta }}_0^2-2R_c(R_c+R)}{(1-k_H^2){\mathrm{\Delta }}_0^2}\left]\frac{\mathit{E}(k_H)}{{\mathrm{\Delta }}_0}\right\}$
	$=$	$-\frac{2}{\pi {\mathrm{\Delta }}_0}\{\mathit{K}([k_H{]}_0)-\frac{e^2}{2^4}\cdot \mathit{K}(k_H)\}-\frac{2}{\pi {\mathrm{\Delta }}_0}\cdot \frac{e^2}{2^4}\{\left[\frac{{\mathrm{\Delta }}_0^2-2R_c(R_c+R)}{(1-k_H^2){\mathrm{\Delta }}_0^2}\right]\}\mathit{E}(k_H)$
$\Rightarrow -\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Psi }}_0+{\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2}{GM}\right]$	$=$	$\mathit{K}([k_H{]}_0)-\frac{e^2}{2^4}\cdot \mathit{K}(k_H)+\frac{e^2}{2^4}[\frac{{\mathrm{\Delta }}_0^2-2R_c(R_c+R)}{(1-k_H^2){\mathrm{\Delta }}_0^2}]\mathit{E}(k_H).$

Compare First Terms

Rewriting the first term in the Huré, et al. (2020) series expression for the potential, we have,

$\frac{{\mathrm{\Psi }}_0}{GM}$

$=$

$-\frac{2}{\pi }\left\{\frac{\mathit{K}([k_H{]}_0)}{[({\varpi }_W+R_c{)}^2+z_W^2{]}^{1/2}}\right\},$

where,

$[k_H{]}_0$

$=$

$[\frac{4{\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}{]}^{1/2}=\{\frac{4{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}{\}}^{1/2}.$

For comparison, the first term in Wong's expression is,

$\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}$

$=$

$-\left(\frac{2^3}{3{\pi }^3}\right){\mathrm{{\rm Y}}}_{W0}({\eta }_0)\left\{\frac{\mathit{K}(k)}{r_1}\right\},$

where,

$a^2$	$\equiv$	$R^2-d^2\Rightarrow a=R_c(1-e^2{)}^{1/2},$
$r_1^2$	$\equiv$	$[\varpi +R_c(1-e^2{)}^{1/2}{]}^2+[z-Z_0{]}^2,$
$k$	$\equiv$	$\{\frac{4\varpi R_c(1-e^2{)}^{1/2}}{[\varpi +R_c(1-e^2{)}^{1/2}{]}^2+[z-Z_0{]}^2}{\}}^{1/2},$
${\mathrm{{\rm Y}}}_{W0}({\eta }_0)$	$=$	$\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[{\cosh }^{}{\eta }_0+1]-E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\},$
	$=$	$\frac{(1-e^2{)}^{1/2}}{e^2}\{-K(k_0)\cdot K(k_0)(1-e)+2K(k_0)\cdot E(k_0)(1+e^2)-E(k_0)\cdot E(k_0)(1+e)\},$
$k_0$	$=$	$[\frac{2}{1+1/e}{]}^{1/2}=[\frac{2e}{1+e}{]}^{1/2}.$

This expression is correct for any value of the aspect ratio, $e$. But let's set $Z_0=0$ — as Huré, et al. (2020) have done — then see how the expression simplifies for an infinitesimally thin hoop, that is, if we let $e\to 0$. First we note that,

$k{|}_{e\to 0}$

$=$

$\{\frac{4\varpi R_c}{[\varpi +R_c{]}^2+z^2}{\}}^{1/2},$

so in this limit the modulus of the complete elliptic integral of the first kind becomes identical to the modulus used by Huré, et al. (2020), $[k_H{]}_0$. Next, we note that,

$r_1{|}_{e\to 0}$

$=$

$[(\varpi +R_c{)}^2+z^2{]}^{1/2}.$

As a result, we can write,

$\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}{|}_{e\to 0}$

$=$

$\frac{{\mathrm{\Psi }}_0}{GM}\cdot \left[\right(\frac{2^2}{3{\pi }^2}){\mathrm{{\rm Y}}}_{W0}({\eta }_0){]}_{e\to 0}.$

Now let's evaluate the coefficient, ${\mathrm{{\rm Y}}}_{W0}$, in the limit of $e\to 0$.

${\mathrm{{\rm Y}}}_{W0}$, in the limit of $e\to 0$. First, drawing from our separate examination of the behavior of complete elliptic integral functions, we appreciate that, $\frac{2^2}{{\pi }^2}[K(k_0)\cdot E(k_0)]$ $=$ $1+\frac{1}{2^5}{k}_0^4+\frac{1}{2^5}{k}_0^6+{𝒪}(k_0^8),$ $\frac{2^2}{{\pi }^2}[K(k_0)\cdot K(k_0)]$ $=$ $1+\frac{1}{2}{k}_0^2+\frac{11}{2^5}{k}_0^4+\frac{17}{2^6}{k}_0^6+{𝒪}(k_0^8),$ $\frac{2^2}{{\pi }^2}[E(k_0)\cdot E(k_0)]$ $=$ $1-\frac{1}{2}{k}_0^2-\frac{1}{2^5}{k}_0^4-\frac{1}{2^6}{k}_0^6+{𝒪}(k_0^8).$ Next, employing the binomial expansion, we find that, $k_0^2$ $=$ $2e(1+e{)}^{-1}$   $=$ $2e(1-e+e^2-e^3+e^4-e^5+\cdots );$ and, $k_0^4$ $=$ $4e^2(1+e{)}^{-2}$   $=$ $4e^2(1-2e+3e^2-4e^3+5e^4-6e^5+\cdots ).$ Hence, we have, $\frac{2^2}{{\pi }^2}\left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W0}({\eta }_0)$ $=$ $-(1-e)[1+\frac{1}{2}{k}_0^2+\frac{11}{2^5}{k}_0^4+{𝒪}(k_0^6)]$     $+2(1+e^2)[1+\frac{1}{2^5}{k}_0^4+{𝒪}(k_0^6)]$     $-(1+e)[1-\frac{1}{2}{k}_0^2-\frac{1}{2^5}{k}_0^4+{𝒪}(k_0^6)]$   $=$ $-(1-e)[1+e(1-e+e^2-e^3+e^4-e^5+\cdots )+\frac{11}{2^3}\cdot e^2(1-2e+3e^2-4e^3+5e^4-6e^5+\cdots )+{𝒪}(e^3)]$     $+2(1+e^2)[1+\frac{1}{2^3}\cdot e^2(1-2e+3e^2-4e^3+5e^4-6e^5+\cdots )+{𝒪}(e^3)]$     $-(1+e)[1-e(1-e+e^2-e^3+e^4-e^5+\cdots )-\frac{1}{2^3}\cdot e^2(1-2e+3e^2-4e^3+5e^4-6e^5+\cdots )+{𝒪}(e^3)]$   $=$ $-(1-e)[1+e(1-e)+\frac{11}{2^3}\cdot e^2]+2(1+e^2)[1+\frac{1}{2^3}\cdot e^2]-(1+e)[1-e(1-e)-\frac{1}{2^3}\cdot e^2]+{𝒪}(e^3)$   $=$ $-[1+e(1-e)+\frac{11}{2^3}\cdot e^2]+e[1+e]+2[1+\frac{1}{2^3}\cdot e^2]+2e^2-[1-e(1-e)-\frac{1}{2^3}\cdot e^2]-e[1-e]+{𝒪}(e^3)$   $=$ $-1-e+e^2-\frac{11}{2^3}\cdot e^2+e+e^2+2+\frac{1}{2^2}\cdot e^2+2e^2-1+e-e^2+\frac{1}{2^3}\cdot e^2-e+e^2+{𝒪}(e^3)$   $=$ $\frac{e^2}{2^3}[2^5-11+3]+{𝒪}(e^3)$   $=$ $3e^2+{𝒪}(e^3)$ $\Rightarrow \frac{2^2}{3{\pi }^2}{\mathrm{{\rm Y}}}_{W0}({\eta }_0)$ $=$ $[1+{𝒪}(e^1)]\cdot (1-e^2{)}^{1/2}$ $\Rightarrow [\frac{2^2}{3{\pi }^2}{\mathrm{{\rm Y}}}_{W0}({\eta }_0){]}_{e\to 0}$ $=$ $1.$

Given that,

$[\frac{2^2}{3{\pi }^2}{\mathrm{{\rm Y}}}_{W0}({\eta }_0){]}_{e\to 0}$

$=$

$1,$

we conclude that,

$\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}{|}_{e\to 0}$

$=$

$\frac{{\mathrm{\Psi }}_0}{GM},$

that is, we conclude that ${\mathrm{\Psi }}_0$ matches ${\mathrm{\Phi }}_{W0}$ in the limit of, $e\to 0$.

Go to Higher Order

Let's keep higher order terms in Wong's n = 0 component, and let's examine contributions to the same order that come from Wong's n = 1 and (if necessary) n = 2 components.

First, note that,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$

$=$

${\mathrm{\Delta }}_0\left(\frac{2^2}{3{\pi }^2}\right){\mathrm{{\rm Y}}}_{W0}({\eta }_0)\left\{\frac{\mathit{K}(k)}{r_1}\right\}.$

Keeping Higher Order in Wong's First Component

$\frac{2^2}{{\pi }^2}[K(k_0)\cdot E(k_0)]$	$=$	$1+\frac{1}{2^5}{k}_0^4+\frac{1}{2^5}{k}_0^6+\frac{231}{2^{13}}{k}_0^8+{𝒪}(k_0^{10}),$
$\frac{2^2}{{\pi }^2}[K(k_0)\cdot K(k_0)]$	$=$	$1+\frac{1}{2}{k}_0^2+\frac{11}{2^5}{k}_0^4+\frac{17}{2^6}{k}_0^6+\frac{1787}{2^{13}}{k}^8+{𝒪}(k_0^{10}),$
$\frac{2^2}{{\pi }^2}[E(k_0)\cdot E(k_0)]$	$=$	$1-\frac{1}{2}{k}_0^2-\frac{1}{2^5}{k}_0^4-\frac{1}{2^6}{k}_0^6-\frac{77}{2^{13}}{k}_0^8+{𝒪}(k_0^{10}).$

Add One Additional Term

$\frac{2^2}{{\pi }^2}[K(k)\cdot E(k)]$ $=$ $\{1+(\frac{1}{2}{)}^2{k}^2+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6+(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2{k}^8+\cdots +[\frac{(2n-1)!!}{2^{n}n!}{]}^2{k}^{2n}+\cdots \}$     $\times \{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4-(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2\frac{k^6}{5}-(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2\frac{k^8}{7}-\cdots [\frac{(2n-1)!!}{2^{n}n!}{]}^2\frac{k^{2n}}{2n-1}-\cdots \}$   $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4-(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2\frac{k^6}{5}-\left(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2\frac{k^8}{7}\right\}+\left\{\right(\frac{1}{2}{)}^2{k}^2\}\times \{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4-\left(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2\frac{k^6}{5}\right\}$     $+\left\{\right(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4\}\times \{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4\}+\{\left(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6\right\}\times \{1-\frac{1}{2^2}{k}^2\}+(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2{k}^8+{𝒪}(k^{10})$   $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-(\frac{5}{2^8})k^6-(\frac{5^2\cdot 7}{2^{14}}\left)k^8\right\}+\left\{\right(\frac{1}{2^2})k^2-\frac{1}{2^4}{k}^4-\frac{3}{2^8}{k}^6-\frac{5}{2^{10}}{k}^8\}$     $+\left(\frac{3^2}{2^6}\right)k^4-\left(\frac{3^2}{2^8}\right)k^6-\left(\frac{3^3}{2^{12}}\right)k^8+\left(\frac{5^2}{2^8}\right)k^6-\left(\frac{5^2}{2^{10}}\right)k^8+\left(\frac{5^2\cdot 7^2}{2^{14}}\right)k^8+{𝒪}(k^{10})$   $=$ $1+[\frac{1}{2^2}-\frac{1}{2^2}]k^2+[\frac{3^2}{2^6}-\frac{3}{2^6}-\frac{1}{2^4}]k^4+[\frac{5^2}{2^8}-\frac{5}{2^8}-\frac{3}{2^8}-\frac{3^2}{2^8}]k^6+\left[\right(\frac{5^2\cdot 7^2}{2^{14}})-(\frac{5^2}{2^{10}})-(\frac{3^3}{2^{12}})-\frac{5}{2^{10}}-(\frac{5^2\cdot 7}{2^{14}}\left)\right]k^8+{𝒪}(k^{10})$   $=$ $1+\frac{1}{2^5}{k}^4+\frac{1}{2^5}{k}^6+\frac{231}{2^{13}}{k}^8+{𝒪}(k^{10})$ $\frac{2^2}{{\pi }^2}[K(k)\cdot K(k)]$ $=$ $\{1+(\frac{1}{2}{)}^2{k}^2+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6+(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2{k}^8+\cdots +[\frac{(2n-1)!!}{2^{n}n!}{]}^2{k}^{2n}+\cdots \}$     $\times \{1+(\frac{1}{2}{)}^2{k}^2+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6+(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2{k}^8+\cdots +[\frac{(2n-1)!!}{2^{n}n!}{]}^2{k}^{2n}+\cdots \}$   $=$ $\{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4+\frac{5^2}{2^8}{k}^6+\frac{5^2\cdot 7^2}{2^{14}}{k}^8\}\times \{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4+\frac{5^2}{2^8}{k}^6+\frac{5^2\cdot 7^2}{2^{14}}{k}^8\}+{𝒪}(k^{10})$   $=$ $\{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4+\frac{5^2}{2^8}{k}^6+\frac{5^2\cdot 7^2}{2^{14}}{k}^8\}+\frac{1}{2^2}{k}^2\{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4+\frac{5^2}{2^8}{k}^6\}+\frac{3^2}{2^6}{k}^4\{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4\}+\left\{\frac{5^2}{2^8}{k}^6\right\}\{1+\frac{1}{2^2}{k}^2\}+\left\{\frac{5^2\cdot 7^2}{2^{14}}{k}^8\right\}+{𝒪}(k^{10})$   $=$ $\{1+\frac{1}{2^2}{k}^2+\frac{3^2}{2^6}{k}^4+\frac{5^2}{2^8}{k}^6+\frac{5^2\cdot 7^2}{2^{14}}{k}^8\}+\{\frac{1}{2^2}{k}^2+\frac{1}{2^4}{k}^4+\frac{3^2}{2^8}{k}^6+\frac{5^2}{2^{10}}{k}^8\}+\{\frac{3^2}{2^6}{k}^4+\frac{3^2}{2^8}{k}^6+\frac{3^4}{2^{12}}{k}^8\}+\{\frac{5^2}{2^8}{k}^6+\frac{5^2}{2^{10}}{k}^8\}+\left\{\frac{5^2\cdot 7^2}{2^{14}}{k}^8\right\}+{𝒪}(k^{10})$   $=$ $1+[\frac{1}{2^2}+\frac{1}{2^2}]k^2+[\frac{3^2}{2^6}+\frac{1}{2^4}+\frac{3^2}{2^6}]k^4+[\frac{5^2}{2^8}+\frac{3^2}{2^8}+\frac{3^2}{2^8}+\frac{5^2}{2^8}]k^6+[\frac{5^2\cdot 7^2}{2^{14}}+\frac{5^2}{2^{10}}+\frac{3^4}{2^{12}}+\frac{5^2}{2^{10}}+\frac{5^2\cdot 7^2}{2^{14}}]k^8+{𝒪}(k^{10})$   $=$ $1+\frac{1}{2}{k}^2+\frac{11}{2^5}{k}^4+\frac{17}{2^6}{k}^6+\frac{1787}{2^{13}}{k}^8+{𝒪}(k^{10})$ $\frac{2^2}{{\pi }^2}[E(k)\cdot E(k)]$ $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4-(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2\frac{k^6}{5}-(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2\frac{k^8}{7}-\cdots [\frac{(2n-1)!!}{2^{n}n!}{]}^2\frac{k^{2n}}{2n-1}-\cdots \}$     $\times \{1-\frac{1}{2^2}{k}^2-\frac{1^2\cdot 3}{2^2\cdot 4^2}{k}^4-(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2\frac{k^6}{5}-(\frac{1\cdot 3\cdot 5\cdot 7}{2^7\cdot 3}{)}^2\frac{k^8}{7}-\cdots [\frac{(2n-1)!!}{2^{n}n!}{]}^2\frac{k^{2n}}{2n-1}-\cdots \}$   $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-\frac{5}{2^8}{k}^6-\frac{5^2\cdot 7}{2^{14}}{k}^8\}\times \{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-\frac{5}{2^8}{k}^6-\frac{5^2\cdot 7}{2^{14}}{k}^8\}+{𝒪}(k^{10})$   $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-\frac{5}{2^8}{k}^6-\frac{5^2\cdot 7}{2^{14}}{k}^8\}+\{-\frac{1}{2^2}{k}^2\}\times \{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-\frac{5}{2^8}{k}^6\}$     $+\{-\frac{3}{2^6}{k}^4\}\times \{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4\}+\{-\frac{5}{2^8}{k}^6\}\times \{1-\frac{1}{2^2}{k}^2\}+\{-\frac{5^2\cdot 7}{2^{14}}{k}^8\}+{𝒪}(k^{10})$   $=$ $\{1-\frac{1}{2^2}{k}^2-\frac{3}{2^6}{k}^4-\frac{5}{2^8}{k}^6-\frac{5^2\cdot 7}{2^{14}}{k}^8\}+\{-\frac{1}{2^2}{k}^2+\frac{1}{2^4}{k}^4+\frac{3}{2^8}{k}^6+\frac{5}{2^{10}}{k}^8\}$     $+\{-\frac{3}{2^6}{k}^4+\frac{3}{2^8}{k}^6+\frac{3^2}{2^{12}}{k}^8\}+\{-\frac{5}{2^8}{k}^6+\frac{5}{2^{10}}{k}^8\}+\{-\frac{5^2\cdot 7}{2^{14}}{k}^8\}+{𝒪}(k^{10})$ $$ $=$ $1+[-\frac{1}{2^2}-\frac{1}{2^2}]k^2+[-\frac{3}{2^6}+\frac{1}{2^4}-\frac{3}{2^6}]k^4+[-\frac{5}{2^8}+\frac{3}{2^8}+\frac{3}{2^8}-\frac{5}{2^8}]k^6+[-\frac{5^2\cdot 7}{2^{14}}+\frac{5}{2^{10}}+\frac{3^2}{2^{12}}+\frac{5}{2^{10}}-\frac{5^2\cdot 7}{2^{14}}]k^8+{𝒪}(k^{10})$ $$ $=$ $1+[-\frac{2}{2^2}]k^2+[-\frac{3}{2^5}+\frac{2}{2^5}]k^4+[-\frac{5}{2^7}+\frac{3}{2^7}]k^6+[-\frac{5^2\cdot 7}{2^{13}}+\frac{5\cdot 2^4}{2^{13}}+\frac{2\cdot 3^2}{2^{13}}]k^8+{𝒪}(k^{10})$   $=$ $1-\frac{1}{2}{k}^2-\frac{1}{2^5}{k}^4-\frac{1}{2^6}{k}^6-\frac{77}{2^{13}}{k}^8+{𝒪}(k^{10})$

Next, employing the binomial expansion, we find that,

$k_0^2$	$=$	$2e(1+e{)}^{-1}$
	$=$	$2e(1-e+e^2-e^3+e^4-e^5+\cdots );$
$k_0^4$	$=$	$4e^2(1+e{)}^{-2}$
	$=$	$4e^2(1-2e+3e^2-4e^3+5e^4-6e^5+\cdots );$
$k_0^6$	$=$	$2^3{e}^3(1+e{)}^{-3}$
	$=$	$2^3{e}^3(1-3e+6e^2-10e^3+15e^4-21e^5+\cdots );$
$k_0^8$	$=$	$2^4{e}^4(1+e{)}^{-4}$
	$=$	$2^4{e}^4(1-4e+10e^2-20e^3+35e^4-56e^5+\cdots ).$

Hence, we have,

$\frac{2^2}{{\pi }^2}\left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W0}({\eta }_0)$	$=$	$-(1-e)[1+\frac{1}{2}{k}_0^2+\frac{11}{2^5}{k}_0^4+\frac{17}{2^6}{k}_0^6+\frac{1787}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
		$+2(1+e^2)[1+\frac{1}{2^5}{k}_0^4+\frac{1}{2^5}{k}_0^6+\frac{231}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
		$-(1+e)[1-\frac{1}{2}{k}_0^2-\frac{1}{2^5}{k}_0^4-\frac{1}{2^6}{k}_0^6-\frac{77}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
	$=$	$-(1-e)[1+e(1-e+e^2-e^3)+\frac{11}{2^3}\cdot e^2(1-2e+3e^2)+\frac{17}{2^6}\cdot 2^3{e}^3(1-3e)+\frac{1787}{2^{13}}\cdot 2^4{e}^4+{𝒪}(e^5)]$
		$+2(1+e^2)[1+\frac{1}{2^5}\cdot 4e^2(1-2e+3e^2)+\frac{1}{2^5}{2}^3{e}^3(1-3e)+\frac{231}{2^{13}}\cdot 2^4{e}^4+{𝒪}(e^5)]$
		$-(1+e)[1-e(1-e+e^2-e^3)-\frac{1}{2^5}\cdot 4e^2(1-2e+3e^2)-\frac{1}{2^6}{2}^3{e}^3(1-3e)-\frac{77}{2^{13}}{2}^4{e}^4+{𝒪}(e^5)]$
	$=$	$-2^{-9}(1-e)[512(1+e-e^2+e^3-e^4)+704\cdot (e^2-2e^3+3e^4)+1088\cdot (e^3-3e^4)+1787\cdot e^4+{𝒪}(e^5)]$
		$+(1+e^2)2^{-9}[1024+128\cdot (e^2-2e^3+3e^4)+256(e^3-3e^4)+462\cdot e^4+{𝒪}(e^5)]$
		$-2^{-9}(1+e)[512(1-e+e^2-e^3+e^4)-64\cdot (e^2-2e^3+3e^4)-64(e^3-3e^4)-77e^4+{𝒪}(e^5)]$
	$=$	$-2^{-9}(1-e)[512+512e+192e^2+e^3(512-1408+1088)+e^4(704-512-3264+1787)+{𝒪}(e^5)]$
		$+2^{-9}(1+e^2)[1024+128e^2+e^3(-256+256)+e^4(384-768+462)+{𝒪}(e^5)]$
		$-2^{-9}(1+e)[512-512e+e^2(512-64)+e^3(-512+128-64)+e^4(512-192-192-77)+{𝒪}(e^5)]$
	$=$	$-2^{-9}(1-e)[512+512e+192e^2+192e^3-1285e^4+{𝒪}(e^5)]$
		$+2^{-9}(1+e^2)[1024+128e^2+78e^4+{𝒪}(e^5)]$
		$-2^{-9}(1+e)[512-512e+448e^2-448e^3+51e^4+{𝒪}(e^5)]$
	$=$	$2^{-9}[(-192+128-448)e^2+(-192+448)e^3+(1285+78-51)e^4+{𝒪}(e^5)]$
		$+2^{-9}e[1024e+(192-448)e^2+(192+448)e^3+{𝒪}(e^4)]+2^{-9}{e}^2[1024+128e^2+{𝒪}(e^3)]$
	$=$	$2^{-9}[(-192+128-448)e^2+2048e^2+(1285+78-51)e^4+(192+448)e^4+128e^4+{𝒪}(e^5)]$
	$=$	$2^{-9}[1536e^2+2080e^4+{𝒪}(e^5)]$
	$=$	$3e^2+\left[\frac{5\cdot 13}{2^4}\right]e^4+{𝒪}(e^5)$
$\Rightarrow \frac{2^2}{3{\pi }^2}{\mathrm{{\rm Y}}}_{W0}({\eta }_0)$	$=$	$[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot (1-e^2{)}^{1/2}.$

Hence,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$

$=$

$[1+{𝒪}(e^2)]\cdot (1-e^2{)}^{1/2}\{\frac{\mathit{K}(k)}{r_1}\}{\mathrm{\Delta }}_0,$

or, more precisely,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$

$=$

$[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot (1-e^2{)}^{1/2}\{\frac{\mathit{K}(k)}{r_1}\}{\mathrm{\Delta }}_0.$

Next Factors

Now,

${\mathrm{\Delta }}_0^2$	$=$	$({\varpi }_W+R_c{)}^2+z_W^2,$
$r_1^2$	$=$	$[{\varpi }_W+R_c(1-e^2{)}^{1/2}{]}^2+z_W^2,$
$\Rightarrow r_1^2-{\mathrm{\Delta }}_0^2$	$=$	$[{\varpi }_W+R_c(1-e^2{)}^{1/2}{]}^2+z_W^2-[({\varpi }_W+R_c{)}^2+z_W^2]$
	$=$	${\varpi }_W^2+2{\varpi }_W{R}_c(1-e^2{)}^{1/2}+R_c^2(1-e^2)+z_W^2-[{\varpi }_W^2+2{\varpi }_W{R}_c+R_c^2+z_W^2]$
	$=$	$2{\varpi }_W{R}_c[(1-e^2{)}^{1/2}-1]-e^2{R}_c^2.$

---

Again, drawing from the binomial theorem, we have,

$(1-e^2{)}^{1/2}$	$=$	$1-\frac{1}{2}{e}^2+\left[\frac{\frac{1}{2}(-\frac{1}{2})}{2}\right]e^4-\left[\frac{\frac{1}{2}(-\frac{1}{2})(-\frac{3}{2})}{3!}\right]e^6+\left[\frac{\frac{1}{2}(-\frac{1}{2})(-\frac{3}{2})(-\frac{5}{2})}{4!}\right]e^8+\cdots$
	$=$	$1-\frac{1}{2}{e}^2-\frac{1}{2^3}{e}^4-\frac{1}{2^4}{e}^6-\frac{5}{2^7}{e}^8-{𝒪}(e^{10}).$

---

$\Rightarrow r_1^2-{\mathrm{\Delta }}_0^2$	$=$	${\varpi }_W{R}_c[-e^2-\frac{1}{2^2}{e}^4-\frac{1}{2^3}{e}^6-\frac{5}{2^6}{e}^8-{𝒪}(e^{10})]-e^2{R}_c^2$
$\Rightarrow \frac{r_1^2}{{\mathrm{\Delta }}_0^2}$	$=$	$1-e^2\left[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}\right]-\frac{{\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}[\frac{1}{2^2}{e}^4+\frac{1}{2^3}{e}^6+\frac{5}{2^6}{e}^8+{𝒪}(e^{10})]$

Now Work on Elliptic Integral Expressions

From a separate discussion we can draw the series expansion of $\mathit{K}(k)$, specifically,

$\frac{2K(k)}{\pi }$

$=$

$1+(\frac{1}{2}{)}^2{k}^2+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6+\cdots$

where,

$\frac{k^2}{4}$	$\equiv$	$\frac{1}{4}\left[\frac{r_1^2-r_2^2}{r_1^2}\right]=\left[\frac{a\varpi }{r_1^2}\right]=\left[\frac{a\varpi }{(\varpi +a{)}^2+(z-Z_0{)}^2}\right]$
	$=$	$\frac{{\varpi }_W{R}_c(1-e^2{)}^{1/2}}{[{\varpi }_W+R_c(1-e^2{)}^{1/2}{]}^2+z_W^2}.$

Also,

$\frac{2K(k_H)}{\pi }$

$=$

$1+(\frac{1}{2}{)}^2{k}_H^2+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}_H^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}_H^6+\cdots$

where,

$\frac{k_H^2}{4}$	$=$	$\frac{1}{{\mathrm{\Delta }}^2}[R(R_c+b\cos {\theta }_H)]$
	$=$	$\frac{{\varpi }_W(R_c+b\cos {\theta }_H)}{[{\varpi }_W+(R_c+b\cos {\theta }_H){]}^2+[z_W-b\sin {\theta }_H{]}^2}.$

What we want to do is write $K(k)$ in terms of $K(k_H)$. Let's try …

$\frac{2K(k)}{\pi }$

$=$

$\frac{2K(k_H)}{\pi }+{\delta }_K,$

where,

${\delta }_K$	$\equiv$	$\frac{2K(k)}{\pi }-\frac{2K(k_H)}{\pi }$
	$=$	$\{1+\frac{k^2}{4}+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}^6+\cdots \}-\{1+\frac{k_H^2}{4}+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}_H^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}_H^6+\cdots \}$
	$\approx$	$\{1+\frac{k^2}{4}\}-\{1+\frac{k_H^2}{4}\}=\frac{k^2}{4}-\frac{k_H^2}{4}$
	$=$	$\frac{{\varpi }_W{R}_c(1-e^2{)}^{1/2}}{[{\varpi }_W+R_c(1-e^2{)}^{1/2}{]}^2+z_W^2}-\frac{{\varpi }_W(R_c+b\cos {\theta }_H)}{[{\varpi }_W+(R_c+b\cos {\theta }_H){]}^2+[z_W-b\sin {\theta }_H{]}^2}$
	$=$	$\{{\varpi }_W{R}_c(1-e^2{)}^{1/2}\left\}\right\{{\varpi }_W^2+R_c^2+z_W^2+2{\varpi }_W{R}_c(1-e^2{)}^{1/2}-R_c^2{e}^2{\}}^{-1}$
		$-\{{\varpi }_W{R}_c(1+e\cos {\theta }_H)\}\{[{\varpi }_W+R_c(1+e\cos {\theta }_H){]}^2+[z_W-R_{c}e\sin {\theta }_H{]}^2{\}}^{-1}$
	$\approx$	$\left\{{\varpi }_W{R}_c\right[1-\frac{1}{2}{e}^2+{𝒪}(e^4)\left]\right\}\{{\varpi }_W^2+R_c^2+z_W^2+2{\varpi }_W{R}_c[1-\frac{1}{2}{e}^2+{𝒪}(e^4)]-R_c^2{e}^2{\}}^{-1}$
		$-\left\{{\varpi }_W{R}_c\right[1+e\cos {\theta }_H\left]\right\}\{{\varpi }_W^2+2{\varpi }_W{R}_c(1+e\cos {\theta }_H)+R_c^2(1+2e\cos {\theta }_H+e^2{\cos }^{}{\theta }_H)+z_W^2-2z_W{R}_{c}e\sin {\theta }_H+R_c^2{e}^2\sin {\theta }_H^2{\}}^{-1}$
	$\approx$	$\{{\varpi }_W{R}_c-e^2[\frac{{\varpi }_W{R}_c}{2}]+{𝒪}(e^4)\}\{[({\varpi }_W+R_c{)}^2+z_W^2]-e^2[{\varpi }_W{R}_c+R_c^2]+{𝒪}(e^4){\}}^{-1}$
		$-\{{\varpi }_W{R}_c+e[{\varpi }_W{R}_c\cos {\theta }_H\left]\right\}\{[({\varpi }_W^2+R_c{)}^2+z_W^2]+2R_{c}e(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)+R_c^2{e}^2{\}}^{-1}$
	$\approx$	$\left\{\frac{{\varpi }_W{R}_c}{[({\varpi }_W+R_c{)}^2+z_W^2]}\right[1-\frac{1}{2}{e}^2\left]\right\}\{1-e^2[\frac{{\varpi }_W{R}_c+R_c^2}{({\varpi }_W+R_c{)}^2+z_W^2}]{\}}^{-1}$
		$-\left\{\frac{{\varpi }_W{R}_c}{[({\varpi }_W^2+R_c{)}^2+z_W^2]}\right[1+e\cos {\theta }_H\left]\right\}\{1+e[\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{({\varpi }_W^2+R_c{)}^2+z_W^2}]+e^2[\frac{R_c^2}{({\varpi }_W^2+R_c{)}^2+z_W^2}]{\}}^{-1}$
	$\approx$	$\frac{{\varpi }_W{R}_c}{[({\varpi }_W+R_c{)}^2+z_W^2]}[1-\frac{1}{2}{e}^2]\{1+e^2[\frac{{\varpi }_W{R}_c+R_c^2}{({\varpi }_W+R_c{)}^2+z_W^2}\left]\right\}$
		$-\frac{{\varpi }_W{R}_c}{[({\varpi }_W^2+R_c{)}^2+z_W^2]}[1+e\cos {\theta }_H]\{1-e[\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{({\varpi }_W^2+R_c{)}^2+z_W^2}]-e^2[\frac{R_c^2}{({\varpi }_W^2+R_c{)}^2+z_W^2}\left]\right\}$
	$\approx$	$-\frac{e\cdot {\varpi }_W{R}_c}{[({\varpi }_W^2+R_c{)}^2+z_W^2]}[\cos {\theta }_H-\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{({\varpi }_W^2+R_c{)}^2+z_W^2}]$
		$+\frac{e^2\cdot {\varpi }_W{R}_c}{[({\varpi }_W+R_c{)}^2+z_W^2]}\left\{\right[\frac{{\varpi }_W{R}_c+R_c^2}{({\varpi }_W+R_c{)}^2+z_W^2}-\frac{1}{2}]-[\frac{R_c^2}{({\varpi }_W^2+R_c{)}^2+z_W^2}]-\cos {\theta }_H[\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{({\varpi }_W^2+R_c{)}^2+z_W^2}\left]\right\}$
	$\approx$	$-\frac{e\cdot {\varpi }_W{R}_c}{[({\varpi }_W^2+R_c{)}^2+z_W^2][({\varpi }_W^2+R_c{)}^2+z_W^2]}[\cos {\theta }_H[({\varpi }_W^2+R_c{)}^2+z_W^2]-2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)]$
		$+\frac{\frac{1}{2}{e}^2\cdot {\varpi }_W{R}_c}{[({\varpi }_W+R_c{)}^2+z_W^2][({\varpi }_W^2+R_c{)}^2+z_W^2]}\{2({\varpi }_W{R}_c+R_c^2)-[({\varpi }_W^2+R_c{)}^2+z_W^2]-2R_c^2-4R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)\}$
	$\approx$	$-\frac{e\cdot {\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}[\cos {\theta }_H-\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{{\mathrm{\Delta }}_0^2}]$
		$+\frac{e^2\cdot {\varpi }_W{R}_c}{2{\mathrm{\Delta }}_0^4}\{2({\varpi }_W{R}_c+R_c^2)-[({\varpi }_W^2+R_c{)}^2+z_W^2]-2R_c^2-4R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)\}.$

Let's subtract $K([k_H{]}_0)$ from the potential expression. But first, let's adopt the shorthand notation …

Given that, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$ $=$ $\{\mathit{K}(k)\}\frac{{\mathrm{\Delta }}_0}{r_1}\left(\frac{2^2}{3{\pi }^2}\right){\mathrm{{\rm Y}}}_{W0}({\eta }_0)$   $=$ $\{\mathit{K}(k)\}\left[\frac{r_1^2}{{\mathrm{\Delta }}_0^2}{]}^{-1/2}\right[1+\left(\frac{5\cdot 13}{2^4\cdot 3}\right)e^2+{𝒪}(e^3)]\cdot (1-e^2{)}^{1/2}$   $=$ $\{\mathit{K}(k)\}\{1+\frac{e^2}{2}[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}]+{𝒪}(e^4)\}[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot [1-\frac{1}{2}{e}^2-{𝒪}(e^4)]$ let's define the variable, ${𝒜}$, such that, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$ $=$ $\{\mathit{K}(k)\}\{1+e^2{𝒜}\}$ $\Rightarrow \{1+e^2\cdot {𝒜}\}$ $=$ $\{1+\frac{e^2}{2}[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}]+{𝒪}(e^4)\}[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot [1-\frac{1}{2}{e}^2-{𝒪}(e^4)]$   $\approx$ $\{1+\frac{e^2}{2}[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}\left]\right\}[1+(\frac{5\cdot 13}{2^4\cdot 3}\left)e^2\right][1-\frac{1}{2}{e}^2]$   $\approx$ $\{1+\frac{e^2}{2}[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}]+(\frac{5\cdot 13}{2^4\cdot 3})e^2-\frac{1}{2}{e}^2\}$ $\Rightarrow 2{𝒜}$ $\approx$ $\left[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}\right]+\left(\frac{41}{2^3\cdot 3}\right).$

We can therefore write,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]-K([k_H{]}_0)$	$\approx$	$-K([k_H{]}_0)+\{K(k_H)+\frac{\pi }{2}\cdot {\delta }_K\}\{1+e^2\cdot {𝒜}\}$
	$\approx$	$-K([k_H{]}_0)+K(k_H)\{1+e^2\cdot {𝒜}\}+\frac{\pi }{2}\cdot {\delta }_K,$

where we should keep in mind that ${\delta }_k$ is ${𝒪}(e^1)$. So, let's examine the piece,

$\frac{2}{\pi }[K(k_H)-K([k_H{]}_0)]$	$=$	$\{1+\frac{k_H^2}{4}+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}_H^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}_H^6+\cdots \}-\{1+\frac{k_H^2}{4}+(\frac{1\cdot 3}{2\cdot 4}{)}^2{k}_H^4+(\frac{1\cdot 3\cdot 5}{2^4\cdot 3}{)}^2{k}_H^6+\cdots {\}}_{e\to 0}$
	$\approx$	$\frac{k_H^2}{4}-[\frac{k_H^2}{4}{]}_{e\to 0}$
	$=$	$\left[\frac{{\varpi }_W{R}_c(1+e\cos {\theta }_H)}{[{\varpi }_W+R_c(1+e\cos {\theta }_H){]}^2+[z_W-R_{c}e\sin {\theta }_H{]}^2}\right]-[\frac{{\varpi }_W{R}_c(1+e\cos {\theta }_H)}{[{\varpi }_W+R_c(1+e\cos {\theta }_H){]}^2+[z_W-R_{c}e\sin {\theta }_H{]}^2}{]}_{e\to 0}$
	$=$	${\varpi }_W{R}_c(1+e\cos {\theta }_H)\{[{\varpi }_W+R_c(1+e\cos {\theta }_H){]}^2+[z_W-R_{c}e\sin {\theta }_H{]}^2{\}}^{-1}-[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}]$
	$=$	${\varpi }_W{R}_c(1+e\cos {\theta }_H)\{{\varpi }_W^2+2{\varpi }_W{R}_c(1+e\cos {\theta }_H)+R_c^2(1+2e\cos {\theta }_H+e^2{\cos }^{}{\theta }_H)+z_W^2-2z_W{R}_{c}e\sin {\theta }_H+R_c^2{e}^2{\sin }^{}{\theta }_H{\}}^{-1}$
		$-\left[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}\right]$
	$=$	$-\left[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}\right]+\left[\frac{{\varpi }_W{R}_c(1+e\cos {\theta }_H)}{({\varpi }_W+R_c{)}^2+z_W^2}\right]\{1+\frac{e[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}+\frac{e^2[R_c^2{\cos }^2{\theta }_H+R_c^2{\sin }^2{\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}{\}}^{-1}$
	$\approx$	$-\left[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}\right]+\left[\frac{{\varpi }_W{R}_c(1+e\cos {\theta }_H)}{({\varpi }_W+R_c{)}^2+z_W^2}\right]\{1-\frac{e[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}$
		$-\frac{e^2[R_c^2{\cos }^2{\theta }_H+R_c^2{\sin }^2{\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}+\frac{e^2[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{[({\varpi }_W+R_c{)}^2+z_W^2{]}^2}\}$
	$\approx$	$\left[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}\right]\{-\frac{e[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}$
		$-\frac{e^2[R_c^2{\cos }^2{\theta }_H+R_c^2{\sin }^2{\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}+\frac{e^2[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{[({\varpi }_W+R_c{)}^2+z_W^2{]}^2}\}$
		$+\left[\frac{{\varpi }_W{R}_c}{({\varpi }_W+R_c{)}^2+z_W^2}\right]\left\{\frac{e\cos {\theta }_H[({\varpi }_W+R_c{)}^2+z_W^2]-e^2\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{({\varpi }_W+R_c{)}^2+z_W^2}\right\}$
	$\approx$	$\left[\frac{{\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}\right]\{e\cos {\theta }_H-\frac{e{R}_c\cos {\theta }_H[2{\varpi }_W+2R_c-2z_W\tan {\theta }_H]}{{\mathrm{\Delta }}_0^2}-\frac{e^2{R}_c^2}{{\mathrm{\Delta }}_0^2}$
		$-\frac{e^2\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{{\mathrm{\Delta }}_0^2}+\frac{e^2[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{{\mathrm{\Delta }}_0^4}\}.$

Now we have,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]-K([k_H{]}_0)$	$\approx$	$-K([k_H{]}_0)+K(k_H)\{1+e^2\cdot {𝒜}\}+\frac{\pi }{2}\cdot {\delta }_K$
$\Rightarrow -{\mathrm{\Delta }}_0\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$	$\approx$	$\frac{2}{\pi }K([k_H{]}_0)+\frac{2}{\pi }[K(k_H)-K([k_H{]}_0)]+{\delta }_K+\frac{2}{\pi }K(k_H)e^2\cdot {𝒜}.$

But, as we have just demonstrated,

$\frac{2}{\pi }[K(k_H)-K([k_H{]}_0)]+{\delta }_K$	$\approx$	$\left[\frac{{\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}\right]\{e\cos {\theta }_H-\frac{e{R}_c\cos {\theta }_H[2{\varpi }_W+2R_c-2z_W\tan {\theta }_H]}{{\mathrm{\Delta }}_0^2}-\frac{e^2{R}_c^2}{{\mathrm{\Delta }}_0^2}$
		$-\frac{e^2\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]}{{\mathrm{\Delta }}_0^2}+\frac{e^2[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{{\mathrm{\Delta }}_0^4}\}$
		$-\frac{e\cdot {\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}[\cos {\theta }_H-\frac{2R_c(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)}{{\mathrm{\Delta }}_0^2}]$
		$+\frac{e^2\cdot {\varpi }_W{R}_c}{2{\mathrm{\Delta }}_0^4}\{2({\varpi }_W{R}_c+R_c^2)-[({\varpi }_W^2+R_c{)}^2+z_W^2]-2R_c^2-4R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)\}.$
	$\approx$	$\left[\frac{e^2\cdot {\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^4}\right]\{-R_c^2-\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]$
		$+({\varpi }_W{R}_c+R_c^2)-\frac{1}{2}[({\varpi }_W^2+R_c{)}^2+z_W^2]-R_c^2-2R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)$
		$+\frac{[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{{\mathrm{\Delta }}_0^2}\}$
TEMPORARY BREAK HERE

Hence,

$-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$	$\approx$	$K([k_H{]}_0)+K(k_H)e^2\cdot {𝒜}+\frac{\pi }{2}[\frac{e^2\cdot {\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^4}\left]\right\{-R_c^2-\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]$
		$+({\varpi }_W{R}_c+R_c^2)-\frac{1}{2}[({\varpi }_W^2+R_c{)}^2+z_W^2]-R_c^2-2R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)$
		$+\frac{[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{{\mathrm{\Delta }}_0^2}\}$

Include Second Wong Term

$\left(\frac{a}{GM}\right){\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}(\varpi ,z){|}_{\mathrm{e}\mathrm{x}\mathrm{t}\mathrm{e}\mathrm{r}\mathrm{i}\mathrm{o}\mathrm{r}}$	$=$	$-\left(\frac{2^3}{3{\pi }^3}\right){\mathrm{{\rm Y}}}_{W1}({\eta }_0)\times \cos \theta \{\frac{a}{r_2}\cdot \mathit{E}(k)\}$
$\Rightarrow -\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}}{GM}\right]$	$=$	$\left(\frac{2^2}{3{\pi }^2}\right){\mathrm{{\rm Y}}}_{W1}({\eta }_0)\times \cos \theta \{\frac{{\mathrm{\Delta }}_0}{r_2}\cdot \mathit{E}(k)\};$

Leading (Upsilon) Coefficient

${\mathrm{{\rm Y}}}_{W1}({\eta }_0)$	$\equiv$	$\left[\frac{\sinh {\eta }_0}{\cosh {\eta }_0}\right]\{K(k_0)\cdot K(k_0)[\cosh {\eta }_0(1-\cosh {\eta }_0)]+2K(k_0)\cdot E(k_0)[(3{\cosh }^{}{\eta }_0-1)]-5E(k_0)\cdot E(k_0)[\cosh {\eta }_0(1+\cosh {\eta }_0)]\}$
$\Rightarrow \left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W1}({\eta }_0)$	$=$	$-(1-e)K(k_0)\cdot K(k_0)+2(3-e^2)K(k_0)\cdot E(k_0)-5(1+e)E(k_0)\cdot E(k_0).$
$\Rightarrow \frac{2^2}{{\pi }^2}\left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W1}({\eta }_0)$	$=$	$-(1-e)[1+\frac{1}{2}{k}_0^2+\frac{11}{2^5}{k}_0^4+\frac{17}{2^6}{k}_0^6+\frac{1787}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
		$+2(3-e^2)[1+\frac{1}{2^5}{k}_0^4+\frac{1}{2^5}{k}_0^6+\frac{231}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
		$-5(1+e)[1-\frac{1}{2}{k}_0^2-\frac{1}{2^5}{k}_0^4-\frac{1}{2^6}{k}_0^6-\frac{77}{2^{13}}{k}_0^8+{𝒪}(k_0^{10})]$
	$=$	$-(1-e)[1+\frac{1}{2}\cdot 2e(1-e+e^2-e^3)+\frac{11}{2^5}\cdot 4e^2(1-2e+3e^2)+\frac{17}{2^6}\cdot 2^3{e}^3(1-3e)+\frac{1787}{2^{13}}\cdot 2^4{e}^4+{𝒪}(e^5)]$
		$+2(3-e^2)[1+\frac{1}{2^5}\cdot 4e^2(1-2e+3e^2)+\frac{1}{2^5}\cdot 2^3{e}^3(1-3e)+\frac{231}{2^{13}}\cdot 2^4{e}^4+{𝒪}(e^5)]$
		$-5(1+e)[1-\frac{1}{2}\cdot 2e(1-e+e^2-e^3)-\frac{1}{2^5}\cdot 4e^2(1-2e+3e^2)-\frac{1}{2^6}\cdot 2^3{e}^3(1-3e)-\frac{77}{2^{13}}\cdot 2^4{e}^4+{𝒪}(e^5)]$
	$=$	$-2^{-9}(1-e)[2^9+2^{9}e(1-e+e^2-e^3)+2^6\cdot 11\cdot e^2(1-2e+3e^2)+2^6\cdot 17\cdot e^3(1-3e)+1787\cdot e^4+{𝒪}(e^5)]$
		$+2^{-9}(6-2e^2)[2^9+2^6\cdot e^2(1-2e+3e^2)+2^7\cdot e^3(1-3e)+231\cdot e^4+{𝒪}(e^5)]$
		$-2^{-9}(5+5e)[2^9-2^9\cdot e(1-e+e^2-e^3)-2^6\cdot e^2(1-2e+3e^2)-2^6\cdot e^3(1-3e)-77\cdot e^4+{𝒪}(e^5)]$
$\Rightarrow \frac{2^{11}}{{\pi }^2}\left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W1}({\eta }_0)$	$=$	$-2^9-2^{9}e(1-e+e^2-e^3)-2^6\cdot 11\cdot e^2(1-2e+3e^2)-2^6\cdot 17\cdot e^3(1-3e)-1787\cdot e^4$
		$+2^{9}e+2^9{e}^2(1-e+e^2)+2^6\cdot 11\cdot e^3(1-2e)+2^6\cdot 17\cdot e^4$
		$+6[2^9+2^6\cdot e^2(1-2e+3e^2)+2^7\cdot e^3(1-3e)+231\cdot e^4]-2^{10}{e}^2-2^7\cdot e^4$
		$+5[-2^9+2^9\cdot e(1-e+e^2-e^3)+2^6\cdot e^2(1-2e+3e^2)+2^6\cdot e^3(1-3e)+77\cdot e^4]$
		$+5[-2^{9}e+2^9\cdot e^2(1-e+e^2)+2^6\cdot e^3(1-2e)+2^6\cdot e^4]+{𝒪}(e^5)$
	$=$	$2^9(e^2-e^3+e^4)-2^6\cdot 11(e^2-2e^3+3e^4)-2^6\cdot 17(e^3-3e^4)-1787e^4$
		$+2^9(e^2-e^3+e^4)+2^6\cdot 11(e^3-2e^4)+2^6\cdot 17e^4$
		$+3\cdot 2^7(e^2-2e^3+3e^4)+3\cdot 2^8(e^3-3e^4)+2\cdot 3^2\cdot 7\cdot 11\cdot e^4-2^{10}{e}^2-2^7\cdot e^4$
		$+5\cdot 2^9(-e^2+e^3-e^4)+5\cdot 2^6(e^2-2e^3+3e^4)+5\cdot 2^6(e^3-3e^4)+5\cdot 77e^4$
		$+5\cdot 2^9(e^2-e^3+e^4)+5\cdot 2^6(e^3-2e^4)+5\cdot 2^6{e}^4+{𝒪}(e^5)$
	$=$	$e^2[2^9-2^6\cdot 11+2^9+3\cdot 2^7-2^{10}-5\cdot 2^9+5\cdot 2^6+5\cdot 2^9]$
		$+e^3[-2^9-2^7\cdot 11-2^6\cdot 17-2^9+2^6\cdot 11-3\cdot 2^8+3\cdot 2^8+5\cdot 2^9-5\cdot 2^7+5\cdot 2^6-5\cdot 2^9+5\cdot 2^6]$
		$+e^4[2^9-3\cdot 11\cdot 2^6+3\cdot 17\cdot 2^6-1787+2^9-11\cdot 2^7+17\cdot 2^6+3^2\cdot 2^7-3^2\cdot 2^8+2\cdot 3^2\cdot 7\cdot 11-2^7-5\cdot 2^9+3\cdot 5\cdot 2^6-3\cdot 5\cdot 2^6+5\cdot 7\cdot 11+5\cdot 2^9-5\cdot 2^7+5\cdot 2^6]+{𝒪}(e^5)$
	$=$	$e^2[-2^6\cdot 11+3\cdot 2^7+5\cdot 2^6]+e^3[-2^{10}-2^7\cdot 11-2^6\cdot 17+2^6\cdot 11]$
		$+e^4[2^{10}-1787+2\cdot 3^2\cdot 7\cdot 11+5\cdot 7\cdot 11+2^6\cdot (3\cdot 17-3\cdot 11-22+17+18-36-2-10+5)]+{𝒪}(e^5)$
	$=$	$2^6{e}^2[0]-2^8\cdot 11e^3+e^4[2^{10}-1787+7\cdot 11\cdot 23-2^8\cdot 3]+{𝒪}(e^5)$
	$=$	$-2^8\cdot 11e^3+2^4\cdot 3\cdot 5e^4+{𝒪}(e^5)$
$\Rightarrow \frac{2^2}{{\pi }^2}\left[\frac{e^2}{(1-e^2{)}^{1/2}}\right]{\mathrm{{\rm Y}}}_{W1}({\eta }_0)$	$=$	$-\frac{11}{2}{e}^3+\frac{3\cdot 5}{2^5}{e}^4+{𝒪}(e^5).$

Floating Comparison Summary

As shown above, the first three terms of the Huré, et al. (2020) series expression may be written as, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Psi }}_0+{\mathrm{\Psi }}_1+{\mathrm{\Psi }}_2}{GM}\right]$ $=$ $\mathit{K}([k_H{]}_0)-\frac{e^2}{2^4}\cdot \mathit{K}(k_H)+\frac{e^2}{2^4}[\frac{{\mathrm{\Delta }}_0^2-2R_c(R_c+R)}{(1-k_H^2){\mathrm{\Delta }}_0^2}]\mathit{E}(k_H).$ Let's see how it compares to the first term of Wong's (1973) expression which, as shown separately above, can be written in the form, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$ $=$ ${\mathrm{\Delta }}_0\left(\frac{2^2}{3{\pi }^2}\right){\mathrm{{\rm Y}}}_{W0}({\eta }_0)\left\{\frac{\mathit{K}(k)}{r_1}\right\}.$ First, as shown above, $\frac{2^2}{3{\pi }^2}{\mathrm{{\rm Y}}}_{W0}({\eta }_0)$ $=$ $[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot (1-e^2{)}^{1/2}.$ Note that, in order to determine the functional form of the ${𝒪}(e^2)$ term in this expression, we will have to include $k_0^8$ terms in the various expressions for products of elliptic integrals. Second, we have shown that, $\frac{r_1^2}{{\mathrm{\Delta }}_0^2}$ $=$ $1-e^2\left[\frac{R_c(R_c+{\varpi }_W)}{{\mathrm{\Delta }}_0^2}\right]-\frac{{\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^2}[\frac{1}{2^2}{e}^4+\frac{1}{2^3}{e}^6+\frac{5}{2^6}{e}^8+{𝒪}(e^{10})]$ $\Rightarrow \frac{{\mathrm{\Delta }}_0}{r_1}$ $\approx$ $1+e^2\left[\frac{R_c(R_c+{\varpi }_W)}{2{\mathrm{\Delta }}_0^2}\right],$       and we are defining ${\delta }_K$ such that, $K(k)$ $=$ $K(k_H)+\frac{\pi }{2}\cdot {\delta }_K.$ Hence, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{0}}}{GM}\right]$ $\approx$ $\{K(k_H)+\frac{\pi }{2}\cdot {\delta }_K\}\{1+e^2[\frac{R_c(R_c+{\varpi }_W)}{2{\mathrm{\Delta }}_0^2}\left]\right\}[1+(\frac{5\cdot 13}{2^4\cdot 3})e^2+{𝒪}(e^3)]\cdot (1-e^2{)}^{1/2}$   $\approx$ $\mathit{K}([k_H{]}_0)+\mathit{K}(k_H)e^2\cdot {𝒜}+\frac{\pi }{2}[\frac{e^2\cdot {\varpi }_W{R}_c}{{\mathrm{\Delta }}_0^4}\left]\right\{-R_c^2-\cos {\theta }_H[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H]$     $+({\varpi }_W{R}_c+R_c^2)-\frac{1}{2}[({\varpi }_W^2+R_c{)}^2+z_W^2]-R_c^2-2R_c\cos {\theta }_H(R_c\cos {\theta }_H+{\varpi }_W\cos {\theta }_H-z_W\sin {\theta }_H)$     $+\frac{[2{\varpi }_W{R}_c\cos {\theta }_H+2R_c^2\cos {\theta }_H-2z_W{R}_c\sin {\theta }_H{]}^2}{{\mathrm{\Delta }}_0^2}\},$ and, ${𝒜}$ $\approx$ $\left[\frac{R_c(R_c+{\varpi }_W)}{2{\mathrm{\Delta }}_0^2}\right]+\left(\frac{41}{2^4\cdot 3}\right).$ Second, $-\frac{\pi {\mathrm{\Delta }}_0}{2}\left[\frac{{\mathrm{\Phi }}_{\mathrm{W}\mathrm{1}}}{GM}\right]$ $=$ $\mathit{E}(k)\left\{\frac{{\mathrm{\Delta }}_0\cdot \cos \theta }{r_2}\right\}\left[\right(\frac{2^2}{3{\pi }^2}){\mathrm{{\rm Y}}}_{W1}({\eta }_0)]$   $=$ $\{\mathit{E}(k_H)+\frac{\pi }{2}\cdot {\delta }_E\}\left\{\frac{{\mathrm{\Delta }}_0\cdot \cos \theta }{r_2}\right\}[-(\frac{11}{2\cdot 3})e+(\frac{5}{2^5})e^2+{𝒪}(e^5)]\cdot (1-e^2{)}^{1/2}.$

Geometric Factor

By definition,

${\mathrm{\Delta }}_0^2$	$=$	$({\varpi }_W+R_c{)}^2+z_W^2,$
$r_1^2$	$=$	$[{\varpi }_W+R_c(1-e^2{)}^{1/2}{]}^2+z_W^2,$
$r_2^2$	$=$	$[{\varpi }_W-R_c(1-e^2{)}^{1/2}{]}^2+z_W^2,$
${\cos }^2\theta$	$=$	$[\frac{r_1^2+r_2^2-4R_c^2(1-e^2)}{2r_1{r}_2}{]}^2.$

Hence,

$r_2^2-{\mathrm{\Delta }}_0^2\cdot {\cos }^2\theta$

$=$

$$

See Also

• Université de Bordeaux (Part 2):  Spheroid-Ring Sequences
• Université de Bordeaux (Part 3):  Discussions Following Dissertation Defense