{"id":"7cab59e7-b136-4908-90e5-c65956794f56","arxiv_id":"2506.03540","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Photoionization models of Mrk 1018 show that broad-line flux changes track disc emission changes, not X-ray changes alone, and that highly ionized Fe K alpha traces the inner coronal region.","lead":"This paper uses photoionization simulations to test how changing an AGN's X-ray and disc emission affects its broad emission lines and the highly ionized iron line. It finds that X-ray changes alone barely affect the optical/UV broad lines, while changes in the disc emission can dim them, supporting the idea that changing-look AGN are powered by intrinsic changes in the accretion flow.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The coronal-origin conclusion for highly ionized Fe Kα is not uniquely determined: it rests on a degenerate radius–density choice and is never compared to observed X-ray line fluxes.","rationale":"The strongest part of the paper is the demonstration that BLR line fluxes are nearly insensitive to α_ox variations alone (Fig. 7): the grid is clean, the trend is monotonic and small, and different X-ray slopes (Fig. 19) and disk temperatures (Fig. 20) do not change the qualitative picture. The Mrk 1018 application correctly reproduces the direction of the observed dimming, and the authors candidly note the overestimation of the magnitude. The most load-bearing weakness is the coronal-origin claim for highly ionized Fe Kα. The model is deliberately set up with log U=3 and n_H=10^12 cm^-3, which forces the emitting region to be extremely close to the black hole (R_FeK=10^-2.5 R_BLR). But because U depends on r and n_H only through the product r^2 n_H, the same high-U state can be produced at a much larger radius if the density is lower. The paper's own Fig. 21 demonstrates this degeneracy: at n_H=10^9 cm^-3, H-like and He-like Fe K emission is already present at R_FeK=10^-1 R_BLR. The abstract's statement that the line 'originates in the coronal part of the accretion disk' is therefore not an output of the model but a restatement of the assumed radius/density combination. No observational X-ray data on Fe Kα in Mrk 1018 (or any CL AGN) are used to validate the predicted luminosity or line ratio, so the conclusion is currently untested. The BLR single-zone assumption is a real limitation, but the authors acknowledge it and it mainly affects the magnitude of the predicted line drops, not the qualitative conclusion. The coronal-origin claim is more consequential because it is a headline result and rests entirely on an unverified parameter choice. A concrete, feasible check is to compare the predicted Fe Kα luminosity and line ratio to X-ray spectra of Mrk 1018 in its two states, with n_H and R_FeK as free parameters. If the data prefer R_FeK > 10^-1 R_BLR or lower densities, the coronal designation would need to be revised. This check would settle whether the concern actually lands.","tokens_in":26083,"tokens_out":11312,"duration_ms":120737,"concrete_test":"Fit the 6–7 keV Fe Kα complex in XMM-Newton/NuSTAR spectra of Mrk 1018 in its bright (Sy1) and faint (Sy1.9) states, measuring the line luminosity or upper limit and the Fe XXV/Fe XXVI line ratio. Run Cloudy with n_H and R_FeK as free parameters (covering n_H=10^9–10^12 cm^-3 and R_FeK from 10^-3 to 1 R_BLR) and compare the predicted line luminosities and ratios to the observed values. If the data allow R_FeK > 10^-1 R_BLR, or require n_H < 10^12 cm^-3, the coronal-origin claim in the abstract is not supported. A simpler analytical check: compute the range of (r, n_H) pairs that yield log U=3 for the Mrk 1018 SEDs in Table 1 and show that R_FeK spans several orders of magnitude.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.2 and the abstract claim that highly ionized Fe Kα originates in the coronal region. The adopted parameters (n_H=10^12 cm^-3, R_FeK=10^-2.5 R_BLR, log U=3) are chosen to make H-like and He-like Fe ions dominant in Fig. 17. However, U=Q(H)/(4πr^2 n_H c) is degenerate in r and n_H: the paper's own Fig. 21 shows that at n_H=10^9 cm^-3, comparable H-like and He-like Fe K emission already appears at R_FeK=10^-1 R_BLR, roughly 30 times farther out. Thus, without an independent density constraint, the model does not force the emitting region to be at a coronal distance; a lower-density medium at a larger radius produces the same qualitative result. Moreover, no predicted Fe Kα luminosity, equivalent width, or line ratio is compared to observed X-ray spectra of Mrk 1018 or any CL AGN. The coronal-origin statement is therefore an input assumption, not an empirically tested output.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents Cloudy photoionization simulations of the response of broad-line region (BLR) optical/UV lines (Hα, Hβ, Mg II, He II) and the Fe Kα X-ray line complex to changes in the AGN continuum. Two scenarios are considered: (i) variation of the X-ray power-law only, with fixed disk emission, and (ii) broadband continuum variations as observed in the changing-look AGN Mrk 1018 (using SEDs from Noda & Done 2018). The main findings are that BLR line fluxes are largely insensitive to X-ray-only changes; for Mrk 1018 the simulations reproduce the observed brightening/dimming trends of the BLR lines, although with predicted flux drops that are considerably larger than observed; and that the highly ionized Fe Kα emission from H-like and He-like Fe ions responds strongly to the X-ray strength of the SED, which the authors interpret as evidence that this emission originates in the coronal region of the accretion disk. The paper also concludes that changing-look behavior in AGN is likely driven by intrinsic changes in the accretion rate.","tokens_in":26286,"tokens_out":4618,"duration_ms":50969,"significance":"If the main conclusions hold, the paper supports the view that changing-look AGN transitions reflect genuine changes in the accretion state rather than purely geometrical effects, and it identifies the highly ionized Fe Kα line as a potentially useful coronal diagnostic. The BLR insensitivity to X-ray-only continuum changes is convincingly demonstrated across a wide parameter grid (Figs. 3-7) and is a useful, robust result. The Fe K part, however, is more model-dependent: the conclusion that the emitting region is coronal rests on a chosen radius and density rather than on an independent constraint, and the predicted Fe K luminosities are never compared to observed X-ray line measurements. With appropriate qualification and a comparison to observed Fe K data, the paper would make a solid contribution to the changing-look AGN literature.","major_comments":[{"comment":"The conclusion that the highly ionized Fe Kα line originates in the coronal region is not uniquely determined by the modeling because the ionization parameter U = Q(H)/(4πr² n_H c) is degenerate in radius and density. The paper's own Fig. 21 shows that at n_H = 10^9 cm^-3, comparable H-like and He-like Fe K emission already appears at R_FeK = 10^-1 R_BLR, which is about 30 times farther out than the adopted R_FeK = 10^-2.5 R_BLR. Thus the adopted combination (n_H = 10^12 cm^-3, R_FeK = 10^-2.5 R_BLR) is not forced by the data; a lower-density medium at larger radius produces the same qualitative result. Without an independent density constraint, the model does not demonstrate that the emitting region is actually at a coronal distance. I recommend either adding such a constraint or explicitly softening the abstract and conclusion claims from 'origin is in the coronal part' to 'consistent with a coronal origin under the adopted parameters.'","section":"Section 4.2, Eq. (3), Fig. 21"},{"comment":"The predicted Fe Kα luminosities, equivalent widths, and line ratios are never compared to observed X-ray spectra of Mrk 1018 or any other changing-look AGN. This is a load-bearing gap because the paper's main novelty for the Fe K line is the claim that its variability traces coronal changes. The simulations currently show that a region with the chosen parameters would respond to X-ray flux variations, but they do not demonstrate that such a region exists or that its predicted line strength is consistent with X-ray observations. I recommend adding a quantitative comparison to available X-ray data (e.g., XMM-Newton or Suzaku spectra of Mrk 1018 or similar CL AGN) or, if that is not possible, clearly labeling the Fe K results as a parameter study rather than an empirical reproduction.","section":"Section 4.2 and Section 5"}],"minor_comments":[{"comment":"The text states 'we assumed that the BLR gas density is nH = 12 cm^-3', which appears to be a typo for n_H = 10^12 cm^-3. Please correct this.","section":"Section 4.1"},{"comment":"The legend labels contain 'Fluorscent' instead of 'Fluorescent' in Figures 17, 21, 22, and 23. Please correct the spelling.","section":"Figure 17 and Appendix Figures 21-23"},{"comment":"The text mentions 'Fe XXV and Fe XVI' as high-ionization lines; this is likely a typo, as Fe XVI is not a high-ionization line in this context and the intended reference is probably Fe XXV and Fe XXVI. Please verify and correct.","section":"Section 3.2 (discussion of Iwasawa et al.)"},{"comment":"The radius-luminosity relation of Bentz et al. (2009) has a reported scatter; including this uncertainty in the propagated Mrk 1018 ionization parameters (Table 1) would make the BLR line predictions more robust and would help assess the significance of the predicted line flux changes.","section":"Section 4.1, Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is already marked as accepted by ApJ, but my assessment is based solely on the content. The BLR response part is solid and the grid study is useful. The primary weakness is the Fe K coronal-origin claim, which currently rests on an assumed radius and density rather than on an empirically tested prediction. I would support publication after the authors either add a comparison to observed Fe K line measurements or substantially qualify the coronal-origin conclusion. This is a fixable issue, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a workmanlike Cloudy study with one clearly demonstrated result and one over-interpreted conclusion. Worth a serious referee, but the coronal-origin claim needs to be softened.\n\nThe systematic grid comparing X-ray-only SED changes to broadband changes is the real content. The demonstration that BLR lines stay nearly flat when the X-ray power-law is varied by orders of magnitude (Figs. 3–7) is solid and will be useful to people arguing against X-ray-only drivers of changing-look AGN. The Mrk 1018 application reproduces the observed dimming direction, and the authors are honest that the magnitude is off: line drops are factors of 10–26 while observed drops are smaller, and they attribute this to the single-zone assumption. That is fair.\n\nThe Fe K part is where I part ways with the abstract. The claim that highly ionized Fe Kα originates in the coronal region, closer in than the BLR, is presented as a finding, but it is partly built into the model. The density (10^12 cm^-3) and radius (10^-2.5 R_BLR) are chosen to make H-like and He-like ions dominate, and Fig. 21 shows that at lower density (10^9 cm^-3) the same qualitative emission appears at 10^-1 R_BLR. The U = Q/(4πr^2 n_H c) degeneracy means the result does not uniquely locate the region. What is missing is any comparison to observed X-ray line fluxes or equivalent widths for Mrk 1018 or any CL AGN. Without that, 'coronal origin' is an interpretive statement, not a measured output.\n\nThe single-zone BLR assumption is the other soft spot, but the authors flag it themselves and it is a standard simplification. No code or data artifacts are provided beyond the public Cloudy code; that is a minor issue for reproducibility.\n\nCitations look fine. The self-citations are to the authors' own previous Cloudy modeling, which is relevant context, not padding.\n\nBottom line: for someone modeling BLR variability and planning XRISM/Athena observations of Fe K, this paper is a useful reference and worth discussing in the reading group. It deserves peer review. I would not cite it as the last word on Fe K location, but I would cite it for the BLR insensitivity result. Recommend engage, with the Fe K conclusion tempered.","headline":"A workmanlike Cloudy study with a solid BLR-insensitivity result and an over-interpreted Fe K coronal-origin claim; worth refereeing, but the coronal location needs to be treated as model input, not output.","tokens_in":26903,"tokens_out":2196,"would_cite":true,"duration_ms":23461,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that changing-look AGN transitions are driven by intrinsic accretion-rate changes, and that the highly ionized Fe Kα line is a coronal diagnostic that tracks them.","keywords":["changing-look AGN","photoionization","broad emission lines","Fe K alpha","accretion disk corona","Eddington ratio","Mrk 1018","Seyfert transitions"],"falsifier":"Take a changing-look AGN such as Mrk 1018 through its bright and faint states with simultaneous optical monitoring and a high-resolution X-ray spectrum, measuring the 6.7 keV He-like Fe Kα line. The model predicts the ionized Fe Kα luminosity falls by about an order of magnitude as $\\log(L/L_{\\rm Edd})$ drops from $-1.14$ to $-2.27$; if the line flux instead stays constant or rises while the continuum dims, the coronal-origin and accretion-rate-driven picture is wrong. A second check: find a changing-look source with a large Hα disappearance but a constant 5100 Å continuum, which the paper's X-ray-only models say should not happen.","tokens_in":25824,"feed_emoji":"🔭","tokens_out":10221,"duration_ms":107877,"temperature":0.7,"pith_summary":"Changing-look AGN flip between spectral types on timescales far shorter than viscous timescales, and this paper asks what actually drives the flip. The authors run photoionization simulations of the broad-line region and the iron Kα region under two continuum-change scenarios: X-ray power-law variation alone, and realistic broadband changes in both disk and X-ray emission. They find that the optical/UV broad lines (Hα, Hβ, Mg II, He II) are nearly indifferent to X-ray-only changes, but respond strongly when the whole spectral energy distribution brightens or dims, as in the well-studied changing-look source Mrk 1018. The highly ionized Fe Kα emission, from H-like and He-like Fe ions, instead tracks the X-ray strength tightly and requires a high-ionization region closer to the black hole than the broad-line region. If correct, the result says changing-look events are genuine accretion-state changes rather than orientation or obscuration effects, and the ionized Fe Kα line is a probe of the corona.","feed_headline":"Accretion changes, not X-ray flares, drive changing-look AGN","feed_subtitle":"X-ray-only variation leaves broad lines almost untouched; matching Mrk 1018's dimming needs broadband disk change.","key_machinery":"The load-bearing tool is the Cloudy photoionization code (C23.01) applied to plane-parallel slabs of solar-composition gas with column density $10^{23}\\ \\mathrm{cm}^{-2}$. The input SED is a big blue bump plus an X-ray power law whose normalization is set by the X-ray-to-optical index $\\alpha_{\\rm ox}$; the ionization parameter $U = Q(H)/(4\\pi r_0^2 n_H c)$ sets the radiation strength. For the realistic Mrk 1018 case, the BLR radius is fixed by a radius-luminosity relation based on the 5100 Å luminosity, gas density is fixed at $n_H = 10^{12}\\ \\mathrm{cm}^{-3}$, and the Fe K emitting radius is scaled as $R_{\\rm FeK} = 10^{-1}$ to $10^{-3} R_{\\rm BLR}$, with $10^{-2.5} R_{\\rm BLR}$ giving the H-like and He-like dominance. The output quantities that carry the argument are the four broad-line luminosities and the decomposition of the Fe Kα flux into H-like, He-like, hot fluorescent (Fe XVIII–XXIII), and cold fluorescent ($\\le$Fe XVII) components.","core_discovery":"On the paper's own terms, the central discovery is a clear separation of drivers in AGN line variability. In a grid of photoionization models with the disk emission held fixed and only the X-ray power-law steepening (larger $\\alpha_{\\rm ox}$), the Hα, Hβ, Mg II, and He II fluxes change by factors of only about 1.2–2.4, far too little to make broad lines appear or disappear. When instead the full spectral energy distribution changes as in Mrk 1018's 2008–2016 transition—bolometric luminosity falling by a factor of 16, with the X-ray luminosity falling about two orders of magnitude—the same lines drop by factors of about 21 (Hα), 26 (Hβ), 24 (Mg II), and 10 (He II) from the highest to the lowest Eddington-ratio state, reproducing the observed dimming from Seyfert 1 to 1.9. For the Fe K region, the simulations show the H-like and He-like Fe components only become dominant when the emitting gas sits at radii near $10^{-2.5}\\,R_{\\rm BLR}$ with a high ionization parameter; at that location the ionized Fe Kα flux tracks the X-ray strength and falls by about an order of magnitude in the faint state. The paper therefore concludes that highly ionized Fe Kα is produced in the coronal region of the accretion disk, and that a change in the mass accretion rate is the likely trigger of the changing-look phenomenon.","pith_inferences":["Beyond the paper, the single-zone assumption is the place to look first if future data disagree: a BLR with a spread in densities and radii would dilute the predicted factor-of-20 line drops and could reconcile the simulations with changing-look sources that show milder line changes.","Beyond the paper, the result implies that a changing-look source caught in its faint state should still show a weak but detectable 6.7 keV Fe Kα component if the accretion flow has not fully shut off; searching for this residual emission is a direct test of the coronal-origin claim.","Beyond the paper, the proposed driver could be tested statistically: in a sample of changing-look AGN, the amplitude of the Hβ change should correlate with the amplitude of the 5100 Å continuum change, not with the X-ray change alone.","Beyond the paper, if the BLR radius lags the luminosity by a dynamical or recombination time, using a radius-luminosity relation at each epoch could misplace the clouds; time-resolved reverberation mapping during a changing-look event would separate this lag effect from intrinsic accretion change."],"forward_implications":["If broad-line disappearance in changing-look AGN requires broadband continuum change, then X-ray-only variability campaigns should not be expected to correlate one-to-one with Hα/Hβ type flips; a strong line transition without any 5100 Å change would point to an additional mechanism.","The ionized Fe Kα line becomes a practical coronal tracer: its 6.7–6.97 keV flux should rise and fall with the Eddington ratio, giving X-ray spectrometers a way to watch the corona evolve across a changing-look event.","The ordering of line responses found here—Hβ dropping most (about 26×) and He II least (about 10×)—gives a quantitative prediction for multi-line monitoring samples of changing-look AGN.","If changing-look phenomena are driven by mass accretion changes, Type 1 and Type 2 are not fixed orientation classes but states a single source can cycle through, so AGN demography needs to be treated as a time-dependent process."],"supporting_citations":[{"why":"Supplies the four Mrk 1018 SEDs and Eddington ratios that drive the broadband-variation models, plus the interpretation of the transition as a change in accretion flow.","marker":"Noda & Done (2018)"},{"why":"Provides the radius-luminosity relation (Eq. 4) used to move the BLR radius as the 5100 Å luminosity changes.","marker":"Bentz et al. (2009)"},{"why":"Cloudy C23.01 is the code with which all line fluxes are computed.","marker":"Gunasekera et al. (2023)"},{"why":"Background code release that Cloudy C23.01 builds upon; the models inherit its atomic data and photoionization treatment.","marker":"Ferland et al. (2017)"},{"why":"Disk-atmosphere density profile used to justify the standard n_H = 10^12 cm^-3 BLR density and gas parameters.","marker":"Adhikari et al. (2016)"},{"why":"Establishes the Fe Kα behavior the simulations reproduce: cold fluorescence at low X-ray flux, H-like/He-like dominance at high illuminating flux.","marker":"Ballantyne & Ross (2002)"},{"why":"Provides the observational Fe Kα-X-ray tandem variability in NGC 1566 used as a consistency check for coronal origin.","marker":"Liang et al. (2022)"},{"why":"Documents Mrk 1018's Hα disappearance together with declining X-ray flux, the changing-look case the broadband models target.","marker":"LaMassa et al. (2017)"},{"why":"Sample result that Hβ is the most variable changing-look line, compared against the ordering predicted by the simulated line drops.","marker":"Zeltyn et al. (2024)"}],"fun_headline_variants":["Accretion shifts, not X-ray flares, flip AGN types","Broadband disk dimming explains changing-look AGN","X-ray-only changes leave broad lines almost intact","Coronal iron line tracks X-ray, CL from accretion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the broad-line region can be represented by a single gas slab of fixed density $n_H = 10^{12}\\ \\mathrm{cm}^{-3}$ whose radius is set by the 5100 Å luminosity through a radius-luminosity relation, so that one ionization parameter describes all lines and the Fe K region; if the real BLR spans a range of densities and radii, or is not in photoionization equilibrium during a fast transition, the computed line drops and the inferred coronal radius would change.","fun_headline_variants_meta":{"raw":{"variants":["Accretion shifts, not X-ray flares, flip AGN types","Broadband disk dimming explains changing-look AGN","X-ray-only changes leave broad lines almost intact","Coronal iron line tracks X-ray, CL from accretion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1768,"prompt_tokens":1172,"completion_tokens":596,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":788,"completion_tokens_details":{"reasoning_tokens":528}},"tokens_in":788,"tokens_out":596,"duration_ms":7312,"temperature":1.0,"reasoning_tokens":528,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:01:11.096484+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a changing-look AGN such as Mrk 1018 through its bright and faint states with simultaneous optical monitoring and a high-resolution X-ray spectrum, measuring the 6.7 keV He-like Fe Kα line. The model predicts the ionized Fe Kα luminosity falls by about an order of magnitude as $\\log(L/L_{\\rm Edd})$ drops from $-1.14$ to $-2.27$; if the line flux instead stays constant or rises while the continuum dims, the coronal-origin and accretion-rate-driven picture is wrong. A second check: find a changing-look source with a large Hα disappearance but a constant 5100 Å continuum, which the paper's X-ray-only models say should not happen.","supporting_citations":[{"cited_title":"M., Yaqoob , T., & Kilgard , R","cited_arxiv_id":null,"evidence_quote":"Documents Mrk 1018's Hα disappearance together with declining X-ray flux, the changing-look case the broadband models target."}],"review_version":1}