We highlight that the anomalous orbits of transNeptunian objects TNOs and an excess in microlensing events in the 5year Optical Gravitational Lensing Experiment data set can be simultaneously explained by a new population of astrophysical bodies with mass several times that of the Earth MSUBSUB. We take these objects to be primordial black holes PBHs and point out the orbits of TNOs would be altered if one of these PBHs was captured by the Solar System inline with the Planet 9 hypothesis. Capture of a free floating planet is a leading explanation for the origin of Planet 9 and we show that the probability of capturing a PBH instead is comparable. The observational constraints on a PBH in the outer Solar System significantly differ from the case of a new ninth planet. This scenario could be confirmed through annihilation signals from the dark matter microhalo around the PBH.
['High Energy Physics - Phenomenology', 'Astrophysics - Cosmology and Nongalactic Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics']
What if Planet 9 is a Primordial Black Hole
{'What if Planet 9 is a Primordial Black Hole?': "Jakub Scholtz 1 and James Unwin 2  \n1 Institute for Particle Physics Phenomenology, Durham University, Durham, DH1 3LE, United Kingdom 2 Department of Physics, University of Illinois at Chicago, Chicago, IL 60607, USA; & Department of Physics, University of California, Berkeley & Theoretical Physics Group, LBNL & Mathematics Sciences Research Institute, Berkeley, CA 94720, USA  \nWe highlight that the anomalous orbits of Trans-Neptunian Objects (TNOs) and an excess in microlensing events in the 5-year OGLE dataset can be simultaneously explained by a new population of astrophysical bodies with mass several times that of Earth ( M ⊕ ). We take these objects to be primordial black holes (PBHs) and point out the orbits of TNOs would be altered if one of these PBHs was captured by the Solar System, inline with the Planet 9 hypothesis. Capture of a free floating planet is a leading explanation for the origin of Planet 9 and we show that the probability of capturing a PBH instead is comparable. The observational constraints on a PBH in the outer Solar System significantly differ from the case of a new ninth planet. This scenario could be confirmed through annihilation signals from the dark matter microhalo around the PBH.  \n1. Introduction. As of this year, two gravitational anomalies of similar mass but very different origins remain to be explained. First, there is a growing body of observational anomalies connected to the orbits of trans-Neptunian Objects (TNOs) [1-3]. These observations have been taken as evidence of a new ninth planet in our solar system, called Planet 9 (P9), with mass M 9 ∼ 5 -15 M ⊕ and orbiting around the Sun at a distance of 300 -1000 AU [4]. Second, there is set of gravitational anomalies recently observed by the Optical Gravitational Lensing Experiment (OGLE). OGLE reported an excess of six ultrashort microlensing events with crossing times of 0 . 1 -0 . 3 days [5]. The lensing objects are located towards the galactic bulge, roughly 8kpc away. These events correspond to lensing by objects of mass M ∼ 0 . 5 M ⊕ -20 M ⊕ [6] and could be interpreted as an unexpected population of free floating planets (FFPs) or as Primordial Black Holes (PBHs).  \nIt is remarkable that these two anomalies correspond to a similar mass scale. Perhaps the most natural explanation is that they are caused by the existence of an unknown population of planets, i.e. the OGLE anomaly should be interpreted as due to FFPs denser than the local star population [7] and P9 might be one of those planets that have been captured by the Solar System. This would imply that our models for planet formation may need to be updated to account for this new population of FFPs, but the current program to hunt for P9 would go unchanged.  \nWe focus on a more exciting possibility: if the OGLE events are due to a population of PBHs then it is possible that the orbital anomalies of TNOs are also due to one of these PBHs that was captured by the Solar System. In this paper, we argue that this scenario is not unreasonable and discuss the observable implications; we estimate the probability of capture of a PBH by the Solar System, highlight that the observational constraints differ significantly between planets and PBHs, and point out that the dark matter (DM) microhalo, which generically forms around such a PBH, can lead to its discovery.  \n2. Two Anomalies. While the structure of the Solar System to semi-major axis a ∼ 100 AU is well explained, for a > 250 AU there are TNO populations whose orbits cannot be readily understood. Observations of TNOs, objects with a > 30 AU, and extreme TNOs (eTNOs) with a > 250 AU exhibit the following anomalies:  \n- i) Unexpected clustering in eTNO orbits [2, 3];\n- ii) The existence of high perihelia ( q ∼ 70 AU) TNOs, such as Sedna, collectively called Sednoids [1, 2, 8];\n- iii) TNOs moving roughly perpendicularly to the planetary plane (with inclination i glyph[greaterorsimilar] 50 · ) [9-11].  \nAn excellent review of TNO anomalies is given in [4].  \nSimulations and analytic arguments indicate that these observations are at odds with the predicted dynamics assuming only the known giant planets. For instance, any coincidence of initial orbits of eTNOs will disperse under evolution (on a timescale of a 10-100 million years [12]), and Solar System simulations imply inclinations are typically bounded by i glyph[lessorsimilar] 40 · [13, 14].  \nNotably, all of the TNO anomalies can be simultaneously explained by a new gravitational source in the outer Solar System. From observations of TNO dynamics one can infer the likely mass and orbit of this hypothetical source, referred to as Planet 9. Simulations of TNOs have identified a number of characteristic benchmark scenarios for a new gravitational source of mass 5 M ⊕ and 10 M ⊕ , which we summarise in Table I.  \nTABLE I. This table contains some benchmark P9 scenarios from [4] stating central values. Here a is the central major axis, e is the ellipticity, i is the inclination. Also given are the perihelion r p = r a (1 -e ) and the aphelion r a = 2 ra -r p . Deviations of O (10%) in parameters provide comparable fits.  \nThe second set of anomalies has been observed by OGLE. The OGLE observations, when interpreted as a population of PBHs, are consistent with a range of masses and density fractions [6] with  \nM ∈ [0 . 5 M ⊕ , 20 M ⊕ ] , f PBH ∈ [0 . 005 , 0 . 1] , (1)  \nwhere f PBH ≡ Ω PBH / Ω DM , in terms of the DM relic density Ω DM h 2 ≈ 0 . 1 with the PBH population forming a subcomponent with relic density Ω PBH . The masses and density fractions are correlated with larger masses corresponding to smaller f PBH .  \nWhile an M ⊕ object is too light to be an astrophysical black hole formed by stellar collapse, PBHs arise from over densities in the early universe [15, 16] and as a result can be substantially lighter than M glyph[circledot] . Formation of PBHs inferred from OGLE has been discussed in [17-19]. We note that PBHs arise from O (1) density fluctuations during radiation domination, due to an increase in the primordial power spectrum. An intriguing coincidence is that PBHs formed during radiation domination via a strong first order phase transition around the electroweak scale are expected to have mass of the same order as P9 with M BH ∼ M ⊕ (125 GeV /T ) 2 [20].  \nSince [6] prefers f PBH < 1, we assume DM particles account for the remaining fraction: f DM = 1 -f PBH ≈ 1. In that case, a PBH will accrete DM and form a microhalo [21-25]. Since the DM densities in these microhalos are typically very high, DM annihilations can be significantly enhanced, leading to potentially detectable signals as we will discuss in Sections 4 and 5.  \n3. Capture Probability. There are three alternative hypotheses for the origin of P9: a) P9 formed on its current orbit ('in situ'); b) P9 formed in the inner Solar System and has been up-scattered into its current orbit; or c) P9 has formed outside of the Solar System and has been captured. While all three scenarios are unlikely, they are still favourable compared to the chance alignment of TNOs [3]. In case of the in situ formation, at a ∼ 500 AU there is typically insufficient time and material to build an Earth mass planet [26-28]. The prospect of a planet forming near Uranus and Neptune before being scattered to its present orbit is low since in order to fall into a stable orbit the planet would need to be appropriately influenced by a passing star (or another mechanism) [4, 29]. The probability of capturing a free floating planet (FFP) is estimated to be similarly improbable, with estimates differing by orders of magnitude depending on assumptions [29-31].  \nWe will argue that while there is a low probability of capturing an Earth mass PBH, it is no more improbable than capturing an FFP of similar mass. The Solar System capture rate can be expressed as follows  \nΓ = 〈 σnv 〉 = ∫ n 0 F ( v + v glyph[circledot] ) d σ d v v d v , (2)  \nwhere F ( v ) and n 0 are the velocity distribution and ambient density of the objects to be captured, d σ/ d v is the differential capture cross-section and v glyph[circledot] ,r is the velocity of the Sun with respect to the rest frame of the objects to be captured.  \nThe differential cross-section (identical for PBHs and free floaters) is significantly suppressed for relative velocities larger than 0 . 25km / s [32], which is much smaller than other velocities in the calculation. As a result, the velocity dispersions in the integrand can be approximated by the zero-order value F ( v glyph[circledot] ). This allows us to cancel the common factor of d σ/ d v in the ratio of PBH to FFP capture rate, which is then well approximated by  \nΓ BH Γ FFP glyph[similarequal] n BH n FFP F PBH ( v glyph[circledot] , PBH ) F FFP ( v glyph[circledot] , FFP ) . (3)  \nWeassume that the PBH velocity distribution is the same as the DM velocity distribution given by the Standard Halo Model [33], with v glyph[circledot] , DM = 220km/s and velocity dispersion σ PBH = v glyph[circledot] / √ 2. The local density of PBHs is related to the local DM density ( ρ DM = 0 . 4 GeV / cm 3 ) and the fact that PBHs comprise a fraction f PBH of this local density:  \nn BH = f PBH ( ρ DM M BH ) ∼ 35pc -3 ( f BH 0 . 05 )( 5 M ⊕ M BH ) .  \nAs for the FFPs, there are two possibilities: While the Solar System could capture a planet when inside a starforming region, for which the FFP density may be as high as 200pc -3 [32], such stellar nurseries are highly disruptive environments. Hence a planet captured in this manner is quite likely subsequently stripped by interactions with nearby stars [27-29]. Instead we consider capture in the field, away from star forming regions. In the field, the FFP density (which we take to the be similar to the local star density) is much lower: n FFP ∼ 0 . 2pc -3 . However, the available time for capture is much longer and the survival probability of a captured object is effectively unity. We assume the FFP velocity dispersion is inherited from the stars in the thin disk: σ ∗ ∼ 40km / s [32]. Remarkably, with these parameters we arrive at:  \nΓ BH Γ FFP ∼ 1 × ( 0 . 2pc -3 n FFP )( 40km / s σ FFP ) 3 ( f BH 0 . 05 )( 5 M ⊕ M BH ) .  \nWe find that the rates are comparable and thus conclude that the probability that an FFP is gravitationally captured by the Solar System in ambient space is roughly comparable to capturing a 5 M ⊕ PBH with f PBH ∼ 0 . 05. Therefore, if one is willing to entertain that the possibility that the TNO orbits indicate a captured planet, it is plausible that the gravitational source in the outer Solar System could instead be a PBH (once we establish evidence for such a PBH population).  \nFinally, we note that gravitational capture normally occurs due to multi-body interactions or drag through  \nthe local environment which leads to energy dissipation. Interestingly, capture of an object with an extended halo presents a new mechanism for dissipation since DM particles will be shed during the encounter. This possibly improves the chance of capture, but understanding the complicated dynamics would require a dedicated study.  \n4. Dark Matter. OGLE [6] indicates f PBH glyph[lessmuch] 1. If the rest is taken up by the DM component, PBHs accrete dense DM microhalos. This is fortunate since DM annihilation provide a potential detection route and in the absence of DM it would be likely impossible to detect an M ⊕ PBH in the Solar System. To describe the DM profile around the PBH one needs to consider its initial configuration and subsequent evolution. We will discuss several characteristic scales that control the properties of the halo: the influence radius r in , truncation distances r t due to striping, and (if applicable) the radius at which DM annihilations shape the halo r max .  \nThe radius of influence r in corresponds to the region in which the PBH dominates the local gravitational potential (assuming the uniform background density approximation) and effectively appears as an O (1) density perturbation. As a result, r in corresponds to the radius that contains DM mass equal to the PBH mass:  \nM BH = 4 π 3 ρ ( t ) r 3 in ( t ) . (4)  \nThe density profile of a halo is relatively simple if the DM kinetic energy can be neglected relative to its potential energy. 1 This is typically the case for M PBH glyph[greaterorsimilar] M ⊕ for DM with mass m glyph[greaterorsimilar] 100 GeV [25], in which case the DM profile is of the form [34] (this profile is typical of self-similar secondary infall):  \nρ ( r ) = ρ eq 2 ( r eq r ) 9 / 4 , (5)  \nin terms of r eq ≡ r in ( t eq ) ∼ 220 AU × ( M BH / 5 M ⊕ ) 1 / 3 the radius of influence at matter-radiation equality, and ρ eq ≡ ρ ( t eq ) glyph[similarequal] 2 . 1 × 10 -19 g / cm 3 the density of Universe at matter-radiation equality.  \nIf the DM can annihilate, the inner DM density may be depleted implying a cross-section dependent region of constant density ρ max = m/ 〈 σv 〉 τ within a radius r max = r eq [ 〈 σv 〉 τρ eq / (2 m )] 4 / 9 where 〈 σv 〉 is the thermally averaged annihilation cross-section and τ is the age of the universe [23]. For the DM models we consider in this work r max is smaller than the Schwarzchild radius r BH and density plateau does not exist. This happens for  \n〈 σv 〉 < 1 . 4 × 10 -54 cm 3 / s ( m 100 GeV ) ( M BH 5 M ⊕ ) 3 / 2 . (6)  \nFurthermore, the DM halo will not have indefinite extent. It will be truncated due to stripping by the Milky Way (MW), by encounters with passing stars, and by the Sun (for a PBH bound to the Solar System). Since we consider a captured PBH, it should have a similar (galactic) orbit to the Sun, thus it will not make any especially close approaches to the center of the Galaxy. This implies very little truncation due to the Galaxy, for an initial halo of M initial ∼ 100 M ⊕ then r t, MW ∼ 10kpc.  \nTidal disruptions from meeting individual stars can be more significant. The tidal radius r t,glyph[star] is dominated by the closest approach to a star. Since the PBH is travelling near the peculiar velocity of the Sun (10km/s) with respect to the local rest, it has travelled around 100 kpc over 10 10 years and has passed O (10 5 ) stars. The typical spacing of stars is 1 pc, and so the distance of closest approach is r ∗ ∼ √ 10 -5 pc = 650AU. The tidal stripping radius due to star encounters is then  \nr t,glyph[star] ∼ r ∗ ( M initial 2 M glyph[circledot] ) 1 3 ∼ 34AU ( r ∗ 650 AU ) ( M initial 100 M ⊕ ) 1 3 .  \nInterestingly, r ∗ is of the same order as the inferred P9 semi-major axis. The halo mass, obtained by integrating to r t,glyph[star] ∼ 34 AU, for the profile in eq. (5) reveals that the total halo mass is typically O (50%) of a PBH of mass 5 M ⊕ glyph[lessorsimilar] M BH glyph[lessorsimilar] 10 M ⊕ , and so any further tidal stripping is controlled by the PBH mass. Once the PBH settles into an orbit around the Sun, tidal radius cuts off the DM halo at  \nr t, glyph[circledot] ∼ r p ( M BH 2 M glyph[circledot] ) 1 3 ∼ 8AU ( r p 400AU ) ( M BH 5 M ⊕ ) 1 3 ,  \nwhich contains DM mass of the order O (15%) of the mass of the PBH.  \n5. Annihilation Signals. On its own, a PBH of mass 5 M ⊕ has a Hawking temperature of 0 . 004 K, making it colder than the CMB, and since it's radius is r BH ∼ 5 cm, the power radiated by the PBH alone is minuscule. However, the DM halo around this PBH can, if annihilating, provide a powerful signal. Annihilations in the PBH halo at the position of P9 would make for a potential FERMILAT source. Moreover, the whole PBH population contributes to gamma ray diffuse emissions [21]. Indeed, for a non-negligible f PBH and DM with a thermal crosssection 〈 σv 〉 0 ∼ 3 × 10 -26 cm 3 / s (i.e. classic WIMP DM), the other PBHs surrounded by DM give a diffuse gamma ray flux that is strongly excluded [24, 25].  \nIn what follows, we take a DM model that generates observable signals, but is not yet excluded, and consider a 'freeze-in' DM candidate [35]. In freeze-in scenarios the DM abundance is initially negligible and is subsequently generated by highly suppressed interactions between the Standard Model and the DM, controlled by a tiny coupling λ glyph[lessmuch] 1, leading to a relic density Ω DM ∝ mY FI ∝ λ 2 .  \nFor a specific model we consider a DM particle χ with mass m coupled to a mediator state φ with mass M φ via a Lagrangian term g ¯ χχφ , and in addition we couple φ to some Standard Model operator with coupling λ . The DM relic density generated by freeze-in is parametrically  \nΩ DM glyph[similarequal] 0 . 2 ( m 100 GeV ) ( λ 6 × 10 -12 ) 2 ( 10 TeV M φ ) . (7)  \nWe assume an annihilation cross-section to Standard Model particles of the form (cf. Supplementary Material)  \n〈 σv 〉 glyph[similarequal] λ 2 g 2 m 2 πM 4 φ , (8)  \nwith the characteristic values in eq. (7) this implies a characteristic cross-section of order  \n〈 σv 〉 ch glyph[similarequal] 1 . 3 × 10 -56 cm 3 / s × ( g 10 -2 ) 2 . (9)  \nFor fixed λ , decreasing g reduces the annihilation crosssection while maintaining the DM abundance. More model details appear in the Supplementary Material and related literature, see e.g. [36-44]. We stress, however, that we present freeze-in as an example of a potentially discoverable scenario that is not currently ruled out.  \n5.a. Point Source Limits. The flux for this freeze-in scenario can be found from the DM annihilation rate  \nΓ = 4 π ∫ r 2 d r ( ρ ( r ) m ) 2 〈 σv 〉 . (10)  \nWe take 〈 σv 〉 glyph[similarequal] 〈 σv 〉 ch , then r max < r BH from eq. (6) and thus the profile of eq. (5) is appropriate. Cutting off the integral at the tidal stripping radius r t, glyph[circledot] we obtain  \nΓ = √ 3 ρ eq 8 πG 3 〈 σv 〉 m 2 = 10 20 s -1 ( 〈 σv 〉 〈 σv 〉 ch )( 100GeV m ) 2 . (11)  \nNote that the above result is independent of the PBH mass (this is not the case for r max > r BH ). Given the annihilation rate, the flux Φ γ and the energy flux Φ E of photons for an observer on Earth is then  \nΦ γ = κ 1 Γ 4 πr 2 9 , and Φ E = κ 2 Γ m 4 πr 2 9 ,  \nwhere κ 1 is the average number of photons per DM annihilation within the observable band of the experiment, while κ 2 is the fraction of energy of the DM annihilation converted into photons within the observable energy band. Both κ 1 and κ 2 are dependent on the mass and branching annihilation channel of the DM candidate and can be determined for example from [45]. To obtain characteristic bounds we set κ 1 ∼ 10 and κ 2 ∼ 1.  \nThe smallest detectable photon flux in the 8 year FERMI-LAT source catalog [46] was J2143.0-5501 with  \nΦ γ = 8 . 8 × 10 -12 photons / cm 2 / s in the band 100 MeV to 100 GeV, while the smallest energy flux in the same band was Φ E = 5 . 96 × 10 -13 erg / cm 2 /s due to 4FGL J1014.6+6126. Using these fluxes as upper bounds we calculate the maximum DM annihilation cross-section allowed by the photon flux limit, taking the rate of eq. (11) for r 9 glyph[similarequal] 400 AU we obtain the bound  \n〈 σv 〉 < 5 . 1 × 10 -56 cm 3 / s ( m 100GeV ) 2 ( γ flux) .  \nThe energy flux bound implies a comparable limit  \n〈 σv 〉 < 2 . 2 × 10 -55 cm 3 / s ( m 100GeV ) ( E flux) .  \nThat the two approaches give similar bounds is a coincidence and for different mass DM, these bounds will differ. Moreover, for more massive DM candidates these constraints weaken as the photons from annihilations become too energetic and different instruments, such as CTA [47], will play an important role.  \nWe note that during the 8 years of data taking the source would have moved by several degrees, making source identification much harder. However, the DM annihilation signals may be easier to spot because the space-time correlation in the signal would enhance the sensitivity of the search. As a result, it is necessary to run a dedicated study in order to determine if there are any candidates that match the P9 trajectory in the FERMILAT dataset, which we will return to in future work [48].  \n5.b. Diffuse Photons Limit. In addition to looking for the P9 PBH there is also an observational bound from the full population of PBHs, this has previously been explored in [21-25] for thermal cross-section WIMP DM.  \nFollowing [24, 25] we bound the diffuse gamma ray flux due to PBH Φ γ by translating the limits on decaying DM of mass m and decay rate Γ DM with the identification  \ndΦ γ d E ∝ f (1 -f 2 )Γ M BH = Γ DM m . (12)  \nFor DM with mass 10 GeV glyph[lessorsimilar] m glyph[lessorsimilar] 10 6 GeV the observational limit is  \nΓ DM glyph[lessorsimilar] 10 -28 s -1 (13)  \nassuming 100% decays to bb , varying by only O (1) over the mass range [49]. Requiring that the differential flux due to Γ is less than the above limit for f = 0 . 05 implies a bound on the rate of  \nΓ glyph[lessorsimilar] 3 . 7 × 10 23 s -1 ( M BH 5 M ⊕ ) ( m 100 GeV ) . (14)  \nComparing with the rate in eq. (11) with the characteristic cross-section of 〈 σv 〉 ch , it is seen that the diffuse bound is readily satisfied in the DM models we consider.  \n6. Exotica. Before closing, we note that alternative 'exotic' compact astrophysical bodies may explain these observations such as DM microhalos (without PBH) e.g. [50-54], bose stars [55], or DM stars [56, 57]. OGLE cannot distinguish between PBH, exotic stars, and planets, however DM microhalos are unlikely to produce OGLE's lensing events. Another possibility is a sizeable DM halo could be shredded during the capture leading a toroidal DM mass distribution around the Sun at ∼ 500 AU with total mass ∼ 10 M ⊕ , this realises the secular approximation (phase space averaged) for a compact object, and is similar to the proposed toroidal baryonic distribution of [58]. Notably, each of these scenarios implies different experimental signatures, distinct from those of a rocky or gas planet.  \n7. Conclusion. This letter highlights that anomalous orbits of TNOs and OGLE's short microlensing events could have the same origin and explores the intriguing scenario that they both arise due to a population of 5 M ⊕ PBHs. While the principal search strategies for a planet is to employ optical [59, 60] and infrared/microwave surveys [61], the signals could be very different for a PBH (or another exotic object). Thus, the PBH hypothesis expands the required experimental program to search for the body responsible for TNO shepherding and motivates dedicated searches for moving sources in x-rays, gamma rays and other high energy cosmic rays. Conversely, if conventional searches fail to find Planet 9 and the evidence for TNO anomalies continues to grow, the PBH P9 hypothesis will become a compelling explanation.  \nAcknowledgements. We thank Martin Bauer for comments on the draft. We are grateful for the hospitality and support of the University of Oxford and the Simons Center for Geometry and Physics (Program: Geometry & Physics of Hitchin Systems). Part of this work was performed at the Aspen Center for Physics (ACP), which is supported by National Science Foundation grant PHY-1607611; the participation of JS at the ACP was supported by the Simons Foundation. JS is also very grateful for the support from the COFUND Fellowship. JU gratefully acknowledges support from the National Science Foundation grant DMS-1440140 while in residence at MSRI during Fall 2019.", 'A. SIZE OF THE PBH': 'The Schwarzschild radius of a black hole is given by  \nr BH = 2 GM BH c 2 glyph[similarequal] 4 . 5cm ( M BH 5 M ⊕ ) . (15)  \nIn Figure 1 we provide an exact scale image of a 5 M ⊕ PBH. The associated DM halo however extends to the stripping radius r t, glyph[circledot] ∼ 8AU, this would imply a DM halo which extends roughly the distance from Earth to Saturn (both in real life and relative to the image).  \nFIG. 1. Exact scale (1:1) illustration of a 5 M ⊕ PBH. Note that a 10 M ⊕ PBH is roughly the size of a ten pin bowling ball.  \n<!-- image -->', "B. DARK MATTER FREEZE-IN VIA A Z '": "To demonstrate a viable freeze-in model with potentially detectable indirect detection signals we consider a specific implementation in which the Standard Model is supplemented by two new 'hidden sector' states, the DM particle χ with mass m and a mediator Z ' with mass M Z ' . We assume a sizeable interaction gZ ' χχ for g ∼ O (10 -2 ), with mass ordering M Z ' > 2 m , and that the hidden sector only couples to the Standard Model fermions f (with mass m f ) via λZ ' ff , which involves a tiny coupling λ ∼ 10 -12 .  \nHere we relabel the mediator Z ' (from φ in the main text) and highlight that it can naturally arise from a kinetic mixing portal [36] involving a new U(1) ' gauge symmetry via the Lagrangian term ∆ L = δB µν X µν where X is the Z ' field strength and B is the hypercharge field strength, where δ parameterises the mixing induced by states (with masses M and M ' ) charged under both U(1) ' and hypercharge δ ∝ g Y g X 32 π 2 log ( M 2 /M ' 2 ) . In motivated models this can lead to couplings appropriate for freezein [37, 38] and one can identify the freeze-in coupling λ glyph[similarequal] δg Y where g Y are Standard Model fermion hypercharge couplings. For further details on Z ' freeze-in portals see e.g. [39-43].  \nThe initial Z ' and χ abundances are negligible but populations are generated via the portal interaction λZ ' ff . The Z ' abundance evolves according to [35]  \n˙ n Z ' +3 n Z ' H glyph[similarequal] ∫ dΠ Z ' dΠ f dΠ f (2 π ) 4 δ 4 ( p Z ' -p f -p f ) |M| 2 Z ' → ff f eq Z ' (16)  \n- [1] M. E. Brown, C. Trujillo, and D. Rabinowitz, Discovery of a candidate inner Oort cloud planetoid , ApJ 617.1 (2004): 645. [astro-ph/0404456].\n- [2] C. Trujillo, and S. S. Sheppard, A Sedna-like body with a perihelion of 80 astronomical units , Nature 507.7493 (2014): 471.\n- [3] K. Batygin, and M. E. Brown, Evidence for a distant giant planet in the solar system , AJ 151.2 (2016): 22. [1601.05438].\n- [4] K. Batygin, F. C. Adams, M. E. Brown, J. C. Becker, The planet nine hypothesis , Phys. Rep. (2019). [1902.10103].\n- [5] P. Mr'oz, A. Udalski, J. Skowron, et al. No large population of unbound or wide-orbit Jupiter-mass planets Nature 548 (2017) 183 [1707.07634].\n- [6] H. Niikura, et al. Constraints on Earth-mass primordial black holes from OGLE 5-year microlensing events, Phys. Rev. D 99 (2019) no.8, 083503 [1901.07120].\n- [7] P. Mr'oz, et al. [The OGLE Collaboration] Two new freefloating or wide-orbit planets from microlensing , A&A, 622, A201 (2019) [1811.00441].  \nwhere f eq Z ' = exp( -E χ /T ) is the equilibrium distribution. The generated abundance of Z ' is [35]  \nY Z ' = 135 8 π 3 g S ∗ √ g ρ ∗ ( M Pl Γ Z ' M 2 Z ' ) (17)  \nwhere Γ Z ' is the partial width involving bath states.  \nSince we take M Z ' > 2 m and g ∼ O (10 -2 ) any Z ' produced decays promptly to DM pairs and thus the resulting DM relic density is Y χ ≈ 2 Y Z ' for λ glyph[lessmuch] g . This can be related to the observed quantity  \nΩ DM glyph[similarequal] 0 . 2 ( m 100 GeV ) ( λ 6 × 10 -12 ) 2 ( 10 TeV M Z ' ) , (18)  \nwhere we have parameterised Γ Z ' glyph[similarequal] λ 2 M Z ' / 8 π .  \nThe feeble coupling λ glyph[lessmuch] 1 is characteristic of freeze-in models and is actually required to satisfy the consistency condition Y DM ( T ) glyph[lessmuch] Y eq DM (T). If λ is too large the hidden sector enters equilibrium ( Y DM ∼ Y eq DM ) and then subsequently freeze-out dynamics sets the DM relic density. However, since λ is very small, direct and indirect detection is typically challenging. For the Z ' model the s-wave annihilation cross-section for m f glyph[lessmuch] m<M Z ' is [44]  \nσ 0 glyph[similarequal] ∑ f N c λ 2 2 π ( 2 g 2 χv m 2 ( M 2 Z ' -4 m 2 ) 2 + g 2 χa m 2 f M 4 Z ' ) , (19)  \nwhere the sum runs over Standard Model fermions and g χa and g χv are the axial and vector couplings to χ . Then eq. (8) is obtained from the above in the limit m glyph[lessmuch] M Z ' . In the main text we highlight that if f PBH is relatively large, indirect detection may be a viable route to search for this class of models in certain parameter ranges.  \n- [8] S. Sheppard, et al. A New High Perihelion Inner Oort Cloud Object , ApJ 157.4 (2019). [1810.00013].\n- [9] R. S. Gomes, J. S. Soares, and R. Brasser, The observation of large semi-major axis Centaurs: Testing for the signature of a planetary-mass solar companion , Icarus 258 (2015): 37-49.\n- [10] Gladman, B., et al. Discovery of the first retrograde transneptunian object , ApJ 697.2 (2009): L91.\n- [11] Y. T. Chen, et al. Discovery of a new retrograde transneptunian object: hint of a common orbital plane for low semimajor axis, high-inclination TNOs and centaurs , ApJ 827.2 (2016): L24. [1608.01808].\n- [12] M. E. Brown, Observational bias and the clustering of distant eccentric Kuiper belt objects , AJ 154.2 (2017): 65. [1706.04175].\n- [13] H. Levison, et al. Origin of the structure of the Kuiper belt during a dynamical instability in the orbits of Uranus and Neptune , Icarus 196.1 (2008): 258-273. [0712.0553].\n- [14] D. Nesvorny, Evidence for slow migration of Neptune from the inclination distribution of Kuiper belt objects.\n- ApJ 150.3 (2015): 73. [1504.06021].\n- [15] S. Hawking, Gravitationally collapsed objects of very low mass, MNRAS 152 (1971) 75.\n- [16] B. J. Carr and S. W. Hawking, Black holes in the early Universe, MNRAS 168 (1974) 399.\n- [17] Y. Tada and S. Yokoyama, Primordial black hole tower: Dark matter, earth-mass, and LIGO black holes, Phys. Rev. D 100 (2019) no.2, 023537 [1904.10298].\n- [18] C. Fu, P. Wu and H. Yu, Primordial Black Holes from Inflation with Nonminimal Derivative Coupling, [1907.05042].\n- [19] B. Carr, S. Clesse, J. Garcia-Bellido and F. Kuhnel, Cosmic Conundra Explained by Thermal History and Primordial Black Holes, [1906.08217].\n- [20] K. Jedamzik, Primordial black hole formation during the QCD epoch, Phys. Rev. D 55 (1997) 5871 [astroph/9605152].\n- [21] G. Bertone, A. R. Zentner and J. Silk, A new signature of DM annihilations: gamma-rays from intermediatemass black holes, Phys. Rev. D 72 (2005) 103517 [astroph/0509565].\n- [22] B. C. Lacki and J. F. Beacom, Primordial Black Holes as DM: Almost All or Almost Nothing, ApJ 720 (2010) L67 [1003.3466].\n- [23] Y. N. Eroshenko, DM density spikes around primordial black holes, Astron. Lett. 42 (2016) no.6, 347 [Pisma Astron. Zh. 42 (2016) no.6, 359] [1607.00612].\n- [24] S. M. Boucenna, F. Kuhnel, T. Ohlsson and L. Visinelli, Novel Constraints on Mixed Dark-Matter Scenarios of Primordial Black Holes and WIMPs, JCAP 1807 (2018) no.07, 003 [1712.06383].\n- [25] J. Adamek, C. T. Byrnes, M. Gosenca and S. Hotchkiss, WIMPs and stellar-mass primordial black holes are incompatible, Phys. Rev. D 100 (2019) no.2, 023506 [1901.08528].\n- [26] S.J. Kenyon, B.C. Bromley, Making planet nine: Pebble accretion at 250-750 AU in a gravitationally unstable ring , ApJ. 825, 33 (2016). [1603.08008].\n- [27] C. H. Heller, Encounters with Protostellar Disks. II. Disruption and Binary Formation , ApJ 455, 252 (1995).\n- [28] E. C. Ostriker, Capture and Induced Disk Accretion in Young Star Encounters , ApJ, 424, 292 (1994).\n- [29] G. Li, and F. C. Adams Interaction cross-sections and Survival Rates for Proposed Solar System Member Planet Nine ApJ 823.1 (2016): L3. [1602.08496].\n- [30] A. Mustill, S. Raymond, and M. Davies, Is there an exoplanet in the Solar system? MNRAS Lett. 460.1 (2016): L109-L113. [1603.07247].\n- [31] R. Parker, T. Lichtenberg, and S. Quanz, Was Planet 9 captured in the Sun's natal star-forming region , MNRAS Lett. 472.1 (2017): L75-L79. [1709.00418].\n- [32] N. Goulinski, and E. Ribak, Capture of free-floating planets by planetary systems. MNRAS 473.2 (2017): 15891595. [1705.10332].\n- [33] A. K. Drukier, K. Freese and D. N. Spergel, Detecting Cold DM Candidates, Phys. Rev. D 33 (1986) 3495.\n- [34] E. Bertschinger, Self-similar secondary infall and accretion in an Einstein-de Sitter universe , ApJ, Suppl. S. 58, 39 (1985).\n- [35] L. J. Hall, K. Jedamzik, J. March-Russell and S. M. West, Freeze-In Production of FIMP DM, JHEP 1003 (2010) 080 [0911.1120].\n- [36] B. Holdom, Two U(1)'s and Epsilon Charge Shifts, Phys. Lett. 166B (1986) 196.  \n- [37] K. R. Dienes, C. F. Kolda and J. March-Russell, Kinetic mixing and the supersymmetric gauge hierarchy, Nucl. Phys. B 492 (1997) 104 [hep-ph/9610479].\n- [38] T. Gherghetta, J. Kersten, K. Olive and M. Pospelov, The Price of Tiny Kinetic Mixing, [1909.00696].\n- [39] X. Chu, T. Hambye and M. H. G. Tytgat, The Four Basic Ways of Creating DM Through a Portal, JCAP 1205 (2012) 034 [1112.0493].\n- [40] X. Chu, Y. Mambrini, J. Quevillon and B. Zaldivar, Thermal and non-thermal production of DM via Z'portal(s), JCAP 1401 (2014) 034 [1306.4677].\n- [41] F. Elahi, C. Kolda and J. Unwin, UltraViolet Freeze-in, JHEP 1503 (2015) 048 [1410.6157].\n- [42] S. Heeba and F. Kahlhoefer, Probing the freeze-in mechanism in dark matter models with U (1) ' gauge extensions, [1908.09834].\n- [43] P. Foldenauer, Light dark matter in a gauged U (1) L µ -L τ model Phys. Rev. D 99 (2019) 035007 [1808.03647].\n- [44] A. Alves, A. Berlin, S. Profumo and F. S. Queiroz, DM Complementarity and the Z ' Portal, Phys. Rev. D 92 (2015) no.8, 083004 [1501.03490].\n- [45] M. Cirelli et al. , PPPC 4 DM ID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection, JCAP 1103 (2011) 051 [1012.4515].\n- [46] [Fermi-LAT Collaboration], Fermi Large Area Telescope Fourth Source Catalog, [1902.10045].\n- [47] M. Doro et al. [CTA Consortium], Dark Matter and Fundamental Physics with the Cherenkov Telescope Array, Astropart. Phys. 43 (2013) 189 [1208.5356].\n- [48] J. Scholtz and J. Unwin, work in progress.\n- [49] T. Cohen, et. al. γ -ray Constraints on Decaying Dark Matter and Implications for IceCube, Phys. Rev. Lett. 119 (2017) no.2, 021102 [1612.05638].\n- [50] J. Diemand, B. Moore and J. Stadel, Earth-mass darkmatter haloes as the first structures in the early Universe, Nature 433 (2005) 389 [astro-ph/0501589].\n- [51] V. Berezinsky, V. Dokuchaev and Y. Eroshenko, Remnants of dark matter clumps, Phys. Rev. D 77 (2008) 083519 [0712.3499].\n- [52] A. M. Green, S. Hofmann and D. J. Schwarz, The First wimpy halos, JCAP 0508 (2005) 003 [astro-ph/0503387].\n- [53] D. Zhang, Impact of Primordial Ultracompact Minihaloes on the Intergalactic Medium and First Structure Formation, MNRAS 418 (2011) 1850 [1011.1935].\n- [54] C. Blanco, M. S. Delos, A. L. Erickcek and D. Hooper, Annihilation Signatures of Hidden Sector Dark Matter Within Early-Forming Microhalos, [1906.00010].\n- [55] P. Jetzer, Boson stars, Phys. Rept. 220 (1992) 163.\n- [56] C. Kouvaris and N. G. Nielsen, Asymmetric Dark Matter Stars, Phys. Rev. D 92 (2015) no.6, 063526 [1507.00959].\n- [57] D. Curtin and J. Setford, How To Discover Mirror Stars, [1909.04071].\n- [58] A. Sefilian and J. Touma, Shepherding in a Selfgravitating Disk of Trans-Neptunian Objects , 2019, AJ 157, 59 [1804.06859].\n- [59] E. F. Linder, and C. Mordasini, Evolution and magnitudes of candidate Planet Nine , A&A 589 (2016): A134. [1602.07465].\n- [60] S. Ginzburg, R. Sari, and A. Loeb, Blackbody Radiation from Isolated Neptunes , ApJ 822.1 (2016): L11. [1603.02876].\n- [61] A. M. Meisner, et al. Searching for Planet Nine with coadded wise and neowise-reactivation images. ApJ 153.2 (2017): 65. [1611.00015]."}