{"id":"1c0eea3c-0fec-4849-8866-574c2fd4291c","arxiv_id":"2507.18687","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In TNG50 Milky Way-like galaxies, ionized envelopes around high-velocity clouds carry about six times more mass than the neutral clouds and their infall rate is comparable to the galactic star formation rate.","lead":"Using the TNG50 cosmological simulation, the authors find that the ionized gas surrounding high-velocity clouds around Milky Way-like galaxies contains about six times more mass than the neutral gas and can fall in fast enough to balance up to 81% of the galaxies' star formation. The result suggests that the diffuse, ionized halo gas that is invisible in 21-cm maps may be the main fuel supply for star formation in disk galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 81% figure depends on Ṁ_ion from half-size radial envelopes without the ambient-CGM step-down used in the 3D envelopes, so random alignment may inflate the headline.","rationale":"I read the paper in good faith. The qualitative picture—ionized envelopes more massive than neutral HVCs, mostly prolate morphologies, and accretion rates that scale with SFR—is internally consistent and builds on the public TNG50 simulations with a clear method and mock-spectra validation. The authors honestly flag the cooling/condensation caveat, which means the headline 'can sustain star formation' is an upper-bound statement rather than a demonstrated pathway. The most load-bearing and checkable weakness is not the cooling step itself but the envelope definition used for the quantitative claim. The half-size method was designed to capture comoving gas, yet the paper's own step-down procedure acknowledges that dvi≤0.5 captures randomly aligned ambient CGM gas; because that step-down is applied to the morphological 3D envelopes and not obviously to the accretion-rate calculation, the population-level 81% and 5.6x values may include an uncontrolled ambient component. The Section 4 validation of half-size masses against 3D envelopes covers only the fiducial galaxy and leaves 20% of HVCs as outliers, so this is not a settled calibration. I therefore agree with the reader's conditional verdict but would sharpen the condition: recompute the population statistics with step-down 3D envelopes before relying on the 81% number. Minor numeric inconsistencies (47 vs 46 galaxies, and the left-panel fit equation printed with Ṁ_ion on the wrong side in the text) are real but do not change the central argument.","tokens_in":18474,"tokens_out":5057,"duration_ms":57423,"concrete_test":"Recompute per-galaxy Ṁ_ion for all 46 galaxies using the connected 3D envelopes with the step-down criterion (dvi starting at 0.5, decrement 0.05 until fewer than 2,000 cells) instead of the half-size radial weighting, and regenerate the right panel of Figure 7. If the peak or median Ṁ_ion/SFR shifts by more than about 20%, or if the neutral-versus-ionized ordering changes, the 81% headline is not robust to the envelope definition. As a null test, randomize cell velocities within the half-size spheres and recompute the alignment-weighted mass; the residual signal measures chance contamination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.2 computes Ṁ_ion by summing cells within twice the half-size radius, weighted by max(0,1−dvi), with dvi≤0.5 defining alignment. This is the basis for Figure 7 and the 81% claim. The same section states that dvi≤0.5 can produce 'unrealistically large ionized envelopes due to random alignment of ambient CGM gas cells,' and therefore the connected 3D envelopes used for morphology are stepped down in dvi by 0.05 until fewer than 2,000 cells. The accretion-rate calculation does not appear to apply that step-down, and the Section 4 validation of half-size masses against the 3D envelopes covers only the fiducial galaxy, with 20% of HVCs disagreeing by more than 15%. If even a fraction of the aligned hot gas in the half-size spheres is ambient CGM moving at similar velocity by chance, both the ~5.6x ionized-to-neutral mass ratio and the 81% Ṁ_ion/SFR peak are biased high. Separately, the paper explicitly defers cooling/condensation and cloud survival to future work, so 'sustain star formation' is conditional; the envelope-contamination issue is more load-bearing because it can be checked directly and could change the headline number rather than just its interpretation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper identifies high-velocity cloud (HVC) analogs in TNG50, defines ionized envelopes around them using velocity alignment (dvi) and temperature thresholds, and computes neutral and ionized masses and accretion rates for one fiducial galaxy and a population of ~47 Milky Way-like galaxies. The central claims are that the ionized envelopes contain ~5.6 times more mass than the neutral HVCs, that the ionized accretion rate can balance ~81% of the galactic star formation rate (with neutral HVCs balancing only ~11%), and that these rates scale with SFR and with each other. The paper also presents mock absorption spectra and covering fractions for comparison with observations, and characterizes the envelopes as mostly prolate.","tokens_in":18730,"tokens_out":6250,"duration_ms":71549,"significance":"If the quantitative claims hold, the paper provides a direct simulation-based argument that ionized CGM gas comoving with HVCs is the dominant accretion reservoir in Milky Way-like galaxies at z=0, strengthening and extending observational estimates such as those of Fox et al. (2019). Strengths include the direct use of simulation cells rather than fitted models, a population of dozens of MW-like galaxies, mock spectra generated with Trident for observational comparison, and explicit caveats about cooling, condensation, and resolution. The main risk is that the headline numbers (5.6x mass ratio, 81% SFR fraction, 6.8 vs 1.3 Msun/yr) depend on the definition of the ionized envelope, and the paper itself notes that the dvi<=0.5 criterion can select ambient CGM gas by chance. Because the accretion-rate estimator does not apply the same step-down decontamination used for the 3D morphological envelopes, the quantitative results require a robustness check before the central claim is fully established.","major_comments":[{"comment":"The accretion rates that feed the headline numbers are computed from the half-size radial envelopes using dvi<=0.5, without the step-down procedure that the same section applies to the 3D Voronoi envelopes. The text explicitly states that dvi<=0.5 can produce \"unrealistically large ionized envelopes due to random alignment of ambient CGM gas cells\", yet the step-down (decreasing dvi by 0.05 until fewer than 2,000 cells) is applied only to the connected 3D envelopes used for morphology, not to the mass and accretion-rate sums that underlie Figure 7 and the 81% and 5.6x numbers. Section 4's validation of the half-size method against the 3D envelopes covers masses only, only in the fiducial galaxy, and reports that 20% of HVCs disagree by more than 15%; it does not compare accretion rates. Please quantify the ambient-CGM contamination by recomputing Mdot_ion with the stepped-down 3D envelopes (or by subtracting a baseline alignment fraction for ambient gas) and report how the 81% peak, the 5.6x mass ratio, and the Figure 7 scalings change with that choice.","section":"Section 2.2"},{"comment":"The title \"Invisible Accretion: Ionized Envelopes of TNG50 HVCs can Sustain Star Formation\" and the concluding sentence assert a causal capacity, while the analysis does not model cooling, condensation, cloud survival, or mixing; the paper explicitly postpones the fate of these HVC analogs to future work. The abstract's parenthetical caveat about sufficient cooling and condensation is appropriate, but the title overstates the result. Please soften the title or add a sentence in the abstract and conclusions clarifying that the 81% figure is an upper-bound-style accretion-rate comparison, not a demonstrated star-formation pathway.","section":"Title and Abstract"}],"minor_comments":[{"comment":"The abstract states that the study covers 47 Milky Way-like galaxies, but Section 3.3 gives 52 total galaxies minus 6 excluded major mergers, yielding 46 galaxies, and the text says \"In these 46 galaxies\". Please reconcile this discrepancy.","section":"Section 3.3"},{"comment":"The text near the left panel of Figure 7 writes the fit as Mdot_neut = 3.16 x Mdot_ion + 2.02, which contradicts the figure's equation Mdot_ion = 3.16 x Mdot_neut + 2.02 and the surrounding claim that ionized accretion is about three times the neutral rate. The text should be corrected to match the figure.","section":"Section 3.3 and Figure 7"},{"comment":"Equation (3) has four free parameters (fmax, rc, alpha, tau), and the half-size is derived from the fit, but the paper does not report uncertainties on the fitted half-sizes or a sensitivity test of the alignment-fraction fit. A brief statement of the fitting uncertainties and their effect on the derived envelope masses would strengthen the robustness discussion.","section":"Section 2.2"},{"comment":"The sentence \"The HVC envelopes are infalling (by selection)\" is imprecise, because the envelope selection is based on velocity alignment rather than on radial infall; the accretion-rate calculation later restricts to infalling material. Please clarify whether the morphology statistics use all envelopes or only those of infalling HVCs.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The core technical risk is the ambient-CGM contamination in the accretion-rate estimator: the half-size dvi<=0.5 criterion is known from the paper's own morphology analysis to pick up randomly aligned CGM gas, and the step-down fix is not applied to the accretion-rate sums that produce the 81% and 5.6x headlines. This is directly checkable with the existing simulation data, and the authors are likely in a position to run the comparison against the stepped-down 3D envelopes. If that test preserves the qualitative and quantitative results, I would support acceptance; if the contamination is substantial, the central claim weakens to a qualitative statement that ionized material dominates neutral material, which is a lesser but still useful result. The 46/47 sample-size discrepancy and the reversed fit equation in Section 3.3 also need a careful proofing pass. No concerns about novelty or scope: the paper fits an astrophysics journal and engages properly with the observational literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core claim—that ionized high-velocity gas carries several times more accreting mass than the neutral HVCs, and that this inflowing component scales with SFR—is likely correct in its qualitative form. The quantitative headline, however, is not yet solid: the 81% figure depends on an envelope definition that the authors themselves flag as vulnerable to random alignment of ambient CGM, and there are internal numerical inconsistencies that need fixing.\n\nWhat's new here is the population study. Paper I did one galaxy; this does 46-47 with a consistent HVC identification pipeline. The mock absorption spectra and the comparison to Fox et al. (2019) are a genuine attempt to connect simulation to observations. The prolate morphology result (73%) is new, and the scaling relations between accretion rates and SFR are useful. The qualitative statement that ionized envelopes contain roughly 5-6 times the neutral mass is robust to reasonable changes in how you define the envelope.\n\nThe soft spots are real. Section 2.2 defines the accretion-rate envelope as all gas within twice the half-size radius, weighted by max(0,1-dvi), with dvi<=0.5 counted as aligned. The same section says dvi<=0.5 can produce 'unrealistically large ionized envelopes due to random alignment of ambient CGM gas cells'—so they step down the cutoff for the connected 3D envelope morphology. That step-down is not applied to the accretion calculation. The validation of half-size masses against 3D envelopes covers only the fiducial galaxy, and 20% of the clouds disagree by more than 15%. This could inflate both the 5.6x mass ratio and the 81% peak. This is checkable and should be checked.\n\nSecond, the cooling and condensation step is explicitly deferred, so 'sustain star formation' is conditional. That's honestly stated, so it's not a fatal flaw, but the headline should be framed as 'would be sufficient if the gas cools,' which the abstract does say.\n\nThird, there are numerical inconsistencies. The text in Section 3.3 writes the fit as Mneut = 3.16 x Mion + 2.02, but Figure 7 shows Mion = 3.16 x Mneut + 2.02. The fit slopes in the center panel imply ~94% and ~18% of SFR for the average galaxy, while the right panel's distribution peaks at 81% and 11%. The abstract and text swing between 81% and 94%. These need to be reconciled.\n\nWho should read this: anyone working on CGM accretion, HVCs, or mock absorption studies. It's a serious paper with a reproducible method on a public simulation. The right response is to send it to peer review and ask the authors to test the sensitivity of the accretion rates to the envelope definition and to fix the internal numbers. I'd engage with it.","headline":"Qualitatively sound population study of ionized HVC envelopes, but the headline 81% needs firmer envelope definitions and internal number fixes before it can be trusted.","tokens_in":19281,"tokens_out":4418,"would_cite":true,"duration_ms":43248,"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":"In simulations of Milky Way-like galaxies, the warm and hot ionized envelopes around high-velocity clouds carry enough inflowing gas to sustain about 81% of star formation, while neutral clouds supply only 11%.","keywords":["high-velocity clouds","circumgalactic medium","ionized gas","TNG50","IllustrisTNG","gas accretion","star formation","Milky Way analogs"],"falsifier":"Trace the identified ionized envelopes forward in a simulation that resolves radiative cooling, turbulent mixing, and thermal conduction. If fewer than roughly half of the envelope cells reach the star-forming disk as cold gas, the claim that ionized HVCs can sustain 81% of the star formation rate fails. Observationally, a robust census of ionized HVC distances and inflow velocities in the Milky Way yielding a total ionized inflow well below the Galactic star formation rate would also contradict it.","tokens_in":18258,"feed_emoji":"🌌","tokens_out":5405,"duration_ms":55643,"temperature":0.7,"pith_summary":"This paper argues that the ionized gas surrounding high-velocity clouds (HVCs) is the main carrier of gas accretion onto Milky Way-like galaxies, not the neutral gas seen in 21-cm emission. Using 47 galaxies from the TNG50 simulation, it finds that the ionized envelopes contain about 5.6 times more mass than the neutral clouds, and that their infall rate can balance 81% of a galaxy's star formation rate, while neutral HVCs supply only 11%. The authors conclude that the diffuse, ionized circumgalactic medium is what actually sustains star formation at low redshift. The result matters because it reframes the dominant fuel source for galaxy growth as the 'invisible' ionized component rather than the easily observed neutral component.","feed_headline":"Ionized HVC envelopes can fuel 81% of star formation","feed_subtitle":"Simulations of 47 Milky Way-like galaxies show the warm-hot envelopes of high-velocity clouds dominate infalling gas.","key_machinery":"The central tool is a velocity-alignment criterion. For each gas cell, the normalized velocity difference dvi = |v_i - v_HVC|/|v_HVC| measures whether the cell is comoving with the HVC. An alignment fraction f(r) is fit with a cored exponential profile to define a 'half-size' radius where f = 0.5, and the ionized envelope is taken as gas with temperature above $10^{4}$.5 K, weighted by max(0, 1 - dvi), within twice the half-size. Accretion rates sum cell masses times the radial velocity divided by radius, with the ionized calculation using alignment weights to avoid double counting. This machinery lets the paper move from neutral H I emission analogs to the otherwise invisible ionized gas in three dimensions.","core_discovery":"In TNG50 Milky Way-like galaxies, every neutral HVC is surrounded by a much more massive warm and hot ionized envelope that is comoving with it, and this envelope, not the neutral cloud, dominates the accreted mass. Across 47 galaxies, the ionized infall rate scales with the star formation rate with slope 0.94, and the distribution of ionized accretion divided by star formation rate peaks at 81%, versus 11% for neutral HVCs. The paper claims that the ionized high-velocity circumgalactic medium can sustain the observed star formation, provided the material can cool and condense into star-forming gas.","pith_inferences":["If the 81% result holds, then the rate at which ionized envelope gas cools and condenses is the true bottleneck for star formation, so simulations resolving turbulent mixing layers should show envelope gas converting to neutral and molecular phases within a dynamical time.","The velocity-alignment technique could be applied to other multiphase circumgalactic structures, such as satellite streams and thermally unstable gas, to test whether similar ionized envelopes are a universal feature of accreting gas.","A testable extension is the prediction that envelope size decreases with galactocentric radius because of increasing shear and pressure; higher-resolution circumgalactic simulations or deep H I plus ultraviolet observations could check this trend.","Because the paper excludes Magellanic-Stream-like interactions, including such major satellite interactions would likely raise the ionized inflow budget further, potentially overshooting the star formation rate balance and indicating that some ionized gas must be ejected or fail to condense."],"forward_implications":["Ionized envelopes are about six times more massive than neutral HVCs and are prolate in 73% of cases, linking them to filamentary infall rather than spherical accretion.","Mock absorption spectra through simulated sightlines show multiphase gas at HVC velocities, so HVCs should be detectable in ultraviolet ions such as C IV, Si III, and O VI even when H I emission is too faint.","Neutral HVC accretion alone accounts for roughly 11% of the star formation rate, implying that H I surveys substantially underestimate the mass accreting onto galaxies.","The linear scaling of ionized accretion rate with star formation rate indicates that low-redshift galaxy growth is tightly connected to diffuse circumgalactic medium inflow.","Combining neutral and ionized HVC accretion gives a total that peaks at 93% of the star formation rate, nearly closing the fuel budget for star formation.","The ionized envelopes are generally hotter than the ambient circumgalactic medium, suggesting that their high anomalous velocities generate shear heating."],"supporting_citations":[{"why":"Paper I; defines the HVC analog selection by cold gas clouds with velocity deviating more than 70 km/s from galactic rotation, which this paper builds on.","marker":"Lucchini et al. 2024"},{"why":"Presents the TNG50 simulation whose Milky Way-like galaxies are analyzed here.","marker":"Nelson et al. 2019b"},{"why":"Describes the IllustrisTNG project and TNG50 data products used for the galaxy sample and simulation outputs.","marker":"Pillepich et al. 2019"},{"why":"Provides the observational covering fractions and mass flow rates that the mock spectra are compared against.","marker":"Fox et al. 2019"},{"why":"Establishes that ionized HVCs with N(H II) greater than N(H I) exist within 15 kpc, the observational anchor for ionized mass dominance.","marker":"Lehner & Howk 2011"},{"why":"Supplies the sample of Milky Way-like galaxies with criteria on total mass, star formation rate, and diskiness.","marker":"Semenov et al. 2024"},{"why":"Gives the neutral-only HVC inflow estimate that motivates the need for an ionized component to explain galactic accretion.","marker":"Putman et al. 2012"}],"fun_headline_variants":["Ionized HVC envelopes outweigh neutral gas 6-to-1","Invisible ionized HVCs may supply 81% of star fuel","Warm ionized halos dominate high-velocity cloud accretion","TNG50 reveals ionized HVC gas scales with star formation","Ionized envelopes of HVCs could sustain Milky Way growth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the ionized envelope gas, once it falls inward, can actually cool and condense into star-forming material; the paper measures inflow rates but does not simulate whether the clouds survive, cool, or mix before reaching the disk.","fun_headline_variants_meta":{"raw":{"variants":["Ionized HVC envelopes outweigh neutral gas 6-to-1","Invisible ionized HVCs may supply 81% of star fuel","Warm ionized halos dominate high-velocity cloud accretion","TNG50 reveals ionized HVC gas scales with star formation","Ionized envelopes of HVCs could sustain Milky Way growth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000255,"raw_usage":{"total_tokens":1593,"prompt_tokens":986,"completion_tokens":607,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":514}},"tokens_in":602,"tokens_out":607,"duration_ms":7191,"temperature":1.0,"reasoning_tokens":514,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:30:46.748103+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Trace the identified ionized envelopes forward in a simulation that resolves radiative cooling, turbulent mixing, and thermal conduction. If fewer than roughly half of the envelope cells reach the star-forming disk as cold gas, the claim that ionized HVCs can sustain 81% of the star formation rate fails. Observationally, a robust census of ionized HVC distances and inflow velocities in the Milky Way yielding a total ionized inflow well below the Galactic star formation rate would also contradict it.","supporting_citations":[],"review_version":1}