🔭 astrophysics

25 Years of Mrk 421 with XMM-Newton: Unveiling Structured Jets and Energy-Dependent Escape

Based on a quarter-century of XMM-Newton observations, this study reveals that the TeV blazar Mrk 421 exhibits complex, non-Gaussian variability and spectral transitions between broken power-law and log-parabola models, providing strong evidence for structured jets or energy-dependent electron escape mechanisms rather than standard one-zone cooling scenarios.

Original authors: Tek P. Adhikari (CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei, China, School of Astronomy and Space Science

Published 2026-10-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Tek P. Adhikari (CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei, China, School of Astronomy and Space Science, University of Science and Technology of China, Hefei, China), Gopal Bhatta (Janusz Gil Institute of Astronomy, University of Zielona Góra, Zielona Góra, Poland), Sangeetha Kizhakkekalam (Janusz Gil Institute of Astronomy, University of Zielona Góra, Zielona Góra, Poland, Doctoral School of Exact and Technical Sciences, University of Zielona Góra, Zielona Góra, Poland), Navaneeth P K (Janusz Gil Institute of Astronomy, University of Zielona Góra, Zielona Góra, Poland, Doctoral School of Exact and Technical Sciences, University of Zielona Góra, Zielona Góra, Poland), Alex Markowitz (Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Warsaw, Poland), Zhicheng He (CAS Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei, China, School of Astronomy and Space Science, University of Science and Technology of China, Hefei, China), Santanu Mondal (Department of Physics and Electronics, Christ University, Bangalore, India)

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

Deep in the vast, dark ocean of space, there exist galaxies that act as cosmic lighthouses, beaming intense radiation across the entire electromagnetic spectrum. Among these, a special class known as blazars stands out because their powerful jets of particles are aimed almost directly at Earth. When we look at these objects, we are seeing the result of matter being accelerated to nearly the speed of light, creating a brilliant, shifting display of energy that changes from minute to minute and year to year. One of the closest and brightest of these cosmic beacons is a galaxy called Markarian 421. For decades, astronomers have watched it, trying to understand the invisible machinery inside its jet that heats up particles to such extreme temperatures. The central mystery has been how these particles gain their energy and why the light they emit changes color so dramatically as the galaxy flares up and calms down.

To solve this puzzle, a team of researchers turned to a massive archive of data collected over a quarter of a century. They focused on the X-ray light from Markarian 421, which is produced by the most energetic electrons in the jet. Using a space telescope called XMM-Newton, they gathered observations from 76 different moments in time, spanning from the year 2000 to 2025. This long timeline allowed them to see the full story of the galaxy's behavior, rather than just a single snapshot. They measured how bright the galaxy was at different times and analyzed the specific "color" or energy of the X-rays it emitted. By comparing the soft, lower-energy X-rays with the hard, high-energy ones, they could trace how the particles inside the jet were moving, accelerating, and losing energy over decades.

The researchers found that Markarian 421 is a creature of extremes. Its brightness fluctuates wildly, sometimes becoming ten times brighter than its quietest states. When the galaxy flares, it does not just get brighter; it also gets "harder," meaning the light shifts toward higher, more energetic X-rays. This pattern, where the source becomes more energetic as it brightens, suggests that the same process is responsible for both the increase in light and the shift in energy. However, the relationship is not perfectly simple. The team noticed that the galaxy does not always follow a single, predictable path. Sometimes, at the same level of brightness, the X-rays can have different colors, implying that the physical conditions inside the jet are changing in complex ways that depend on the history of the flare, not just its current intensity.

When the scientists looked closely at the shape of the X-ray spectrum, they discovered that a simple, straight-line description of the energy distribution was not enough to explain what they saw. In most cases, the data required a more complex model that showed a distinct change in slope, or a smooth curve, rather than a uniform line. This indicates that the electrons producing the light are not behaving as a single, uniform group. Instead, the data suggests a scenario where the jet has a structured environment. The particles seem to be accelerated at a specific front, perhaps a shockwave, and then move into a turbulent region where they lose energy. The fact that the point where the energy distribution changes remains clustered around a specific value of 2 keV across all these years suggests that the physical size and conditions of this emission region are remarkably stable, even as the brightness varies.

Perhaps the most significant finding challenges a long-held assumption about how these jets work. Standard theories predicted that as high-energy electrons radiate their energy away, the spectrum should steepen by a specific, large amount. However, the measurements from this 25-year study showed a much gentler slope change. This small difference is a crucial clue. It suggests that the electrons are not just cooling down in a simple, uniform zone. Instead, the researchers propose that the particles are escaping from the emission region at different speeds depending on their energy. High-energy particles seem to leave the zone faster than low-energy ones, which softens the expected drop in the spectrum. This energy-dependent escape, combined with the acceleration and cooling processes, paints a picture of a dynamic, multi-layered jet where particles are constantly being injected, shuffled, and lost.

By piecing together two and a half decades of observations, the study provides a detailed map of how a relativistic jet evolves over time. It confirms that the jet is a complex system where acceleration and cooling compete, and where the escape of particles plays a vital role in shaping the light we see. The results do not just describe a single event but reveal the underlying rules that govern the behavior of these powerful cosmic engines. The work suggests that to truly understand blazars, we must move beyond simple models and consider how particles move through structured, changing environments. This deeper understanding brings us closer to unraveling the physics of the most energetic processes in the universe.

Technical Summary: 25 Years of Mrk 421 with XMM-Newton: Unveiling Structured Jets and Energy-Dependent Escape

Problem Statement
The high-synchrotron-peaked (HSP) blazar Mrk 421 serves as a critical laboratory for understanding particle acceleration and radiation mechanisms in relativistic jets. While previous studies have characterized its variability on timescales ranging from hours to years, a comprehensive, long-term spectral and timing analysis spanning a quarter-century has been lacking. Specifically, the physical mechanisms governing the source's spectral evolution—whether driven by simple one-zone synchrotron cooling, structured jet scenarios, or energy-dependent particle escape—remain debated. Standard models often predict a spectral steepening ( $\Delta\Gamma$ ) of 0.5 above a cooling break, yet observations frequently show deviations. This work addresses the need to resolve long-term jet evolution, constrain acceleration/cooling models, and link spectral parameters to flux states using the most extensive XMM-Newton dataset to date.

Methodology
The authors analyzed nearly 25 years (2000–2025) of archival XMM-Newton/pn observations, comprising 76 distinct epochs in both IMAGING and TIMING modes.

Data Reduction: Standard SAS (v21.0) procedures were employed, including background flare filtering (excluding intervals >0.4 counts/s), pile-up correction via iterative core excision, and extraction of light curves in soft (0.3–2 keV) and hard (2–10 keV) bands with 100 s time bins.
Timing Analysis: The study calculated fractional variability ( $F_{var}$ ) to quantify intrinsic variability amplitudes and examined correlations between soft and hard bands. Flux distributions were fitted with lognormal functions to characterize the source's brightness states.
Spectral Analysis: Spectra were fitted using XSPEC (v12.13.1) over the 0.6–10 keV range, corrected for Galactic absorption ( $N_H = 1.92 \times 10^{20}$ cm $^{-2}$ ). Three models were tested: Power-Law (PL), Broken Power-Law (BPL), and Log-Parabola (LP). Model selection was rigorously determined using the F-test, corrected Akaike Information Criterion (AICc), and Bayesian Information Criterion (BIC).
Parameter Correlations: The authors investigated the "harder-when-brighter" trend via Hardness Ratios (HR) and analyzed correlations between spectral indices ( $\Gamma_1, \Gamma_2$ ), break energy ( $E_b$ ), and flux states.

Key Results

Variability and Flux Distribution: Mrk 421 exhibits flux variations of approximately one order of magnitude, ranging from $\sim$ 100 to $>900$ counts/s. The flux distribution is non-Gaussian and lognormal-like, dominated by low-to-intermediate states with a tail of rare, intense flares.
Harder-When-Brighter Trend: A significant "harder-when-brighter" trend was confirmed. The fractional variability in the hard band (2–10 keV) correlates tightly with the soft band (0.3–2 keV) but with a steeper slope (1.41), indicating that high-energy electrons undergo more rapid changes in their energy distribution than low-energy electrons.
Spectral Modeling: The simple Power-Law model was insufficient for any epoch. The best-fit models were predominantly BPL (59 epochs) and LP (17 epochs). This fluctuation between models suggests the emission region transitions between a distinct acceleration shockfront and downstream stochastic turbulence.
Spectral Breaks and Steepening: The spectral break energy ( $E_b$ ) for the BPL model remains clustered around 2 keV across all flux states. Crucially, the spectral steepening above the break ( $\Delta\Gamma = \Gamma_2 - \Gamma_1$ ) is mild, averaging $\approx 0.19$ , which is significantly smaller than the canonical 0.5 predicted by standard one-zone synchrotron cooling models.
State-Dependent Evolution: The correlation between spectral indices and break energy reverses polarity depending on the flux state. Low-flux states show negative correlations between indices and $E_b$ , while high-flux states show weak positive correlations, implying distinct physical regimes govern quiescent versus flaring phases.
Hysteresis: Intra-observation HR analysis revealed complex hysteresis loops (both clockwise and counter-clockwise), indicating time-dependent interplays between particle acceleration, cooling, and escape timescales.

Significance and Conclusions
The paper concludes that the standard one-zone cooling approximation is insufficient to explain the long-term X-ray behavior of Mrk 421. The observed mild spectral steepening ( $\Delta\Gamma \approx 0.19$ ) supports scenarios involving structured jets or energy-dependent particle escape. By modeling the electron continuity equation with an energy-dependent escape timescale ( $t_{esc} \propto \gamma^{-k}$ ), the authors derive a diffusion index of $k \approx 0.62$ , which naturally produces the observed weak steepening.

The study provides the longest temporal baseline (25 years) for Mrk 421, demonstrating that while the leptonic framework broadly applies, the jet emission is governed by a complex interplay of particle acceleration, radiative cooling, and energy-dependent transport. The stability of the break energy around 2 keV suggests a relatively stable physical scale in the emitting region, with variability driven primarily by changes in particle injection rates rather than major structural changes in the jet. These findings offer key constraints for future models of relativistic jets in high-synchrotron-peaked blazars.

More like this