{"id":"b68ef645-6ed3-4c8d-b552-1a908752c15f","arxiv_id":"2607.24969","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Jet mechanical luminosity estimated from cm radio continuum implies 1–10% hadronic acceleration efficiency in Gamma-Loud Protostars, consistent with shock-acceleration theory.","lead":"The paper derives a formula that turns radio continuum brightness of young stars into the mechanical power of their jets, then compares that power to gamma-ray output. It finds that protostellar jets can supply the observed cosmic rays at 1–10% efficiency, supporting jets as the accelerators in Gamma-Loud Protostars.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The 1–10% efficiency rests on Eq. 12's fixed τ = 10^5 yr accumulation assumption plus the Ṗ_jet = Ṗ_out substitution; both are order-unity-to-order-of-magnitude levers on ε that are assumed, not measured.","rationale":"The reader identified the Ṗ_jet ≈ Ṗ_out substitution as the weakest assumption; I agree it is load-bearing (ε is sensitive to it as f^{3/2}) but find it shares top billing with the τ = 10^5 yr accumulation assumption in Eq. 12, which is a linear, potentially order-of-magnitude lever on ε that the paper adopts without justification — hence \"partial\" rather than \"agree.\" The internal algebra of Eqs. 1–6 checks out: I verified the substitutions (Eq. 2→3→5→6) reproduce the stated coefficients and exponents (0.71→2.5 with exponent 0.87; 50 L⊙ with exponent 1.10, including the M⊙ yr⁻¹ (km/s)² → L⊙ conversion), so there is no internal inconsistency. The concern is purely about unconstrained astrophysical normalizations, and the paper is unusually candid about its own limitations — it labels the efficiencies \"rough estimates,\" acknowledges the survey bias toward photoionized regions, and calls for dedicated radio follow-up. That candor, plus the fact that the claimed 1–10% is explicitly provisional, means the concern does not warrant moving the verdict off CONDITIONAL; it sharpens what the condition should be: an independent L_jet determination and a defensible confinement time before the efficiency range is quoted as consistent with DSA. The validation in Fig. 1 is partly a one-parameter rescaling (η is fit to make Eq. 10 match the L_out data), so only the exponent agreement (0.65 vs 0.72) constitutes an independent check — worth noting but not itself load-bearing.","tokens_in":10836,"tokens_out":5059,"duration_ms":160910,"concrete_test":"For the GLPs with spatially resolved jet data (e.g., HH 80–81, S255 NIRS 3, where V_jet and geometry are measured independently), compute L_jet directly from Eq. 1/2 using measured V_jet, θ0, and flux — bypassing Eq. 4 and the Ṗ_jet = Ṗ_out step entirely — then recompute ε = W_p/(τ L_jet) for τ spanning 10^4–10^6 yr with t_esc estimated from source size and a standard diffusion coefficient. If ε stays within 1–10% across this range, the claim is robust; if it exits the band (expected if t_esc ≪ 10^5 yr), the 1–10% figure is an artifact of the assumed accumulation time and phenomenological normalization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim — ε ≈ 1–10%, \"consistent with diffusive shock acceleration\" — comes from Fig. 5 via Eq. 12, E_CR ≈ ετL_jet. Two structural assumptions set the absolute value of ε, and neither is constrained by the data presented.\n\n(1) Eq. 12 equates the instantaneous CR reservoir energy inferred from the gamma-ray luminosity (W_p, set by the current L_γ/ρ and the pp-loss or confinement time used in Méndez-Gallego et al. 2026) with the time-integrated jet output over a fixed τ = 10^5 yr. Physically, W_p = ε L_jet × min(τ, t_esc, t_loss). If CRs escape the gamma-ray production region on ~10^4 yr — plausible given the source sizes and typical diffusion coefficients — the same W_p requires ε ten times larger, pushing the range to 10–100% and into tension with the DSA benchmark the paper invokes. ε scales linearly with this unconstrained ratio, and no escape/loss argument is given to justify τ = 10^5 yr as the relevant accumulation time.\n\n(2) L_jet enters through Eq. 6, whose normalization is anchored to the Ṗ_jet ≈ Ṗ_out substitution (Sect. 2, after Eq. 4). Because Ṁ_jet ∝ Ṗ_jet^{1/2} at fixed Sνd² (Eq. 3), L_jet = 0.5 Ṗ_jet²/Ṁ_jet ∝ Ṗ_jet^{3/2}. Any momentum boost or deficit between jet and outflow (energy-driven phases, inclination/opacity corrections to Ṗ_out) shifts L_jet — and hence ε — by f^{3/2}; f ~ 2–3 moves ε by a factor ~3–5. The paper's own asymmetric normalization error on Eq. 6 (50 +120/−40 L⊙) already spans a factor ~2.4.\n\nThe product of these two levers means the quoted 1–10% band is largely a statement about the adopted τ and the phenomenological normalization, not a measurement that discriminates 1% from 30%. The near-parallel L_bol scalings of E_CR (∝ L_bol^0.79) and L_jet (∝ L_bol^0.65) also guarantee a tight \"correlation\" regardless of the underlying physics, so the correlation itself adds little evidential weight beyond the normalizations. The authors do flag the radio-data biases and call the values rough, which tempers but","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript constructs an estimator for the mechanical luminosity of protostellar jets from cm radio continuum luminosity. Starting from the Reynolds (1986) free-free model for a conical ionized jet (Eq. 1), the authors adopt a fixed 10% ionization fraction (Eq. 2), eliminate the jet velocity in favour of the momentum rate (Eq. 3), and then substitute the Anglada et al. (2018) phenomenological Ṗ_out–Sνd² correlation (Eq. 4) under the assumption Ṗ_jet ≈ Ṗ_out, obtaining Ṁ_jet(Sνd²) (Eq. 5) and L_jet(Sνd²) (Eq. 6). The derivation is checked against the Maud et al. (2015) and Li et al. (2018) outflow mechanical-luminosity data (Fig. 1) and found consistent in slope, with a fitted ratio η = L_jet/L_out ≈ 70. The authors then cross-match the Gamma-Loud Protostar (GLP) sample of Méndez-Gallego et al. (2026) with VLASS and SARAO, classify counterparts into protostellar-jet and H II-region candidates (Fig. 2), and compare the hadronic CR energy E_CR inferred from the gamma-ray emission with ετL_jet for τ = 10^5 yr (Eq. 12, Fig. 5), concluding that an acceleration efficiency ε ≈ 1–10% suffices, consistent with diffusive shock acceleration.","tokens_in":11332,"tokens_out":4632,"duration_ms":185059,"significance":"If the result holds, this is the first systematic framework connecting a purely radio-observable quantity to the jet mechanical power of a sample of Galactic hadronic accelerator candidates, and it delivers a falsifiable, testable prediction: dedicated sensitive radio observations of GLPs should recover jets with L_jet such that ε lands in the DSA range. The derivation itself is transparent — every algebraic step from Eqs. (1)–(6) is shown with propagated uncertainties — and the internal consistency check against independent IR/CO outflow data (Fig. 1, slopes 0.65 vs 0.72) is a genuine, non-trivial test of the slope. The use of public all-sky surveys (VLASS, SARAO) makes the method immediately reproducible and applicable to larger samples. The authors are also commendably explicit about the survey-sensitivity bias toward H II regions. The main weakness is that the headline number (ε ~ 1–10%) carries normalization systematics — the fixed accumulation time in Eq. (12), the Ṗ_jet = Ṗ_out substitution, and the documented radio-selection bias — that are currently asserted rather than quantified.","major_comments":[{"comment":"The central quantitative result (ε ≈ 1–10%) follows from equating E_CR — the instantaneous CR reservoir inferred from the current L_γ/ρ under a loss/confinement model in Méndez-Gallego et al. (2026) — with the time-integrated jet output ετL_jet at a fixed τ = 10^5 yr. Physically the reservoir satisfies E_CR ≈ ε L_jet × min(τ, t_esc, t_loss). No argument is given that CR escape and pp-loss times exceed 10^5 yr in the gamma-ray production region; for the ambient densities typical of massive star-forming cores, t_pp can be ~10^4 yr, in which case the same E_CR requires ε an order of magnitude larger (10–100%), in tension with the DSA benchmark the paper invokes. Since ε scales linearly with this unconstrained ratio, the quoted range is not a measurement but a consequence of the adopted τ. The authors should either justify the accumulation time with explicit escape/loss estimates for the GLP","section":"§4.2, Eq. (12) and Fig. 5"},{"comment":"The substitution Ṗ_jet ≈ Ṗ_out is the step that anchors the normalization of Eq. (6), and it is load-bearing: because Ṁ_jet ∝ Ṗ_jet^{1/2} at fixed Sνd² (Eq. 3), L_jet = 0.5 Ṗ_jet²/Ṁ_jet ∝ Ṗ_jet^{3/2}. Any momentum boost or deficit between jet and outflow — energy-driven phases (Ṗ_out > Ṗ_jet), or inclination/opacity corrections to the CO-derived Ṗ_out — shifts L_jet, and hence ε, by f^{3/2}; f ~ 2–3 moves ε by a factor ~3–5. The paper's own asymmetric normalization error on Eq. (6) (50 +120/−40 L⊙) already spans a factor ~2.4, yet Fig. 5 presents ε bands without propagating this normalization systematic into the efficiency range. The momentum-conservation assumption is standard in the field and defensible, but its quantitative leverage on the final ε must be shown explicitly (e.g., a shaded systematic band in Fig. 5), and the authors should state whether the momentum-driven ass","section":"§2, after Eq. (4)"},{"comment":"The efficiency is derived 'ignoring the observational biases' that the paper itself documents: VLASS/SARAO detect only ~30% of RMS YSOs and no close SARAO counterparts, many GLP counterparts are upper limits or H II-region dominated (Fig. 2), and Fig. 4 shows a systematic rightward shift of the jet-classified points relative to the expected relation, which the authors interpret as possible radio-flux excess. Because ε ∝ E_CR/L_jet and L_jet ∝ (Sνd²)^1.1, a radio excess (H II contamination or non-jet emission) inflates L_jet and biases ε systematically low — i.e., directly toward the quoted 1–10% range. The efficiency estimate should either be restricted to the clean protostellar-jet subsample (filled markers in Figs. 2/4) with the H II candidates excluded or shown separately, or the bias should be propagated as a systematic direction on ε (stating whether 1–10% is a lower limit). As writ","section":"§4.2, Fig. 5"}],"minor_comments":[{"comment":"All VLASS (3 GHz) and SARAO (1.3 GHz) fluxes are scaled to 3.6 cm with a single spectral index α = 0.6. For SARAO this extrapolates over a factor ~6.4 in frequency, so Sνd² varies by a factor ~3 for α between 0 and 0.6, and the introduction itself cites negative spectral indices from synchrotron emission in jets. A sentence quantifying the sensitivity of L_jet (∝ (Sνd²)^1.1) to the assumed α, per survey, would be useful.","section":"§4.1"},{"comment":"The statement that Eqs. (7) and (10) 'are consistent to one another' should be sharpened: the test is of the slope only (0.65 ± 0.07 vs 0.72 ± 0.15), since the normalization offset η = 10^1.84 is fitted to the data. It would strengthen the paper to note that η ≈ 70 = V_jet/V_out implies V_out ~ 3 km s^-1 for V_jet = 200 km s^-1, which is a reasonable molecular-outflow velocity and thus a physical consistency check rather than a tautology.","section":"§3, Fig. 1"},{"comment":"The phrase 'assuming V̇_jet = 0' is unclear notation; presumably it means a time-steady jet velocity. Please reword.","section":"§2, Eq. (3)"},{"comment":"The 5 arcmin cross-match radius is large compared with typical protostellar-jet angular sizes and the Fermi localization regions; please justify the choice and comment on the expected false-association rate, since misassociated counterparts feed directly into Fig. 5.","section":"§4.1"},{"comment":"Please state explicitly whether the H II-candidate sources (empty markers in Figs. 2–4) are included in the efficiency comparison, and give the number of sources entering the 1–10% estimate.","section":"Fig. 5"},{"comment":"Minor language/typo items: 'an unique opportunity' (Abstract); 'we are interested in compare' (§2); 'kinetical power' (§4); 'hinting the presence' (Conclusions). Also, the dispersion in the Sνd²–L_bol relation (Eq. 8) is acknowledged in §3 but should be carried through to the uncertainty discussion of Figs. 4–5.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The gamma-ray inputs (E_CR values and the L_γ/ρ–L_bol correlation) are taken entirely from the companion paper Méndez-Gallego et al. (2026, Nature Astronomy, arXiv:2606.19445), which shares the first author and much of the author list. This is not improper, but the present paper's central claim is only as strong as that companion's hadronic interpretation and E_CR estimates, which are not independently assessable here; the editor may wish to confirm the companion paper's status. The derivation paper itself is a legitimate, separable contribution."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful bit here is a transparent L_jet(Sνd²) scaling (Eq. 6) and its L_bol form (Eq. 10), built by folding Reynolds free-free mass-loss into the Anglada radio–momentum and radio–bolometric relations, then checking consistency against Maud/Li outflow luminosities. That check works: they recover η≈L_jet/L_out~70 and mass-loss rates in the right ballpark. Applying it to the GLP sample and getting ε~1–10% is the first systematic efficiency census for these objects. The algebra is clear, circularity is low (the formula does not embed the gamma-ray data), and they are honest that VLASS/SARAO mostly catch H II regions or upper limits, so many points are not clean jet detections.\n\nThe soft spot is real but proportional. Absolute ε comes from E_CR≈ετL_jet with fixed τ=10^5 yr and the Ṗ_jet≈Ṗ_out substitution that normalizes Eq. 6. Escape or loss times shorter than τ move ε linearly; a factor-of-few mismatch in momentum rates moves L_jet (and ε) by f^{3/2}. Their own asymmetric error bar on the normalization already spans ~2–3. The near-parallel L_bol slopes of E_CR and L_jet also mean a correlation is almost guaranteed once normalizations are set, so Fig. 5 is more a consistency check than an independent measurement of efficiency. They flag the radio bias and call the numbers rough; that is fair. Free parameters (ionization fraction, θ0, τ, η) are standard order-unity choices, not hidden knobs.\n\nNo code or source tables, so reproducibility is only moderate. Citations look appropriate; mild self-citation to the GLP catalog is expected and not load-bearing for the radio derivation.\n\nThis is for people working on YSO jets, Galactic CRs, or multiwavelength follow-up of the GLP sample. It is a practical incremental tool, not a reorganization of the field. I would send it to referees: the derivation is solid enough and the caveats are stated. Engage if you need a radio-to-jet-power conversion or are planning deeper cm observations of these sources; treat the 1–10% band as order-of-magnitude until better radio associations and a justified accumulation time exist.","headline":"Clean algebraic tool for L_jet from cm radio, applied to GLPs for a provisional 1–10% efficiency; the absolute ε is set by assumed τ and Ṗ_jet≈Ṗ_out, not measured.","tokens_in":12269,"tokens_out":613,"would_cite":true,"duration_ms":18030,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Protostellar jets can supply the cosmic-ray energy seen in gamma-loud young stars at 1–10% efficiency.","keywords":["protostellar jets","gamma-loud protostars","mechanical luminosity","cosmic-ray acceleration","radio continuum","hadronic emission","acceleration efficiency"],"falsifier":"Deep, high-resolution centimetre maps of the same gamma-loud sources that cleanly separate jet free-free emission from H II-region emission and yield radio luminosities systematically lower (or higher) than the present upper limits would push the inferred efficiencies outside the 1–10 percent band.","tokens_in":11800,"feed_emoji":"📡","tokens_out":822,"duration_ms":14677,"temperature":0.7,"pith_summary":"Gamma-loud protostars are young stars whose jets appear to accelerate protons that produce detectable gamma rays. This paper builds a practical bridge from ordinary centimetre radio brightness to the mechanical power carried by those jets, then compares that power with the cosmic-ray energy required by the gamma-ray data. The derived jet powers match independent infrared and radio trends, and the cosmic-ray energy tracks the injected mechanical energy. The implied conversion efficiency sits between 1 and 10 percent, the range expected for shock acceleration. The result matters because it turns an all-sky radio survey into a census of jet kinetic power and lets observers test whether jets alone can explain the non-thermal emission without extra accelerators.","feed_headline":"Jets power gamma-loud protostars at 1–10% efficiency","feed_subtitle":"Radio brightness yields jet kinetic power; cosmic-ray energy tracks it as shock theory predicts","key_machinery":"The jet mechanical luminosity scaling L_jet / L_⊙ ≈ 50 (S_ν d² / mJy kpc²)^1.10, obtained by substituting the phenomenological momentum-rate–radio relation into the theoretical mass-loss expression under momentum conservation.","core_discovery":"Combining the theoretical free-free mass-loss formula for conical jets with the observed radio–outflow-momentum correlation yields a closed expression for jet mechanical luminosity in terms of centimetre luminosity alone. When this luminosity is integrated over a typical jet lifetime and compared with the cosmic-ray energy needed to produce the observed gamma rays, the two quantities correlate and the required acceleration efficiency falls in the 1–10 percent range.","pith_inferences":["If the efficiency remains near 10 percent across a larger, bias-free sample, protostellar jets could supply a non-negligible fraction of the low-energy cosmic-ray budget inside star-forming regions.","The same radio-to-power conversion can be applied to extragalactic jet-driven systems once analogous free-free and momentum correlations are established.","Systematic under-estimates of jet radio flux would raise the true efficiency ceiling and might require additional acceleration sites after all."],"forward_implications":["Centimetre continuum becomes a direct estimator of jet kinetic power for any young stellar object, not only the gamma-loud subset.","A measured efficiency of a few percent supports diffusive shock acceleration inside protostellar jets as the origin of the gamma rays.","Sources whose radio emission lies far above the jet locus are flagged as likely H II contaminants that may hide fainter true jets.","The same scaling predicts which presently gamma-quiet jets should become detectable once deeper gamma-ray or radio data arrive."],"fun_headline_variants":["Radio cm glow unlocks jet power in gamma-loud protostars","Jet kinetic luminosity from radio alone matches gamma energetics","Protostellar jets accelerate cosmic rays at 1–10% efficiency","cm radio traces mechanical power feeding GLP gamma rays","Free-free radio yields jet L_mech; CR energy tracks at few percent"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The jet and the slower molecular outflow are assumed to carry the same momentum rate, so the radio–outflow correlation can be inserted directly into the jet mass-loss formula.","fun_headline_variants_meta":{"raw":{"variants":["Radio cm glow unlocks jet power in gamma-loud protostars","Jet kinetic luminosity from radio alone matches gamma energetics","Protostellar jets accelerate cosmic rays at 1–10% efficiency","cm radio traces mechanical power feeding GLP gamma rays","Free-free radio yields jet L_mech; CR energy tracks at few percent"]},"model":"grok-4.5","effort":"low","cost_usd":0.001613,"raw_usage":{"total_tokens":832,"prompt_tokens":736,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":16128000,"prompt_tokens_details":{"text_tokens":736,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":20,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":736,"tokens_out":76,"duration_ms":2604,"temperature":1.0,"reasoning_tokens":20,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T04:40:26.133087+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Deep, high-resolution centimetre maps of the same gamma-loud sources that cleanly separate jet free-free emission from H II-region emission and yield radio luminosities systematically lower (or higher) than the present upper limits would push the inferred efficiencies outside the 1–10 percent band.","supporting_citations":[],"review_version":1}