🔭 astrophysics

Pulsar based modeling of point spread function of Fermi Large Area Telescope

By deriving improved Point Spread Function parameters from bright pulsar data that contradict the standard Fermi/LAT models, this study demonstrates that the previously reported extended emission around the blazar Mrk 501 is likely an artifact of PSF inaccuracies rather than a genuine astrophysical signal.

Original authors: J. Blunier, A. Neronov, D. Semikoz

Published 2026-06-26

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: J. Blunier, A. Neronov, D. Semikoz

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to take a picture of a single, bright star in the night sky. You want to know if there is a faint, fuzzy cloud of light surrounding that star. To do this, you need to know exactly how your camera blurs the image of a perfect point of light. In astronomy, this blur is called the Point Spread Function (PSF). If you don't know exactly how your camera blurs things, you might mistake a natural blur for a real cloud, or miss a real cloud because you think it's just a blur.

This paper is about fixing the "blur map" for the Fermi Large Area Telescope (LAT), a space telescope that looks at high-energy light (gamma rays).

Here is the story of what the authors found, explained simply:

1. The Problem: A Flawed Map

The scientists who built the Fermi telescope created a mathematical model (a set of rules) to predict how the telescope blurs light. They used computer simulations to build this map.

However, the authors of this paper suspected the map might be slightly wrong. They compared the official "blur map" against real observations of pulsars. Pulsars are like cosmic lighthouses—neutron stars that spin incredibly fast and flash beams of light toward Earth. Because they are so far away, they should look like perfect, tiny dots. If the telescope's blur map is accurate, the light from a pulsar should fit the map perfectly.

2. The Detective Work: Using Cosmic Lighthouses

To test the map, the authors used a clever trick. Pulsars flash on and off very quickly.

On-pulse: When the lighthouse beam is pointing at us.
Off-pulse: When the beam is pointing away.

The "off-pulse" time still has background noise (like static on a radio), but no light from the pulsar itself. By subtracting the "off-pulse" noise from the "on-pulse" signal, the scientists isolated the pure light of the pulsar. They then looked at how that pure light was spread out in the sky.

The Discovery: They found that the official computer-generated map did not match the real pulsar data. The real light was spread out differently than the map predicted. In some areas, the map was too sharp; in others, it was too fuzzy. It was like trying to use a map of a city that was drawn 20 years ago to navigate a city that has since been rebuilt.

3. The Solution: A New, Better Map

Because the official map was inaccurate, the authors created a new, "data-driven" map. Instead of relying on computer simulations, they used the actual measurements from the pulsars to build a new set of rules for how the telescope blurs light.

They tested this new map and found it fit the pulsar data perfectly, whereas the old map left big gaps and errors.

4. The Big Test: The Case of Mrk 501

Why does this matter? Recently, another group of scientists claimed to have found a giant, fuzzy cloud of light around a distant galaxy called Mrk 501. They thought this cloud was evidence of invisible magnetic fields floating in the empty space between galaxies.

The authors of this paper decided to re-analyze the data for Mrk 501, but this time, they used their new, corrected blur map instead of the old one.

The Result: When they used the new map, the "fuzzy cloud" disappeared.

With the old map: The data looked like there was a cloud.
With the new map: The data looked like just a single point of light with no cloud.

The authors concluded that the "detection" of the cloud was likely a false alarm. It wasn't a real cloud; it was just an artifact caused by using a slightly wrong blur map. The telescope's imperfections had tricked the previous analysis into seeing something that wasn't there.

Summary

The Analogy: The Fermi telescope is a camera. The "PSF" is the instruction manual for how that camera blurs images.
The Issue: The official manual was slightly wrong.
The Fix: The authors used real "lighthouses" (pulsars) to write a corrected manual.
The Outcome: When they used the corrected manual to look at a famous galaxy (Mrk 501), a previously reported "glow" around it vanished. The glow was never real; it was just a misunderstanding of how the camera works.

This paper serves as a reminder that in science, even the tools we use to measure the universe (like our "blur maps") need to be constantly checked against reality to ensure we aren't seeing things that aren't there.

Technical Summary: Pulsar-based Modeling of the Point Spread Function of the Fermi Large Area Telescope

Problem Statement
The sensitivity of searches for extended gamma-ray emission around extragalactic sources is fundamentally limited by the precision of the telescope's Point Spread Function (PSF). Inaccuracies in the PSF models provided by the Fermi Large Area Telescope (LAT) Instrument Response Functions (IRFs) can lead to false positive detections of source extension. Previous studies have demonstrated that early PSF modeling based on Monte-Carlo simulations failed to reproduce observed event distributions, leading to spurious claims of extended emission around blazars (e.g., Ando & Kusenko 2010). While the Fermi/LAT collaboration has updated their PSF models using improved Monte-Carlo simulations, validation studies currently rely primarily on the 68% containment radius. The authors argue that this metric is insufficient for assessing uncertainties in the context of extended emission searches, as different PSF profiles (e.g., Gaussian vs. exponential with power-law tails) can share the same containment radius but differ significantly in shape.

Methodology
The authors propose a data-driven approach to assess and improve the Fermi/LAT PSF models by comparing the analytical models in the IRFs against signals from bright pulsars (Vela, Geminga, and Crab).

Data Selection: The study utilizes the P8R3_SOURCEVETO_V3 event class, selected for its low residual charged cosmic-ray background. Data spans 17 years (March 2008 – January 2026) in the 0.1–300 GeV energy range.
Signal Isolation: To isolate the point-source signal from extended pulsar wind nebulae and diffuse Galactic emission, the authors employ a pulse-phase analysis. They subtract the "off-pulse" flux from the "on-pulse" flux using pulsar ephemerides from the 3rd Fermi/LAT pulsar catalog. This subtraction removes non-periodic background components while preserving the morphology of the point source.
Model Comparison: The authors fit the analytical PSF model (a sum of two King functions describing "core" and "tail" components) to the stacked pulsar data across various energy bins, off-axis angles, and event types (PSF0–PSF3, FRONT/BACK).
Parameter Optimization: They compare the fit quality ( $\chi^2$ ) of the standard IRF parameters against a re-fitted set where the PSF parameters ( $f_c, \sigma_c, \gamma_c, \sigma_t, \gamma_t$ ) are left free to vary.

Key Contributions and Results

Identification of PSF Inconsistencies: The analysis reveals that the PSF parameters provided in the Fermi/LAT IRF files are inconsistent with the observed pulsar data. In more than half of the energy and off-axis angle bins, the IRF parameters lie outside the 90% confidence range of the measurements derived from pulsar data. The discrepancies manifest as mismatches in both the PSF core and the PSF tail, with residuals reaching beyond 10% over wide angular ranges.
Development of an Improved PSF Model: By fitting the King function parameters directly to the pulsar data, the authors derive an improved set of PSF parameters that are consistent with observations. While the 68% containment radii of the original and revised models agree within a few percent, the shape of the revised model differs significantly, particularly in the tails.
Re-analysis of Mrk 501: The authors apply the revised PSF to a re-analysis of the blazar Mrk 501, a source previously reported (Webar et al. 2025) to exhibit extended emission indicative of an Intergalactic Magnetic Field (IGMF) with strength $B \sim 10^{-15}$ $B \sim 1 0^{- 15}$ G.
- Using the standard IRF PSF, the analysis reproduces a detection of extended emission with a significance of $\sim 4.57\sigma$ ($TS = 20.92$).
- Using the revised, pulsar-based PSF, the evidence for extended emission disappears. The test statistic drops to $TS = 4.08 $($ \epsilon = 0.40 \pm 0.12$), which is not statistically significant.
- The revised PSF yields a better fit for the point source signal itself (increasing the point source $TS$ from 49,375 to 49,579).
- Consequently, the extended emission model corresponding to the IGMF parameters claimed by Webar et al. (2025) is ruled out at the 95% confidence level when using the revised PSF.

Significance and Limitations
The paper demonstrates that precise PSF modeling is critical for distinguishing between true extended emission and artifacts caused by PSF model mismatches. The authors conclude that the previously reported detection of extended emission around Mrk 501 is likely a false positive resulting from the use of PSF parameters that do not accurately describe the real data.

The authors note a specific limitation of their data-driven approach: it relies on high-statistics pulsed signals from pulsars, which are only available below $\sim 10$ GeV. Above this energy, pulsar fluxes are suppressed, necessitating a return to Monte-Carlo-based PSF determination for high-energy searches. However, for the 0.1–10 GeV range, the study establishes that the current IRF parameterizations require revision to avoid systematic errors in extended emission searches.

1. The Problem: A Flawed Map

2. The Detective Work: Using Cosmic Lighthouses

3. The Solution: A New, Better Map

4. The Big Test: The Case of Mrk 501

Summary

More like this