{"id":"85178158-3f0f-4aa6-91c6-17eaf2c21e9c","arxiv_id":"2507.19957","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"VLBA images of 12 FLASH HI absorbers at z=0.44-0.68 show that the radio core can account for the full absorption in 7 sources, raising the estimated HI optical depths by up to a factor of three after covering-factor corrections.","lead":"Astronomers used the VLBA network to image 12 distant galaxies where neutral hydrogen absorbs radio light, and found that for 7 of them the compact core alone can explain the full absorption. The work provides the first parsec-scale maps for these intermediate-redshift HI absorbers and suggests their true gas column densities may be up to three times higher than ASKAP measured.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Frequency mismatch between 1.4 GHz continuum ratios and ~0.7 GHz absorption fractions can change core-sufficiency classifications and the factor-three VOD correction.","rationale":"The reader's weakest_assumption is the frequency mismatch between the 1.4 GHz VLBA/RACS flux ratios and the ~0.7–1.0 GHz ASKAP absorption fractions. I examined the derivation in Appendix 1 and the application in §3.2.2 and agree this is the most load-bearing concern. The inequality in Eq. (9) is exactly Smin/SRACS > PF, where the left side is a 1.4 GHz ratio and the right side is a line-frequency ratio; this requires equal spectral indices for the compact and total emission. The sample's strong preference for peaked or inverted SEDs (§4.1) makes a positive α_c − α_t plausible for many sources, so the bias is not second-order: it can flip the 'core sufficient' classification and materially change the VOD upper limits from which the median factor-3 correction is derived. The paper's own statement that 'a direct comparison cannot be made' is an explicit admission of missing support, which the reviewing rule requires flagging. I considered whether the logical issue that Score/SRACS is only a lower limit on f if the gas covers the entire core is even more fundamental, but the frequency mismatch enters first and contaminates every ratio in the sample, so I treat it as the single most load-bearing concern. The VLBA data and source classifications remain valuable, and the paper is transparent about the caveat, so the reader's CONDITIONAL verdict is appropriate; my read does not change it.","tokens_in":25381,"tokens_out":12960,"duration_ms":170128,"concrete_test":"For the four targets with astrogeo images (0023+010, 0141-231, 0518-245, 2007-245), fit the same components at 1.4 GHz and at 4–8 GHz to measure α_c and α_t; for all 12 targets, use RACS 1.4 GHz and ASKAP 0.7–1.0 GHz total fluxes to derive α_t, and adopt α_c = 0, +0.5, and +1.0 as a sensitivity range. Rescale each core fraction as R(ν21) = R(1.4 GHz) × (ν21/1.4 GHz)^(α_c − α_t), re-run the §3.2.2 classification, and recompute the Table 3 VOD upper limits. If any 'core sufficient' target changes class, or if the median VOD correction ratio changes by more than about 30%, the frequency-mismatch bias is material and the central factor-three claim is not yet established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix 1 (Eq. 9) and §3.2.2 compare Smin/SRACS, measured with VLBA and RACS at ~1.4 GHz, against PF = |SHI|/SASKAP, measured at the redshifted HI line near 0.7–1.0 GHz. This comparison is valid only if the compact-core to total flux-density ratio is the same at both frequencies. Writing R(ν) for the core fraction and S∝ν^α for each component, R(ν21)/R(1.4 GHz) = (ν21/1.4 GHz)^(α_c − α_t). The sample is dominated by peaked-spectrum sources (9/12, §4.1); if the compact core is self-absorbed below turnover (α_c > 0) while extended emission is optically thin (α_t < 0), then α_c − α_t > 0 and the core fraction at 0.7 GHz can be less than half its 1.4 GHz value (e.g., 0.5^1.2 ≈ 0.44). Consequently, the tabulated Smin/SRACS values can overestimate the covering factor at the line frequency, and sources such as 2236-251 (Score/SRACS = 0.33 vs PF = 0.32) or 1136+004 (0.27 vs 0.18) could switch from 'core sufficient' to 'total VLBA flux required'. That changes Smin and the Table 3 VOD upper limits, and the median factor-3 VOD correction is computed from exactly those limits. The paper acknowledges that a direct comparison cannot be made but does not quantify the resulting bias; that omission is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents VLBA 1.4 GHz continuum observations of twelve H I 21-cm absorbers detected in the FLASH pilot surveys at redshifts 0.4 < z < 0.7. The authors classify the parsec-scale source morphologies, identify the radio core in each source, and develop a methodology to test whether the core alone, or the total emission detected in the VLBA image, has sufficient flux density to cause the entire H I absorption measured by ASKAP. Using the ratio of VLBA to RACS flux densities as a lower limit on the gas covering factor, they derive upper limits on the velocity-integrated optical depth (VOD). They report that for seven of twelve sources the core is sufficient to explain the absorption, that three compact sources likely have most of their VLBA-scale emission occulted, and that the distribution of H I VODs at 0.4 < z < 1.0 could increase by up to a factor of three after covering-factor correction. The paper also highlights two high-VOD candidates, 0023+010 and 0903+010.","tokens_in":25702,"tokens_out":7524,"duration_ms":81664,"significance":"If the central claims hold, this work provides a practical method for correcting H I 21-cm absorption measurements for covering factors at intermediate redshifts, where spectroscopic VLBI at the redshifted line frequency is not currently feasible. The methodology is clearly described, the per-source tables are internally consistent, and the arithmetic in Appendix 1 is correct. The paper makes good use of public astrogeo images and is transparent about its assumptions and small sample size. The conclusion that the VOD distribution may be underestimated by a median factor of roughly 2.8 is important for evolutionary studies of cold gas, and the identification of two absorbers with VOD limits above 100 km/s is of significant interest. However, the central quantitative results depend on comparing flux ratios at 1.4 GHz with absorption fractions at ~0.7 GHz, and this frequency mismatch is not quantified. Because the sample is dominated by peaked-spectrum sources, this is a load-bearing issue that affects the 'seven of twelve' claim and the factor-three correction.","major_comments":[{"comment":"The test Smin/SRACS > PF compares the 1.4 GHz VLBA-to-RACS flux ratio with the ~0.7 GHz ASKAP absorption peak fraction. This is valid only if the core-to-total flux ratio is the same at both frequencies. Writing R(ν) for the core fraction and S ∝ ν^α for each component, R(ν21)/R(1.4 GHz) = (ν21/1.4 GHz)^(α_c − α_t). For a self-absorbed core (α_c > 0) with optically thin extended emission (α_t < 0), this ratio can be about 0.4–0.5, so the tabulated Score/SRACS values can overestimate the covering factor at the line frequency. In this sample, nine of twelve sources have peaked SEDs (§4.1), and borderline cases such as 2236-251 (Score/SRACS = 0.33 vs PF = 0.32) and 1136+004 (0.27 vs 0.18) could change classification if the bias is at this level. The paper acknowledges in Appendix 1 that a direct comparison cannot be made, but it does not quantify the resulting systematic error; this is load-bearing for both the 'seven of twelve' claim and the median factor-2.8 VOD correction.","section":"Appendix 1, Eq. (9) and §3.2.2"},{"comment":"The source 2007-245 is included among the seven sources for which 'the core has sufficient flux density' (§4.1), but in §3.3 the authors deliberately adopt Smin = Stot instead of Score because the absorber may be intervening. The methodology in §3.2.2 defines core sufficiency by the condition PF < Score/SRACS, which is satisfied for this source, yet the adopted covering-factor limit is based on Stot. Similarly, 0023+010 is excluded from the seven in §4.1 even though §3.3 states that its core is sufficient based on Score/Stot. The counting rule and the treatment of these two exceptions should be stated explicitly so that the reader can verify the 'seven of twelve' claim.","section":"§3.3 and §4.1"},{"comment":"The VLBA observations were made in February–April 2024, while the RACS-mid and FLASH/ASKAP measurements are from different epochs. For 0023+010 the authors note that Stot exceeds SRACS and interpret this as variability, but no variability check is presented for the remaining sources. Because the inequality tests in §3.2.2 are sensitive to the flux ratios (e.g., 2236-251, where Score/SRACS = 0.33 and PF = 0.32 are nearly equal), variability at the ten to twenty percent level, which is common in compact AGN, could change the classification of individual sources and hence the upper-limit VODs and the median correction factor.","section":"§2 and Table 3"}],"minor_comments":[{"comment":"Table 1 lists source '2236-351' while all other tables and the text use '2236-251'; the coordinates match 2236-251, so this appears to be a typo.","section":"Table 1"},{"comment":"The Figure 3 caption gives zHI = 0.5159 for 0023+010, while Table 2 lists z = 0.6745; one of these is incorrect and should be corrected.","section":"Figure 3 / Table 2"},{"comment":"The text states Smin/Stot = 0.36, but Table 3 reports Smin/SRACS = 0.36; since the authors use Stot as the total flux for this source, the table entry should be labeled consistently with the text.","section":"§3.3, 0023+010"},{"comment":"The text includes 'NHI > 1022 cm–1', but the unit should be cm^-2 for a column density; please correct this typographical error.","section":"§4.3"},{"comment":"The text contains formatting artifacts such as 'Y oon et al.' and 'V ol.' in the references; these should be fixed to 'Yoon et al.' and 'Vol.'.","section":"Throughout"},{"comment":"Since PF is defined as the peak absorption fraction, it would be helpful to restate explicitly that the inequality test in §3.2.2 applies to the line peak, while the covering-factor correction is applied per channel when computing the VOD upper limit.","section":"§3.4"}],"recommendation":"major_revision","confidential_remarks":"This is a solid observational paper from an experienced team, and the underlying data are well presented. The frequency-mismatch issue identified in the major comments is the key obstacle: the authors acknowledge it but do not quantify it. If they can show, using the multi-frequency astrogeo images and the RACS/ASKAP total flux densities already in hand, that the core-to-total ratio does not change dramatically between 1.4 GHz and ~0.7 GHz for this sample, the quantitative conclusions would be robust. Otherwise, the claims should be reframed as conditional on spectral-index assumptions, which would substantially weaken the headline results. The paper is not circular and the arithmetic is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the first parsec-scale continuum imaging of FLASH HI 21-cm absorbers at z~0.5-0.7, and the data are a solid contribution. The per-source morphologies, core identifications, and the Smin covering-factor framework are clearly presented, and the tables are internally consistent. The two high-VOD candidates (0023+010 and 0903+010) are interesting targets for follow-up.\n\nThe soft spot is the frequency mismatch in the covering-factor comparison. Equation 9 compares Smin/SRACS measured at 1.4 GHz with PF = |SHI|/SASKAP measured at the redshifted line near 0.7-1.0 GHz. The paper acknowledges this in Appendix 1 but does not quantify the bias. If the compact core has a rising spectrum below turnover (common in the 9/12 peaked-spectrum sources) while extended emission steepens, the core fraction at 0.7 GHz can be well below the 1.4 GHz value. That would overestimate covering factors and the corrected VOD upper limits. For sources with PF close to Score/SRACS (2236-251: 0.32 vs 0.33; 1136+004: 0.18 vs 0.27), the classification could flip from core-sufficient to total-VLBA-flux required, which changes Smin and the Table 3 limits. Since the median factor-3 VOD increase is computed from exactly those limits, this is a load-bearing caveat. It does not kill the paper; the high-margin cases (e.g., 0141-231, PF=0.09 vs core fraction 0.57) are robust, and 0903+010's argument likely strengthens if the 0.7 GHz core fraction is lower. But the headline number needs a sensitivity analysis using plausible spectral indices or multi-frequency VLBI.\n\nTwo smaller concerns: the 2007-245 exception (using Stot instead of Score because of associated/intervening ambiguity) is ad hoc and should be flagged more explicitly, and the factor-3 statement is based on upper limits, not measured VODs, on a sample of 12 - the paper does caution about this, which is good. Variability is discussed for one source only; ASKAP and RACS are different epochs, so variability could matter more than they say.\n\nOverall, this is a well-executed observational paper with a clear method and new data. It deserves a serious referee. I would ask the authors to quantify the spectral-index bias in Eq. 9, ideally using the multi-frequency astrogeo images they already have, and to present the factor-3 claim with error bars or a range.\n\nRecommendation: send to peer review. It should be published after revisions.","headline":"Useful first VLBI continuum look at FLASH HI absorbers, but the factor-of-three VOD correction rests on an unquantified 1.4-to-0.7 GHz comparison that deserves scrutiny.","tokens_in":26333,"tokens_out":4063,"would_cite":true,"duration_ms":47444,"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":"Parsec-scale radio cores contain enough flux to account for all of the detected H I 21-cm absorption in seven of twelve FLASH absorbers, and covering-factor corrections could raise the population's velocity-integrated optical depths by up…","keywords":["galaxies: active","galaxies: ISM","methods: observational","radio lines: galaxies","radio continuum: general","surveys","H I 21-cm absorption","very long baseline interferometry"],"falsifier":"Direct spectroscopic VLBI observations of the H I 21-cm line at the redshifted frequency (~700-1000 MHz) toward these twelve sources would settle the claim: if the line is detected in absorption against an unresolved component whose flux density is smaller than the adopted core, or is not detected against the core, then the core-based covering factors and the factor-of-three correction are wrong. A cheaper first test is comparing the 1.4 GHz and 5 GHz VLBA core fractions of the same targets, since a systematic difference would show that the assumed core-to-total ratio is frequency dependent.","tokens_in":25165,"feed_emoji":"📡","tokens_out":11401,"duration_ms":116475,"temperature":0.7,"pith_summary":"This paper uses VLBA 1.4 GHz images of twelve neutral-hydrogen (H I) 21-cm absorbers discovered by FLASH, an untargeted ASKAP survey at redshifts 0.4-0.7, to ask whether the absorption is produced against the compact radio core or against more extended emission. It finds that in seven of the twelve sources the parsec-scale core alone has enough flux density to produce the entire absorption seen in the ASKAP spectra, and that for three compact sources with projected sizes of 305-409 pc most of the milliarcsecond-scale emission is probably occulted by gas. Using the measured core and total VLBA flux densities as lower limits on the gas covering factor, the paper derives corrected velocity-integrated optical depths (VODs) that are a median factor of 2.8 higher than the uncorrected ASKAP values, so the VOD distribution at 0.4