{"id":"7b7b1149-5321-4e3f-8aba-0cf06f78bd46","arxiv_id":"2411.15357","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"The Lyα1/Lyα2 doublet ratio of H-like ions is proposed and demonstrated as a plasma column density diagnostic, decreasing from 2 to 1 as column density increases from about 1e20 to 1e24 cm^-2.","lead":"This paper shows that the brightness ratio of two closely spaced X-ray lines from hydrogen-like ions changes with the amount of gas the light passes through, from 2:1 in thin gas to 1:1 in thick gas. The authors test this as a new way to measure column densities in X-ray plasmas and demonstrate it on XRISM data from the binary Cen X-3.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The column-density inference from the Lyα ratio is sensitive to the assumed microturbulent velocity; Cen X-3 uses zero turbulence, so a plausible wind/turbulence of a few hundred km/s could shift the inferred NH by a factor of several, which the paper does not quantify.","rationale":"The paper's qualitative central claim is physically reasonable: the resolved Lyα fine-structure ratio should respond to line optical depth, and Cloudy is a well-established code with publicly available model data. The XRISM observation genuinely resolves the Fe XXVI doublet, and the observed ratio 1.35 ± 0.11 is a real measurement. The concern is not with the code or the observation but with the unquantified sensitivity of the column-density mapping to the assumed microturbulent velocity, which directly affects line trapping. The reader's weakest-assumption analysis already identified environment dependence, especially turbulence and geometry; my stress-test sharpens this into a specific mechanism that can be tested by a small grid of Cloudy runs. If the turbulence sensitivity is large, the diagnostic remains useful as an order-of-magnitude column indicator but cannot support a precise NH without an independent turbulence measurement. Because the reader already issued a CONDITIONAL verdict based on missing uncertainty propagation, my concern reinforces that verdict rather than moving it, hence UNCHANGED.","tokens_in":5747,"tokens_out":3695,"duration_ms":39450,"concrete_test":"Rerun the Cen X-3 Cloudy grid of Figure 4 with microturbulent velocity vturb = 0, 50, 150, and 300 km/s, holding density, ionization parameter, continuum shape, and plane-parallel geometry fixed. For each grid, read off the NH where Lyα1/(Lyα2+M1) = 1.35 ± 0.11. If the inferred NH changes by more than a factor of about 2 between vturb = 0 and 150 km/s, then the paper must either constrain the turbulence from other observed line widths or explicitly present the diagnostic as order-of-magnitude only.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim for Cen X-3 is the mapping from the observed Lyα1/Lyα2 = 1.35 ± 0.11 to NH ≈ 2×10^22 cm^-2 through a Cloudy model with v_turb = 0. In line radiative transfer, the line-center optical depth is inversely proportional to the Doppler b parameter; for a fixed ionic column, larger b means less photon trapping, so the ratio-versus-NH curve moves to higher NH as turbulence increases. The paper itself shows this environment sensitivity by adopting 150 km/s microturbulence for the Perseus model and zero for Cen X-3, but it never computes how the Cen X-3 inference changes under a plausible wind turbulence of 50–300 km/s. Since high-mass X-ray binary winds are supersonic and turbulent, zero turbulence is likely an underestimate. If the curve shifts by a factor of 2–3 in NH, the diagnostic remains only order-of-magnitude without an independent velocity or turbulence constraint. The cross-check against the Heα z/w ratio agreeing 'within a factor of a few' is consistent with this uncertainty and does not close it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new column-density diagnostic for X-ray plasmas based on the intensity ratio of the fine-structure components Lyα1 and Lyα2 of the Lyman-α doublet of H-like ions, which can be resolved for the first time with XRISM's microcalorimeter. Cloudy simulations predict that the ratio is about 2 at low column density and approaches 1 at high column density as resonance-line trapping increases. The authors demonstrate the method on a 48 ks XRISM Resolve eclipse observation of Cen X-3, where the Fe XXVI Lyα doublet is resolved and fitted with a ratio of 1.35 ± 0.11. A customized plane-parallel photoionized Cloudy model, assuming log ξ ≈ 4, n = 10^9 cm^-3, and zero turbulence, yields NH ≈ 2 × 10^22 cm^-2, which is said to agree with the Heα z/w method within a factor of a few. The paper argues that the diagnostic is useful for a wide range of collisionally and radiatively ionized plasmas.","tokens_in":6085,"tokens_out":5976,"duration_ms":59402,"significance":"If the diagnostic is robust, it is a genuinely new and simple tool that leverages a capability unique to microcalorimeters, and it would be applicable to many X-ray binaries, active galactic nuclei, and clusters. The paper has concrete strengths: the XRISM spectral fit is clean with reduced chi-squared below 1.2, the Cloudy models and figures are publicly available, and the inferred column density is checked against an independent line-ratio method. The physical explanation via line optical depth and escape probability is clear. However, the quantitative demonstration is only as strong as the assumed plasma model, and the paper does not yet quantify how the inferred NH depends on the most uncertain model inputs, particularly the turbulent velocity.","major_comments":[{"comment":"The Cen X-3 column-density inference is made with the turbulent velocity fixed to zero, yet no test of this assumption is presented. In the line-transfer regime relevant here, the line-center optical depth scales inversely with the Doppler b parameter, so for a fixed ionic column a larger microturbulent velocity reduces photon trapping and shifts the Lyα1/Lyα2-versus-NH curve toward higher NH. Since the stellar wind in a high-mass X-ray binary is expected to be turbulent, the authors should quantify how the inferred NH changes for, e.g., v_turb = 50, 150, and 300 km/s. Without this, the stated NH ≈ 2 × 10^22 cm^-2 has an unquantified systematic uncertainty that could be as large as the claimed factor-of-a-few consistency with the Heα z/w method.","section":"Section 3, Figure 4"},{"comment":"The model parameters used for the Cen X-3 Cloudy simulation (photon index -1.8, LX = 10^37 erg/s, rin = 10^12 cm, log ξ ≈ 4, n = 10^9 cm^-3, plane-parallel geometry) are adopted without uncertainties, and the reported NH carries only the statistical uncertainty from the observed ratio 1.35 ± 0.11. The diagnostic is model-dependent, as the Figure 2 caption itself acknowledges, so the paper should provide a systematic error budget or an explicit sensitivity study over these parameters. The agreement 'within a factor of a few' with the Heα z/w ratio is a useful sanity check but does not substitute for propagating the model uncertainty.","section":"Section 3"},{"comment":"The paper presents Figure 2 as a general diagnostic for selected elements, but the curves are computed for a specific Perseus-like collisionally ionized plasma with fixed temperature, density, and 150 km/s turbulence, and the caption explicitly warns that the absolute NH scale changes with environment. The abstract and Section 4 claim usefulness for a wide range of plasmas, but the manuscript does not show how the calibration changes with temperature, density, abundances, or geometry. To support the general claim, the authors should either provide a small parameter study or explicitly restrict the claim to a model-dependent mapping that must be recomputed for each source.","section":"Section 2, Figure 2; Section 4"}],"minor_comments":[{"comment":"The code name appears as 'C loudy' with an inserted space in several places (e.g., the abstract and Section 2), which appears to be a formatting artifact; please correct.","section":"Abstract, Section 2"},{"comment":"The photon index is written as a negative value ('-1.8'); in standard X-ray spectroscopy the photon index is usually quoted as a positive number (Γ = 1.8), so please clarify the sign convention or correct the value.","section":"Section 3"},{"comment":"The figure caption states that the M1 transition intensity is added to Lyα2, but the text in Section 3 does not mention this when describing the two-Gaussian fit; please state in the text whether the fitted Lyα2 component includes the unresolved M1 contribution and how this affects the comparison with the Cloudy ratio.","section":"Figure 4 caption"},{"comment":"The acceptance date is given as 'Accepted January 18, 2024', which precedes the received date and appears to be a typo for 2025; please correct.","section":"Manuscript header"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript is within scope for A&A Letters and the central idea is promising. My recommendation of major revision is driven by the unquantified systematic uncertainty in the turbulence and model-parameter choices, not by any concern about circularity. The Cloudy calibration is independent of the observed Cen X-3 ratio, so the circularity concern in the stress-test note does not land. If the authors add a sensitivity analysis for v_turb and other model parameters, I would be happy to see the revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The punchline: this is a short, honest paper that proposes a genuinely new use of the resolved Lyα1/Lyα2 doublet ratio as a column density diagnostic, and demonstrates it on a XRISM eclipse observation of Cen X-3. The underlying physics is not new—line ratios going from 2:1 to 1:1 with optical depth is textbook—but the application is: XRISM is the first instrument that can resolve this doublet beyond the Sun, and the Cen X-3 result is the first time it's been used as a column gauge. The fit is clean (reduced χ² < 1.2), the Cloudy models and figures are public, and the paper is transparent that the calibration is environment-dependent (Figure 2 caption). That's real credit.\n\nThe soft spots are the ones the reader flagged. The central quantitative claim—NH ≈ 2×10²² cm⁻² for Cen X-3—comes from a single Cloudy model with zero microturbulence, and no uncertainty is propagated from either the spectral fit or the model assumptions. For a high-mass X-ray binary, a turbulent velocity of a few hundred km/s is plausible, and since the line optical depth goes as the inverse Doppler b, that can shift the inferred NH by a factor of a few. The paper does not quantify this, nor does it argue that zero turbulence is the right choice. The cross-check against the Heα z/w ratio agreeing 'within a factor of a few' is consistent but does not close the issue, because that method is vulnerable to the same environmental assumptions.\n\nI don't think this is a fatal flaw. The paper's claim is that the ratio is a useful column density indicator, not that it gives precise NH values in all environments. For a two-page letter the analysis is proportionate, and the missing systematics are a natural follow-up rather than an error that invalidates the demonstration. Still, a referee should ask for an explicit statement of how the inferred NH changes under plausible variations in turbulence, geometry, and ionization parameter, and for an error budget on the quoted column density.\n\nBottom line: X-ray spectroscopists working with microcalorimeter data and Cloudy users will want to read this. It deserves a serious referee. I'd send it out, and ask for the caveat on systematics to be added before acceptance.","headline":"A solid short letter proposing a useful but environment-sensitive Lyα doublet column density diagnostic, with a clean Cen X-3 demonstration that would benefit from an explicit turbulence/error analysis.","tokens_in":6633,"tokens_out":3373,"would_cite":true,"duration_ms":30943,"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 Ly-alpha1/Ly-alpha2 ratio of hydrogen-like ions changes from about 2 to about 1 with column density, making it a practical X-ray plasma column density indicator.","keywords":["X-ray spectroscopy","column density","line intensity ratios","Lyman-alpha doublet","hydrogen-like ions","radiative transfer","X-ray binaries","microcalorimeters"],"falsifier":"In a source with an independently known column density, for example from photoelectric absorption edges, a measured $\\mathrm{Ly}\\alpha_1/\\mathrm{Ly}\\alpha_2$ ratio that is inconsistent with the Cloudy calibration curve for that element and plasma state would falsify the claim; so would ratios from different H-like ions in one spectrum that imply incompatible $N_{\\mathrm{H}}$ values.","tokens_in":5578,"feed_emoji":"🔭","tokens_out":18530,"duration_ms":134748,"temperature":0.7,"pith_summary":"With the spectral resolution of X-ray microcalorimeters, the two fine-structure components of the Lyman-$\\alpha$ line of hydrogen-like ions, $\\mathrm{Ly}\\alpha_1$ and $\\mathrm{Ly}\\alpha_2$, are now resolvable in astrophysical sources. The paper argues that their intensity ratio is a plasma column density indicator: it is about 2 when the plasma is optically thin and approaches 1 as the line becomes optically thick. This is demonstrated using the spectral simulation code Cloudy and applied to an XRISM observation of the high-mass X-ray binary Centaurus X-3, where the Fe XXVI doublet ratio of $1.35\\pm0.11$ implies a hydrogen-equivalent column density of roughly $2\\times10^{22}\\,\\mathrm{cm}^{-2}$. If the claim holds, any X-ray plasma with a resolvable H-like Lyman-$\\alpha$ doublet can have its column density estimated from a single line ratio.","feed_headline":"2-to-1 Ly-alpha ratio measures plasma column density","feed_subtitle":"Resolved by XRISM, the Ly-alpha doublet ratio drops from 2 to 1 as column thickness grows, tested on Cen X-3.","key_machinery":"The central object is the fine-structure doublet of the Lyman-$\\alpha$ line of one-electron ions: the $2p\\,{}^2P_{3/2}\\to1s\\,{}^2S_{1/2}$ transition ($\\mathrm{Ly}\\alpha_1$) and the $2p\\,{}^2P_{1/2}\\to1s\\,{}^2S_{1/2}$ transition ($\\mathrm{Ly}\\alpha_2$). The mechanism is line radiative transfer. At low optical depth the line photons escape freely and the intensities follow the statistical weights of the upper levels, producing a ratio near 2; at high optical depth photons are repeatedly absorbed and re-emitted until they escape from near the $\\tau=2/3$ surface, driving the ratio toward 1. Cloudy computes this transfer with an escape probability formalism and partial redistribution for resonance lines, and converts the $\\tau=2/3$ location into an $N_{\\mathrm{H}}$ value using the modeled abundance, charge, and level populations.","core_discovery":"The paper's central claim is that the ratio $\\mathrm{Ly}\\alpha_1/\\mathrm{Ly}\\alpha_2$ of hydrogen-like ions changes monotonically with the hydrogen-equivalent column density $N_{\\mathrm{H}}$, from about 2 at low column density, where the upper fine-structure levels are populated in proportion to their statistical weights, to about 1 at high column density, where the emergent flux is determined near the depth at which the line optical depth reaches $2/3$. Cloudy simulations provide the calibration curve for a given plasma environment. The authors demonstrate the diagnostic with the Fe XXVI doublet in an XRISM eclipse observation of Centaurus X-3, obtaining $N_{\\mathrm{H}}\\sim2\\times10^{22}\\,\\mathrm{cm}^{-2}$, consistent with an independent He$\\alpha$ z/w ratio estimate within a factor of a few.","pith_inferences":["The paper leaves implicit that the ratio could be inverted to constrain turbulence or geometry: in systems where $N_{\\mathrm{H}}$ is known from continuum absorption or other line ratios, the observed $\\mathrm{Ly}\\alpha_1/\\mathrm{Ly}\\alpha_2$ would measure the Doppler broadening or the geometrical arrangement of the emitting plasma.","A multi-element extension is natural: if the Cloudy curves are correct, ratios from several H-like ions in one spectrum must converge on a single $N_{\\mathrm{H}}$, giving an internal consistency test that the single-line version does not provide.","Observing Cen X-3 across orbital phases or at successive eclipses would test the predicted ratio-$N_{\\mathrm{H}}$ relation directly, since the absorbing column along the line of sight changes with the binary geometry."],"forward_implications":["The Ly$\\alpha$ doublet ratio becomes a column density diagnostic for any X-ray plasma with a resolvable H-like Ly$\\alpha$ line, whether the plasma is collisionally or radiatively ionized.","Different elements are sensitive to different column density ranges, roughly $10^{20}$ to $10^{24}\\,\\mathrm{cm}^{-2}$, so ratios from several H-like ions can probe different depths in the same plasma.","X-ray microcalorimeters such as XRISM resolve the Ly$\\alpha$ doublet for $Z\\ge18$ (argon) and the Ly$\\beta$ doublet for $Z\\ge26$ (iron), which grating spectrometers could not do.","In the Cen X-3 demonstration, the observed ratio of $1.35\\pm0.11$ maps through the Cloudy model to $N_{\\mathrm{H}}\\sim2\\times10^{22}\\,\\mathrm{cm}^{-2}$, in line with an independent He$\\alpha$ z/w estimate."],"supporting_citations":[{"why":"Provides the Cloudy one-electron microphysics and line-profile treatment that make the doublet resolvable in simulations at microcalorimeter resolution.","marker":"Gunasekera et al. (2024)"},{"why":"Documents the Cloudy code whose atomic and radiative transfer framework is used for all ratio predictions.","marker":"Chatzikos et al. (2023)"},{"why":"Sets out the four line-transfer regimes used to choose escape probability and redistribution behavior in Cloudy.","marker":"Hummer (1962)"},{"why":"Supplies the Eddington-Barbier interpretation that explains why the ratio saturates at 1:1 for optically thick lines.","marker":"Rutten (2003)"},{"why":"Provides the solar Fe XXVI measurement of the optically thin doublet ratio near 2, anchoring the low-column limit.","marker":"Tanaka (1986)"},{"why":"Supplies the APED line-energy data used to fix the Ly alpha line centers in the Cen X-3 spectral fit.","marker":"Smith et al. (2001)"},{"why":"Gives the independent He alpha z/w column density estimate for Cen X-3 used to check the new diagnostic.","marker":"Chakraborty et al. (2021)"},{"why":"Reports the XRISM Resolve observation of Cen X-3 whose eclipse spectra contain the resolved Fe XXVI doublet.","marker":"Mochizuki et al. (2024)"}],"fun_headline_variants":["Ly-alpha doublet ratio reveals plasma column density","XRISM sees Ly-alpha ratio drop from 2 to 1 with column depth","New column density probe: Ly-alpha1 to Ly-alpha2 ratio","Cloudy-calibrated Ly-alpha doublet ratio measures NH in plasmas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calibration is environment-dependent: the same observed ratio maps to different column densities if the plasma's temperature, density, turbulence, abundances, or geometry differs from the assumed Cloudy model.","fun_headline_variants_meta":{"raw":{"variants":["Ly-alpha doublet ratio reveals plasma column density","XRISM sees Ly-alpha ratio drop from 2 to 1 with column depth","New column density probe: Ly-alpha1 to Ly-alpha2 ratio","Cloudy-calibrated Ly-alpha doublet ratio measures NH in plasmas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1238,"prompt_tokens":866,"completion_tokens":372,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":482,"completion_tokens_details":{"reasoning_tokens":296}},"tokens_in":482,"tokens_out":372,"duration_ms":3599,"temperature":1.0,"reasoning_tokens":296,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:23:26.386453+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a source with an independently known column density, for example from photoelectric absorption edges, a measured $\\mathrm{Ly}\\alpha_1/\\mathrm{Ly}\\alpha_2$ ratio that is inconsistent with the Cloudy calibration curve for that element and plasma state would falsify the claim; so would ratios from different H-like ions in one spectrum that imply incompatible $N_{\\mathrm{H}}$ values.","supporting_citations":[{"cited_title":"Cloudy and the High-Resolution Microcalorimeter Revolution: Optical, UV, and X-ray Spectra of One-electron Systems","cited_arxiv_id":"2412.01606","evidence_quote":"Provides the Cloudy one-electron microphysics and line-profile treatment that make the doublet resolvable in simulations at microcalorimeter resolution."},{"cited_title":"2023, Rev","cited_arxiv_id":null,"evidence_quote":"Documents the Cloudy code whose atomic and radiative transfer framework is used for all ratio predictions."},{"cited_title":"J., Chatzikos, M., Guzmán, F., & Su, Y","cited_arxiv_id":null,"evidence_quote":"Gives the independent He alpha z/w column density estimate for Cen X-3 used to check the new diagnostic."},{"cited_title":"L., et al","cited_arxiv_id":null,"evidence_quote":"Reports the XRISM Resolve observation of Cen X-3 whose eclipse spectra contain the resolved Fe XXVI doublet."}],"review_version":1}