<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article">
 <front>
  <journal-meta>
   <journal-id journal-id-type="publisher-id">
    jhepgc
   </journal-id>
   <journal-title-group>
    <journal-title>
     Journal of High Energy Physics, Gravitation and Cosmology
    </journal-title>
   </journal-title-group>
   <issn pub-type="epub">
    2380-4327
   </issn>
   <issn publication-format="print">
    2380-4335
   </issn>
   <publisher>
    <publisher-name>
     Scientific Research Publishing
    </publisher-name>
   </publisher>
  </journal-meta>
  <article-meta>
   <article-id pub-id-type="doi">
    10.4236/jhepgc.2025.111003
   </article-id>
   <article-id pub-id-type="publisher-id">
    jhepgc-138767
   </article-id>
   <article-categories>
    <subj-group subj-group-type="heading">
     <subject>
      Articles
     </subject>
    </subj-group>
    <subj-group subj-group-type="Discipline-v2">
     <subject>
      Physics 
     </subject>
     <subject>
       Mathematics
     </subject>
    </subj-group>
   </article-categories>
   <title-group>
    The Origin of Rings
   </title-group>
   <contrib-group>
    <contrib contrib-type="author" xlink:type="simple">
     <name name-style="western">
      <surname>
       Mark
      </surname>
      <given-names>
       Linick
      </given-names>
     </name>
    </contrib>
   </contrib-group> 
   <aff id="affnull">
    <addr-line>
     aLong Beach, CA, USA
    </addr-line> 
   </aff> 
   <pub-date pub-type="epub">
    <day>
     02
    </day> 
    <month>
     01
    </month>
    <year>
     2025
    </year>
   </pub-date> 
   <volume>
    11
   </volume> 
   <issue>
    01
   </issue>
   <fpage>
    19
   </fpage>
   <lpage>
    27
   </lpage>
   <history>
    <date date-type="received">
     <day>
      7,
     </day>
     <month>
      October
     </month>
     <year>
      2024
     </year>
    </date>
    <date date-type="published">
     <day>
      5,
     </day>
     <month>
      October
     </month>
     <year>
      2024
     </year> 
    </date> 
    <date date-type="accepted">
     <day>
      5,
     </day>
     <month>
      January
     </month>
     <year>
      2025
     </year> 
    </date>
   </history>
   <permissions>
    <copyright-statement>
     © Copyright 2014 by authors and Scientific Research Publishing Inc. 
    </copyright-statement>
    <copyright-year>
     2014
    </copyright-year>
    <license>
     <license-p>
      This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
     </license-p>
    </license>
   </permissions>
   <abstract>
    The origin of rings is one of the final arguments in a series explaining solar system formation. Liquid moons placed in circular orbits close to the planet (inside the Roche limit) will deposit ring material as they migrate. Once outside the limit, moons retain this material. Finally, the shape of the moon, Pan, is explained, using the same mechanism that forms rings.
   </abstract>
   <kwd-group> 
    <kwd>
     Liquid
    </kwd> 
    <kwd>
      Moons
    </kwd> 
    <kwd>
      Roche Limit
    </kwd>
   </kwd-group>
  </article-meta>
 </front>
 <body>
  <sec id="s1">
   <title>1. Introduction</title>
   <p>Planetary rings are fascinating, where their beauty only adds to the mystery of their formation. However, the mechanism of formation is one of the least understood processes related to planet formation and evolution. For a recent explanation of the formation process, see Charnoz <xref ref-type="bibr" rid="scirp.138767-1">
     [1]
    </xref>. A variety of mechanisms leading to ring formation are considered.</p>
   <p>The origin of rings is one of the last arguments in a series explaining solar system formation (craters and tidal locking are the others). It is based on formations beginning as a liquid. As explained earlier, solar systems are believed to have a liquid phase early in development. This difficult concept is made more believable by demonstrating the origin of rings. It starts with liquid moons placed in circular orbits about the planet. Liquids far from the planet stay intact (moons)—liquids close to a planet form rings (inside the Roche limit).</p>
   <p>The shape of the liquid (moon) is approximated as an ellipsoid and the gravity is determined on the surface at the semi-major axis. The tidal force is calculated and the distance is found where free objects are no longer held to the surface—the Roche limit. Particles free of the (liquid) surface are immediately placed in orbit, forming rings. Finally, to demonstrate the mechanism, we show how the moon, Pan, acquired its strange shape.</p>
  </sec><sec id="s2">
   <title>2. Motivation</title>
   <p>The distance from the planet, where the tidal force becomes dominant over the self gravity of the moon, is called the Roche limit. For scientists, it seems there are only two Roche limits of interest-one for rigid moons (moons held together by chemical bonds) and one for fluid moons (moons held together by self gravity-the so called “rubble” moon).</p>
   <p>I’m not sure what a rubble moon is, a comet, maybe? If these moons exist, they are assumed to enter the Roche limit from the outside. It would be difficult for the debris of any moon entering from the outside, to obtain the orbit we seek for rings (circular). An example of this type of moon may be the Comet Shoemaker-Levy 9. Generally, the path the particles take, is determined by the orbital character of the moons that produce them. We will dismiss the “rubble” moon for this argument in favor of the rigid moon. That is, the rigid moon holds the clue, if interpreted correctly.</p>
   <p>Most rigid moons (the moons we actually see) are in circular orbits. If inside the Roche limit, free particles on the surface of these moons will obtain a similar orbit, with initial conditions inherited from the moon. In <xref ref-type="fig" rid="fig1(a)">
     Figure 1(a)
    </xref>, two objects released from opposite sides follow sleight elliptical orbits. In <xref ref-type="fig" rid="fig1(b)">
     Figure 1(b)
    </xref>, objects released at intervals show the pattern of rings. This is what we want. This would form rings.</p>
   <fig id="fig1" position="float">
    <label>Figure 1</label>
    <caption>
     <title>Figure 1. Ring formation of moon inside Roche limit: (a) Released particles follow elliptical paths shown as dashed lines; (b) Particles released at intervals form rings.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId16.jpeg?20250113050637" />
   </fig>
   <p>The problem is, rigid moons have no free particles (surface particles). If they did, the Roche limit for a rigid moon is 100,000 Km, well inside the distance to the ring’s edge at 137,000 Km—for Saturn, that is. How did rings form outside this distance?</p>
   <p>The mistake, I believe, is using a rigid moon (a solid) to calculate the Roche limit. We want to assume the moon was in a molten state (a liquid) when first placed in orbit. Then the tidal force will influence the shape (ellipsoidal) and the surface gravity of the moon (the liquid is held together by chemical bonds). The rigid moon (a sphere) is a later stage product, formed after cooling and migration.</p>
   <p>The molten moon placed in obit must possess two properties: 1) a Roche limit that extends to the current edge of the rings (~137,000 Km) and 2) a mechanism that converts the moon to the free particles found in rings. (How the liquid moon is placed in orbit was explained earlier—The Origin of Moons.)</p>
  </sec><sec id="s3">
   <title>3. Tidal Shape of Liquids</title>
   <p>Since the shape of the liquid will change with distance in orbit, we need to generate ellipsoids [Roche ellipsoid?] to examine the tidal behavior. Roche ellipsoids are unnecessary if we can demonstrate with a simple approximation.</p>
   <p>For a liquid, there is little to no tensile strength. If the change in the radius is simply proportional to the forces at the surface, particularly at the axis, then the equations of the major ( 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math>) and minor ( 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       b 
     </mi> 
    </math>) axis are as follows:</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        a 
      </mi> 
      <mn>
        1 
      </mn> 
      <mo>
        = 
      </mo> 
      <msub> 
       <mi>
         r 
       </mi> 
       <mrow> 
        <mi>
          m 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          n 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        ∗ 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mn>
          1 
        </mn> 
        <mo>
          + 
        </mo> 
        <mrow> 
         <mrow> 
          <mi>
            F 
          </mi> 
          <mn>
            1 
          </mn> 
         </mrow> 
         <mo>
           / 
         </mo> 
         <mrow> 
          <msub> 
           <mi>
             F 
           </mi> 
           <mrow> 
            <mi>
              m 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              n 
            </mi> 
           </mrow> 
          </msub> 
         </mrow> 
        </mrow> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(1)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        b 
      </mi> 
      <mn>
        1 
      </mn> 
      <mo>
        = 
      </mo> 
      <msub> 
       <mi>
         r 
       </mi> 
       <mrow> 
        <mi>
          m 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          n 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        ∗ 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mn>
          1 
        </mn> 
        <mo>
          − 
        </mo> 
        <mrow> 
         <mrow> 
          <mi>
            F 
          </mi> 
          <mn>
            2 
          </mn> 
         </mrow> 
         <mo>
           / 
         </mo> 
         <mrow> 
          <msub> 
           <mi>
             F 
           </mi> 
           <mrow> 
            <mi>
              m 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              n 
            </mi> 
           </mrow> 
          </msub> 
         </mrow> 
        </mrow> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(2)</p>
   <p>where</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         F 
       </mi> 
       <mrow> 
        <mi>
          m 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          n 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          m 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msubsup> 
         <mi>
           r 
         </mi> 
         <mrow> 
          <mi>
            m 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            n 
          </mi> 
         </mrow> 
         <mn>
           2 
         </mn> 
        </msubsup> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math>(3)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        F 
      </mi> 
      <mn>
        1 
      </mn> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mrow> 
          <mrow> 
           <mo>
             ( 
           </mo> 
           <mrow> 
            <mi>
              d 
            </mi> 
            <mo>
              − 
            </mo> 
            <msub> 
             <mi>
               r 
             </mi> 
             <mrow> 
              <mi>
                m 
              </mi> 
              <mi>
                o 
              </mi> 
              <mi>
                o 
              </mi> 
              <mi>
                n 
              </mi> 
             </mrow> 
            </msub> 
           </mrow> 
           <mo>
             ) 
           </mo> 
          </mrow> 
         </mrow> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        − 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mi>
          d 
        </mi> 
        <mo>
          − 
        </mo> 
        <msub> 
         <mi>
           r 
         </mi> 
         <mrow> 
          <mi>
            m 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            n 
          </mi> 
         </mrow> 
        </msub> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
      <msup> 
       <mi>
         w 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
     </mrow> 
    </math>(4)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        F 
      </mi> 
      <mn>
        2 
      </mn> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mrow> 
         <mrow> 
          <msub> 
           <mi>
             r 
           </mi> 
           <mrow> 
            <mi>
              m 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              o 
            </mi> 
            <mi>
              n 
            </mi> 
           </mrow> 
          </msub> 
         </mrow> 
         <mo>
           / 
         </mo> 
         <mi>
           R 
         </mi> 
        </mrow> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
      <mo>
        ∗ 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           R 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mn>
        , 
      </mn> 
      <mtext>
          
      </mtext> 
      <mtext>
          
      </mtext> 
      <mtext>
        where 
      </mtext> 
      <mtext>
          
      </mtext> 
      <mtext>
          
      </mtext> 
      <mi>
        R 
      </mi> 
      <mo>
        = 
      </mo> 
      <msqrt> 
       <mrow> 
        <msup> 
         <mi>
           d 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
        <mo>
          + 
        </mo> 
        <msubsup> 
         <mi>
           r 
         </mi> 
         <mrow> 
          <mi>
            m 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            n 
          </mi> 
         </mrow> 
         <mn>
           2 
         </mn> 
        </msubsup> 
       </mrow> 
      </msqrt> 
     </mrow> 
    </math>(5)</p>
   <p>Then adjusting for constant volume</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mtext>
        ratio 
      </mtext> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          a 
        </mi> 
        <mn>
          1 
        </mn> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <mi>
          b 
        </mi> 
        <mn>
          1 
        </mn> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math></p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        b 
      </mi> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <msub> 
         <mi>
           r 
         </mi> 
         <mrow> 
          <mi>
            m 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            o 
          </mi> 
          <mi>
            n 
          </mi> 
         </mrow> 
        </msub> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mrow> 
          <mtext>
            ratio 
          </mtext> 
         </mrow> 
         <mrow> 
          <mrow> 
           <mn>
             1 
           </mn> 
           <mo>
             / 
           </mo> 
           <mn>
             3 
           </mn> 
          </mrow> 
         </mrow> 
        </msup> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math></p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        a 
      </mi> 
      <mo>
        = 
      </mo> 
      <mi>
        b 
      </mi> 
      <mo>
        ∗ 
      </mo> 
      <mtext>
        ratio 
      </mtext> 
     </mrow> 
    </math></p>
   <table-wrap id="table1">
    <label>
     <xref ref-type="table" rid="table1">
      Table 1
     </xref></label>
    <caption>
     <title>
      <xref ref-type="bibr" rid="scirp.138767-"></xref>Table 1. Calculated and measured dimensions for some moons of Saturn.</title>
    </caption>
    <table class="MsoTableGrid custom-table" border="0" cellspacing="0" cellpadding="0"> 
     <tr> 
      <td class="custom-bottom-td acenter" width="29.39%"><p style="text-align:center">Saturn Moons</p></td> 
      <td class="custom-bottom-td acenter" width="29.39%"><p style="text-align:center">Calculated</p><p style="text-align:center">a × b (×b) (Km)</p></td> 
      <td class="custom-bottom-td acenter" width="41.21%"><p style="text-align:center">Measured</p><p style="text-align:center">Dimensions (Km)</p></td> 
     </tr> 
     <tr> 
      <td class="custom-top-td acenter" width="29.39%"><p style="text-align:center">Titan</p></td> 
      <td class="custom-top-td acenter" width="29.39%"><p style="text-align:center">5152 × 5150</p></td> 
      <td class="custom-top-td acenter" width="41.21%"><p style="text-align:center">5149 × 5149 × 5150</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Rhea</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">1530 × 1526</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">1530 × 1526 × 1525</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Dione</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">1128 × 1120</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">1128 × 1123 × 1119</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Tethys</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">1077 × 1055</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">1077 × 1057 × 1053</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Enceladus</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">513 × 500</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">513 × 503 × 497</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Mimas</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">416 × 387</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">416 × 393 × 381</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Janus</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">209 × 166</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">203 × 185 × 153</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Epimetheus</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">135 × 107</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">130 × 114 × 106</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Pandora</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">103 × 72</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">104 × 81 × 64</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Prometheus</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">111 × 75</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">136 × 79 × 59</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Atlas</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">40 × 26</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">41 × 35 × 19</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="29.39%"><p style="text-align:center">Pan</p></td> 
      <td class="acenter" width="29.39%"><p style="text-align:center">39 × 24</p></td> 
      <td class="acenter" width="41.21%"><p style="text-align:center">34 × 31 × 20</p></td> 
     </tr> 
    </table>
   </table-wrap>
   <p>The only comparison for a homogeneous liquid was an approximation ( 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        e 
      </mi> 
      <mo>
        ≪ 
      </mo> 
      <mn>
        1 
      </mn> 
     </mrow> 
    </math>) between the Earth and the moon <xref ref-type="bibr" rid="scirp.138767-2">
     [2]
    </xref>, where we both calculated a tidal distortion of 50 meters—between 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math> and 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       b 
     </mi> 
    </math>. Calculations for some moons of Saturn are listed in <xref ref-type="table" rid="table1">
     Table 1
    </xref> along with the measured values <xref ref-type="bibr" rid="scirp.138767-3">
     [3]
    </xref> <xref ref-type="bibr" rid="scirp.138767-4">
     [4]
    </xref>.</p>
   <p>Most moons match well. We conclude the method we use to generate ellipsoids is adequate for what we wish to demonstrate.</p>
  </sec><sec id="s4">
   <title>4. Gravitational Force of Ellipsoid</title>
   <p>We calculate the self gravity of a liquid ellipsoidal moon at points on the surface at the semi-major axis. These are the points closest and furthest from the planet. We demonstrate with a sphere first, then move to the ellipsoid.</p>
   <p>Usually, the gravitational force of a sphere to an external point is accomplished by summing the potential of shells – integrating them from 0 to 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       R 
     </mi> 
    </math>. But this can also be accomplished by summing the force of disks between 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       R 
     </mi> 
    </math> and 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        − 
      </mo> 
      <mi>
        R 
      </mi> 
     </mrow> 
    </math>.</p>
   <p>Start with the gravitational force from a disk centered at the origin, a distance 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       d 
     </mi> 
    </math> from the point 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       P 
     </mi> 
    </math>, as shown in <xref ref-type="fig" rid="fig2">
     Figure 2
    </xref>. The integration is standard and the equations of force are Equation (6) or Equation (7):</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mi>
        t 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mn>
          1 
        </mn> 
        <mo>
          − 
        </mo> 
        <mi>
          cos 
        </mi> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mi>
           α 
         </mi> 
         <mo>
           ) 
         </mo> 
        </mrow> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(6)</p>
   <p>or</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mi>
        t 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mn>
          1 
        </mn> 
        <mo>
          − 
        </mo> 
        <mfrac> 
         <mi>
           d 
         </mi> 
         <mrow> 
          <msqrt> 
           <mrow> 
            <msup> 
             <mi>
               R 
             </mi> 
             <mn>
               2 
             </mn> 
            </msup> 
            <mo>
              + 
            </mo> 
            <msup> 
             <mi>
               d 
             </mi> 
             <mn>
               2 
             </mn> 
            </msup> 
           </mrow> 
          </msqrt> 
         </mrow> 
        </mfrac> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(7)</p>
   <p>where 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       G 
     </mi> 
    </math> is the gravitational constant, 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       m 
     </mi> 
    </math> the mass at 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       P 
     </mi> 
    </math>, 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       σ 
     </mi> 
    </math> the volume density of the disk, 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       R 
     </mi> 
    </math> the disk radius, 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       d 
     </mi> 
    </math> the distance of 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       P 
     </mi> 
    </math> from the disk and 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       t 
     </mi> 
    </math> the disk thickness.</p>
   <fig id="fig2" position="float">
    <label>Figure 2</label>
    <caption>
     <title>Figure 2. Gravitational force from a disk acting on a point P perpendicular to the face.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId73.jpeg?20250113050638" />
   </fig>
   <p>With the radius of the disk and the distance to 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       P 
     </mi> 
    </math> a function of 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       x 
     </mi> 
    </math>, we make the following substitutions in Equation (7) to get,</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mn>
          1 
        </mn> 
        <mo>
          − 
        </mo> 
        <mfrac> 
         <mrow> 
          <mi>
            d 
          </mi> 
          <mo>
            − 
          </mo> 
          <mi>
            x 
          </mi> 
         </mrow> 
         <mrow> 
          <msqrt> 
           <mrow> 
            <mrow> 
             <mo>
               ( 
             </mo> 
             <mrow> 
              <msup> 
               <mi>
                 R 
               </mi> 
               <mn>
                 2 
               </mn> 
              </msup> 
              <mo>
                − 
              </mo> 
              <msup> 
               <mi>
                 x 
               </mi> 
               <mn>
                 2 
               </mn> 
              </msup> 
             </mrow> 
             <mo>
               ) 
             </mo> 
            </mrow> 
            <mo>
              + 
            </mo> 
            <msup> 
             <mrow> 
              <mrow> 
               <mo>
                 ( 
               </mo> 
               <mrow> 
                <mi>
                  d 
                </mi> 
                <mo>
                  − 
                </mo> 
                <mi>
                  x 
                </mi> 
               </mrow> 
               <mo>
                 ) 
               </mo> 
              </mrow> 
             </mrow> 
             <mn>
               2 
             </mn> 
            </msup> 
           </mrow> 
          </msqrt> 
         </mrow> 
        </mfrac> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
      <mi>
        d 
      </mi> 
      <mi>
        x 
      </mi> 
     </mrow> 
    </math>(8)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        = 
      </mo> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mtext>
             d 
           </mtext> 
           <mi>
             x 
           </mi> 
          </mrow> 
         </mrow> 
        </mstyle> 
        <mo>
          − 
        </mo> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mfrac> 
            <mrow> 
             <mrow> 
              <mo>
                ( 
              </mo> 
              <mrow> 
               <mi>
                 d 
               </mi> 
               <mo>
                 − 
               </mo> 
               <mi>
                 x 
               </mi> 
              </mrow> 
              <mo>
                ) 
              </mo> 
             </mrow> 
             <mtext>
               d 
             </mtext> 
             <mi>
               x 
             </mi> 
            </mrow> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mrow> 
                <mo>
                  ( 
                </mo> 
                <mrow> 
                 <msup> 
                  <mi>
                    R 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                 <mo>
                   − 
                 </mo> 
                 <msup> 
                  <mi>
                    x 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                </mrow> 
                <mo>
                  ) 
                </mo> 
               </mrow> 
               <mo>
                 + 
               </mo> 
               <msup> 
                <mrow> 
                 <mrow> 
                  <mo>
                    ( 
                  </mo> 
                  <mrow> 
                   <mi>
                     d 
                   </mi> 
                   <mo>
                     − 
                   </mo> 
                   <mi>
                     x 
                   </mi> 
                  </mrow> 
                  <mo>
                    ) 
                  </mo> 
                 </mrow> 
                </mrow> 
                <mn>
                  2 
                </mn> 
               </msup> 
              </mrow> 
             </msqrt> 
            </mrow> 
           </mfrac> 
          </mrow> 
         </mrow> 
        </mstyle> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(9)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        = 
      </mo> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mn>
          2 
        </mn> 
        <mi>
          R 
        </mi> 
        <mo>
          − 
        </mo> 
        <mn>
          2 
        </mn> 
        <mi>
          R 
        </mi> 
        <mo>
          + 
        </mo> 
        <mrow> 
         <mrow> 
          <mn>
            2 
          </mn> 
          <msup> 
           <mi>
             R 
           </mi> 
           <mn>
             3 
           </mn> 
          </msup> 
         </mrow> 
         <mo>
           / 
         </mo> 
         <mrow> 
          <mn>
            3 
          </mn> 
          <msup> 
           <mi>
             d 
           </mi> 
           <mn>
             2 
           </mn> 
          </msup> 
         </mrow> 
        </mrow> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(10)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
        <mi>
          m 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           d 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math> (11)</p>
   <p>Substituting 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       R 
     </mi> 
    </math> for 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       d 
     </mi> 
    </math>, the force on the surface of the sphere can be factored as follows:</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mtext>
             d 
           </mtext> 
           <mi>
             x 
           </mi> 
          </mrow> 
         </mrow> 
        </mstyle> 
        <mo>
          − 
        </mo> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mfrac> 
            <mrow> 
             <mrow> 
              <mo>
                ( 
              </mo> 
              <mrow> 
               <mi>
                 R 
               </mi> 
               <mo>
                 − 
               </mo> 
               <mi>
                 x 
               </mi> 
              </mrow> 
              <mo>
                ) 
              </mo> 
             </mrow> 
             <mtext>
               d 
             </mtext> 
             <mi>
               x 
             </mi> 
            </mrow> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mrow> 
                <mo>
                  ( 
                </mo> 
                <mrow> 
                 <msup> 
                  <mi>
                    R 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                 <mo>
                   − 
                 </mo> 
                 <msup> 
                  <mi>
                    x 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                </mrow> 
                <mo>
                  ) 
                </mo> 
               </mrow> 
               <mo>
                 + 
               </mo> 
               <msup> 
                <mrow> 
                 <mrow> 
                  <mo>
                    ( 
                  </mo> 
                  <mrow> 
                   <mi>
                     R 
                   </mi> 
                   <mo>
                     − 
                   </mo> 
                   <mi>
                     x 
                   </mi> 
                  </mrow> 
                  <mo>
                    ) 
                  </mo> 
                 </mrow> 
                </mrow> 
                <mn>
                  2 
                </mn> 
               </msup> 
              </mrow> 
             </msqrt> 
            </mrow> 
           </mfrac> 
          </mrow> 
         </mrow> 
        </mstyle> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(12)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        = 
      </mo> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mrow> 
       <mo>
         [ 
       </mo> 
       <mrow> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mtext>
             d 
           </mtext> 
           <mi>
             x 
           </mi> 
          </mrow> 
         </mrow> 
        </mstyle> 
        <mo>
          − 
        </mo> 
        <mstyle displaystyle="true"> 
         <mrow> 
          <msubsup> 
           <mo>
             ∫ 
           </mo> 
           <mrow> 
            <mo>
              − 
            </mo> 
            <mi>
              R 
            </mi> 
           </mrow> 
           <mi>
             R 
           </mi> 
          </msubsup> 
          <mrow> 
           <mfrac> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mi>
                 R 
               </mi> 
               <mo>
                 − 
               </mo> 
               <mi>
                 x 
               </mi> 
              </mrow> 
             </msqrt> 
             <mtext>
                 
             </mtext> 
             <mtext>
               d 
             </mtext> 
             <mi>
               x 
             </mi> 
            </mrow> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mn>
                 2 
               </mn> 
               <mi>
                 R 
               </mi> 
              </mrow> 
             </msqrt> 
            </mrow> 
           </mfrac> 
          </mrow> 
         </mrow> 
        </mstyle> 
       </mrow> 
       <mo>
         ] 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(13)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
        <mi>
          m 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           R 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math> (14)</p>
   <p>The equation for the ellipse is 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mrow> 
       <mrow> 
        <msup> 
         <mi>
           x 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           a 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        + 
      </mo> 
      <mrow> 
       <mrow> 
        <msup> 
         <mi>
           y 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           b 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        = 
      </mo> 
      <mn>
        1 
      </mn> 
     </mrow> 
    </math>, 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        a 
      </mi> 
      <mo>
        &gt; 
      </mo> 
      <mi>
        b 
      </mi> 
     </mrow> 
    </math>. Substituting 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       y 
     </mi> 
    </math> for 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       R 
     </mi> 
    </math> and 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math> for 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       d 
     </mi> 
    </math>, Equation (7) for the disks becomes,</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mstyle displaystyle="true"> 
       <mrow> 
        <msubsup> 
         <mo>
           ∫ 
         </mo> 
         <mrow> 
          <mo>
            − 
          </mo> 
          <mi>
            a 
          </mi> 
         </mrow> 
         <mi>
           a 
         </mi> 
        </msubsup> 
        <mrow> 
         <mrow> 
          <mo>
            [ 
          </mo> 
          <mrow> 
           <mn>
             1 
           </mn> 
           <mo>
             − 
           </mo> 
           <mfrac> 
            <mrow> 
             <mi>
               a 
             </mi> 
             <mo>
               − 
             </mo> 
             <mi>
               x 
             </mi> 
            </mrow> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mrow> 
                <mo>
                  ( 
                </mo> 
                <mrow> 
                 <msup> 
                  <mi>
                    b 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                 <mo>
                   − 
                 </mo> 
                 <mrow> 
                  <mo>
                    ( 
                  </mo> 
                  <mrow> 
                   <mrow> 
                    <mrow> 
                     <msup> 
                      <mi>
                        b 
                      </mi> 
                      <mn>
                        2 
                      </mn> 
                     </msup> 
                    </mrow> 
                    <mo>
                      / 
                    </mo> 
                    <mrow> 
                     <msup> 
                      <mi>
                        a 
                      </mi> 
                      <mn>
                        2 
                      </mn> 
                     </msup> 
                    </mrow> 
                   </mrow> 
                  </mrow> 
                  <mo>
                    ) 
                  </mo> 
                 </mrow> 
                 <msup> 
                  <mi>
                    x 
                  </mi> 
                  <mn>
                    2 
                  </mn> 
                 </msup> 
                </mrow> 
                <mo>
                  ) 
                </mo> 
               </mrow> 
               <mo>
                 + 
               </mo> 
               <msup> 
                <mrow> 
                 <mrow> 
                  <mo>
                    ( 
                  </mo> 
                  <mrow> 
                   <mi>
                     a 
                   </mi> 
                   <mo>
                     − 
                   </mo> 
                   <mi>
                     x 
                   </mi> 
                  </mrow> 
                  <mo>
                    ) 
                  </mo> 
                 </mrow> 
                </mrow> 
                <mn>
                  2 
                </mn> 
               </msup> 
              </mrow> 
             </msqrt> 
            </mrow> 
           </mfrac> 
          </mrow> 
          <mo>
            ] 
          </mo> 
         </mrow> 
         <mtext>
           d 
         </mtext> 
         <mi>
           x 
         </mi> 
        </mrow> 
       </mrow> 
      </mstyle> 
     </mrow> 
    </math>(15)</p>
   <p>or</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <mstyle displaystyle="true"> 
       <mrow> 
        <msubsup> 
         <mo>
           ∫ 
         </mo> 
         <mrow> 
          <mo>
            − 
          </mo> 
          <mi>
            a 
          </mi> 
         </mrow> 
         <mi>
           a 
         </mi> 
        </msubsup> 
        <mrow> 
         <mrow> 
          <mo>
            [ 
          </mo> 
          <mrow> 
           <mn>
             1 
           </mn> 
           <mo>
             − 
           </mo> 
           <mfrac> 
            <mrow> 
             <mi>
               a 
             </mi> 
             <mo>
               − 
             </mo> 
             <mi>
               x 
             </mi> 
             <mn> 
              <mo stretchy="false">
                ) 
              </mo> 
             </mn> 
            </mrow> 
            <mrow> 
             <msqrt> 
              <mrow> 
               <mi>
                 A 
               </mi> 
               <msup> 
                <mi>
                  x 
                </mi> 
                <mn>
                  2 
                </mn> 
               </msup> 
               <mo>
                 + 
               </mo> 
               <mi>
                 B 
               </mi> 
               <mi>
                 x 
               </mi> 
               <mo>
                 + 
               </mo> 
               <mi>
                 C 
               </mi> 
              </mrow> 
             </msqrt> 
            </mrow> 
           </mfrac> 
          </mrow> 
          <mo>
            ] 
          </mo> 
         </mrow> 
         <mtext>
           d 
         </mtext> 
         <mi>
           x 
         </mi> 
        </mrow> 
       </mrow> 
      </mstyle> 
     </mrow> 
    </math>(16)</p>
   <p>with</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        A 
      </mi> 
      <mo>
        = 
      </mo> 
      <mn>
        1 
      </mn> 
      <mo>
        − 
      </mo> 
      <mrow> 
       <mrow> 
        <msup> 
         <mi>
           b 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           a 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        , 
      </mo> 
      <mi>
        B 
      </mi> 
      <mo>
        = 
      </mo> 
      <mo>
        − 
      </mo> 
      <mn>
        2 
      </mn> 
      <mi>
        a 
      </mi> 
      <mo>
        , 
      </mo> 
      <mi>
        C 
      </mi> 
      <mo>
        = 
      </mo> 
      <msup> 
       <mi>
         a 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
      <mo>
        + 
      </mo> 
      <msup> 
       <mi>
         b 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
     </mrow> 
    </math></p>
   <p>The gravitational attraction on the surface of the ellipsoid at point 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math> becomes</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        G 
      </mi> 
      <mi>
        m 
      </mi> 
      <mi>
        σ 
      </mi> 
      <msubsup> 
       <mrow> 
        <mrow> 
         <mo>
           [ 
         </mo> 
         <mrow> 
          <mi>
            x 
          </mi> 
          <mo>
            − 
          </mo> 
          <mrow> 
           <mo>
             ( 
           </mo> 
           <mrow> 
            <mi>
              a 
            </mi> 
            <mo>
              + 
            </mo> 
            <mrow> 
             <mi>
               B 
             </mi> 
             <mo>
               / 
             </mo> 
             <mi>
               A 
             </mi> 
            </mrow> 
           </mrow> 
           <mo>
             ) 
           </mo> 
          </mrow> 
          <mi>
            E 
          </mi> 
          <mo>
            + 
          </mo> 
          <mrow> 
           <mi>
             D 
           </mi> 
           <mo>
             / 
           </mo> 
           <mi>
             A 
           </mi> 
          </mrow> 
         </mrow> 
         <mo>
           ] 
         </mo> 
        </mrow> 
       </mrow> 
       <mrow> 
        <mo>
          − 
        </mo> 
        <mi>
          a 
        </mi> 
       </mrow> 
       <mrow> 
        <mo>
          + 
        </mo> 
        <mi>
          a 
        </mi> 
       </mrow> 
      </msubsup> 
     </mrow> 
    </math>(17)</p>
   <p>where</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        D 
      </mi> 
      <mo>
        = 
      </mo> 
      <msqrt> 
       <mrow> 
        <mi>
          A 
        </mi> 
        <msup> 
         <mi>
           x 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
        <mo>
          + 
        </mo> 
        <mi>
          B 
        </mi> 
        <mi>
          x 
        </mi> 
        <mo>
          + 
        </mo> 
        <mi>
          C 
        </mi> 
       </mrow> 
      </msqrt> 
     </mrow> 
    </math></p>
   <p>and</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        E 
      </mi> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          ln 
        </mi> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mrow> 
          <mn>
            2 
          </mn> 
          <msqrt> 
           <mi>
             A 
           </mi> 
          </msqrt> 
          <mi>
            D 
          </mi> 
          <mo>
            + 
          </mo> 
          <mn>
            2 
          </mn> 
          <mi>
            A 
          </mi> 
          <mi>
            x 
          </mi> 
          <mo>
            + 
          </mo> 
          <mi>
            B 
          </mi> 
         </mrow> 
         <mo>
           ) 
         </mo> 
        </mrow> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msqrt> 
         <mi>
           A 
         </mi> 
        </msqrt> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math></p>
   <p>which can be programmed.</p>
  </sec><sec id="s5">
   <title>5. Roche Limit for Liquids</title>
   <p>The tidal force for a homogeneous liquid ellipsoid can now be calculated. An example using the parameters for the moon Pan, the closest moon that orbits Saturn, is plotted in <xref ref-type="fig" rid="fig3">
     Figure 3
    </xref>. The liquid ellipsoid is represented as a solid green line using Equation (18), where 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math> is the semi-major axis of the ellipsoid. The solid ellipsoid is represented as a solid blue line using Equation (18). The force for a rigid sphere is plotted as a dashed blue line using Equation (19), where 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       r 
     </mi> 
    </math> is the radius of the sphere.</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         F 
       </mi> 
       <mrow> 
        <mi>
          T 
        </mi> 
        <mi>
          i 
        </mi> 
        <mi>
          d 
        </mi> 
        <mi>
          a 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mrow> 
          <mrow> 
           <mo>
             ( 
           </mo> 
           <mrow> 
            <mi>
              d 
            </mi> 
            <mo>
              − 
            </mo> 
            <mi>
              a 
            </mi> 
           </mrow> 
           <mo>
             ) 
           </mo> 
          </mrow> 
         </mrow> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        − 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mi>
          d 
        </mi> 
        <mo>
          − 
        </mo> 
        <mi>
          a 
        </mi> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
      <msup> 
       <mi>
         w 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
      <mo>
        − 
      </mo> 
      <msub> 
       <mi>
         F 
       </mi> 
       <mrow> 
        <mi>
          e 
        </mi> 
        <mi>
          l 
        </mi> 
        <mi>
          l 
        </mi> 
        <mi>
          i 
        </mi> 
        <mi>
          p 
        </mi> 
        <mi>
          s 
        </mi> 
        <mi>
          o 
        </mi> 
        <mi>
          i 
        </mi> 
        <mi>
          d 
        </mi> 
       </mrow> 
      </msub> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         a 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(18)</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         F 
       </mi> 
       <mrow> 
        <mi>
          T 
        </mi> 
        <mi>
          i 
        </mi> 
        <mi>
          d 
        </mi> 
        <mi>
          a 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          M 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mrow> 
          <mrow> 
           <mo>
             ( 
           </mo> 
           <mrow> 
            <mi>
              d 
            </mi> 
            <mo>
              − 
            </mo> 
            <mi>
              r 
            </mi> 
           </mrow> 
           <mo>
             ) 
           </mo> 
          </mrow> 
         </mrow> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
      <mo>
        − 
      </mo> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mi>
          d 
        </mi> 
        <mo>
          − 
        </mo> 
        <mi>
          r 
        </mi> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
      <msup> 
       <mi>
         w 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
      <mo>
        − 
      </mo> 
      <mrow> 
       <mrow> 
        <mi>
          G 
        </mi> 
        <mi>
          m 
        </mi> 
       </mrow> 
       <mo>
         / 
       </mo> 
       <mrow> 
        <msup> 
         <mi>
           r 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mrow> 
     </mrow> 
    </math>(19)</p>
   <p>The plots below show the Roche limit is influenced by density and shape. The liquid ellipsoid (solid green) shows an almost exact match to the edge of the rings (~137,000 Km). This remarkable find shows the character of the moon when first placed in orbit inside the Roche limit as a liquid.</p>
   <fig id="fig3" position="float">
    <label>Figure 3</label>
    <caption>
     <title>Figure 3. Tidal force for different moon types using the parameters of the moon Pan.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId130.jpeg?20250113050638" />
   </fig>
   <p>The solid ellipsoid (solid blue) shows a Roche limit closer to the planet (~122,000 Km). This is probably unimportant. Larger mass moons did not cool inside the Roche limit, as none have large distortions.</p>
   <p>The Roche limit for the solid sphere is much closer to the planet. This would never cause the rings we see today. Only liquids match the ring pattern we see today.</p>
   <p>We now need a mechanism that creates free particles.</p>
  </sec><sec id="s6">
   <title>6. Free Particles</title>
   <p>A liquid ellipsoid and a spherically rigid moon are both held together by forces other than self gravity, i.e. chemical bonds. Once each cross their respective Roche limits, neither can retain objects on their surface. But neither moon has objects on the surface and as long as the moon itself is stable, nothing much should happen but a small change in shape. Mechanisms are needed creating free particles to participate in ring formation.</p>
   <p>The tidal force would be an instinctive choice. A distorted liquid may drip (see Tate’s Law). Or crust formed from cooling may drift off the ends, if bouyant. However, both methods reduce the mass of the moon, which reduces the tidal force at the surface. Neither mechanism participates to completion. This always results in a shepard moon and shepard moons are seldom seen in rings. Further, the type of particles produced do not match the particles seen in rings.</p>
   <p>Remembering the moon is a liquid, it becomes obvious the mechanism causing free particles has nothing to do with force, and everything to do with temperature hot liquids boil. Then rings are formed by the evaporation of a liquid moon in a circular orbit within the Roche limit. The efficiency of evaporation and the composition of the liquid determines the size of the particles placed in orbit.</p>
   <p>We generalize the idea that all objects in orbit around a planet begin as a liquid. Larger mass moons will migrate—losing particles inside the Roche limit, retaining particles once outside. Inside the limit, smaller mass moons evaporate completely.</p>
   <p>Near the limit, moons take an intermediate form. To demonstrate and support our argument, we explain the shape of the moon Pan (and the similar shape of Atlas).</p>
  </sec><sec id="s7">
   <title>7. The Shape of Pan</title>
   <p>The moon Pan has a curious shape. In <xref ref-type="fig" rid="fig4">
     Figure 4
    </xref>, Pan appears flattened with a ridge, or ring, around the equator. In contrast, most moons we see are only slightly distorted from spherical. We can explain the difference if we allow moons to begin as a liquid near the Roche limit. Note: as explained earlier, lower mass moons like Pan lose the ability to migrate.</p>
   <fig id="fig4" position="float">
    <label>Figure 4</label>
    <caption>
     <title>Figure 4. Pan: The innermost moon of Saturn.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId131.jpeg?20250113050639" />
   </fig>
   <p>Moons produce particles through evaporation (or boiling), but only lose particles near the poles along the semi-major axis (point a in our graphs). Tidal forces are greatest at the poles of the semi-major axis.</p>
   <p>Particles are retained along the circumference of the semi-minor axis (points at b in our graphs and our future ring). Points on the semi-minor axis are the same distance as the orbital path. Tidal forces are absent along the orbital path.</p>
   <p>Tidal forces effect the state of the matter. The lower pressure at 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math> prolongs evaporation and the higher pressure at 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       b 
     </mi> 
    </math> promotes condensation. When a crust begins to form, it produces a ring of circumference 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mn>
        2 
      </mn> 
      <mi>
        π 
      </mi> 
      <mi>
        b 
      </mi> 
     </mrow> 
    </math>. As mass (volume) loss continues (at 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       a 
     </mi> 
    </math>), crust develops on the inside of the ring (at 
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       b 
     </mi> 
    </math>). The crust is buoyant and should become self supportive (structurally) at some point. The change to the ellipsoid is diagrammed in <xref ref-type="fig" rid="fig5">
     Figure 5
    </xref>.</p>
   <fig id="fig5" position="float">
    <label>Figure 5</label>
    <caption>
     <title>Figure 5. (Left) Ring of crust starts to form, 

      <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
  
        <mn>
         
   2
  
        </mn>
  
        <mi>
         
   π
  
        </mi>
  
        <mi>
         
   b
  
        </mi>
 
       </mrow>

      </math>. (Right) At a later time, the ring/disc increases from the inside as the volume of ellipsoid decreases.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId142.jpeg?20250113050639" />
   </fig>
   <p>The next progression of changes is shown in <xref ref-type="fig" rid="fig6">
     Figure 6
    </xref>. The volume of the ellipsoid decreases and the size of the disk increases until a change in moment of inertia produces a 90˚ rotation about some axis in the plane of the disk. This rotation is produced by a force [moment] from the planet acting on some line in the plane of the disk, perpendicular to the axis of rotation. We chose a line through B. The result, seen in the center, is followed by a second rotation about B, aligning the disk in the plane of the orbit. This torque is represented by the arrows at the ends of the disk. Finally, the compression of the semi-major axis, still liquid, is shown on the right (and see <xref ref-type="fig" rid="fig4">
     Figure 4
    </xref>).</p>
   <fig id="fig6" position="float">
    <label>Figure 6</label>
    <caption>
     <title>Figure 6. Left) Ellipsoid showing disk before rotation. Center) 90˚ rotation of B about center. Arrows designate torque on disk from the planet about B. Right) Disk lies in plane of orbit. Semi-major axis compresses.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/2181202-rId145.jpeg?20250113050639" />
   </fig>
  </sec><sec id="s8">
   <title>8. Conclusion</title>
   <p>We suggest liquid moons boil when placed in orbit—remaining a moon outside the limit, evaporate into a ring inside the limit and take an indeterminate shape at the limit. Pan is just inside the Roche limit, but far enough to lose mass. Slow mass loss while cooling guides the shape.</p>
   <p>In space, liquids are a forgotten state of matter. It is hoped the arguments we present show the possibilities and encourage participation.</p>
  </sec>
 </body><back>
  <ref-list>
   <title>References</title>
   <ref id="scirp.138767-ref1">
    <label>1</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Charnoz, S., Canup, R.M., Crida, A. and Dones, L. (2017) The Origin of Planetary Ring Systems. &gt;https://arxiv.org/pdf/1703.09741
    </mixed-citation>
   </ref>
   <ref id="scirp.138767-ref2">
    <label>2</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Fitzpatrick, R. (2016) Roche Ellipsoids. &gt;https://farside.ph.utexas.edu/teaching/336L/Fluid/node37.html
    </mixed-citation>
   </ref>
   <ref id="scirp.138767-ref3">
    <label>3</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Thomas, P.C. and Helfenstein, P. (2020) The Small Inner Satellites of Saturn: Shapes, Structures and Some Implications. Icarus, 344, Article ID: 113355. &gt;https://doi.org/10.1016/j.icarus.2019.06.016
    </mixed-citation>
   </ref>
   <ref id="scirp.138767-ref4">
    <label>4</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Thomas, P.C. (2010) Sizes, Shapes, and Derived Properties of the Saturnian Satellites after the Cassini Nominal Mission. Icarus, 208, 395-401. &gt;https://doi.org/10.1016/j.icarus.2010.01.025
    </mixed-citation>
   </ref>
  </ref-list>
 </back>
</article>