🔭 astrophysics

Characterizing the Gamma-ray Emission from Low-Luminosity AGN

Using 14.4 years of Fermi-LAT data and a newly released Python stacking library, this study characterizes the gamma-ray emission of low-luminosity AGN, finding that individually detected sources are likely dominated by jet emission explainable by synchrotron self-Compton radiation, while the subthreshold population is consistent with star-formation activity.

Original authors: Chris Karwin, Nikita Khatiya, Margot Boughelilba, Xiurui Zhao, Anita Reimer, Marco Ajello

Published 2026-04-08

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Chris Karwin, Nikita Khatiya, Margot Boughelilba, Xiurui Zhao, Anita Reimer, Marco Ajello

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 the universe is a vast, bustling city. Most of the "power plants" in this city are Active Galactic Nuclei (AGN)—supermassive black holes at the center of galaxies that are feasting on gas and dust, shooting out massive beams of energy like cosmic lighthouses.

For a long time, astronomers have been looking at the brightest, loudest power plants in the city. But this new paper focuses on the quiet, low-power ones: the Low-Luminosity AGN (LLAGN). These are the black holes that are barely eating, glowing dimly, and moving slowly.

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

1. The Mystery of the Quiet Giants

The team, led by Christopher Karwin, wanted to know: Do these quiet black holes shoot out high-energy gamma rays (the most energetic light in the universe), and if so, why?

Previously, only four of these quiet giants had been spotted by the Fermi-LAT telescope (a giant gamma-ray camera in space). The researchers decided to look at 186 of these quiet sources at once. Since they are too faint to see individually, they used a clever trick called "stacking."

The Analogy: Imagine trying to hear a single whisper in a noisy stadium. You can't hear one person. But if you ask 186 people to whisper the same thing at the exact same time, their voices combine into a roar that you can hear. That's what "stacking" does: it combines the faint signals of many galaxies to reveal a pattern that was previously hidden.

2. The Big Discovery: Two Different Voices

When they "listened" to the combined signal of these 186 quiet galaxies, they found a clear answer, but it had two distinct voices:

Voice A: The Star-Forming Noise.
The researchers found that the signal from the quiet galaxies (especially the spiral-shaped ones like our Milky Way) looks exactly like the signal from star formation.
- The Metaphor: Think of a busy construction site. The gamma rays aren't coming from the black hole engine itself, but from the "construction crew" (new stars) building around it. The more stars being born, the more gamma rays are produced.
- The Result: For most of these quiet galaxies, the gamma rays are likely just a side effect of new stars being born, not the black hole doing anything special.
Voice B: The Jet Whisper.
However, there was also a faint hint of a connection to radio waves coming from the center of the galaxy.
- The Metaphor: This is like hearing a faint hum from a jet engine. It suggests that in some of these galaxies, the black hole is still shooting out a tiny, compact jet of particles, even if it's not a massive firehose like the loud AGN.
- The Result: It's a mix. The signal is mostly star formation, but a little bit of it might be the black hole's jet.

3. A New Star in the Sky

While stacking the faint signals, the team also found something new. One specific galaxy, NGC 4374, was finally loud enough to be seen on its own.

This brings the total number of confirmed "quiet" gamma-ray black holes to five. It's like finding a new friend in a crowd of strangers who finally spoke up.

4. How the "Quiet" Jets Work (The Physics Part)

For the few galaxies they could see clearly (NGC 315, NGC 4261, and the new NGC 4374), the team built a computer model to figure out how they are making gamma rays.

They found that these jets are very different from the famous, powerful ones (like in M87).

The Old Idea: Jets are like high-pressure fire hoses, moving near the speed of light, with strong magnetic fields.
The New Reality: These quiet jets are like a slow-moving, lazy river.
- They are moving slowly (not near light speed).
- They are "weakly magnetized" (the magnetic field is weak).
- They are "particle-dominated" (they are full of heavy particles, not just energy).

The Metaphor: Imagine a powerful sports car (the loud AGN) vs. a slow, heavy delivery truck (the quiet AGN). The sports car zooms by with a loud engine (strong magnetic fields). The delivery truck moves slowly, but it's so heavy with cargo (particles) that it still manages to make a little bit of noise (gamma rays) just by bumping into things.

5. Why This Matters

This paper is important because it solves a puzzle: Why are most nearby black holes so quiet?

It tells us that for the "quiet" ones, the gamma rays we see are mostly just the background noise of new stars being born.
But for the few that do have jets, those jets are slow, heavy, and weak.
The team also released their "stacking" software code to the public, so other scientists can use this "whisper-amplifier" tool to study other faint objects in the universe.

The Bottom Line

The universe is full of "quiet" black holes. They aren't the roaring monsters we often imagine. Instead, they are mostly just sitting there while new stars are born around them. But occasionally, they do let out a slow, heavy burp of energy from a tiny jet, proving that even the quietest giants have a little bit of fire left in them.

1. Problem Statement

Low-Luminosity Active Galactic Nuclei (LLAGN) constitute the majority of AGN in the local Universe ( $d_L \le 100$ Mpc), characterized by bolometric luminosities $\lesssim 10^{42}$ erg s $^{-1}$ and radiatively inefficient accretion flows (RIAFs). While theoretical models predict high-energy $\gamma$ -ray emission from both the accretion disks and relativistic jets of LLAGN, observational evidence has been scarce. Prior to this study, the Fermi Large Area Telescope (Fermi-LAT) had only significantly detected four LLAGN (NGC 1275, NGC 4486/M87, NGC 315, and NGC 4261).

The central scientific challenge addressed in this work is to characterize the $\gamma$ -ray emission mechanisms of the broader, undetected LLAGN population. Specifically, the authors aim to determine whether the $\gamma$ -ray signal in these subthreshold sources is dominated by:

Star-Formation Activity (SFA): Common in spiral galaxies, producing $\gamma$ -rays via cosmic ray interactions.
Compact Jets: Relativistic jets associated with the AGN, potentially producing $\gamma$ -rays via Synchrotron Self-Compton (SSC) processes.
A Mixed Origin: A combination of both mechanisms.

2. Methodology

Sample Selection

The study utilizes the Palomar spectroscopic survey of bright, northern galaxies (486 galaxies).

Target Sample: The authors selected Seyferts, LINERs, and transition nuclei, resulting in an initial sample of 186 subthreshold sources.
Exclusions: Sources already detected by Fermi-LAT (4 sources), a new detection (NGC 4374), sources matching blazar catalogs (Roma-BZCAT, CRATES, Yuan & Wang), and NGC 4151 (due to confusion with ultra-fast outflows) were removed or analyzed separately.
Control Sample: 65 galaxies from the Palomar survey with no visible emission lines were used as a control to verify the signal is not instrumental or background-related.

Data Analysis

Dataset: 14.4 years of Fermi-LAT data (August 2008 – January 2023).
Energy Range: 1–800 GeV for the main analysis; 0.3–1 GeV for low-energy checks.
Technique: A stacking analysis was performed using the publicly released fermi-stacking Python library. This technique aggregates the likelihood signals of many sub-threshold sources to detect an average population signal.
- Joint Likelihood Fit: Combined four Point Spread Function (PSF) event types (PSF0–PSF3) to optimize sensitivity.
- Background Modeling: Modeled Galactic diffuse emission (gll_iem_v07), isotropic background, and point sources from the 4FGL-DR3 catalog.
- Validation: The method was validated using Monte Carlo simulations and by analyzing "blank sky" regions, confirming no spurious detections.

Spectral Energy Distribution (SED) Modeling

For the individually detected sources (NGC 315, NGC 4261, and the new detection NGC 4374), the authors performed detailed broadband SED modeling:

Model: One-zone steady-state Synchrotron Self-Compton (SSC) leptonic model.
Constraints: Used archival X-ray data (Chandra) and radio VLBI kinematics to constrain jet viewing angles, bulk Lorentz factors, and magnetic field strengths.
Assumptions: The models assume a charge-neutralized jet with cold protons and relativistic electrons, focusing on the inner jet regions ( $\sim 10^2 - 10^4 r_g$ ).

3. Key Contributions

New Detection: The first significant detection of NGC 4374 (M84) in $\gamma$ -rays (TS = 31.4), increasing the number of firmly detected LLAGN to five.
Subthreshold Population Detection: A statistically significant detection of the stacked signal from 186 subthreshold LLAGN sources (TS = 31.2, $\sim 5.2\sigma$ ).
Decomposition of Signal Sources:
- Demonstrated that the stacked signal is dominated by spiral galaxies (TS = 39.7, $6.0\sigma$ ), while non-spirals (ellipticals/lenticulars) showed no signal.
- Established a robust correlation between the $\gamma$ -ray luminosity ( $L_\gamma$ ) and the infrared luminosity ( $L_{IR}$ ), consistent with the scaling relation for Star-Forming Galaxies (SFGs).
- Identified a correlation between $L_\gamma$ and core radio luminosity ( $L_{15\text{GHz}}$ ), suggesting a potential contribution from compact jets, though this remains ambiguous due to the radio-IR correlation in SFGs.
Public Release: The authors released the fermi-stacking Python library to facilitate future subthreshold source analyses.

4. Results

Stacking Analysis Results

Overall Signal: The full sample of 186 sources yielded a best-fit photon index of $\Gamma \approx 2.3$ and an integrated flux of $5.5 \times 10^{-12}$ ph cm $^{-2}$ s $^{-1}$ .
Morphology Dependence: The signal is almost exclusively driven by spiral galaxies (6.0 $\sigma$ significance). Non-spiral galaxies showed no significant signal ( $\sim 0.8\sigma$ ).
Correlations:
- $L_\gamma$ vs. $L_{IR}$ : A strong correlation was found ( $\alpha \approx 1.4-1.5$ ), matching the SFG scaling relation. This implies SFA is the dominant driver for the subthreshold population.
- $L_\gamma$ vs. $L_{radio}$ : A correlation exists ( $\alpha \approx 0.5-0.55$ ). While this could indicate jet emission, the authors note it may also be a secondary effect of the radio-IR correlation in SFGs.
- Three-Component Test: A joint fit including both IR and radio luminosities showed a slight improvement over the IR-only fit, suggesting a possible mixed contribution, but the significance was insufficient for a definitive conclusion.

Individual Source Modeling (NGC 315, NGC 4261, NGC 4374)

The broadband SEDs of these three FR I radio galaxies were successfully modeled using SSC emission.
Physical Parameters: The models require the inner jet emission regions to be:
- Weakly magnetized: Magnetic energy density is negligible compared to particle energy density ( $u_B / u_{particles} \sim 10^{-3}$ ).
- Particle-dominated: The total energy density is strongly dominated by particles (electrons and protons).
- Slowly moving: Bulk Lorentz factors are low ( $\Gamma \sim 1.0 - 1.5$ , corresponding to $\beta < 0.3$ ).
Implication: These sources are predicted to be faint at TeV energies due to the low bulk motion and weak magnetic fields.

5. Significance and Conclusion

This work provides the most comprehensive characterization of $\gamma$ -ray emission from the local LLAGN population to date. The key findings are:

Heterogeneity of LLAGN: The LLAGN population is not homogeneous. The $\gamma$ -ray signal in the subthreshold sample is primarily driven by star-formation activity in spiral galaxies, rather than AGN jet activity.
Jet vs. Disk: While the subthreshold population is SFA-dominated, the individually detected LLAGN (which are typically early-type galaxies with strong jets) are dominated by jet emission (SSC).
Physical Constraints: For the detected jetted LLAGN, the $\gamma$ -ray emission can only be explained by SSC models if the jets are slow, weakly magnetized, and particle-dominated. This challenges standard models that often assume highly relativistic, magnetically dominated jets for high-energy emission.
Methodological Advance: The successful application of stacking techniques to LLAGN demonstrates the power of deep Fermi-LAT observations combined with statistical methods to study faint astrophysical populations that are individually undetectable.

In summary, the study concludes that while compact jets are present in some LLAGN and contribute to the $\gamma$ -ray sky (especially in early-type galaxies), the bulk of the $\gamma$ -ray emission from the local LLAGN population arises from star-formation processes in spiral host galaxies.