Pith. sign in

REVIEW 4 major objections 5 minor 7 references

Spectral State Switching in Mrk 421: Results from the AstroSat LAXPC/SXT Observations

T0 review · 4 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The blazar Mrk 421 alternates between two discrete X-ray spectral states, according to flux- and time-resolved spectroscopy with AstroSat data.

desk verdict The two-state spectral claim doesn't survive contact with the paper's own time-resolved table — Γ1 drifts continuously with flux, so the flux-resolved 'clusters' are likely binning artifacts. read the letter →

arxiv 2512.21885 v2 pith:L24M6DTK submitted 2025-12-26 astro-ph.HE

Sikandar Akbar , Zahir Shah , Ranjeev Misra , Naseer Iqbal This is my paper
classification astro-ph.HE
keywords blazarMrk421X-rayspectroscopyspectralstatesbrokenpower-lawdoublelognormalflux-resolvedparticleacceleration
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper aims to show that the archetypal high-synchrotron-peaked blazar Mrk 421 routinely occupies two dominant X-ray spectral states. Using simultaneous soft and hard X-ray observations from AstroSat, the authors find that the low-energy particle index clusters around two discrete values (~2.5 and ~2.1) across ten flux states, and that time-resolved spectroscopy yields a double lognormal flux distribution and a bimodal index distribution. The same behaviour appears in two additional long observations spanning 2017–2019. If correct, this indicates that the jet's particle acceleration switches between two regimes, connecting observable spectral states to distinct physical conditions in the jet.

What carries the argument

The synchrotron-convolved broken power-law model (synconv ⊗ n(ξ)) used to fit joint SXT and LAXPC20 spectra, which provides a better fit than a log-parabola in all flux states. The analysis divides the 100-second binned light curve into ten flux states of 8 counts s⁻¹ width for flux-resolved spectroscopy, and uses 10-ks segments for time-resolved spectroscopy; the distributions of flux and index are then fit with two-component (double lognormal and double normal) functions to characterise bimodality.

What would settle it

Extract flux-resolved spectra using narrower flux bins (e.g., 2 counts s⁻¹) over the same dataset and check whether Γ1 varies smoothly with mid-bin flux; if it does, the apparent two-state pattern disappears. Additionally, test whether a single skewed distribution (e.g., a gamma or lognormal) fits the time-resolved index histogram as well as the double normal; if yes, the bimodality is not required.

Watch

Extended reading notes

Core claim

The central claim is that flux-resolved spectroscopy of Mrk 421, using a synchrotron-convolved broken power-law model, reveals that the low-energy particle index is not a smooth function of flux but clusters around two values across flux states—indicating two spectral states in the source. This is corroborated by time-resolved spectroscopy: the flux histogram is best described by a double lognormal distribution (χ²ν = 1.05 for one observation, 0.91 for three combined) and the index histogram by a double normal distribution, and the behaviour persists across observations from 2017–2019. The authors conclude that Mrk 421 routinely occupies two dominant spectral states and interpret this within

Load-bearing premise

The conclusion rests on reading the two clusters of low-energy index values as genuinely discrete spectral states; this could collapse if the clustering is an artifact of the chosen flux binning and the index actually hardens continuously with flux.

Editorial extensions

If this is right

  • If the two states are real, jet models for Mrk 421 must account for discrete switches in particle acceleration efficiency rather than a continuum.
  • The double lognormal flux distribution supports multiplicative variability, i.e., Gaussian fluctuations in acceleration timescales, as proposed in earlier theoretical work.
  • The persistent two-state behaviour over multiple epochs suggests it is an intrinsic property of the source, not a methodological artifact.
  • The increasing break energy with flux (harder-when-brighter) links the two-state switching to shifts of the synchrotron peak to higher energies.
  • The low-energy index states being stable while the high-energy index varies implies that the variability is driven by changes in the high-energy electron population.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The two-cluster pattern in flux-resolved data could partly reflect the coarse flux binning (8 counts/s, with the lowest two states merged); a monotonic hardening of Γ1 with flux would appear step-like if bins are wide. A finer binning test would separate a continuous trend from true two-state switching.
  • The time-resolved index values in Table 4 appear to vary gradually (from ~2.9 to ~2.2) within the observation; a double-normal fit might be capturing the ends of a continuous distribution rather than two distinct states, and a single skewed distribution could be tested.
  • If the two states are genuine, similar bimodality might be expected in optical/UV and GeV/TeV light curves of the same source, offering a multi-wavelength falsification, and possibly in other HBL blazars with long monitoring.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper presents a time- and flux-resolved X-ray spectral analysis of Mrk 421 using simultaneous AstroSat SXT and LAXPC20 observations from 3–8 January 2017, plus two additional long observations spanning 2017–2019. The authors divide the 100-s binned light curve into ten count-rate flux states (S1–S10) and fit a synchrotron-convolved broken power law to each state, reporting that the low-energy index Γ1 clusters around two discrete values (≈2.5 for S1–S6 and ≈2.1 for S7–S10). They also perform time-resolved spectroscopy in 10-ks segments and fit double-lognormal mixtures to the flux distribution and double-normal mixtures to the index distribution. The paper concludes that Mrk 421 'routinely occupies two dominant spectral states', implying discrete switching between two particle-acceleration regimes.

Significance. If the two-state claim were firmly established, it would be an interesting result for jet-acceleration physics, adding to earlier reports of bimodal distributions in blazar variability. The paper also applies standard reduction pipelines and gives new flux-resolved results with a physically motivated emission model. However, the central claim is not supported by the evidence presented: the paper's own time-resolved table shows a continuous evolution of Γ1 with flux, and no statistical test is performed to distinguish discrete states from a continuous Γ1–flux trend. The apparent bimodality in the index histogram can be a consequence of the double-lognormal flux distribution under a monotonic index–flux relation, making the two 'independent' routes redundant. The significance of the work therefore hinges on analysis that is currently missing.

major comments (4)
  1. [§2.3, Table 4] Table 4 directly contradicts the two-state interpretation. The time-resolved Γ1 values drift continuously from 2.90 at 5000 s to 2.20 at 335000 s, with intermediate values 2.68, 2.57, 2.49, 2.39, etc., and a flux-dip segment at 355000 s returning to 2.78. The flux-resolved Table 1 averages over count-rate bins; the S6–S7 jump is simply a bin boundary in this smooth trend. The paper never tests a continuous Γ1–flux relation against a step model, nor does it examine the robustness of the apparent two clusters to the bin width (8 counts/s), bin edges, or the merging of the top two states. Given the significant Γ1–normalization correlation reported in Fig. 4 (r_s = −0.70, p = 0.025), the two-cluster appearance is not evidence for discrete states without such a test.
  2. [§2.3–§2.4] The claimed independence of the flux-resolved and time-resolved routes is not valid. Both are derived from the identical set of joint SXT+LAXPC20 spectral fits and from the same underlying count-rate data; the time-resolved parameters are simply taken from the authors' previous work (Akbar et al. 2025), and the flux-resolved analysis is a binning of the same light curve. The double-lognormal flux distribution, when mapped through a monotonic Γ1(log F) relation, will produce a bimodal-looking index histogram even if the underlying relation is continuous. The paper does not test whether the index distribution is bimodal after removing a smooth flux–index correlation, so the 'agreement' between the two methods is not independent confirmation.
  3. [§2.4.1, Table 2] The statistical evidence for bimodality is weak. For the single long observation, the double-normal fit to the index histogram has χ²ν = 2.04, which is a poor fit (χ²ν ≈ 2 with the number of bins used); the Anderson–Darling statistic against a single Gaussian is 0.74, only marginally above the 5% critical value of 0.726. The paper does not report the single-component χ²ν for comparison, nor does it carry out a formal significance test (e.g., F-test or likelihood-ratio test) of the improvement from adding a second component. The same criticism applies to the combined dataset in Table 3. A poor fit by a two-component model does not by itself establish bimodality.
  4. [§2.2, Fig. 2] The claim that the broken-power-law model is 'consistently a better fit' than the log-parabola is based on a visual comparison of reduced χ² values only, without any statistical model-comparison test (e.g., F-test, AIC, or BIC). The reduced χ² values in Fig. 2 appear close to unity for both models in all states. This point is secondary to the main conclusion, but it underpins the choice of the BPL parameters and should be quantified properly.
minor comments (5)
  1. [Abstract] The abstract contains an incomplete phrase: 'two dominant spectral;' should read 'two dominant spectral states'. Also, the abstract gives F_rms values for both SXT and LAXPC20 but the text in §2.1 says 'the obtained F_var value shows that the source had significant varibility' with typos ('varibility').
  2. [Fig. 6] The left panel of Fig. 6 is labeled 'Index' but the x-axis values (≈0.34–0.48) appear to be log10(Γ1) rather than Γ1 itself; the caption and Table 2 are inconsistent about whether the double-lognormal fit for the index is applied to the linear index or its logarithm. This should be clarified.
  3. [§1.2] The text says 'we restricted the spectral analysis to the 3–20 keV energy range' but the light curve is generated in the 3–30 keV band; the abstract states the LAXPC energy range as 3–80 keV. Please make the energy ranges consistent throughout.
  4. [Eq. (7)] Equation (7) contains a factor 10^x ln 10 that is appropriate for converting a normal distribution in linear space to a variable x = log10(F), but this factor is not appropriate when x is a photon index. The definition of the 'double normal' model used for the index distribution needs to be stated without this factor or otherwise clarified.
  5. [§2.2, Table 1] In Table 1, the normalization uncertainties are quoted as '+0.00/−0.00' for S1–S10, which appears to be a typographical artifact; the normalizations must have non-zero uncertainties. Please check the quoted errors.

Circularity Check

2 steps flagged · score 5.0 of 10

Central 'two-state' claim reduces to a binned continuous hardening trend; the 'independent' time-resolved confirmation is inherited from the same authors' previous broken power-law fits.

  1. renaming known result [Section 2.2 / Fig. 3-4 and Table 1 (flux-resolved) vs. Table 4 and Fig. 4; see also Section 3]
    "The variation of the broken power-law parameters across flux states (Fig. 3) reveals that the low-energy index (Γ1) is distributed over two distinct levels, suggesting the presence of two electron populations or emission zones contributing to the observed X-ray spectrum."

    The 'states' S1–S10 are defined in §2.1 as contiguous 8 counts/s intervals of the LAXPC count rate ('each having a width of 8 counts s−1'), so the state variable is just a binned coordinate of flux. The same paper's Table 4 shows Γ1 declining continuously from 2.90 to 2.20 as logF rises from −9.13 to −8.74, with intermediate values (2.68, 2.57, 2.49, 2.39, …) throughout; Fig. 4 independently reports a significant Γ1–normalization correlation (rs = −0.70, p = 0.025). The S1–S6/S7–S10 'clusters' are therefore the binned endpoints of the known harder-when-brighter trend sliced at one bin boundary. No continuous Γ1(logF) relation is fitted or rejected, so the discrete-state conclusion is a relabeling of the binned trend rather than a new, independent result.

  2. self citation load bearing [Section 2.3 and Section 3 (Discussion)]
    "Here, we reproduce only the spectral parameters derived from the broken power-law model for each 10 ks segment, which are directly taken from that analysis and used in the present study for comparison with flux-resolved results. ... The agreement between these approaches, which probe variability in fundamentally different ways, indicates that the observed state-switching is not an artifact of a particular analysis method or timescale selection."

    The claimed independent confirmation of two-state behaviour is not independent: the time-resolved spectral parameters are taken from the same authors' prior paper (Akbar et al. 2025), and the time-resolved fluxes are computed with the same spectral model ('using the same spectral model (synchrotron convolved broken power law model)'). Both the flux-resolved and time-resolved routes therefore share the same observation, the same BPL model, and the same fitting procedure. The 'independence' is a self-consistency check within the authors' own analyses, not an external or methodologically distinct verification; the known flux–index correlation further ensures the two histograms are linked rather than constituting two separate lines of evidence.

full rationale

The paper's central claim — that Mrk 421 routinely occupies two discrete spectral states — rests on (a) an apparent clustering of Γ1 into two levels in flux-resolved spectroscopy and (b) a bimodal index distribution plus double-lognormal flux distribution in time-resolved data. Neither strand is derived independently of the binning or of the authors' prior fits. Flux states are defined as fixed-width count-rate bins, so parameter values per state are conditional averages over flux intervals; the paper's own Table 4 shows a continuous, roughly monotonic Γ1 decline with rising flux, and Fig. 4 reports a significant Γ1–normalization correlation. The step-like S1–S6/S7–S10 pattern is thus a binned representation of a known continuous harder-when-brighter relation, not a newly derived discrete-state result. The time-resolved 'confirmation' is not independent: the parameters are inherited from Akbar et al. (2025), the fluxes are recomputed with the same synchrotron-convolved broken power-law model, and the flux and index histograms are linked by the same flux–index correlation. These issues do not make the spectral fits themselves circular, and much of the paper's data reduction and model comparison is self-contained; but the specific two-state interpretation is substantially over-determined by the chosen binning and by self-referential reuse of prior fits. Hence a partial circularity score of 5 is appropriate.

Assumptions & free parameters 3 free parameters · 5 assumptions · 1 invented entities

The central claim rests on standard single-zone synchrotron modeling, a set of ad hoc flux-state binning choices, mixture-model parameters fitted to the X-ray histograms, and a large inherited block of time-resolved fits from the authors' prior paper. No new entities are introduced with independent evidence; the 'two states' construct is an interpretation of the same data.

free parameters (3)
  • Flux-state binning choices (width/threshold/merging) = 8 counts/s; <20 counts/s excluded; top two states merged
    Ad hoc in Section 2.1; the apparent Γ1 clustering is shaped by these choices and no robustness check is reported.
  • Histogram mixture parameters (a, µ, σ) = e.g., flux double-lognormal: a=0.47±0.11, µ1=−9.06±0.03, µ2=−8.81±0.01 (Table 2)
    Fitted parameters used to claim double-lognormal flux / bimodal index; uncertainties in the index second component are large (a=0.17±0.10).
  • Relative normalization and gain offset between SXT and LAXPC20 = free in fits (values not tabulated)
    Nuisance calibration parameters that affect the joint spectral fits and hence the quoted Γ1 values and errors.
assumptions (5)
  • domain assumption X-ray emission is single-zone leptonic synchrotron radiation from a homogeneous spherical region
    Section 2.2: the synconv model assumes a homogeneous sphere with magnetic field B and isotropic electron distribution n(γ); standard for blazar X-ray spectral fits but not independently tested here.
  • domain assumption The synchrotron convolution model (synconv, Eq. 4) correctly implements the pitch-angle-averaged emissivity
    The results depend entirely on the convolved broken power-law fit in XSPEC; no independent validation of the model implementation is given.
  • domain assumption A lognormal flux distribution implies multiplicative variability, and a double lognormal implies two states
    Section 2.4 and Discussion invoke Sinha et al. (2018) and Khatoon et al. (2020); the interpretation that Gaussian index fluctuations produce lognormal flux is adopted, not derived here.
  • domain assumption The time-resolved spectral parameters in Table 4 (from Akbar et al. 2025) are correct and model-consistent
    Section 2.3 takes parameters 'directly from that analysis'; the combined-distribution claim inherits those fits without re-derivation.
  • standard math Two-component Gaussian/lognormal mixture models (Eqs. 5-7) adequately represent the histograms
    Mixture fitting is standard, but parameter identifiability is marginal given the small number of segments and large component uncertainties.
invented entities (1)
  • Two discrete particle-index states / two electron populations or emission zones
    purpose: Explains the bimodal Γ1 and flux distributions
    In Section 3 the paper suggests two electron populations/emission zones, but presents no independent observational handle (e.g., a spectral signature in another band); the entity is an interpretation of the same X-ray data, not a separately constrained physical object.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spectral State Switching in Mrk 421: Results from the AstroSat LAXPC/SXT Observations." pith.science (2026). https://pith.science/paper/L24M6DTK

@misc{pith2026251221885,
  author       = {Pith},
  title        = {Pith review of: Spectral State Switching in Mrk 421: Results from the AstroSat LAXPC/SXT Observations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/L24M6DTK}},
  note         = {Machine review of arXiv:2512.21885}
}
abstract

We carried a detailed time and flux resolved X-ray spectral analysis of the high-synchrotron-peaked blazar Mrk\,421 using simultaneous LAXPC20 and SXT observations. The 100\,s binned LAXPC20 light curve obtained during 3--8 January 2017 reveals pronounced flux variability. The source exhibits a fractional variability amplitude of $F_{\mathrm{rms}} = 0.210 \pm 0.005$ in the SXT band and $F_{\mathrm{rms}} = 0.316 \pm 0.006$ in the LAXPC20 band. During this interval, the source reached a peak LAXPC20 count rate of 122.94\,counts\,s$^{-1}$, while the peak count rate in the SXT light curve is 26.79\,counts\,s$^{-1}$. This enabled us to carry out flux-resolved spectroscopy by dividing the 100\,s binned LAXPC20 light curve into ten flux states (S1--S10), each spanning a width of 8\,counts\,s$^{-1}$. For each flux state, simultaneous SXT and LAXPC20 spectra were extracted and fitted jointly. We find that the spectra in these states are well described by a synchrotron-convolved broken power-law, which provides a better fit than a log-parabola model. The low-energy particle index (index before the break) is found to cluster around two discrete values across flux states indicating two spectra states in the source. The break energy consistently moves to high energy with increase in flux level in these states. Time-resolved spectroscopy (10-ks segments) confirms that the flux histogram is best modelled as a double lognormal distribution and the index histogram is double normal. Inclusion of two additional long observations spanning 2017-2019 shows the same double-state behaviour on longer timescales. Together, the results indicate that Mrk\,421 routinely occupies two dominant spectral; in a leptonic synchrotron framework this can be explained by Gaussian-like fluctuations in acceleration conditions producing lognormal flux states.

Figures

Figures reproduced from arXiv: 2512.21885 by the authors.

Figure 1
Figure 1. One-hour binned SXT and LAXPC light curves of Mrk 421 ob [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Comparison of the reduced chi-square (χ 2 ν ) values obtained from the broken power-law and log-parabola fits for different spectral states of Mrk 421. The hydrogen column density was fixed at NH = 1.33 × 1020 cm−2 , as obtained from the LAB survey (Kalberla et al., 2005). Galactic absorption was modeled using the TBABS routine. A relative normalization constant was ap￾plied to the SXT spectrum to account for inter-… view at source ↗
Figure 3
Figure 3. Variation of the broken power-law spectral parameters of Mrk 421 with flux state (S1–S10). [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Correlations among the spectral parameters of the broken power-law model for Mrk 421. [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Variation of the break energy with flux state (S1–S10) obtained from flux resolved spectroscopy. [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: Comparison of the index distribution (left panel) and the flux(erg cm [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: Comparison of the index distribution (left panel) and the flux(erg cm [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

7 extracted references · 5 linked inside Pith

  1. [1996]

    nature 383, 319–320

    Extremely rapid bursts of TeV photons from the active galaxy Markarian 421. nature 383, 319–320. doi:10.1038/383319a0. Gaskell, C.M., 2004. Lognormal X-Ray Flux Variations in an Extreme Narrow-Line Seyfert 1 Galaxy. Astrophys. J. Lett. 612, L21–L24. doi:10.1086/424565. Ghisellini, G., Tavecchio, F., Foschini, L., Ghirlanda, G., Maraschi, L., Celotti, A., ...

    arXiv 2004

  2. [2001]

    The Astrophys- ical Journal 559, 187

    Simultaneous x-ray and tev gamma-ray observation of the tev blazar markarian 421 during 2000 february and may. The Astrophys- ical Journal 559, 187. Kushwaha, P ., Chandra, S., Misra, R., Sahayanathan, S., Singh, K.P ., Baliyan, K.S., 2016a. Evidence for Two Lognormal States in Multi- wavelength Flux Variation of FSRQ PKS 1510-089. Astrophys. J. Lett. 822...

    arXiv 2000

  3. [2010]

    VHEγ-ray emission of PKS 2155-304: spectral and temporal variability. Astron. Astrophys. 520, A83. doi:10.1051/0004-6361/ 201014484. Hota, J., Shah, Z., Khatoon, R., Misra, R., Pradhan, A.C., Gogoi, R., 2021. Understanding the X-ray spectral curvature of Mkn 421 using broad- band AstroSat observations. Mon. Not. R. Astron. Soc. 508, 5921–

    doi:10.1051/0004-6361/ 2021

  4. [2016]

    The Astrophysical Journal 831, 92

    Multifrequency photo-polarimetric webt observation campaign on the blazar s5 0716+ 714: Source microvariability and search for characteristic timescales. The Astrophysical Journal 831, 92. Biteau, J., Giebels, B., 2012. The minijets-in-a-jet statistical model and the rms-flux correlation. Astron. Astrophys. 548, A123. doi:10.1051/ 0004-6361/201220056,arXi...

    arXiv 2012

  5. [2019]

    Variability studies and modelling of the blazar PKS 2155-304 in the light of a decade of multi-wavelength observations. Mon. Not. R. Astron. Soc. 484, 749–759. doi:10.1093/mnras/stz027, arXiv:1901.01743. Falomo, R., Pian, E., Treves, A., 2014. An optical view of bl lacertae objects. The Astronomy and Astrophysics Review 22, 1–38. Gaidos, J.A., Akerlof, C....

    arXiv 1901

  6. [2201]

    McHardy, I., 2010

    doi:10.1093/mnras/staf620,arXiv:2505.23338. McHardy, I., 2010. X-Ray Variability of AGN and Relationship to Galactic Black Hole Binary Systems, in: Belloni, T. (Ed.), Lecture Notes in Physics, Berlin Springer Verlag, p. 203. doi:10.1007/ 978-3-540-76937-8_8. Misra, R., Roy, J., Yadav, J., 2021. An alternative scheme to estimate astrosat/laxpc background f...

    arXiv 2010

  7. [5934]

    Kalberla, P .M.W., Burton, W.B., Hartmann, D., Arnal, E.M., Bajaja, E., Morras, R., P ¨oppel, W.G.L., 2005

    doi:10.1093/mnras/stab2903,arXiv:2110.03344. Kalberla, P .M.W., Burton, W.B., Hartmann, D., Arnal, E.M., Bajaja, E., Morras, R., P ¨oppel, W.G.L., 2005. The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI. Final data release of the combined LDS and IAR surveys with improved stray-radiation corrections. As- tron. Astrophys. 440, 775–782. doi:10.1051/0004...

    arXiv 2005

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.