MAGIC observations of NGC 4278. The first low-luminosity radio galaxy with compact jets detected at TeV energies

1. Problem and Scientific Context

Low-luminosity Active Galactic Nuclei (LLAGN) with compact radio jets are the most common class of radio-loud AGN in the local universe ($z \lesssim 0.1$). Despite their prevalence, they had remained undetected in the Very High Energy (VHE; $E > 100$ GeV) band until recently. The Large High Altitude Air Shower Observatory (LHAASO) recently reported the first TeV detection of such a source, NGC 4278, a nearby ($D_L \approx 16.4$ Mpc) radio-loud AGN powered by a radiatively inefficient accretion flow (RIAF).

The scientific challenge addressed in this paper is to:

• Investigate the high-energy properties of NGC 4278 during both its flaring and subsequent quasi-quiescent states using independent data.
• Determine if the TeV emission can be explained by standard leptonic jet models (Synchrotron Self-Compton, SSC).
• Understand the physical mechanisms driving the transition between flaring and quiescent states and estimate the jet power to see if it aligns with previous time-averaged estimates.

2. Methodology

The authors employed a multi-wavelength approach combining archival data from ground-based Cherenkov telescopes, space-based gamma-ray observatories, and X-ray telescopes.

• MAGIC Observations (VHE):

  • NGC 4278 lies $\sim 1^\circ$ from two known TeV blazars (1ES 1215+303 and 1ES 1218+304). The authors re-analyzed archival MAGIC data (2010–2024) taken while observing these blazars.
  • Data Processing: Used the MAGIC Analysis and Reconstruction Software (MARS). A Random Forest algorithm separated hadronic background from gamma-ray signals.
  • Analysis Strategy: Since NGC 4278 was not the primary target, the team performed an off-axis analysis, masking the two blazars to estimate background. They calculated significance and derived 95% confidence upper limits (ULs) for fluxes above 150 GeV.
  • Note: The MAGIC observation window did not overlap with the specific LHAASO flaring period (Aug 2021 – Jan 2022).

• Fermi-LAT (GeV):

  • Analyzed data from Oct 2022 to July 2024 (post-flare) using the Fermipy package and P8R3_SOURCE_V3 instrument response functions.
  • Performed a binned likelihood fit including Galactic diffuse emission and 4FGL-DR4 catalog sources.

• Swift-XRT (X-ray):

  • Analyzed three new snapshots (May–June 2024) to characterize the post-flare X-ray state.
  • Modeled the nuclear emission as an absorbed power law with fixed Galactic and intrinsic absorption ($N_H = 4 \times 10^{20} \text{ cm}^{-2}$).

• Spectral Energy Distribution (SED) Modeling:

  • Used the JetSeT code to model the broadband SED (Radio to TeV) for both the Flaring State (LHAASO active phase) and the Quasi-Quiescent State (post-flare).
  • Model Assumptions: A one-zone SSC model (spherical blob, relativistic plasma). External Compton (EC) was neglected due to the RIAF nature of the source (lack of strong external photon fields).
  • Optimization: Markov Chain Monte Carlo (MCMC) methods were used to find best-fit parameters, assuming a power-law electron distribution $n(\gamma) \propto \gamma^{-p}$.

3. Key Results

A. Observational Findings

• MAGIC Detection: No statistically significant gamma-ray emission was detected from NGC 4278. The significance was $0.7\sigma$.
• Upper Limits: The derived integral flux upper limit for the entire dataset ($>150$ GeV) is $F_{UL} = 1.5 \times 10^{-12} \text{ ph s}^{-1} \text{ cm}^{-2}$. These limits are consistent with LHAASO results, effectively constraining the spectrum down to 100 GeV.
• Fermi-LAT: No significant detection in the post-flare period ($TS = 2.5$, $\sim 1.2\sigma$).
• Swift-XRT: Detected X-ray emission with a photon index $\Gamma \approx 2.1$ and an unabsorbed flux of $(8 \pm 2) \times 10^{-13} \text{ erg s}^{-1} \text{ cm}^{-2}$ in the 0.5–8 keV band.

B. SED Modeling and Physical Interpretation

• Flaring State:

  • The X-ray emission is interpreted as synchrotron radiation, while the TeV emission is Inverse Compton (SSC).
  • Particle Dominance: The emitting region is extremely particle-dominated ($U_B/U_e \approx 2 \times 10^{-5}$), far exceeding typical equipartition values found in BL Lac objects.
  • Peak Frequencies: The synchrotron peak is inferred at extreme high frequencies ($\nu_s \sim 5 \times 10^{19}$ Hz), and the Compton peak at $\nu_C \sim 4 \times 10^{26}$ Hz.
  • Acceleration: An efficient acceleration mechanism is required to reach maximum electron Lorentz factors of $\gamma_{max} \sim 4 \times 10^7$.

• Quasi-Quiescent State:

  • The source remains detectable at TeV energies but with lower flux.
  • Transition Mechanism: A simple radiative cooling model (without particle re-injection) failed to reproduce the observed flux levels after 2 years; the predicted cooling was too rapid. This implies continuous or episodic particle injection is necessary to maintain the emission.
  • Parameter Shifts: The quiescent state suggests a lower magnetic field ($B \sim 0.7$ mG vs. $8$ mG) and a larger emitting region radius ($R \sim 4.7 \times 10^{16}$ cm vs. $1.3 \times 10^{16}$ cm). The synchrotron peak shifts to lower frequencies ($\nu_s \sim 8 \times 10^{17}$ Hz), suggesting transitional BL Lac behavior.

C. Jet Power Estimates

• The total jet power ($L_{jet}$) is estimated at $\sim 10^{42}$ erg s$^{-1}$ for both states.
• Composition: >99% of the power is kinetic, shared roughly equally between relativistic electrons and cold protons. Radiative losses account for <1%.
• Consistency: These estimates align remarkably well with previous time-averaged observational constraints derived from radio core luminosity and cavity power ($L_{jet} \sim 1-2 \times 10^{42}$ erg s$^{-1}$).

4. Significance and Conclusions

• Confirmation of LLAGN Physics: The study confirms that low-luminosity AGN with compact jets can accelerate particles to TeV energies, challenging the notion that such power is exclusive to high-luminosity blazars.
• Recurrent Flaring: The agreement between the modeled jet power during flares and the time-averaged jet power supports the hypothesis that NGC 4278 undergoes recurrent, episodic flaring events. These events are sufficient to explain the long-term energy budget without requiring a continuous, high-power jet.
• Jet Structure: The results support the view that these jets are composed of discrete, compact components (mini-jets) rather than a continuous flow, consistent with VLBI observations.
• Confinement: The jet power is modest ($L_{jet} \lesssim 10^{43}$ erg s$^{-1}$), suggesting the jets are likely confined by the dense nuclear environment of the host galaxy and cannot break out to form large-scale radio structures.
• Future Prospects: NGC 4278 represents a new class of targets for future VHE observatories like the Cherenkov Telescope Array (CTA), which will have the sensitivity to detect similar low-luminosity, compact radio galaxies.

In summary, this paper provides the first independent VHE constraints on NGC 4278, validates the LHAASO detection, and uses multi-wavelength modeling to characterize the extreme particle acceleration and episodic nature of jets in low-luminosity AGN.