{"id":"2b1fc9dc-e60e-44e4-a177-2dc87b83f607","arxiv_id":"2506.15110","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using a double power-law source spectrum with an extended cutoff, the authors predict that Virgo A dominates the highest-energy UHECR flux with a lighter composition than Centaurus A and Fornax A, implying hemispheric differences.","lead":"This paper models how ultra-high-energy cosmic rays from four nearby radio galaxies travel to Earth. It predicts that Virgo A should produce a brighter, lighter-composition signal at the highest energies than Centaurus A or Fornax A, so northern and southern sky observations could differ.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The North-South asymmetry rests on M87's adopted spine Lorentz factor of 7; the paper's own VirA3 run with Gamma=3 removes the high-energy tail, so the claim is not robust to the observationally allowed range.","rationale":"The reader's weakest_assumption already identifies Eq. (2) and the Table 1 values of Eb and <Gamma>_spine as the fragile premise. My reading agrees and sharpens the point: the contrast between VirA1 and CenA1/ForA1 is dominated by the squared Lorentz factor in the exponential cutoff, and the paper's own VirA3 model demonstrates that lowering M87's <Gamma>_spine from 7 to 3 erases the distinctive high-energy tail and light composition. This is not an external or exotic objection; it is an internal sensitivity already present in the manuscript. The concern is load-bearing because the central claim is the predicted North-South asymmetry at E > 10^20.2 eV, and that asymmetry disappears for a plausible and frequently adopted value of the M87 jet Lorentz factor. The paper's robustness claim covers only +/-20% parameter variations, which is far narrower than the factor-of-two uncertainty in the M87 spine Lorentz factor. Despite this, the paper remains a clean model study with a falsifiable prediction, so the appropriate verdict is still CONDITIONAL rather than rejection. Since the reader's verdict is already CONDITIONAL, no verdict change is needed.","tokens_in":20801,"tokens_out":7402,"duration_ms":84342,"concrete_test":"Recompute the Figure 6 comparison with VirA3 (or a new VirA model with <Gamma>_spine = 3-4, all other Table 2 settings unchanged) overlaid on CenA1/ForA1, and evaluate the E > 10^20.2 eV flux and <lnA>. If the Virgo A tail and lighter composition no longer exceed CenA/ForA, the hemispheric asymmetry claim depends on the chosen upper-end M87 Lorentz factor. A secondary check would replace the apparent-speed-derived Gamma with a flow Lorentz factor from a jet-dynamics or MHD model of M87's kpc jet and repeat the calculation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive input is the extended cutoff Z0 Eb <Gamma>_spine^2 in Eq. (2). Because the cutoff enters exponentially, the Virgo A versus Centaurus A/Fornax A contrast is set almost entirely by <Gamma>_spine^2: 49 for VirA1 versus 1.44/2.25 for CenA1/ForA1. Table 1 adopts <Gamma>_spine = 7.0 for M87 from kpc-scale superluminal motions, but those are pattern speeds; flow Lorentz factors inferred for M87's kpc jet have substantial systematic uncertainty, and values near 3 are commonly quoted. The authors' own VirA3 model, identical to VirA1 except <Gamma>_spine = 3, shows a steep dN/dE drop above 10^20.2 eV and a heavier <lnA>, i.e., the behavior the paper attributes to Centaurus A/Fornax A rather than to Virgo A. Thus the headline asymmetry is not a robust consequence of the DPL model; it is a reflection of one end of the allowed <Gamma>_spine range for M87. The stated robustness to +/-20% parameter variation does not cover this much larger, observationally allowed excursion, and the paper acknowledges that further refinement of such source parameters would be speculative.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript adopts a double power-law (DPL) source spectrum with an 'extended' exponential cutoff, previously derived by the authors from relativistic hydrodynamic plus Monte Carlo simulations of FR-type jets, and uses CRPropa to propagate ultra-high-energy cosmic rays from Virgo A, Centaurus A, Fornax A, and Cygnus A to Earth. The authors compute arrival energy spectra and mean logarithmic mass for a grid of jet models (Table 2), varying break energy, spine Lorentz factor, host metallicity, EBL model, and propagation distance range. The central finding is that, because Virgo A is assigned a high spine Lorentz factor, its DPL source has a high-energy tail that survives the short 16 Mpc propagation, producing a brighter and proton-dominated flux above about 10^20.2 eV, while Centaurus A and Fornax A, with lower Lorentz factors, produce heavier and steeper spectra. Cygnus A, despite its high power, contributes little because of its 250 Mpc distance. The authors conclude that, if radio galaxies are the dominant UHECR sources, the energy spectrum and mass composition observed in the Northern and Southern skies should differ at the highest energies.","tokens_in":21086,"tokens_out":7334,"duration_ms":82916,"significance":"If the prediction is correct, it supplies a concrete, testable link between jet kinematics and UHECR spectrum and composition, potentially bearing on the reported PAO/TA differences. The paper's strengths are its systematic model grid, its use of the public CRPropa code with standard EBL and photodisintegration processes, its direct comparison with the conventional single power-law injection, and its unusually explicit statements about parameter uncertainty. The physical mechanism, an extended cutoff scaling as Z0 Eb <Gamma>^2_spine, is well motivated by the authors' earlier simulations. However, the quantitative predictions are substantially less secure than the qualitative trends, because the headline Virgo A versus Centaurus A/Fornax A contrast is controlled by two poorly pinned parameters, Eb and <Gamma>_spine. As a scenario study the paper is valuable; as a robust prediction it currently overreaches.","major_comments":[{"comment":"The load-bearing asymmetry between Virgo A and the southern sources is not robust to the observationally allowed range of M87's spine Lorentz factor. Equation (2) places the exponential cutoff at Z0 Eb <Gamma>^2_spine, so the contrast is controlled by the adopted values: 49 for VirA1 versus 1.44 and 2.25 for CenA1 and ForA1. Table 1 sets <Gamma>_spine = 7 for Virgo A based on kpc-scale superluminal pattern speeds, but pattern speeds are not necessarily flow speeds, and published estimates for M87's kpc jet include values near 3. The authors' own VirA3 run, identical to VirA1 except <Gamma>_spine = 3, shows a steep dN/dE decline above about 10^20.2 eV and a heavier <lnA>, qualitatively matching the behavior the paper attributes to Centaurus A and Fornax A (Figure 4). The +/-20% shaded bands in Figure 6 do not cover this factor-of-two excursion, so the statements in the abstract and Section 4 that the qualitative trends are robust are not supported. Please quantify the threshold in <Gamma>_spine (and Eb) above which the high-energy tail survives, add an observationally motivated lower-Lorentz-factor run for Virgo A, and present the hemispheric asymmetry explicitly as conditional on <Gamma>_spine being near the upper end of the allowed range.","section":"Section 3, Figures 4 and 6; Table 1; Eq. (2)"},{"comment":"The break energy Eb is equally load-bearing, and its acknowledged uncertainty is not propagated into the robustness claim. The adopted relation Eb ~ 45 EeV x phi*xi (Qj/Qn)^alpha relies on phi*xi chosen from a representative range 0.15-0.7, with the text stating that further refinement would be speculative. Because the exponential cutoff scale is Z0 Eb <Gamma>^2_spine, the difference between VirA1 (Eb = 10 EeV) and VirA2 (Eb = 30 EeV) changes both the propagated spectrum and the composition substantially (Figure 4). The +/-20% variations in Figure 6 are smaller than the factor-of-three to five spread admitted for phi*xi in Table 1 and the surrounding text. The authors should either perform a joint sensitivity study over the full adopted ranges of phi*xi and <Gamma>_spine, or explicitly restrict all conclusions to the chosen fiducial values rather than claiming robustness to 'plausible parameter values'.","section":"Section 2.1, Eq. (3); Table 1"}],"minor_comments":[{"comment":"The propagation function is computed for a power-law injection with gamma = -2 and then applied to all DPL and SPL models. Figure 3 shows weak gamma dependence for proton-only and iron-only injections over the distance ranges considered, so the approximation is defensible, but a short validation using the actual mixed-composition VirA1 and CenA1 spectra at the highest energies would remove residual doubt.","section":"Section 2.4, Eq. (6)"},{"comment":"The captions state that the <lnA> values at Earth are artificially flattened above about 10^20.3 eV owing to limited statistics, but the flattening is not marked in the plots; adding a grey band or a note in each panel would prevent readers from interpreting those flat segments as physical.","section":"Figures 4 and 5"},{"comment":"The notation using negative values of s1 and s2 inside exponents with a minus sign is confusing; the authors should define the spectral slopes explicitly (for example, dN/dE0 ~ E^-0.6 below Eb and ~ E^-2.6 above Eb) to avoid ambiguity about the sign convention.","section":"Eq. (2) and Section 2.1"},{"comment":"The phrase 'Northern and Southern Hemispheres' should be 'northern and southern sky exposures,' since Virgo A is visible from both hemispheres and the relevant difference is in the exposure-weighted contributions from specific sources.","section":"Section 3 and Figure 6"},{"comment":"The statement that accounting for magnetic deflections 'up to twice the baseline level' corresponds to the VirA6 model with dp = (1-3)d is imprecise; the maximum distance is doubled but the range shape also changes, so the wording should be adjusted.","section":"Section 4, item 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is transparent, technically careful, and within the scope of the journal. The main issue is the gap between the acknowledged parameter uncertainty and the robust-conclusion language attached to the central hemispheric-asymmetry claim. A threshold analysis over the observationally allowed ranges of <Gamma>_spine and Eb, or a clearly conditional framing of the conclusions, would resolve the concern. I saw no citation or novelty concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read it carefully, and I think the reader's conditional verdict is about right. What's genuinely new here is applying the authors' double power-law source spectrum with the extended exponential cutoff to CRPropa propagation for the four nearest radio galaxies. The headline prediction—that if radio galaxies dominate, the northern sky should be brighter and proton-dominated above ~10^20.2 eV while the southern sky should be heavier and drop more steeply—is concrete, falsifiable, and not present in the prior SPL-based studies. The paper also does some things well: it varies Eb, spine Lorentz factor, metallicity, EBL model, and propagation distance, uses the public CRPropa code, and is honest about the arbitrariness in phi*xi and about the fact that <lnA> at Earth is only reliably computed up to 10^20.3 eV. The propagation-function treatment is a reasonable approximation, and the EBL and distance variations show the qualitative results don't hinge on those choices.\n\nThe soft spot is the one the stress-test note flags, and it lands. The entire North-South contrast between Virgo A and Centaurus A/Fornax A is set by the exponential cutoff at Z0 Eb <Gamma>_spine^2. For VirA1, <Gamma>_spine^2 = 49; for CenA1 and ForA1 it's 1.44 and 2.25. That contrast is the paper. But the adopted 7.0 for M87 comes from kpc-scale superluminal pattern speeds, and flow Lorentz factors near 3 are commonly quoted for that jet. The authors' own VirA3 run, with Gamma=3, shows Virgo A behaving like Centaurus/Fornax: steep drop above 10^20.2 eV and heavier composition. So the headline asymmetry is a prediction conditional on the high end of the allowed Gamma range, not a robust consequence of the DPL model. The paper's robustness claim based on plus/minus 20 percent variations in parameters doesn't cover this much larger, observationally plausible excursion, and the paper itself acknowledges that further refinement of source parameters would be speculative. That should be conceded more directly. The authors don't hide the sensitivity—they discuss it in the VirA1 versus VirA3 comparison—but the summary overstates robustness.\n\nThe self-referential element is real but not a fatal flaw: the source spectrum comes from their own previous simulations rather than being independently derived, but the propagation and composition calculation is separate and standard. Missing code or data artifacts and the lack of statistical/systematic uncertainties on the high-energy composition plots are minor but fair criticisms.\n\nBottom line: this deserves a serious referee. It's a coherent model study with a testable prediction, and the community would benefit from having the conditional nature of that prediction spelled out. I'd sent it out, with the explicit request that the authors reframe the asymmetry as dependent on M87's spine Lorentz factor and add a sensitivity run at a lower, still-plausible Gamma rather than claiming robustness.","headline":"Worth engaging: a clean model study whose new, testable North-South asymmetry prediction is not as robust as claimed—it sits on M87's assumed spine Lorentz factor of 7, and their own run with Gamma=3 erases it.","tokens_in":829,"tokens_out":1074,"would_cite":true,"duration_ms":38003,"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":"Virgo A's faster jet gives the northern sky a brighter, proton-dominated flux above $10^{20.2}$ eV, while the southern sky falls off steeply and stays heavy.","keywords":["ultra-high-energy cosmic rays","Fanaroff-Riley radio galaxies","relativistic jets","double power-law source spectrum","extended exponential cutoff","shear acceleration","cosmic-ray mass composition","cosmic-ray propagation"],"falsifier":"Measure Virgo A's kpc-scale jet spine speed with very long baseline interferometry or proper-motion monitoring: a Lorentz factor below about 3 would push the extended cutoff below $10^{20.2}$ eV and erase the predicted northern proton excess. Alternatively, if a northern-sky observatory above $10^{20.2}$ eV sees a composition as heavy as the southern sky, or no excess toward Virgo A, the central claim fails.","tokens_in":20619,"feed_emoji":"🌌","tokens_out":12982,"duration_ms":116160,"temperature":0.7,"pith_summary":"This paper tries to establish that if nearby radio galaxies are the main sources of ultra-high-energy cosmic rays (particles above about $10^{18}$ eV), the highest-energy sky should look different in the two hemispheres: brighter and more proton-dominated in the north, fainter and heavier in the south. The mechanism is a newly adopted source spectrum, a double power law with an extended exponential cutoff whose characteristic energy grows as the square of the jet spine's Lorentz factor, which gives the fast jet of Virgo A a long high-energy tail that slower jets lack. Propagating this spectrum to Earth with a Monte Carlo code, the authors find that Virgo A's arriving cosmic rays above about $10^{20.2}$ eV are mostly protons, while Centaurus A and Fornax A contribute a steeper, heavier flux. The consequence, if true, is a concrete astrophysical explanation for the apparent north-south differences in cosmic-ray spectra and composition at the highest energies.","feed_headline":"Virgo A's jet would brighten the north's highest-energy cosmic-ray sky","feed_subtitle":"Protons from Virgo A, heavy nuclei from Centaurus A and Fornax A would split the sky at the highest energies.","key_machinery":"The load-bearing object is the double power-law source spectrum with an 'extended' exponential cutoff, Equation (2): $$dN/dE_0\\,dt = S_n f_r(A_0)\\left[(E_0/Z_0E_b)^{-s_1} + (E_0/Z_0E_b)^{-s_2}\\right]^{-1}\\exp[-E_0/(Z_0E_b\\langle\\Gamma\\rangle_{\\rm spine}^2)],$$ with $s_1\\approx -0.5$ to $-0.6$, $s_2\\approx -2.6$, and a cutoff extended by the square of the mean jet-spine Lorentz factor. This form encodes the authors' earlier finding that relativistic shear acceleration pushes particles beyond the break energy $E_b$ up to a cutoff at $Z_0 E_b \\langle\\Gamma\\rangle_{\\rm spine}^2$. The break energy follows $E_b\\approx 45\\,{\\rm EeV}\\,\\phi\\xi\\,(Q_j/Q_n)^\\alpha$ with $\\alpha=1/4$ for Fanaroff-Riley type I (FR-I) jets and $1/3$ for type II (FR-II) jets. A Monte Carlo propagation code then maps the spectrum to Earth, including pair production, photopion losses, and photo-disintegration; the high value $\\langle\\Gamma\\rangle_{\\rm spine}\\approx 7$ for Virgo A is what produces its distinctive surviving proton tail.","core_discovery":"The central discovery is a predicted hemispheric asymmetry in the arriving ultra-high-energy cosmic-ray population, driven by the jet-spine Lorentz factor of the nearest radio galaxies. Using a source spectrum with an extended exponential cutoff at rigidity-scaled energy $Z_0 E_b \\langle\\Gamma\\rangle_{\\rm spine}^2$, the authors find that Virgo A ($\\langle\\Gamma\\rangle_{\\rm spine}\\approx 7$) injects a substantial flux beyond $10^{20.2}$ eV, and that photo-disintegration during the 16 Mpc journey converts much of that heavy-nucleus tail into protons. Centaurus A and Fornax A, with spine Lorentz factors of about $1.2$ and $1.5$, lack this extended tail; their arriving spectra drop steeply above $10^{20.2}$ eV and their composition stays heavier. Cygnus A, despite the highest Lorentz factor, contributes little because its 250 Mpc distance strips the flux. The authors therefore state that if radio galaxies are major sources of ultra-high-energy cosmic rays, the northern sky (dominated by Virgo A) and the southern sky (dominated by Centaurus A and Fornax A) should differ in both flux and composition at the highest energies.","pith_inferences":["If this source spectrum is right, part of the observed discrepancy between northern and southern ultra-high-energy measurements could be astrophysical rather than instrumental; the paper presents the asymmetry as a prediction of the radio-galaxy origin, not as a resolution of that discrepancy.","A full-sky observatory at the highest energies could distinguish this model from a uniform source population: the predicted pattern is a proton excess localized toward Virgo A, with a corresponding deficit elsewhere.","The mechanism should generalize to any nearby jet with a spine Lorentz factor above about 5, so measuring the spine speeds of the full local radio-galaxy sample would turn the model into a quantitative sky map of the highest-energy cosmic rays.","The main confounder is magnetic deflection; the paper brackets path-length increases with one-dimensional ensembles, but a three-dimensional magnetohydrodynamic propagation study would test whether the Virgo A protons remain associated with the source direction at $10^{20.2}$ eV."],"forward_implications":["Above roughly $10^{20.2}$ eV, a northern-sky observatory should see a proton-dominated excess associated with Virgo A, while southern-sky data should show a steeper cutoff and a heavier composition.","The north-south difference should appear mainly at the highest energies; below about $10^{19.5}$ eV the model does not demand a strong hemispheric split in either spectrum or composition.","In this model, the jet's spine Lorentz factor is a first-order control on the arriving composition: lowering it from 7 to 3 steepens the cutoff and markedly increases the mean mass at Earth.","Cygnus A, despite being the most powerful jet considered, should not appear as a high-energy source at Earth, because propagation over 250 Mpc removes nearly all of its flux.","Accurate measurements of Virgo A's kpc-scale jet speed would directly sharpen or weaken the predicted northern high-energy tail."],"supporting_citations":[{"why":"Supply the double power-law source spectrum with extended exponential cutoff and the Monte Carlo transport results that justify it.","marker":"J. Seo et al. (2023, 2024)"},{"why":"Provide the Monte Carlo propagation code used to simulate energy losses and photo-disintegration from source to Earth.","marker":"R. Alves Batista et al. (2016, 2022)"},{"why":"Gives the previous single power-law treatment of the same nearby radio galaxies against which the DPL predictions are compared.","marker":"B. Eichmann et al. (2022)"},{"why":"Supplies the catalog of 42 local radio galaxies used to select the four prominent sources.","marker":"S. van Velzen et al. (2012)"},{"why":"Supports the superluminal-motion estimate of Virgo A's spine Lorentz factor near 7 that drives the extended high-energy tail.","marker":"J. A. Biretta et al. (1999)"},{"why":"Supports the low spine Lorentz factor of Centaurus A that makes its predicted spectrum steep and heavy.","marker":"S. Wykes et al. (2019)"},{"why":"Provides the jet-power scaling used to set the break energy and connect jet dynamics to the source spectrum.","marker":"A. Bhattacharjee et al. (2024)"}],"fun_headline_variants":["Nearest radio jets split cosmic-ray sky by hemisphere","Virgo A's fast jet predicts northern cosmic-ray surplus","Cosmic-ray composition may differ between north and south skies","Radio galaxies could engender asymmetric ultra-high-energy sky","Northern and southern ultra-high-energy cosmic-ray skies diverge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the adopted source spectrum and on Virgo A's jet being as fast as claimed, with a spine Lorentz factor near 7; the paper itself notes that the factor $\\phi\\xi$ setting the break energy is chosen from a representative range and that further refinement would be speculative.","fun_headline_variants_meta":{"raw":{"variants":["Nearest radio jets split cosmic-ray sky by hemisphere","Virgo A's fast jet predicts northern cosmic-ray surplus","Cosmic-ray composition may differ between north and south skies","Radio galaxies could engender asymmetric ultra-high-energy sky","Northern and southern ultra-high-energy cosmic-ray skies diverge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00064,"raw_usage":{"total_tokens":3039,"prompt_tokens":1130,"completion_tokens":1909,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":746,"completion_tokens_details":{"reasoning_tokens":1830}},"tokens_in":746,"tokens_out":1909,"duration_ms":14187,"temperature":1.0,"reasoning_tokens":1830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:43:41.867645+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure Virgo A's kpc-scale jet spine speed with very long baseline interferometry or proper-motion monitoring: a Lorentz factor below about 3 would push the extended cutoff below $10^{20.2}$ eV and erase the predicted northern proton excess. Alternatively, if a northern-sky observatory above $10^{20.2}$ eV sees a composition as heavy as the southern sky, or no excess toward Virgo A, the central claim fails.","supporting_citations":[{"cited_title":"T., Nulsen, P","cited_arxiv_id":null,"evidence_quote":"Supports the low spine Lorentz factor of Centaurus A that makes its predicted spectrum steep and heavy."}],"review_version":1}