Oort Cloud Bombardment by Dark Matter

The realization that primordial black holes (PBHs) might be some fraction of the dark matter begged the question, how often do PBHs enter the solar system? For a Neptune radius solar system the answer is, rarely. For an Oort cloud sized system the answer is different. Simulations of bombardment of the Oort cloud by dark matter suggest that dislodgement of protocomets and their entry into the inner solar system can match the observed frequency of comets, if that PBH fraction is high enough. Comets were traditionally considered as messengers, usually omens. After 50 years of puzzlement regarding dark matter, we need a hint from the dark universe about the size and nature of dark matter particles.

💡 Research Summary

The paper “Oort Cloud Bombardment by Dark Matter” investigates whether macroscopic dark‑matter (DM) objects—specifically lunar‑mass (≈10⁻⁷ M☉) primordial black holes (PBHs) or similar free‑floating moons (FFM)—can perturb the distant Oort cloud enough to send comets into the inner Solar System at rates comparable to those observed. The authors begin by reviewing the conventional picture: the Galactic tidal field (especially the vertical component) and passing stars are the dominant external forces shaping the Oort cloud and injecting long‑period comets. They then argue that if a non‑negligible fraction of the Galactic dark halo consists of compact, massive objects rather than elementary particles (WIMPs, axions), an additional perturbation channel should exist.

To quantify this, the authors construct a “toy model”. They adopt a total Oort‑cloud comet population of N≈4×10¹³ objects larger than 10 km, based on size‑frequency distributions (α≈3.6 for D>2.8 km, α≈0.5 for smaller bodies). The spatial distribution follows the Dehnen (1993) density law n(r)∝a r²/(r+a)² with a scale radius a=10⁵ AU, but they also compare alternative power‑law profiles. In the simulations they populate the cloud with 250 000 test comets (≈6×10⁻⁹ of the full population) and let a stream of DM objects pass through. Each DM object has a mass M (varied from 10⁻¹⁰ to 10⁻⁵ M☉) and a velocity drawn from a combination of the Sun’s orbital speed (220 km s⁻¹) and an isotropic halo dispersion (110 km s⁻¹). The impact parameter b is assumed uniform within a maximum value b_max≈M AU (i.e., larger DM masses have proportionally larger cross‑sections). The momentum transferred to a comet during a close encounter is approximated by Δp=2GMm/(b v), where m is the comet mass and v the relative speed.

The integration scheme is simple: at each timestep (0.1–0.3 yr) the comet’s velocity is updated by the instantaneous acceleration, and its position is advanced accordingly. The authors record any comet whose perihelion falls within 300 AU, deeming it “delivered” to the inner Solar System. Because the simulation neglects planetary perturbations, inter‑comet gravity, and encounters with b>b_max, they introduce two correction factors when converting the raw delivery count (n) into a physical delivery rate (rate₂). The first factor accounts for the fraction of DM trajectories omitted due to the finite b_max (∝4a²/b_max²/N). The second rescales the simulated comet sample to the full Oort‑cloud population. The raw DM flux through the cloud is given by rate₁=4πρa²v/m, where ρ≈0.01 M☉ pc⁻³ is the local dark‑matter density and m the DM particle mass. Substituting the numbers yields a DM encounter rate roughly ten million times larger than the stellar encounter rate, because the number density of lunar‑mass objects is enormous even though each carries little mass.

Figure 5 (referenced in the text) shows the resulting delivery rate as a function of DM mass and b_max. For DM masses between 10⁻⁹ and 10⁻⁵ M☉, and for plausible b_max values, the model predicts delivery rates ranging from ≲0.1 yr⁻¹ up to ≳10 yr⁻¹. If only a fraction f_PBH of the total dark halo consists of such objects, the rates scale linearly with f_PBH. The authors find that f_PBH≈0.1 (i.e., 10 % of the halo mass in lunar‑mass PBHs) yields a delivery rate of order 1–2 comets per year, which matches the observed influx of long‑period comets (≈1.5 yr⁻¹ for hyperbolic orbits, potentially double that when near‑parabolic comets are included).

The paper discusses several caveats. First, the omission of planetary scattering likely underestimates the true delivery efficiency, as giant planets can both eject comets and funnel them inward. Second, the choice of b_max is somewhat arbitrary; a more realistic treatment would integrate over the full range of impact parameters, which would increase computational cost but could change the rate substantially. Third, the simulation’s statistical sample (250 k comets) is tiny compared with the real cloud, leading to Poisson noise in the delivery count. Fourth, the velocity distribution of DM is simplified; a full Galactic model would include anisotropies, substructure, and possible streams. Finally, the model assumes all DM objects have the same mass, whereas a realistic mass spectrum could alter both the encounter rate and the momentum transfer.

In the discussion, the authors turn to observational prospects. They argue that upcoming wide‑field surveys—Rubin’s LSST and the Roman Space Telescope—could detect microlensing events caused by lunar‑mass PBHs if the cadence is optimized for short‑timescale events (≈15 min). Current LSST cadence (low cadence, hemispheric coverage) is not ideal; a dedicated high‑cadence, moderate‑field program would improve sensitivity. They note that the limiting detectable mass for microlensing is ≈10⁻¹¹ M☉, and that moving to the ultraviolet could lower this limit by a factor of two due to reduced wave‑optics effects. The authors also point out that their code (comet6.f) and data are publicly available, facilitating reproducibility.

The conclusion emphasizes that macroscopic dark matter in the lunar‑mass range could plausibly contribute a non‑negligible fraction of the observed long‑period comet flux, provided that about 10 % of the halo mass resides in such objects. This would not only link dark‑matter physics to cometary dynamics but also hint at a broader role for PBHs (or similar compact objects) in the early Galaxy, possibly influencing the delivery of water and other volatiles to the inner Solar System. The paper calls for more sophisticated N‑body simulations that include planetary perturbations, a realistic spectrum of DM masses, and larger impact‑parameter coverage, as well as targeted microlensing campaigns to either detect or rule out a substantial lunar‑mass PBH component.