{"id":"d6c6b15b-2ab9-40aa-a032-2b5598c45378","arxiv_id":"2411.09903","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Combining HIPASS, ALFALFA, and FASHI over 76% of the sky yields the most complete local HI mass function, with alpha=-1.30, log(M*/h70^-2 Msun)=9.86, and Omega_HI=4.54e-4 h70^-1.","lead":"Astronomers combined three radio surveys covering 76% of the sky to measure the local hydrogen mass function, the count of galaxies by their atomic hydrogen content. The result includes the first such measurement from the FAST FASHI survey and pins down the cosmic hydrogen abundance with lower cosmic variance than any single survey.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"HIPASS effective volumes are extended to z=0.05 even though HIPASS observed only z<0.042, inflating the largest sky-area component and biasing the combined HIMF and Omega_HI.","rationale":"The reader's verdict was CONDITIONAL, with the primary caveat being the universality of the HIMF across the surveyed sky. I agree that universality is a real limitation, but it is a physical assumption shared by every HIMF measurement and the paper explicitly discusses large-scale structure effects. A more load-bearing and more easily settled problem is the HIPASS Vmax extrapolation in Section 3.3. The paper explicitly assumes no evolution between z=0.042 and z=0.05 to let HIPASS effective volumes reach z=0.05, but Vmax is defined by the survey's actual selection function, not by cosmic evolution. HIPASS observed only z<0.042, and the southern sky is not observed by the other surveys beyond that redshift, so counting that volume undercounts the HIPASS density. Because HIPASS carries 58% of the combined sky area, the effect on the weighted total is sizeable—up to roughly 70% at the knee—and likely shifts Omega_HI by more than the quoted statistical error. The reader's concern about the 'Vmax using total sky area' wording is actually not an inconsistency: using the total sky area for every galaxy is algebraically equivalent to a sky-area-weighted average of the per-survey HIMFs. The issue is specifically the missing redshift truncation for HIPASS. The proposed check, recomputing with Vmax capped at z=0.042, is straightforward and would settle whether the headline numbers survive. If the shift is large, the claim of an unbiased Omega_HI needs revision; if small, the method is validated. In either case the verdict should remain CONDITIONAL pending the corrected volume calculation, so I keep the reader's conditional verdict unchanged.","tokens_in":18914,"tokens_out":14081,"duration_ms":145935,"concrete_test":"Recompute the total HIMF via Eq. 5 with HIPASS Vmax truncated at the survey redshift limit z=0.042 (i.e., Vmax_i = A_HIPASS * integral from 0 to min(z_lim, z_flux_max) of dV/dz dz), keeping ALFALFA and FASHI at z=0.05. Compare the new best-fit Schechter parameters and the integral of Eq. 8 to the quoted values; if Omega_HI shifts by more than about 10%, or phi_s by more than about 15%, the reported central values are not robust to the correct HIPASS effective volume.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 computes Vmax for HIPASS galaxies out to z=0.05, stating that no evolution is assumed between z=0.042 and 0.05. In the 1/Vmax estimator, Vmax is the volume actually surveyed in which a galaxy could have been detected; HIPASS did not observe z>0.042, and the southern footprint is not covered by ALFALFA or FASHI beyond that redshift. The no-evolution assumption cannot turn unobserved volume into usable volume. With HIPASS contributing 18,291 deg^2 (58% of the combined sky area), an un-truncated Vmax directly suppresses the HIPASS contribution in the sky-area-weighted total: for galaxies bright enough to reach the survey redshift limit, the correct Vmax is smaller by a factor (0.042/0.05)^3 ≈ 0.59, meaning the HIPASS HIMF is underweighted by roughly 70% at and above the knee. This propagates into the total HIMF, the best-fit phi_s, and the headline Omega_HI = 4.54e-4, and the quoted uncertainties do not include this systematic. This is a concrete selection-function error, unlike the unavoidable question of whether the HIMF is universal across the sky.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first HIMF measurement for the FASHI survey and combines HIPASS, ALFALFA, and FASHI (split into north and south) to derive a 'total' HIMF over 31,528 deg^2 at 0<z<0.05, claiming this is the most complete local HIMF to date. The authors apply a uniform distance framework (Cosmicflows-4), a common completeness treatment based on S21 and W50, and the 1/Vmax estimator. They fit the total HIMF with a single Schechter function (alpha=-1.30, log Ms=9.86, phi_s=6.58e-3, Omega_HI=4.54e-4) and also present a double-Schechter fit. The central claim is that combining the three surveys suppresses cosmic variance and yields an unbiased local HIMF.","tokens_in":19203,"tokens_out":21261,"duration_ms":201320,"significance":"If the central measurement is sound, this would be a valuable reference product: the first FASHI HIMF, a uniform re-analysis of three major HI surveys, and a nearly all-sky local HIMF. The completeness analysis is careful, and the agreement of the ALFALFA and FASHI plateau levels in Fig. 3 is a genuine cross-validation. The use of a consistent distance and completeness pipeline across all samples is a real strength. However, the current combination scheme in Sect. 3.3 is not the standard volume-weighted 1/Vmax combination, and this affects the headline HIMF and Omega_HI. The significance therefore depends on whether the Vmax treatment is corrected and the results re-derived.","major_comments":[{"comment":"The sentence 'When combining the four samples, we calculated Vmax for each galaxy using the total sky area of 31528 deg2' is not a valid application of the 1/Vmax estimator for a union of surveys with disjoint footprints. A galaxy detected in HIPASS could not have been detected over the ALFALFA or FASHI footprints, so its maximum volume should be (Omega_HIPASS/4pi)*(4pi/3) D_max^3, not (Omega_total/4pi)*(4pi/3) D_max^3. Using the total sky area for every galaxy is algebraically equivalent to forming the sky-area-weighted average of the four sample HIMFs; indeed, the total values in Table 3 (e.g., log M_HI=9.5: 1.279e-2) match the area-weighted sum of the four sample entries. This is not the volume-weighted combination that gives the number density in the combined survey volume. Because HIPASS has the largest area (18,291 deg^2) but the shallowest depth, area-weighting overweights HIPASS relative to its volume share, pulling the combined HIMF down where HIPASS amplitudes are low. The authors should recompute each galaxy's Vmax using the sky area of its own survey and then combine the samples with volume weighting, or provide an explicit statistical justification for why the area-weighted average is the appropriate estimator for the cosmic HIMF. The fitted parameters in Table 4 and Omega_HI in Section 4 will likely change, and the quoted uncertainties do not capture this systematic.","section":"Sect. 3.3, Eq. (5) and Table 3"},{"comment":"The statement 'Since the HIPASS sample is limited to a slightly smaller redshift of z < 0.042, we can assume that there is no evolution within the redshift range of 0.042 < z < 0.05' is used to extend HIPASS Vmax to z=0.05. In the 1/Vmax method, Vmax is the volume actually surveyed in which a source could have been detected; HIPASS did not observe z>0.042. The no-evolution assumption only says the space density is the same in that interval; it cannot make unobserved volume available for detection. HIPASS galaxies with D_max > D(z=0.042) should have Vmax truncated at the survey redshift limit. The numerical impact is modest because HIPASS's flux limit gives D_max ~ 72 Mpc (z ~ 0.017) for a knee-mass galaxy, so only the highest masses (log M_HI >~ 10.65) are affected, but the reasoning is incorrect and the effect on the high-mass bin should be quantified explicitly.","section":"Sect. 3.3, HIPASS redshift limit"}],"minor_comments":[{"comment":"There is a typo: 'ALAFLFA' should be 'ALFALFA'.","section":"Sect. 1, last paragraph"},{"comment":"The definition of Vmax_i is not stated explicitly. Please add an equation or sentence defining Vmax_i = (Omega_survey/4pi) * (4pi/3) D_max^3 for each survey, and clarify the units and sky area used.","section":"Sect. 3.3, Eq. (5)"},{"comment":"The table caption should state that the 'Total' column is the area-weighted combination of the four samples; if the method is revised after the major comment above, the caption and values will need to be updated.","section":"Table 3 caption"},{"comment":"The phrase 'increasing MH i by 0.06 dex and decreasing phi by -0.09 dex' has a double negative; it should read 'decreasing phi by 0.09 dex'.","section":"Sect. 5.2"},{"comment":"The improvement in reduced chi-square from 56/17 to 24/15 is reported without a statistical significance test. A simple delta-chi-square or F-test probability would help support the claim that the double Schechter is preferred.","section":"Sect. 4, double Schechter fit"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely question and the FASHI HIMF is a useful new product. The completeness work is careful and the plateau agreement in Fig. 3 is a nice validation. The main concern is the Vmax combination scheme: the use of the total sky area for every galaxy is not the standard 1/Vmax estimator and biases the combined HIMF and Omega_HI. This is fixable by recomputing per-survey Vmax and volume-weighting, and should be done before publication. The HIPASS redshift-limit extension is a secondary, more easily fixed issue. I would not reject the paper on the basis of the redshift-limit point alone, but the area-weighting issue is load-bearing for the headline numbers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The punchline: the first FASHI HIMF and a three-survey combination covering 76% of the sky are genuinely new and useful. But the 1/Vmax computation for the combined sample as described in Section 3.3 is wrong, and it biases the headline Omega_HI in a way the quoted errors do not reflect.\n\nWhat is actually good: the FASHI measurement alone is worth the paper, and the uniform re-processing of HIPASS, ALFALFA, and FASHI with the same distance model, completeness fitting, and Schechter fitting is the right approach. The completeness analysis is careful; the agreement of the ALFALFA and FASHI plateau levels in Fig. 3 is a real validation. The distance-robustness test in Fig. 6 is also useful. The resulting parameters (alpha = -1.30, log Ms = 9.86, Omega_HI = 4.54e-4) are plausible and consistent with Guo et al. (2023).\n\nThe soft spot is load-bearing. Section 3.3 says Vmax is computed for all galaxies using the total sky area of 31,528 deg², and for HIPASS galaxies it 'assumes no evolution' out to z=0.05 even though HIPASS observed only z<0.042. That is not a no-evolution assumption; it is an extrapolation into volume HIPASS never observed. For a HIPASS galaxy, Vmax should be based on the HIPASS footprint and redshift limit. Using the pooled sky area and the extended redshift limit lowers HIPASS's weight by roughly a factor of three at and above the knee, since HIPASS is 58% of the combined sky. That changes the combined HIMF, the Schechter fits, and Omega_HI. The fix is straightforward: compute each survey's contribution with its own Vmax and combine, or use a joint likelihood that respects each survey's footprint. The individual HIMFs in Table 3 appear to use per-survey volumes, so the total should be their area-weighted average, not a pooled Vmax.\n\nThe second concern is less concrete but real: 'unbiased' relies on the HIMF being universal across the surveyed sky. Fig. 5 shows factor-level north/south differences at M_HI < 1e9, attributed to cosmic variance. Environmental dependence is another possibility, though it does not invalidate the measurement as an average. The citation pattern is appropriate; the FASHI-related self-citations are justified by the dataset. The paper's own final note about limited low-mass statistics is fair, but it does not cover the HIPASS Vmax problem.\n\nWho this is for: anyone working on local HI surveys, WALLABY, MIGHTEE, or FASHI follow-ups. It deserves a serious referee. I would not desk-reject it, but I would send it back for major revision to correct the Vmax treatment and re-state the 'unbiased' claim.","headline":"First FASHI HIMF and an ambitious three-survey combination, but the pooled 1/Vmax treatment underweights HIPASS and the quoted Omega_HI is biased.","tokens_in":19784,"tokens_out":4783,"would_cite":true,"duration_ms":46029,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By combining HIPASS, ALFALFA, and FASHI into one 76%-sky sample processed with identical methods, this paper establishes the most complete local-universe HI mass function to date, including the first FASHI measurement.","keywords":["HI mass function","local universe","HIPASS","ALFALFA","FASHI","Schechter function","cosmic HI abundance","1/Vmax method"],"falsifier":"A uniform deep HI survey covering both hemispheres could measure the HI mass function separately in the north and south: if the two disagree by more than the quoted uncertainties below $10^9\\,M_\\odot$, the universal-HIMF premise fails. The no-evolution assumption could be tested directly by measuring the HIMF in the $0.042<z<0.05$ slice with a deep survey; a measured slope or amplitude change there would invalidate the HIPASS volume extension.","tokens_in":18666,"feed_emoji":"📡","tokens_out":10311,"duration_ms":87868,"temperature":0.7,"pith_summary":"This paper sets out to measure how many galaxies of each neutral-hydrogen mass fill the local universe by pooling three blind HI surveys—HIPASS, ALFALFA, and the new FASHI catalogue—into a single $31{,}528$ square-degree sample covering 76% of the sky. To make the combination meaningful, the authors recalculate distances, completeness, and the mass function in exactly the same way for every survey. The payoff is a local HI mass function with cosmic variance strongly suppressed: a Schechter fit with low-mass slope $\\alpha=-1.30$, knee mass $\\log(M_s)=9.86$, and cosmic HI abundance $\\Omega_{\\rm HI}=4.54\\times10^{-4}$. The paper also reports the first HI mass function from FASHI and finds that a double Schechter function with two knee masses fits the data better than a single one. If correct, this gives future HI surveys a local benchmark for the gas content of galaxies and how it evolves.","feed_headline":"Three surveys, 76% of sky, one local HI mass function","feed_subtitle":"Combining HIPASS, ALFALFA, and FASHI pins Omega_HI at 4.54e-4 and finds a double-Schechter shape.","key_machinery":"The load-bearing machinery is a unified $1/V_{\\max}$ estimator run on the merged sample. Each galaxy contributes weight $f_{\\rm rms,i}/[C(S_{21,i}|W_{50,i})V_{\\max,i}]$, where $C$ is the completeness from error-function fits to the $S_{21}^{3/2}\\,dn/d\\log S_{21}$ plateau in each line-width bin, $V_{\\max}$ is the maximum volume in which the galaxy would still clear the 50% completeness limit, and $f_{\\rm rms}$ is a pixel-level weight that corrects FASHI's uneven exposure. Distances come from the same Cosmicflows-4 model for all three surveys, and $V_{\\max}$ is computed over the full $31{,}528\\,\\mathrm{deg}^2$ area. This uniformity is what lets the authors compare and then combine surveys with very different selection functions.","core_discovery":"Combining HIPASS, ALFALFA, FASHI north, and FASHI south at $0<z<0.05$ with all four samples reduced through one pipeline gives the most complete HI mass function of the local universe measured so far. A single Schechter fit to the total yields $\\alpha=-1.30\\pm0.01$, $\\log(M_s/h_{70}^{-2}\\,M_\\odot)=9.86\\pm0.01$, and $\\phi_s=(6.58\\pm0.23)\\times10^{-3}\\,h_{70}^3\\,\\mathrm{Mpc}^{-3}\\,\\mathrm{dex}^{-1}$, which integrates to $\\Omega_{\\rm HI}=(4.54\\pm0.20)\\times10^{-4}\\,h_{70}^{-1}$. A double Schechter function with a shared slope and two knees at $\\log(M_{s1})=9.96$ and $\\log(M_{s2})=9.65$ fits better, with the reduced $\\chi^2$ dropping from $56/17$ to $24/15$, especially improving the high-mass end. The authors interpret the differences between individual survey HIMFs at low masses as large-scale structure effects and argue that averaging over 76% of the sky suppresses cosmic variance enough to give an unbiased local estimate.","pith_inferences":["If the universal-HIMF premise holds, the north/south differences seen in the individual surveys are pure sample variance and should shrink as deeper surveys cover more volume; if they persist, the combined function is an average of environmentally different HIMFs and its low-mass end is not a single universal curve.","The two Schechter knees around $10^{9.6}$ and $10^{10}\\,M_\\odot$ map naturally onto populations split by halo mass; a direct test would be to recompute the HIMF using only central or only satellite galaxies from a group catalogue.","The pixel-weighting correction for FASHI's uneven exposure is a template for other time-filler surveys whose depth varies across the sky.","Applying the same pipeline to future FASHI data releases, or to other all-sky HI surveys, would show how much of the remaining scatter is methodological rather than cosmic."],"forward_implications":["Future HI surveys can compare their deeper or higher-redshift measurements against a local benchmark that is not tied to one patch of sky.","The local cosmic HI abundance $\\Omega_{\\rm HI}\\approx4.5\\times10^{-4}$ provides an anchor for studies of HI evolution towards $z\\sim0.2$.","The better-fitting double Schechter function implies that the local HI mass function may be shaped by two distinct galaxy populations, possibly centrals and satellites.","The individual survey HIMFs provide a direct measure of how much large-scale structure biases single-survey estimates at low HI masses.","The FASHI HIMF presented here is the first from that survey and can be updated as future FASHI releases cover more sky."],"supporting_citations":[{"why":"Supplies the HIPASS survey catalogue and its completeness and redshift characteristics used as one of the three combined samples.","marker":"Meyer et al. (2004)"},{"why":"Provides the completeness method, error-function fits to the flux-density plateau, that is adopted uniformly for all three surveys.","marker":"Haynes et al. (2011)"},{"why":"Establishes the 1/Veff method, the comparison with 1/Vmax, and the formula used to integrate the Schechter fit for Omega_HI.","marker":"Martin et al. (2010)"},{"why":"Provides the final ALFALFA HIMF, the variance formula used for the total, and evidence that separate ALFALFA regions differ because of cosmic variance.","marker":"Jones et al. (2018)"},{"why":"Releases the first FASHI catalogue and documents its depth, noise, and sky coverage, the new survey whose HIMF is measured here.","marker":"Zhang et al. (2024)"},{"why":"Provides the Cosmicflows-4 Wiener-filter distance estimates applied uniformly to all three surveys.","marker":"Valade et al. (2024)"},{"why":"Gives the ALFALFA 100% HIMF with 90% completeness cuts that the combined result is compared against and found consistent with.","marker":"Guo et al. (2023)"},{"why":"The original HIPASS HIMF measurement, used as the southern-sky baseline and comparison at the high-mass end.","marker":"Zwaan et al. (2005)"}],"fun_headline_variants":["Most complete local HI mass function from 76% of sky","Combining three surveys reveals double-Schechter HI mass function","FASHI plus HIPASS and ALFALFA: unbiased HI mass function","76% sky coverage yields definitive HI mass function shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The measurement assumes the intrinsic HI mass function is the same across the whole surveyed sky, so that averaging over regions with different large-scale structure recovers the true cosmic function; it also assumes no HI evolution between $z=0.042$ and $z=0.05$, which lets the HIPASS volume be extended to the common redshift limit.","fun_headline_variants_meta":{"raw":{"variants":["Most complete local HI mass function from 76% of sky","Combining three surveys reveals double-Schechter HI mass function","FASHI plus HIPASS and ALFALFA: unbiased HI mass function","76% sky coverage yields definitive HI mass function shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1545,"prompt_tokens":1097,"completion_tokens":448,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":713,"completion_tokens_details":{"reasoning_tokens":375}},"tokens_in":713,"tokens_out":448,"duration_ms":5241,"temperature":1.0,"reasoning_tokens":375,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:11:22.921274+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A uniform deep HI survey covering both hemispheres could measure the HI mass function separately in the north and south: if the two disagree by more than the quoted uncertainties below $10^9\\,M_\\odot$, the universal-HIMF premise fails. The no-evolution assumption could be tested directly by measuring the HIMF in the $0.042<z<0.05$ slice with a deep survey; a measured slope or amplitude change there would invalidate the HIPASS volume extension.","supporting_citations":[{"cited_title":"J., Zwaan, M","cited_arxiv_id":null,"evidence_quote":"Supplies the HIPASS survey catalogue and its completeness and redshift characteristics used as one of the three combined samples."},{"cited_title":"G., Haynes, M","cited_arxiv_id":null,"evidence_quote":"Provides the final ALFALFA HIMF, the variance formula used for the total, and evidence that separate ALFALFA regions differ because of cosmic variance."},{"cited_title":"I., Pomarède, D., et al","cited_arxiv_id":null,"evidence_quote":"Provides the Cosmicflows-4 Wiener-filter distance estimates applied uniformly to all three surveys."},{"cited_title":"A., Meyer, M","cited_arxiv_id":null,"evidence_quote":"The original HIPASS HIMF measurement, used as the southern-sky baseline and comparison at the high-mass end."}],"review_version":1}