<h1> Short answer </h1>
<p>Short answer is no, it can't prove or disprove the existence of Planet 9. Even if there is a significant difference between the barycenter of the solar system with and without Planet 9 we wouldn't be able to tell without hundreds if not thousands of years of precise data. If we don't have Planet 9 as a reference and are limited to the Sun and the eight known planets for reference then without another reference in the same frame we have no way to observe the motion of our reference.</p>
<p>It is perhaps useful to consider the solar system (i.e. the Sun plus the eight known planets) as orbiting <em>together</em> around the barycenter of the solar system plus Planet 9 system. Without seeing Planet 9 as a reference how would we have any notion that this other barycenter even exists?</p>
<p>Of course Planet 9 would still be exerting gravitational influence and over time this would cause a deviation between where the barycenter <em>should</em> be and where it actually is -- that is the only way we could use the motion of the barycenter to infer anything about individual bodies. Essentially, the barycenter is a convolution of the gravitational influence of all the mass,<sup>5</sup> its deconvolution requires analysis of temporal variation.</p>
<hr />
<h1> Unnecessarily long answer </h1>
<p><strong>Note</strong>   -   When I discuss variations in the barycenter I will be referring to the barycenter of the solar system in heliocentric, not barycentric, coordinates (in barycentric coordinates of course there is no variation in the barycenter, by definition). This simplifies some of the qualitative analysis.</p>
<hr />
<p>Rigorously answering this requires answering several separate questions (some which aren't necessarily relevant in answering the broader question but do provide useful background):</p>
<ol>
<li>Does the location of the barycenter over time contain enough information to recover information about the contributing masses?</li>
<li>How precisely can we predict the location of the barycenter?</li>
<li>How much deviation would a <span class="math-container">$5$</span>-<span class="math-container">$10~M$</span><sub>⊕</sub> planet at <span class="math-container">$400$</span>-<span class="math-container">$800$</span> AU cause to the center of mass of the solar system?</li>
<li>If no such deviation exists, does that rule out the existence of Planet 9? Conversely, if such a deviation does exist, does it rule <strong>in</strong> the existence of Planet 9?</li>
</ol>
<hr />
<h3> 1 - Location of the barycenter over time </h3>
<p>The short answer to question <strong>1</strong> is yes (sort of). There is information on the orbits of the contributing masses which can be backed out from simple analysis of the location of the barycenter <strong>over time</strong>.</p>
<p>The emphasis here is because at any given point in time the location of the barycenter contains no information about the location of individual bodies, only the entire mass distribution. However, since most of the mass in the solar system is discretized into large bodies (as opposed to a gas or dust cloud, protoplanetary disk, etc), the influence of each massive body on the location of the barycenter will vary with a frequency proportional to its orbital period -- and, critically, this is fairly measurable (but only if we can also know the relative positions of the bodies). Consequently, the influence of the individual masses can be retrieved (this is not necessarily the case for gas or dust clouds and protoplanetary disks).</p>
<p>In a simple system, retrieving information about the orbits is fairly straightforward (although I will neglect the rigorous demonstration for brevity). The issues arise when the system becomes more complicated.</p>
<p>One source of confound is resonant orbits. Certain resonances would defy deconvolution because the frequencies of their respective influences on the barycenter are synchronous. Furthermore, multiple resonances can't necessarily be disambiguated: the frequency pattern in the motion of the barycenter caused by 2 bodies in a resonance could be replicated by 3 bodies in a resonance. Pluto and Neptune are in a <a href="https://en.wikipedia.org/wiki/Orbital_resonance#Mean-motion_resonances_in_the_Solar_System" rel="nofollow noreferrer">2:3 resonance</a>, so I don't think we can immediately dismiss the possibility of a theoretical Planet 9 being in a resonance of some kind and that this could hinder deconvolving its effect on the barycenter.</p>
<p>Another significant challenge is the qualification that the barycenter must be observed <strong>over time</strong>, and more specifically how much time is needed to deconvolve the various influences. As I mentioned, each body's influence on the barycenter is cyclic based on their orbital period and eccentricity. Using the Bulirsch-Stoer integrator from <a href="https://github.com/WMiller256/Vulcan" rel="nofollow noreferrer">Vulcan</a> and tracking an energy analog of the barycenter's motion<sup>1</sup> as a function of time we can analyze the effect of the different orbital frequencies.</p>
<p>Let's consider a three body system of the Sun, Saturn, and Jupiter. In this system the influence of Saturn and Jupiter should clearly be visible in the resulting sinusoid. After integrating this system for <span class="math-container">$2$</span> Saturn years (<span class="math-container">$60$</span> Earth years) the energy-analogy over time looks like this:</p>
<p><a href="https://i.sstatic.net/CweUg.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/CweUg.png" alt="enter image description here" /></a></p>
<p>It should be pretty clear that the above behavior can be created by two sinusoids with different frequencies, to any who find this unclear I recommend having a look at <a href="https://demonstrations.wolfram.com/CombiningSineWaves/" rel="nofollow noreferrer">this tool</a>. If wanted to we could use Fourier analysis to deconvolve those two frequencies and recover the orbital periods of Jupiter and Saturn (I will also neglect this for brevity).</p>
<p>Let's consider a slightly worse case: a three body system comprised of the Sun, Jupiter, and a hypothetical <a href="https://en.wikipedia.org/wiki/Hot_Jupiter#:%7E:text=Hot%20Jupiters%20are%20a%20class,the%20moniker%20%22hot%20Jupiters%22." rel="nofollow noreferrer">hot jupiter</a> with a semi-major axis of <span class="math-container">$1.55^{}$</span>x<span class="math-container">$10^{10}$</span> m on a relatively circular orbit. Our energy-analog after <span class="math-container">$3$</span> hot-jupiter years (<span class="math-container">$12$</span> Earth days) looks like this:</p>
<p><a href="https://i.sstatic.net/Dabr5.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/Dabr5.png" alt="enter image description here" /></a></p>
<p><strong>Where did the other sinusoid go?</strong></p>
<p>Its effect is still there, but we haven't sampled sufficient time to see it. We would need to look at <span class="math-container">$1,095$</span> hot-jupiter years (<span class="math-container">$12$</span> Earth years) in order to see one orbit of our toy Jupiter<sup>2</sup> Fortunately we don't have to integrate quite that long to see the effect of the toy Jupiter on the energy-analog. After <span class="math-container">$320$</span> hot-jupiter years (<span class="math-container">$3.5$</span> Earth years):</p>
<p><a href="https://i.sstatic.net/aYz3A.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/aYz3A.png" alt="enter image description here" /></a></p>
<p>This becomes even more complicated with non-zero eccentricities and unstable orbits. For a similar three body, hot-jupiter system but with the toy Jupiter having an eccentricity of <span class="math-container">$0.65$</span> over a period of <span class="math-container">$7$</span> Earth years the energy-analog is:<sup>3</sup></p>
<p><a href="https://i.sstatic.net/dbdUa.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/dbdUa.png" alt="enter image description here" /></a></p>
<p>There are other ways to draw inferences about the orbital parameters from the relative motion of the barycenter (e.g. the shape of individual cycles), but that is for a different question.</p>
<p><strong>In sum, the location of the barycenter over time could certainly be used to find a mass imbalance, but but not necessarily argue for the existence of Planet 9 -- that would depend on the specifics of the observations.</strong></p>
<hr />
<h2> 2 - Precision of barycenter predictions </h2>
<p>There are number of significant barriers to accurately calculating the barycenter of the solar system but the most troublesome are the uncertainties regarding the interiors of Saturn and Jupiter. Specifically, the behavior of liquid metallic hydrogen at such massive pressures and (more so for Saturn) understanding their gravitational moments (Fortney <a href="https://zenodo.org/record/1230846#.Xw02XihKiiN" rel="nofollow noreferrer">2004</a>).</p>
<p>The issue this raises is that without knowing the centers of mass for Jupiter and Saturn (which contain 92% of planetary mass in the solar system) with sufficient precision we don't know well enough what the solar system barycenter <em>should</em> be to determine whether or not the true barycenter differs enough to indicate the existence of Planet 9.</p>
<p>With information from the Cassini spacecraft combined with data from the <a href="https://public.nrao.edu/telescopes/vlba/" rel="nofollow noreferrer">VLBA radio telescope</a> the Saturn ephemerides was constrained to <span class="math-container">$4$</span> km in <a href="http://nasa.gov/jpl/cassini/scientists-pinpoint-saturn-with-exquisite-accuracy" rel="nofollow noreferrer">2015</a>. This represented an improvement of about an order of magnitude. The Juno spacecraft provided a similar improvement to the Jupiter ephemerides, constraining it to within <span class="math-container">$10$</span> km in <a href="https://ui.adsabs.harvard.edu/abs/2019AAS...23330201J/abstract" rel="nofollow noreferrer">2019</a>.</p>
<p>What is less well constrained is how constant these parameters are, so it bears noting that an inaccuracy of <span class="math-container">$350$</span> km in the location of the center of mass of Saturn's core within the planet corresponds to an uncertainty of <span class="math-container">$100$</span> m in the position of the solar system barycenter. Jupiter lacks a dense core like Saturn's, so it is more difficult to quantify a direct relationship between inaccuracies, but suffice it to say that an inaccuracy of <span class="math-container">$100$</span> km in Jupiter's center of mass corresponds to an uncertainty of <span class="math-container">$100$</span> m in the barycenter.<sup>4</sup></p>
<p>So assuming Saturn's center of mass has not drifted much since September of 2018, <strong>I would estimate the precision of the calculated solar system barycenter to be roughly <span class="math-container">$\pm15$</span> m.</strong></p>
<hr />
<h2> 3 - Deviation from predicted barycenter </h2>
<p>Calculating the effect Planet 9 would have on the barycenter is actually quite trivial thanks to superposition (see <a href="https://stackoverflow.com/questions/59768976/how-can-i-work-out-the-gravitational-force-of-any-2d-polygon/59797041#59797041">this answer</a> for a more rigorous treatment of the limitations of this concept when it comes to gravitational interactions). The equation for calculating the distance from the primary body to the barycenter in a two-body system,</p>
<p><span class="math-container">$$
r_1 = \frac{a}{1 + \frac{m_1}{m_2}}
$$</span></p>
<p>can also be applied to calculating the barycenter of two barycenters where <span class="math-container">$a$</span> is the distance between the two barycenters, <span class="math-container">$m_1$</span> is the mass which contributes to the primary barycenter, <span class="math-container">$m_2$</span> is the mass which contributes to the secondary barycenter, and <span class="math-container">$r_1$</span> is the distance from the primary barycenter to the shared barycenter. And since <span class="math-container">$99.86$</span>% of the mass of the solar system is contained in the sun it is a rigorous approximation to use</p>
<p><span class="math-container">$$
r_1 = \frac{200 \cdot 1.496 \text{x} 10^{11}}
           {1 + \frac{1.989e30}{5 \cdot 5.972 \text{x} 10^{24}}} 
    = 4 \text{x} 10^8 \text{ m}
$$</span></p>
<p>as a lower bound, and</p>
<p><span class="math-container">$$
r_1 = \frac{1200 \cdot 1.496 \text{x} 10^{11}}
           {1 + \frac{1.989e30}{10 \cdot 5.972 \text{x} 10^{24}}} 
    = 5 \text{x} 10^9 \text{ m}
$$</span></p>
<p>as an upper bound for the distance from the barycenter of the solar system without Planet 9 to the barycenter with it.</p>
<hr />
<h2> 4 - Interpretation </h2>
<p>Even the lower bound of <span class="math-container">$4\text{x}10^{8}$</span> m is a huge deviation, so did we just demonstrate that Planet 9 can't exist? Unfortunately no. This goes back to number <strong>1</strong> - the problem is time. There are a number of ways to think about this, but the salient point is that the influence on the barycenter of the solar system would have been persistent since well before the earliest observational records. Planet 9 would have to have an orbital period on the order of <span class="math-container">$10,000$</span> years. Over periods of thousands of years it would have obvious and readily measurable effects on the motion of the other planets (and that motion would be reflected in the motion of the heliocentric barycenter), but otherwise the effects are simply too small to disambiguate from noise.</p>
<hr />
<p><sup>1</sup> I create this analog by calculating a hypothetical energy of the system if the center of the sun was the barycenter of the system and comparing it to the initial value. It is not tremendously physically meaningful but it does reduce the parameter-space of the temporal position of the barycenter to two dimensions, thereby making it easier to show the effect I discuss graphically.</p>
<p><sup>2</sup> This is one of the reasons simulating hot jupiters is so difficult (and one of the things the Vulcan project is intended to remedy): the limiting factor for orbital integration is always the fastest-orbiting body (although multi-step integrators mitigate this somewhat) so systems with large differences between orbital periods take proportionally longer to simulate. Using the <code>Mercury</code> code, a gigayear solar system simulation with a hot-jupiter on an 8 day orbit would require roughly 10 months of simulation time on a high-end desktop computer.</p>
<p><sup>3</sup> Note that the energy analog is degenerate here because the hot-jupiter will be ejected after approximately 10 Earth years.</p>
<p><sup>4</sup> These uncertainty relationships are from my own calculations and may be flawed as they have not been rigorously tested or reviewed.</p>
<p><sup>5</sup> And I do literally mean <strong>all</strong> the mass. As in all the mass in the universe, gravity has an infinite range after all....</p>
<hr />

CC BY-SA 4.0
2020-07-14T08:52:14.570Z
I realize the plots are lacking axis titles, I intend to remedy that but have to modify the Vulcan source to do so.
William Miller
CC BY-SA 4.0
2020-10-07T17:20:08.707Z
I've taken the liberty of backing up the current state of Vulcan to archive.org to have some host redundancy to prevent link rot (I credited it to you): https://archive.org/details/vulcan-master
called2voyage
CC BY-SA 4.0
2020-07-14T08:22:24.650Z
2022-07-11T03:08:53.923Z
2022-07-11T03:08:53.923Z
William Miller
<p>I've just read this recent news article, <a href="https://www.sciencealert.com/astronomers-have-located-the-centre-of-the-solar-system-to-within-100-metres" rel="noreferrer">Astronomers Have Located The Centre of The Solar System to Within 100 Metres</a>, reporting on a paper[1] that claims narrowing the Solar System barycentre to within 100 meters, just outside the sun, and I wonder, could this result help prove or disprove the hypothesis of a massive <a href="https://en.wikipedia.org/wiki/Planet_Nine" rel="noreferrer">Planet Nine</a> as proposed by Batygin and Brown[2] in 2016?</p>
<p>I imagine that any unknown mass concentration would change the barycentre and generate a error, and if we could measure this error, it should be possible to pinpoint the position of said mass concentration. Can we just use probe tracking data from the <a href="https://en.wikipedia.org/wiki/Voyager_program" rel="noreferrer">Voyagers</a> or <a href="https://en.wikipedia.org/wiki/Pioneer_10" rel="noreferrer">Pioneers</a> to see if their hyperbolic trajectories match the barycentre reported in the paper, for example?</p>
<p>[1]. <a href="https://arxiv.org/pdf/2001.00595.pdf" rel="noreferrer">Vallisneri, M., et al. “Modeling the Uncertainties of Solar System Ephemerides for Robust Gravitational-wave Searches with Pulsar-timing Arrays”. The Astrophysical Journal, vol. 893, no 2, abril de 2020, p. 112. DOI.org (Crossref), doi:10.3847/1538-4357/ab7b67.</a></p>
<p>[2]. <a href="https://arxiv.org/abs/1601.05438" rel="noreferrer">Batygin, Konstantin, e Michael E. Brown. “Evidence for a Distant Giant Planet in the Solar System”. The Astronomical Journal, vol. 151, no 2, janeiro de 2016, p. 22. arXiv.org, doi:10.3847/0004-6256/151/2/22.</a></p>

CC BY-SA 4.0
2020-07-12T12:05:40.853Z
Here's a recent press release about an alternative explanation for the oddities of the outer Solar System that doesn't require Planet Nine: https://www.colorado.edu/today/2020/07/07/outer-solar-system
PM 2Ring
CC BY-SA 4.0
2020-07-12T15:28:20.847Z
FWIW, a rough calculation of $rm/M_\odot$ for Planet 9, using $r=400$ AU & $m=5M_\oplus=15×10^{-6}M_\odot$, gives around 900,000 km. That's a pretty large contribution to the barycentre. But of course Planet Nine (if it exists) isn't the only mass in the outer Solar System.
PM 2Ring
CC BY-SA 4.0
2020-07-12T19:22:10.120Z
@PM2Ring, wouldn't the smaller masses like comets be scattered at random, and so cancel each other out, on average? That would make easier to filter the signal of large mass objects.
ksousa
CC BY-SA 4.0
2020-07-12T20:09:01.227Z
@ksousa Yes, the stuff that's randomly distributed *should* mostly cancel out, but a few non-random lumps of mass could counter-balance Planet 9. OTOH, it's unlikely that there's a lot of mass out there, regardless of how it's distributed. Hopefully, someone who knows more about this stuff will write an actual answer. ;)
PM 2Ring
CC BY-SA 4.0
2020-07-13T21:03:41.767Z
@PM2Ring the estimated total mass of the Kuiper belt per this Wikipedia source arXiv:1811.05191. doi:10.1134/S1063773718090050. is only about 2% of the mass of the Earth. This was based on gravitational influence on the planets and spacecraft. That being said, I would think that if Planet 9 exists, it's effects would have been seen in this study.
Jack R. Woods
<planet><solar-system><orbital-mechanics><9th-planet><space-probe>
Can the paper narrowing Solar System's barycentre to within 100m help find Planet Nine?
ksousa
# Can the paper narrowing Solar System's barycentre to within 100m help find Planet Nine?

(Asked by: username_0 on 2020-07-12)

I've just read this recent news article, [Astronomers Have Located The Centre of The Solar System to Within 100 Metres](https://www.sciencealert.com/astronomers-have-located-the-centre-of-the-solar-system-to-within-100-metres), reporting on a paper[1] that claims narrowing the Solar System barycentre to within 100 meters, just outside the sun, and I wonder, could this result help prove or disprove the hypothesis of a massive [Planet Nine](https://en.wikipedia.org/wiki/Planet_Nine) as proposed by Batygin and Brown[2] in 2016?

I imagine that any unknown mass concentration would change the barycentre and generate a error, and if we could measure this error, it should be possible to pinpoint the position of said mass concentration. Can we just use probe tracking data from the [Voyagers](https://en.wikipedia.org/wiki/Voyager_program) or [Pioneers](https://en.wikipedia.org/wiki/Pioneer_10) to see if their hyperbolic trajectories match the barycentre reported in the paper, for example?

[1]. [Vallisneri, M., et al. “Modeling the Uncertainties of Solar System Ephemerides for Robust Gravitational-wave Searches with Pulsar-timing Arrays”. The Astrophysical Journal, vol. 893, no 2, abril de 2020, p. 112. DOI.org (Crossref), doi:10.3847/1538-4357/ab7b67.](https://arxiv.org/pdf/2001.00595.pdf)

[2]. [Batygin, Konstantin, e Michael E. Brown. “Evidence for a Distant Giant Planet in the Solar System”. The Astronomical Journal, vol. 151, no 2, janeiro de 2016, p. 22. arXiv.org, doi:10.3847/0004-6256/151/2/22.](https://arxiv.org/abs/1601.05438)

**Tags:** planet, solar-system, orbital-mechanics, 9th-planet, space-probe

> **Comments:**
> **username_1**: Here's a recent press release about an alternative explanation for the oddities of the outer Solar System that doesn't require Planet Nine: https://www.colorado.edu/today/2020/07/07/outer-solar-system  
> **username_1**: FWIW, a rough calculation of $rm/M_\odot$ for Planet 9, using $r=400$ AU & $m=5M_\oplus=15×10^{-6}M_\odot$, gives around 900,000 km. That's a pretty large contribution to the barycentre. But of course Planet Nine (if it exists) isn't the only mass in the outer Solar System.  
> **username_0**: @PM2Ring, wouldn't the smaller masses like comets be scattered at random, and so cancel each other out, on average? That would make easier to filter the signal of large mass objects.  
> **username_1**: @username_0 Yes, the stuff that's randomly distributed *should* mostly cancel out, but a few non-random lumps of mass could counter-balance Planet 9. OTOH, it's unlikely that there's a lot of mass out there, regardless of how it's distributed. Hopefully, someone who knows more about this stuff will write an actual answer. ;)  
> **username_2**: @PM2Ring the estimated total mass of the Kuiper belt per this Wikipedia source arXiv:1811.05191. doi:10.1134/S1063773718090050. is only about 2% of the mass of the Earth. This was based on gravitational influence on the planets and spacecraft. That being said, I would think that if Planet 9 exists, it's effects would have been seen in this study.  

---

## 1 answers

### Answer by username_3 on 2020-07-14. Score: 20 | Accepted answer

# Short answer

Short answer is no, it can't prove or disprove the existence of Planet 9. Even if there is a significant difference between the barycenter of the solar system with and without Planet 9 we wouldn't be able to tell without hundreds if not thousands of years of precise data. If we don't have Planet 9 as a reference and are limited to the Sun and the eight known planets for reference then without another reference in the same frame we have no way to observe the motion of our reference.

It is perhaps useful to consider the solar system (i.e. the Sun plus the eight known planets) as orbiting *together* around the barycenter of the solar system plus Planet 9 system. Without seeing Planet 9 as a reference how would we have any notion that this other barycenter even exists?

Of course Planet 9 would still be exerting gravitational influence and over time this would cause a deviation between where the barycenter *should* be and where it actually is -- that is the only way we could use the motion of the barycenter to infer anything about individual bodies. Essentially, the barycenter is a convolution of the gravitational influence of all the mass,5 its deconvolution requires analysis of temporal variation.

---

# Unnecessarily long answer

**Note**   -   When I discuss variations in the barycenter I will be referring to the barycenter of the solar system in heliocentric, not barycentric, coordinates (in barycentric coordinates of course there is no variation in the barycenter, by definition). This simplifies some of the qualitative analysis.

---

Rigorously answering this requires answering several separate questions (some which aren't necessarily relevant in answering the broader question but do provide useful background):

1. Does the location of the barycenter over time contain enough information to recover information about the contributing masses?
2. How precisely can we predict the location of the barycenter?
3. How much deviation would a $5$-$10~M$⊕ planet at $400$-$800$ AU cause to the center of mass of the solar system?
4. If no such deviation exists, does that rule out the existence of Planet 9? Conversely, if such a deviation does exist, does it rule **in** the existence of Planet 9?

---

### 1 - Location of the barycenter over time

The short answer to question **1** is yes (sort of). There is information on the orbits of the contributing masses which can be backed out from simple analysis of the location of the barycenter **over time**.

The emphasis here is because at any given point in time the location of the barycenter contains no information about the location of individual bodies, only the entire mass distribution. However, since most of the mass in the solar system is discretized into large bodies (as opposed to a gas or dust cloud, protoplanetary disk, etc), the influence of each massive body on the location of the barycenter will vary with a frequency proportional to its orbital period -- and, critically, this is fairly measurable (but only if we can also know the relative positions of the bodies). Consequently, the influence of the individual masses can be retrieved (this is not necessarily the case for gas or dust clouds and protoplanetary disks).

In a simple system, retrieving information about the orbits is fairly straightforward (although I will neglect the rigorous demonstration for brevity). The issues arise when the system becomes more complicated.

One source of confound is resonant orbits. Certain resonances would defy deconvolution because the frequencies of their respective influences on the barycenter are synchronous. Furthermore, multiple resonances can't necessarily be disambiguated: the frequency pattern in the motion of the barycenter caused by 2 bodies in a resonance could be replicated by 3 bodies in a resonance. Pluto and Neptune are in a [2:3 resonance](https://en.wikipedia.org/wiki/Orbital_resonance#Mean-motion_resonances_in_the_Solar_System), so I don't think we can immediately dismiss the possibility of a theoretical Planet 9 being in a resonance of some kind and that this could hinder deconvolving its effect on the barycenter.

Another significant challenge is the qualification that the barycenter must be observed **over time**, and more specifically how much time is needed to deconvolve the various influences. As I mentioned, each body's influence on the barycenter is cyclic based on their orbital period and eccentricity. Using the Bulirsch-Stoer integrator from [Vulcan](https://github.com/WMiller256/Vulcan) and tracking an energy analog of the barycenter's motion1 as a function of time we can analyze the effect of the different orbital frequencies.

Let's consider a three body system of the Sun, Saturn, and Jupiter. In this system the influence of Saturn and Jupiter should clearly be visible in the resulting sinusoid. After integrating this system for $2$ Saturn years ($60$ Earth years) the energy-analogy over time looks like this:

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

It should be pretty clear that the above behavior can be created by two sinusoids with different frequencies, to any who find this unclear I recommend having a look at [this tool](https://demonstrations.wolfram.com/CombiningSineWaves/). If wanted to we could use Fourier analysis to deconvolve those two frequencies and recover the orbital periods of Jupiter and Saturn (I will also neglect this for brevity).

Let's consider a slightly worse case: a three body system comprised of the Sun, Jupiter, and a hypothetical [hot jupiter](https://en.wikipedia.org/wiki/Hot_Jupiter#:%7E:text=Hot%20Jupiters%20are%20a%20class,the%20moniker%20%22hot%20Jupiters%22.) with a semi-major axis of $1.55^{}$x$10^{10}$ m on a relatively circular orbit. Our energy-analog after $3$ hot-jupiter years ($12$ Earth days) looks like this:

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

**Where did the other sinusoid go?**

Its effect is still there, but we haven't sampled sufficient time to see it. We would need to look at $1,095$ hot-jupiter years ($12$ Earth years) in order to see one orbit of our toy Jupiter2 Fortunately we don't have to integrate quite that long to see the effect of the toy Jupiter on the energy-analog. After $320$ hot-jupiter years ($3.5$ Earth years):

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

This becomes even more complicated with non-zero eccentricities and unstable orbits. For a similar three body, hot-jupiter system but with the toy Jupiter having an eccentricity of $0.65$ over a period of $7$ Earth years the energy-analog is:3

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

There are other ways to draw inferences about the orbital parameters from the relative motion of the barycenter (e.g. the shape of individual cycles), but that is for a different question.

**In sum, the location of the barycenter over time could certainly be used to find a mass imbalance, but but not necessarily argue for the existence of Planet 9 -- that would depend on the specifics of the observations.**

---

## 2 - Precision of barycenter predictions

There are number of significant barriers to accurately calculating the barycenter of the solar system but the most troublesome are the uncertainties regarding the interiors of Saturn and Jupiter. Specifically, the behavior of liquid metallic hydrogen at such massive pressures and (more so for Saturn) understanding their gravitational moments (Fortney [2004](https://zenodo.org/record/1230846#.Xw02XihKiiN)).

The issue this raises is that without knowing the centers of mass for Jupiter and Saturn (which contain 92% of planetary mass in the solar system) with sufficient precision we don't know well enough what the solar system barycenter *should* be to determine whether or not the true barycenter differs enough to indicate the existence of Planet 9.

With information from the Cassini spacecraft combined with data from the [VLBA radio telescope](https://public.nrao.edu/telescopes/vlba/) the Saturn ephemerides was constrained to $4$ km in [2015](http://nasa.gov/jpl/cassini/scientists-pinpoint-saturn-with-exquisite-accuracy). This represented an improvement of about an order of magnitude. The Juno spacecraft provided a similar improvement to the Jupiter ephemerides, constraining it to within $10$ km in [2019](https://ui.adsabs.harvard.edu/abs/2019AAS...23330201J/abstract).

What is less well constrained is how constant these parameters are, so it bears noting that an inaccuracy of $350$ km in the location of the center of mass of Saturn's core within the planet corresponds to an uncertainty of $100$ m in the position of the solar system barycenter. Jupiter lacks a dense core like Saturn's, so it is more difficult to quantify a direct relationship between inaccuracies, but suffice it to say that an inaccuracy of $100$ km in Jupiter's center of mass corresponds to an uncertainty of $100$ m in the barycenter.4

So assuming Saturn's center of mass has not drifted much since September of 2018, **I would estimate the precision of the calculated solar system barycenter to be roughly $\pm15$ m.**

---

## 3 - Deviation from predicted barycenter

Calculating the effect Planet 9 would have on the barycenter is actually quite trivial thanks to superposition (see [this answer](https://stackoverflow.com/questions/59768976/how-can-i-work-out-the-gravitational-force-of-any-2d-polygon/59797041#59797041) for a more rigorous treatment of the limitations of this concept when it comes to gravitational interactions). The equation for calculating the distance from the primary body to the barycenter in a two-body system,

$$
r\_1 = \frac{a}{1 + \frac{m\_1}{m\_2}}
$$

can also be applied to calculating the barycenter of two barycenters where $a$ is the distance between the two barycenters, $m\_1$ is the mass which contributes to the primary barycenter, $m\_2$ is the mass which contributes to the secondary barycenter, and $r\_1$ is the distance from the primary barycenter to the shared barycenter. And since $99.86$% of the mass of the solar system is contained in the sun it is a rigorous approximation to use

$$
r\_1 = \frac{200 \cdot 1.496 \text{x} 10^{11}}
{1 + \frac{1.989e30}{5 \cdot 5.972 \text{x} 10^{24}}}
= 4 \text{x} 10^8 \text{ m}
$$

as a lower bound, and

$$
r\_1 = \frac{1200 \cdot 1.496 \text{x} 10^{11}}
{1 + \frac{1.989e30}{10 \cdot 5.972 \text{x} 10^{24}}}
= 5 \text{x} 10^9 \text{ m}
$$

as an upper bound for the distance from the barycenter of the solar system without Planet 9 to the barycenter with it.

---

## 4 - Interpretation

Even the lower bound of $4\text{x}10^{8}$ m is a huge deviation, so did we just demonstrate that Planet 9 can't exist? Unfortunately no. This goes back to number **1** - the problem is time. There are a number of ways to think about this, but the salient point is that the influence on the barycenter of the solar system would have been persistent since well before the earliest observational records. Planet 9 would have to have an orbital period on the order of $10,000$ years. Over periods of thousands of years it would have obvious and readily measurable effects on the motion of the other planets (and that motion would be reflected in the motion of the heliocentric barycenter), but otherwise the effects are simply too small to disambiguate from noise.

---

1 I create this analog by calculating a hypothetical energy of the system if the center of the sun was the barycenter of the system and comparing it to the initial value. It is not tremendously physically meaningful but it does reduce the parameter-space of the temporal position of the barycenter to two dimensions, thereby making it easier to show the effect I discuss graphically.

2 This is one of the reasons simulating hot jupiters is so difficult (and one of the things the Vulcan project is intended to remedy): the limiting factor for orbital integration is always the fastest-orbiting body (although multi-step integrators mitigate this somewhat) so systems with large differences between orbital periods take proportionally longer to simulate. Using the `Mercury` code, a gigayear solar system simulation with a hot-jupiter on an 8 day orbit would require roughly 10 months of simulation time on a high-end desktop computer.

3 Note that the energy analog is degenerate here because the hot-jupiter will be ejected after approximately 10 Earth years.

4 These uncertainty relationships are from my own calculations and may be flawed as they have not been rigorously tested or reviewed.

5 And I do literally mean **all** the mass. As in all the mass in the universe, gravity has an infinite range after all....

---

> **Comments:**
> **username_3**: I realize the plots are lacking axis titles, I intend to remedy that but have to modify the Vulcan source to do so.  
> **username_4**: I've taken the liberty of backing up the current state of Vulcan to archive.org to have some host redundancy to prevent link rot (I credited it to you): https://archive.org/details/vulcan-master  

---

