{"id":"b6e65edf-1fc5-45a0-9dd3-da5c01db1153","arxiv_id":"2501.02056","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Hydrogen absorption limits the cool circumgalactic gas in dwarf galaxies to 5e8-2e9 solar masses, below 10% of the halo baryon budget.","lead":"Dwarf galaxies' cool gas halos are measured to contain at most 5e8 to 2e9 solar masses, under 10% of each halo's expected baryon budget. The rest of the baryons must be hidden in hotter gas or ejected into intergalactic space.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported 'upper limits' are not formal upper limits: the highest NHI values are lower limits treated as outliers, and the by-eye fits assume true columns near the lower bounds; the paper's own sensitivity test allows a factor ~1.4-2 higher mass.","rationale":"I read the paper as making a carefully bounded claim: for volume-filling, photoionized cool gas, a power-law density model fitted to HI absorption yields an upper limit on the cool CGM mass in dwarf halos. The analytic argument in Sec. 2.2 and Appendix B is sound for a specified column-density profile: given the observed NHI profile, pointwise minimization of density at fixed emissivity maximizes mass, and that minimum is achieved by fV = 1. The Cloudy treatment, the redshift and metallicity uncertainty estimates, and the public data table all support reproducibility and lend credibility to the direction of the result. My concern is not with the analytic upper-limit theorem but with the empirical step that turns censored observations into a headline number. The sample contains numerous lower limits on NHI, including some very high values, and these are handled by fitting curves that are 'close to many lower limits' and by discounting the highest ones. Since for fV = 1 the derived mass is an increasing function of the assumed column, treating a lower limit as if the true column were near the limit produces a mass that is better described as a lower envelope of the mass needed to explain the minimum observed columns, not a formal upper limit on the mass that could be present. The paper's own sensitivity test shows that this choice changes the result by factors of 1.3-2.0, which is larger than the claimed metallicity and redshift uncertainties and large enough to alter the '10% of baryons' and '15% total baryons' statements. This does not overturn the qualitative conclusion that dwarf CGM cool gas is a small baryon reservoir, but it does mean the headline upper limit should be presented with a censored-data treatment or clearly relabeled as an upper limit for the representative/average sightline population. The reader's verdict is already CONDITIONAL and explicitly asks for a transparent treatment of the high lower limits, so my read does not change the verdict; it sharpens the condition that should be met before the numerical upper limits are quoted as formal bounds.","tokens_in":20820,"tokens_out":14584,"duration_ms":172308,"concrete_test":"Recompute the three mass-bin results with a formal censored likelihood: detections enter as log-normal measurements, upper limits as survival probabilities P(N_model < N_UL), and lower limits as P(N_model > N_LL). Also compute a strict upper-limit variant that maximizes M_cCGM over (a_n, n_H,0) subject to N_model being at or above every lower-limit value and at or below every upper-limit value. Compare the resulting 95% upper bounds on M_cCGM to Table 1; if any bin's bound exceeds the quoted value by more than a factor of 1.5, the headline 'upper limit' is not robust to the adopted lower-limit treatment.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantity is an upper limit on M_cCGM, but the procedure in Sec. 2.4 does not actually construct one. The binned sample contains detections, upper limits, and lower limits; model curves are chosen to be 'consistent with most measurements' and 'close to many lower limits', while the few high lower limits are set aside as unrepresentative. Because the fV = 1 inferred mass increases monotonically with the assumed NHI, a genuine upper limit requires maximizing the mass over all column values allowed by the data, or at least a well-defined censored-data fit. Instead, lower limits are effectively treated as point values near their bounds. Sec. 2.5 quantifies the consequence: replacing lower limits by 3x (5x) larger values raises M_cCGM by factors of 1.3-1.6 (1.5-2.0) per bin, so the headline range 5e8-2e9 M_sun and the '<10% of baryons' statement can shift upward by a factor of two, potentially exceeding the abstract's 15% total baryon budget once stars are included. Thus the headline result is contingent on an outlier-rejection and censoring convention that is not formalized, rather than a model-independent upper limit. The analytic fV = 1 proof for a specified column profile is not the weak point; the weak point is applying that proof to censored, heterogeneous sightlines without a statistical treatment of the censoring.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a phenomenological model of cool (T≈10^4 K), photoionized circumgalactic gas in dwarf galaxies and applies it to archival HI absorption measurements from the Z24 and M24 samples. The model assumes a single power-law density profile n_H ∝ r^{-a_n} between 0.1 R200m and R200m, a constant volume filling fraction fV, and gas in photoionization and heating equilibrium with a uniform metagalactic UV background. The authors show analytically that, for a fixed HI column density profile, the required cool-gas mass scales as M ∝ N_HI^{1/2} fV^{1/2}, so volume-filling (fV=1) gas maximizes the inferred mass. Fitting the model to HI columns binned into three halo-mass bins yields M_cCGM ≈ 5×10^8–2×10^9 M_sun for fV=1, corresponding to ≈5–10% of the cosmological halo baryon budget; clumpy models with fV=0.01 give masses ≈11 times lower. The paper also estimates metallicity and redshift uncertainties, discusses gas densities, pressures, baryon and metal budgets, accretion rates, and compares with more massive halos.","tokens_in":21219,"tokens_out":6374,"duration_ms":60787,"significance":"If the upper-limit construction were made formal, this would be a valuable first empirical constraint on the cool CGM mass in 10^10–3×10^11 M_sun halos, directly relevant to the missing-baryon problem and to feedback models in dwarf galaxies. The analytic scaling argument in §2.2 and Appendix B is clean, and the conservative direction of assigning all measured HI to the cool, volume-filling phase is physically sensible. The paper is also transparent about many model limitations and provides a useful comparison with previous work on more massive halos. However, the headline 'upper limit' is not actually constructed as an upper limit in the data-fitting step: the treatment of lower limits and outliers in §2.4–2.5 is a censoring convention, not a maximization. The significance of the paper therefore hinges on whether that gap can be closed by a formal statistical treatment or by reframing the claims.","major_comments":[{"comment":"The central quantity is reported as an upper limit on M_cCGM, but the fitting procedure in §2.4 does not construct one: the models are chosen to be 'consistent with most measurements' and 'close to many of the lower limits', while the few high lower limits are set aside as 'may not be representative'. Because M_cCGM increases monotonically with the assumed N_HI (Eq. B4/B5), a genuine upper limit requires either maximizing the mass over all column values allowed by the data, including lower limits and outliers, or performing a formal censored-data fit. The paper's own sensitivity test in §2.5 shows that replacing lower-limit columns by factors of 3 and 5 raises M_cCGM by factors of 1.3–1.6 and 1.5–2.0 per bin, respectively; thus the headline range 5×10^8–2×10^9 M_sun and the '<10% of baryons' statement can shift upward by roughly a factor of two, potentially exceeding the abstract's '<15% total baryon budget' once stars are included. The claims need to be reframed as conditional on the censoring/outlier convention, or the upper-limit construction needs to be made formal.","section":"§2.4, §2.5"},{"comment":"No goodness-of-fit or statistical uncertainty is reported for the manual grid search over a_n and n_H,0. Table 1 lists single parameter values for each bin, but there is no measure of how many models are acceptable, what the scatter among data points contributes, or how the inferred M_cCGM varies within the acceptable region of the grid. Reporting a chi-square-like statistic or a likelihood surface over the grid, and propagating the resulting parameter ranges to M_cCGM, would make the comparison with the data quantifiable and would allow the reader to assess whether the quoted factor-of-two sensitivity is dominated by the censoring or by the model fit itself.","section":"§2.4, Table 1"},{"comment":"The redshift treatment is not conservative for an upper limit. The nominal models adopt the z=0 MGRF, while the sample has median z=0.1 and 90% of galaxies at z<0.2; §2.5 states that the required masses are about 30% higher at z=0.1 and about 60% higher at z=0.2. Since a higher MGRF intensity implies more ionization and hence a larger inferred gas mass for fixed N_HI, modeling all sightlines at z=0 pushes the inferred masses systematically low. The abstract's statement that the redshift-range uncertainty is ≈15% is based on a symmetric 0≤z≤0.1 interval (the z=0.05 case), not on the actual sample redshift distribution. For an upper limit, one should either evaluate each sightline at its own redshift or adopt a high-redshift MGRF for the limit, and the reported uncertainty should reflect the sample distribution.","section":"§2.5, §2.1, Abstract"}],"minor_comments":[{"comment":"In the sentence 'producing a column of NHI ∼ 10^14 cm^-3', the units should be cm^-2, not cm^-3.","section":"§3.5"},{"comment":"The statement 'the mass uncertainty corresponding to 0 ≤ z ≤ 0.1 is ±15%' is inconsistent with the earlier statement that z=0.1 gives a +30% mass increase; clarify that the ±15% is the symmetric uncertainty around z=0.05, not around the nominal z=0 model.","section":"§2.5"},{"comment":"The dotted curves are described both as power-law fits to guide the eye and as clumpy fV=0.01 models; distinguish these two uses clearly in the caption or in the text.","section":"Figure 2 caption"},{"comment":"State explicitly that R200m and M200m are median values of each halo-mass bin, since the binning procedure uses medians rather than the full halo-mass distributions.","section":"Table 1"},{"comment":"The phrase 'taking L = 2 R_CGM (a good approximation out to ≈0.5 R_CGM)' is unclear; the path length through a sphere at impact parameter b is 2(R_CGM^2 - b^2)^1/2, so specify the approximation and its range of validity more precisely.","section":"Appendix B, Eq. (B5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the analytic derivation is sound; the main issue is the gap between the 'upper limit' language and the actual censoring/outlier procedure in §2.4–2.5. I would support publication after the upper-limit construction is formalized or the claims are reframed accordingly. The authors' reliance on shared-author papers (Z24, F24) for data and model choices is disclosed and does not by itself constitute a problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first attempt to map cool CGM mass against halo mass for dwarfs, and the core argument—that volume-filling gas maximizes the mass needed to reproduce a given HI column—is clean and correct. The result, M_cCGM ≲ 5e8–2e9 M_sun (5–10% of the halo baryon budget), is a genuinely new empirical constraint and will be cited.\n\nWhat it does well: the analytic derivation in §2.2/Appendix B is transparent and the f_V^{1/2} scaling is derived, not assumed. The paper is honest about its uncertainties: metallicity (±10%), redshift (30–60% at z=0.1–0.2), and it runs a sensitivity test replacing lower limits with 3× and 5× higher values. The comparison to Z24 and more massive halos is balanced. Data table is included, which helps reproducibility.\n\nThe soft spots are real but not fatal. The headline 'upper limits' are not formal: the fits are by-eye, with no goodness-of-fit or statistical treatment of censored data. Lower limits are effectively treated as point values near their bounds, and the few high lower limits are set aside as unrepresentative. The paper's own sensitivity test shows that raising lower limits by 3× increases the mass by factors of 1.3–1.6 per bin, and 5× by 1.5–2.0. That means the '≲10% of baryons' statement could become ~20% if the high columns are real. The authors acknowledge this, but the presentation still leans on the optimistic reading. A proper censored-data fit—or at least a conservative fit that maximizes over allowed columns—would make the headline claim rigorous. Also, the binning uses median halo masses without propagating scatter, which is minor.\n\nThe stress-test note is on target. The analytic proof for a specified column profile is fine; the weakness is applying it to heterogeneous, censored sightlines. I don't think this undermines the direction of the result—dwarfs clearly do not have most of their baryons in the cool CGM—but the numerical upper limits carry an unquantified systematic that a reviewer should push on.\n\nBottom line: this deserves a serious referee. It's a solid, useful paper that needs moderate revision (quantitative fitting, clearer censoring treatment, possibly a more conservative headline). I'd cite it. Reading group: maybe—it's a good example of how to present upper limits with acknowledged caveats, and the discussion of baryon budget is thought-provoking.","headline":"A useful first empirical constraint on cool CGM mass in dwarf halos, with an honest but informal handling of censored data that should be tightened in review.","tokens_in":21670,"tokens_out":1565,"would_cite":true,"duration_ms":15379,"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":"By modeling HI absorption around 40 dwarf galaxies, this paper proves volume-filling gas maximizes the inferred cool-gas mass and places upper limits of $5\\times10^8$–$2\\times10^9\\,M_\\odot$ on that mass.","keywords":["circumgalactic medium","dwarf galaxies","neutral hydrogen absorption","cool gas mass","photoionization equilibrium","baryon budget","volume filling fraction","HI column density"],"falsifier":"A measurement of $N_{\\rm HI}>10^{16}\\,{\\rm cm}^{-2}$ at $r_\\perp/R_{200m}>0.9$ in a dwarf halo would exceed what the volume-filling power-law model can produce, forcing either a steeper density profile, higher clumpiness, or a substantial contribution from gas phases outside the cool CGM.","tokens_in":20667,"feed_emoji":"🌌","tokens_out":6815,"duration_ms":66279,"temperature":0.7,"pith_summary":"The paper asks how much cool ($\\approx10^4$ K), photoionized gas surrounds dwarf galaxies, a phase that is hard to observe and has been poorly measured. Using archival HI absorption measurements toward 40 dwarf galaxies, the authors model the gas as a single power-law density profile in photoionization equilibrium with the metagalactic UV background. They show analytically that, for a fixed observed HI column, smooth volume-filling gas requires the largest total mass, so their mass estimates are conservative upper limits. For halo masses $M_{200m}=10^{10}$–$3\\times10^{11}\\,M_\\odot$ they infer $M_{\\rm cCGM}=5\\times10^8$–$2\\times10^9\\,M_\\odot$, less than 10% of each halo's baryon budget, and they argue that stars plus cool gas account for less than 15% of the baryons in dwarfs. This matters because it limits the fuel available for star formation and sharpens the question of where the missing baryons reside.","feed_headline":"Dwarf halos cap cool gas below 10% of baryons","feed_subtitle":"HI absorption limits the cool circumgalactic medium of dwarfs to 5e8–2e9 solar masses.","key_machinery":"The central object is a phenomenological density profile $n_H(r)=n_{H,0}(r/R_{\\rm CGM})^{-a_n}$ for gas between 0.1 and 1 $R_{200m}$, with a constant volume filling fraction $f_V$ and the gas in heating, cooling, and ionization equilibrium with the metagalactic UV background. The argument is carried by the analytic identity $M_{\\rm cCGM}\\propto N_{\\rm HI}^{1/2}f_V^{1/2}$, derived from $N_{\\rm HI}\\approx10\\,L\\,n_H^2 f_V$: at the low densities of this gas the neutral fraction is roughly proportional to density, so the column depends quadratically on density, and since $f_V\\le1$ the mass is maximized when the gas fills the volume. A standard photoionization code supplies the equilibrium temperature and neutral fractions, and the model converts observed HI columns into a spherical cool-gas mass and its scaling with clumpiness.","core_discovery":"For a given HI column density, the mass of cool photoionized gas needed to reproduce it grows with the volume filling fraction as $M_{\\rm cCGM}\\propto f_V^{1/2}$, so $f_V=1$ gives an upper limit on the mass. The authors fit power-law density profiles $n_H(r)=n_{H,0}(r/R_{200m})^{-a_n}$ with constant $f_V$ to observed HI columns in three halo-mass bins and find $M_{\\rm cCGM}=5\\times10^8$–$2\\times10^9\\,M_\\odot$, corresponding to 5–10% of the cosmological baryon budget. Clumpy models with $f_V=0.01$ reproduce the same columns with about 11 times less mass, matching the analytic scaling. The paper concludes that dwarf galaxies hold at most $\\lesssim15\\%$ of their baryons in stars and cool CGM, with the remainder in a warm/hot phase or ejected into the intergalactic medium.","pith_inferences":["Inference: Because the mass scales only as $f_V^{1/2}$, even an order-of-magnitude error in clumpiness changes the inferred mass by only a factor of about three, so the conclusion that cool gas is a minority of the baryon budget is robust unless the filling fraction is extremely small.","Inference: The paper's upper-limit framework could be sharpened by measuring metal-line absorption in the same sightlines; a metallicity estimate would fix the equilibrium temperature and break the remaining degeneracy between density and $f_V$.","Inference: If future observations stack O VI or X-ray data around dwarf halos, they could distinguish whether the missing ~85% of baryons remain as warm/hot CGM or have been ejected into the IGM, a test the paper does not perform.","Inference: The analytic relation suggests that HI column profiles alone cannot uniquely determine cool-gas mass, so progress on this question will likely require kinematic information from resolved line profiles to estimate the gas distribution along the line of sight."],"forward_implications":["For the three halo-mass bins, the inferred cool CGM mass rises with halo mass from about $5\\times10^8$ to $2\\times10^9\\,M_\\odot$ while staying below 10% of each halo's cosmological baryon budget.","Stars plus cool CGM together make up less than 15% of the baryon budget in dwarf halos, so the majority of baryons must be warm/hot or already ejected.","If the cool gas is clumpy with $f_V=0.01$, the same HI columns require roughly 11 times less mass, so clumpiness cannot hide a large cool reservoir unless $f_V$ is far smaller.","Assuming infall on a dynamical timescale, the cool CGM can supply accretion rates at or above the current star formation rates for about a gigayear.","A constant cool-gas baryon fraction of a few to ten percent may extend from dwarf to $L_*$ halos, though survey and method differences still prevent a firm trend."],"supporting_citations":[{"why":"Supplies the base archival sample of dwarf galaxies with HI absorption measurements, stellar masses, and the initial cool-gas mass estimate that this work extends.","marker":"Z24"},{"why":"Adds the CUBS survey HI absorption measurements for low-mass galaxies, extending the sample and enabling halo-mass binning.","marker":"M24"},{"why":"Provides the equilibrium temperature relation $T_{\\rm eq}(n_H,Z)$ used to set gas temperature and ionization state in the models.","marker":"F24"},{"why":"Supplies the metagalactic UV background spectrum at $z=0$ used in the photoionization calculations.","marker":"Khaire & Srianand (2019)"},{"why":"The photoionization code used to compute neutral hydrogen fractions and equilibrium temperatures.","marker":"Ferland et al. 2017"},{"why":"Converts stellar masses to halo masses for the sample, setting the $M_{200m}$ bins.","marker":"Behroozi et al. 2019"},{"why":"Provides the cosmological conversion to $M_{200m}$ and $R_{200m}$ used to normalize impact parameters.","marker":"Diemer 2018"},{"why":"Supplies the cosmological baryon fraction used to compute the halo baryon budget and the reported percentages.","marker":"Planck Collaboration et al. 2016"}],"fun_headline_variants":["Dwarf galaxies cap cool gas at 10% of baryon budget","Cool gas in dwarf halos: maximum 2e9 solar masses","Volume-filling gas sets upper limit on dwarf cool mass","Dwarf baryons: at most 15% in stars and cool CGM","Clumpy gas would slash dwarf cool mass by 11x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the observed HI column densities come entirely from cool gas that fills space with a constant volume filling fraction and follows a single power-law density profile in photoionization and thermal equilibrium with a uniform UV background; if the gas is clumpy, multi-phase, or out of equilibrium, the quoted mass limits shift.","fun_headline_variants_meta":{"raw":{"variants":["Dwarf galaxies cap cool gas at 10% of baryon budget","Cool gas in dwarf halos: maximum 2e9 solar masses","Volume-filling gas sets upper limit on dwarf cool mass","Dwarf baryons: at most 15% in stars and cool CGM","Clumpy gas would slash dwarf cool mass by 11x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001677,"raw_usage":{"total_tokens":6730,"prompt_tokens":1105,"completion_tokens":5625,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":721,"completion_tokens_details":{"reasoning_tokens":5532}},"tokens_in":721,"tokens_out":5625,"duration_ms":42887,"temperature":1.0,"reasoning_tokens":5532,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:14:51.325877+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of $N_{\\rm HI}>10^{16}\\,{\\rm cm}^{-2}$ at $r_\\perp/R_{200m}>0.9$ in a dwarf halo would exceed what the volume-filling power-law model can produce, forcing either a steeper density profile, higher clumpiness, or a substantial contribution from gas phases outside the cool CGM.","supporting_citations":[],"review_version":1}