REVIEW 3 major objections 3 minor 36 references
FAST H I 21 cm study of blueberry galaxies
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Blueberry galaxies appear nearly devoid of atomic hydrogen: FAST detects H I in only 2 of 28 targets, and the remaining upper limits imply depletion times around 0.5 Gyr.
desk verdict First real H I census of blueberries finds a low detection rate, but the size of the gas-poor offset rests on an unvalidated 50 km/s line-width assumption. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the H I 21 cm line as a tracer of the cold atomic gas reservoir, observed with the Five-hundred-meter Aperture Spherical radio Telescope. The load-bearing step is the conversion of nondetection noise into a $3\sigma$ H I mass upper limit, Equation (6), which assumes a 50 km/s line width for every source; that limit, combined with stellar masses and SFRs from the parent sample, gives the gas fraction $f_{\rm HI}=M_{\rm HI}/M_\ast$ and depletion time $\tau_{\rm HI}=M_{\rm HI}/{\rm SFR}$ plotted against main-sequence relations. The O32 ratio, defined as O[III] $\lambda$5007 / O[II] $\lambda$3727, organizes the sample and connects the low detection rate to the high-O32 Green Pea population, where earlier work also found few H I detections.
What would settle it
A targeted interferometric H I survey of the 26 nondetections reaching $M_{\rm HI}<2\times10^{8}\,M_\odot$ would settle it: recovering emission in more than a few percent of targets, or showing that the assumed 50 km/s line width is too narrow by a factor of several, would overturn the conclusion that blueberries are systematically gas-poor.
Extended reading notes
Core claim
The central claim is that blueberry galaxies at $z\lesssim0.05$ have H I contents well below what main-sequence scaling relations predict for their stellar masses and star-formation rates. Only two of 28 targets (J1026+0426 and J1132+0809) show significant H I emission, and those detections are entangled with neighboring galaxies in FAST's 3-arcmin beam; if the H I is apportioned by stellar-mass ratio, the blueberries themselves hold roughly $5.4\times10^{7}\,M_\odot$ and $1.7\times10^{8}\,M_\odot$, respectively. For the 26 nondetections, the median $3\sigma$ upper limit on H I mass is $2.0\times10^{8}\,M_\odot$, the gas fraction $f_{\rm HI}=M_{\rm HI}/M_\ast$ lies below the $f_{\rm HI}$--$M_\ast$ and $f_{\rm HI}$--sSFR relations for main-sequence galaxies, and the median $3\sigma$ upper limit on H I depletion time $\tau_{\rm HI}=M_{\rm HI}/{\rm SFR}$ is 0.53 Gyr, about an order of magnitude shorter than typical local main-sequence values. The resulting H I detection rate is only 7.1%, and the paper interprets this as evidence that blueberries, like high-O32 Green Peas, exhaust their atomic gas on short timescales and are likely seen in a brief, gas-hungry starburst phase.
Load-bearing premise
The quantitative conclusions hinge on treating the $3\sigma$ upper limits, computed with an assumed 50 km/s line width, as accurate proxies for the true H I content of the nondetections, and on attributing the two detected profiles to the blueberries rather than to neighboring galaxies inside FAST's large beam.
Editorial extensions
If this is right
- Blueberry starbursts cannot be sustained by their present atomic reservoirs for more than about 0.5 Gyr; either the star formation is short-lived or the gas must be accreted or recycled from the surroundings.
- In these low-mass, high-sSFR systems atomic gas is not the dominant baryonic component, breaking the usual expectation that dwarf galaxies are H I-dominated.
- The low detection rate at high O32 means that selecting galaxies by extreme ionization parameter efficiently picks out galaxies with little neutral gas, strengthening O32-based searches for Lyman-continuum leaker candidates.
- The H I content must be mapped at higher angular resolution before the detected systems can be used to constrain where LyC photons escape, since the FAST beam blends the blueberries with neighbors.
- Molecular gas follow-up is needed because the standard H I-to-H2 conversion for main-sequence galaxies may not apply in these low-metallicity, high-density starbursts.
Reading between the lines
- Going beyond the paper, the same upper-limit method could be tested statistically: if the true H I line widths of nondetected dwarfs follow the roughly 20-30 km/s widths typical of dwarf galaxies, the 50 km/s assumption is conservative and the low gas fractions are robust; if they are broader, the offset from main-sequence relations shrinks.
- The paper's logic implies that the fraction of O32 greater than or similar to 10 dwarfs with detectable H I should drop sharply with decreasing stellar mass; future wide-field H I surveys could map this boundary and predict which systems are most likely to leak Lyman-continuum radiation.
- If these galaxies are really depleting their gas in about 0.5 Gyr, they may be caught in the final, gas-poor phase of a merger-driven starburst, connecting the short depletion times to the disturbed H I morphologies seen in interferometric images of individual Green Peas and blueberries.
- A direct measurement of molecular gas in a handful of blueberries would test whether the atomic reservoir is small because it has already been converted into H2, which would make the low H I content a sign of rapid fuelling rather than of genuinely little gas.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This Letter reports FAST H I 21 cm observations of 28 blueberry galaxies, a class of low-redshift, compact, high-sSFR starbursts considered local analogs of high-redshift Lyα emitters. The authors detect H I toward two objects, J1026+0426 and J1132+0809, and derive 3σ upper limits for the remaining 26 galaxies. Using Equation (6), which assumes a 50 km/s line width for nondetections, they obtain a median 3σ upper limit on H I mass of about 2.0×10^8 M_sun, a median f_HI upper limit of 10.2, and a median τ_HI upper limit of 0.53 Gyr. They compare these values with the f_HI-sSFR and f_HI-M* relations for main-sequence galaxies and conclude that blueberries have low atomic gas fractions and depletion times roughly an order of magnitude below local main-sequence values, with a low detection rate of 2/28 (7.1%).
Significance. If the quantitative conclusions hold, this paper would strengthen the emerging picture that extremely compact, high-sSFR dwarf galaxies such as blueberries are gas-poor, depleting their atomic gas reservoirs on sub-Gyr timescales. The main strengths are the large sample of blueberries observed in H I, the use of the standard H I mass formula, explicit discussion of neighbor contamination, and comparison with external relations from Catinella et al. (2018) and Kanekar et al. (2021). The low detection rate itself is likely robust to many caveats. However, the quantitative offsets in f_HI and τ_HI rest on an assumed line width that is not supported by the detected profiles, and on uncertain attribution of the two detections; these need to be addressed before the central quantitative claims can be accepted.
major comments (3)
- [§3, Eq. (6) and Table 1] The upper limits on M_HI, and hence the median f_HI and τ_HI claims, assume a fixed line width of 50 km/s for all nondetections. This assumption is not supported by the two detected profiles, whose FWHMs are 148.8±4.4 and 124.4±7.3 km/s (Table 1 notes). Because the 3σ integrated-flux limit scales as sqrt(W/50), a realistic W≈150 km/s raises every upper limit by a factor of 1.73: the median M_HI upper limit becomes ~3.5×10^8 M_sun, the median log f_HI upper limit moves from ~1.00 to ~1.25, above the extrapolated main-sequence value of ~1.16 at log M*=7.3, and the median τ_HI upper limit rises from 0.53 to ~0.92 Gyr. The claimed order-of-magnitude deficit in τ_HI is therefore not robust to the line-width prior. Please justify the 50 km/s choice using measured widths from BBs, GPs, or low-mass BCDs, or present the upper limits as a function of assumed line width.
- [§3, Table 1 notes, Fig. 3] Both detections carry substantial attribution uncertainty. For J1026+0426, the quoted log M_HI=7.73 is obtained by assuming that the BB-to-neighbor H I mass ratio equals their stellar mass ratio, while the neighboring Sbc galaxy's optical velocity is closer to the detected V_HI than the BB's optical velocity is. For J1132+0809, the lower value log M_HI=8.24 assumes the 12″ companion is at the H I velocity, and the alternative with no companion at that redshift gives log M_HI=9.56 (Table 1 note d). Because only two detections are reported, both scenarios should be propagated through the f_HI and τ_HI comparisons in Fig. 3 and the summary, not only in table notes; otherwise the reader cannot determine how much of the 'gas-poor blueberries' conclusion depends on a specific assumption about the neighbors.
- [§3.1, Fig. 3(a)] The conclusion that BBs have lower f_HI than main-sequence relations relies on comparing the median 3σ upper limit (10.2) with the extrapolated expected value (14.6). A 3σ upper limit is not a measurement, and the difference of only ~0.16 dex is small; many individual upper limits in Table 1 are above the expected M_HI at their stellar mass, e.g., J0820+5431 has log M_HI<7.94 against an expected log M_HI≈7.76. The claim that BBs 'tend to have lower f_HI' needs a proper upper-limit treatment, such as a Kaplan-Meier or survival-analysis estimate of the median, or a binomial test of the fraction of limits below the expected relation, before it can be regarded as established.
minor comments (3)
- [§3.1] In the paragraph discussing f_HI for the detected sources, 'J1123+0809' should be 'J1132+0809'; J1123+2050 is a nondetection and the text is about the detected source with a neighboring blue galaxy.
- [§3.3] The statement that 'all 26 BBs have an O32 ratio greater than 8.3' would be clearer if it specified that this refers to the 26 sources with available O[II] fluxes; Table 1 lists two sources (J1136+3427 and J1323−0132) with no O32 value.
- [Eq. (6)] The line-width dependence in Equation (6) is presented as a fixed numerical factor with Δv; please state explicitly whether the 50 km/s width is intended as a boxcar, Gaussian FWHM, or other profile shape, since the conversion from per-channel noise to integrated-flux limit depends on this choice.
Circularity Check
No significant circularity: the H I mass and upper-limit calculations use standard formulae with an explicit line-width prior, comparisons rely on external galaxy relations, and the only self-citation is a comparison sample, not load-bearing.
full rationale
The paper's central claims are observational: a 2/28 detection rate, median 3-sigma H I mass upper limits, and gas fractions and depletion times derived from those limits. The H I mass formula (Eq. 4) is the standard relation. The nondetection upper limits (Eq. 6) adopt an explicit 50 km/s line-width assumption; this is an analysis prior, not a parameter fitted to the target result. Changing the assumed width scales the upper limits, but that is a robustness concern about the sensitivity estimate, not a circular reduction of the conclusion to its input. The comparisons to f_HI-sSFR and f_HI-M* relations use external benchmarks from Catinella et al. (2018) and Kanekar et al. (2021), and the depletion-time comparison to main-sequence galaxies uses Saintonge et al. (2016) and Catinella et al. (2018). The Chandola et al. (2024) self-citation appears as a comparison sample for BCDs and for the method of estimating stellar masses of neighbors; it does not supply the central 'low gas fraction' or 'short depletion time' result. The paper does not invoke any uniqueness theorem, smuggle in an ansatz by citation, or rename a known result. The two detections are deblended from neighboring galaxies using an explicit stellar-mass-ratio assumption, which is a stated model choice rather than a forced equivalence. Overall, the derivation chain is self-contained against external data, and no step reduces to its own input by construction. The single point of caution is the 50 km/s line-width prior for nondetections, but this is an assumption with stated sensitivity implications, not circularity.
Assumptions & free parameters
free parameters (1)
- Assumed H I line width =
50 km/s
assumptions (5)
- domain assumption The 50 km/s line width approximates the true H I profiles of blueberries.
- ad hoc to paper For J1026+0426, the H I mass ratio of the blueberry to its companion equals their stellar mass ratio.
- ad hoc to paper The photometric redshift of the companion to J1132+0809 is not at the H I velocity, so most of the detected H I belongs to the blueberry.
- domain assumption The FAST 3-arcmin beam cannot spatially separate the blueberry from nearby galaxies, so all line flux within the beam is attributed using the mass-ratio and redshift assumptions.
- standard math The standard H I mass formula and the ON-OFF calibration equations are valid.
Cite this review
Pith. "Pith review of FAST H I 21 cm study of blueberry galaxies." pith.science (2026). https://pith.science/paper/GVE4BR67
@misc{pith2026241113527,
author = {Pith},
title = {Pith review of: FAST H I 21 cm study of blueberry galaxies},
year = {2026},
howpublished = {\url{https://pith.science/paper/GVE4BR67}},
note = {Machine review of arXiv:2411.13527}
}
abstract
Green Peas (GPs) and blueberry galaxies (BBs) are thought to be local analogs ($z<$0.1) of high redshift Ly$\alpha$ emitters. H I study of these can help us understand the star formation in the primordial Universe. In this Letter, we present the results of H I 21 cm study of 28 high specific star formation rate (sSFR $\gtrsim$10$^{-8}$ yr$^{-1}$) BBs at $z\lesssim$0.05 with the Five-hundred-meter Aperture Spherical radio Telescope. We report significant H I detection towards two BBs namely J1026+0426 and J1132+0809, and discuss possible H I contribution from neighboring galaxies. The median 3$\sigma$ upper limit of $\sim$2.0$\times$10$^{8}$ M$_{\odot}$ was obtained on H I mass for galaxies with nondetections. We find BBs tend to have lower H I-to-stellar mass ratio or gas fraction ($f_{\rm HI}$) than expected from $f_{\rm HI}$-sSFR and $f_{\rm HI}$-$M_{\ast}$ relations for main-sequence galaxies. The BBs also have a median 3$\sigma$ upper limit on H I gas depletion time scale ($\tau_{\rm HI}$) $\sim$0.5 Gyr, about 1 order of magnitude lower than $\tau_{\rm HI}$ for local main-sequence galaxies. We find a significantly low H I detection rate of 2/28 (7.1$^{+9.4}_{-4.6}$ \%) towards these galaxies, which is similar to previous H I studies of low redshift GPs of high ionization parameter indicator, O32 $\equiv$O[{\sc iii}]$\lambda$5007/O[{\sc ii}]$\lambda$3727 ratios $\gtrsim$10.
Figures
Reference graph
Works this paper leans on
-
[1]
Abazajian, K. N., Adelman-McCarthy, J. K., Ag¨ ueros, M. A., et al. 2009, ApJS, 182, 543, doi: 10.1088/0067-0049/182/2/543
-
[2]
2020, ApJS, 249, 3, doi: 10.3847/1538-4365/ab929e Astropy Collaboration, Robitaille, T
Ahumada, R., Allende Prieto, C., Almeida, A., et al. 2020, ApJS, 249, 3, doi: 10.3847/1538-4365/ab929e Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Collabor...
-
[3]
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5, doi: 10.1086/130766
doi:10.1086/130766 1981
-
[4]
2009, MNRAS, 399, 1191, doi: 10.1111/j.1365-2966.2009.15383.x
Cardamone, C., Schawinski, K., Sarzi, M., et al. 2009, MNRAS, 399, 1191, doi: 10.1111/j.1365-2966.2009.15383.x
arXiv 2009
-
[5]
2018, MNRAS, 476, 875, doi: 10.1093/mnras/sty089
Catinella, B., Saintonge, A., Janowiecki, S., et al. 2018, MNRAS, 476, 875, doi: 10.1093/mnras/sty089
-
[6]
2024, MNRAS, 527, 603, doi: 10.1093/mnras/stad3018
Chandola, Y., Li, D., Tsai, C.-W., et al. 2024, MNRAS, 527, 603, doi: 10.1093/mnras/stad3018
-
[7]
2024, MNRAS, 531, 5140, doi: 10.1093/mnras/stae1490
Dutta, S., Bera, A., Bait, O., et al. 2024, MNRAS, 531, 5140, doi: 10.1093/mnras/stae1490
-
[8]
Flury, S. R., Jaskot, A. E., Ferguson, H. C., et al. 2022, ApJ, 930, 126, doi: 10.3847/1538-4357/ac61e4
Show all 36 references
-
[9]
2007, ApJ, 671, 278, doi: 10.1086/522955
Gawiser, E., Francke, H., Lai, K., et al. 2007, ApJ, 671, 278, doi: 10.1086/522955
2007 doi
-
[10]
1986, ApJ, 303, 336, doi: 10.1086/164079 Gil de Paz, A., Madore, B
Gehrels, N. 1986, ApJ, 303, 336, doi: 10.1086/164079 Gil de Paz, A., Madore, B. F., & Pevunova, O. 2003, ApJS, 147, 29, doi: 10.1086/374737
1986 doi
-
[11]
R., Millman, K
Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2
2020 doi
-
[12]
L., & Erb, D
Henry, A., Scarlata, C., Martin, C. L., & Erb, D. 2015, ApJ, 809, 19, doi: 10.1088/0004-637X/809/1/19 10
2015 doi
-
[13]
P., Giovanelli, R., & Brinchmann, J
Huang, S., Haynes, M. P., Giovanelli, R., & Brinchmann, J. 2012, ApJ, 756, 113, doi: 10.1088/0004-637X/756/2/113
2012 doi
-
[14]
Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55
2007 doi
-
[15]
I., Schaerer, D., Thuan, T
Izotov, Y. I., Schaerer, D., Thuan, T. X., et al. 2016, MNRAS, 461, 3683, doi: 10.1093/mnras/stw1205
2016 doi
-
[16]
I., Schaerer, D., Worseck, G., et al
Izotov, Y. I., Schaerer, D., Worseck, G., et al. 2018a, MNRAS, 474, 4514, doi: 10.1093/mnras/stx3115
-
[17]
I., Worseck, G., Schaerer, D., et al
Izotov, Y. I., Worseck, G., Schaerer, D., et al. 2021, MNRAS, 503, 1734, doi: 10.1093/mnras/stab612 —. 2018b, MNRAS, 478, 4851, doi: 10.1093/mnras/sty1378
2021 doi
-
[18]
E., & Oey, M
Jaskot, A. E., & Oey, M. S. 2013, ApJ, 766, 91, doi: 10.1088/0004-637X/766/2/91
2013 doi
-
[19]
2020, Research in Astronomy and Astrophysics, 20, 064, doi: 10.1088/1674-4527/20/5/64
Jiang, P., Tang, N.-Y., Hou, L.-G., et al. 2020, Research in Astronomy and Astrophysics, 20, 064, doi: 10.1088/1674-4527/20/5/64
2020 doi
-
[20]
E., & Yang, H
Jiang, T., Malhotra, S., Rhoads, J. E., & Yang, H. 2019, ApJ, 872, 145, doi: 10.3847/1538-4357/aaee8a
2019 doi
-
[21]
2021, ApJL, 913, L15, doi: 10.3847/2041-8213/abfb76
Kanekar, N., Ghosh, T., Rhoads, J., et al. 2021, ApJL, 913, L15, doi: 10.3847/2041-8213/abfb76
2021 doi
- [22]
-
[23]
S., Staveley-Smith, L., Kilborn, V
Koribalski, B. S., Staveley-Smith, L., Kilborn, V. A., et al. 2004, AJ, 128, 16, doi: 10.1086/421744 Le Reste, A., Cannon, J. M., Hayes, M. J., et al. 2024, MNRAS, 528, 757, doi: 10.1093/mnras/stad3910
2004 doi
-
[24]
H., Jaskot, A
McKinney, J. H., Jaskot, A. E., Oey, M. S., et al. 2019, ApJ, 874, 52, doi: 10.3847/1538-4357/ab08eb
2019 doi
-
[25]
2014, MNRAS, 442, 900, doi: 10.1093/mnras/stu902
Nakajima, K., & Ouchi, M. 2014, MNRAS, 442, 900, doi: 10.1093/mnras/stu902
2014 doi
-
[26]
2011, International Journal of Modern Physics D, 20, 989, doi: 10.1142/S0218271811019335
Nan, R., Li, D., Jin, C., et al. 2011, International Journal of Modern Physics D, 20, 989, doi: 10.1142/S0218271811019335
2011 doi
-
[27]
A., et al
Paalvast, M., Verhamme, A., Straka, L. A., et al. 2018, A&A, 618, A40, doi: 10.1051/0004-6361/201832866
2018 doi
-
[28]
N., et al
Purkayastha, S., Kanekar, N., Chengalur, J. N., et al. 2022, ApJL, 933, L11, doi: 10.3847/2041-8213/ac7522
2022 doi
-
[29]
2016, MNRAS, 462, 1749, doi: 10.1093/mnras/stw1715
Saintonge, A., Catinella, B., Cortese, L., et al. 2016, MNRAS, 462, 1749, doi: 10.1093/mnras/stw1715
2016 doi
- [30]
-
[31]
Taylor, M. B. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 347, Astronomical Data Analysis Software and Systems XIV, ed. P. Shopbell, M. Britton, & R. Ebert, 29
2005
-
[32]
X., Goehring, K
Thuan, T. X., Goehring, K. M., Hibbard, J. E., Izotov, Y. I., & Hunt, L. K. 2016, MNRAS, 463, 4268, doi: 10.1093/mnras/stw2259
2016 doi
-
[33]
E., & Wang, J
Yang, H., Malhotra, S., Rhoads, J. E., & Wang, J. 2017a, ApJ, 847, 38, doi: 10.3847/1538-4357/aa8809
-
[34]
2017b, ApJ, 844, 171, doi: 10.3847/1538-4357/aa7d4d
Yang, H., Malhotra, S., Gronke, M., et al. 2017b, ApJ, 844, 171, doi: 10.3847/1538-4357/aa7d4d
-
[35]
2019, Science China
Zhang, K., Wu, J., Li, D., et al. 2019, Science China
2019
-
[36]
Physics, Mechanics, and Astronomy, 62, 959506, doi: 10.1007/s11433-019-9383-y
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.