How does blazar-boosted dark matter produce elastic and inelastic scattering signatures?

Why inelastic scattering changes the BBDM picture

Blazars are active galactic nuclei with a pair of extremely relativistic jets, one closely aligned with our line of sight. The paper works within the idea that dark matter particles near the blazar center can be accelerated to ultrarelativistic energies by scattering against high-energy particles in the jet, producing a boosted dark matter population. Earlier BBDM studies assumed a constant dark matter-proton cross section and considered only elastic scattering. This paper argues that those simplifications are not sufficient: blazar jets can reach luminosities around 10^48 erg/s and carry particle spectra from GeV to above PeV, so the energy dependence of the cross section matters, and inelastic processes become significant at such energies.

The distinction is not cosmetic. Elastic scattering keeps the scattering partners intact and dominates the low-energy part of the boosted flux. Inelastic scattering, treated here as DIS between dark matter and protons, becomes important at higher energies and lets the boosted dark matter spectrum extend into regimes that a pure elastic calculation does not reach. The paper's central claim is therefore methodological as much as phenomenological: the shape of the predicted BBDM flux depends on which scattering processes are included.

A vector mediator sets the scattering strength

The calculation is performed in a vector portal model in which dark matter interacts with Standard Model quarks through a massive vector mediator V. The differential cross section for dark matter-proton scattering depends on the mediator mass m_V and the dark matter mass m_chi. The paper focuses on light dark matter, m_chi less than or about 0.1 GeV, which is the regime where boosted dark matter scenarios are most attractive, and adopts two benchmark mediator masses, m_V = 1 GeV and m_V = 10 GeV.

For the inelastic channel, the DIS differential cross section is written in terms of the energy transfer and the Bjorken variables x and y, with the proton structure function F2 built from quark and antiquark parton distribution functions; the study uses the MSTW 2008 NNLO PDFs. The elastic and inelastic contributions are then combined into a differential cross section with respect to the dark matter kinetic energy after scattering. This is the object that feeds the boosted flux calculation, so the choice of mediator mass and the inclusion of DIS both propagate directly into the predicted event rates and the resulting limits.

Two blazars as the source benchmarks

The paper applies the framework to two sources, TXS 0506+056 and BL Lacertae, using source parameters such as redshift, luminosity distance, black hole mass, jet bulk Lorentz factor, viewing angle, proton spectral index, and jet proton luminosity. TXS 0506+056 is treated with a lepto-hadronic model and BL Lacertae with a hadronic model. The dark matter density near the central black hole is described by a spike profile, normalized so that the mass within a specified radial range matches the black hole mass, and the line-of-sight integral of that profile is factored into the flux calculation.

The boosted flux at Earth is then obtained by integrating the dark matter-proton differential cross section against the jet particle spectrum over the jet emission region, weighted by the dark matter line-of-sight integral and the source distance. The paper notes that the DIS treatment in this flux expression is an approximation because the scattering angle is not fixed by kinematics in inelastic scattering; it reports that for m_chi below about 100 MeV, the mass range of main interest, the difference from the exact result becomes negligible as m_chi decreases, and at m_chi = 10 MeV the two results are nearly identical.

DIS extends the flux to higher energies

For both sources, the paper compares the full calculation including elastic and inelastic scattering with a pure elastic scenario, for dark matter masses of 1 MeV, 10 MeV, and 100 MeV, with g_chi = g_q = 0.1 and m_V = 1 GeV. In the low-energy region the elastic contribution dominates. At higher energies the DIS contribution becomes significant and allows the BBDM flux to extend to much higher energies than the elastic-only case. This is the direct numerical basis for the paper's argument that inelastic effects cannot be omitted when describing the boosted dark matter flux.

The practical consequence is that the high-energy end of the spectrum, which is where the neutrino and gamma-ray signatures live, is populated by the inelastic channel. A calculation that stops at elastic scattering would therefore miss the part of the flux most relevant to indirect detection with high-energy observatories.

Pions turn DIS into neutrinos and gamma rays

In the inelastic channel, dark matter scattering on protons produces mesons and baryons, and the paper treats the process as chi + p -> chi + mesons + baryons. Because pions are the lightest mesons, they dominate the fragmentation, with production rates significantly exceeding those of heavier mesons and baryons, so the calculation includes only pion contributions. Charged pions decay through muons into neutrinos and antineutrinos, while neutral pions decay into two photons.

These decay chains are what connect the scattering calculation to observables. The paper uses a specially modified version of GENIE to compute the DIS process between boosted dark matter and protons, and compares the analytical DIS cross section with the numerical result. The resulting high-energy neutrino and gamma-ray fluxes are the signatures that the paper confronts with data from neutrino and gamma-ray observatories.

Combined limits surpass earlier BBDM constraints

The paper derives constraints on the dark matter-proton cross section sigma_chi_p by comparing its predictions with multi-messenger data from ground- and space-based observatories. The constraint figures combine limits from Super-Kamiokande, IceCube, and Fermi-LAT/MAGIC, with the IceCube curves shown for m_V = 10 GeV and m_V = 1 GeV, alongside limits from PandaX-4T, SENSEI, and big bang nucleosynthesis, and a grey region excluded by LEP searches for Z -> gamma V. The paper states that its results improve on previous limits from the constant cross section and pure elastic scattering scenarios by several orders of magnitude, and that the constraints from the high-energy neutrino and gamma-ray signatures of DIS significantly surpass those from traditional direct detection experiments.

The paper also uses LZ and Borexino as effective direct detectors for the boosted dark matter flux, and reports that the resulting constraints reach sensitivity comparable to the LEP Z -> gamma V search limits. The comparison that matters for a reader is therefore between three treatments: constant cross section with pure elastic scattering, energy-dependent elastic scattering, and the full elastic-plus-DIS calculation. The paper's claim is that the last of these produces the strongest and most complete constraints.

A possible dark matter origin for IceCube neutrinos

The paper highlights that dark matter-proton DIS is a novel source of high-energy astrophysical neutrinos, and that for energetic sources such as blazars this channel becomes especially important. It connects this to IceCube's 2017 observation of high-energy neutrinos from the direction of TXS 0506+056 and suggests that the event may come from dark matter. The paper presents this as a motivating possibility rather than a demonstrated identification, and it argues that a combined analysis of neutrino observations, the gamma-ray spectrum, and direct detection signatures offers a way to constrain blazar emission models and probe the particle nature of dark matter.

This is the paper's most speculative but also most distinctive claim. It does not assert that the IceCube neutrinos are dark matter in origin; it shows that the BBDM framework with DIS can produce high-energy neutrinos from blazars and that this interpretation can be tested against the same multi-messenger data used to set the cross section limits.

Where the conclusions stop

The paper lists several factors that could affect its results and that are left for future work: the complex environment at the blazar center may absorb high-energy gamma rays, the actual jet geometry may deviate from the idealized conical structure assumed, proton trajectories may be bent by strong magnetic fields near the central engine rather than following straight lines, and the precise launching point of the jet remains uncertain. These caveats qualify the flux predictions and the derived limits, and they are stated alongside the results rather than as a separate verdict on the study.

The paper also notes a closely related study that explored neutrino signals from inelastic dark matter-proton scattering in TXS 0506+056 and AP Librae. It describes its own contribution as more comprehensive because it includes neutrino signals together with constraints from direct detection experiments and high-energy photon observations. The priority a reader should take away is therefore the combination of energy-dependent scattering, DIS, and multi-messenger constraints in one framework, applied to two specific blazars, rather than any single number in isolation.

Sources used in this answer