Pith. sign in

REVIEW 3 major objections 6 minor 27 references

A radio view on Gamma-Loud Protostars: Derivation of jet mechanical luminosity

T0 review · 3 major / 6 minor · reviewed 2026-07-31 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Protostellar jets can supply the cosmic-ray energy seen in gamma-loud young stars at 1–10% efficiency.

desk verdict 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. read the letter →

arxiv 2607.24969 v1 pith:EPUZ5EKR submitted 2026-07-27 astro-ph.HE astro-ph.SR

Javier Méndez-Gallego , Carlos Carrasco-González , Rubén López-Coto , Iván Agudo , Stefano Menchiari , Rubén Fedriani , Emma de Oña Wilhelmi , Guillem Anglada This is my paper · ORCID
classification astro-ph.HEastro-ph.SR
keywords protostellarjetsgamma-loudprotostarsmechanicalluminositycosmic-rayaccelerationradiocontinuumhadronicemissionefficiency
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

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.

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 (3)
  1. [§4.2, Eq. (12) and Fig. 5] 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
  2. [§2, after Eq. (4)] 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
  3. [§4.2, Fig. 5] 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
minor comments (6)
  1. [§4.1] 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.
  2. [§3, Fig. 1] 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.
  3. [§2, Eq. (3)] The phrase 'assuming V̇_jet = 0' is unclear notation; presumably it means a time-steady jet velocity. Please reword.
  4. [§4.1] 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.
  5. [Fig. 5] 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.
  6. [General] 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.

Circularity Check

1 steps flagged · score 2.0 of 10

Mild author-overlap self-citation supplies the GLP γ-ray sample; the L_jet derivation and ε ratio are not forced by construction.

  1. self citation load bearing [Sect. 4.2, Eq. (12) and Fig. 5; also Abstract/Introduction citing Méndez-Gallego et al. 2026]
    "Using archival data from VLASS and SARAO, together with the relation derived in Eq. (6), we estimate the ε required to account for the observed γ-ray emission. ... Figure 5 shows the hadronic non-thermal energy estimated to reproduce the observed γ-ray emission under a purely hadronic scenario (Méndez-Gallego et al. 2026)."

    The numerical E_CR (and thus the 1–10% ε band) is taken entirely from the authors’ own prior GLP analysis. That is load-bearing for the efficiency claim, but it is not circular: L_jet is built from radio relations that do not embed L_γ, and ε is their ratio under an external accumulation model. Standard follow-up self-citation, scored mildly.

full rationale

The load-bearing chain for jet mechanical luminosity (Eqs. 1–6) combines an external free-free jet model (Reynolds 1986) with the phenomenological Ṗ_out–radio and radio–L_bol relations compiled in Anglada et al. (2018). Substituting Ṗ_jet ≈ Ṗ_out is an explicit physical assumption, not a quantity defined from the γ-ray data. Validation against Maud et al. / Li et al. L_out–L_bol data (Fig. 1) checks slope consistency and introduces a single normalization factor η ≡ L_jet/L_out; that fit does not enter the efficiency formula. E_CR and L_γ/ρ come from the authors’ prior GLP catalog (Méndez-Gallego et al. 2026) and are independent observables from the radio side. Efficiency then follows from the model relation E_CR ≈ ε τ L_jet with an assumed τ. None of these steps makes ε or the CR–mechanical correlation true by definition of its inputs. The only circularity-adjacent element is ordinary follow-up self-citation of the GLP sample and co-author phenomenology, which is not load-bearing in the self-definitional sense. Absolute ε remains assumption-sensitive (τ, Ṗ_jet/Ṗ_out), but that is a correctness issue, not circularity.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central efficiency claim rests on a short chain: free-free jet theory + empirical radio–momentum correlation + momentum conservation + fixed ionization fraction + fixed jet lifetime + hadronic gamma-ray interpretation. Most pieces are standard domain assumptions; the largest free choices are the 10% ionization fraction, au = 10^5 yr, and the heta_0 average that propagate into the prefactor of Eq. 6.

free parameters (4)
  • ionization fraction Ṁ_ion / Ṁ_jet = 0.1
    Fixed at 10% following Cabrit (2007) and Anglada et al. (2018); directly scales the mass-loss and therefore L_jet.
  • jet lifetime π = 10^5 yr
    Adopted representative value τ = 10^5 yr (Frank et al. 2014) to convert L_jet into total kinetic energy E_CR/ε.
  • jet opening angle θ_0 = 33 ± 16 deg
    Set to 33 ± 16° from Anglada et al. (2018); enters the prefactor and uncertainty of Eqs. 2–6.
  • η = L_jet / L_out scaling = 10^1.84 ± 0.06
    Best-fit factor 10^1.84 ± 0.06 used to overlay Eq. 10 on Maud/Li outflow data in Fig. 1; not used for the efficiency points themselves but validates the framework.
assumptions (5)
  • domain assumption Reynolds (1986) free-free conical-jet model gives Ṁ_ion from S_ u d² (Eq. 1)
    Standard radio-jet theory adopted without re-derivation; Section 2.
  • domain assumption Ṗ_jet ≈ Ṗ_out (momentum conservation between jet and molecular outflow)
    Explicitly invoked after Eq. 4 to substitute the Anglada phenomenological relation into the theoretical mass-loss expression.
  • domain assumption Anglada et al. (2018) empirical Ṗ_out–S_ u d² and S_ u d²–L_bol power laws hold for the GLP population
    Used as Eqs. 4 and 8; the paper checks consistency but does not re-fit.
  • domain assumption Gamma-ray emission of GLPs is hadronic and traces jet-accelerated protons
    Taken from Méndez-Gallego et al. (2026); required to interpret L_\gamma/ ho as CR output.
  • domain assumption Constant spectral index 0.6 for scaling VLASS/SARAO fluxes to 3.6 cm
    Applied uniformly in Section 4.1; standard for partially optically thick thermal jets.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A radio view on Gamma-Loud Protostars: Derivation of jet mechanical luminosity." pith.science (2026). https://pith.science/paper/EPUZ5EKR

@misc{pith2026260724969,
  author       = {Pith},
  title        = {Pith review of: A radio view on Gamma-Loud Protostars: Derivation of jet mechanical luminosity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EPUZ5EKR}},
  note         = {Machine review of arXiv:2607.24969}
}
read the original abstract

Context. Gamma-Loud Protostars (GLPs) have been recently reported as Galactic hadronic accelerators whose acceleration site is situated in their protostellar jets. Theoretical and observational analysis situate radio cm luminosity as a thermal tracer of jet activity, presenting an unique opportunity to study jet properties as accelerators. Aims. We aim to develop a way to estimate the kinetic power of protostellar jets and compare their energetics with those extracted from the non-thermal gamma-ray side of GLPs. Methods. We combine theoretical and phenomenological relations to estimate the jet mechanical luminosity based on the radio cm luminosity. We relate the resulting values to the non-thermal contribution of GLPs, studying its behaviour and efficiency. Results. The derivation of the jet power successfully reproduces infrared and radio observations. The cosmic-ray energy correlates to the injected mechanical energy, implying an acceleration efficiency of 1--10%. Future radio observations are needed to find the definite accelerators and obtain reliable efficiencies.

Figures

Figures reproduced from arXiv: 2607.24969 by the authors.

Figure 1
Figure 1. Mechanical luminosity of molecular outflows versus [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Radio luminosity at 3.6 cm versus bolometric luminosity [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. CR output versus cm luminosity at 3.6 cm. Marker styles [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: Non-thermal CR energy calculated for a hadronic sce [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 1 linked inside Pith

  1. [1]

    F., & Carrasco-González, C

    Anglada, G., Rodríguez, L. F., & Carrasco-González, C. 2018, A&A Rev., 26, 3

    2018

  2. [2]

    T., Padovani, M., & Marcowith, A

    Araudo, A. T., Padovani, M., & Marcowith, A. 2021, MNRAS, 504, 2405

    2021

  3. [3]

    2022, MNRAS, 510, 2277

    Araya, M., Gutiérrez, L., & Kerby, S. 2022, MNRAS, 510, 2277

    2022

  4. [4]

    K., et al

    Beuther, H., Schilke, P., Sridharan, T. K., et al. 2002, A&A, 383, 892

    2002

  5. [5]

    Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883

    1982

  6. [6]

    1996, A&A, 311, 858

    Bontemps, S., Andre, P., Terebey, S., & Cabrit, S. 1996, A&A, 311, 858

    1996

  7. [7]

    2007, in Lecture Notes in Physics, Berlin Springer Verlag, ed

    Cabrit, S. 2007, in Lecture Notes in Physics, Berlin Springer Verlag, ed. J. Fer- reira, C. Dougados, & E. Whelan, V ol. 723, 21

    2007

  8. [8]

    & Bertout, C

    Cabrit, S. & Bertout, C. 1992, A&A, 261, 274

    1992

Show all 27 references
  1. [9]

    & Spitkovsky, A

    Caprioli, D. & Spitkovsky, A. 2014, ApJ, 783, 91 Caratti o Garatti, A., Giannini, T., Nisini, B., & Lorenzetti, D. 2006, A&A, 449, 1077 Caratti o Garatti, A., Stecklum, B., Linz, H., Garcia Lopez, R., & Sanna, A. 2015, A&A, 573, A82 Carrasco-González, C., Rodríguez, L. F., Ang...

    2014

  2. [10]

    2022, ApJ, 930, 28 de Oña Wilhelmi, E., López-Coto, R., Aharonian, F., et al

    Cristofari, P., Blasi, P., & Caprioli, D. 2022, ApJ, 930, 28 de Oña Wilhelmi, E., López-Coto, R., Aharonian, F., et al. 2024, Nature Astron- omy, 8, 425 de Oña Wilhelmi, E., López-Coto, R., & Su, Y . 2023, MNRAS, 523, 105

    2022

  3. [11]

    M., Arce, H

    Dunham, M. M., Arce, H. G., Mardones, D., et al. 2014, ApJ, 783, 29

    2014

  4. [12]

    Fedriani, R., Caratti o Garatti, A., Purser, S. J. D., et al. 2019, Nature Communi- cations, 10, 3630

    2019

  5. [13]

    P., Cabrit, S., et al

    Frank, A., Ray, T. P., Cabrit, S., et al. 2014, in Protostars and Planets VI, ed. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 451–474

    2014

  6. [14]

    D., Camilo, F., et al

    Goedhart, S., Cotton, W. D., Camilo, F., et al. 2024, MNRAS, 531, 649

    2024

  7. [15]

    A., Chandler, C

    Lacy, M., Baum, S. A., Chandler, C. J., et al. 2020, PASP, 132, 035001

    2020

  8. [16]

    2018, ApJ, 867, 167

    Li, Q., Zhou, J., Esimbek, J., et al. 2018, ApJ, 867, 167

    2018

  9. [17]

    L., Hoare, M

    Lumsden, S. L., Hoare, M. G., Urquhart, J. S., et al. 2013, ApJS, 208, 11

    2013

  10. [18]

    T., Moore, T

    Maud, L. T., Moore, T. J. T., Lumsden, S. L., et al. 2015, MNRAS, 453, 645 Méndez-Gallego, J., López-Coto, R., de Oña Wilhelmi, E., et al. 2025, A&A, 695, A11 Méndez-Gallego, J., López-Coto, R., de Oña Wilhelmi, E., et al. 2026, Nature Astronomy [arXiv:2606.19445]

    2015 arXiv

  11. [19]

    O., Hoare, M

    Obonyo, W. O., Hoare, M. G., Lumsden, S. L., et al. 2024, MNRAS, 533, 3862 Article number, page 5 A&A proofs:manuscript no. main

    2024

  12. [20]

    2015, A&A, 582, L13

    Padovani, M., Hennebelle, P., Marcowith, A., & Ferrière, K. 2015, A&A, 582, L13

    2015

  13. [21]

    2016, A&A, 590, A8

    Padovani, M., Marcowith, A., Hennebelle, P., & Ferrière, K. 2016, A&A, 590, A8

    2016

  14. [22]

    S., et al

    Peretti, E., Amato, E., Cerri, S. S., et al. 2026, arXiv e-prints, arXiv:2603.16647

    2026

  15. [23]

    2024, Nature Astronomy, 8, 530

    Peron, G., Casanova, S., Gabici, S., Baghmanyan, V ., & Aharonian, F. 2024, Nature Astronomy, 8, 530

    2024

  16. [24]

    Pudritz, R. E. & Norman, C. A. 1983, ApJ, 274, 677

    1983

  17. [25]

    Reynolds, S. P. 1986, ApJ, 304, 713 Rodríguez-Kamenetzky, A., Carrasco-González, C., Araudo, A., et al. 2016, ApJ, 818, 27 Rodríguez-Kamenetzky, A., Pasetto, A., Carrasco-González, C., et al. 2025, ApJ, 978, L31

    1986

  18. [26]

    A., van der Marel, N., Garufi, A., et al

    Rota, A. A., van der Marel, N., Garufi, A., et al. 2025, A&A, 700, A32

    2025

  19. [27]

    Yang, L.-N., Tang, S., & Tam, P.-H. T. 2026, MNRAS, 549, stag896 Article number, page 6

    2026

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.