{"id":"1035e2bf-e102-4edc-99ee-b11be2486697","arxiv_id":"1909.00815","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Supernovae-heated superbubbles with high entropy are buoyantly accelerated through the circumgalactic medium, producing slow, long-lived outflows that match observed CGM kinematics and simulation recycling times.","lead":"This paper proposes that hot bubbles from supernova explosions rise through a galaxy's gaseous halo by buoyancy, not by initial momentum, and can remain there for over a billion years despite moving slowly. If correct, it changes how galactic winds, gas recycling, and star formation regulation are modeled and explains puzzling slow outflows seen in observations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Long Gyr recycling times hinge on bubbles evolving as coherent constant-entropy parcels; realistic entrainment and mixing would erase the overshoot and oscillations shown in Fig 3, so the model's central claim is not yet established.","rationale":"The paper's analytic framework is internally consistent and the buoyancy mechanism is physically plausible; the equations are explicit and the Hill et al. comparison provides independent support for the production of high-entropy gas. The difficulty is that the quantitative headline result, multi-Gyr recycling at ~100 km/s, requires the bubble to remain a closed, non-mixing object. The reader's weakest assumption correctly identifies this. My own reading sharpens it: entrainment is not just a possible instability but an inevitable feature of buoyant thermals in stratified fluids. The conservative oscillator in Eq 21-24 has no entropic mixing channel; any real mixing converts the predicted oscillations into a single damped plume rise, reducing both the maximum height and the residence time. The paper's own caveats in Sections 3.3 and 6.5 concede this. Existing evidence, such as AGN bubbles in clusters, shows survival to large radii in a much lower-density medium, so it does not resolve the L* CGM regime. A targeted 3D simulation with the paper's own initial conditions would settle whether the long-lived coherent bubble exists. Until that check, the central claim should remain conditional rather than accepted.","tokens_in":37486,"tokens_out":13259,"duration_ms":141662,"concrete_test":"Run a 3D hydro simulation of a single pressure-equilibrated bubble with the Fig 3 initial conditions (KSB=5.84 keV cm^2, mSB from Eq 10, h=200 pc, namb=1 cm^-3) rising through a hydrostatic CGM with K(r)=K200(r/R200)^1.1 and the density normalization of Eq 29, including radiative cooling and a resolution of <10 pc. Track the maximum radius reached by the connected region whose mass-weighted entropy remains within a factor of 2 of KSB, and measure the time until that entropy contrast is lost. If the coherent region fails to reach the Eq 24 neutral-buoyancy radius, or mixes before 1 Gyr, then the long recycling times claimed in Fig 3 are not robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim, that slow superbubbles persist in the CGM for >Gyr, is carried by Eq 21-24 and the oscillatory trajectories of Fig 3. Those equations model a single Lagrangian bubble with fixed entropy KSB and fixed mass mSB, with drag computed for a rigid sphere (Eq 25). This is exactly the regime in which a buoyant blob conserves its entropy contrast and oscillates about the neutral-buoyancy radius. A real buoyant thermal in a stratified medium entrains ambient gas through Kelvin-Helmholtz and Rayleigh-Taylor vorticity at its surface; the entropy contrast is diluted, the effective cross-section grows, and the motion becomes a one-shot plume rise rather than a damped oscillator. The authors explicitly concede in Section 3.3 that long-term oscillations \"will not occur if the bubble is mixed into the CGM,\" and in Section 6.5 that the coherent-bubble treatment is a rough approximation. Because the >Gyr recycling times and high apocenters in Fig 3 are the product of the conservative oscillator, the coherence assumption is load-bearing, not a detail. The cited cluster observations show bubbles survive to ~100 kpc in a much more rarified, magnetized ICM; they do not establish Gyr survival in an L* CGM. The MUGS2 comparison is suggestive but partly circular, since those runs use the authors' own superbubble feedback model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes that supernova-driven superbubbles with entropy exceeding the local CGM entropy are buoyantly accelerated after breaking out of the disc, producing outflows and fountains that reach high galactocentric radii and re-accrete on >Gyr timescales despite modest velocities. The authors derive an equation of motion (Eq. 32) combining buoyancy, gravity, and drag for a spherical bubble in a hydrostatic, power-law CGM, with the superbubble initial entropy from Weaver et al. (1977) and Mac Low & McCray (1988) (Eq. 15). Numerical integrations show bubbles reaching ~100 kpc, overshooting their neutral-buoyancy radius, and oscillating with damping by drag. The framework is compared with MUGS2 cosmological simulations, Hill et al. (2018) ISM simulations, and COS-Halos/Stocke et al. observations, and the authors report qualitative and quantitative agreement.","tokens_in":37702,"tokens_out":4470,"duration_ms":49518,"significance":"If correct, this framework resolves a central tension in galactic outflow studies: observed CGM kinematics are slow (typically ~100-200 km/s) while recycling times inferred from simulations and metal distributions require >Gyr residence. The analytic derivation is a clear strength: it is not fitted to the observations it explains, it builds on classical buoyancy arguments, and it makes falsifiable predictions for phase-space distributions of CGM absorbers. The use of public high-resolution ISM simulations to test the breakout entropy normalization (Eq. 15) is particularly valuable. The main limitation is that the quantitative central claim of long recycling times and high apocenters rests on the assumption that bubbles remain coherent, constant-entropy objects for Gyr; the authors themselves flag this in Sec. 3.3 and Sec. 6.5. Because this load-bearing assumption is acknowledged but not quantitatively tested, the paper needs revision rather than acceptance as-is.","major_comments":[{"comment":"The central result, that slow superbubbles persist in the CGM for >Gyr, is produced by treating the bubble as a single Lagrangian parcel with fixed entropy KSB and fixed mass mSB, with drag computed for a rigid sphere (Eq. 25). Under this assumption the bubble conserves its entropy contrast and oscillates about the neutral-buoyancy radius, giving the long recycling times shown in Fig. 3. Real buoyant thermals in a stratified medium entrain ambient gas through Kelvin-Helmholtz and Rayleigh-Taylor instabilities, which dilutes the entropy contrast and increases the effective cross-section; the motion becomes a one-shot plume rise rather than a damped oscillator. The authors explicitly concede in Sec. 3.3 that long-term oscillations 'will not occur if the bubble is mixed into the CGM' and in Sec. 6.5 that the coherent-bubble treatment is a rough approximation. Because the >Gyr recycling times and high apocenters in Figs. 3, 6, and 7 are products of this conservative oscillator, the coherence assumption is load-bearing. The cited cluster observations (Sec. 3.3) show bubbles surviving to ~100 kpc in a much more rarified, magnetized ICM, and do not establish Gyr survival in an L* CGM. To support the central claim, the authors should provide a quantitative entrainment/mixing model or high-resolution 3D simulations of bubble rise in an L* CGM demonstrating that the coherent-bubble regime is realized.","section":"Sec. 2.3-2.4, Eqs. (21)-(25), Fig. 3"},{"comment":"The cooling analysis is decoupled from the dynamics: trajectories are integrated with constant KSB, and the cooling time is then compared with the flight time. If the bubble loses entropy gradually, its buoyancy decreases continuously, and the trajectory, apocenter, and recycling time change before the gas reaches the quoted cooling time. The binary distinction between 'adiabatic' and 'cooled' used in Figs. 7-9 is therefore not sufficient to establish the cooling-limited heights in Fig. 9. A coupled treatment, or a demonstration that cooling does not appreciably alter the trajectories for the cases shown, is needed to support the claim that entropy-driven winds can reach ~100 kpc before radiative losses become important.","section":"Sec. 3.1, Figs. 5, 8, 9"},{"comment":"The comparison to MUGS2 is suggestive but not an independent test of the model. MUGS2 uses the Keller et al. (2014) superbubble feedback implementation, which is based on the same Weaver et al. (1977) and Mac Low & McCray (1988) evaporation physics used to derive Eq. (15). The correlation between initial entropy and re-accretion time shown in Fig. 13 is therefore expected in part by construction. The authors should state this limitation explicitly and, ideally, compare with simulations using a different feedback implementation to break the circularity.","section":"Sec. 4.1, Fig. 13"}],"minor_comments":[{"comment":"There are several typographical errors: 'viral entropy' should be 'virial entropy' in Sec. 5 and elsewhere, and 'As we will will see later' appears in Sec. 3.2.","section":"Throughout"},{"comment":"The axis label 'Height of Blob (kpc)' is unusual; 'Galactocentric radius' or 'Height above disc plane' would be clearer.","section":"Fig. 3"},{"comment":"The sentence 'This height is always the first turnover point in the flight (as subsequent oscillations are damped by drag), except in the case where the apoapsis time is > 10 Gyr' is slightly ambiguous and should be rephrased to clarify whether the maximum height is defined by the first turning point or by the end of the integration.","section":"Sec. 3.2"},{"comment":"The discussion of entrainment and multiphase gas would benefit from a more concrete statement of which of the model's predictions survive if the bubble fragments; currently the text says mixing 'will not completely halt' the outflow, but the quantitative impact on recycling times is not assessed.","section":"Sec. 3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and the core idea is interesting and likely publishable after revision. The main issue is that the headline quantitative claim depends on the coherent-bubble assumption, which the authors themselves concede is rough; the revision should either provide a quantitative justification for this assumption or substantially soften the abstract and conclusion. The MUGS2 comparison should be reframed as a consistency check rather than an independent validation, given the shared feedback physics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The reader's conditional verdict is about right. The core idea is simple: supernova-heated bubbles with entropy above the local CGM entropy keep rising after break-out, so outflows can be slow yet reach far out and recycle slowly. Eq. 15 is a genuinely new closed-form expression for bubble entropy at break-out from ISM density, scale height, and cluster luminosity; Eq. 32 is a clean equation of motion with buoyancy, gravity, and drag. I also give them credit for not hiding the rough edges: Sections 3.3 and 6.5 explicitly concede that coherent, non-mixing bubbles are an idealization and that the Gyr oscillations in Fig. 3 would not occur if the bubble mixes.\n\nThe stress-test note is correct that the long recycling time is carried by the conservative oscillator. If a bubble shreds or entrains, its entropy contrast drops, drag rises, and it becomes a one-shot plume, not a damped oscillator. That would shorten the recycling time and lower the apoapsis. The authors argue that mixing may not halt the outflow, and cited cluster bubbles do rise over 100 kpc, but that is a different regime (rarer gas, higher entropy contrast). So the headline claim of >Gyr residence is not yet established quantitatively. The MUGS2 comparison is also partly circular because the simulations use the same feedback model, though the entropy-vs-lifetime correlation in Fig. 13 is genuinely interesting and independent of the analytic bubble model. The COS-Halos comparison is qualitative; the contours in Fig. 19 are suggestive, not a fit.\n\nNone of this kills the paper. The central physical mechanism, buoyancy as a slow, continuous accelerator, should survive moderate mixing; what gets lost is the precise oscillation timescale. The paper reads as an honest extension of Bower et al. (2017) and Lochhaas et al. (2018), and its parameter-free derivation deserves a serious referee. I would send it to review, with a request that the authors either simulate or estimate a mixing timescale for the bubbles and show how the recycling times degrade when the coherent-bubble assumption is relaxed. That seems like the one revision that would make the strongest claims defensible.","headline":"A clean analytic framework for buoyancy-driven outflows whose headline Gyr recycling times rest on a coherent-bubble assumption the authors themselves flag; worth serious review, but the strongest claims need a mixing timescale.","tokens_in":38298,"tokens_out":2578,"would_cite":true,"duration_ms":28168,"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":"Superbubbles with entropy above the surrounding halo gas keep accelerating upward after leaving the disc, reaching ~100 kpc and recycling in >1 Gyr even when launched slowly.","keywords":["galactic winds","circumgalactic medium","superbubbles","buoyancy","entropy","outflows","fountains","star formation feedback"],"falsifier":"A high-resolution simulation of a superbubble rising through an entropy-stratified halo that shows the bubble fully mixing into the background within ~100 Myr, before it reaches the radius where $K_{\\rm SB}=K(r)$, would falsify the long-recycling prediction. Likewise, an observation of cool, low-entropy CGM gas at ~100 kpc moving at ~100 km/s with no coexisting high-entropy phase would contradict the entropy-driven picture.","tokens_in":37203,"feed_emoji":"🌬️","tokens_out":6557,"duration_ms":56768,"temperature":0.7,"pith_summary":"This paper establishes that the hot bubbles blown by clustered supernovae keep accelerating after they break out of a galaxy's disc, lifted by buoyancy in the entropy-stratified circumgalactic medium (CGM). A bubble whose entropy exceeds the local CGM entropy feels an upward force that drives it to high galactocentric radii and keeps it in the CGM for more than a gigayear, even when its initial velocity is far below the escape velocity. If correct, this resolves the tension between slow observed outflow velocities and the long recycling times implied by simulations and by the presence of metals far outside galaxies.","feed_headline":"Slow galactic winds float for a billion years on buoyancy","feed_subtitle":"Hot superbubbles keep rising through the galaxy's halo, explaining slow outflows and billion-year recycling times.","key_machinery":"The central object is the entropy $K = k_B T n^{1-\\gamma}$, the adiabatic invariant that replaces thermodynamic entropy in astrophysical haloes. The load-bearing identity is the buoyant acceleration $\\ddot{r} = \\nabla\\phi\\,[(K_{\\rm SB}/K(r))^{3/5}-1]$, which converts the classical Archimedean buoyancy force into a statement about entropy contrast. Combined with a power-law CGM entropy profile $K(r)=K_{200}(r/R_{200})^{\\alpha}$, a hydrostatic density profile, and a drag term for a pressure-confined sphere of fixed mass, this yields the full equation of motion that is integrated to produce bubble trajectories. The initial superbubble entropy is derived from the standard luminosity-driven superbubble solution, evaluated at breakout when the bubble radius equals the ISM scale height.","core_discovery":"The central discovery is that a superbubble's fate after breaking out of a star-forming disc is set by its entropy relative to the entropy-stratified CGM, not by its launch speed. The paper derives an equation of motion from Archimedes' principle in entropy form: a bubble with entropy $K_{\\rm SB}$ greater than the local CGM entropy $K(r)$ feels an upward acceleration proportional to $[(K_{\\rm SB}/K(r))^{3/5}-1]$. The bubble rises until it reaches the radius where $K_{\\rm SB}=K(r)$, overshoots due to momentum conservation, and then oscillates about the buoyant equilibrium, damped by drag. For a Milky Way-like halo this trajectory reaches roughly 100 kpc and returns in more than a gigayear while keeping velocities near 100 km/s, well below escape speed. The paper shows that these bubbles carry significant mass loading and that the predicted kinematics match both self-consistent cosmological simulations and CGM absorption-line observations.","pith_inferences":["If the coherent-bubble assumption fails, much of the material might still be transported outward as mixed high-entropy gas, but the sharp >1 Gyr recycling times and clean oscillations would not occur; this is testable with 3D simulations of bubbles rising through stratified haloes.","The same entropy-buoyancy argument naturally extends to AGN-heated bubbles in cluster haloes, where buoyancy is already invoked, and could provide a unified description of supernova- and AGN-driven outflows.","A distinguishing observable is that entropy-driven winds rise as discrete, anisotropic bubbles, whereas global wind models produce smooth shells; spatially resolved absorption-line kinematics could separate the two pictures.","The model predicts a monotonic relation between the entropy of ejected gas and its re-accretion time, so correlating outflow entropy with recycling times in simulations would provide a direct quantitative test."],"forward_implications":["Slow, mass-loaded outflows can persist in the CGM for more than a gigayear, so observed velocities near 100 km/s do not imply rapid re-accretion onto the galaxy.","Entropy-driven fountains recycle gas on gigayear timescales, meaning a gas parcel only needs to be ejected a few times to spend most of its life outside the star-forming disc.","In haloes above roughly $10^{12}\\,M_\\odot$, the virial entropy exceeds typical superbubble entropy, suppressing buoyant uplift and explaining why supernova feedback loses effectiveness at that mass scale.","Bubbles launched from thicker, more diffuse ISM reach higher altitudes and cool more slowly, and lower metallicity extends their cooling times, making the mechanism stronger at high redshift.","The predicted bubble trajectories pass through the temperature range where O VI absorption is strong, directly connecting the model to ultraviolet CGM absorption surveys."],"supporting_citations":[{"why":"Supplies the self-similar luminosity-driven superbubble solution used to compute interior density, temperature, and entropy.","marker":"Weaver et al. (1977)"},{"why":"Provides the superbubble interior density and mass scalings used to derive the breakout entropy K_SB.","marker":"Mac Low & McCray (1988)"},{"why":"Gives the pressure-equilibrium approximation that lets the bubble's density be written in terms of entropy and background pressure.","marker":"Kompaneets (1960)"},{"why":"Establishes the buoyancy criterion for a fluid parcel in an entropy-stratified medium, the physical basis of the acceleration.","marker":"Chandrasekhar (1961)"},{"why":"Provides COS-Halos data on CGM density, temperature, and kinematics that serve as the primary observational comparison.","marker":"Werk et al. (2014)"},{"why":"Prior entropy-threshold argument for feedback effectiveness that this paper extends into a full kinematic model.","marker":"Bower et al. (2017)"},{"why":"High-resolution ISM simulations used to test whether realistic superbubbles generate the predicted high entropies.","marker":"Hill et al. (2018b)"},{"why":"Superbubble feedback model implemented in the cosmological simulations whose fountain recycling times are compared to the predictions.","marker":"Keller et al. (2014)"},{"why":"Provides the analytic framework for entropy-stratified CGM atmospheres and buoyancy in haloes.","marker":"Voit et al. (2017)"}],"fun_headline_variants":["Entropy lifts galactic winds gently for a billion years","Buoyant bubbles drive slow, long-lived galactic winds","Superbubbles float on entropy, explaining slow outflows","Entropy, not speed, controls how galactic winds rise","Hot bubbles rise for gigayears via entropy buoyancy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The rising bubble remains a single coherent, pressure-confined object with constant entropy and mass for gigayear timescales, with no mixing, fragmentation, or disruption; if bubbles break apart, the long recycling times and high apocenters do not occur.","fun_headline_variants_meta":{"raw":{"variants":["Entropy lifts galactic winds gently for a billion years","Buoyant bubbles drive slow, long-lived galactic winds","Superbubbles float on entropy, explaining slow outflows","Entropy, not speed, controls how galactic winds rise","Hot bubbles rise for gigayears via entropy buoyancy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001129,"raw_usage":{"total_tokens":4730,"prompt_tokens":1022,"completion_tokens":3708,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":3627}},"tokens_in":638,"tokens_out":3708,"duration_ms":117479,"temperature":1.0,"reasoning_tokens":3627,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:36:00.447505+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution simulation of a superbubble rising through an entropy-stratified halo that shows the bubble fully mixing into the background within ~100 Myr, before it reaches the radius where $K_{\\rm SB}=K(r)$, would falsify the long-recycling prediction. Likewise, an observation of cool, low-entropy CGM gas at ~100 kpc moving at ~100 km/s with no coexisting high-entropy phase would contradict the entropy-driven picture.","supporting_citations":[],"review_version":1}