🔭 astrophysics

Not so Swift: 20 years of multiwavelength observations of Mrk 421 and Mrk 501

This study analyzes 20 years of multiwavelength Swift observations of blazars Mrk 421 and Mrk 501 to characterize their long-term variability, revealing energy-dependent flux distributions, a general harder-when-brighter trend in X-rays, and a potential 390-day quasi-periodicity in Mrk 501's X-ray emission, while finding no correlation between optical/UV and X-ray bands.

Original authors: Gabrielle L. Taylor, Stefan J. Wagner, Alicja Wierzcholska, Michael Zacharias

Published 2026-02-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Gabrielle L. Taylor, Stefan J. Wagner, Alicja Wierzcholska, Michael Zacharias

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 two cosmic lighthouses, Mrk 421 and Mrk 501, shining billions of light-years away. These aren't ordinary lighthouses; they are blazars, which are supermassive black holes at the centers of galaxies, shooting powerful beams of energy (jets) directly at Earth. Because we are looking straight down the barrel of these jets, they appear incredibly bright and change their brightness rapidly.

For 20 years (from 2005 to 2025), astronomers have been watching these two "lighthouses" with a space telescope called Swift. They watched them glow in UV light (like a blacklight) and X-rays (the kind that see through your skin). The goal was to figure out how these black holes are behaving and to test our theories about how they work.

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

1. The "Harder-When-Brighter" Rule

Think of the light from these blazars like a car engine.

When the engine is idling (low brightness), it runs on a "soft" fuel mix (lower energy light).
When you step on the gas and the engine roars (high brightness), it switches to a "hard" fuel mix (high energy light).

The astronomers found that both Mrk 421 and Mrk 501 follow this rule perfectly. When they get brighter, their light becomes "harder" (more energetic). This suggests that the particles inside the jet are being accelerated very quickly, much faster than they can cool down.

2. The "Two-Engine" Mystery (Multi-Zone Model)

For a long time, scientists thought these blazars were powered by a single, giant engine (a "one-zone" model). If this were true, the UV light and the X-ray light should change brightness at the exact same time, like two headlights on the same car turning on and off together.

The Surprise: The data showed that the UV light and X-ray light do not dance together. Sometimes the X-rays flare up while the UV stays calm, and vice versa.

The Analogy: Imagine a house with two separate rooms. In one room, a party is raging (X-rays flaring), while in the other room, everyone is sleeping (UV staying calm).
The Conclusion: This proves that these blazars aren't powered by one simple engine. They are likely powered by multiple zones or engines working independently within the jet. Sometimes they sync up, but often they are doing their own thing.

3. The "Lognormal" Lottery

When the astronomers looked at the history of the blazars' brightness, they noticed a specific pattern. The brightness didn't just go up and down randomly; it followed a specific mathematical curve called a lognormal distribution.

The Analogy: Think of a snowball rolling down a hill. As it rolls, it picks up more snow. The bigger it gets, the more snow it picks up in the next second. The size of the snowball at the bottom isn't random; it's the result of a chain reaction where small changes get multiplied into huge changes.
What it means: This suggests that the "fuel" for these blazars comes from the black hole's accretion disk (the swirling matter around it). Small, random fluctuations in the disk get multiplied as they travel out into the jet, creating the huge flares we see.

4. The "Ghost Rhythm" (Quasi-Periodicity)

Scientists love finding rhythms in nature, like a heartbeat. They wondered: Do these blazars pulse with a regular beat?

Mrk 421: No rhythm found. It's like a drummer playing completely randomly.
Mrk 501: There might be a heartbeat! In the X-ray data, they found a potential rhythm of about 390 days (roughly a year).
The Analogy: Imagine a lighthouse that spins. If the light beam hits a specific spot on a rotating platform, you might see a flash every time the platform completes a turn.
The Cause: This rhythm could be caused by a "hot spot" of gas orbiting the black hole, or perhaps the black hole itself is wobbling (precessing) like a spinning top. It's a clue that something physical and massive is moving in a circle near the black hole.

5. The "Silent Partner" (No Correlation)

One of the biggest questions was: Does the UV light tell us anything about the X-ray light?

The Result: No. Even though they come from the same object, they don't seem to talk to each other on long timescales.
The Analogy: It's like having two neighbors. One is a drummer (X-rays) and the other is a singer (UV). Sometimes they practice at the same time, but usually, they are on completely different schedules. You can't predict when the singer will start just by listening to the drummer.

The Big Picture

This 20-year study is like watching a movie in slow motion. By looking at such a long time span, the astronomers realized that these cosmic monsters are complex, multi-layered systems.

They aren't simple, single-beam machines. They are chaotic, multi-zone environments where different parts of the jet act independently, fueled by a "snowball effect" from the black hole's disk. While Mrk 501 might have a mysterious yearly rhythm, Mrk 421 remains a chaotic free spirit.

In short: The universe is messier and more complex than our simple models suggested. These black holes are not just one big engine; they are a whole orchestra of engines playing different tunes, sometimes in sync, but often not.

1. Problem Statement

Blazars, specifically High-frequency peaked BL Lac objects (HBLs) like Mrk 421 and Mrk 501, exhibit complex variability across the electromagnetic spectrum. While short-term variability (days to weeks) and gamma-ray behavior are well-studied, the long-term temporal behavior (decades) at synchrotron energies (optical, UV, X-ray) remains less understood.

A central challenge in blazar physics is distinguishing between one-zone and multi-zone emission models.

One-zone models (e.g., Synchrotron Self-Compton) predict strong correlations between all bands within the synchrotron peak (optical/UV and X-ray).
Multi-zone models suggest that different emission regions or mechanisms operate independently, potentially explaining the frequent lack of correlation observed between optical/UV and X-ray bands.

Previous studies often focused on short timescales (<1 week) or lacked simultaneous coverage across sufficient wavelengths. This paper aims to characterize the long-term (20-year) temporal behavior of Mrk 421 and Mrk 501 to assess current emission models and investigate the presence of Quasi-Periodic Oscillations (QPOs).

2. Methodology

The authors analyzed a unique dataset comprising 20 years of quasi-simultaneous observations (2005–2025) from the Neil Gehrels Swift Observatory.

Data Sources:
- UV/Optical: Swift-UVOT telescope using six filters (V, B, U, UVW1, UVM2, UVW2). Narrow-band UV filters were combined into a single weighted "UV" light curve to maximize data density.
- X-ray: Swift-XRT telescope (0.3–10 keV).
Data Processing:
- Fluxes were corrected for Galactic dust extinction.
- X-ray spectra were fitted using two models: a Single Power Law (SPL) and a Log-Parabola (LP). The F-test was used to determine the best fit.
- Spectral indices ( $\Gamma$ or $\alpha$ ) were derived to classify spectral states as "hard" ( $\Gamma < 2$ ) or "soft" ( $\Gamma > 2$ ).
Analytical Techniques:
- Fractional Variability ( $F_{var}$ ): Calculated to quantify variability amplitude across bands.
- Flux Distributions: Tested against lognormal distributions to identify underlying stochastic processes (additive vs. multiplicative).
- Time Series Analysis: Employed three methods to search for periodicity:
  1. Structure Functions (SF): To identify characteristic timescales and break points.
  2. Lomb-Scargle Periodograms (LSP): Generalized to handle uneven sampling and measurement errors.
  3. Discrete Correlation/Autocorrelation Functions (DCF/DACF): To detect time lags and periodicities, validated against synthetic light curves (Monte Carlo simulations) to establish significance ( $\ge 3\sigma$ ).
- Cross-Correlation: Discrete Cross-Correlation Functions (DCCF) were used to test for time lags between UV and X-ray bands.

3. Key Contributions

Dataset Scale: This study utilizes the longest continuous light curves (20 years) for these sources at synchrotron energies, allowing for the probing of timescales previously inaccessible.
Spectral Modeling: Confirmed that curved Log-Parabola models are statistically superior to Power Laws for describing the X-ray spectra of both sources over the long term.
Correlation Analysis: Systematically tested correlations not just between bands, but between specific spectral states (hard vs. soft X-rays) and UV flux.
QPO Search: Conducted a rigorous search for QPOs using multiple statistical methods, distinguishing between physical signals and windowing artifacts.

4. Key Results

A. Variability and Spectral Behavior

Energy Dependence: Both sources exhibit energy-dependent variability. X-ray flux shows significantly higher fractional variability ( $F_{var} \approx 0.6-0.7$ ) compared to UV/Optical ( $F_{var} \approx 0.1-0.4$ ). This supports the hypothesis that X-rays are produced by higher-energy electrons with shorter cooling times.
Harder-When-Brighter: Both sources display a clear "harder-when-brighter" trend. As X-ray flux increases, the photon index hardens (spectral slope decreases). This suggests rapid electron acceleration relative to cooling rates.
Flux Distributions: Flux distributions follow a lognormal shape, but the degree of skewness (lognormality) is energy-dependent. X-ray distributions are more strongly lognormal than UV distributions. This energy dependence suggests that photons originate from different emission zones, supporting a multi-zone model.

B. Correlations

Optical/UV vs. X-ray: No significant simultaneous correlation was found between the optical/UV and X-ray bands over the full 20-year dataset.
Spectral State Dependence: When X-ray data were split into "hard" and "soft" subpopulations, a weak increase in correlation with UV was observed, but it remained low.
- Mrk 421: UV correlated slightly better with hard X-rays.
- Mrk 501: UV correlated slightly better with soft X-rays.
Conclusion on Correlations: The lack of consistent correlation contradicts the one-zone SSC model and strongly favors a multi-zone emission model, where different zones may flare independently or with varying time delays.

C. Quasi-Periodic Oscillations (QPOs)

Mrk 421: No significant QPOs were detected in UV or X-ray bands.
Mrk 501:
- X-ray: A potential QPO was detected at a period of $\sim 400$ days (observed frame) / $\sim 390$ days (host frame) with $>3\sigma$ significance in the X-ray LSP.
- UV: No QPOs detected in the UV band.
- DACF Consistency: The DACF showed a trough at $\sim 1400$ days (consistent with $3.5 \times 400$ days) but lacked a clear peak at the fundamental period, suggesting the signal is weak or the sampling is insufficient to confirm it definitively.
Interpretation: If real, the 390-day period in Mrk 501 could be attributed to accretion disc instabilities (e.g., hot spots, Lense-Thirring precession) imprinting on the jet, rather than a supermassive binary black hole system (which is constrained by gravitational limits for such periods).

5. Significance and Conclusions

This paper provides a comprehensive 20-year baseline analysis of two of the nearest HBLs, offering critical constraints on blazar physics:

Model Validation: The results strongly disfavor a simple one-zone model. The lack of consistent optical/UV-X-ray correlation, combined with energy-dependent lognormality and variability, necessitates a multi-zone emission model.
Physical Mechanisms: The "harder-when-brighter" trend and lognormal flux distributions point to stochastic multiplicative processes, likely driven by fluctuations in the accretion disc that propagate to the jet.
QPO Implications: The tentative detection of a $\sim 390$ -day QPO in Mrk 501 X-rays offers a potential window into the dynamics of the accretion disc-jet connection, though further monitoring is required to confirm the periodicity.
Future Outlook: The study highlights the limitations of current sampling rates in detecting short-term QPOs and the need for denser, continuous monitoring to resolve the complex interplay between different emission zones in blazars.

In summary, the 20-year dataset reveals that Mrk 421 and Mrk 501 are complex, multi-component systems where long-term variability is governed by stochastic processes likely rooted in the accretion disc, rather than a single, coherent emission zone.