{"id":"afc757b5-eaed-4210-a43c-520b0725cbaf","arxiv_id":"1909.02329","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"In a Bootes I-like ultrafaint dwarf simulated at 0.95 pc resolution, stellar winds and supernovae expel most newly made metals but retain over 80% of the cold gas, so stellar feedback alone cannot quench these galaxies.","lead":"This study uses high-resolution 3D simulations to show that supernova feedback alone cannot blow away the cold gas of an ultrafaint dwarf galaxy like Bootes I. Instead, only the hot metal-rich ejecta escapes, pointing to environmental effects as the likely cause of gas removal.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'unavoidable' environmental conclusion overreaches: only thermal winds plus SNe are modeled, while radiation pressure (with a momentum budget comparable to the binding momentum) is omitted, so 'stellar feedback alone' is not actually tested.","rationale":"The paper's simulation result is internally consistent: under the stated assumptions (smooth hydrostatic gas, continuous thermal injection from static OB associations), the cold gas is retained and the metals escape. The problem is the generalization. The central claim, as phrased in the abstract and conclusions, is that stellar feedback alone cannot quench UFDs and that environmental effects are unavoidable. But the model only includes thermal feedback from winds and SNe. Radiation pressure, photoionization, and cosmic rays are also products of massive stars and are known to affect gas in dwarf galaxies. A simple estimate shows that the cumulative radiation momentum is of the same order as the gas binding momentum, so this omission is not obviously negligible. Thus the conclusion overreaches the model. This is an internal scope mismatch rather than a disagreement with consensus. The reader's weakest_assumption correctly notes the idealized initial conditions, but the more decisive issue is the definition of 'stellar feedback' in the conclusion. The proposed test, adding a radiation-pressure source term to the fiducial run, would settle whether the omitted channel changes the gas-loss fraction materially. If it does, the conclusion needs to be narrowed to 'thermal feedback alone' or supplemented with a treatment of other channels. The verdict remains CONDITIONAL because the underlying simulation is a legitimate, though idealized, study; the condition is that the conclusions must be reworded to match the modeled physics, or the missing physics must be shown to be negligible.","tokens_in":20557,"tokens_out":13729,"duration_ms":151321,"concrete_test":"Rerun the fiducial high-resolution adiabatic simulation (A-10-04-01-P) with an added radiation-pressure term: inject momentum at rate L/c into the same ten 4-pc OB-association spheres, using the Leitherer et al. (2014) luminosity of the 10^5 M_sun population over 30 Myr, keeping all other physics identical. Compare the gas mass retained inside the 2-kpc box at t = 30 Myr with the ~80% shown in Fig. 6. If the retained fraction drops below ~60%, the central claim that stellar feedback alone cannot remove the gas is not supported; if it stays above ~80%, the omission is quantitatively unimportant for this conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 defines the feedback as thermal only: 'The massive star feedback is modelled through thermal energy deposition, and no other mechanism (such as, for instance, radiation pressure) is included.' The abstract and Section 5, however, generalize this to 'stellar feedback alone' being unable to remove the gas and conclude that environmental effects are 'unavoidable.' That inference is not supported by the model because a known stellar-feedback channel is absent. A momentum budget for the adopted 10^5 M_sun population using Leitherer et al. (2014) luminosities gives L ~ 10^42 erg/s, so over 30 Myr the radiation momentum is L t / c ~ 3e46 g cm/s, comparable to M_gas v_esc ~ 1.2e40 g x 3e6 cm/s ~ 4e46 g cm/s. Radiation pressure alone could therefore unbind a substantial fraction of the gas, and similar considerations apply to photoionization and cosmic rays. The simulations establish that thermal feedback from winds and SNe is insufficient; they do not establish that stellar feedback as a whole is insufficient. The 'unavoidable' wording should be conditioned on the omitted mechanisms being negligible.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents high-resolution AMR simulations of an isolated, non-rotating Boötes I-like ultrafaint dwarf galaxy, initialized as a smooth, single-phase gas distribution in hydrostatic equilibrium with a static Burkert dark-matter halo. A coeval 10^5 M_sun stellar population, represented by ten static OB associations, injects mass and thermal energy from stellar winds and supernovae for 30 Myr. The authors run one high-resolution adiabatic run to 30 Myr and one radiative run to 20 Myr, plus lower-resolution variants with different numbers and positions of OB associations. They find that less than 20-30 per cent of the initial gas is lost, no complete blow-away occurs, and the bulk of the cold gas remains bound, whereas 70-80 per cent of the SN ejecta is lost from the box by 30 Myr in the adiabatic run. They contrast this with the analytic estimate that about 50 SNe should suffice and conclude that stellar feedback alone cannot remove the gas and that environmental mechanisms are unavoidable.","tokens_in":20811,"tokens_out":15236,"duration_ms":151572,"significance":"If the result holds for the modeled feedback channels, it is a useful counterpoint to simple energy-coupling arguments and supports previous hydrodynamical studies that find SN feedback in low-mass haloes to be inefficient at driving global outflows. The paper's strengths are the ~1 pc resolution in the high-resolution runs, the explicit resolution and configuration checks in Appendix A, the separate tracking of gas and metals, and the direct comparison between a binding-energy estimate and the simulated outcome. However, because the model includes only thermal energy deposition from winds and SNe, the significance for the broader question of 'stellar feedback alone' is limited; radiation pressure, photoionization and cosmic rays are not tested, and the paper's own limitation statements in Sections 4.3-4.4 and Appendix A acknowledge several of the missing ingredients.","major_comments":[{"comment":"The simulations model stellar feedback only as thermal energy and mass deposition from winds and SNe: §2 states that 'no other mechanism (such as, for instance, radiation pressure) is included,' and §4.4 lists radiation, cosmic rays, and magnetic fields as missing ingredients. The abstract and §5, however, conclude that 'stellar feedback alone' cannot remove the gas and that environmental effects are 'unavoidable.' That inference is not supported by the model. A simple momentum estimate for the adopted 10^5 M_sun population, using L ~ 10^42 erg/s from Leitherer et al. (2014), gives L t / c ~ 3 x 10^46 g cm/s over 30 Myr, comparable to the binding momentum M_gas v_esc ~ 4 x 10^46 g cm/s. Radiation pressure could therefore contribute to unbinding a substantial fraction of the gas, and similar order-of-magnitude considerations apply to photoionization and cosmic rays. I request that the conclusions be restricted to thermal feedback from winds and SNe, or that the model be extended, or that a quantitative argument be given for why the omitted channels are negligible.","section":"Abstract; §2; §5"},{"comment":"The quantitative headline numbers (20-30 per cent gas loss and 70-80 per cent ejecta loss by 30 Myr) come from a single high-resolution realization for the fiducial configuration (A-10-04-01-P and R-10-04-01-P in Table A.1), while the lower-resolution suite explores variations in the number and placement of OB associations. The paper itself states that the exact gas-loss fraction depends more on association location than on number, and Fig. A.3 shows appreciable run-to-run spread. I would like the abstract and the §5 bullet points to present these ranges as illustrative of the sampled configurations rather than as robust predictions, or the authors to add at least one additional high-resolution realization to quantify the scatter. The qualitative conclusion that no complete blow-away occurs is supported by the existing suite; this comment concerns the precision of the quantitative claims.","section":"§3 and Appendix A"}],"minor_comments":[{"comment":"The definition of M0 in Eq. (4) is confusing: M0 is used in Eq. (3) as a mass in units of M_sun, but writing 'M0 = MDM/5.8 M_sun' mixes a dimensionless ratio with the M_sun unit. Please clarify the intended dimensional statement.","section":"§4.1, Eqs. (3)-(4)"},{"comment":"There is a typo in the sentence comparing the adiabatic and cooling metallicity profiles: 'predited' should be 'predicted'.","section":"§3.3"},{"comment":"Several typos occur in these sections: 'feeback' should be 'feedback', 'simultations' should be 'simulations', and the reference 'V orobyov et al. 2015' should have the correct spacing 'Vorobyov et al. 2015'.","section":"§4.4 and §5"},{"comment":"The x-axis labels of Fig. A.2 appear incomplete: bins such as 10-100 pc, 200-300 pc, 500-600 pc, and 600-700 pc are not visible in the printed labels. Please ensure the full bin ranges are shown or described in the caption.","section":"Fig. A.2"},{"comment":"The abstract does not mention that the radiative high-resolution run is truncated at 20 Myr rather than 30 Myr; given that the conclusions are drawn partly from the 30 Myr adiabatic run, a brief statement of the radiative run's duration would help readers assess the comparison.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is carefully written and the simulations are described in detail; the main issue is the gap between the thermal-only feedback modeled in Section 2 and the broad 'stellar feedback alone' and 'unavoidable environmental effects' language in the abstract and Section 5. If the authors are willing to condition those conclusions on the modeled channels and to frame the quantitative loss fractions as configuration-dependent, the paper would be a solid contribution to the series. The single-realization issue at high resolution is secondary but should be acknowledged in the quantitative claims. No concerns about novelty or scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis one is worth reading if you work on ultrafaint dwarfs, but the headline claim overreaches a bit. What the paper does well: it runs very high-resolution (0.95 pc) AMR simulations of a Bootes I-like UFD with thermal feedback from winds and SNe, and shows—robustly, I think—that this feedback does not expel the bulk of the cold gas, even in the adiabatic limit. The separation between retained cold gas and escaping metal-rich ejecta is a genuinely useful quantitative result, and the paper is unusually honest about its limitations: the radiative run stops at 20 Myr, only one random OB association configuration is done at high resolution, and the ISM is smooth and single-phase.\n\nThe soft spot is real: the paper only includes thermal energy deposition. It explicitly says radiation pressure is not modeled. Yet the abstract and conclusions say 'stellar feedback alone' cannot remove the gas and that environmental effects are 'unavoidable.' That is a step too far. A quick momentum estimate—L t / c for a 10^5 Msun population over 30 Myr is within a factor of a few of M_gas v_esc—suggests radiation pressure could matter, and photoionization and cosmic rays are also left out. So the paper actually demonstrates that thermal feedback from winds and SNe is insufficient, not that stellar feedback in general is. That needs rewording, and probably some discussion of why radiation pressure would or wouldn't help in a low-metallicity UFD.\n\nOther soft spots: the high-res result is a single realization, so we don't know the stochastic spread; the appendix tests different OB association counts and concentrations at lower resolution, but not the full random distribution at high res. And no code or data are public, so independent checking is limited. These are reproducibility niggles, not signs of error.\n\nIn sum: the core result is solid and the paper is a useful addition. It deserves a serious referee, and with a toned-down conclusion it would be acceptable as is. I'd send it to review if I were the editor, but I'd ask the authors to fix the overreach and add a sentence or two on the missing feedback channels.\n\nRecommendation: engage with it, cite it if you work on dwarf feedback, and if you review it, push on the radiation-pressure caveat.","headline":"Solid high-resolution simulation showing thermal feedback alone can't remove cold gas from a UFD, but the paper overreaches by concluding stellar feedback in general is insufficient when it never modeled radiation pressure.","tokens_in":21372,"tokens_out":2850,"would_cite":false,"duration_ms":29653,"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":"Stellar feedback alone cannot remove the cold gas from an ultrafaint dwarf galaxy like Boötes I, even when the total supernova energy exceeds the gas binding energy.","keywords":["ultrafaint dwarf galaxies","stellar feedback","supernova-driven outflows","gas removal","hydrodynamic simulations","Boötes I","galactic winds","chemical enrichment"],"falsifier":"Run the same Boötes I-like galaxy with a clumpy, multiphase initial gas distribution and OB associations that move with the flow, keeping the total injected energy fixed; if more than half of the cold gas leaves the simulation box within 30 Myr, the claim that stellar feedback alone cannot remove the bulk of the gas would be falsified.","tokens_in":20329,"feed_emoji":"💥","tokens_out":7524,"duration_ms":77800,"temperature":0.7,"pith_summary":"This paper asks whether stellar feedback from a burst of star formation can, by itself, blow the cold gas out of an ultrafaint dwarf galaxy like Boötes I. Using high-resolution three-dimensional hydrodynamic simulations with an instantaneous $10^5\\,M_\\odot$ stellar population injecting winds and supernova energy for 30 Myr, it finds the answer is no: even though the total supernova energy exceeds the gas binding energy, no galactic-scale outflow develops and less than 20–30 percent of the cold gas is lost. The same simulations find the opposite for metal-rich ejecta: by 30 Myr, 70–80 percent of the supernova ejecta has left the simulation box. If this is right, the gas-free state of ultrafaint dwarfs cannot be explained by internal feedback alone, and environmental mechanisms such as ram-pressure stripping or tidal interactions must finish the cleaning. This matters because it challenges the common analytic assumption that a galaxy's gas is ejected once cumulative supernova energy exceeds binding energy.","feed_headline":"Supernova energy alone cannot sweep out an ultrafaint dwarf","feed_subtitle":"Simulations of Boötes I keep over 70% of cold gas while leaking metals; environment must finish the job.","key_machinery":"The load-bearing machinery is the spatially resolved, thermal injection of stellar winds and supernova energy from ten fixed OB associations inside a high-resolution, adaptive-mesh grid, with the galaxy initialized as a smooth, single-phase, non-rotating gas in hydrostatic equilibrium inside a static dark-matter potential. The comparison that carries the argument is between the integrated energy budget and the actual coupling: only a small filling factor of hot superbubbles forms, shocks from inner associations compress the central gas, and energy vents along low-density chimneys instead of accelerating the cold gas coherently. A passive metallicity tracer follows the separate fate of the enriched ejecta, showing that it leaves through the same channels while the ambient medium stays bound.","core_discovery":"On the paper's own terms, the central discovery is that stellar feedback in an ultrafaint dwarf is a selective leak, not a wind. With ten OB associations randomly scattered through a smooth, single-phase gas distribution in a dark-matter-dominated halo modeled on Boötes I, supernovae carve hot bubbles and chimneys through which metal-enriched ejecta escapes, while the cold ambient gas remains trapped. In the adiabatic run, less than 20–30 percent of the initial gas mass is lost after 30 Myr; in the run with radiative cooling, the cold gas is even less affected, and the final density profile is nearly indistinguishable from the initial one. Meanwhile, by 30 Myr the adiabatic simulation has lost 70–80 percent of the supernova ejecta, with roughly half the ejecta already gone by 20 Myr in the cooling run. The authors conclude that the simple energetic criterion—enough supernova energy to exceed the gas binding energy—does not guarantee gas removal, because the injected energy is channeled along low-density paths rather than coupled to the bulk of the gas.","pith_inferences":["Inference: a two-stage cleaning—stellar feedback leaks the metals, then an environmental process strips the remaining cold gas—would naturally produce the old, very metal-poor stellar populations of ultrafaint dwarfs, though the paper does not simulate that sequence.","Inference: because the authors use fixed, non-moving OB associations and note that this maximizes feedback effectiveness, allowing the sources to move with the flow would likely spread and dilute the energy deposition and make gas removal even harder.","Inference: a testable extension is to repeat the runs with a clumpy, multiphase initial interstellar medium; the likely outcome is that cold clumps survive while hot ejecta vents even more efficiently, preserving the qualitative conclusion while changing the quoted percentages.","Inference: the 30 Myr window covers one instantaneous burst; a longer, lower-intensity star-formation history consistent with the integrated-galactic-initial-mass-function argument would produce fewer massive stars per unit time, so the feedback failure would likely be more severe, not less."],"forward_implications":["The analytic rule of thumb that gas is ejected once cumulative supernova energy exceeds binding energy is not a reliable predictor for ultrafaint dwarfs; simulations must resolve where the energy actually goes.","If stellar feedback cannot remove the cold gas, explaining gas-free ultrafaint dwarfs requires external agents such as ram-pressure stripping, tidal interactions, or reionization to do the work, which the paper identifies as the next step.","Gas loss and metal loss decouple: a dwarf can shed most of its freshly produced metals while retaining most of its original cold gas, so chemical-evolution models that tie the two together need revision.","Even in the adiabatic limit, with radiative cooling switched off and an instantaneous coeval starburst, no galactic-scale outflow develops, so the failure of feedback is not a cooling artifact."],"supporting_citations":[{"why":"Supplies the mass and energy injection rates for stellar winds and supernovae that drive the feedback in the simulations.","marker":"Leitherer et al. (2014)"},{"why":"Provides the resolution criterion showing that the bubble cooling radius is resolved so injected momentum is captured correctly.","marker":"Kim & Ostriker (2015)"},{"why":"Justifies treating discrete supernovae in an OB association as a continuous energy source once the blast waves become subsonic.","marker":"Mac Low & McCray (1988)"},{"why":"Supplies the explanation that scattered explosion sites vent energy through low-density channels and compress central gas, preventing blow-away.","marker":"Mori et al. (2002)"},{"why":"Earlier three-dimensional hydrodynamical study of single supernovae in ultrafaint-sized haloes that this work builds upon.","marker":"Bland-Hawthorn et al. (2015)"},{"why":"Previous simulations of isolated low-mass haloes with extended star formation whose conclusions about ultrafaint dwarf formation are compared with the present results.","marker":"Webster et al. (2015)"},{"why":"Provides the OB-association sampling procedure and the analogous hydrodynamic setup whose concentrated sources drive a wind, the contrast case for this paper.","marker":"Calura et al. (2015)"},{"why":"Supports the claim that at high resolution stellar feedback results are insensitive to subgrid cooling details.","marker":"Read et al. (2016)"}],"fun_headline_variants":["Supernovae leak metals, trap gas in ultrafaint dwarfs","Stellar feedback in ultrafaint dwarfs: leak, not wind","Gas stays, metals go: supernova feedback in ultrafaint dwarfs","Supernova energy leaks away: gas retains, metals escape","Selective leak: supernovae strip metals, not gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations assume the gas begins as a smooth, single-phase, non-rotating cloud in hydrostatic equilibrium, with all star-forming associations fixed in place for the whole 30 Myr; if the real interstellar medium was clumpy or multiphase, or the sources moved, the retained-gas and metal-loss fractions could change.","fun_headline_variants_meta":{"raw":{"variants":["Supernovae leak metals, trap gas in ultrafaint dwarfs","Stellar feedback in ultrafaint dwarfs: leak, not wind","Gas stays, metals go: supernova feedback in ultrafaint dwarfs","Supernova energy leaks away: gas retains, metals escape","Selective leak: supernovae strip metals, not gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000937,"raw_usage":{"total_tokens":4079,"prompt_tokens":1089,"completion_tokens":2990,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2911}},"tokens_in":705,"tokens_out":2990,"duration_ms":21030,"temperature":1.0,"reasoning_tokens":2911,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:53:27.677738+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Boötes I-like galaxy with a clumpy, multiphase initial gas distribution and OB associations that move with the flow, keeping the total injected energy fixed; if more than half of the cold gas leaves the simulation box within 30 Myr, the claim that stellar feedback alone cannot remove the bulk of the gas would be falsified.","supporting_citations":[{"cited_title":"2014, , 212, 14","cited_arxiv_id":null,"evidence_quote":"Supplies the mass and energy injection rates for stellar winds and supernovae that drive the feedback in the simulations."},{"cited_title":"& Ostriker , E","cited_arxiv_id":null,"evidence_quote":"Provides the resolution criterion showing that the bubble cooling radius is resolved so injected momentum is captured correctly."},{"cited_title":"2002, , 571, 40","cited_arxiv_id":null,"evidence_quote":"Supplies the explanation that scattered explosion sites vent energy through low-density channels and compress central gas, preventing blow-away."},{"cited_title":"2015, , 807, 154","cited_arxiv_id":null,"evidence_quote":"Earlier three-dimensional hydrodynamical study of single supernovae in ultrafaint-sized haloes that this work builds upon."},{"cited_title":"2015, , 799, L21","cited_arxiv_id":null,"evidence_quote":"Previous simulations of isolated low-mass haloes with extended star formation whose conclusions about ultrafaint dwarf formation are compared with the present results."},{"cited_title":"G., Romano , D., & D'Ercole , A","cited_arxiv_id":null,"evidence_quote":"Provides the OB-association sampling procedure and the analogous hydrodynamic setup whose concentrated sources drive a wind, the contrast case for this paper."},{"cited_title":"I., Agertz , O., & Collins , M","cited_arxiv_id":null,"evidence_quote":"Supports the claim that at high resolution stellar feedback results are insensitive to subgrid cooling details."}],"review_version":1}