🔭 astrophysics

X-Ray Spectral Variability of the TeV HBL Blazar PG 1553+113 with XMM-Newton

This study analyzes 30 XMM-Newton observations of the TeV HBL blazar PG 1553+113 from 2001 to 2024, revealing that its X-ray spectra are predominantly well-described by Log-Parabola models with significant variability in spectral parameters, suggesting that spectral evolution is driven by changes in particle acceleration or cooling within the jet.

Original authors: P. U. Devanand, Alok C. Gupta, Paul J. Wiita, V. Jithesh, Archana Gupta

Published 2026-05-01

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: P. U. Devanand, Alok C. Gupta, Paul J. Wiita, V. Jithesh, Archana Gupta

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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 filled with cosmic lighthouses called blazars. These aren't ordinary lighthouses; they are powered by supermassive black holes at the centers of distant galaxies, shooting out beams of energy so powerful and fast (close to the speed of light) that they look incredibly bright to us, even though they are billions of light-years away.

One of these cosmic beacons is named PG 1553+113. It's a "High-energy peaked BL Lac object," which is a fancy way of saying it's a blazar that shines brightest in high-energy X-rays.

This paper is like a detective story where astronomers acted as cosmic detectives, using a space telescope called XMM-Newton to take a massive collection of "snapshots" (spectra) of PG 1553+113 over 23 years (from 2001 to 2024). They wanted to understand how the light from this object changes over time and what that tells us about the physics happening inside the jet.

Here is the breakdown of their findings in simple terms:

1. The Shape of the Light Curve

When you look at the light from this blazar, it doesn't look like a straight line. It's more like a curved hill.

The "Log-Parabola" (The Curved Hill): For 14 out of the 30 snapshots, the light curve was clearly curved. The authors found that a mathematical shape called a "Log-Parabola" fit these curves perfectly.
- Analogy: Imagine rolling a ball down a hill. Sometimes the hill is a straight ramp (a simple Power Law), but often, the hill curves, getting steeper or flatter as you go. This curvature tells us that the particles (electrons) inside the jet are being accelerated and cooled in a complex way, like a crowd of runners where the faster ones start to get tired and slow down at different rates.
The "Power Law" (The Straight Ramp): For the other 15 snapshots, the light curve was simple and straight. A "Power Law" model was enough to describe it.
- Analogy: In these moments, the "hill" was so straight that the curvature was too subtle to measure, or the signal was too faint to see the bend. It's like looking at a very long, straight road where the curve is too far away to notice.
The "Broken" Hill: In a few specific cases, the curve wasn't just a smooth hill; it looked like a hill that suddenly changed direction. This suggested that two different types of light were mixing together.

2. Finding the "Peak" of the Energy

Every blazar has a "peak" energy where it shines the brightest. For PG 1553+113, this peak is usually in the Ultraviolet (UV) or soft X-ray range.

The Problem: The X-ray telescope they used (XMM-Newton) couldn't see the very low-energy part of the light where the peak usually sits. It was like trying to find the top of a mountain while standing in a valley; you can see the slope, but not the very top.
The Solution: They combined the X-ray data with data from the telescope's "Optical Monitor" (which sees visible and UV light). By stitching these two views together, they could finally see the top of the mountain.
The Result: They found the peak energy was located in the Ultraviolet range (between 4.6 and 48.6 electron-volts). This confirms that the light is coming from electrons spiraling in magnetic fields (synchrotron radiation).

3. The "Ghost" in the Machine (Inverse Compton)

The most exciting discovery happened in three specific observations taken in February 2017.

The Mystery: In these snapshots, the light curve showed a strange "hardening" or a second bump at the high-energy end. It looked like the smooth hill had a second, smaller hill attached to its side.
The Explanation: The authors believe this is a sign of a second process called Inverse Compton scattering.
- Analogy: Imagine the main light (Synchrotron) is a stream of tennis balls. Occasionally, these tennis balls get hit by a fast-moving bat (high-energy photons) and get kicked up into a much higher energy state. This creates a "ghost" signal—a second, fainter stream of light riding on top of the main one.
The Verdict: Two of these observations (Obs A and Obs B) showed very strong evidence of this "ghost" light. A third one (Obs C) showed a very weak hint of it. This is rare for this type of blazar, which usually only shows the main "tennis ball" stream.

4. What Does This Mean?

The study concludes that PG 1553+113 is a dynamic object that changes its behavior over time.

Sometimes it behaves simply (a straight line).
Sometimes it shows complex curvature (a curved hill), suggesting changes in how particles are accelerated or cooled.
Occasionally, it reveals a hidden second component (the "ghost" light), suggesting that the jet is temporarily mixing two different types of energy production.

The authors also noted that when the blazar gets brighter, its light spectrum gets "harder" (shifts toward higher energies), which is a common behavior for this class of cosmic objects.

In summary: By taking 30 snapshots over two decades, the team mapped out the changing "personality" of PG 1553+113. They confirmed that its light usually follows a curved pattern, but occasionally, it reveals a hidden, complex structure where two different energy processes are happening at once. This helps scientists understand the extreme physics of the jets shooting out from supermassive black holes.

1. Problem Statement and Motivation

The blazar PG 1553+113 is a high-energy-peaked BL Lac object (HBL) and a TeV emitter, known for its intense flux variability and a claimed quasi-periodic oscillation (QPO) of ~2.2 years across multiple bands. While previous studies have analyzed individual observations or small subsets of data from various satellites (BeppoSAX, Swift, Chandra, etc.), there was no comprehensive study utilizing the entire archive of XMM-Newton observations available up to late 2024.

The primary scientific goals were to:

Characterize the X-ray spectral variability of PG 1553+113 over a 23-year baseline.
Determine whether the spectra are best described by simple Power-Law (PL) models or curved Log-Parabolic (LP) models, which imply stochastic particle acceleration.
Investigate the presence of a second emission component (Inverse Compton) manifesting as spectral hardening (concave curvature) in the X-ray band.
Estimate the synchrotron peak frequency ( $\nu_s$ ) by combining X-ray data with simultaneous Optical/UV data from the XMM-Newton Optical Monitor (OM).

2. Methodology

Data Selection and Reduction:

Dataset: The authors analyzed 30 EPIC-PN X-ray spectra derived from 37 available XMM-Newton observations (excluding 5 unpublished and 8 without EPIC science products). One observation (Obs ID: 0790380501) was split into scheduled and unscheduled components.
Energy Range: Spectral analysis was restricted to 0.6 – 7.0 keV. The lower limit was set to 0.6 keV to avoid instrumental uncertainties in timing mode data and the upper limit to 7.0 keV because background counts became comparable to source counts at higher energies in many observations.
Reduction: Data were processed using the Science Analysis System (SAS) v22.1.0. Pile-up corrections were applied for Prime Full Window (PFW) modes, and background regions were carefully selected.
Joint Fitting: For 27 observations with simultaneous Optical Monitor (OM) data, joint spectral fitting was performed using X-ray (EPIC-PN) and Optical/UV (OM) data. A 3% systematic error was added to OM data to account for cross-normalization issues.

Spectral Modeling:
The spectra were fitted using XSPEC with three primary models, all convolved with Galactic absorption (tbabs):

Power Law (PL): $F(E) = N E^{-\Gamma}$
Log-Parabola (LP): $F(E) = N (E/E_1)^{-(\alpha + \beta \log(E/E_1))}$ , where $\alpha$ is the photon index at pivot energy $E_1=1$ keV, and $\beta$ is the curvature parameter.
Broken Power Law (BPL): Used for spectra showing significant negative curvature (spectral hardening) to test for an additional Inverse Compton (IC) component.

Statistical Criteria:

Model Selection: An F-test was used to compare PL and LP models. If the null hypothesis (PL is sufficient) was rejected with $p \le 0.01$ , the LP model was preferred.
IC Signature Search: Candidates for an IC component were identified if they showed significant negative curvature ( $|\beta| \ge \beta_{err}$ ) in LP fits and were better described by a stable BPL model with independent soft ( $\Gamma_1$ ) and hard ( $\Gamma_2$ ) photon indices at the 99% confidence level.

3. Key Contributions and Results

A. Spectral Fitting Statistics (30 Observations)

Log-Parabolic Dominance: 14 out of 30 spectra were best fitted by the LP model with 99% confidence.
- Parameters: $\alpha \simeq 2.13 - 2.80$ and $\beta \simeq 0.04 - 0.18$ .
- One spectrum (Obs ID: 0790381401) showed a preference for LP with negative curvature ( $\beta \simeq -0.08$ ).
Power Law Dominance: 15 spectra were adequately described by simple PL models with photon indices $\Gamma \simeq 2.53 - 2.69$ . This suggests that during these epochs, the observed energy band sampled the linear portion of the synchrotron continuum where curvature was weak or statistically unconstrained.
Correlations:
- No significant correlation was found between curvature ( $\beta$ ) and flux, nor between $\alpha$ and $\beta$ . This implies the spectral curvature is not governed solely by a simple energy-dependent acceleration scenario but likely involves competing effects of stochastic acceleration and radiative cooling.
- A strong negative correlation ( $r_{lin} = -0.744$ ) was found between the photon index $\alpha$ and flux ( $F_{0.6-7}$ ). This "harder-when-brighter" behavior is typical of HBLs and suggests the injection of fresh, harder electrons during high-flux states.

B. Synchrotron Peak Estimation (Joint PN+OM)

X-ray data alone could not constrain the synchrotron peak ( $E_p$ ) as it fell below the 0.2 keV limit.
Joint fitting with OM data (using a reddening-corrected LP model, LPR) allowed the estimation of the rest-frame synchrotron peak energy.
Results: For ~20 observations, $E_p$ was constrained in the range 4.59 – 48.61 eV (UV regime). This aligns with previous SED studies placing the peak between 8–33 eV.
For some observations, the peak appeared near the optical V-band edge, making it difficult to constrain with OM data alone.

C. Search for Inverse Compton (IC) Signatures
The authors identified three candidate observations (Obs IDs: 0790381401, 0790381501, 0790381001) showing potential spectral hardening indicative of an IC component overlapping with the synchrotron tail.

Strong Candidates (Obs A & B):
- Obs A (0790381401, Feb 13, 2017): Showed significant negative curvature in LP fits. BPL fits yielded $\Gamma_1 \approx 2.72$ , $\Gamma_2 \approx 2.55$ , and a break energy $E_{break} \approx 2.87$ keV. The spectral hardening $\Delta\Gamma \approx 0.17$ .
- Obs B (0790381501, Feb 13, 2017, ~8 hrs later): Also showed negative curvature. BPL fit yielded $\Gamma_1 \approx 2.57$ , $\Gamma_2 \approx 2.33$ , $E_{break} \approx 4.26$ keV, and $\Delta\Gamma \approx 0.24$ .
- Significance: Both observations showed stable BPL fits with independent parameters at 99% confidence. The presence of a high-energy tail in PL fits that was rectified by BPL fits, combined with the low flux state of Obs A, strongly suggests an intermixing of the synchrotron and IC components.
Weak Candidate (Obs C): Obs ID 0790381001 showed mild hardening ( $\Delta\Gamma \approx 0.04$ ) but lacked a distinct high-energy tail in PL fits, making the IC signature less certain.

4. Significance and Conclusions

Comprehensive Baseline: This study provides the most extensive spectral analysis of PG 1553+113 using XMM-Newton data to date, covering 30 observations over two decades.
Acceleration Mechanisms: The prevalence of Log-Parabolic spectra supports the stochastic (second-order Fermi) acceleration model, where the probability of particle acceleration decreases with energy. The lack of $\alpha$ - $\beta$ correlation suggests complex interplay between acceleration and cooling timescales.
Transient IC Detection: The detection of spectral hardening in two specific observations (Feb 2017) provides rare evidence of the Inverse Compton component becoming visible in the soft X-ray band of an HBL. This suggests that PG 1553+113 can temporarily exhibit spectral characteristics similar to Intermediate-energy peaked Blazars (IBLs) when the synchrotron and IC components overlap.
Synchrotron Peak Location: The joint analysis confirms the synchrotron peak lies in the UV range (4.6–48.6 eV), consistent with its classification as an HBL, but highlights the difficulty in constraining this peak without broad-band coverage extending into the optical/UV.

Future Outlook: The authors note that five additional XMM-Newton observations are expected to be released in 2026, which will further improve temporal coverage and allow for a more detailed investigation of the long-term spectral evolution and the QPO mechanism in this source.