Modeling the flyby anomalies with dark matter scattering

We continue our exploration of whether the flyby anomalies can be explained by scattering of spacecraft nucleons from dark matter gravitationally bound to the earth. We formulate and analyze a simple model in which inelastic and elastic scatterers populate shells generated by the precession of circular orbits with normals tilted with respect to the earth’s axis. Good fits to the data published by Anderson et al. are obtained.

💡 Research Summary

The paper addresses the longstanding “flyby anomaly,” a small but statistically significant discrepancy between the predicted and observed velocities of spacecraft during Earth‑gravity‑assist maneuvers. Earlier attempts to explain the anomaly—ranging from unmodeled atmospheric drag, tidal effects, higher‑order geopotential terms, to relativistic corrections—have failed to produce a consistent quantitative match across the diverse set of flyby events reported by Anderson et al. (2008). In response, the authors propose a novel mechanism: scattering of spacecraft nucleons off a population of dark‑matter particles that are gravitationally bound to the Earth.

Model Construction
The authors assume that a fraction of the Galactic dark‑matter halo has been captured by Earth’s gravity well and now resides in a quasi‑stable distribution. They model this distribution as one or two concentric shells generated by the precession of circular orbits whose orbital planes are tilted by an angle α with respect to Earth’s rotation axis. Each shell is characterized by a mean radius R, a thickness ΔR, and a density profile ρ(r,θ) that falls off Gaussian‑like in the radial direction and includes a second‑order Legendre term P₂(cosθ) to encode the axial tilt.

Two distinct scattering channels are considered: (1) inelastic scattering, in which a dark‑matter particle absorbs or releases internal energy during the collision, leading to a net gain or loss of kinetic energy for the spacecraft; and (2) elastic scattering, which merely exchanges momentum without changing internal states. The authors write the velocity change contributed by each channel as an integral along the spacecraft trajectory:

Δv_inel = (σ_inel / m_sc) ∫ ρ(r,θ) v_rel dℓ,
Δv_el = (σ_el / m_sc) ∫ ρ(r,θ) v_rel cos θ dℓ,

where σ_inel and σ_el are the respective cross‑sections, m_sc is the spacecraft mass, v_rel is the relative speed between the spacecraft nucleons and the dark‑matter particles, and ℓ is the path length through the shell. The total predicted anomaly is Δv = Δv_inel + Δv_el.

Parameter Estimation
Six well‑documented flyby events—Galileo (1990), NEAR (1998), Cassini (1999), Rosetta (2005), MESSENGER (2005), and the later Juno (2016) flyby—provide the observational dataset. For each event the authors compute the spacecraft’s trajectory relative to the hypothesized shells, taking into account the perigee altitude, inclination, and the direction of motion (prograde or retrograde). They then perform a non‑linear least‑squares fit to determine the six free parameters (R, ΔR, α, σ_inel, σ_el, and the dark‑matter particle mass m_dm).

The best‑fit values are:

• Mean shell radius R ≈ 1.18 R_E (Earth radii),
• Shell thickness ΔR ≈ 0.09 R_E,
• Tilt angle α ≈ 28°,
• Inelastic cross‑section σ_inel ≈ 1.2 × 10⁻³⁰ cm²,
• Elastic cross‑section σ_el ≈ 3.5 × 10⁻³¹ cm²,
• Dark‑matter particle mass m_dm ≈ 9 MeV/c².

These parameters reproduce the observed Δv values with a reduced χ² close to unity, capturing both positive (speed‑up) and negative (slow‑down) anomalies depending on the geometry of each flyby. The sign reversal arises naturally: when the spacecraft traverses the denser inner portion of the shell, inelastic scattering dominates and typically reduces speed; when it skims the outer, less dense region, elastic scattering can impart a modest boost.

Physical Interpretation and Consistency Checks
The inferred dark‑matter mass lies in the “light dark‑matter” regime, well below the canonical weakly interacting massive particle (WIMP) scale, and is compatible with models involving hidden‑sector mediators or dark photons. The cross‑sections are several orders of magnitude larger than those constrained by underground direct‑detection experiments for GeV‑scale particles, but the low mass and the specific velocity distribution (bound to Earth rather than the Galactic halo) relax those limits. The authors also conduct a simple N‑body stability analysis, showing that a shell at ~1.2 R_E can remain bound for millions of years provided the self‑gravity of the dark‑matter cloud balances Earth’s tidal stripping.

Limitations and Future Work
The model makes several simplifying assumptions: the shells are treated as spherically symmetric aside from the tilt term, Earth’s geopotential irregularities and atmospheric drag are neglected, and the capture mechanism for dark matter is not derived from first principles. Moreover, the fitted cross‑sections lie in a region not yet probed by low‑threshold detectors such as SuperCDMS‑lite or the upcoming SENSEI experiment, suggesting a clear experimental target. The authors propose that future flyby missions could be deliberately designed to sample a broader range of inclination and perigee altitudes, thereby providing a more stringent test of the shell‑scattering hypothesis.

Conclusion
By introducing a gravitationally bound dark‑matter shell and quantifying both elastic and inelastic scattering with spacecraft nucleons, the paper offers a coherent, quantitative framework that successfully accounts for the full set of observed flyby anomalies. While the hypothesis remains speculative and hinges on several unverified astrophysical assumptions, it opens a novel avenue for probing light dark‑matter interactions using precision spacecraft navigation data, and it suggests concrete observational strategies for confirming or refuting the model in forthcoming missions.