{"id":"f4ccf126-544a-47bd-a1f3-d4f1676a4189","arxiv_id":"2507.02195","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"XRISM's precision X-ray detector shows the iron K-alpha line of Centaurus A is a narrow doublet plus a broad component, implying emitting gas spans roughly 0.001 to 10 parsecs.","lead":"This paper examines the finest X-ray spectrum yet taken of the iron line in Centaurus A, the nearest active galaxy, and finds that the line contains two narrow peaks plus a broad component. The broad component indicates that the gas producing the line is spread across a huge range of distances, from about a thousandth of a parsec to tens of parsecs from the black hole.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The q≈2 and extended-emission conclusion rests on the internally inconsistent rdblur*MYTorusL convolution; q is likely biased by the model mismatch, and the required k=4.1 shows the reprocessor is not self-consistent.","rationale":"The reader's weakest assumption identifies the rdblur*MYTorusL inconsistency, and I agree that it is the central model weakness. The broad-line detection and the resolved narrow doublet are supported by the data and are not in question. What is load-bearing is the quantitative interpretation: the claim that the line shape requires a single extended emission region with q≈2 spanning four decades in radius. That conclusion depends entirely on the combined rdblur*MYTorusL model, whose internal inconsistency is acknowledged by the authors in Section 3.2. The formal error bar on q (1.99 ± 0.03) is misleading because it does not include the systematic offset from using incompatible emitter geometries. The additional requirement of a large relative normalization k=4.1, only partially explained by variability, strengthens the concern that the model is compensating for its own inconsistency rather than measuring an astrophysical emissivity profile. I do not see a reason to reject the paper: the observational result, including the FWHM measurements and the need for a broad component, is solid and worth publishing. The physical interpretation, however, should remain conditional until a self-consistent forward-modeling test is performed. Therefore the reader's CONDITIONAL verdict is appropriate and my read does not change it.","tokens_in":23392,"tokens_out":8137,"duration_ms":110852,"concrete_test":"Forward-model test: generate a mock XRISM/Resolve spectrum with a Monte Carlo reprocessor (or a public ray-tracing/torus code) using a known emissivity law q_true=2 over 10^-3 to 10 pc at i=24° and NH~1e23 cm^-2, then fit it with rdblur*MYTorusL exactly as in the paper. If the recovered q differs from q_true by more than the reported ±0.03, the q measurement is dominated by the model mismatch. A simpler side check: replace MYTorusL in Model A with a narrow Hölzer doublet and refit; if the fit statistic and q are nearly unchanged, the extended-torus part of the model is irrelevant to the q constraint.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most striking physical result is not the broad Fe Kα detection (which is robust) but the claim that a single emission region with emissivity index q≈2 extending from ~10^-3 pc to ~10 pc is required (Sections 3.2 and 5). This claim is produced by convolving a MYTorusL line profile with rdblur. The two models are geometrically inconsistent: MYTorusL generates a rest-frame line from an extended toroidal reprocessor, while rdblur re-assigns Doppler shifts according to a disk power-law emissivity r^-q between Rin and Rout. The authors explicitly acknowledge that 'the information about the location from which the Fe Kα photons were emitted in the MYTorusL model is lost' (Section 3.2). In the Compton-thin regime the pre-blurred MYTorusL line is nearly a narrow 2:1 doublet, so the broad-to-narrow ratio is essentially set by rdblur's assumed emissivity law. A mismatch between the actual radial emissivity of the torus and the assumed r^-q power law will be absorbed into q, making the quoted ±0.03 uncertainty purely statistical and not representative of model error. Model A also requires k=4.1 (Section 3.2), and the authors' variability argument accounts for only ~2.3 of this factor, so the model cannot simultaneously reproduce the line normalization and shape. Until q is recovered from a model that computes Fe Kα emissivity and kinematics self-consistently, the astrophysical conclusion q≈2 is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the Fe Kα complex in Cen A using the first XRISM/Resolve spectrum, with a Chandra-based estimate of contaminating emission within the XRISM beam. The authors report that the line consists of narrow Fe Kα1/α2 cores with FWHM ≈ 480 km/s plus a broad component with FWHM ≈ 4300 km/s of comparable flux. They then fit the 6.0–6.8 keV band with phenomenological Gaussians and with several combinations of the MYTorusL reprocessor table model blurred by the rdblur disk-blurring kernel. The preferred single-component model (model A) yields an emissivity index q = 1.99 ± 0.03, inner radius ≈ 5.4×10^2 r_g, outer radius ≈ 6.6×10^6 r_g, and inclination 24+13−7 deg, but requires a relative normalization k = 4.1. Alternative models with two or three components, with reverberation-mapping radius priors, and with fixed q = 3 are also explored. The paper concludes that the line shape requires an extended emission region from ~10^−3 pc to ~10 pc and that q ≈ 2.","tokens_in":23743,"tokens_out":4759,"duration_ms":56419,"significance":"If the q ≈ 2 result and the derived radial extent are robust, the paper would provide a rare, high-resolution X-ray constraint on the location of Fe Kα-emitting gas in a nearby AGN, with implications for coronal geometry and the relationship between the X-ray and optical/IR line-emitting structures. The observational work is careful in several respects: the broad-line detection is supported by clear residuals in the two-Gaussian fit (Fig. 3) and by BIC comparisons; the contamination subtraction using simultaneous Chandra data is a strength; and the authors explicitly enumerate the model caveats in Sections 3.2 and 4. However, the central astrophysical inference — that a single emission region with q ≈ 2 extending over four decades in radius is required — rests on the rdblur*MYTorusL convolution, a model combination the paper itself describes as internally inconsistent. In its current form the paper does not establish the q ≈ 2 claim beyond the model's assumptions, so the significance of the headline result is not yet realized.","major_comments":[{"comment":"The central claim that the Fe Kα line requires an emission region extending from ~10^−3 pc to ~10 pc with emissivity index q = 1.99 ± 0.03 is derived from convolving the MYTorusL line profile with rdblur. As the authors themselves note in Section 3.2 and again in Section 4, rdblur assumes the pre-blurred line is emitted at a single radius, whereas MYTorusL produces a line from an extended toroidal reprocessor so that 'the information about the location from which the Fe Kα photons were emitted in the MYTorusL model is lost'. In the Compton-thin regime relevant here, the unblurred MYTorusL line is close to the narrow 2:1 laboratory doublet, so the broad-to-narrow flux ratio and the line wings are essentially set by the assumed r^−q emissivity law in rdblur. The quoted ±0.03 uncertainty is therefore purely statistical and does not include the dominant model systematic. The q ≈ 2 conclusion is not established until it is recovered from a model that computes the Fe Kα line emissivity and Doppler/relativistic kinematics self-consistently (e.g., a Monte Carlo reprocessor with an explicit radial emissivity), or until the authors quantify the bias in q, Rin, and Rout caused by the single-radius assumption, for example by fitting simulated spectra generated with a known emissivity law.","section":"Section 3.2, model A and Table 2"},{"comment":"Model A requires a relative normalization k = 4.1+6.0−1.3, and the authors' own variability estimate based on the Swift/BAT light curve (Section 3.2 and Fig. 5) accounts for only a factor 2.3 ± 0.8 of this offset. This means that at k = 1 the MYTorusL line flux is significantly underpredicted, and the model is not a self-consistent physical description of both the line shape and the line normalization. Because the radial emissivity inference is drawn from the same model, the large k signals a mismatch between the assumed reprocessor and the actual emitting region, further weakening the q ≈ 2 claim. The paper lists possible explanations (past variability, non-solar abundances, geometry differences), but these are not tested; the revision should either remove the k-unity physical interpretation, add an explicit caveat that q is derived from a model whose normalization is ad hoc, or demonstrate with a concrete test that the normalization deficit does not bias q.","section":"Section 3.2, relative normalization k = 4.1"},{"comment":"The statement that 'the Fe Kα line in Cen A has a complex shape, which can only be explained by contribution from emission regions distributed over a wide range of radii, from ∼10^2−3 r_g ∼ 10^−3 pc to >∼10^6 r_g, which corresponds to >∼10 pc' is too strong given the acknowledged model inconsistencies of model A and the fact that the paper itself shows an alternative model (model D, with q = 3) that reproduces the profile with three components, not a single continuous radial distribution. The data robustly require a broad line, but the specific radial extent and especially the continuous extended region are model-dependent. The conclusions should be rephrased to say that the line shape is consistent with, but does not uniquely require, a single extended region with q ≈ 2, and the systematic uncertainty from the rdblur*MYTorusL combination should be reflected in the abstract and conclusions.","section":"Section 5, first paragraph"}],"minor_comments":[{"comment":"The statement that 'Additional Gaussian profiles only provided minimal improvements in the fit statistic and were ruled out as necessary components of the line profile by the BIC' is not quantified; please report the C-statistic and BIC values for the models with one, two, three, and four Gaussian components, or state that they are available as supplementary material.","section":"Section 3.1"},{"comment":"Several entries in Table 2 are ambiguous because a missing lower or upper error is not explicitly flagged; for example, model C lists Rin = 5.7+4.5 × 10^4 r_g with no lower error, and model D lists Rout = 8.3 × 10^8 r_g with no errors at all. Please add a systematic notation, such as a footnote stating that a missing bound means the parameter is unconstrained in that direction.","section":"Table 2"},{"comment":"The paper states in Section 2.3 that all parameters are presented with 1σ uncertainties, but in Section 3.2 the black hole mass error is quoted as a 3σ error and then rescaled; please clarify the convention at the first use of uncertainties (including asymmetric errors) so the reader knows which quantities are 1σ and which are not.","section":"Section 2.3 and Section 3.2"},{"comment":"The notation 'rdblur * atable{mytl_V000HLZnEp000_v01.fits}' is difficult to parse in prose; after the first use, please define a shorthand name (e.g., 'the MYTorusL table model') and then use that shorthand consistently throughout the text.","section":"Section 3.2"},{"comment":"The caption sentence 'The 4.1+6.0−1.3 multiple of that is depicted as a dashed line' is grammatically incomplete and unclear; please specify that the dashed line is the average BAT count rate multiplied by the best-fit relative normalization k, and that the pink region is the corresponding 1σ range.","section":"Fig. 5 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially interesting and the data analysis is generally careful, but the headline result (q ≈ 2 and the extended radial distribution) is built on the rdblur*MYTorusL combination that the authors themselves acknowledge to be internally inconsistent. The revision should either obtain q from a self-consistent model, or substantially soften the abstract/conclusion claims and explicitly present the q and radius measurements as conditional on the model assumptions. The use of MYTorusL is legitimate and the coauthor overlap with one of the model developers is disclosed through the reference list; I do not see a conflict issue that needs action."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper has a solid new observational result and a wobbly physical interpretation. XRISM/Resolve clearly separates the Fe Kα1/α2 doublet at 480 km/s and requires an additional broad component at 4300 km/s with comparable flux. Cen A is the nearest AGN, so this is a reference measurement. The authors do the appropriate checks—Chandra contamination subtraction, two RMFs, BIC comparisons—and the broad line survives. Credit where it is due.\n\nThe soft spot is the step from line shape to geometry. The q≈2 and Rin/Rout constraints come from convolving MYTorusL with rdblur. As the authors state in Section 3.2, rdblur assumes emission at a single radius while MYTorusL's line comes from an extended torus; the radius information is lost. The combination is not self-consistent. I think the stress-test note is right: q is likely absorbing the mismatch between the actual emissivity of the torus and the assumed r^-q law. The quoted ±0.03 is statistical only. Model A also needs k=4.1, and past variability accounts for only ~2.3 of that, so the normalization is not explained.\n\nThat said, the paper does not hide this. Section 4 is unusually explicit about the modeling caveats, and the abstract flags the k issue. The claim that q≈2 is \"required\" is too strong; the data require a broad component, and q≈2 is one way to fit it. Model D shows q=3 works with three components. So the headline should be the resolved doublet and broad-line detection, not q≈2.\n\nWho is this for? Anyone working on AGN X-ray spectroscopy or Cen A specifically. The measurement is a genuine step forward even if the physical model is provisional. It deserves a serious referee because the data are unique and the caveats are clearly spelled out. I would accept it with pressure to soften the interpretation language, not desk-reject.","headline":"Resolved Fe Kα doublet and a broad component in Cen A are real and worth knowing; the q≈2 extended-emission interpretation rests on a model combination the authors themselves admit is inconsistent.","tokens_in":24364,"tokens_out":1484,"would_cite":true,"duration_ms":16665,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The Fe Kα line in Centaurus A demands an emission region spanning roughly 0.001 to 10 parsecs.","keywords":["active galactic nuclei","X-ray spectroscopy","Fe K alpha","Centaurus A","XRISM/Resolve","emissivity index","reverberation mapping","broad-line region"],"falsifier":"A decisive test would be a fully self-consistent model that applies relativistic kinematics and ray tracing in situ to an extended toroidal reprocessor: if such a model fits the same XRISM/Resolve spectrum with $q=3$ or with a different radial extent, the paper's core claims would fail. A simpler check is whether the broad 4300 km/s component survives a re-fit with an independently recalibrated Resolve line-spread function; removing that component would collapse the argument for emission out to tens of parsecs.","tokens_in":23178,"feed_emoji":"🔭","tokens_out":10074,"duration_ms":108077,"temperature":0.7,"pith_summary":"Using the high-resolution XRISM/Resolve spectrum of the nearest active galactic nucleus, the paper shows that the Fe Kα line is not a single narrow feature but a resolved doublet with FWHM of about 480 km/s sitting on a broad component with FWHM of about 4300 km/s and comparable flux. The authors argue that this line shape can only be produced by an emission region extending from about $10^{-3}$ pc to about $10^{1}$ pc, with emissivity falling as $r^{-q}$ and $q\\approx 1.99$. If correct, this is one of the first direct measurements of an AGN line-emitting region spanning four orders of magnitude in radius, from the broad-line-region scale to the torus or molecular-disk scale. The fits also favor a low inclination of about $24^{\\circ}$, consistent with the jet orientation, rather than the high inclination implied by the Seyfert-2 classification.","feed_headline":"Iron line of Cen A spans 0.001 to 10 parsecs","feed_subtitle":"XRISM/Resolve reveals a narrow Fe Kα doublet plus an equally bright broad component, pointing to an extended corona.","key_machinery":"The load-bearing ingredients are the laboratory Fe Kα1/Kα2 line shape (energies, widths, and 2:1 intensity ratio), the roughly 5 eV resolution of the XRISM/Resolve microcalorimeter at 7 keV that separates the doublet and the broad wings, and the composite spectral model rdblur * MYTorusL. In that model, MYTorusL supplies a pre-broadened toroidal line profile including the Compton shoulder, while rdblur applies Doppler and relativistic blurring for a disk surface between an inner and outer radius with radial emissivity $\\epsilon(r)\\propto r^{-q}$. Model comparisons use the C-statistic for fitting and the Bayesian information criterion for penalizing extra components, and a multiplicative constant $k$ absorbs delays, abundances, and geometry differences between the assumed and true reprocessor.","core_discovery":"The paper's central claim is that the Fe Kα complex of Centaurus A contains a narrow doublet and a broad base whose combined profile requires a radially extended, roughly uniform reprocessor. Two narrow peaks at the laboratory energies of Fe Kα1 and Fe Kα2 are resolved with FWHM $(4.8\\pm0.2)\\times10^2$ km/s each, while a broad component with FWHM $(4.3\\pm0.3)\\times10^3$ km/s carries a similar flux. Modelling the line with a blurred toroidal-reprocessor profile (rdblur applied to MYTorusL), the authors find a best fit with inner radius $5.4^{+7.8}_{-4.0}\\times10^2\\,r_g$ ($1.4^{+2.1}_{-1.1}\\times10^{-3}$ pc) and outer radius $6.6^{+13.8}_{-3.2}\\times10^6\\,r_g$ ($17^{+36}_{-9}$ pc), emissivity index $q=1.99\\pm0.03$, and inclination $24^{+13}_{-7}$ degrees. If $q$ is fixed to the commonly assumed value of 3, the same line profile can only be reproduced with three separate emitting components, which the authors take as evidence that the fitted $q\\approx2$ is intrinsic rather than an artifact of one parameterization.","pith_inferences":["A direct testable extension is that other nearby AGNs observed with XRISM/Resolve should show similar narrow-doublet-plus-broad profiles with $q\\approx2$ if an extended corona is the norm rather than a lamppost.","The large normalization factor $k\\approx4$ could be a delayed response of the line to a brighter past continuum; a monitoring campaign tracking line flux and continuum simultaneously would test whether variability alone explains it.","The apparent ring-like inner component described in the appendix, if confirmed, could be a tidally disrupted structure; searching for similar narrow radial annuli in other high-resolution AGN spectra would provide a test.","Fitting the Fe Kβ line and the broad-band continuum, as the authors indicate they plan to do with NuSTAR and XRISM/Xtend, could break the degeneracy among $q$, inclination, and $k$ that limits the present fits."],"forward_implications":["If $q\\approx2$ is real, the ionizing corona must be extended, or the reprocessor must bend so that the illuminating flux does not decay as fast as $r^{-2}$, changing standard assumptions about AGN disk-corona geometry.","The fitted outer radius of roughly 10 pc puts the outer line-emitting region within reach of ALMA or JWST, offering an independent, spatially resolved test of the spectral fit.","The absence of a comparably broad component in optical and infrared lines implies that the broad Fe Kα emission traces gas that is hidden, dust-free, or otherwise invisible at longer wavelengths.","The low inclination of about $24^{\\circ}$ agrees with jet and inner-disk measurements and implies that Cen A's Seyfert-2 appearance comes from line-of-sight obscuration rather than an edge-on disk.","The result extends to a second object the pattern seen in NGC 4151, where the Fe Kα line also required emission from multiple radii, suggesting such complex line profiles may be common in nearby AGNs."],"supporting_citations":[{"why":"Supplies the laboratory energies, widths, and intensity ratio of the Fe Kα1/Kα2 doublet that fix the narrow-line centroids and the broad-line peak energy.","marker":"Hölzer et al. 1997"},{"why":"Introduces the MYTorus model of toroidal reprocessing, the basis of the MYTorusL line profile used in all the fits.","marker":"Murphy & Yaqoob 2009"},{"why":"Provides the updated MYTorusL table with the Hölzer Fe Kα profile at 2 eV resolution used for the XRISM data.","marker":"Yaqoob 2024"},{"why":"Gives the reverberation-mapping radii that define the two outer components of model C and provide the comparison for the fitted radial extents.","marker":"Iwata et al. 2024"},{"why":"Reports the analogous XRISM/Resolve Fe Kα analysis of NGC 4151 with three components, which motivates the multi-component models applied here.","marker":"XRISM Collaboration et al. 2024"},{"why":"Describes the Resolve microcalorimeter whose spectral resolution is what allows the doublet and broad base to be separated.","marker":"Porter et al. 2024"},{"why":"Provides prior Chandra/HETG measurements of the Cen A Fe Kα line, column density, and redshift that this analysis refines.","marker":"Bogensberger et al. 2024b"}],"fun_headline_variants":["Cen A's Fe Kα line covers 0.001–10 pc","Fe Kα doublet and broad base: Cen A's extended reprocessor","q≈2 emissivity fits Cen A's Fe Kα line shape","Cen A's iron line demands reprocessor from 0.001 pc to 10 pc","Broad Fe Kα component not seen in optical: extended corona?"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inferred radial extent, $q\\approx2$, and inclination assume that stacking rdblur (which treats the pre-blurred line as emitted at a single radius) on top of MYTorusL (which computes the line from an extended torus) does not seriously bias the fit, an inconsistency the authors themselves flag as potentially offsetting some best-fit parameters.","fun_headline_variants_meta":{"raw":{"variants":["Cen A's Fe Kα line covers 0.001–10 pc","Fe Kα doublet and broad base: Cen A's extended reprocessor","q≈2 emissivity fits Cen A's Fe Kα line shape","Cen A's iron line demands reprocessor from 0.001 pc to 10 pc","Broad Fe Kα component not seen in optical: extended corona?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001243,"raw_usage":{"total_tokens":5203,"prompt_tokens":1151,"completion_tokens":4052,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":767,"completion_tokens_details":{"reasoning_tokens":3948}},"tokens_in":767,"tokens_out":4052,"duration_ms":30498,"temperature":1.0,"reasoning_tokens":3948,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:35:07.364973+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a fully self-consistent model that applies relativistic kinematics and ray tracing in situ to an extended toroidal reprocessor: if such a model fits the same XRISM/Resolve spectrum with $q=3$ or with a different radial extent, the paper's core claims would fail. A simpler check is whether the broad 4300 km/s component survives a re-fit with an independently recalibrated Resolve line-spread function; removing that component would collapse the argument for emission out to tens of parsecs.","supporting_citations":[{"cited_title":"D., & Yaqoob, T","cited_arxiv_id":null,"evidence_quote":"Introduces the MYTorus model of toroidal reprocessing, the basis of the MYTorusL line profile used in all the fits."},{"cited_title":"2024, PASJ, 76, 923,","cited_arxiv_id":null,"evidence_quote":"Gives the reverberation-mapping radii that define the two outer components of model C and provide the comparison for the fitted radial extents."},{"cited_title":"S., Kilbourne, C","cited_arxiv_id":null,"evidence_quote":"Describes the Resolve microcalorimeter whose spectral resolution is what allows the doublet and broad base to be separated."}],"review_version":1}