{"id":"131413ed-0bdb-4414-beb4-2ec1b6a8ac7a","arxiv_id":"2507.18162","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Numerical solutions of the shear acceleration transport equation with synchrotron losses show mildly relativistic AGN jets can accelerate electrons beyond Lorentz factor 10^8, with an effective cutoff scaling as gamma_eff ~ gamma_e,max / gamma_b0^3.","lead":"This paper models how electrons are accelerated in the fast, shearing outer layers of jets from active galaxies, including energy losses from synchrotron radiation. It finds that mildly relativistic jets can push electrons past 50 TeV, with a new formula for the energy cutoff in spatially averaged spectra.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (13) is stated without derivation; the numerical agreement is the only support for the claimed gamma_eff_max scaling, and a radius-dependent scattering time could alter that scaling.","rationale":"The reader identified the spatially uniform scattering time as the weakest assumption; this is a real concern because Eq. (13) explicitly depends on the radial profile of the acceleration rate via gamma_max(r) ∝ gamma_b(r)^6, and a radius-dependent tau would change that profile. However, the more immediately checkable load-bearing issue is that Eq. (13) itself is stated without derivation, so the numerical agreement in Figs. 1-3 is the only support for the headline scaling. The paper does not provide code or data, and the fits are visual without reported uncertainties. Still, the model is internally consistent under its stated assumptions: the transport equation (4) with tau = tau(p) is well posed, and the qualitative conclusion that larger gamma_b,0 boosts the effective cutoff is plausible. The central quantitative claim, gamma_eff_max ~ 2.4e8, is therefore conditionally acceptable pending an independent derivation or a numerical check of the fitted coefficients. I agree with the reader that the concern does not invalidate the paper, but the omitted derivation of the pivotal Eq. (13) is a more specific and more directly testable vulnerability than the general spatial-variation caveat.","tokens_in":9151,"tokens_out":2524,"duration_ms":21471,"concrete_test":"Independently derive Eq. (13) by performing the radial integral I(p) = 2π ∫ r f(r,p) dr with f(r,p) ∝ exp[-6.5 (p/p_max(r))^{2/3}], p_max(r) = p_* gamma_b(r)^6, for the linear profile (6). Compare the asymptotic large-p slope to the claimed effective cutoff gamma_eff_max = gamma_b,0^3 gamma_*; if the slope is not equivalent to Eq. (13), the central energy estimate shifts. Also re-fit the numerical curves in Figs. 1 and 3 to extract the fitted coefficient a in exp(-a (p/p0)^{2/3}) and check whether a corresponds to gamma_eff_max from Eq. (13) for beta_0 = 0.7 and 0.95.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The manuscript omits the derivation of Eq. (13), the pivotal formula for the headline claim that electrons reach gamma_eff_max ~ 2.4e8. The text states that integrating f(r,p) ∝ exp[-6.5 (gamma/gamma_max(r))^{2/3}] over the cross-section and taking an asymptotic expansion 'reproduces the exponential shape' with gamma_eff_max ≃ gamma_b,0^3 gamma_*. A crucial issue is that this asymptotic approximation assumes the radial integral is dominated by small r, where gamma_max(r) is largest. For the linear profile with beta_0 = 0.7 (gamma_b,0 ~ 1.4), the gamma_b,0^3 factor is only ~2.7, so the effective cutoff is close to gamma_*; yet the numerical result in Fig. 1 shows a cutoff around 600 gamma_0 = 6e7 (with gamma_0 = 1e5), close to Eq. (13) with a modest factor. The larger concern is that the authors infer gamma_max from an exponential fit of form exp(-a (p/p0)^{2/3}) without providing the fitted coefficient values needed to reproduce Eq. (13). If the correct radial weighting in the asymptotic expansion differs from the assumed gamma_max(r) ∝ gamma_b(r)^6 — which is itself only an approximate local balance — the claimed gamma_eff_max scaling with gamma_b,0^3 could be off by an order of magnitude. The paper explicitly defers a non-uniform scattering time to future work (Section 3: 'For convenience... tau = tau0 p^alpha' and 'a more complex situation arises for a non-uniform scattering time; we leave its detailed analysis to future work'), but a spatially varying tau would change the radial profile of gamma_max and hence the gamma_b,0^3 factor in Eq. (13). This is not an internal inconsistency, but the central quantitative claim is not independently supported without the omitted derivation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates synchrotron-limited electron acceleration in relativistic shearing flows, motivated by extended X-ray and VHE emission from large-scale AGN jets. It presents numerical solutions of the steady-state, space-dependent particle transport equation (Eq. 4) for two flow profiles — a linearly decreasing profile (Eq. 6) and a power-law profile (Eq. 7) — assuming a momentum-dependent scattering time tau(p) = tau0 p^alpha with alpha = 1/3, a mono-energetic source at p0, and synchrotron losses. The numerical spectra exhibit sub-exponential cutoffs of the form exp[-(p/p0)^(2/3)], similar to the box-model expectation (Eq. 3). The authors then integrate a locally defined cutoff gamma_max(r) = gamma_* gamma_b(r)^6 over the jet cross-section and claim that the asymptotic expansion of the resulting integral yields an effective cutoff gamma_eff_max = gamma_e,max / gamma_b0^3, quoted in Eq. (13). For parameters representative of mildly relativistic large-scale jets, this gives gamma_eff_max ~ 2.4e8, implying electron Lorentz factors beyond 1e8 and energies above 50 TeV. The paper also compares the effective cutoff with the spatially averaged box-model estimate (Eq. 12) and discusses observational implications, including smoother synchrotron spectra and possible year-scale X-ray variability.","tokens_in":9541,"tokens_out":4796,"duration_ms":46650,"significance":"If the central result holds, the paper offers a falsifiable prediction for the electron cutoff in large-scale AGN jets and a quantitative bridge between the commonly used box model and the space-dependent shear acceleration formalism. The numerical solutions to Eq. (4) provide a useful benchmark for future treatments that include spatial dependence of the scattering time or magnetic field. The strength of the paper lies in its clear statement of the transport equation, boundary conditions, and parameter choices, and in the explicit comparison of numerical cutoffs with analytic expectations for two flow profiles. However, the pivotal effective-cutoff formula (Eq. 13) is stated without a derivation, and the numerical method is described only as a 'finite element method' with no details on mesh, tolerances, or convergence. These omissions make the central quantitative claim difficult to verify and currently limit the paper to a conditional contribution.","major_comments":[{"comment":"The derivation of the effective cutoff Lorentz factor gamma_eff_max is not shown. The text states that integrating f(r,p) proportional to exp[-6.5 (gamma/gamma_max(r))^(2/3)] over the cross-section and taking the asymptotic expansion of the resultant imaginary error function 'reproduces the exponential shape' with gamma_eff_max ~ gamma_b0^3 gamma_*, but no intermediate steps are given. This is the load-bearing formula of the paper: the conclusion that electrons reach gamma_eff_max ~ 2.4e8 depends entirely on this expansion. Please provide the full derivation, including the integral representation, the change of variables, the asymptotic expansion, and the condition under which the small-r region dominates. Without this, Eq. (13) is unsupported and cannot be checked.","section":"Section 3.3, Eq. (13)"},{"comment":"The numerical fits used to infer gamma_max are not quantified. The figure captions quote fit coefficients 0.09 and 0.05 in the exponentials exp(-0.09 [p/p0]^(2/3)) and exp(-0.05 [p/p0]^(2/3)), respectively, but the text does not explain how these coefficients map to the claimed reference values gamma_max = 600 gamma0 = 6e7 (beta0=0.7) and gamma_max = 1.5e3 gamma0 = 1.5e9 (beta0=0.95). Reconstructing the conversion is essential to judge the stated 'very good agreement' with Eq. (13), especially because the fit coefficients are not given for the power-law profile case. Please provide the fitted coefficient values and the explicit relation used to convert them into gamma_max.","section":"Section 3.1, Figures 1 and 3"},{"comment":"The assumption that the scattering time tau depends only on momentum, taup = tau0 p^alpha, is central to the derivation of Eq. (13), because the radial integral assumes gamma_max(r) = gamma_* gamma_b(r)^6 with a spatially constant prefactor. The authors acknowledge in Section 3 that a non-uniform scattering time is left to future work, but they do not assess how sensitive the effective cutoff is to this simplification. Since a spatially varying tau would change the local acceleration and escape rates, it could alter the radial weighting in the integral and hence the gamma_b0^3 scaling in Eq. (13). Please provide a quantitative robustness test, for example a simple power-law tau(r) dependence, to show that the claimed effective cutoff is not an artifact of the constant-tau assumption.","section":"Section 3, Eq. (4), and Section 3.3"},{"comment":"The numerical method is described only as 'using a finite element method' with no details on the mesh resolution, element order, tolerance parameters, or convergence checks. The boundary condition at r=0 is stated, but the treatment of the outer boundary at r=r2 is not specified. Since the paper's central argument relies on comparing numerically extracted cutoffs with Eq. (13), the absence of these details prevents reproduction of the results. Please include a description of the discretization and a convergence test, or refer to a publicly available code/script if one was used.","section":"Section 3, numerical method"}],"minor_comments":[{"comment":"There is a typo in the phrase 'linearly deceasing flow profile' in Section 3.1; it should read 'linearly decreasing'.","section":"Section 3.1 and Figure 1"},{"comment":"The sentence 'for which evidence has been been recently reported' contains a duplicated 'been'.","section":"Section 4"},{"comment":"The parameter set is stated as B = 10 microG at the end of Section 3, but Figure 4 uses B = 20 microG; please state this explicitly in the main text before the figure caption so the reader is not confused.","section":"Section 3.2, Figure 4"},{"comment":"The definition of w in Eq. (11) is unclear: is beta02 the relative velocity between the on-axis flow and the outer radius? The text says beta02 = (beta0 - beta2)/(1 - beta0 beta2), but beta2 = beta(r2) is not defined until later. Please clarify the notation at first use.","section":"Section 3.3, Eq. (11)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a compact paper building on the authors' prior work. The missing derivation of Eq. (13) is the main obstacle; it is a short calculation that should be straightforward to add. The numerical method section also needs more detail. If the derivation turns out to have a subtle flaw (for instance, if the asymptotic expansion is not dominated by the small-r region for the low-beta0 case), the gamma_b0^3 scaling could change, which would weaken the claimed enhancement beyond the box-model estimate. The paper is within scope for ApJ and the topic is relevant, but it currently lacks the technical transparency needed for a solid pubblicazione. I recommend major revision rather than rejection because the issues appear fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a modest but real advance in the shear acceleration program. The new content is the space-dependent numerical solutions with synchrotron losses, plus the effective cutoff formula (Eq. 13) showing that the spatially integrated spectrum cuts off at gamma_e,max / gamma_b0^3 rather than at the local maximum. That is a useful correction to the box model.\n\nThe paper does a decent job: the transport equation is stated, the two flow profiles are sensible, and the numerical solutions are compared with the expected exponential shape for both beta0=0.7 and 0.95. The claimed agreement with Eq. (13) is plausible, and the authors are upfront that a spatially varying scattering time is deferred to future work. The heavy citation of prior work by the same group is normal here; this is a self-consistent extension of that earlier box-model work, not an independent test.\n\nThe soft spots: the derivation of Eq. (13) is a one-liner. The coefficient 6.5 clearly comes from Eq. (3) with alpha=1/3, but the asymptotic expansion of the integral is not shown, so the reader cannot verify the gamma_b0^3 factor without redoing the math. The numerical validation is visual only—no error bars, no convergence study, no released code. And the spatially constant tau assumption, while acknowledged, is potentially the most fragile part because real turbulence will vary across the shear layer. That could change the radial weighting and hence the gamma_b0^3 factor. These are real concerns but they do not break the paper; they make it a conditional accept rather than a clean one.\n\nI would send this to peer review. It is not groundbreaking, but it is a legitimate contribution that gives phenomenologists a handy formula. The referee should ask for a fuller derivation of Eq. (13) and, ideally, a convergence check or the code.","headline":"Solid incremental step in shear acceleration: new numerical solutions with synchrotron losses and a practical effective-cutoff formula, but the key derivation is underreported.","tokens_in":10149,"tokens_out":3659,"would_cite":true,"duration_ms":33164,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Synchrotron cooling still allows shear acceleration to push electrons in mildly relativistic AGN jets to Lorentz factors above 10^8, beyond 50 TeV.","keywords":["Fermi shear acceleration","synchrotron losses","relativistic jets","active galactic nuclei","electron transport","high-energy astrophysics","X-ray emission"],"falsifier":"Measure a spatially resolved X-ray synchrotron spectrum of a large-scale AGN jet and, using independent estimates of the magnetic field, shear width, and spine Lorentz factor, infer the electron cutoff Lorentz factor from the location and shape of the spectral rollover. If the inferred cutoff is more than an order of magnitude below $\\gamma^{\\rm eff}_{\\rm max}$ from Eq. (13) for the same parameters, or if the high-energy rollover is a sharp exponential rather than the predicted sub-exponential with flux index $\\eta/(2+\\eta)=1/4$, the synchrotron-limited shear acceleration picture is falsified.","tokens_in":8899,"feed_emoji":"⚡","tokens_out":12747,"duration_ms":118255,"temperature":0.7,"pith_summary":"Fermi-type shear acceleration—particles gaining energy by scattering off velocity differences in the flow—is a leading candidate for keeping ultra-relativistic electrons in the kilo-parsec-scale jets of active galactic nuclei, but synchrotron losses could in principle cap their energy. This paper argues that the cap is much higher than a homogeneous box model suggests: because the fastest shear acceleration happens where the flow Lorentz factor $\\gamma_b$ (the relativistic boost) is largest, near the inner shear boundary, the local maximum electron energy scales as $\\gamma_b^6$. After averaging over the jet cross-section, the integrated electron spectrum shows a sub-exponential cutoff with an effective maximum Lorentz factor $\\gamma^{\\rm eff}_{\\rm max}\\simeq 2.4\\times10^8\\,\\xi_2^{-3/2}(\\gamma_{b,0}/2)^{3}(100\\,{\\rm pc}/\\Delta r)^{2}(10\\,\\mu{\\rm G}/B)^{7/2}$, where $\\xi_2$ is the turbulence normalization in units of 0.2. If this is correct, mildly relativistic large-scale jets can accelerate electrons beyond 50 TeV, and the radio and X-ray emitting regions need not be cospatial.","feed_headline":"Shear flows can push jet electrons past 50 TeV","feed_subtitle":"Numerical transport solutions put the effective electron cutoff near Lorentz factor 10^9 in mildly relativistic jets.","key_machinery":"The load-bearing object is the generalized, steady-state particle transport equation for a cylindrical relativistic shear flow, in which the acceleration is controlled by the shear coefficient $\\Gamma_s(r)\\propto \\gamma_b(r)^4(\\mathrm{d}\\beta/\\mathrm{d}r)^2$ together with a momentum-diffusion operator based on the scattering time $\\tau=\\tau_0 p^\\alpha$ and a spatial diffusion term $\\kappa=c^2\\tau/3$. The argument's key move is to balance the local shear acceleration time $t_{\\rm acc}(r)\\propto[\\gamma_b^4(\\mathrm{d}\\beta/\\mathrm{d}r)^2]^{-1}$ against the synchrotron cooling time, which yields a local maximum Lorentz factor proportional to $\\gamma_b^6$; averaging the resulting cutoff over the jet cross-section then compresses this to the effective scaling $\\gamma^{\\rm eff}_{\\max}\\propto \\gamma_b^3$. For Kolmogorov turbulence ($\\alpha=1/3$) the integrated electron spectrum develops the sub-exponential cutoff $\\exp[-(\\gamma/\\gamma_{\\max})^{2/3}]$, whose shape is the observable signature of the mechanism.","core_discovery":"The paper's central claim is that synchrotron losses do not limit electron acceleration in mildly relativistic shearing flows to the low values implied by a homogeneous box model. Solving the steady-state, spatially dependent transport equation for the phase-space distribution $f(r,p)$ in a cylindrical jet with shear coefficient $\\Gamma_s(r)=c^2\\,\\gamma_b(r)^4(\\mathrm{d}\\beta/\\mathrm{d}r)^2/15$ and scattering time $\\tau=\\tau_0 p^\\alpha$ (with $\\alpha=1/3$ for Kolmogorov turbulence), the authors find that the highest momenta are produced where $\\gamma_b(r)$ is largest, near the inner shear boundary. The local balance of acceleration and synchrotron cooling gives $\\gamma_{e,\\max}\\propto \\gamma_b^6$, and integrating over the jet cross-section turns this into the effective cutoff $\\gamma^{\\rm eff}_{\\rm max}\\simeq 2.4\\times10^8\\,\\xi_2^{-3/2}(\\gamma_{b,0}/2)^{3}(100\\,{\\rm pc}/\\Delta r)^{2}(10\\,\\mu{\\rm G}/B)^{7/2}$. Numerical solutions for linear and power-law velocity profiles confirm this estimate and reproduce a sub-exponential cutoff of the form $\\exp[-(\\gamma/\\gamma_{\\max})^{2/3}]$ in the spatially integrated spectrum. The result implies electron energies beyond 50 TeV in mildly relativistic large-scale AGN jets and suggests a multizone picture in which X-ray-emitting particles are concentrated near the jet axis or inner sheath.","pith_inferences":["A direct observational test of Eq. (13) is to measure, for a sample of large-scale jets with independent estimates of $B$, $\\Delta r$, and $\\gamma_{b,0}$, whether the inferred electron cutoff energy follows the predicted $B^{-7/2}(\\Delta r)^{-2}(\\gamma_{b,0}/2)^3$ scaling.","Because the model assumes a scattering time that depends only on momentum, the most immediate extension is to let $\\tau$ vary with radius and test whether the effective cutoff remains as high as Eq. (13) predicts when turbulence strength changes across the shear layer.","The same transport framework could be applied to other relativistic shear flows, such as gamma-ray-burst jets or pulsar wind nebulae, where the generic signatures would be a local cutoff scaling as $\\gamma_b^6$ and a sub-exponential integrated cutoff whose index depends on the turbulence spectrum.","If the multizone picture is right, spatially resolved X-ray observations should find the X-ray-emitting particles concentrated toward the jet axis or inner sheath, with shorter variability timescales than the radio-emitting outer sheath."],"forward_implications":["Electrons in mildly relativistic large-scale AGN jets can be accelerated to Lorentz factors of $10^8$–$10^9$, corresponding to energies beyond 50 TeV, even when synchrotron losses are included.","The spatially integrated electron spectrum has a sub-exponential cutoff with index $\\eta=2/3$, which translates into a smooth synchrotron spectrum with flux cutoff index $\\eta/(2+\\eta)=1/4$; observed X-ray rollovers should therefore be gradual rather than sharp.","The highest-energy electrons are concentrated near the inner shear boundary, so radio and X-ray emitting regions need not be cospatial, favoring a multizone interpretation of large-scale jet emission.","The spatially averaged box-model acceleration timescale gives a conservative lower bound on the cutoff; the effective cutoff is higher by roughly a factor $\\gamma_b^3$, reinforcing earlier conclusions about the cosmic-ray potential of large-scale jets.","For face-on viewing, differential Doppler boosting can push the observed cutoff frequency beyond the local $\\gamma_{e,\\max}$ value, which may explain reported year-timescale, few-percent X-ray flux variations in large-scale jets."],"supporting_citations":[{"why":"Supplies the analytical estimate for the local maximum electron Lorentz factor scaling as $\\gamma_b^6$ and the Kolmogorov-scaling turbulence parameters used for the jet (Eq. 9).","marker":"Rieger & Duffy 2019"},{"why":"Provides the leaky-box treatment with synchrotron losses and the exponential cutoff form (Eq. 3) that the numerical solutions are compared against.","marker":"Liu et al. 2017"},{"why":"Establishes the box model and parameter set for large-scale AGN jets, including the magnetic field and coherence length, and the cosmic-ray implications that the effective-timescale result reinforces.","marker":"Wang et al. 2021"},{"why":"Derives the generalized space-dependent transport equation (Eq. 4) with the relativistic shear coefficient, which the numerical solutions solve.","marker":"Webb et al. 2018"},{"why":"Shows how the accelerated spectrum depends on the flow profile and provides the weighting factor for the mean acceleration timescale in the linear profile case.","marker":"Rieger & Duffy 2022"},{"why":"Gives the relation converting an exponential particle cutoff into a smoother synchrotron spectrum, used for the predicted observational cutoff shape.","marker":"Zirakashvili & Aharonian 2007"},{"why":"Motivates the parameter range by reporting very-high-energy emission from a large-scale AGN jet that requires ultra-relativistic electrons.","marker":"H.E.S.S. Collaboration et al. 2020"}],"fun_headline_variants":["Shear flows push jet electrons past 50 TeV","Jet shear boosts electron cutoff beyond 50 TeV","Mildly relativistic shear lifts electron limit to 50 TeV","Electron acceleration in shearing jets tops 50 TeV","Synchrotron limit pushed back by jet shear to 50 TeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that the average time between magnetic scatterings depends only on a particle's momentum, not on where it sits in the jet; if the turbulence strength or its coherence length varies across the shear layer, the local acceleration and escape rates change and the derived cutoff energies would shift.","fun_headline_variants_meta":{"raw":{"variants":["Shear flows push jet electrons past 50 TeV","Jet shear boosts electron cutoff beyond 50 TeV","Mildly relativistic shear lifts electron limit to 50 TeV","Electron acceleration in shearing jets tops 50 TeV","Synchrotron limit pushed back by jet shear to 50 TeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000426,"raw_usage":{"total_tokens":2189,"prompt_tokens":958,"completion_tokens":1231,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":1144}},"tokens_in":574,"tokens_out":1231,"duration_ms":12409,"temperature":1.0,"reasoning_tokens":1144,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:39:10.521515+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a spatially resolved X-ray synchrotron spectrum of a large-scale AGN jet and, using independent estimates of the magnetic field, shear width, and spine Lorentz factor, infer the electron cutoff Lorentz factor from the location and shape of the spectral rollover. If the inferred cutoff is more than an order of magnitude below $\\gamma^{\\rm eff}_{\\rm max}$ from Eq. (13) for the same parameters, or if the high-energy rollover is a sharp exponential rather than the predicted sub-exponential with flux index $\\eta/(2+\\eta)=1/4$, the synchrotron-limited shear acceleration picture is falsified.","supporting_citations":[{"cited_title":"An Introduction to Particle Acceleration in Shearing Flows","cited_arxiv_id":"1909.07237","evidence_quote":"Supplies the analytical estimate for the local maximum electron Lorentz factor scaling as $\\gamma_b^6$ and the Kolmogorov-scaling turbulence parameters used for the jet (Eq. 9)."}],"review_version":1}