Apps/DysonPotential: Difference between revisions

From jetwiki
Jump to navigation Jump to search
No edit summary
Line 27: Line 27:
<table border="1" cellpadding="5" align="center">
<table border="1" cellpadding="5" align="center">
<tr><td align="center" colspan="1">
<tr><td align="center" colspan="1">
'''Equation image extracted without modification from p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D Dyson (1893a)]'''<br />
'''Equation &amp; text copied without modification from <br />p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D Dyson (1893a)]'''<br />
''The Potential of an Anchor Ring'', Phil. Trans. Royal Soc. London. A., Vol. 184
''The Potential of an Anchor Ring''<br /> Phil. Trans. Royal Soc. London. A., Vol. 184
</td></tr>
</td></tr>
<tr>
<tr>
<td>
<td align="left">
[[File:DysonExternalPotentialEquation.png|550px|center|To be inserted:  the Potential Exterior to an Anchor Ring]]
<!-- [[File:DysonExternalPotentialEquation.png|550px|center|To be inserted:  the Potential Exterior to an Anchor Ring]] -->
Therefore, at any external point,
<table border="0" align="center" width="100%" cellpadding="5">
<tr>
  <td align="right"><math>V</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{4F}{R + R_1} \biggl\{1 - \frac{1}{8} \frac{a^2}{c^2} \cos^2 \frac{\psi}{2}
-
\frac{1}{768} \frac{a^4}{c^4} \biggl[ 5 + 8\cos\psi - \cos^2\psi - 4\cos^3\psi - \frac{4c^2}{RR_1} \cos 2\psi \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">&nbsp;</td>
  <td align="center">&nbsp;</td>
  <td align="left">
&nbsp; &nbsp; &nbsp;&nbsp; &nbsp; &nbsp;<math>+</math>&nbsp; &amp;c.<math>\biggr\}</math>
  </td>
</tr>
 
<tr>
  <td align="right">&nbsp;</td>
  <td align="center">&nbsp;</td>
  <td align="left">
<math>+~\frac{E(R+R_1)}{RR_1} \biggl\{\frac{a^2}{8c^2} \cos\psi
-
\frac{1}{192} \frac{a^4}{c^4} \biggl[ 2\cos^2\psi - 4\cos\psi + \frac{2c^2}{RR_1} \cos 2\psi \biggr] + </math>
&nbsp; &amp;c.<math>\biggr\} \, ,</math>
  </td>
</tr>
</table>
where <math>\psi</math> is the angle between <math>R</math> and <math>R_1</math>, and the modulus of the elliptic functions is
<table border="0" align="center" width="100%" cellpadding="5">
<tr>
  <td align="center"><math>\frac{R_1 - R}{R_1 + R} \, .</math></td>
</tr>
</table>
</td>
</td>
</tr>
</tr>

Revision as of 13:21, 15 September 2021

Dyson (1893)

Dyson
(1893)

Our focus, here, is on the pioneering work of F. W. Dyson (1893a, Philosophical Transactions of the Royal Society of London. A., 184, 43 - 95) and (1893b, Philosophical Transactions of the Royal Society of London. A., 184, 1041 - 1106). He used analytic techniques to determine the approximate equilibrium structure of axisymmetric, uniformly rotating, incompressible tori. C.-Y. Wong (1974, ApJ, 190, 675 - 694) extended Dyson's work, using numerical techniques to obtain more accurate equilibrium structures for incompressible tori having solid body rotation. Since then, Y. Eriguchi & D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870 - 1875) and I. Hachisu, J. E. Tohline & Y. Eriguchi (1987, ApJ, 323, 592 - 613) have mapped out the full sequence of Dyson-Wong tori, beginning from a bifurcation point on the Maclaurin spheroid sequence.

The most challenging aspect of each of these studies has been the development of an analytic and/or computational technique that can be used to accurately determine the gravitational potential of toroidal-shaped configurations. With this in mind, it should be appreciated that, in a paper that preceded his 1974 work, C.-Y. Wong (1973, Annals of Physics, 77, 279) derived an analytic expression for the exact potential (inside as well as outside) of axisymmetric, uniform-density tori having an arbitrarily specified ratio of the major to minor (cross-sectional) radii, R/d. This is an outstanding accomplishment that has received little attention in the astrophysics literature and, therefore, has heretofore been under-appreciated. In a separate, accompanying discussion, we detail how Wong accomplished this task.

External Potential

His Derived Expression

(See an accompanying Ramblings Chapter for additional derivation details.) On p. 62, in §8 of Dyson (1893a), we find the following approximate expression for the potential at point "P", anywhere exterior to an anchor ring:

Anchor Ring Schematic
Anchor Ring Schematic

Caption: Anchor ring schematic, adapted from figure near the top of §2 (on p. 47) of Dyson (1893a)

Equation & text copied without modification from
p. 62 of Dyson (1893a)

The Potential of an Anchor Ring
Phil. Trans. Royal Soc. London. A., Vol. 184

Therefore, at any external point,

V =

4FR+R1{118a2c2cos2ψ21768a4c4[5+8cosψcos2ψ4cos3ψ4c2RR1cos2ψ]

   

          +  &c.}

   

+E(R+R1)RR1{a28c2cosψ1192a4c4[2cos2ψ4cosψ+2c2RR1cos2ψ]+   &c.},

where ψ is the angle between R and R1, and the modulus of the elliptic functions is

R1RR1+R.

In Dyson's expression, the leading factor of F is the complete elliptic integral of the first kind, namely,

F=F(μ)

0π/2dϕ1μ2sin2ϕ,

where, μ(R1R)/(R1+R). Similarly, E=E(μ) is the complete elliptic integral of the second kind.

Comparison With Thin Ring Approximation

In the limit of a/c0, Dyson's expression gives,

VDyson

=

4K(μ)R+R1,

where we have acknowledged that, in the twenty-first century, the complete elliptic integral of the first kind is more customarily represented by the letter, K. In a separate discussion, we have shown that the gravitational potential of an infinitesimally thin ring is given precisely by the expression,

[πGM]ΦTR

=

2K(k)R1,

where, k[1(R/R1)2]1/2. Is Dyson's expression identical to this one when a/c=0 ?

Proof

Taking a queue from our accompanying discussion of toroidal coordinates, if we adopt the variable notation,

ηln(R1R),

then we can write,

coshη=12[eη+eη]

=

R2+R122RR1,

which implies that,

[2cothη+1]1/2=[1e2η]1/2

=

[1(RR1)2]1/2.

This is the definition of the parameter, k, in the expression for ΦTR. Now, if we employ the Descending Landen Transformation for the complete elliptic integral of the first kind, we can make the substitution,

K(k)

=

(1+k1)K(k1),

      where,      

k1

11k21+1k2.

But notice that, 1k2=eη, in which case,

k1

=

1eη1+eη

=

1R/R11+R/R1

=

R1RR1+R,

which is the definition of the parameter, μ, in the expression for VDyson. Hence, we can write,

[πGM]ΦTR

=

2R1[(1+k1)K(k1)]

 

=

2K(μ)R1[1+R1RR1+R]

 

=

4K(μ)R1+R.

Aside from the adopted sign convention, this is indeed precisely the expression given by VDyson when a/c=0 .

Evaluation

Dyson's Figures

In his effort to illustrate the behavior of equipotential contours in the space exterior to various anchor rings, Dyson evaluated his expression for the potential up through 𝒪(a2c2); that is, he evaluated the function,

V2VDyson|𝒪(a2/c2)

=

4K(μ)R+R1[118(a2c2)cos2(ψ2)]+(R+R1)E(μ)RR1[18(a2c2)cosψ].

Figures 1 - 6 from Dyson (1893a) — replicated immediately below — show his resulting set of contours for six cases: Tori (anchor rings) having aspect ratios of a/c=0,1/5,2/5,3/5,4/5,1. Click on an image to view the contour plot at higher resolution. In what follows we present results from our own evaluation of this "V2" function for the single case of an anchor ring having a/c=2/5.


Figures 1 - 6 extracted without modification from pp. 63-66 of F. W. Dyson (1893)

The Potential of an Anchor Ring, Phil. Trans. Royal Soc. London. A., Vol. 184

The Potential Exterior to an Anchor Ring; R/d = infinity
The Potential Exterior to an Anchor Ring; R/d = infinity
The Potential Exterior to an Anchor Ring; R/d = 5
The Potential Exterior to an Anchor Ring; R/d = 5
The Potential Exterior to an Anchor Ring; R/d = 2.5
The Potential Exterior to an Anchor Ring; R/d = 2.5
The Potential Exterior to an Anchor Ring; R/d = 1.667
The Potential Exterior to an Anchor Ring; R/d = 1.667
The Potential Exterior to an Anchor Ring; R/d = 1.25
The Potential Exterior to an Anchor Ring; R/d = 1.25
The Potential Exterior to an Anchor Ring; R/d = 1
The Potential Exterior to an Anchor Ring; R/d = 1

Our Attempt to Replicate

First, let's test the accuracy of Dyson's (1893a) "series expansion" expression for the elliptic integrals, K(μ) and E(μ); in the following table, the high-precision evaluations labeled "Numerical Recipes" have been drawn from the tabulated data that is provided in our accompanying discussion of incomplete elliptic integrals. Drawing from our accompanying set of Key mathematical relations — in which k, rather than μ, represents the function modulus — the relevant series-expansion expressions are:

2K(k)π

=

1+(12)2k2+(1324)2k4+(135243)2k6+(1357273)2k8++[(2n1)!!2nn!]2k2n+

Gradshteyn & Ryzhik (1965), §8.113.1


2E(k)π

=

1122k21232242k4(135243)2k65(1357273)2k87[(2n1)!!2nn!]2k2n2n1

Gradshteyn & Ryzhik (1965), §8.114.1

These expressions — up through 𝒪(μ4) — can be found in the middle of p. 58 of Dyson (1893a). We strongly suspect that, in constructing the equipotential contours shown in his figures 1-6, Dyson used expressions for K(μ) and E(μ) that were more accurate than this. For example, we found it necessary to include terms up through 𝒪(μ10) in order to match to three digits accuracy the potential contour values and coordinate locations reported by Dyson.

μ Numerical Recipes Series expansion up through 𝒪(μ4) Series expansion up through 𝒪(μ10)
K(μ) E(μ) K(μ) E(μ) K(μ) E(μ)
0.34202014 1.62002589 1.52379921 1.6198 1.5239 1.6200263 1.5237989
0.57357644 1.73124518 1.43229097 1.7239 1.4336 1.73124518 1.43230
0.76604444 1.93558110 1.30553909 1.8773 1.3150 1.93558109 1.3061
0.90630779 2.30878680 1.16382796 2.042 1.199 2.308784 1.1700
0.98480775 3.15338525 1.04011440 2.16 1.12 3.150 1.069

We actually used the "descending Landen transformation" to evaluate K(μ) through 𝒪(μ10).


For c=1 and a specification of the ratio, a/c, take the following steps to map out an equipotential curve that has V2=V0:

  • Choose a value of Ra
    • Guess a value of (cR)R1(c+R)ϖ=(R12R2)/(4c)     and,     z=±R12(c+ϖ)2
    • Set cosψ=(R12+R24c2)/(2RR1)
    • Evaluate the function, V2
    • If V2V0 to the desired accuracy, loop back up and guess another value of R1
  • If V2=V0 to the desired accuracy, save the coordinate location, (ϖ,z), and loop back up to pick another value of R


The Potential Exterior to an Anchor Ring; R/d = 2.5
The Potential Exterior to an Anchor Ring; R/d = 2.5
   Compare with Dyson

Tabulated Data

As the data in the following table documents, we have been able to construct equipotential contours that agree with Dyson, not only qualitatively, but quantitatively. For example:

  • The dark green contour has been designed to touch the surface of the torus precisely where its outermost edge cuts through the equatorial plane (ϖ,z)=(1.4,0). This means that R=0.4 and R1=2.4. (These four coordinate values are highlighted in pink in the second major column of the table.) When we plugged these values of R and R1 into Dyson's expression for V2, we determined that the value of the potential at this point on the torus surface is 0.8551 — see the yellow-highlighted heading of the second major table column. Compare this to the value of 0.855 that Dyson has printed just below the Figure 3 x-axis where a fiducial identifies the coordinate, ϖ=1.4. As has been catalogued at the bottom of table column #2, we have found that this dark-green contour touches the vertical axis at the coordinate location, (ϖ,z)=(0,0.572), for which, R1=R=1.1518.
  • By design — see the coordinate values highlighted in pink in table column #1 — our outermost (pink) contour touches the equatorial plane at (ϖ,z)=(1.5,0)(R,R1)=(0.5,2.5). When we plugged these values of R and R1 into Dyson's expression for V2, we determined that the value of the potential at this point outside the torus is 0.7737 — see the yellow-highlighted heading of table column #1. Compare this to the value of 0.777 that Dyson has printed just below the Figure 3 x-axis where a fiducial identifies the coordinate, ϖ=1.5. As has been catalogued at the bottom of table column #1, we have found that this pink contour touches the vertical axis at the coordinate location, (ϖ,z)=(0,0.794), for which, R1=R=1.2766.
  • Similarly, we have constructed contours that intersect the equatorial plane at the fiducials marking ϖ=0.0 (red curve & table column #5), ϖ=0.2 (light-green curve & table column #6), and ϖ=0.4 (light-blue curve & table column #7). According to our calculations, they correspond, respectively, to values of the potential, V2=0.9800 (Dyson's corresponding fiducial label is 0.980), V2=0.9896 (Dyson's corresponding fiducial label is 0.990), and V2=1.0212 (Dyson's corresponding fiducial label is 1.021).
  • Finally, we constructed two contours (blue and orange) by initially specifying the value of the potential, rather than specifying the coordinate values (R,R1). We used the values of the potential that Dyson associated with the fiducials along the vertical axis at (ϖ,z)=(0.0,0.4) and at (ϖ,z)=(0.0,0.2): Respectively, V2=0.912 — blue contour detailed in our table column #3 — and V2=0.961— orange contour detailed in our table column #4. We determined that these two contour curves intersected the vertical axis at, respectively, (ϖ,z)=(0.0,0.402) and (ϖ,z)=(0.0,0.204), that is, at coordinate locations that were nearly identical to the locations labeled by Dyson.


Coordinates of Points that Trace Seven Different Equipotential Contours External to the Anchor Ring With c/a=5/2

Column #1 Column #2 Column #3 Column #4 Column #5 Column #6 Column #7
V2 = 0.7737

R

R1

ϖ

z

0.5000

2.5000

1.500

0.000

0.5005

2.4990

1.499

0.043

0.504

2.4889

1.485

0.137

0.510

2.4720

1.463

0.215

0.520

2.4445

1.426

0.298

0.530

2.4177

1.391

0.358

0.550

2.3665

1.324

0.444

0.580

2.2940

1.232

0.532

0.610

2.2265

1.146

0.592

0.640

2.1632

1.067

0.636

0.700

2.0465

0.925

0.696

0.800

1.8745

0.718

0.749

0.9000

1.7240

0.541

0.774

1.000

1.5890

0.381

0.786

1.100

1.4670

0.236

0.791

1.2000

1.3558

0.100

0.793

1.277

1.2766

0.000

0.794



V2 = 0.8551

R

R1

ϖ

z

0.400

2.4000

1.400

0.000

0.405

2.3830

1.379

0.144

0.410

2.3668

1.358

0.199

0.425

2.3190

1.299

0.302

0.450

2.2458

1.210

0.398

0.480

2.1655

1.115

0.466

0.520

2.0690

1.003

0.520

0.570

1.9610

0.880

0.557

0.620

1.8635

0.772

0.577

0.700

1.7240

0.621

0.588

0.800

1.5712

0.457

0.588

0.900

1.4360

0.313

0.581

1.000

1.3147

0.182

0.575

1.100

1.2050

0.061

0.572

1.1518

1.1518

0.000

0.572


V2 = 0.9120

R

R1

ϖ

z

1.0776

1.0776

0.000

0.402

1.000

1.1582

0.085

0.404

0.950

1.2135

0.143

0.409

0.900

1.2715

0.202

0.416

0.800

1.3979

0.328

0.435

0.700

1.5401

0.470

0.458

0.600

1.7040

0.636

0.477

0.550

1.7970

0.732

0.480

0.500

1.8998

0.840

0.474

0.475

1.9560

0.900

0.464

0.440

2.0410

0.993

0.440

0.400

2.1510

1.117

0.383


V2 = 0.9610

R

R1

ϖ

z

1.0206

1.0206

0.000

0.204

0.9500

1.0937

0.073

0.210

0.900

1.1488

0.127

0.221

0.800

1.2685

0.242

0.257

0.700

1.4030

0.370

0.304

0.600

1.5572

0.516

0.355

0.550

1.6440

0.600

0.378

0.500

1.7395

0.694

0.395

0.450

1.8462

0.801

0.404

0.410

1.9690

0.929

0.394



V2 = 0.9800

R

R1

ϖ

z

1.0000

1.0000

0.000

0.000

0.900

1.1053

0.103

0.072

0.800

1.2225

0.214

0.147

0.700

1.3543

0.336

0.222

0.600

1.5050

0.476

0.293

0.550

1.5897

0.556

0.325

0.500

1.6827

0.645

0.352

0.450

1.7865

0.747

0.372

0.400

1.9050

0.867

0.377


V2 = 0.9896

R

R1

ϖ

z

0.8000

1.2000

0.200

0.000

0.7950

1.2062

0.206

0.034

0.780

1.2248

0.223

0.068

0.760

1.2503

0.246

0.099

0.730

1.2895

0.282

0.134

0.700

1.3305

0.320

0.166

0.650

1.4022

0.386

0.213

0.600

1.4796

0.457

0.256

0.550

1.5633

0.535

0.294

0.500

1.6552

0.622

0.328

0.450

1.7573

0.721

0.353

0.400

1.8737

0.838

0.366



V2 = 1.0212

R

R1

ϖ

z

0.6000

1.4000

0.400

0.000

0.5950

1.4078

0.407

0.048

0.580

1.4315

0.428

0.097

0.570

1.4477

0.443

0.120

0.540

1.4978

0.488

0.171

0.500

1.5688

0.553

0.224

0.450

1.6663

0.644

0.275

0.400

1.7767

0.749

0.312

Intermediate Step

Objective

As has been reprinted above, on p. 62 of Dyson's Part I we find his power-series expression for the external potential, namely,

πVDysonGM|𝒪(a4/c4)

=

4K(μ)R+R1{118(a2c2)cos2(ψ2)1768(ac)4[5+8cosψcos2ψ4cos3ψ4c2RR1cos2ψ]}

 

 

+(R+R1)E(μ)RR1{18(ac)2cosψ1192(ac)4[2cos2ψ4cosψ+2c2RR1cos2ψ]},

where — as in the context of toroidal coordinates — we occasionally will make the substitution, eη=R1/R, and therefore,

μ

R1RR1+R=eη1eη+1.

In order to facilitate matching boundary conditions at the surface of the torus, between the exterior and interior expressions for the gravitational potential, Dyson rewrites this Part I expression for the external potential and — explicitly evaluating it on the torus surface — sets, R=a. Specifically, on p. 1049 of Dyson's Part II we find equation (6), which reads,

V2πa2

=

ln(8ca)+12(ac)[ln(8ca)54]cosχ+{116[ln(8ca)52]+316[ln(8ca)+17367236]cos2χ}(a2c2)

 

 

+{332[ln(8ca)2512]cosχ+564[ln(8ca)+7244824]cos3χ}(a3c3)

 

 

+{9256[ln(8ca)2]+7128[ln(8ca)191682]cos2χ+351024[ln(8ca)2+19120]cos4χ}(a4c4)+

In order to obtain this alternate power-series expression, Dyson …

  • Expresses angular variations in terms of the angle, χ, instead of the angle, ψ; these two angles are identified in the above schematic.
  • Employs power-series expansions of both elliptic integral functions, K(μ) and E(μ).
  • Uses the binomial theorem to develop a number of other power-series expressions.

In what follows we will attempt to demonstrate that this second (Part II, equation 6) expression is identical to the first.

The Ratio R1/c

Note that, via the law of cosines,

R12

=

(2c)2+R24Rccosχ

(R1c)2

=

4+(Rc)24(Rc)cosχ

At the surface of the torus, where R=a, we therefore have,

R1c

=

2[1(ac)cosχ+14(ac)2]1/2.

Low Order

Employing the binomial theorem, we can write,

(R1c)1

=

12[1(ac)cosχ+14(ac)2]1/2

 

12{112[(ac)cosχ+14(ac)2]+38[(ac)cosχ+14(ac)20]2}

 

12{1+12(ac)cosχ18(ac)2+38(ac)2cos2χ}

 

12{1+12(ac)cosχ+(ac)2[38cos2χ18]}

1+ac(R1c)1

1+12(ac)+14(ac)2cosχ

[1+ac(R1c)1]1

1[12(ac)+14(ac)2cosχ]+[12(ac)+14(ac)2cosχ]2

 

112(ac)+14(ac)2(1cosχ)

Higher Order

Adopting the shorthand notation,

γ12(R1c),      and       b(ac)cosχ+14(ac)2,

and employing the binomial theorem, we can write,

γ=[1+b]1/2

=

1+12b123b2+124b335273b4+𝒪(a5c5)

 

=

1+12[(ac)cosχ+14(ac)2]123[(ac)cosχ+14(ac)2]2

 

 

+124[(ac)cosχ+14(ac)2]335273[(ac)cosχ+14(ac)2]4+𝒪(a5c5)

 

=

112(ac)cosχ+123(ac)2123[(ac)2cos2χ12(ac)3cosχ+124(ac)4]

 

 

+124[(ac)cosχ+14(ac)2][(ac)2cos2χ12(ac)3cosχ]35273[(ac)4cos4χ]+𝒪(a5c5)

 

=

112(ac)cosχ+123(ac)2123(ac)2cos2χ+124(ac)3cosχ127(ac)4

 

 

124(ac)3cos3χ+125(ac)4cos2χ+126(ac)4cos2χ35273(ac)4cos4χ+𝒪(a5c5)

 

=

112(ac)cosχ+123(ac)2(1cos2χ)+124(ac)3(cosχcos3χ)+127(ac)4[1+6cos2χ5cos4χ]+𝒪(a5c5).

Also, we have,

1γ=2(R1c)1=[1+b]1/2

=

112b+323b235243b3+357273b4+𝒪(a5c5)

 

=

112[(ac)cosχ+14(ac)2]+323[(ac)cosχ+14(ac)2]2

 

 

35243[(ac)cosχ+14(ac)2]3+357273[(ac)cosχ+14(ac)2]4+𝒪(a5c5)

 

=

1+12(ac)cosχ123(ac)2+323[(ac)2cos2χ12(ac)3cosχ+124(ac)4]

 

 

35243[(ac)cosχ+14(ac)2][(ac)2cos2χ12(ac)3cosχ]+357273[(ac)4cos4χ]+𝒪(a5c5)

 

=

1+12(ac)cosχ123(ac)2+323(ac)2cos2χ324(ac)3cosχ+327(ac)4

 

 

+35243[(ac)3cos3χ12(ac)4cos2χ]35263[(ac)4cos2χ]+357273[(ac)4cos4χ]+𝒪(a5c5)

 

=

1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]

 

 

+127(ac)4[330cos2χ+35cos4χ]+𝒪(a5c5).

Hence,

1+(ac)(R1c)1

=

1+12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]+125(ac)4[5cos3χ3cosχ]+𝒪(a5c5).

And, adopting the shorthand notation,

d12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]+125(ac)4[5cos3χ3cosχ],

we have,

[1+(ac)(R1c)1]1

=

1d+d2d3+d4+𝒪(a5c5)

 

=

112(ac)122(ac)2cosχ124(ac)3[3cos2χ1]125(ac)4[5cos3χ3cosχ]

 

 

+{12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]}{12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]}

 

 

[12(ac)+122(ac)2cosχ][12(ac)+122(ac)2cosχ]2+124(ac)4+𝒪(a5c5)

 

=

112(ac)122(ac)2cosχ124(ac)3[3cos2χ1]125(ac)4[5cos3χ3cosχ]

 

 

+12(ac){12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]}

 

 

+{122(ac)2cosχ}{12(ac)+122(ac)2cosχ}+124(ac)3[3cos2χ1]{12(ac)}

 

 

[12(ac)+122(ac)2cosχ][12(ac)+122(ac)2cosχ]2+124(ac)4+𝒪(a5c5)

 

=

112(ac)122(ac)2cosχ124(ac)3[3cos2χ1]125(ac)4[5cos3χ3cosχ]

 

 

+122(ac)2+123(ac)3cosχ+125(ac)4[3cos2χ1]+123(ac)3cosχ+124(ac)4cos2χ+125(ac)4[3cos2χ1]

 

 

[122(ac)2+122(ac)3cosχ+124(ac)4cos2χ][12(ac)+122(ac)2cosχ]+124(ac)4+𝒪(a5c5)

 

=

112(ac)+122(ac)2(1cosχ)+124(ac)3[2cosχ+2(cosχ1)2(3cos2χ1)]

 

 

+125(ac)4[(3cos2χ1)+2cos2χ+(3cos2χ1)2cosχ4cosχ+2(5cos3χ3cosχ)]+𝒪(a5c5)

 

=

112(ac)+122(ac)2(1cosχ)+123(ac)3(2cosχ3cos2χ)

 

 

+125(ac)4(9cosχ+8cos2χ5cos3χ)+𝒪(a5c5).

Relationship Between Angles

Drawing on the Law of Cosines, as above, we can state that on the torus surface,

R12

=

(2c)2+a24accosχ

Alternatively, applying the Law of Cosines to the angle, ψ, we have,

(2c)2

=

R12+a22aR1cosψ

cosψ

=

R12+a24c22aR1.

Therefore, anywhere along the surface of the torus, we can switch from one of these angles to the other via the relation,

cosψ

=

12a[4c2+a24accosχ+a24c2][4c2+a24accosχ]1/2

 

=

[cosχ12(ac)][1(ac)+14(ac)2cosχ]1/2.

Cosine ψ Expansion

Employing the binomial theorem, we therefore can write,

cosψ

=

[cosχ12(ac)]{112[(ac)+14(ac)2cosχ]+38[(ac)+14(ac)2cosχ]2

 

 

524[(ac)]3+5727[(ac)]4+𝒪(a5c5)}

 

=

[cosχ12(ac)]{1+12(ac)123(ac)2cosχ+38[(ac)214(ac)3cosχ+124(ac)4cos2χ]

 

 

+524(ac)3+5727(ac)4+𝒪(a5c5)}

 

=

[cosχ12(ac)]{1+12(ac)+123(ac)2[3cosχ]+125(ac)3[103cosχ]+127(ac)4[3cos2χ+57]+𝒪(a5c5)}

cosψ|𝒪(a2/c2)

=

[cosχ12(ac)]{1+12(ac)+123(ac)2[3cosχ]}

 

=

cosχ{1+12(ac)+123(ac)2[3cosχ]}+12(ac){1+12(ac)}

 

=

cosχ+12(ac)(1cosχ)+123(ac)2[23cosχ+cos2χ].

Cosine-Squared Expansion

Letting,

b[(ac)+14(ac)2cosχ],

via the binomial theorem we have,

cos2ψ

=

[cosχ12(ac)]2[1(ac)+14(ac)2cosχ]1

 

=

[cosχ12(ac)]2{1b+b2b3+b4𝒪(b5)}

 

=

[cosχ12(ac)]2{1[(ac)+14(ac)2cosχ]+[(ac)+14(ac)2cosχ]2

 

 

[(ac)+14(ac)2cosχ]3+[(ac)+14(ac)2cosχ]4+𝒪(a5c5)}

 

=

[cosχ12(ac)]2{1+[(ac)14(ac)2cosχ]+[(ac)212(ac)3cosχ+124(ac)4cos2χ]

 

 

+[(ac)14(ac)2cosχ][(ac)212(ac)3cosχ+124(ac)4cos2χ]

 

 

+[(ac)212(ac)3cosχ+124(ac)4cos2χ]2+𝒪(a5c5)}

 

=

[cosχ12(ac)]2{1+(ac)14(ac)2cosχ+(ac)212(ac)3cosχ+124(ac)4cos2χ

 

 

+(ac)312(ac)4cosχ14(ac)4cosχ+(ac)4+𝒪(a5c5)}

 

=

[cos2χ(ac)cosχ+122(ac)2]{1+(ac)+(ac)2[114cosχ]+(ac)3[112cosχ]

 

 

+(ac)4[1+124cos2χ34cosχ]+𝒪(a5c5)}

 

=

cos2χ{1+(ac)+(ac)2[114cosχ]+(ac)3[112cosχ]+(ac)4[1+124cos2χ34cosχ]}

 

 

(ac)cosχ{1+(ac)+(ac)2[114cosχ]+(ac)3[112cosχ]}

 

 

+122(ac)2{1+(ac)+(ac)2[114cosχ]}+𝒪(a5c5)

 

=

cos2χ+(ac)[cos2χcosχ]+(ac)2[cos2χ14cos3χcosχ+122]

 

 

+(ac)3{cos2χ[112cosχ]cosχ[114cosχ]+122}

 

 

+(ac)4{cos2χ[1+124cos2χ34cosχ]cosχ[112cosχ]+122[114cosχ]}+𝒪(a5c5)

 

=

cos2χ+(ac)[cos2χcosχ]+(ac)2[122cosχ+cos2χ14cos3χ]

 

 

+(ac)3[122cosχ+54cos2χ12cos3χ]+(ac)4[1221724cosχ+32cos2χ34cos3χ+124cos4χ]+𝒪(a5c5)

Cosine-Cubed Expansion

Again, letting,

b[(ac)+14(ac)2cosχ],

via the binomial theorem we have,

cos3ψ

=

[12(ac)cosχ]3[1(ac)+14(ac)2cosχ]3/2

 

=

[12(ac)cosχ]3{132[b]+3523[b]2357243[b]3+3579273[b]4+𝒪(b5)}

 

=

[12(ac)cosχ]3{132[(ac)+14(ac)2cosχ]+3523[(ac)+14(ac)2cosχ]2

 

 

5724[(ac)+14(ac)2cosχ]3+57927[(ac)+14(ac)2cosχ]4+𝒪(b5)}

 

=

[12(ac)cosχ]3{1+32(ac)323(ac)2cosχ+3523[(ac)212(ac)3cosχ+124(ac)4cos2χ]

 

 

5724[(ac)+14(ac)2cosχ][(ac)212(ac)3cosχ+124(ac)4cos2χ]

 

 

+57927[(ac)212(ac)3cosχ+124(ac)4cos2χ]2+𝒪(b5)}

 

=

[12(ac)cosχ]3{1+32(ac)323(ac)2cosχ+3523[(ac)212(ac)3cosχ+124(ac)4cos2χ]

 

 

+5724(ac)[(ac)212(ac)3cosχ]5726(ac)4cosχ+325727(ac)4+𝒪(a5c5)}

 

=

[cos2χ(ac)cosχ+122(ac)2][12(ac)cosχ]{1+32(ac)323(ac)2cosχ+3523(ac)2

 

 

3524(ac)3cosχ+3527(ac)4cos2χ+5724(ac)35725(ac)4cosχ5726(ac)4cosχ+325727(ac)4+𝒪(a5c5)}

 

=

{cos3χ+32(ac)cos2χ322(ac)2cosχ+123(ac)3}

 

 

×{1+32(ac)+323(ac)2[5cosχ]+524(ac)3[73cosχ]+527(ac)4[3cos2χ227cosχ27cosχ+327]+𝒪(a5c5)}

 

=

cos3χ{1+32(ac)+323(ac)2[5cosχ]+524(ac)3[73cosχ]+527(ac)4[3cos2χ227cosχ27cosχ+327]}

 

 

+32(ac)cos2χ{1+32(ac)+323(ac)2[5cosχ]+524(ac)3[73cosχ]}

 

 

322(ac)2cosχ{1+32(ac)+323(ac)2[5cosχ]}+123(ac)3{1+32(ac)}+𝒪(a5c5)

 

=

cos3χ{1+32(ac)+323(ac)2[5cosχ]+524(ac)3[73cosχ]+527(ac)4[3cos2χ227cosχ27cosχ+327]}

 

 

+32cos2χ{(ac)+32(ac)2+323(ac)3[5cosχ]+524(ac)4[73cosχ]}

 

 

322cosχ{(ac)2+32(ac)3+323(ac)4[5cosχ]}+123(ac)3+324(ac)4

 

 

+𝒪(a5c5)

Coefficients of Elliptic Integrals

Rewriting the external potential, as provided in the above-stated objective, and evaluating it at the torus surface,

πVDysonGM|𝒪(a4/c4)

=

4K(μ)a+R1{tK}+(a+R1)E(μ)aR1{tE}

 

=

4K(μ)c(R1c)1[1+(ac)(R1c)1]1{tK}+E(μ)a[1+(ac)(R1c)1]{tE},

where,

tK

118(a2c2)cos2(ψ2)1768(ac)4[5+8cosψcos2ψ4cos3ψ4c2RR1cos2ψ],

and,

tE

18(ac)2cosψ1192(ac)4[2cos2ψ4cosψ+2c2RR1cos2ψ].

Given our derived power-series expressions for various trigonometric functions, these coefficients can be rewritten as,

tK

=

1124(ac)2(1+cosψ)+1263(ac)3(R1c)1(2cos2ψ1)1283(ac)4[5+8cosψcos2ψ4cos3ψ]

 

=

1124(ac)2{1cosχ+12(ac)(1cosχ)+123(ac)2[23cosχ+cos2χ]}

 

 

+1263(ac)3(R1c)1{2[cos2χ+(ac)(cos2χcosχ)]1}

 

 

1283(ac)4{58cosχcos2χ+4cos3χ}+𝒪(a5c5)

 

=

1124(ac)2[1cosχ]125(ac)3[1cosχ]127(ac)4[23cosχ+cos2χ]

 

 

+1263(R1c)1{(ac)3(2cos2χ1)+2(ac)4(cos2χcosχ)}

 

 

1283(ac)4{58cosχcos2χ+4cos3χ}+𝒪(a5c5)

 

=

1124(ac)2[1cosχ]+1263(ac)3{(R1c)1(2cos2χ1)23(1cosχ)}

 

 

+1283(ac)4{23(R1c)1(cos2χcosχ)[58cosχcos2χ+4cos3χ]23[23cosχ+cos2χ]}+𝒪(a5c5)

 

=

1124(ac)2(1cosχ)+1263(ac)3[(R1c)1(2cos2χ1)6+6cosχ]

 

 

+1283(ac)4[23(R1c)1(cos2χcosχ)17+26cosχ5cos2χ4cos3χ]+𝒪(a5c5);

and,

tE

=

18(ac)2cosψ1253(ac)3(R1c)1{2[cos2ψ]1}1253(ac)4[cos2ψ2cosψ]

 

=

18(ac)2{cosχ+12(ac)(1cosχ)+123(ac)2[23cosχ+cos2χ]}

 

 

1253(ac)3(R1c)1{2[cos2χ+(ac)(cos2χcosχ)]1}1253(ac)4[cos2χ+2cosχ]+𝒪(a5c5)

 

=

123(ac)2cosχ+124(ac)3(1cosχ)+126(ac)4[23cosχ+cos2χ]

 

 

1253(ac)3(R1c)1[2cos2χ+2(ac)(cos2χcosχ)1]1253(ac)4(cos2χ+2cosχ)+𝒪(a5c5)

 

=

123(ac)2cosχ+124(ac)3(1cosχ)1253(ac)3(R1c)1(2cos2χ1)

 

 

1243(ac)4(R1c)1(cos2χcosχ)+1263(ac)4[613cosχ+cos2χ]+𝒪(a5c5)

 

=

123(ac)2cosχ+1253(ac)3[6(1cosχ)(R1c)1(2cos2χ1)]

 

 

+1263(ac)4[(613cosχ+cos2χ)4(R1c)1(cos2χcosχ)]+𝒪(a5c5).

Now, inserting to the appropriate order the above expression for the ratio, R1/c — namely,

(R1c)1

=

12+122(ac)cosχ+𝒪(a2c2),

we have,

tK

=

1124(ac)2(1cosχ)+1263(ac)3{[12+122(ac)cosχ](2cos2χ1)6+6cosχ}

 

 

+1283(ac)4{22(cos2χcosχ)17+26cosχ5cos2χ4cos3χ}+𝒪(a5c5)

 

=

1124(ac)2(1cosχ)+1273(ac)3(13+12cosχ+2cos2χ)+1283(ac)4(2cos3χcosχ)

 

 

+1283(ac)4(17+22cosχcos2χ4cos3χ)+𝒪(a5c5)

 

=

1124(ac)2(1cosχ)+1273(ac)3(13+12cosχ+2cos2χ)

 

 

+1283(ac)4(17+21cosχcos2χ2cos3χ)+𝒪(a5c5);

and,

tE

=

123(ac)2cosχ+1253(ac)3{6(1cosχ)+[12+122(ac)cosχ](12cos2χ)}

 

 

+1263(ac)4{(613cosχ+cos2χ)+2(cosχcos2χ)}+𝒪(a5c5)

 

=

123(ac)2cosχ+1263(ac)3(1312cosχ2cos2χ)+1273(ac)4(cosχ2cos3χ)

 

 

+1263(ac)4(611cosχcos2χ)+𝒪(a5c5)

 

=

123(ac)2cosχ+1263(ac)3(1312cosχ2cos2χ)

 

 

+1273(ac)4(1221cosχ2cos2χ2cos3χ)+𝒪(a5c5).

Alternate "Small" Argument of Elliptic Integrals

Defining the "small parameter,"

k

1μ2

 

=

[1(eη1eη+1)2]1/2

 

=

[(eη+1)2(eη1)2(eη+1)2]1/2

 

=

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

 

=

[4eη(eη+1)2]1/2=[4eη(1+eη)2]1/2

 

=

2(RR1)1/2[1+RR1]1.

At the surface of the torus, where R=a, we therefore have,

k

=

2(ac)1/2(R1c)1/2[1+accR1]1

4k

=

2(ac)1/2(R1c)1/2[1+accR1]

 

=

2(ca)1/2[4+(ac)24(ac)cosχ]1/4{1+ac[4+(ac)24(ac)cosχ]1/2}

 

=

(23ca)1/2[1(ac)cosχ+14(ac)2]1/4{1+a2c[1(ac)cosχ+(a2c)2]1/2}

ln4k

=

12ln(23ca)+14ln[1(ac)cosχ+14(ac)2]+ln{1+a2c[1(ac)cosχ+(a2c)2]1/2}

 

12ln(23ca)+14ln[1(ac)cosχ]+ln[1+a2c]

 

12ln(23ca)+a2c[112cosχ]

Alternatively, if — as above — we adopt the shorthand notation, γR1/(2c), we can write,

Summary

γ

=

112(ac)cosχ+123(ac)2(1cos2χ)+124(ac)3(cosχcos3χ)+127(ac)4[1+6cos2χ5cos4χ]+𝒪(a5c5);

1γ

=

1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]+𝒪(a5c5).

k

=

(2ac)1/2γ1/2[1+(ac)12γ]1

 

=

(2ac)1/2γ1/2[1+(ac)1γ+(ac)214γ2]1/2

 

=

(2ac)1/2[γ+(ac)+14(ac)21γ]1/2

(k)2m

=

(2ac)m[γ+(ac)+14(ac)21γ]m=(2ac)mΓm

4k

=

(23ca)1/2[γ+(ac)+14(ac)21γ]1/2=(23ca)1/2Γ1/2

ln4k

=

12ln(23ca)+12lnΓ,

where,

Γ

γ+(ac)+14(ac)21γ

 

=

112(ac)cosχ+123(ac)2(1cos2χ)+124(ac)3(cosχcos3χ)+127(ac)4[1+6cos2χ5cos4χ]+(ac)

 

 

+122(ac)2{1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]}+𝒪(a5c5)

 

=

1+(ac)(112cosχ)+123(ac)2(1cos2χ)+124(ac)3(cosχcos3χ)+127(ac)4(1+6cos2χ5cos4χ)

 

 

+122(ac)2+123(ac)3cosχ+125(ac)4(3cos2χ1)+𝒪(a5c5)

 

=

1+(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)+127(ac)4(5+18cos2χ5cos4χ)+𝒪(a5c5).

Now, if we adopt the shorthand notation,

g(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)+127(ac)4(5+18cos2χ5cos4χ)+𝒪(a5c5),

we also have,

lnΓ=ln(1+g)

=

g12g2+13g314g4+𝒪(a5c5)

 

=

(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)+127(ac)4(5+18cos2χ5cos4χ)

 

 

12[(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)]2

 

 

+13[(ac)(112cosχ)+123(ac)2(3cos2χ)]314[(ac)(112cosχ)]4+𝒪(a5c5)

 

=

(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)+127(ac)4(5+18cos2χ5cos4χ)

 

 

12{122(ac)2(2cosχ)2+124(ac)3(3cos2χ)(2cosχ)+125(ac)4(3cosχcos3χ)(2cosχ)

 

 

+123(ac)3(3cos2χ)(112cosχ)+126(ac)4(3cos2χ)(3cos2χ)+124(ac)4(3cosχcos3χ)(112cosχ)}

 

 

+13[(ac)(112cosχ)+123(ac)2(3cos2χ)][(ac)2(112cosχ)2+122(ac)3(112cosχ)(3cos2χ)]126(ac)4(2cosχ)4+𝒪(a5c5)

 

=

12(ac)(2cosχ)+123(ac)2(3cos2χ)123(ac)2(2cosχ)2

 

 

+124(ac)3(3cosχcos3χ)125(ac)3(3cos2χ)(2cosχ)125(ac)3(3cos2χ)(2cosχ)+1233(ac)3(2cosχ)2(2cosχ)

 

 

126(ac)4(3cosχcos3χ)(2cosχ)127(ac)4(3cos2χ)(3cos2χ)126(ac)4(3cosχcos3χ)(2cosχ)

 

 

+1243(ac)4(2cosχ)(3cos2χ)(2cosχ)+1253(ac)4(3cos2χ)(2cosχ)2+127(ac)4(5+18cos2χ5cos4χ)126(ac)4(2cosχ)4+𝒪(a5c5)

And,

(k)2

=

(2ac)Γ1

 

=

(2ac)(1+g)1

 

=

(2ac){1g+g2g3}+𝒪(a5c5)

 

=

(2ac){1[(ac)(112cosχ)+123(ac)2(3cos2χ)+124(ac)3(3cosχcos3χ)]

 

 

+[(ac)(112cosχ)+123(ac)2(3cos2χ)]2[(ac)(112cosχ)]3}+𝒪(a5c5)

 

=

(2ac){112(ac)(2cosχ)+123(ac)2[(3cos2χ)+2(2cosχ)2]

 

 

+124(ac)3[(3cosχcos3χ)+(2cosχ)(3cos2χ)(2cosχ)3]}+𝒪(a5c5)

And,

(k)4

=

(2ac)2Γ2

 

=

(2ac)2(1+g)2

 

=

(2ac)2{12g+3g2}+𝒪(a5c5)

 

=

(2ac)2{12[(ac)(112cosχ)+123(ac)2(3cos2χ)]+3[(ac)(112cosχ)]2}+𝒪(a5c5)

 

=

(2ac)2{1(ac)(2cosχ)+122(ac)2[3(2cosχ)2(3cos2χ)]}+𝒪(a5c5)

 

=

4(ac)2(ac)3(84cosχ)+(ac)4(912cosχ+4cos2χ)+𝒪(a5c5)

Elliptic Integral Expressions

Hence, drawing from our set of Key Expressions for the complete elliptic integral of the first kind, specifically,

K(μ)

=

ln4k'+122(ln4k'212)(k)2+(1324)2(ln4k'212234)(k)4+(135246)2(ln4k'212234256)(k)6+

Gradshteyn & Ryzhik (1965), §8.113.3

where:   k'(1μ2)1/2

we can write,

K(μ)

=

ln4k'+122(ln4k'1)k2+3226(ln4k'76)k4+5228(ln4k'3730)k6+

Now, we recognize that,

(a2)4K(μ)R1+a

=

(a2c)4K(μ)[R1c+ac]1

 

=

(ac)2K(μ)(R1c)1[1+ac(R1c)1]1

 

(ac)(R1c)1[1+ac(R1c)1]1{2ln4k'+12(ln4k'1)k2}.


Also, drawing from our set of Key Expressions for the complete elliptic integral of the second kind, specifically,

E(μ)

=

1+12(ln4k112)(k)2+123224(ln4k212134)(k)4+1232522426(ln4k212234156)(k)6+

Gradshteyn & Ryzhik (1965), §8.114.3

where:   k'(1μ2)1/2

we have,

E(μ)

=

1+12(ln4k12)(k)2+324(ln4k1134)(k)4+325273(ln4k11231235)(k)6+

(a2)(R1+R)E(μ)RR1

12[1+ac(R1c)1]{1+12(ln4k12)(k)2}.

External Potential at Torus Surface

Initial Low Resolution

Hence,

(a2)V2

[118(a2c2)cos2(ψ2)](ac)(R1c)1[1+ac(R1c)1]1{2ln4k'+12(ln4k'1)k2}

 

 

+[18(a2c2)cosψ]12[1+ac(R1c)1]{1+12(ln4k12)(k)2}

cV2

[118(a2c2)cos2(ψ2)0](R1c)1[1+ac(R1c)1]1{4ln4k'+(ln4k'1)k2}

 

 

+[18(ac)cosψ][1+ac(R1c)1]{1+12(ln4k12)(k)2}.

Hence,

cV2

(R1c)1[1+ac(R1c)1]1{4ln4k'+(ln4k'1)k2}

 

 

+[18(ac)cosψ][1+ac(R1c)1]{1+12(ln4k12)(k)2}

 

{1+12(ac)cosχ+(ac)2[38cos2χ18]}[112(ac)+14(ac)2(1cosχ)]{2ln4k'+12(ln4k'1)k2}

 

 

+18(ac)cosψ[1+12(ac)+14(ac)2cosχ]{1+12(ln4k12)(k)2}.

To order (a/c)1, this gives,

cV2

2ln4k'+12(ln4k'1)k2+(ac)(cosχ1)ln4k'+18(ac)cosψ

 

2ln4k'+(ac){1+cosχln4k'+18cosψ}

 

2ln4k'+(ac){1+cosχ[ln4k'18]}

 

ln(23ca)+ac[112cosχ]+(ac){12cosχ[ln(23ca)14]1}

 

ln(23ca)+ac{112cosχ+12cosχ[ln(23ca)14]1}

 

ln(23ca)+12(ac)[ln(23ca)54]cosχ.



We are trying to match equation (6) in Dyson's (1893b) "Part II", that is,

V2πa2

=

ln(8ca)+12(ac)[ln(8ca)54]cosχ+{116[ln(8ca)52]+316[ln(8ca)+17367236]cos2χ}(a2c2)

 

 

+{332[ln(8ca)2512]cosχ+564[ln(8ca)+7244824]cos3χ}(a3c3)

 

 

+{9256[ln(8ca)2]+7128[ln(8ca)191682]cos2χ+351024[ln(8ca)2+19120]cos4χ}(a4c4)+

High Resolution

Summary

2(R1c)1

=

1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]+𝒪(a5c5).

[1+(ac)(R1c)1]1

=

112(ac)+122(ac)2(1cosχ)+123(ac)3(2cosχ3cos2χ)+125(ac)4(9cosχ+8cos2χ5cos3χ)+𝒪(a5c5).

tK

=

1124(ac)2(1cosχ)+1273(ac)3(13+12cosχ+2cos2χ)+1283(ac)4(17+21cosχcos2χ2cos3χ)+𝒪(a5c5).

1+(ac)(R1c)1

=

1+12(ac)+122(ac)2cosχ+124(ac)3(3cos2χ1)+125(ac)4(5cos3χ3cosχ)+𝒪(a5c5).

tE

=

123(ac)2cosχ+1263(ac)3(1312cosχ2cos2χ)+1273(ac)4(1221cosχ2cos2χ2cos3χ)+fE5(ac)5+𝒪(a6c6).


πVDysonGM/c

=

2K(μ)[2(R1c)1][1+(ac)(R1c)1]1{tK}+E(μ)[1+(ac)(R1c)1]{(ca)tE}

 

=

2K(μ){2(R1c)1}

 

 

×{[1+(ac)(R1c)1]1}

 

 

×{tK}

 

 

+E(μ){1+(ac)(R1c)1}

 

 

×{(ca)tE}

 

=

2K(μ){1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]}

 

 

×{112(ac)+122(ac)2(1cosχ)+123(ac)3(2cosχ3cos2χ)+125(ac)4(9cosχ+8cos2χ5cos3χ)}

 

 

×{1124(ac)2(1cosχ)+1273(ac)3(13+12cosχ+2cos2χ)+1283(ac)4(17+21cosχcos2χ2cos3χ)}

 

 

+E(μ){1+12(ac)+122(ac)2cosχ+124(ac)3[3cos2χ1]+125(ac)4[5cos3χ3cosχ]}

 

 

×{123(ac)cosχ+1263(ac)2(1312cosχ2cos2χ)+1273(ac)3(1221cosχ2cos2χ2cos3χ)+fE5(ac)4}+𝒪(a5c5)

 

=

2K(μ){1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]}

 

 

×{1124(ac)2(1cosχ)+1273(ac)3(13+12cosχ+2cos2χ)+1283(ac)4(17+21cosχcos2χ2cos3χ)

 

 

12(ac)+125(ac)3(1cosχ)1283(ac)4(13+12cosχ+2cos2χ)+122(ac)2(1cosχ)126(ac)4(1cosχ)2

 

 

+123(ac)3(2cosχ3cos2χ)+125(ac)4(9cosχ+8cos2χ5cos3χ)}

 

 

+E(μ){123(ac)cosχ[1+12(ac)+122(ac)2cosχ+124(ac)3(3cos2χ1)]

 

 

+1263(ac)2(1312cosχ2cos2χ)[1+12(ac)+122(ac)2cosχ]

 

 

+1273(ac)3(1221cosχ2cos2χ2cos3χ)[1+12(ac)]+fE5(ac)4}+𝒪(a5c5)

 

=

2K(μ){1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]+127(ac)4[330cos2χ+35cos4χ]}

 

 

×{112(ac)+324(ac)2(1cosχ)+1273(ac)3[(13+12cosχ+2cos2χ)+12(1cosχ)+48(2cosχ3cos2χ)]

 

 

+1283(ac)4[(17+21cosχcos2χ2cos3χ)(13+12cosχ+2cos2χ)12(1cosχ)2+24(9cosχ+8cos2χ5cos3χ)]}

 

 

+E(μ){123(ac)cosχ124(ac)2cosχ125(ac)3cos2χ127(ac)4(3cos3χcosχ)

 

 

+1263(ac)2(1312cosχ2cos2χ)+1273(ac)3(1312cosχ2cos2χ)+1283(ac)4(1312cosχ2cos2χ)cosχ

 

 

+1273(ac)3(1221cosχ2cos2χ2cos3χ)+1283(ac)4(1221cosχ2cos2χ2cos3χ)+fE5(ac)4}+𝒪(a5c5)

 

=

2K(μ){1+12(ac)cosχ+123(ac)2(3cos2χ1)+124(ac)3(5cos3χ3cosχ)+127(ac)4(330cos2χ+35cos4χ)}

 

 

×{112(ac)+324(ac)2(1cosχ)+1273(ac)3[(13+12cosχ+2cos2χ)+(1212cosχ)+(96cosχ144cos2χ)]

 

 

+1283(ac)4[(17+21cosχcos2χ2cos3χ)+(1312cosχ2cos2χ)+(12+24cosχ12cos2χ)+(216cosχ+192cos2χ120cos3χ)]}

 

 

+E(μ){123(ac)cosχ+1263(ac)2[(1312cosχ2cos2χ)12cosχ]

 

 

+1273(ac)3[(1312cosχ2cos2χ)+(1221cosχ2cos2χ2cos3χ)12cos2χ]

 

 

+1283(ac)4[(13cosχ12cos2χ2cos3χ)+(1221cosχ2cos2χ2cos3χ)6(3cos3χcosχ)+283fE5]}+𝒪(a5c5)

That is,

πVDysonGM/c

=

2K(μ){1+12(ac)cosχ+123(ac)2(3cos2χ1)+124(ac)3(5cos3χ3cosχ)+127(ac)4(330cos2χ+35cos4χ)}

 

 

×{112(ac)+324(ac)2(1cosχ)+1273(ac)3(1+96cosχ142cos2χ)+1283(ac)4(16183cosχ+177cos2χ122cos3χ)}

 

 

+E(μ){123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)+1273(ac)3(2533cosχ16cos2χ2cos3χ)

 

 

+1283(ac)4(122cosχ14cos2χ22cos3χ+283fE5)}+𝒪(a5c5)

 

=

2K(μ){112(ac)+324(ac)2(1cosχ)+1273(ac)3(1+96cosχ142cos2χ)+1283(ac)4(16183cosχ+177cos2χ122cos3χ)

 

 

+12(ac)cosχ122(ac)2cosχ+325(ac)3(cosχcos2χ)+1283(ac)4(cosχ+96cos2χ142cos3χ)

 

 

+123(ac)2(3cos2χ1)+124(ac)3(13cos2χ)+327(ac)4(1cosχ)(3cos2χ1)

 

 

+124(ac)3(5cos3χ3cosχ)+125(ac)4(3cosχ5cos3χ)+127(ac)4(330cos2χ+35cos4χ)}

 

 

+E(μ){123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)+1273(ac)3(2533cosχ16cos2χ2cos3χ)

 

 

+1283(ac)4(122cosχ14cos2χ22cos3χ+283fE5)}+𝒪(a5c5)

 

=

2K(μ){1+12(ac)(cosχ1)+124(ac)2[3(1cosχ)4cosχ+2(3cos2χ1)]

 

 

+1273(ac)3[36(cosχcos2χ)+24(13cosχ3cos2χ+5cos3χ)+(1+96cosχ142cos2χ)]

 

 

+1283(ac)4[(16183cosχ+177cos2χ122cos3χ)+18(1cosχ)(3cos2χ1)+(cosχ+96cos2χ142cos3χ)+24(3cosχ5cos3χ)+6(330cos2χ+35cos4χ)]}

 

 

+E(μ){123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)+1273(ac)3(2533cosχ16cos2χ2cos3χ)

 

 

+1283(ac)4(122cosχ14cos2χ22cos3χ+283fE5)}+𝒪(a5c5)

 

=

2K(μ){1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)+1273(ac)3(23+60cosχ250cos2χ+120cos3χ)

 

 

+1283(ac)4(1694cosχ+147cos2χ194cos3χ+210cos4χ)}

 

 

+E(μ){123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)+1273(ac)3(2533cosχ16cos2χ2cos3χ)

 

 

+1283(ac)4(122cosχ14cos2χ22cos3χ+283fE5)}+𝒪(a5c5)

Insert Expressions for K and E

Summary

k

=

2(ac)1/2(R1c)1/2[1+ac(R1c)1]1

4k

=

2(ac)1/2(R1c)1/2[1+ac(R1c)1]

 

=

2(ca)1/2[4+(ac)24(ac)cosχ]1/4{1+ac[4+(ac)24(ac)cosχ]1/2}

 

=

(23ca)1/2[1(ac)cosχ+14(ac)2]1/4{1+a2c[1(ac)cosχ+(a2c)2]1/2}

ln4k

=

12ln(23ca)+14ln[1(ac)cosχ+14(ac)2]+ln{1+a2c[1(ac)cosχ+(a2c)2]1/2}

 

12ln(23ca)+14ln[1(ac)cosχ]+ln[1+a2c]

 

12ln(23ca)+a2c[112cosχ]


Remember that (see above),

K(μ)

=

ln4k'+122(ln4k'1)k2+3226(ln4k'76)k4+5228(ln4k'3730)k6+

And (see above),

E(μ)

=

1+12(ln4k12)(k)2+324(ln4k1134)(k)4+325273(ln4k11231235)(k)6+


We are trying to match equation (6) in Dyson's (1893b) "Part II", that is,

V2πa2

=

ln(8ca)+12(ac)[ln(8ca)54]cosχ+{116[ln(8ca)52]+316[ln(8ca)+17367236]cos2χ}(a2c2)

 

 

+{332[ln(8ca)2512]cosχ+564[ln(8ca)+7244824]cos3χ}(a3c3)

 

 

+{9256[ln(8ca)2]+7128[ln(8ca)191682]cos2χ+351024[ln(8ca)2+19120]cos4χ}(a4c4)+


To First Order

(k)2

=

(2ac){112(ac)(2cosχ)+123(ac)2[(3cos2χ)+2(2cosχ)2]

 

 

+124(ac)3[(3cosχcos3χ)+(2cosχ)(3cos2χ)(2cosχ)3]}+𝒪(a5c5)

ln4k

=

12ln(23ca)+12lnΓ

 

=

12ln(23ca)+12{12(ac)(2cosχ)+123(ac)2(3cos2χ)123(ac)2(2cosχ)2+𝒪(a3c3)}

Hence,

πVDysonGM/c

=

2K(μ){1+12(ac)(cosχ1)}+E(μ){123(ac)cosχ}+𝒪(a2c2)

 

[1+12(ac)(cosχ1)]{2K(μ)}123(ac)cosχ{E(μ)}

 

[1+12(ac)(cosχ1)]{2ln4k+12[ln4k1]k'2}123(ac)cosχ{1+12[ln4k12]k'20}

 

{2ln4k+12[ln4k1](2ac)}+12(ac)(cosχ1){2ln4k+12[ln4k1]k'20}123(ac)cosχ{1+12[ln4k12]k'20}

 

2ln4k+(ac)[ln4k1]+(ac)(cosχ1){ln4k}123(ac)cosχ

 

2ln4k+(ac)[ln4k1+(cosχ1)(ln4k)123cosχ]

 

2ln4k+(ac)[1+cosχ(ln4k)123cosχ]

 

ln(23ca)+12(ac){(2cosχ)2+cosχ[ln(23ca)]14cosχ}

 

ln(23ca)+12(ac)[ln(23ca)54]cosχ

To Second Order

πVDysonGM/c

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{2K(μ)}

 

 

+{123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)}{E(μ)}+𝒪(a3c3)

 

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{2ln4k'+12(ln4k'1)k2+3225(ln4k'76)k4}

 

 

+{123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)}{1+12(ln4k12)(k)2}+𝒪(a3c3)

 

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{2[12ln(23ca)+12lnΓ]+12[12ln(23ca)+12lnΓ1]k2+3225[12ln(23ca)+12lnΓ76]k4}

 

 

+{123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)}{1+12[12ln(23ca)+12lnΓ12](k)2}+𝒪(a3c3)

 

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{[ln(23ca)+lnΓ]+122[ln(23ca)+lnΓ2](2ac)[112(ac)(2cosχ)]+326[3ln(23ca)+3lnΓ07]4(ac)2}

 

 

+{123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)}{1+122[ln(23ca)+lnΓ01](2ac)[112(ac)0(2cosχ)]}+𝒪(a3c3),

where,

lnΓ

=

12(ac)(2cosχ)+123(ac)2(1+4cosχ2cos2χ)+𝒪(a3c3).

Hence,

πVDysonGM/c

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{[ln(23ca)+12(ac)(2cosχ)+123(ac)2(1+4cosχ2cos2χ)]

 

 

+12(ac)[ln(23ca)2+12(ac)(2cosχ)][112(ac)(2cosχ)]+324(ac)2[3ln(23ca)7]}

 

 

+{123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)}{1+12(ac)[ln(23ca)1]}+𝒪(a3c3)

 

=

{1+12(ac)(cosχ1)+124(ac)2(17cosχ+6cos2χ)}{[ln(23ca)+12(ac)(2cosχ)+123(ac)2(1+4cosχ2cos2χ)]

 

 

+12(ac)[ln(23ca)2]+122(ac)2(2cosχ)122(ac)2[ln(23ca)2](2cosχ)+324(ac)2[3ln(23ca)7]}

 

 

123(ac)cosχ+1263(ac)2(1324cosχ2cos2χ)124(ac)2cosχ[ln(23ca)1]+𝒪(a3c3)

 

=

ln(23ca)+12(ac)(2cosχ)123(ac)cosχ+12(ac)[ln(23ca)2]+12(ac)(cosχ1)[ln(23ca)]

 

 

+123(ac)2(1+4cosχ2cos2χ)+122(ac)2(2cosχ)122(ac)2[ln(23ca)2](2cosχ)+324(ac)2[3ln(23ca)7]

 

 

+122(ac)2(cosχ1)(2cosχ)+124(ac)2(17cosχ+6cos2χ)ln(23ca)

 

 

+1263(ac)2(1324cosχ2cos2χ)124(ac)2cosχ[ln(23ca)1]+𝒪(a3c3)

 

=

ln(23ca)+123(ac){4(2cosχ)cosχ+[4ln(23ca)8]+[4ln(23ca)](cosχ1)}

 

 

+1263(ac)2{24(1+4cosχ2cos2χ)+48(2cosχ)48[ln(23ca)2](2cosχ)+36[3ln(23ca)7]

 

 

+48(cosχ1)(2cosχ)+12(17cosχ+6cos2χ)ln(23ca)+(1324cosχ2cos2χ)12cosχ[ln(23ca)1]}+𝒪(a3c3)

 

=

ln(23ca)+123(ac){4(2cosχ)cosχ+[4ln(23ca)8]+[4ln(23ca)](cosχ1)}

 

 

+1263(ac)2{(24+96cosχ48cos2χ)+(9648cosχ)+[2ln(23ca)](9648cosχ)+[108ln(23ca)252]

 

 

+48(3cosχcos2χ2)+(1284cosχ+72cos2χ)ln(23ca)+(1324cosχ2cos2χ)+12cosχ12cosχ[ln(23ca)]}+𝒪(a3c3)

 

=

ln(23ca)+12(ac)[ln(23ca)54]cosχ

 

 

+1263(ac)2{24+96cosχ48cos2χ+9648cosχ+19225296cosχ+(1312cosχ2cos2χ)+(144cosχ48cos2χ96)

 

 

+36[ln(23ca)]cosχ+(1284cosχ+72cos2χ96+108)[ln(23ca)]}+𝒪(a3c3)

 

=

ln(23ca)+12(ac)[ln(23ca)54]cosχ

 

 

+1263(ac)2{71+84cosχ98cos2χ+(2448cosχ+72cos2χ)[ln(23ca)]}+𝒪(a3c3)

 

=

ln(23ca)+12(ac)[ln(23ca)54]cosχ+1263(ac)2{71+84cosχ98cos2χ+24ln(23ca)(12cosχ+3cos2χ)}+𝒪(a3c3)

In an effort to compare this expression with equation (6) from Dyson's (1893b) "Part II", we should make the substitutions,

ln(23ca)(λ+2)       and       2cos2χ1+cos2χ.

This means,

πVDysonGM/c|𝒪(a2/c2)

=

+1263(ac)2{71+84cosχ49(1+cos2χ)+24(λ+2)(12cosχ)+36(λ+2)(1+cos2χ)}

 

=

+1263(ac)2{71+84cosχ4949cos2χ+24(λ+22λcosχ4cosχ)+36(λ+2+λcos2χ+2cos2χ)}

 

=

+1263(ac)2{71+84cosχ4949cos2χ+24λ+4848λcosχ96cosχ+36λ+72+36λcos2χ+72cos2χ}

 

=

+1263(ac)2{60λ48λcosχ12cosχ+36λcos2χ+23cos2χ}

 

=

+(ac)2{5λ16(4λ+1)16cosχ+3(λ+2336)16cos2χ}.

This expression differs from the 2nd-order term in Dyson's equation (6) by the amount,

Δ|𝒪(a2/c2)

=

(ac)2{5λ16(4λ+1)16cosχ+3(λ+2336)16cos2χ}(ac)2{λ1216+3(λ+1736)16cos2χ}

 

=

11612(ac)2{60λ(48λ+12)cosχ+(36λ+23)cos2χ}11612(ac)2{12λ6+(36λ+17)cos2χ}

 

=

11612(ac)2{48λ+6(48λ+12)cosχ+(6)cos2χ}

 

=

125(ac)2{(cos2χ1)(8λ+2)(1+cosχ)}

 

=

124(ac)2{(cosχ1)(4λ+1)}(1+cosχ)

 

=

124(ac)2(cosχ24λ)(1+cosχ)

Interior Potential

In equation (9) on p. 1050 of his "Part II", Dyson (1893b) presents the following power-series expression for the gravitational potential at points inside the torus:

V

=

2πa2{L+12(1R2a2)+ac[(L1)2(Ra)R38a3]cosχ

 

 

+a2c2[(L14)16+(L1)8(R2a2)364(R4a4)+3(L54)16(R2a2)cos2χ596(R4a4)cos2χ]+}

where,

L

ln(8ca).

Note that, for the example illustrated above, a/c=2/5 and, hence, L=ln(20)=2.99573. Therefore, at any point on the surface of this example torus,

V2πa2

ln(20)+25[ln(20)1218]cosχ

 

 

+2252[4ln(20)164+ln(20)18364+12ln(20)1564cos2χ10364cos2χ]

 

=

ln(20)+240[4ln(20)5]cosχ+222652{[4ln(20)1]+[8ln(20)8]3+[12ln(20)15]cos2χ(103)cos2χ}

 

=

ln(20)+120[4ln(20)5]cosχ+12452{4ln(20)10+13[36ln(20)55]cos2χ}


Red Contour

Let's construct an equipotential contour that extends the red contour into the interior region. Let's begin by evaluating Dyson's interior potential expression at the coordinate location where the red contour touches the surface of the torus. According to Column #5, this point on the surface has coordinates, (ϖ,z)=(0.867,0.377); or, equivalently, R=a=0.4 and,

cosχ

=

1ϖR=0.3325

χ

=

cos1(0.3325)=1.2318.

Hence, for this specific point on the torus surface, we find,

V2πa2

3+[720]cosχ+12452{2+[533]cos2χ}=3.0870.

Contour
(column #)
V2 (external) Surface ϖ χ (radians)
blue (3) 0.9120 1.117 1.8676
orange (4) 0.961 0.929 1.392
red (5) 0.9800 0.867 1.2318
light green (6) 0.9800 0.838 1.1538
light blue (7) 1.0212 0.749 0.8925

See Also

The following quotes have been taken from Petroff & Horatschek (2008):

§1:   "The problem of the self-gravitating ring captured the interest of such renowned scientists as Kowalewsky (1885), Poincaré (1885a,b,c) and Dyson (1892, 1893). Each of them tackled the problem of an axially symmetric, homogeneous ring in equilibrium by expanding it about the thin ring limit. In particular, Dyson provided a solution to fourth order in the parameter σ=a/b, where a=rt provides a measure for the radius of the cross-section of the ring and b=ϖt the distance of the cross-section's centre of mass from the axis of rotation."

§7:   "In their work on homogeneous rings, Poincaré and Kowalewsky, whose results disagreed to first order, both had made mistakes as Dyson has shown. His result to fourth order is also erroneous as we point out in Appendix B."

  1. Shortly after their equation (3.2), Marcus, Press & Teukolsky make the following statement: "… we know that an equilibrium incompressible configuration must rotate uniformly on cylinders (the famous "Poincaré-Wavre" theorem, cf. Tassoul 1977, &Sect;4.3) …"
Tiled Menu

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