Issue 
A&A
Volume 646, February 2021



Article Number  L14  
Number of page(s)  9  
Section  Letters to the Editor  
DOI  https://doi.org/10.1051/00046361/202140311  
Published online  18 February 2021 
Letter to the Editor
Memories of past close encounters in extreme transNeptunian space: Finding unseen planets using pure random searches
^{1}
Universidad Complutense de Madrid, Ciudad Universitaria, 28040 Madrid, Spain
email: nbplanet@ucm.es
^{2}
AEGORA Research Group, Facultad de Ciencias Matemáticas, Universidad Complutense de Madrid, Ciudad Universitaria, 28040 Madrid, Spain
Received:
11
January
2021
Accepted:
3
February
2021
Context. The paths followed by the known extreme transNeptunian objects (ETNOs) effectively avoid direct gravitational perturbations from the four giant planets, yet their orbital eccentricities are in the range between 0.69−0.97. Solar system dynamics studies show that such high values of the eccentricity can be produced via close encounters or secular perturbations. In both cases, the presence of yettobediscovered transPlutonian planets is required. Recent observational evidence cannot exclude the existence, at 600 AU from the Sun, of a planet of five Earth masses.
Aims. If the high eccentricities of the known ETNOs are the result of relatively recent close encounters with putative planets, the mutual nodal distances of sizeable groups of ETNOs with their assumed perturber may still be small enough to be identifiable geometrically. In order to confirm or reject this possibility, we used Monte Carlo random search techniques.
Methods. Two arbitrary orbits may lead to close encounters when their mutual nodal distance is sufficiently small. We generated billions of random planetary orbits with parameters within the relevant ranges and computed the mutual nodal distances with a set of randomly generated orbits with parameters consistent with those of the known ETNOs and their uncertainties. We monitored which planetary orbits had the maximum number of potential close encounters with synthetic ETNOs and we studied the resulting distributions.
Results. We provide narrow ranges for the orbital parameters of putative planets that may have experienced orbitchanging encounters with known ETNOs. Some sections of the available orbital parameter space are strongly disfavored by our analysis.
Conclusions. Our calculations suggest that more than one perturber is required if scattering is the main source of orbital modification for the known ETNOs. Perturbers might not be located farther than 600 AU and they have to follow moderately eccentric and inclined orbits to be capable of experiencing close encounters with multiple known ETNOs.
Key words: methods: data analysis / methods: numerical / celestial mechanics / planets and satellites: detection / minor planets / asteroids: general / Kuiper belt: general
© ESO 2021
1. Introduction
Extreme transNeptunian objects (ETNOs) serve as unique probes into the gravity perturbations shaping the outer solar system beyond the classical transNeptunian or Kuiper belt (see e.g., Kaib et al. 2019). The trajectories followed by the known ETNOs effectively avoid direct gravitational perturbations from the four giant planets, yet their orbital eccentricities are in the range 0.69−0.97. Efficient drivers for the eccentricity excitation of small bodies include close encounters with planets (see e.g., Carusi et al. 1990) and the von ZeipelLidovKozai mechanism (von Zeipel 1910; Kozai 1962; Lidov 1962; Ito & Ohtsuka 2019). In both cases above and in the case of ETNOs, the presence of yettobediscovered transPlutonian planets is required.
Fienga et al. (2020) used the INPOP19a planetary ephemerides that include Jupiterupdated positions by the Juno mission and a reanalysis of Cassini observations to show that there is no clear evidence for the existence of the socalled Planet 9 predicted by Batygin & Brown (2016) as an explanation for the orbital architecture of the known ETNOs. However, Fienga et al. (2020) concluded that if Planet 9 exists, it cannot be closer than 500 AU, if it has a mass of 5 M_{⊕}, and no closer than 650 AU, if it has a mass of 10 M_{⊕}. The latest version of the Planet 9 hypothesis (Batygin et al. 2019) predicts the existence of a planet with a mass in the range 5−10 M_{⊕}, following an orbit with a value of the semimajor axis in the range of 400−800 AU, eccentricity in the range of 0.2−0.5, and inclination in the interval between (15°, 25°). A number of exoplanets have already been observed orbiting at hundreds of AUs from their host stars (see e.g., Bailey et al. 2014; Naud et al. 2014; Nguyen et al. 2021) and theoretical calculations have confirmed plausible pathways for their formation (see e.g., Kenyon & Bromley 2015, 2016).
If the high eccentricities of known ETNOs are the result of relatively recent close encounters with putative planets, the mutual nodal distances of sizeable groups of ETNOs with their assumed perturber may still be small enough to be identifiable geometrically. Here, we use Monte Carlo random search techniques to identify orbits that may lead to the maximum number of potential close encounters with synthetic ETNOs whose orbital parameters are consistent with those of the real ETNOs and their uncertainties. This Letter is organized as follows. In Sect. 2, we review our methodology and present the data used in our analyses. In Sect. 3, we apply our methodology and discuss its results. Our conclusions are summarized in Sect. 4.
2. Methods and data description
This work explores a “what if” scenario in which the starting hypothesis states that a sizeable number of known ETNOs have experienced relatively recent close encounters with putative planets; the timescale comes constrained by their orbital periods that are in the range 1867−50116 yr, so encounters may have taken place during the last 10^{3} − 10^{5} yr. In this work, therefore, we are testing the hypothesis statistically. If this hypothesis is plausible, a statistically significant number of compatible planetary perturber orbits should emerge from the analysis of a very large sample of orbits. If the distributions of one or more of the orbital parameters of the candidate are flat, the starting hypothesis must be rejected as this would show that there is no favored orbital solution for the perturber; conversely, if all the distributions produce consistent intervals that are statistically significant, the plausibility of the starting hypothesis can be considered as confirmed. Plausibility concerns the likelihood of acceptance, not the likelihood of being true or better than competing scenarios. The problem under investigation here is equivalent to a nondifferenciable optimization that is wellsuited for a uniform Monte Carlo random search (Metropolis & Ulam 1949). Our methodology brings together geometry and statistics in our attempt to find the confocal ellipse that passes the closest to the maximum number of known orbits of a certain dynamical class; this approach is fundamentally different from those involving Nbody calculations and statistics (see e.g., de la Fuente Marcos et al. 2016, 2017; de León et al. 2017). Kalinicheva & Chernetenko (2020) recently applied a geometric approach within the context of the Planet 9 hypothesis.
2.1. Methodology
Two arbitrary orbits may experience close encounters when their mutual nodal distance is sufficiently small. Recurrent (or even single) encounters within 1 Hill radius (Hamilton & Burns 1992) of a massive body may change the orbit of a small body significantly. The mutual nodal distance between the orbits of a small body (an ETNO in our case) and an arbitrary planet can be computed as described in Appendix A. Orbits are defined by the values of semimajor axis, a (that controls size), eccentricity, e (that controls shape), and those of the angular elements – inclination, i, longitude of the ascending node, Ω, and argument of perihelion, ω – that control the orientation in space of the orbit; the perihelion distance or pericenter, q, is given by the expression q = a (1 − e). The actual position of an object in its orbit is controled by a fourth angle, the mean anomaly, M. Our geometric approach leaves this angle out of the analysis and, therefore, it is not capable of predicting the current location of the perturber, if it is, in fact, real.
We generated 2 × 10^{10} random planetary orbits with uniformly distributed relevant parameters: Ω_{p} and ω_{p} ∈(0°, 360° ), i_{p} ∈ (0°, 80° ), e_{p} ∈ (0, 0.9), and q_{p} ∈ (x, 1000 − x) AU, with x = 300, 400, 500, and 600 AU so a_{p} = q_{p}/(1 − e_{p}). For each random planetary orbit, we computed the mutual nodal distances between a set of synthetic ETNOs and the planet. Each set of synthetic ETNOs was generated using the mean values and standard deviations of the orbit determinations of the known ETNOs as pointed out in Appendix B.
For each combination of random planetary orbit and synthetic (but compatible with the observations) set of ETNOs, we have counted how many synthetic ETNOs had at least one mutual nodal distance with the planet under 5 AU (for x = 300 AU), 7.5 AU (for x = 400 AU), and 10 AU (for x = 500 and 600 AU). We then proceeded to record the random planetary orbit if the count was ≥5, in order to maximize the number of potential close approaches between planet and set of ETNOs. We then studied the resulting distributions. In order to analyze the results, we produced histograms using the Matplotlib library (Hunter 2007) with sets of bins computed using NumPy (van der Walt et al. 2011) by applying the Freedman and Diaconis rule (Freedman & Diaconis 1981). Instead of using frequencybased histograms, we considered counts to form a probability density so the area under the histogram will sum to one.
The nodal distance separation criteria for selective counting are not arbitrary but motivated by the results in Fienga et al. (2020), a 2 M_{⊕} has a Hill radius of 4.5 AU (if a_{p} = 400 AU and e_{p} = 0.1), a 5 M_{⊕} has a Hill radius of 9.2 AU (if a_{p} = 542 AU and e_{p} = 0.01). Therefore, we implicitly assume that the farther away the pericenter of the planet is, the more massive it may be.
2.2. Data sources
Here, we work with publicly available data from Jet Propulsion Laboratory’s (JPL) SmallBody Database (SBDB)^{1} and HORIZONS online solar system data and ephemeris computation service^{2}, both provided by the Solar System Dynamics Group (Giorgini 2015). Assuming the definition in Trujillo & Sheppard (2014), that the ETNOs have q > 30 AU and a > 150 AU, the sample of known ETNOs now includes 39 objects with reliable orbits (see Appendix B) whose data have been retrieved from JPL’s SBDB and HORIZONS using tools provided by the Python package Astroquery (Ginsburg et al. 2019). In the following section, we use barycentric elements because within the context of the ETNOs, barycentric orbit determinations better account for their changing nature as Jupiter follows its 12 yr orbit around the Sun; Appendix C shows results based on heliocentric orbits that are consistent with those of the barycentric ones.
3. Results and discussion
As pointed out above, the motivation behind this study is the belief that a fossil record of planetary encounters might be preserved in the distribution of the orbital elements of ETNOs. With this hypothesis in mind, we monitored which planetary orbits had the maximum number of potential close encounters with synthetic ETNOs and analyzed the resulting distributions that are shown in Fig. 1. The first (leftmost) column of panels shows the results based on the assumption that q_{p} > 300 AU, including 20 304 orbits with a number of potential close approaches in the range between 5−7; the next one, shows results for q_{p} > 400 AU also with a range of potential close approaches of 5−7 for 5671 orbits; the following column of panels displays results for q_{p} > 500 AU with a range of 5−6 for 2635 orbits; the right column shows panels with results for q_{p} > 600 AU and the number of potential close approaches between planet and set of ETNOs is just 5 for 56 orbits. From these values and the plots, it is increasingly difficult to find consistent perturbers as we move farther away from the Sun.
Fig. 1.
Barycentric orbital elements of putative perturber planets. Distributions of barycentric orbital elements of planetary orbits that may result in close encounters (under 5 AU for q_{p} > 300 AU, under 7.5 AU for q_{p} > 400 AU, and under 10 AU for q_{p} > 500 AU and q_{p} > 600 AU) with five or more presentday extreme transNeptunian objects (ETNOs). Each column of panels shows the cumulative results of 2 × 10^{10} experiments. From left to right: results of imposing q_{p} > 300 AU (left column, 20 304 orbits), q_{p} > 400 AU (5671 orbits), q_{p} > 500 AU (2635 orbits), and q_{p} > 600 AU (right column, 56 orbits). Median values are shown as vertical blue lines. Based on the available data on known ETNOs, the presence of massive perturbers well beyond 600 AU is strongly excluded within the context of the hypotheses centered in this study. 
Table 1 shows a summary of the central values and dispersions of the orbital parameters of the sample of orbits that may experience close encounters with multiple known ETNOs. We consider that results for q_{p} > 300 AU and q_{p} > 400 AU are statistically consistent and are compatible with a_{p} ∈ (331, 489) AU, e_{p} < 0.1, i_{p} ∈ (10° , 20° ), Ω_{p} ∼ 180°, and ω_{p} ∈ ( − 60° , 60° ). The value of a_{p} is probably ∼400 AU but the uncertainty is significant, the value of the eccentricity is well constrained and it has to be low, the inclination is perhaps the best constrained value and it is ∼14°, the value of Ω_{p} is very likely ∼180°, but the value of ω_{p} is poorly constrained, perhaps ∼50°. If scattering is the main source of orbital modification for the group of ETNOs that mainly move at (300, 400) AU from the Sun, the orbit of the putative perturber is relatively wellconstrained and according to Fienga et al. (2020), it must have a mass < 5 M_{⊕}.
Summary of central values and dispersions of optimal barycentric orbits.
Table 1 shows that our approach is far less conclusive in the cases of q_{p} > 500 AU and q_{p} > 600 AU as these values produce distributions of the orbital parameters that may not be statistically compatible. In any case, we must emphasize that > 50% of the known ETNOs cannot interact directly with a perturber with q_{p} > 500 AU because their aphelion distances, Q = a (1 + e), are below 500 AU: 9 ETNOs have Q < 300 AU, 18 have Q < 400 AU, and 23 have Q < 500 AU. Therefore, the distribution in Q and our calculations suggest that more than one perturber is required if scattering is the main source of orbital modification for the known ETNOs. Perturbers might not be located farther than 600 AU and they have to follow moderately eccentric and inclined orbits to be capable of experiencing presentday close encounters with multiple known ETNOs.
At this point, it can be argued that the results in Fig. 1 may be affected by a statistical artifact. In order to test the statistical significance of our results, we repeated the Monte Carlo random search on an input sample of scrambled data as explained in Appendix D. By randomly reassigning the values of the orbital elements of the ETNOs, we preserve the original distributions of the parameters, but destroy any possible correlations that may have been induced by close encounters with massive perturbers (or the effects of hypothetical meanmotion or secular resonances). The results of this significance test are shown in Fig. D.1: the distribution of i_{p} becomes flat and almost the same happens to the distribution of Ω_{p}. In other words, for the scrambled data, it is not possible to find a statistically significant orbital solution that is compatible with the starting hypothesis and we conclude that the results in Fig. 1 are unlikely to be due to statistical artifacts.
4. Summary and conclusions
When considering the wellstudied case of Neptune and the regular transNeptunian objects, we observe that these objects are not part of a single population, but they are organized into several dynamical classes. Some objects never experience close encounters with Neptune due to the existence of protection mechanisms such as meanmotion or secular resonances, as in the case of Neptune’s Trojans or Pluto (see e.g., Milani et al. 1989; Wan et al. 2001), others undergo close encounters that may send them towards the inner solar system (centaurs) or outwards to become scattered objects (see the recent review by Saillenfest 2020). If transPlutonian planets exist, their perturbations may shape the extreme transNeptunian space in a similar fashion and the current sample of ETNOs may include the signatures of their presence (see e.g., Saillenfest et al. 2017).
Our results show that it is possible to find suitable orbits for which the mutual nodal distances of sizeable (≥5, comprising at least 13% of the sample) groups of ETNOs with their assumed perturber could be small enough for a close encounter to occur, at least in theory (assuming that no protection mechanisms, such as meanmotion or secular resonances, are in place to avoid the flyby). This was our original aim. In addition, the results presented in Sect. 3 clearly indicate that the most probable planetary orbit consistent with the starting hypothesis is the one obtained for the experiment with q_{p} > 300 AU. The number of consistent orbits in this case is 3.6 times higher than the one found for q_{p} > 400 AU. Our results seem to be incompatible with those attributed to a statistical artifact (see Appendix D).
Prior to the announcement of the Planet 9 hypothesis (Batygin & Brown 2016), Trujillo & Sheppard (2014) had already argued for the existence of a planet at 250 AU within the context of the ETNOs – driving von ZeipelLidovKozai librations on 2012 VP_{113} – and de la Fuente Marcos & de la Fuente Marcos (2014) suggested that the limited data available at the time were more compatible with the presence of two massive perturbers, one of them close to 300 AU. These massive perturbers were initially proposed based on data corresponding to a sample of 13 objects, whereas the current sample has tripled this number. If we repeat the experiment discussed in Sect. 3 for q_{p} > 200 AU (see Appendix E), we find that the number of consistent orbits, although larger than the one generated in the experiment for q_{p} > 400 AU, is still lower than that of the most statistically significant case, namely, 8234 versus 20 304 for q_{p} > 300 AU. The most likely orbit is still similar, in terms of shape and orientation in space, to the most probable one in the q_{p} > 300 AU case.
Our results are consistent with the conclusions of the study presented in de la Fuente Marcos & de la Fuente Marcos (2017). Although our approach has not been able to single out a statistically significant, unique planetary orbit that may be responsible for the orbital architecture observed in extreme transNeptunian space, we provide narrow ranges for the orbital parameters of putative planets that may have experienced orbitchanging encounters with known ETNOs. Some sections of the available orbital parameter space are strongly excluded by the findings of our analysis.
Acknowledgments
We thank the anonymous referee for a constructive and timely report, S. J. Aarseth, J. de León, J. Licandro, A. CabreraLavers, J.M. Petit, M. T. Bannister, D. P. Whitmire, G. Carraro, E. Costa, D. Fabrycky, A. V. Tutukov, S. Mashchenko, S. Deen and J. Higley for comments on ETNOs and A. I. Gómez de Castro for providing access to computing facilities. This work was partially supported by the Spanish ‘Ministerio de Economía y Competitividad’ (MINECO) under grant ESP201787813R. In preparation of this Letter, we made use of the NASA Astrophysics Data System and the MPC data server.
References
 Bailey, V., Meshkat, T., Reiter, M., et al. 2014, ApJ, 780, L4 [Google Scholar]
 Batygin, K., & Brown, M. E. 2016, AJ, 151, 22 [Google Scholar]
 Batygin, K., Adams, F. C., Brown, M. E., et al. 2019, Phys. Rep., 805, 1 [Google Scholar]
 Carusi, A., Valsecchi, G. B., & Greenberg, R. 1990, Celest. Mech. Dyn. Astron., 49, 111 [Google Scholar]
 de la Fuente Marcos, C., & de la Fuente Marcos, R. 2014, MNRAS, 443, L59 [Google Scholar]
 de la Fuente Marcos, C., & de la Fuente Marcos, R. 2017, MNRAS, 471, L61 [Google Scholar]
 de la Fuente Marcos, C., de la Fuente Marcos, R., & Aarseth, S. J. 2016, MNRAS, 460, L123 [Google Scholar]
 de la Fuente Marcos, C., de la Fuente Marcos, R., & Aarseth, S. J. 2017, Ap&SS, 362, 198 [Google Scholar]
 de León, J., de la Fuente Marcos, C., & de la Fuente Marcos, R. 2017, MNRAS, 467, L66 [Google Scholar]
 Fienga, A., Di Ruscio, A., Bernus, L., et al. 2020, A&A, 640, A6 [CrossRef] [EDP Sciences] [Google Scholar]
 Freedman, D., & Diaconis, P. 1981, Prob. Theor. Relat. Fields, 57, 453 [Google Scholar]
 Ginsburg, A., Sipőcz, B. M., Brasseur, C. E., et al. 2019, AJ, 157, 98 [Google Scholar]
 Giorgini, J. D. 2015, IAU Gen. Assem., 22, 2256293 [NASA ADS] [Google Scholar]
 Hamilton, D. P., & Burns, J. A. 1992, Icarus, 96, 43 [Google Scholar]
 Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
 Ito, T., & Ohtsuka, K. 2019, Monogr. Environ. Earth Planets, 7, 1 [Google Scholar]
 Kaib, N. A., Pike, R., Lawler, S., et al. 2019, AJ, 158, 43 [Google Scholar]
 Kalinicheva, O. V., & Chernetenko, Y. A. 2020, Astrophys. Bull., 75, 459 [Google Scholar]
 Kenyon, S. J., & Bromley, B. C. 2015, ApJ, 806, 42 [Google Scholar]
 Kenyon, S. J., & Bromley, B. C. 2016, ApJ, 825, 33 [Google Scholar]
 Kozai, Y. 1962, AJ, 67, 591 [Google Scholar]
 Lidov, M. L. 1962, Planet. Space Sci., 9, 719 [Google Scholar]
 Metropolis, N., & Ulam, S. 1949, J. Am. Stat. Assoc., 44, 335 [Google Scholar]
 Milani, A., Nobili, A. M., & Carpino, M. 1989, Icarus, 82, 200 [Google Scholar]
 Naud, M.E., Artigau, É., Malo, L., et al. 2014, ApJ, 787, 5 [Google Scholar]
 Nguyen, M. M., De Rosa, R. J., & Kalas, P. 2021, AJ, 161, 22 [Google Scholar]
 Saillenfest, M. 2020, Celest. Mech. Dyn. Astron., 132, 12 [Google Scholar]
 Saillenfest, M., Fouchard, M., Tommei, G., et al. 2017, Celest. Mech. Dyn. Astron., 129, 329 [Google Scholar]
 Trujillo, C. A., & Sheppard, S. S. 2014, Nature, 507, 471 [Google Scholar]
 van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22 [Google Scholar]
 von Zeipel, H. 1910, Astron. Nachr., 183, 345 [Google Scholar]
 Wan, X.S., Huang, T.Y., & Innanen, K. A. 2001, AJ, 121, 1155 [Google Scholar]
Appendix A: Mutual nodal distance: formulae
Our statistical analyses are based on studying the distribution of nodal distances between two Keplerian trajectories (one ETNO and one hypothetical planet) with a common focus; therefore, the core of our approach is purely geometrical and gravitational interactions are not directly taken into account. The mutual nodal distance between the orbits of a small body (an ETNO in our case) and an arbitrary planet can be written as (see Eqs. (16) and (17) in Saillenfest et al. 2017):
where
ΔΩ_{p} = Ω − Ω_{p}, and a, e, i, Ω and ω are the orbital elements of the small body, and a_{p}, e_{p}, i_{p}, Ω_{p}, and ω_{p} those of the arbitrary planet. For each random search experiment (each analysis consists of 2 × 10^{10}), we compute the orbital elements of the putative perturber using the expressions:
where x = 300, 400, 500, and 600 AU (see Sect. 2, or 200 AU for Appendix E), and r_{j} with j = 1, 5, are random numbers in the interval (0, 1) with a uniform distribution.
Appendix B: Extreme transNeptunian objects: Data
The barycentric orbit determinations used as input data in the uniform Monte Carlo random search discussed in Sect. 3 are shown in Table B.1. The data are referred to epoch 2459000.5 Barycentric Dynamical Time (TDB) and they have been retrieved (as of 11January2021) from JPL’s SBDB and HORIZONS using tools provided by the Python package Astroquery (Ginsburg et al. 2019).
Barycentric orbital elements and 1σ uncertainties of known ETNOs with robust orbit determinations.
The orbits of the synthetic ETNOs are obtained using the expressions:
where R_{j} with j = 1, 5, are univariate Gaussian random numbers.
Appendix C: Heliocentric orbit determinations: results
If we repeat the calculations discussed in Sect. 3 using heliocentric orbit determinations instead of barycentric ones as input data, we obtain the distributions in Fig. C.1 with central values and dispersions summarized in Table C.1.
Fig. C.1.
Heliocentric orbital elements of putative perturber planets. Distributions of heliocentric orbital elements of planetary orbits that may result in close encounters (under 5 AU for q_{p} > 300 AU, under 7.5 AU for q_{p} > 400 AU, and under 10 AU for q_{p} > 500 AU and q_{p} > 600 AU) with five or more presentday extreme transNeptunian objects (ETNOs). Each column of panels shows the cumulative results of 2 × 10^{10} experiments. From left to right: results of imposing q_{p} > 300 AU (left column, 18 369 orbits), q_{p} > 400 AU (6594 orbits), q_{p} > 500 AU (3567 orbits), and q_{p} > 600 AU (right column, 105 orbits). Median values are shown as vertical blue lines. Based on the available data on known ETNOs, the presence of massive perturbers well beyond 600 AU is strongly excluded within the context of the hypotheses centered in this study. 
Summary of central values and dispersions of optimal heliocentric orbits.
Appendix D: Statistical significance
In order to verify that our results do not come from statistical artifacts, we randomly scramble the orbit parameters data used as input and repeat the uniform Monte Carlo random searches discussed in Sect. 3 and Appendix C. In these experiments, the set of synthetic ETNOs is such that, for example, the first object may have the value of a of object #5, e of #27, i of #7, Ω of #37, and ω of #13. By randomly rearranging the values of the orbital elements of the ETNOs, we retain the original distributions of the parameters, but destroy any possible correlations existing among them. The results of these statistical significance tests are shown in Figs. D.1 and D.2. The distributions in i_{p} and Ω_{p} are flattened and no statistically significant perturbing orbits are produced.
Fig. D.1.
Barycentric orbital elements of putative perturber planets. As Fig. 1 but using scrambled data as input. 
Fig. D.2.
Heliocentric orbital elements of putative perturber planets. As Fig. C.1 but using scrambled data as input. 
Appendix E: Results at 200 AU
Fienga et al. (2020) focused on testing for the presence of possible planets at 400, 500, 600, 650, 700, 750, and 800 AU with masses of 5 or 10 M_{⊕}. The existence of 5 M_{⊕} planets at 400 or 500 AU is strongly disfavored by their results (see their Fig. 5, top panels). However, a hypothetical Earthlike planet at 200−400 AU from the Sun may still induce significant gravitational effects if close encounters are possible, due to its relatively large value for the Hill radius (e.g., 2.2 AU if a_{p} = 200 AU, e_{p} = 0.1 and 1 M_{⊕}). We repeated the analysis, imposing q_{p} > 200 AU, and we obtained 8234 orbits with a number of potential close approaches in the range between 5−7. Here, we count how many synthetic ETNOs had at least one mutual nodal distance with the planet under 2 AU. The median values and 16th and 84th percentiles (absolute maximum in parentheses) from the Monte Carlo random searches whose distributions are shown in Fig. E.1 are: AU (258 AU), (0.03), (10°), (196°), and (350°).
Fig. E.1.
Barycentric orbital elements of putative perturber planets. As Fig. 1 but for q_{p} > 200 AU. 
All Tables
Barycentric orbital elements and 1σ uncertainties of known ETNOs with robust orbit determinations.
All Figures
Fig. 1.
Barycentric orbital elements of putative perturber planets. Distributions of barycentric orbital elements of planetary orbits that may result in close encounters (under 5 AU for q_{p} > 300 AU, under 7.5 AU for q_{p} > 400 AU, and under 10 AU for q_{p} > 500 AU and q_{p} > 600 AU) with five or more presentday extreme transNeptunian objects (ETNOs). Each column of panels shows the cumulative results of 2 × 10^{10} experiments. From left to right: results of imposing q_{p} > 300 AU (left column, 20 304 orbits), q_{p} > 400 AU (5671 orbits), q_{p} > 500 AU (2635 orbits), and q_{p} > 600 AU (right column, 56 orbits). Median values are shown as vertical blue lines. Based on the available data on known ETNOs, the presence of massive perturbers well beyond 600 AU is strongly excluded within the context of the hypotheses centered in this study. 

In the text 
Fig. C.1.
Heliocentric orbital elements of putative perturber planets. Distributions of heliocentric orbital elements of planetary orbits that may result in close encounters (under 5 AU for q_{p} > 300 AU, under 7.5 AU for q_{p} > 400 AU, and under 10 AU for q_{p} > 500 AU and q_{p} > 600 AU) with five or more presentday extreme transNeptunian objects (ETNOs). Each column of panels shows the cumulative results of 2 × 10^{10} experiments. From left to right: results of imposing q_{p} > 300 AU (left column, 18 369 orbits), q_{p} > 400 AU (6594 orbits), q_{p} > 500 AU (3567 orbits), and q_{p} > 600 AU (right column, 105 orbits). Median values are shown as vertical blue lines. Based on the available data on known ETNOs, the presence of massive perturbers well beyond 600 AU is strongly excluded within the context of the hypotheses centered in this study. 

In the text 
Fig. D.1.
Barycentric orbital elements of putative perturber planets. As Fig. 1 but using scrambled data as input. 

In the text 
Fig. D.2.
Heliocentric orbital elements of putative perturber planets. As Fig. C.1 but using scrambled data as input. 

In the text 
Fig. E.1.
Barycentric orbital elements of putative perturber planets. As Fig. 1 but for q_{p} > 200 AU. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.