{"id":"1ffac6ba-989c-4689-a46c-170bb69f87db","arxiv_id":"1908.08023","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"New NICER X-ray observations of NGC 4388 reveal a Compton-thin obscuring column of about 2.6 x 10^23 cm^-2, roughly half the column seen by NuSTAR in 2013, confirming rapid variability in the absorber.","lead":"Astronomers used NASA's NICER X-ray telescope to measure the thick cloud of gas hiding the black hole at the center of the galaxy NGC 4388. The cloud is dense but not thick enough to block all X-rays, and it appears to have thinned by about half since 2013, showing that such obscuration can change on human timescales.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The NICER column is a 105.6 ks average over ~16 months of a source known to have unveiling events; the coarse hardness check does not rule out a low-column fraction biasing the average, so the factor-of-two variability versus NuSTAR is not yet secured.","rationale":"The paper is a careful spectral analysis with two model families, MCMC uncertainties, and transparent limitations, and the core N_H measurement near 2.6e23 cm^-2 is probably robust to many systematics. The cold-neutral-gas concern identified by the reader is valid, but it is explicitly acknowledged in the final paragraph and partially tested: the authors report that replacing the neutral obscuration with ionized obscuration through zxipcf yields significantly worse fits. That reduces its force as the single most load-bearing issue. The least secure condition is instead temporal: the NICER N_H is fitted to a 105.6 ks exposure accumulated over about 16 months, while the comparison value is a single 2013 snapshot. NGC 4388 is known for short-timescale unveiling events, so a brief low-column interval included in the summed spectrum could bias the single-zone N_H downward, exactly in the direction of the claimed variability. The hardness-ratio check is too coarse, as the authors acknowledge, and no time-resolved spectral analysis is presented. A time-resolved spectral fit would directly settle whether the average is representative. Because the variability claim is central and this check is missing, I would adjust the verdict from ACCEPT to CONDITIONAL: the paper should be accepted only if the authors add such a time-resolved analysis or appropriately qualify the variability conclusion.","tokens_in":16305,"tokens_out":15115,"duration_ms":151928,"concrete_test":"Split the 105.6 ks NICER exposure into segments of roughly 10-20 ks each, fix the spectral shape parameters other than the normalization and N_H,LOS to the best-fit mytorus values from Table 1, and fit N_H,LOS independently in each segment. If all segment values are consistent with about 2.6e23 cm^-2 within 1-sigma, the time-averaged value is representative and the comparison with the 2013 NuSTAR snapshot stands. If the segments show N_H values clustering near the 2013 NuSTAR value or near zero, the average is a mixture of states and the variability claim should be reframed as a change in the distribution or duty cycle of column densities, not a simple factor-of-two decline.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central variability claim compares a single-epoch NuSTAR column (about 5e23 cm^-2 in 2013) to a column fitted from a 105.6 ks time-averaged NICER spectrum accumulated over roughly 16 months. NGC 4388 is known to exhibit short-timescale 'unveiling' events, and the NICER light curve in Figure 1 shows rare high-flux intervals. If even a small fraction of the total exposure sampled a low-column state, a single-zone absorber fitted to the summed spectrum will return an effective N_H biased below the dominant state, which would mimic the reported factor-of-two decline. The authors' only check is a hardness ratio that they themselves describe as having large errors; no time-resolved spectral fits are presented. Therefore the 'genuine variability' conclusion is not yet secured: the NICER column may be an exposure-weighted mixture of different column states rather than a representative 2017-2019 column. The cold-neutral-gas limitation is acknowledged and partially tested (zxipcf-based ionized obscuration is reported to fit significantly worse), so the temporal mixing issue is the least secure condition for the central claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a NICER X-ray spectral analysis of the Seyfert-2 galaxy NGC 4388 using a time-averaged 0.6–10 keV spectrum with 105.6 ks of net exposure accumulated between 2017 December and 2019 March. The authors fit the spectrum with two independent model families, pexmon and mytorus, each combined with an XSTAR ionized absorber, scattered power-law components, and three mekal plasma components. Both models give consistent line-of-sight neutral column densities, N_H ≈ 2.6 × 10^23 cm^-2, which are Compton-thin and broadly consistent with earlier Suzaku and XMM-Newton measurements but about half of the column reported from a 2013 NuSTAR observation; this difference is interpreted as genuine variability. The Fe K-alpha line is nominally resolved with a width corresponding to ~1200 km/s, suggesting an origin in the optical broad line region. The paper also reports ionized absorption in the Fe K band, derives lower limits on the Fe K production radius using rdblur, and presents simulated XRISM and Athena spectra as future tests of the inferred geometry.","tokens_in":16679,"tokens_out":5325,"duration_ms":55387,"significance":"If the central result holds, the paper provides a clean, well-characterized measurement of a Compton-thin, variable absorber in a nearby Seyfert-2 galaxy, with geometric inferences that link the obscuring/reflecting gas to the broad line region. The use of two independent physical model families with statistically consistent column densities is a genuine strength, as are the reported MCMC uncertainties and the explicit comparison with previous Suzaku, XMM-Newton, and NuSTAR analyses. The simulated XRISM/Athena spectra give falsifiable predictions for future high-resolution X-ray observatories. The paper is transparent about the limitations of the cold-neutral-gas assumption and about calibration-related residuals, which is commendable.","major_comments":[{"comment":"The central variability claim—that the obscuring column declined by roughly a factor of two between the 2013 NuSTAR epoch and the NICER epoch—is not yet secured because the NICER column is fitted to a 105.6 ks spectrum averaged over ~16 months. The paper itself notes rare high-flux intervals in the light curve and states that the hardness-ratio errors are large. If even a small fraction of the exposure caught an unveiling low-column state of the kind reported by Elvis et al. (2004), a single-zone absorber fitted to the summed spectrum will return an effective N_H biased below the dominant-state value, mimicking a genuine decline. The hardness-ratio check in Fig. 1 is not a substitute for time-resolved spectral fitting because the hardness ratio is not a direct column-density diagnostic and has large uncertainties. To support the 'genuine variability' conclusion, the authors should either fit spectra extracted from low-flux and high-flux intervals separately, add a two-zone absorption component to the time-averaged model, or present a quantitative demonstration that the observed sporadic high-flux intervals contribute negligibly to the fitted N_H. Absent that, the abstract and Section 4 should be rephrased to describe the difference as apparent or tentative.","section":"§2 and §4, Fig. 1"},{"comment":"The abstract states that the robustness of the variability is 'reinforced by the use of consistent models and procedures,' but the quoted NuSTAR measurements are not model-identical to the present fits. Kamraj et al. (2017) used pexrav, which lacks the self-consistent Fe K lines included in pexmon, and their mytorus fit used a single scattering/line component and a single inclination. Masini et al. (2016) used a mytorus implementation that is similar but not exactly the same as the one in this paper. Because pexmon/pexrav and decoupled mytorus geometries differ, part of the measured N_H difference could reflect model systematics rather than source variability. The wording in Section 4 ('could reflect a genuine reduction') is appropriately cautious, but the abstract and title are not. The authors should either re-fit the 2013 NuSTAR data directly with the exact models used here, or explicitly frame the variability as tentative pending a uniform reanalysis.","section":"§4, Kamraj et al. (2017) and Masini et al. (2016)"},{"comment":"The final paragraph appropriately acknowledges that pexmon and mytorus assume cold, neutral gas and that the gas may not be entirely neutral if it resides in the BLR. This is not a peripheral caveat: the quoted N_H values are neutral-equivalent columns, and a partially ionized absorber or reflector can bias the inferred column. The paper reports that exploratory xillver fits are equally good and that replacing the neutral absorber with zxipcf gives significantly worse fits, but those tests do not quantify how much the best-fit N_H changes under partial ionization, and the xillver and zxipcf geometries are not equivalent to the pexmon/mytorus setups. Given that the Fe K line width places the gas near or within the BLR, the systematic error on N_H from the ionization assumption could be substantially larger than the statistical errors of ±0.03 × 10^23 cm^-2. The authors should either add a quantitative systematic-error estimate for N_H under a partially ionized absorption/reflection model, or explicitly label the measured quantity as a 'neutral-equivalent column density' throughout the abstract and conclusions.","section":"§3.2, §3.3, and final paragraph of §4"}],"minor_comments":[{"comment":"In the model equation, 'contant[6]' should be 'constant[6]'.","section":"§3.2"},{"comment":"The text cites 'Vasudevan & Fabian (2006)' for the bolometric correction range, but the reference list contains Vasudevan & Fabian (2007) and (2009); please correct the citation year.","section":"§3.3"},{"comment":"The units of F_{0.6-10} and F_{0.6-10,unabs} are missing the exponent on cm: they should be erg cm^-2 s^-1.","section":"Table 1"},{"comment":"Please state how the 90% confidence lower limits in Table 2 were computed (e.g., Δχ^2 or MCMC percentile), since the table gives only the χ^2 value at the boundary.","section":"§3.4 and Table 2"},{"comment":"The sentence 'The pexmon model is not formally acceptable in a statistical sense' is important; please quantify the calibration-residual contribution to χ^2 (e.g., the χ^2 contribution of the 2–3 keV band) so readers can judge the impact on parameter uncertainties.","section":"§3.2"}],"recommendation":"major_revision","confidential_remarks":"This is a careful observational paper with transparent modeling and a plausible central result. The main issue is not the quality of the spectral fits but the strength of the variability claim: the time-averaging of the NICER data and the non-identical NuSTAR model comparisons both weaken the factor-of-two conclusion. The cold-neutral-gas assumption is acknowledged but its quantitative effect on N_H is not assessed. I believe the manuscript is publishable after the authors either provide time-resolved spectral evidence or temper the variability claim and add a systematic-error discussion. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The paper delivers a new NICER spectrum of NGC 4388 with a careful dual-model analysis; both pexmon and mytorus give N_H around 2.6e23 cm^-2, consistent with Suzaku/XMM and about half the NuSTAR 2013 value. That consistency is real. The second thing is that the variability claim is less secure than the abstract suggests. The NICER spectrum is a 105.6 ks sum taken over roughly 16 months, and the source is known to have unveiling events. A small amount of low-column exposure can pull a single-zone fit below the dominant state. The authors' only check is a hardness ratio with large errors; no time-resolved spectral fits are shown. So the factor-of-two decline versus NuSTAR could be an average artifact rather than a genuine change in the dominant column.\n\nWhat is new and good: the resolved Fe K line width, the ionized absorber detection, and the radius lower limits are all new. The two model families agree, and the MCMC errors are reported. The paper is transparent about the pexmon fit not being formally acceptable and about the cold-neutral-gas assumption; exploratory fits with ionized obscuration through zxipcf are reported as worse, which is the right way to address that concern. The literature citations to past Suzaku/XMM and NuSTAR measurements look fair, and the simulated XRISM/Athena spectra are clearly labeled as extrapolations from the fitted model.\n\nThe soft spots are proportionate. The time-averaging issue is the most important, because it directly bears on the central variability conclusion. The pexmon chi-square problem is minor and likely calibration-related around 2-3 keV. The F-test for the ionized absorber is simple but not a reason to reject the detection. The mekal components are ad hoc but not central. The authors already acknowledge the cold-neutral limitation in their final paragraph, so that is not a hidden flaw.\n\nWho should read this: AGN observers and anyone working on obscuration geometry or Fe K diagnostics. It is a useful measurement paper even if the variability conclusion needs to be softened. I would send it to peer review, and I would ask the authors to either do a time-resolved spectral analysis or explicitly present the NICER column as an exposure-weighted average and phrase the NuSTAR comparison accordingly.","headline":"NICER gives a clean, well-modeled column density for NGC 4388, but the headline variability versus NuSTAR rests on a 16-month average and is not yet secured.","tokens_in":17278,"tokens_out":2616,"would_cite":true,"duration_ms":25095,"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":"Obscuring gas around NGC 4388 halved in six years, new X-ray spectra show.","keywords":["active galactic nuclei","Seyfert 2","NGC 4388","X-ray spectroscopy","column density variability","Fe K line","broad line region","NICER"],"falsifier":"Re-observe NGC 4388 with NuSTAR and NICER simultaneously: if NuSTAR again measures a column near $5\\times 10^{23}$ cm$^{-2}$ while NICER measures about $2.6\\times 10^{23}$ cm$^{-2}$, the reported variability is an artifact of cross-calibration or model assumptions; if both instruments agree on the lower value, genuine variability is confirmed.","tokens_in":16092,"feed_emoji":"🔭","tokens_out":13074,"duration_ms":114248,"temperature":0.7,"pith_summary":"The paper analyzes a 105.6 ks time-averaged NICER spectrum of the Seyfert-2 active galaxy NGC 4388 to locate and characterize the gas that hides the central engine. Two independent spectral models, pexmon and mytorus, give consistent answers: the line-of-sight absorbing column is $N_H \\simeq 2.6\\times 10^{23}$ cm$^{-2}$, which is high but below the Compton-thick threshold. That value is only about half the column measured by NuSTAR in late 2013, so the paper argues the obscuration genuinely changed rather than being an artifact of different models. The neutral Fe K$\\alpha$ line is resolved at about 1200 km s$^{-1}$, consistent with an origin in the optical broad line region, and an ionized absorber is required at the $4\\sigma$ level. If right, this means the 'torus' is not a static parsec-scale structure but dynamic gas much closer to the black hole.","feed_headline":"NGC 4388's obscuring gas halved in six years","feed_subtitle":"Two model families agree the column is Compton-thin and half the 2013 value, pointing to changing gas.","key_machinery":"The argument rests on two independent spectral model families fitted to the same 105.6 ks NICER spectrum. The first, pexmon, is built on a cold, neutral reflection slab (pexrav) with Fe and Ni K lines whose strengths are tied to the reflection fraction by atomic physics. The second, mytorus, is a self-consistent torus model with linked absorption, scattering, and line-emission components; the version used here decouples the line-of-sight (near-side) and face-on (far-side) views so the obscuring geometry and the reflecting geometry can be constrained separately. Both model constructions include the same XSTAR photoionized absorber in front of the intrinsic cutoff power law, three mekal plasma components for the soft X-ray emission, and a small scattered polar component. Agreement between the two families on $N_H$ is the load-bearing result; a relativistic blurring function (rdblur) is then added to the pexmon model to turn the Fe K line width into lower limits on the inner radius of the neutral reflector.","core_discovery":"Using a summed, time-averaged NICER spectrum of NGC 4388 with a net exposure of 105.6 ks, the paper finds that the intrinsic power-law continuum from the central engine passes through a Compton-thin neutral absorber with column density $N_H = 2.67^{+0.02}_{-0.03}\\times 10^{23}$ cm$^{-2}$ under pexmon and $N_H = 2.64\\pm 0.03\\times 10^{23}$ cm$^{-2}$ under mytorus. These values are formally consistent, and both models find a small reflection fraction relative to the direct continuum, ruling out a Compton-thick line of sight. Compared with $N_H \\simeq 4$--$6.5\\times 10^{23}$ cm$^{-2}$ measured from 2013 NuSTAR data, the NICER value is about half as large, and the paper argues this difference is genuine variability because the same model families and procedures were used. The neutral Fe K$\\alpha$ line is resolved at $\\sigma = 40^{+10}_{-10}$ eV, corresponding to a velocity broadening around $1200^{+400}_{-400}$ km s$^{-1}$; interpreted as Keplerian motion, this places the line-forming region at roughly $10^4\\,GM/c^2$, inside the optical BLR. The data also require a photoionized (XSTAR) absorber at the $4\\sigma$ level, similar to warm absorbers seen in Seyfert-1 galaxies, and three soft thermal plasma components describe the low-energy spectrum.","pith_inferences":["If the obscuring gas is part of the optical BLR, then the Seyfert-1/Seyfert-2 distinction in NGC 4388 is a line-of-sight alignment through the same close-in cloud population, rather than a difference in the presence of a parsec-scale torus; future spectropolarimetry or Fe K reverberation could test this directly.","A natural monitoring experiment is to track $N_H$ with NICER on weeks-to-months timescales; if the gas is clumpy BLR material, the column should show stochastic jumps and possibly brief Compton-thick excursions, which would distinguish cloud crossings from a smooth radial gradient.","The paper's assumed radial ordering, with ionized absorber interior to the neutral gas, could be tested by catching a continuum flare: the warm absorber's ionization parameter should respond before the neutral column changes, yielding a direct distance estimate for each component."],"forward_implications":["With the line of sight Compton-thin, the strong Fe K edge and apparent reflection cannot come from a Compton-thick torus; the mytorus fits favor scattering inside the near-side absorber, though a degenerate far-side contribution remains.","Because the 2013 NuSTAR measurements and the 2018-2019 NICER measurements were made with consistent model families, the factor-of-two difference is evidence that the column density varies on timescales of years, consistent with the much faster variability reported for this source in 2004.","A resolved neutral Fe K$\\alpha$ line with $\\sigma \\simeq 40$ eV and implied velocity $\\sim 1200$ km s$^{-1}$ is consistent with an origin in the optical broad line region, tying the line to gas much closer than a parsec-scale torus.","Relativistic blurring fits place lower limits on the inner radius of the neutral reflector from $r \\geq 270\\,GM/c^2$ to $r \\geq 4500\\,GM/c^2$ depending on inclination, allowing the obscuration to begin well inside the classical torus.","Simulations in the paper show XRISM/Resolve should detect dynamical line broadening in 100 ks if the reflector sits at $1600\\,GM/c^2$, and Athena/X-IFU should detect factor-of-two column changes in 10 ks."],"supporting_citations":[{"why":"The 2013 NuSTAR measurement that found roughly twice the absorbing column; the comparison anchor for the variability claim.","marker":"Kamraj et al. 2017"},{"why":"Independent NuSTAR analysis using the same mytorus implementation that also found a higher column, reinforcing genuine variability.","marker":"Masini et al. 2016"},{"why":"XMM-Newton measurements of $N_H$ in broad agreement with the NICER value, anchoring the historical stability.","marker":"Beckmann et al. 2004"},{"why":"Suzaku measurements of $N_H$ in broad agreement with the NICER value.","marker":"Shirai et al. 2008"},{"why":"Reported hour-scale absorption variability in NGC 4388, the precedent for a rapidly changing absorber.","marker":"Elvis et al. 2004"},{"why":"Provides the pexrav cold-reflection model on which pexmon is built.","marker":"Magdziarz & Zdziarski 1995"},{"why":"Adds Fe and Ni K lines to pexrav with atomic-physics-linked strengths, defining the pexmon model used here.","marker":"Nandra et al. 2007"},{"why":"Introduces the mytorus model for self-consistent torus absorption, scattering, and lines.","marker":"Murphy & Yaqoob 2009"},{"why":"Provides the mytorus model tables and decoupled geometry used for the second family of fits.","marker":"Yaqoob et al. 2010"},{"why":"Shows that narrow Fe K lines in about half of Seyfert-1s are consistent with BLR production, supporting the paper's line-origin interpretation.","marker":"Shu et al. 2010"}],"fun_headline_variants":["NICER shows NGC 4388's gas column halved since 2013","Seyfert-2's veil thins: NGC 4388 column drops to half","Compton-thin confirmed: NGC 4388's column halves vs NuSTAR","NICER sees variable obscuration in NGC 4388: column halves","Two models, one answer: NGC 4388's absorber is Compton-thin and halved"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The models assume the obscuring and reflecting gas is cold and neutral with a single uniform column along the line of sight, so if the gas is partly ionized, as it would likely be if it sits in the broad line region, the reported column density could be systematically biased.","fun_headline_variants_meta":{"raw":{"variants":["NICER shows NGC 4388's gas column halved since 2013","Seyfert-2's veil thins: NGC 4388 column drops to half","Compton-thin confirmed: NGC 4388's column halves vs NuSTAR","NICER sees variable obscuration in NGC 4388: column halves","Two models, one answer: NGC 4388's absorber is Compton-thin and halved"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2936,"prompt_tokens":1182,"completion_tokens":1754,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":798,"completion_tokens_details":{"reasoning_tokens":1642}},"tokens_in":798,"tokens_out":1754,"duration_ms":13196,"temperature":1.0,"reasoning_tokens":1642,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:51:54.877417+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observe NGC 4388 with NuSTAR and NICER simultaneously: if NuSTAR again measures a column near $5\\times 10^{23}$ cm$^{-2}$ while NICER measures about $2.6\\times 10^{23}$ cm$^{-2}$, the reported variability is an artifact of cross-calibration or model assumptions; if both instruments agree on the lower value, genuine variability is confirmed.","supporting_citations":[],"review_version":1}