{"id":"389bdf66-032d-42ee-9a7e-4d836809e551","arxiv_id":"2411.16674","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Simulations show that stochastic shear acceleration in the turbulent sheath of a mildly relativistic kiloparsec-scale jet can accelerate protons to near the Hillas limit, suggesting Centaurus A and FR II jets as UHECR sources.","lead":"This paper uses computer simulations of the jets of radio galaxies to test whether they can accelerate cosmic rays to ultra-high energies. The authors find that the turbulent outer layer of a mildly relativistic jet, like the one in Centaurus A, can push protons to EeV energies and may explain part of the observed cosmic-ray anisotropy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Cen A claim is tied to spectra time-integrated to 1000 R0/c, but the physical jet lifetime is only 110-360 R0/c; without a cap the quoted EeV peak and >1 EeV tail are overestimated.","rationale":"The paper is a serious and internally consistent simulation study. The RMHD setup is documented, multiple magnetizations and velocities are explored, and Appendix C provides a reasonable check that the results are not controlled by the choice of injection energy. Credit is due for these robustness efforts. However, the central astrophysical extrapolation has a soft spot that is more specific than the reader's frozen-field concern: even if the frozen snapshot were a perfect representation of the local shear, the time-integrated spectra used for the Cen A conclusion extend to 1000 R0/c, whereas the paper itself estimates the physical jet propagation time at only 110-360 R0/c. In a steady jet with continuous injection and advection, a particle should not contribute to the spectrum for longer than the residence time in the acceleration region; counting the same particles at every snapshot up to 1000 R0/c over-weights the high-energy tail. The reader's weakest assumption is closely related, but the integration-time cap is a sharper and independently checkable issue. If the authors can show that truncating at 360 R0/c leaves the spectral peak near 0.2-0.7 EeV and a non-negligible >1 EeV population, the Cen A dipole claim would stand. Without that demonstration, the claim is conditional at best. Since the reader already assigned CONDITIONAL, and our concern is a concrete version of the same limitation, the verdict should remain unchanged.","tokens_in":12723,"tokens_out":9085,"duration_ms":89262,"concrete_test":"Recompute the FR Ia/b/c time-integrated spectra shown in Figure 6 using t_max=360 R0/c (and, as a bracketing case, t_max=110 R0/c) instead of 1000 R0/c. If the E^2 dN/dE peak drops below ~0.2 EeV or the fraction of particles above 1 EeV falls by more than a factor of ~2, the abstract's claim that Cen A can account for the UHECR dipole is not supported. This requires no new physics, only re-binning the existing particle histories.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 estimates the Cen A jet propagation time as t_jet,min ~ (110-360) R0/c, while Section 3 and Figure 6 integrate test-particle spectra over time up to 1000 R0/c to mimic continuous injection. These two numbers are in tension. In a real kpc-scale jet, a particle is advected downstream and leaves the acceleration region after at most ~t_jet,min; the frozen snapshot plus periodic y-boundary instead keeps simulated particles in the same accelerating shear layer for the full 1000 R0/c. Since the same particles are counted at every snapshot, the time-integrated spectrum over-weights long-lived, high-energy particles. Capping the integration at 360 R0/c will lower the spectral peak and reduce the number of >1 EeV particles, and the paper does not show that the truncated spectra still reach the energies invoked for the Cen A dipole. The FR II conclusions are less affected, since those jets are longer, but the Cen A claim in the abstract depends on the uncapped FR I integration. This is not a criticism of the acceleration mechanism itself, which is supported by the injection-energy test in Appendix C; it is a question of whether the simulated acceleration time is astrophysically allowed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter reports 3D RMHD simulations of mildly relativistic jets (FR I-like beta0=0.6 and FR II-like beta0=0.9) with test-particle tracking of protons in the turbulent spine-sheath, aiming to show that stochastic shear acceleration can produce ultra-high-energy cosmic rays in kiloparsec-scale jets. With parameters guided by Centaurus A (R0 ~ 0.1 kpc, B ~ 27-87 microG, sigma = 0.02-0.2), the time-integrated proton spectra peak at ~0.2-0.7 EeV and extend beyond 1 EeV; for FR II parameters the rigidity spectra peak around ~2 EV. The authors conclude that Cen A's kpc-scale jet could account for the observed UHECR dipole anisotropy and that FR II jets could account for the observed spectrum with a rigidity cutoff at a few EV. The paper includes robustness checks for injection radius (Appendix A), injection energy (Appendix C), and MHD freeze time (FR Ia vs FR Ib), and the central quantitative claim is that the spectral peak is approximately 0.1 times the Hillas energy Emax = q beta B R_j computed from simulated quantities.","tokens_in":117,"tokens_out":8691,"duration_ms":125722,"significance":"If the acceleration result holds, this is a valuable contribution: it provides a concrete, parameter-based demonstration that shear acceleration in the turbulent sheath of a mildly relativistic jet can approach the Hillas limit, with a clean relation between the spectral peak and the Hillas energy. The simulations are internally consistent, use no sub-grid physics, and include checks on injection radius and injection energy; the Hillas limit is an external benchmark rather than a fitted quantity, so the central acceleration claim is not circular. The main weaknesses are in the astrophysical extrapolation: the quoted Cen A spectra are integrated to 1000 R0/c although the jet crossing time is only 110-360 R0/c, and the step from source spectra to the observed UHECR dipole and rigidity cutoff lacks a propagation and deflection model. These issues affect the astrophysical conclusions while leaving the basic acceleration mechanism supported by the simulation evidence.","major_comments":[{"comment":"The Cen A conclusions rely on spectra integrated to t = 1000 R0/c, but the physical jet crossing time quoted in Section 2 is only tjet,min ~ 110-360 R0/c. In the frozen-field setup with periodic y-boundaries, test particles remain in the same shear layer for the full 1000 R0/c and are counted at every snapshot, so the time-integrated spectra over-weight long-lived, high-energy particles. The authors should either cap the integration at 110-360 R0/c or demonstrate that the truncated spectra still peak at ~0.2-0.7 EeV and retain a >1 EeV tail. Without this, the abstract's claim that Cen A's kpc-scale jet could account for the UHECR dipole is not supported by the presented simulations.","section":"Section 2 and Section 3/Figure 6"},{"comment":"The step from the simulated proton spectra to the observed UHECR dipole and rigidity cutoff is made without a propagation and deflection model. The observed dipole amplitude at 4-8 EeV and >=32 EeV depends on the source location relative to the observer, Galactic and extragalactic magnetic fields, energy losses, and composition; a source spectrum peaking at 0.2-0.7 EeV does not by itself imply a dipole at those energies. The authors should either add quantitative transport estimates (even a simplified deflected-propagation calculation) or soften the claims from 'could account for' to 'is consistent with'. This is a gap in the astrophysical interpretation, not in the acceleration simulation itself.","section":"Section 4"},{"comment":"The justification for freezing the RMHD fields at t = tfrozen and continuing test-particle evolution to 1000 R0/c is not fully established. The paper states that this choice 'reflects the acceleration of particles at various locations along the jet,' but all particles in fact see the same frozen turbulent fields for the entire post-freeze time. A more faithful treatment would advect particles downstream or through a time-evolving flow so that their residence time in the acceleration region is limited by the jet length. At minimum, the authors should show that the spectral results are insensitive to replacing the frozen snapshot with a time-evolving sequence up to t ~ 360 R0/c for the FR I runs.","section":"Section 2, frozen snapshot paragraph"}],"minor_comments":[{"comment":"Typo: 'Futher studies' should be 'Further studies'.","section":"Section 4"},{"comment":"The phrase 'The vertical lines denotes' should be 'The vertical lines denote'.","section":"Figure 6 caption"},{"comment":"The paper uses 'EV', 'EeV', and 'Exavolts' for rigidity/energy units; please define these clearly at first use to avoid confusion between rigidity (volts) and particle energy (eV).","section":"Abstract and Section 4"},{"comment":"The normalization factors are said to be chosen differently for visualization, but the relative normalization across runs is not shown; please state explicitly how the spectra are normalized.","section":"Figure 4 caption"},{"comment":"The sentence 'the results are found to be consistent when injecting at higher energies, as discussed in Appendix C' would be clearer if it explicitly noted that the comparison is between FR Ib at late times and FR Ib-hi at earlier times.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a well-executed Letter whose central acceleration mechanism is supported by the simulations and robustness checks. The key issue is the finite jet propagation time for Cen A versus the 1000 R0/c integration used to quote peak energies; this should be addressed with truncated spectra before the Cen A claim is accepted. The propagation/deflection gap is also important for the dipole and rigidity-cutoff statements. The novelty relative to Wang et al. (2023) is mainly the UHECR parameter application and test-particle spectra; the authors should ensure the new elements are clearly demarcated in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid simulation result, but the headline astrophysical claim goes beyond what the simulation actually shows. The key new result is that stochastic-shear acceleration in the resolved, KHI-driven turbulent sheath of a mildly relativistic kpc-scale jet can push protons to roughly 0.1 of the Hillas limit, with time-integrated spectra peaking around 0.1 Emax. This is an output of the simulation, not a fitted parameter, and it is a genuinely new quantitative demonstration.\n\nThe numerics are careful. The grid is high-resolution, the turbulence spectra show a Kolmogorov range over two decades, and the tests on injection radius (Appendix A), injection energy (Appendix C), and frozen time (FR Ia vs Ib) are all reasonable. The authors also check that injecting at a higher energy, where the Larmor radius resonates with the resolved turbulence, gives the same acceleration rate and spectra. For a simulation letter, that is solid work.\n\nThe soft spot is the time integration. The paper estimates Cen A's de-projected jet length to be 6-22 kpc, giving a minimum propagation time tjet,min ~ 110-360 R0/c. But the particle spectra in Figure 6 are integrated up to 1000 R0/c, with periodic y-boundaries that keep particles in the same frozen turbulent sheath for the entire run. In a real jet, a particle is advected downstream and leaves the acceleration region after roughly tjet,min. Truncating the integration at 360 R0/c will lower the spectral peak and reduce the high-energy tail; the paper does not show that the truncated spectra still reach the energies invoked for the Cen A dipole. This is not a critique of the acceleration mechanism itself - the injection-energy test supports that - but it does undermine the abstract's claim about Cen A. The FR II conclusions are less affected, since those jets are assumed to be hundreds of R0 long.\n\nA second, smaller gap is the jump from proton spectra to the observed UHECR dipole and rigidity cutoff. There is no propagation or deflection model, and no flux normalization. The authors themselves hedge here, saying the results 'suggest' and 'it seems likely.' That part should be treated as a suggestion, not a demonstration. Also, no code or data are released; that is a minor issue for a Letter but worth noting.\n\nWho is this for? Anyone working on UHECR source models or shear acceleration in jets. It deserves a serious referee. My recommendation: send it to peer review, but the referee should press for a time-capped analysis for the Cen A case, or a clear argument for why the long integration is physically justified. Without that, the abstract needs to be toned down.","headline":"Solid simulation showing shear acceleration reaches ~0.1 Hillas limit in kpc-scale jets, but the Cen A dipole claim rests on an integration time longer than the physical jet lifetime.","tokens_in":13542,"tokens_out":6137,"would_cite":true,"duration_ms":57147,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.70.Sa","98.54.Gr","95.30.Qd"],"model":"deepseek-v4-flash","headline":"This paper argues that stochastic shear acceleration in the kiloparsec-scale turbulent sheaths of mildly relativistic jets can accelerate protons to EeV energies, potentially explaining the Centaurus A dipole and the UHECR spectrum.","keywords":["ultra-high-energy cosmic rays","shear acceleration","relativistic jets","Centaurus A","radio galaxies","RMHD simulations","test-particle simulations","Kelvin-Helmholtz instability"],"falsifier":"Run the same particle injection without freezing the RMHD fields, stopping the integration at the physical jet crossing time $t_{\\rm jet,min} \\simeq (110{-}360)\\,R_0/c$ for Cen A; if the resulting proton spectra no longer peak in the $0.2$–$0.7$ EeV range or extend beyond $1$ EeV, the claimed explanation of the UHECR dipole by Cen A's kpc-scale jet would be falsified.","tokens_in":12537,"feed_emoji":"🌌","tokens_out":9427,"duration_ms":78195,"temperature":0.7,"pith_summary":"The paper argues that the turbulent sheath of a mildly relativistic kiloparsec-scale jet—the kind seen in Centaurus A—can stochastically accelerate protons to energies around $0.2$–$0.7$ EeV, with tails beyond $1$ EeV, close to the theoretical Hillas limit. If true, this gives a concrete physical mechanism for the observed dipole anisotropy in ultra-high-energy cosmic rays pointing toward Cen A, and, for faster Fanaroff-Riley II jets, a way to produce the observed spectrum with a rigidity cutoff at a few exavolts. The authors reach this conclusion by combining high-resolution relativistic magnetohydrodynamic simulations with test-particle tracking, without sub-grid physics, and show that the time-integrated spectra peak at roughly one tenth of the maximum energy across a range of magnetizations and velocities. The central claim is that stochastic shear acceleration, rather than shock acceleration, is the operative process in these mildly relativistic kpc-scale jets.","feed_headline":"Jet sheaths can push cosmic rays to EeV energies","feed_subtitle":"Cen A-like jets can produce the observed dipole; more powerful jets reach exavolt rigidities.","key_machinery":"The central object is the turbulent spine-sheath layer formed by the Kelvin-Helmholtz instability at the interface between the fast jet spine and the surrounding cocoon. In this layer the flow speed drops smoothly from $\\beta_0$ to zero, and magnetic turbulence is sustained by the instability. Charged particles gyrating in this layer experience repeated head-on and trailing collisions with magnetic fluctuations embedded in the sheared flow, gaining energy on average—a stochastic Fermi-II process that in a sheared velocity profile acts as shear acceleration. The simulations resolve this layer down to the Larmor radius scale of injected particles, then freeze the MHD snapshot at the saturated KHI stage and continue test-particle integration to $1000\\,R_0/c$.","core_discovery":"The paper's central claim is that stochastic Fermi-type acceleration in the Kelvin-Helmholtz-driven turbulent sheath of a jet with bulk speed $\\beta_0 = 0.6$ and magnetization $\\sigma = 0.02$–$0.2$ produces proton spectra peaking at $\\sim 0.2$–$0.7$ EeV and extending beyond $1$ EeV, consistent with the energies required for Cen A to be the source of the UHECR dipole anisotropy. For a more powerful FR II-type jet ($\\beta_0 = 0.9$, radius $R_0 = 1$ kpc), the rigidity spectrum reaches peak energies around $2$ EV on timescales of $120$–$400\\,R_0/c$, suggesting such sources could account for the isotropic UHECR component with a cutoff at a few EV. The paper also finds that the spectral peak satisfies $E_{\\rm peak} \\approx 0.1 E_{\\max}$, where $E_{\\max} = q\\,\\beta\\,B\\,R_j$ is the Hillas limit evaluated with the simulated jet radius.","pith_inferences":["A testable extension would be to run the same test particles in a live, non-frozen MHD turbulence with a physical jet length limit; the paper's FR Ib run ($t_{\\rm frozen}=60R_0/c$) suggests some acceleration on this timescale, but whether the peak stays near $0.5$ EeV for Cen A's actual $110$–$360\\,R_0/c$ crossing time is not yet demonstrated.","The same stochastic-shear mechanism has been proposed for electron acceleration in kpc-scale jets; if protons and electrons share the same turbulent sheath, the simulated cosmic-ray spectra could be combined with multi-wavelength jet models to predict the neutrino and TeV gamma-ray output from pion production in the jet.","One observable implication is that if Cen A's jet magnetic field is at the lower end of the $10$–$60\\,\\mu$G range simulated, the Hillas energy drops, pushing the predicted proton peak below $0.2$ EeV and weakening the dipole connection; future Faraday rotation or spectral ageing measurements of the sheath field can directly narrow this."],"forward_implications":["Cen A's kiloparsec jet becomes a quantitatively viable UHECR source: the simulated proton peak at $0.2$–$0.7$ EeV and the hard spectrum below it are compatible with the observed dipole amplitude rising from $1.7\\%$ at $4$–$8$ EeV to $17\\%$ at $\\geq 32$ EeV.","More powerful FR II galaxies, with faster jets and larger radii, can provide the isotropic component: rigidity peak energies around $2$ EV on timescales of a few hundred $R_0/c$ correspond to the few-EV rigidity cutoff favored by spectrum and composition data.","Mixed composition arises naturally: if heavy elements are entrained into the jet, the same acceleration mechanism pushes them to higher rigidities, so events at $\\sim 100$ EV would be heavy nuclei from FR II jets with high magnetization or large radius.","The empirical scaling $E_{\\rm peak} \\approx 0.1 E_{\\max}$ gives a simple prediction: the peak energy of the cosmic-ray spectrum from a jet is set by one tenth of its Hillas energy, so jet magnetic field and radius measurements directly predict the spectral peak."],"supporting_citations":[{"why":"Supplies the RMHD simulation framework and the KHI sheath profile that the present runs refine at higher resolution.","marker":"Wang et al. 2023"},{"why":"Introduces shear acceleration as a Fermi-type process and gives the analytic basis for the energization mechanism.","marker":"Rieger & Duffy 2004"},{"why":"Provides the modern kinetic treatment of stochastic shear acceleration in relativistic shear flows.","marker":"Lemoine 2019"},{"why":"Defines the maximum-energy bound $E_{\\max}=q\\beta BR_j$ used to normalize the simulated spectra.","marker":"Hillas 1984"},{"why":"Sets the magnetic field and emission constraints for Cen A that fix the simulation parameters.","marker":"H. E. S. S. Collaboration et al. 2020"},{"why":"Measures the Cen A jet speed $\\beta \\sim 0.5$–$0.7$, motivating $\\beta_0=0.6$.","marker":"Hardcastle et al. 2003"},{"why":"Provides the observed dipole anisotropy amplitude and energy dependence that the Cen A hypothesis must match.","marker":"Pierre Auger Collaboration et al. 2017"},{"why":"Gives the recent combined spectrum and composition analysis linking Cen A to the UHECR dipole.","marker":"Abdul Halim et al. 2024"}],"fun_headline_variants":["Cen A's jet sheath accelerates cosmic rays to EeV","Jet turbulence boosts cosmic rays to UHE energies","Kiloparsec jet sheaths reach EeV cosmic ray energies","FR II jets push cosmic rays to exavolt rigidity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that accelerating test particles for up to $1000\\,R_0/c$ in static, frozen turbulence faithfully represents what happens in a real jet, whose turbulence keeps evolving and whose finite length limits Cen A to roughly $110$–$360\\,R_0/c$ of propagation time.","fun_headline_variants_meta":{"raw":{"variants":["Cen A's jet sheath accelerates cosmic rays to EeV","Jet turbulence boosts cosmic rays to UHE energies","Kiloparsec jet sheaths reach EeV cosmic ray energies","FR II jets push cosmic rays to exavolt rigidity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000221,"raw_usage":{"total_tokens":1459,"prompt_tokens":963,"completion_tokens":496,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":428}},"tokens_in":579,"tokens_out":496,"duration_ms":5229,"temperature":1.0,"reasoning_tokens":428,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:49:34.583384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same particle injection without freezing the RMHD fields, stopping the integration at the physical jet crossing time $t_{\\rm jet,min} \\simeq (110{-}360)\\,R_0/c$ for Cen A; if the resulting proton spectra no longer peak in the $0.2$–$0.7$ EeV range or extend beyond $1$ EeV, the claimed explanation of the UHECR dipole by Cen A's kpc-scale jet would be falsified.","supporting_citations":[],"review_version":1}