{"id":"e33d0e84-dc88-48e4-8f9b-6ab4ce67dc81","arxiv_id":"2506.16624","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Cluster winds from resolved star clusters supply at least 70% of the shock energy and power the galactic superwind in the starburst merger VV 114.","lead":"Astronomers mapped the gas around 66 star clusters in the merging galaxy VV 114 and found that winds from these clusters carry enough energy to drive most of the galaxy's shock-heated gas and a galactic-scale outflow. The result is one of the clearest links yet between individual star clusters and galaxy-wide feedback in a starburst.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 70% cluster-wind share depends on identifying [S II] n_e profiles with the hot CC85 wind fluid; this phase identification is unverified, and a C2-only re-fit would test it.","rationale":"Read in good faith, the paper supplies a rich, internally consistent data set and several independent energetics estimates: the soft X-ray luminosity implies ΣĖ* ~ 2×10^43 erg s^-1 and the SFR-based starburst99 estimate gives ~1.4×10^43 erg s^-1, both near the 1.6×10^43 erg s^-1 cluster-wind sum. So the qualitative conclusion that stellar feedback is energetically important does not rest on the CC85 fits alone. What the CC85 fits add is the spatially resolved identification of the resolved clusters as the engines, and the specific '70% of L_shock' number. The phase mismatch between warm [S II] gas and the hot wind is the least secure link in that chain: it is an assumption about what the observable traces, not a disagreement with the community consensus. The reader's weakest_assumption names exactly this issue; I agree. I would not move to reject because the X-ray and SFR cross-checks independently bracket the same energy scale, and the C2-only refit is a feasible, decisive test. The current CONDITIONAL verdict is the right one pending that test.","tokens_in":30730,"tokens_out":10086,"duration_ms":118339,"concrete_test":"Refit the CC85 model (Eqs. 7-13) to radial electron density profiles constructed only from the broad C2 component of the [S II] λλ6716,6731 doublet - the gas independently classified as shock-ionized in Section 3.4 - rather than the single-component fluxes used in Section 4.2. Compare the summed Ė_T with 1.6×10^43 erg s^-1. If the C2-only fits do not reproduce the r^-2 profiles and the fitted Ė_T values, the [S II] densities are tracing the photoionized/ambient gas rather than the cluster-wind fluid, and the 70% claim is unsupported. If the C2-only fits agree within uncertainties, the phase-identification concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central number ΣĖ_T ≳ 1.6×10^43 erg s^-1 (Section 4.2) is obtained by fitting the CC85 adiabatic wind solution (Eqs. 4-13) to radial electron density profiles measured from the [S II] doublet at T_e = 10^4 K. The model density ρ in Eqs. (9)-(13) is the mass density of the hot wind fluid (T ~ 10^7-10^8 K), not the warm [S II]-emitting gas. The paper never supplies a physical relation between these phases - no pressure equilibrium, filling factor, or cooling/entrainment model - and it uses single-component [S II] fluxes, mixing the narrow photoionized C1 gas with the broad shock-tracing C2 gas. If the observed r^-2 profiles are set by the ambient ISM distribution, photoionization stratification, or unresolved shells rather than by the wind, the fitted Ė_T values and the 70% comparison to L_shock lose their physical meaning. The same fits also fix v_infinity = 2000 km/s; for a fixed density profile the inferred Ė_T scales roughly as v_infinity^3, so an unvalidated velocity choice is a separate amplifier of the phase ambiguity. This is the most load-bearing step in the paper because the abstract's 'at least 70%' and the superwind 'confirmation' both rely on ΣĖ_T.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents KCWI and MUSE IFU spectroscopy of the local luminous infrared merger VV 114, together with HST narrow-band imaging, and argues that resolved young star clusters are the dominant feedback engine. The authors decompose optical emission lines into a narrow (C1) and broad (C2) component, associate C2 with radiative shocks via diagnostic line-ratio mixing sequences, find an anticorrelation between shock fraction and cluster surface density, and fit CC85 adiabatic wind models to [S II]-derived radial electron density profiles around 66 clusters. Summing the fitted energy injection rates gives ΣĖ_T ≳ 1.6×10^43 erg s^-1, quoted as 70% of L_shock, and this value is compared with independent estimates from soft X-ray luminosity, SFR-scaled Starburst99 models, and the kinetic energy of the proposed galactic superwind. The paper concludes that cluster winds and supernovae power most of the shock luminosity and drive the superwind in VV 114.","tokens_in":31059,"tokens_out":6809,"duration_ms":74561,"significance":"If the central claim holds, this is an important empirical demonstration that a resolved young cluster population, rather than an AGN or tidal dynamics alone, can supply most of the mechanical energy driving shocks and a galactic superwind in a LIRG. The study makes good use of a rich dataset: high-resolution KCWI and MUSE spectroscopy, stringent δBIC>1000 selection of the broad component, a transparent shock-mixing framework built from MAPPINGS grids, a published HST cluster catalog, and multiple cross-checks of the energy budget against X-ray and SFR-scaled expectations. The paper is therefore significant for feedback studies and for connecting cluster-scale energy injection to galaxy-scale outflows. The headline 70% number, however, rests on a chain of model-dependent conversions, in particular the identification of [S II] density profiles with the hot CC85 wind fluid and a fixed 2000 km/s velocity used in both the cluster-wind fits and the superwind energy budget; the quantitative claim needs revision or additional tests before it can be accepted at face value.","major_comments":[{"comment":"The central quantity ΣĖ_T is obtained by fitting the CC85 adiabatic wind solution to electron density profiles measured from the [S II] doublet at T_e = 10^4 K, using single-component fluxes that mix the C1 and C2 kinematic components. The model density in Eqs. (9)–(13) is the mass density of the hot (10^7–10^8 K) wind fluid, and no physical relation is supplied between that fluid and the warm [S II]-emitting gas, such as pressure equilibrium, a filling factor, or an entrainment/cooling model. If the observed r^-2 profiles are set by photoionization stratification, the ambient ISM distribution, or unresolved shells rather than by the wind, the fitted Ė_T values and the 70% comparison to L_shock lose their physical meaning. A re-fit using only C2 fluxes, or an explicit multiphase model, is needed to support the headline claim.","section":"§4.2, Eqs. (4)–(13)"},{"comment":"The total shock luminosity is estimated as L_shock = 75 × L_Hα,C2, assuming an average shock velocity of 160 km/s from Rich et al. (2010). No uncertainty or sensitivity is quoted for this conversion, even though the C2 line widths and the inferred shock velocities in Figure 9 span roughly 100–500 km/s. Because L_shock is the denominator of the headline 70% share, the authors should quantify how L_shock and the resulting cluster-wind fraction vary when the shock velocity distribution implied by their own C2 measurements is used instead of a single value.","section":"§4.2 (L_shock conversion)"},{"comment":"The superwind kinetic energy budget uses v_exp = 2000 km/s, identical to the terminal velocity v_inf fixed in all CC85 cluster-wind fits. For a fixed observed density profile, the fitted Ė_T scales approximately as v_inf^3 through Eqs. (9)–(13), and the wind kinetic luminosity Ė_kin = E_kin/t_dyn scales as v_exp^3 for fixed bubble radius and mass. The close agreement between Ė_kin ≈ 1.2×10^43 erg/s and ΣĖ_T ≈ 1.6×10^43 erg/s is therefore partly by construction. The independent comparisons with the soft X-ray and SFR-derived mechanical luminosities are not affected, but the statement that the superwind 'is confirmed' to be driven by the starburst should be softened or supported by fits with different fixed velocities or an independently measured wind speed.","section":"§4.4.1 and §4.2"},{"comment":"The anticorrelation between shock fraction and local cluster surface density may be partly produced by the definition of f_shock in Eq. (3): R_mix = (1−f_shock) R_HII + f_shock R_shock. Near clusters, strong HII-region photoionization will lower the inferred f_shock even if the absolute shock emission per unit volume is spatially constant. The authors should verify the trend using the absolute C2 shock luminosity, or by estimating f_shock from C2-only line ratios, before interpreting the anticorrelation as evidence for shocks concentrated in interstitial regions.","section":"§4.1.2 and Fig. 8 (right panels)"}],"minor_comments":[{"comment":"The text reads 'the infrared bolometric luminosity of VV 114 is L_IR = 4.5×10^11 M⊙'; the units should be L⊙, not M⊙.","section":"§4.2 (last paragraph)"},{"comment":"The shock model grid is described as having 'a preshock density of 1 cm−1'; this appears to be a typo for 1 cm^−3.","section":"§3.3"},{"comment":"The treatment of the lower-limit fits flagged with '≥' in Tables A1 and A2 should be stated explicitly: the summation of ΣĖ_T should give the sum both with lower-limit values treated as point values and with them treated as lower bounds, so the reader can see how the 70% figure depends on this choice.","section":"§4.2 and Tables A1/A2"},{"comment":"The histograms and the analysis use δBIC cut based on Hα for selecting spaxels, but the C2 [OIII] and [OII] measurements are drawn from the same spaxels; please justify that the Hα-based BIC selection does not bias the higher-ionization C2 line ratios and kinematics.","section":"§2.5 and Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for ApJ and does not raise novelty or citation concerns. The main risk is the unvalidated phase identification in §4.2; I would ask the authors to attempt the C2-only density-profile fit or to add an explicit multiphase/filling-factor model, and to propagate the shock-velocity and v_inf uncertainties into the 70% claim. If such tests cannot be performed, the quantitative claims should be substantially weakened and reframed as an upper limit or a model-dependent estimate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The paper is a serious, data-heavy attempt to measure where the energy for feedback in VV 114 comes from, and the best part is empirical: a resolved shock-fraction map built from a newly constructed photoionization/shock mixing sequence, and a clean anticorrelation between shock fraction and local surface density of resolved clusters. That anticorrelation is the real observational result and should survive independent of the energetics modeling.\n\nThe new systematic piece is applying the CC85 adiabatic wind model to radial [S II] electron density profiles around 66 clusters and summing the fitted energy injection rates to get ΣĖ_T ≈ 1.6×10^43 erg s^-1, which they compare to the total shock luminosity and claim at least 70% of the shock energy comes from the clusters. The paper is careful in the line decomposition (BIC thresholds, SNR cuts, C1/C2 distinction) and the C2 broad component is convincingly identified with radiative shocks: the line ratios occupy the shock mixing locus, and σ_C2 correlates with predicted shock velocity.\n\nThe soft spot is load-bearing. The [S II] densities are measured at T_e ~ 10^4 K, i.e. in the warm ionized gas, while the CC85 model describes a hot (10^7-10^8 K) volume-filling wind. The paper never supplies the physical link between those phases — no pressure equilibrium argument, no filling factor or entrainment model. If the density profiles are set by photoionization stratification or unresolved shells rather than by the wind itself, the fitted Ė_T values don't mean what the abstract claims. This is the central unresolved question.\n\nThere is also a degree of circularity in the \"confirmation\" of the superwind: the kinetic energy of the bubble is estimated using the same v_exp = 2000 km/s that was fixed in the cluster wind fits, so the agreement between the two energy budgets is at least partly built in. And the fitted Ė_T values are quoted without uncertainties; the tables give no error bars, which matters because the sum is dominated by a few high-Ė_T clusters.\n\nThe L_shock = 75 L_Hα,C2 conversion from Rich et al. (2010) assumes an average shock velocity of 160 km/s; that's another model-dependent factor. So the quantitative \"70%\" is not robust yet; the qualitative conclusion that the starburst cluster population is a major feedback source is plausible and consistent with the data.\n\nWho it's for: people doing resolved feedback in mergers and cluster wind theory. It deserves a serious referee, with major comments asking for a phase-consistency check, error propagation on Ė_T, and an explicit test of the velocity dependence. I'd bring it to a reading group because the phase mismatch is a great discussion topic.","headline":"A data-rich case that cluster winds matter in VV 114, but the headline 70% share rests on an unverified identification of [S II] densities with the hot CC85 wind.","tokens_in":31783,"tokens_out":3400,"would_cite":true,"duration_ms":33500,"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 the starburst merger VV 114, resolved star clusters inject at least $1.6\\times10^{43}$ erg s$^{-1}$, matching about 70% of the galaxy's shock luminosity and powering its superwind.","keywords":["VV 114","luminous infrared galaxies","galactic feedback","star clusters","cluster winds","interstellar shocks","galactic superwind","integral field spectroscopy"],"falsifier":"Measure the same radial electron-density profiles around star clusters that lack young massive stars and around control positions with no cluster: if identical declining profiles appear, the CC85 attribution is wrong. Alternatively, high-spatial-resolution X-ray spectroscopy that resolves the hot gas around a single cluster and finds a density profile inconsistent with $1/r^2$, or a gas temperature far below the $10^7$-$10^8$ K wind, would refute the model fits.","tokens_in":30508,"feed_emoji":"💨","tokens_out":11248,"duration_ms":98518,"temperature":0.7,"pith_summary":"The paper argues that in the merging, infrared-luminous starburst VV 114, the dominant engine of feedback is the population of resolved young star clusters, not primarily the buried AGN or the tidal dynamics of the merger. Using wide-field integral-field spectroscopy, it separates every emission line into a narrow photoionized component and a broad component whose line ratios and kinematics match radiative shocks. Around individual clusters, radial electron-density profiles follow the predictions of adiabatically expanding cluster winds; summing the fitted energy-injection rates gives $\\sum\\dot{E}_T \\gtrsim 1.6\\times10^{43}$ erg s$^{-1}$, about 70% of the total shock luminosity. This budget agrees with independent estimates from soft-band X-rays and from the star-formation rate, and it also matches the kinetic energy needed to drive the observed $\\sim$4 kpc galactic superwind bubble. If correct, this is a resolved, quantitative demonstration that starburst clusters can dominate feedback energetics in a luminous infrared galaxy.","feed_headline":"Star clusters drive 70% of shock energy in a starburst merger","feed_subtitle":"Expanding-wind fits around 66 star clusters reproduce the galaxy's X-ray and outflow energy budgets.","key_machinery":"The load-bearing object is the CC85 model of a steady-state, spherically expanding adiabatic wind: within an injection radius $R_{\\rm sc}$, massive stellar winds and supernovae deposit energy $\\dot{E}_T$ and mass $\\dot{M}_T$ into a hot, overpressurized fluid, and outside $R_{\\rm sc}$ the density falls as $1/r^2$. The paper fits these predicted density profiles to [S II] electron densities measured around each cluster, fixing the terminal wind velocity at 2000 km s$^{-1}$ while leaving $\\dot{E}_T$ and $R_{\\rm sc}$ free; summing the fitted $\\dot{E}_T$ values yields the 70% figure. A second mechanism, the shock-mixing sequence, linearly interpolates between H II-region photoionization models and radiative shock models to assign each spatial pixel a shock fraction $f_{\\rm shock}$ and a shock velocity, which is what connects the broad emission component to cluster-driven shocks.","core_discovery":"The central discovery claim is that the resolved star-cluster population of VV 114 supplies the dominant share of the galaxy's shock energy. The authors fit the CC85 adiabatic spherical wind model to [S II] electron-density profiles around 66 of the 180 cataloged clusters and, after summing the best-fit energy-injection rates, obtain $\\sum\\dot{E}_T \\gtrsim 1.6\\times10^{43}$ erg s$^{-1}$. Comparing this with the total shock luminosity $L_{\\rm shock}=2.3\\times10^{43}$ erg s$^{-1}$ derived from the broad H$\\alpha$ component leads to the conclusion that at least 70% of the shock luminosity is driven by star clusters. The same total agrees with the mechanical energy injection rate estimated from the soft-band X-ray luminosity ($\\sim2\\times10^{43}$ erg s$^{-1}$) and with Starburst99 predictions from the star-formation rate ($1.4\\times10^{43}$ erg s$^{-1}$). A blueshifted broad component across most of the galaxy and a $\\sim$4 kpc bubble in the shock-fraction maps, coincident with the soft X-ray emission, are interpreted as a galactic superwind whose kinetic energy injection rate ($1.2\\times10^{43}$ erg s$^{-1}$ for $n_e=100$ cm$^{-3}$) agrees with the summed cluster winds, confirming the starburst as the wind's power source.","pith_inferences":["If the phase-matching assumption holds, a direct prediction is that the diffuse hot wind around clusters should be detectable in X-ray emission lines with surface-brightness profiles following the same $1/r^2$ falloff, a test the current galaxy-wide X-ray data cannot yet make.","Applying the same CC85 fitting procedure to other LIRGs with HST cluster catalogs would test whether the 70% figure is a general property of starburst mergers or specific to VV 114; the paper's X-ray and star-formation cross-checks provide a template for that survey.","The anticorrelation between shock fraction and cluster surface density suggests wind-wind collisions in the interstitial gas, rather than each cluster independently shocking the ambient medium; high-resolution simulations of clustered winds could discriminate these geometries.","Because the optical IFU misses the most heavily embedded feedback, the 70% may be a lower bound on the starburst's true contribution: near-IR [Fe II] mapping reveals shocked gas hidden behind $A_V>6$ screens whose energy source remains ambiguous."],"forward_implications":["If the 70% figure holds, the resolved star clusters are the primary feedback engine in VV 114, contributing more shock energy than the AGN or the merger dynamics.","The galactic superwind's energy budget is closed by cluster winds: the kinetic energy injection rate of the $\\sim$4 kpc bubble matches $\\sum\\dot{E}_T$, implying the starburst alone drives the outflow.","The estimated shock dissipation timescale for the $10^{10}\\,M_\\odot$ molecular reservoir, roughly 40 Myr against a 140 Myr dynamical time, implies shocks can drain the gas of rotational support and drive it inward, potentially fueling a future starburst.","The [Fe II]/Pa$\\beta$ ratio is shown to be a dust-penetrating shock diagnostic whose spatial pattern tracks the optical line-ratio maps, extending shock mapping into heavily obscured regions.","For clusters with flat density profiles the fitted injection radii and energy rates are lower limits, so the true total cluster energy input may exceed $1.6\\times10^{43}$ erg s$^{-1}$."],"supporting_citations":[{"why":"Supplies the steady-state adiabatic spherical wind model (Eqs. 4-13) used to fit radial density profiles and derive each cluster's $\\dot{E}_T$.","marker":"Chevalier & Clegg (1985)"},{"why":"Provides the HST-based catalog of 180 resolved star clusters whose positions define the sample and the density-profile groupings.","marker":"Linden et al. (2021)"},{"why":"Supplies the soft-band X-ray luminosity ($L_X=2\\times10^{41}$ erg s$^{-1}$) used as an independent mechanical-energy cross-check and the outflow measurements supporting the superwind.","marker":"Grimes et al. (2006)"},{"why":"Provides the conversion from broad H$\\alpha$ luminosity to total shock luminosity ($L_{\\rm shock}\\approx75\\,L_{{\\rm H}\\alpha,{\\rm C2}}$) used to compute the 70% fraction.","marker":"Rich et al. (2010)"},{"why":"Starburst99 models used to convert the star-formation rate into a mechanical energy injection rate for comparison with the cluster-wind sum.","marker":"Leitherer et al. (1999)"},{"why":"Supplies the photoionization and radiative shock models used to build the shock-mixing sequence and assign shock fractions and shock velocities.","marker":"Dopita & Sutherland (1996)"},{"why":"Establishes the dissipative-infall analysis (40 Myr timescale) that connects shock energy to the molecular gas reservoir in a comparable merger.","marker":"Rich et al. (2011)"}],"fun_headline_variants":["Cluster winds power 70% of VV 114's shock energy","Starburst clusters fuel galaxy-wide superwind in VV 114","VV 114's star clusters dominate shock energy budget","Cluster winds match VV 114's X-ray and outflow energy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the [S II] electron-density profiles around clusters trace the hot, adiabatically expanding cluster-wind fluid rather than warmer photoionized gas; if the profiles instead trace photoionization structure, unresolved shells, or the ambient ISM geometry, the fitted mechanical energy rates and the 70% figure lose their physical meaning.","fun_headline_variants_meta":{"raw":{"variants":["Cluster winds power 70% of VV 114's shock energy","Starburst clusters fuel galaxy-wide superwind in VV 114","VV 114's star clusters dominate shock energy budget","Cluster winds match VV 114's X-ray and outflow energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001029,"raw_usage":{"total_tokens":4446,"prompt_tokens":1164,"completion_tokens":3282,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":780,"completion_tokens_details":{"reasoning_tokens":3208}},"tokens_in":780,"tokens_out":3282,"duration_ms":23880,"temperature":1.0,"reasoning_tokens":3208,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:22:54.830649+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same radial electron-density profiles around star clusters that lack young massive stars and around control positions with no cluster: if identical declining profiles appear, the CC85 attribution is wrong. Alternatively, high-spatial-resolution X-ray spectroscopy that resolves the hot gas around a single cluster and finds a density profile inconsistent with $1/r^2$, or a gas temperature far below the $10^7$-$10^8$ K wind, would refute the model fits.","supporting_citations":[{"cited_title":"T., Evans , A","cited_arxiv_id":null,"evidence_quote":"Provides the HST-based catalog of 180 resolved star clusters whose positions define the sample and the density-profile groupings."},{"cited_title":"P., Heckman , T., Hoopes , C., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the soft-band X-ray luminosity ($L_X=2\\times10^{41}$ erg s$^{-1}$) used as an independent mechanical-energy cross-check and the outflow measurements supporting the superwind."}],"review_version":2}