<p>The <em>big essential fact</em> about momentum is that it is conserved (in the presence of a central force, as is the case here). So the Angular momentum of Pluto today is the same as it was yesterday, and the same as last year and (excepting perturbations) the same as it has ever been.</p>
<p>It is easiest to calculate for a body that is moving perpendicular to its position vector. This is always true for circular motion, about the centre of the circle.  It is not true for elliptical motion, except at apoapsis and periapsis.</p>
<p>At these times <span class="math-container">$L = mvr$</span>.  For Pluto the periapsis speed is <span class="math-container">$v = 6.10\, km/s$</span> the distance is 4.44 billion km and the mass is <span class="math-container">$1.31 × 10^{22}$</span> kg.  To get the angular momentum you've got to multiply them together. If you want SI units, convert those km to m first.</p>
<p>The angular momentum today is the same.</p>
<p>Alternately you can use the relationship</p>
<p><span class="math-container">$L= \sqrt{\mu \ell}$</span> where <span class="math-container">$\mu = GM =1.33×10^{20}$</span> and <span class="math-container">$\ell$</span> is the semi latus rectum or <span class="math-container">$\ell=a(1+e^2)$</span>, and you have to plug in the semimajor axis for Pluto and eccentricity of its orbit.</p>
<p>Or you can look it up, since it is constant <span class="math-container">$L=3.6\times 10^{38}$</span></p>

CC BY-SA 4.0
2021-02-19T02:17:26.113Z
How important are satellites in calculation of angular momentum of a planet? Do the 79 satellites of Jupiter change/affect angular momentum of Jupiter?
Majoris
CC BY-SA 4.0
2021-01-27T21:07:09.650Z
2021-01-27T21:28:47.557Z
2021-01-27T21:28:47.557Z
James K
<blockquote>
<p>Where can I get the actual data to calculate this?</p>
</blockquote>
<p>If you are a programmer, you could use the SPICE library: <a href="https://naif.jpl.nasa.gov/naif/toolkit.html" rel="nofollow noreferrer">https://naif.jpl.nasa.gov/naif/toolkit.html</a> and the data files available here: <a href="https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/" rel="nofollow noreferrer">https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/</a> or you could use <a href="https://wgc.jpl.nasa.gov:8443/webgeocalc/#StateVector" rel="nofollow noreferrer">https://wgc.jpl.nasa.gov:8443/webgeocalc/#StateVector</a> or HORIZONS: <a href="https://ssd.jpl.nasa.gov/horizons.cgi" rel="nofollow noreferrer">https://ssd.jpl.nasa.gov/horizons.cgi</a></p>
<p>The masses are listed here: <a href="https://naif.jpl.nasa.gov/pub/naif/generic_kernels/pck/" rel="nofollow noreferrer">https://naif.jpl.nasa.gov/pub/naif/generic_kernels/pck/</a> file &quot;gm_de431.tpc&quot;.</p>
<p>For example, the state of Pluto (<strong>here I'm using Pluto system: Pluto + satellites</strong>) on 2021-01-27 00:00 UTC in the mean ecliptic and equinox of J2000 reference frame is:<br></p>
<p><strong>r</strong> = (2.11488 · 10<sup>12</sup>, -4.65878 · 10<sup>12</sup>, -1.13045 · 10<sup>11</sup>) m<br>
<strong>v</strong> = (5088.6; 1090.04; -1584.49) m/s<br></p>
<p>GM = 977 · 10<sup>9</sup> m<sup>3</sup>/s<sup>2</sup>. If we use G = 6.6743 · 10<sup>-11</sup> m<sup>3</sup> kg<sup>-1</sup> s<sup>-2</sup> (<a href="https://physics.nist.gov/cgi-bin/cuu/Value?bg" rel="nofollow noreferrer">https://physics.nist.gov/cgi-bin/cuu/Value?bg</a>), M= 1.4638239 · 10<sup>22</sup> kg.</p>
<p><strong>L</strong> = <strong>r</strong> x M<strong>v</strong> = (1.09861 · 10<sup>38</sup>, 4.06324 · 10<sup>37</sup>, 3.80771 · 10<sup>38</sup>) kgm<sup>2</sup>/s, the magnitude is 3.9838 · 10<sup>38</sup> kgm<sup>2</sup>/s.</p>
<p>If we do the calculations for 20 years ago (2001-01-27 00:00 UTC) we get:</p>
<p><strong>L</strong> = (1.09685 · 10<sup>38</sup>, 4.04601 · 10<sup>37</sup>, 3.78485 · 10<sup>38</sup>) kgm<sup>2</sup>/s, the magnitude is 3.9613 · 10<sup>38</sup>  kgm<sup>2</sup>/s.</p>
<p><strong>EDIT</strong></p>
<p>For Mars, 1970-08-25 00:00 UTC (the notation 123e+011 means 123 · 10<sup>11</sup>):</p>
<p><strong>r</strong> = (-1.9051451e+011; 1.58990113e+011; 8.02428218e+009) m<br>
<strong>v</strong> = (-14604.5348; -16539.1841; 13.5952337) m/s<br>
<strong>L</strong> = (8.65493375e+037; -7.35385224e+037; 3.51194063e+039) kgm<sup>2</sup>/s<br>
||<strong>L</strong>|| = 3.51377655e+039 kgm<sup>2</sup>/s.</p>
<p>For today (Jan 28) I get:</p>
<p><strong>r</strong> = (4.11669129e+010; 2.27295989e+011; 3.75334592e+009) m<br>
<strong>v</strong> = (-22924.868; 6375.0766; 695.965252) m/s<br>
<strong>L</strong> = (8.61551162e+037; -7.35994071e+037; 3.51209564e+039) kgm<sup>2</sup>/s<br>
||L|| = 3.51392307e+039 kgm<sup>2</sup>/s.</p>
<p><strong>WARNING!</strong> Do not to assume that all the digits are significant. I wrote the numbers with an unrealistic accuracy just to verify the calculations.</p>
<hr />
<p><strong>EDIT #2</strong></p>
<p>It could be interesting to see the variation of <em>L</em> against the time.<br>
<br>
The following graph shows the value of <em>(L - Lmin) / Lmin</em>, where <em>Lmin</em> is the smallest calculated angular momentum and <em>L</em> is the instantaneous angular momentum.<br>
The calculations are done with the SPICE library and the DE430, jup310 and mar097 ephemerides. The independent variable is the time wrt the orbital period of each planet. The starting epoch is 2021-01-28 00:00 UTC.</p>
<p><a href="https://i.sstatic.net/XknyA.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/XknyA.png" alt="enter image description here" /></a></p>
<p>The graph is as expected. Mercury, well inside the Sun gravity well, is only slightly perturbed by the other bodies.<br>
The Earth is heavily perturbed by the Moon and also Jupiter is significantly perturbed by its moons.</p>

CC BY-SA 4.0
2021-01-28T15:53:47.027Z
A BIG thank you! This helps so much. I loved the State Vector calculator. :) Would you mind editing your answer with one more example for another planet and another random date, let's say for Mars on date 25 Aug 1970.
Majoris
CC BY-SA 4.0
2021-01-28T18:36:22.190Z
Oh, that's fair enough. Even Horizons does that in their generated ephemeris data, and they tell you not to assume all the digits are significant.
PM 2Ring
CC BY-SA 4.0
2021-01-28T19:30:57.750Z
@PM2Ring I wrote the numbers with unrealistic accuracy, so Majoris can verify his calculations. For Pluto, I already wrote "here I'm using Pluto system: Pluto + satellites".
Cristiano
CC BY-SA 4.0
2021-01-30T20:05:38.917Z
@Cristiano Awesome! That graph is what I want to do. We should connect. How do I contact you?
Majoris
CC BY-SA 4.0
2021-01-30T21:03:18.920Z
@Majoris If it's not a marriage proposal, I guess you can write here. :-) What's the problem?
Cristiano
CC BY-SA 4.0
2021-01-28T10:33:45.553Z
2021-01-30T11:15:19.740Z
2021-01-30T11:15:19.740Z
Cristiano
<p>I want to calculate the Angular Momentum of any planet at any point on the orbit around the Sun.
Example: I want to calculate the Angular Momentum of Pluto today, on Jan 27, 2021 00:00:00 Hrs</p>
<p>Where can I get the actual data to calculate this?
I will need Mass, Velocity and Radius to calculate. Do I need the θ too?
Is it readily calculated by NASA?</p>
<p><a href="http://hyperphysics.phy-astr.gsu.edu/hbase/amom.html" rel="nofollow noreferrer">http://hyperphysics.phy-astr.gsu.edu/hbase/amom.html</a></p>
<p><a href="https://i.sstatic.net/iEWYi.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/iEWYi.jpg" alt="enter image description here" /></a>
<a href="https://i.sstatic.net/NTiOC.gif" rel="nofollow noreferrer"><img src="https://i.sstatic.net/NTiOC.gif" alt="enter image description here" /></a></p>

CC BY-SA 4.0
2021-01-27T22:18:53.843Z
Why do you think the angular momentum will change?!
ProfRob
CC BY-SA 4.0
2021-01-28T15:13:27.077Z
One could want to calculate the angular momentum at many points in time to demonstrate that it does _not_ change.
Phil Frost
CC BY-SA 4.0
2021-01-28T17:09:37.983Z
@PhilFrost To demonstrate that it varies, I'd say.
Cristiano
<orbit><planet><solar-system><orbital-mechanics>
Calculating the Angular Momentum of a planet
Majoris
# Calculating the Angular Momentum of a planet

(Asked by: username_0 on 2021-01-27)

I want to calculate the Angular Momentum of any planet at any point on the orbit around the Sun.
Example: I want to calculate the Angular Momentum of Pluto today, on Jan 27, 2021 00:00:00 Hrs

Where can I get the actual data to calculate this?
I will need Mass, Velocity and Radius to calculate. Do I need the θ too?
Is it readily calculated by NASA?

<http://hyperphysics.phy-astr.gsu.edu/hbase/amom.html>

[![enter image description here](https://i.sstatic.net/iEWYi.jpg)](https://i.sstatic.net/iEWYi.jpg)
[![enter image description here](https://i.sstatic.net/NTiOC.gif)](https://i.sstatic.net/NTiOC.gif)

**Tags:** orbit, planet, solar-system, orbital-mechanics
> **Comments:**
> **username_1**: Why do you think the angular momentum will change?!  
> **username_2**: One could want to calculate the angular momentum at many points in time to demonstrate that it does _not_ change.  
> **username_3**: @PhilFrost To demonstrate that it varies, I'd say.  

---

## 2 answers

### Answer by username_3 on 2021-01-28. Score: 5 | Accepted answer

> Where can I get the actual data to calculate this?

If you are a programmer, you could use the SPICE library: <https://naif.jpl.nasa.gov/naif/toolkit.html> and the data files available here: <https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/> or you could use <https://wgc.jpl.nasa.gov:8443/webgeocalc/#StateVector> or HORIZONS: <https://ssd.jpl.nasa.gov/horizons.cgi>

The masses are listed here: <https://naif.jpl.nasa.gov/pub/naif/generic_kernels/pck/> file "gm\_de431.tpc".

For example, the state of Pluto (**here I'm using Pluto system: Pluto + satellites**) on 2021-01-27 00:00 UTC in the mean ecliptic and equinox of J2000 reference frame is:

**r** = (2.11488 · 1012, -4.65878 · 1012, -1.13045 · 1011) m  
**v** = (5088.6; 1090.04; -1584.49) m/s

GM = 977 · 109 m3/s2. If we use G = 6.6743 · 10-11 m3 kg-1 s-2 (<https://physics.nist.gov/cgi-bin/cuu/Value?bg>), M= 1.4638239 · 1022 kg.

**L** = **r** x M**v** = (1.09861 · 1038, 4.06324 · 1037, 3.80771 · 1038) kgm2/s, the magnitude is 3.9838 · 1038 kgm2/s.

If we do the calculations for 20 years ago (2001-01-27 00:00 UTC) we get:

**L** = (1.09685 · 1038, 4.04601 · 1037, 3.78485 · 1038) kgm2/s, the magnitude is 3.9613 · 1038 kgm2/s.

**EDIT**

For Mars, 1970-08-25 00:00 UTC (the notation 123e+011 means 123 · 1011):

**r** = (-1.9051451e+011; 1.58990113e+011; 8.02428218e+009) m  
**v** = (-14604.5348; -16539.1841; 13.5952337) m/s  
**L** = (8.65493375e+037; -7.35385224e+037; 3.51194063e+039) kgm2/s  
||**L**|| = 3.51377655e+039 kgm2/s.

For today (Jan 28) I get:

**r** = (4.11669129e+010; 2.27295989e+011; 3.75334592e+009) m  
**v** = (-22924.868; 6375.0766; 695.965252) m/s  
**L** = (8.61551162e+037; -7.35994071e+037; 3.51209564e+039) kgm2/s  
||L|| = 3.51392307e+039 kgm2/s.

**WARNING!** Do not to assume that all the digits are significant. I wrote the numbers with an unrealistic accuracy just to verify the calculations.

---

**EDIT #2**

It could be interesting to see the variation of *L* against the time.  
  
The following graph shows the value of *(L - Lmin) / Lmin*, where *Lmin* is the smallest calculated angular momentum and *L* is the instantaneous angular momentum.  
The calculations are done with the SPICE library and the DE430, jup310 and mar097 ephemerides. The independent variable is the time wrt the orbital period of each planet. The starting epoch is 2021-01-28 00:00 UTC.

[![enter image description here](https://i.sstatic.net/XknyA.png)](https://i.sstatic.net/XknyA.png)

The graph is as expected. Mercury, well inside the Sun gravity well, is only slightly perturbed by the other bodies.  
The Earth is heavily perturbed by the Moon and also Jupiter is significantly perturbed by its moons.

> **Comments:**
> **username_0**: A BIG thank you! This helps so much. I loved the State Vector calculator. :) Would you mind editing your answer with one more example for another planet and another random date, let's say for Mars on date 25 Aug 1970.  
> **username_5**: Oh, that's fair enough. Even Horizons does that in their generated ephemeris data, and they tell you not to assume all the digits are significant.  
> **username_3**: @PM2Ring I wrote the numbers with unrealistic accuracy, so Majoris can verify his calculations. For Pluto, I already wrote "here I'm using Pluto system: Pluto + satellites".  
> **username_0**: @Cristiano Awesome! That graph is what I want to do. We should connect. How do I contact you?  
> **username_3**: @Majoris If it's not a marriage proposal, I guess you can write here. :-) What's the problem?  

---

### Answer by username_4 on 2021-01-27. Score: 11

The *big essential fact* about momentum is that it is conserved (in the presence of a central force, as is the case here). So the Angular momentum of Pluto today is the same as it was yesterday, and the same as last year and (excepting perturbations) the same as it has ever been.

It is easiest to calculate for a body that is moving perpendicular to its position vector. This is always true for circular motion, about the centre of the circle. It is not true for elliptical motion, except at apoapsis and periapsis.

At these times $L = mvr$. For Pluto the periapsis speed is $v = 6.10\, km/s$ the distance is 4.44 billion km and the mass is $1.31 × 10^{22}$ kg. To get the angular momentum you've got to multiply them together. If you want SI units, convert those km to m first.

The angular momentum today is the same.

Alternately you can use the relationship

$L= \sqrt{\mu \ell}$ where $\mu = GM =1.33×10^{20}$ and $\ell$ is the semi latus rectum or $\ell=a(1+e^2)$, and you have to plug in the semimajor axis for Pluto and eccentricity of its orbit.

Or you can look it up, since it is constant $L=3.6\times 10^{38}$

> **Comments:**
> **username_0**: How important are satellites in calculation of angular momentum of a planet? Do the 79 satellites of Jupiter change/affect angular momentum of Jupiter?  

---

