The Positron Puzzle

The paper addresses why the Milky Way exhibits a prominent 511 keV line with a bulge-to-disk luminosity ratio ~1, unlike most astrophysical tracers, and what sources and processes produce and annihilate the positrons responsible. Using five decades of measurements—especially INTEGRAL/SPI since 2002—the work frames three central questions: What do we see (morphology and spectrum), where do the positrons come from (source populations and particle processes), and why does it look as it does (propagation and annihilation environments). It highlights: (1) latest fluxes for bulge, disk, and Galactic center; (2) spectral evidence that e+ annihilate predominantly in warm neutral and warm ionized ISM phases with high positronium fraction; (3) the need to reconcile source distributions, propagation over 0.1–10 Myr and 0.1–10 kpc, and annihilation conditions to explain the observed image; and (4) the limitations of MeV gamma-ray imaging and spectroscopy that complicate interpretation.

The review compiles and critiques prior observational, theoretical, and modeling efforts: early balloon and OSSE detections; INTEGRAL/SPI’s refined bulge/disk measurements and high-resolution spectroscopy; studies of annihilation conditions (e.g., Jean et al. 2006; Churazov et al. 2005, 2011; Siegert et al. 2019a) indicating warm ISM phases and high positronium fractions; propagation models (Jean et al. 2009; Martin et al. 2012; Alexis et al. 2014; Higdon et al. 2009) that struggle to reproduce SPI morphology; tracer map correlations (Knoedlseder et al. 2005; Siegert 2017) suggesting a correlation with old stellar light rather than gas; and numerous source assessments: nucleosynthesis (26Al, 44Ti, 56Ni), SNe Ia escape fractions, novae, X-ray binaries (pair plasma), pulsars, Sgr A*, cosmic-ray secondaries, and dark matter. It also reviews critiques of previous constraints on injection energies from annihilation-in-flight based on simplistic spectral fits, and highlights instrument limitations (coded-mask insensitivity to isotropic emission, strong backgrounds).

As a synthesis/review, the paper employs: (1) maximum likelihood modeling of INTEGRAL/SPI count data with bulge, disk, and central components, assessing asymmetries and point-like emission near the Galactic center; (2) high-resolution spectral decomposition of the 511 keV line into narrow and broad components plus ortho-positronium continuum to infer annihilation media (warm neutral/ionized ISM) and positronium fraction; (3) template-fitting with a relative likelihood improvement metric (Δχ2) to compare 511 keV and ortho-positronium bands against multi-wavelength tracer maps (IR starlight, Hα, HI 21 cm, CO, dust, synchrotron, free-free, GeV π0 emission, 1.809 MeV 26Al), identifying which distributions best match SPI data; (4) back-of-the-envelope and population-based production rate calculations for candidate sources using radioactive decay laws (branching ratios, lifetimes, ejecta masses), XRB pair production and duty cycles, and pulsar pair production scaling; (5) transport considerations via collisional and collisionless regimes, and a general Fokker–Planck (cosmic-ray transport) framework including diffusion, advection, re-acceleration, energy losses, and sinks (annihilation, escape); (6) reassessment of annihilation in flight by embedding it within flexible diffuse continuum models (e.g., GALPROP-based components) rather than fixed power-laws; (7) forward-looking analysis strategies (spherical harmonic decompositions of future images, hierarchical multi-component modeling, handling time-variable Solar System foregrounds) to robustly disentangle astrophysical components.

• Imaging and fluxes: Using INTEGRAL/SPI, recent fluxes are FB ≈ (8.9–10.1)×10−4, FD ≈ (13.1–20.1)×10−4, and FGCS ≈ (0.6–1.2)×10−4 ph cm−2 s−1, with bulge-to-disk flux ratio ≈0.6±0.1 translating to a luminosity ratio LB/LD = 1.0±0.1 (assuming effective distances 8.2 kpc bulge, 6.5 kpc disk). The bulge shows an asymmetry (peaking near l ≈ −1°) and a point-like central component within SPI’s ~2.7° PSF (~400 pc).
• Spectroscopy and ISM conditions: The 511 keV line profile and ortho-positronium continuum imply dominant annihilation in warm neutral and warm ionized ISM phases (temperatures ~7000–40000 K; ionization 2–25%), with total positronium fraction ~92–97% and contributions from in-flight formation and thermal recombination; some broadening may arise from Galactic rotation.
• Tracer map comparison: Only near-IR starlight (1.25–4.9 μm; COBE/DIRBE) consistently improves fits for both 511 keV and ortho-positronium bands, indicating old stellar light correlates with annihilation sites more than gas tracers (HI, CO), dust, synchrotron, free-free, or GeV π0 maps.
• Source budgets (order-of-magnitude): 26Al (~2.8 M⊙ in Galaxy) yields ~3.2×10^42 e+ s−1; 44Ti from ccSNe (~10−4 M⊙ per SN; ccSN rate ~2/century) gives ~1.6×10^42 e+ s−1 (potentially doubled including SNe Ia); SNe Ia 56Ni escape is uncertain: canonical 0.5 M⊙ with ~3% escape yields ~5.1×10^42 e+ s−1, but constraints from SN2014J allow up to ~2×10^43 e+ s−1; classical novae contribute <~3×10^41 e+ s−1 by theory (1275 keV 22Na upper limits allow <5×10^42 e+ s−1); X-ray binaries/microquasars could contribute up to ~(3±3)×10^43 e+ s−1 assuming ~10^4 sources and duty cycle 10−3–10−2; pulsars ~5×10^42 e+ s−1; cosmic-ray secondaries ~(1–2)×10^42 e+ s−1; Sgr A* episodic activity could supply bulge-scale flux in bursts; DM and PBH scenarios face strong constraints from 511 keV and MeV continuum (Reticulum II limits; isotropic background bounds).
• Propagation: MeV positron transport remains uncertain; collisional (nearly ballistic) and collisionless (wave-particle) regimes predict different propagation lengths (up to kpc) and timescales (0.1–10 Myr). Thermalization in typical warm ISM occurs on ~0.1–1 Myr, after which advection with gas flow dominates annihilation sites.
• Annihilation in flight: Prior upper limits on injection energy (3–7 MeV) are not robust when diffuse continuum is modeled flexibly (e.g., GALPROP components). Depending on component choices and possible unobserved isotropic/halo annihilation, in-flight contributions could be larger without violating MeV data.
• Foregrounds and isotropic components: Time-variable Solar System albedo (asteroids, trojans, Moon) and potential isotropic halo/IGM annihilation are plausible and largely invisible to coded-mask instruments, potentially biasing morphology inferences.
• Overall: Multiple source classes likely contribute; the best spatial tracers are old stellar populations; instrument systematics and modeling degeneracies are key to resolve. Future Compton telescopes (e.g., COSI) and advanced analyses can disentangle components and test scenarios (e.g., globular cluster 511 keV from flares or MSP-related activity).

The findings collectively indicate that no single source class can uniquely explain the 511 keV morphology and flux. Nucleosynthesis (26Al, 44Ti, 56Ni) plausibly supplies a non-negligible baseline, but cannot alone reproduce the bulge prominence and total rate; compact objects (XRBs, pulsars, Sgr A*) could add significant bulge-concentrated positrons, yet suffer from large uncertainties in population statistics, duty cycles, escape fractions, and spectral states. The strongest spatial correlation with near-IR starlight suggests that annihilation sites follow the old stellar distribution, implying either in-situ production (e.g., intermittent stellar flares, some XRBs) or limited propagation from nearby sources into warm ISM phases. Transport physics remains a central uncertainty: kpc-scale propagation can smear initial source distributions into gas-like morphologies, but current data do not favor spiral arm signatures. Prior constraints on positron injection energies based on fixed power-law continua are model-dependent; allowing for physically motivated diffuse components and potential isotropic contributions relaxes these bounds. Instrumental limitations (coded-mask insensitivity to isotropic/slow gradients) and time-variable Solar System foregrounds can mimic or dilute astrophysical structures, necessitating careful treatment. The review argues for integrated, hierarchical modeling that ties source populations, injection spectra, transport, and annihilation physics to multi-wavelength and cosmic-ray data, enabling direct fits to raw gamma-ray data and robust discrimination among scenarios (including DM).

This review reframes the Positron Puzzle as a multi-component, system-level problem: the 511 keV emission likely arises from a combination of nucleosynthesis products, compact object pair production, and possibly intermittent stellar flares, with annihilation predominantly in warm ISM phases and morphology correlated with old stellar light. It highlights methodological pitfalls (overly rigid spectral models, inadequate handling of isotropic/foreground components) and emphasizes instrument-driven limitations that can bias conclusions. Key directions include: deploying wide-field Compton telescopes (e.g., COSI) to detect isotropic and extended components; constructing hierarchical, physically grounded Galactic models (sources, injection, transport, annihilation) constrained by radio-to-gamma data and cosmic-ray measurements; applying structural analyses (e.g., spherical harmonics) to future images; quantifying and modeling Solar System foreground variability; and targeted searches (e.g., globular clusters) to test flare or MSP-related scenarios. Solving the positron puzzle will, in turn, inform dark matter searches and stellar/ISM physics (e.g., coronal heating, superbubbles).

• Instrumental: MeV gamma-ray astronomy suffers from strong instrumental backgrounds, limited angular resolution, and coded-mask insensitivity to isotropic emission or shallow gradients; SPI exposure is non-uniform.
• Modeling: Significant uncertainties in source population sizes, spatial distributions, duty cycles, and escape fractions; degeneracies among diffuse emission components in the MeV band; reliance on template fits and reconstructed images that can bias morphology.
• Propagation physics: Poorly constrained low-energy positron transport (collisional vs collisionless, diffusion coefficients, advection), thermalization timescales vs gas dynamics, and annihilation in flight probability in realistic ISM conditions.
• Data gaps: No direct MeV detections of novae lines to calibrate yields; sparse or contentious annihilation features from XRBs; limited sensitivity to high-latitude/halo and isotropic components; uncertain time-variable Solar System foreground contributions.