{"id":"05a4b2de-f594-476f-a36a-ef8d3a8b0830","arxiv_id":"2504.16019","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":20,"one_line_summary":"A leptohadronic jet model of PKS 1717+177 can accommodate the 244 EeV Amaterasu event as a proton, but only with an assumed Lorentz invariance violation and a surprisingly strong extragalactic magnetic field.","lead":"This paper explores whether the blazar PKS 1717+177 could have produced the record-energy Amaterasu cosmic ray, using a jet model that includes both normal matter and protons. The authors conclude this is plausible only if a speculative violation of Lorentz symmetry stretches the proton's travel distance and a strong magnetic field bends its path by 2.5 degrees.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The propagation leg hinges on an input LIV coefficient δhad,0=1e-21 that is neither derived nor independently supported; the paper's own Rmax consistency argument suggests the preferred value may differ, so this unverified BSM parameter is the least secure load-bearing assumption.","rationale":"I read the paper as a motivated scenario, not a detection claim. The SED modeling is transparent, uses public codes (GAMERA, CRPropa), and the authors flag the large Rext, the strong EGMF, the low neutrino rate, and limited statistics. The GMF backtracking result is a genuine computational check and supports the claim of negligible Galactic deflection. However, the bridge from the source to the event is the LIV coefficient; the paper's own discussion makes clear this is a chosen value within an upper bound, not a measured or derived one. The additional coupling to Rmax through the single-event normalization (Fig. 2) adds a quantitative tension that is not resolved. Because the paper explicitly frames the conclusion as a possible origin and provides a detailed multimessenger consistency check, I do not see an internal contradiction that would justify rejection. The conditional verdict is appropriate; the check proposed above would determine whether the LIV assumption is actually sufficient and allowed.","tokens_in":16060,"tokens_out":18656,"duration_ms":191113,"concrete_test":"Run a public UHECR propagation code (e.g., SimProp or CRPropa with a Coleman–Glashow LIV term) for protons injected at 277 EeV (244 EeV at Earth) from z = 0.137 with δhad,0 = 1e-21, including adiabatic, pair-production, and photopion losses, and compare the surviving flux at Earth against the TA exposure used in Sec. 3.2. If the survival probability is below ~10% or the expected event count is not O(1), the association fails. As a cross-check, evaluate the Pierre Auger 2022 public likelihood for δhad,0 = 1e-21; if this value lies outside the 95% CL allowed band, the propagation premise is not viable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The propagation of a 244 EeV proton from PKS 1717+177 at d_c = 590 Mpc is only possible because of the assumed Lorentz invariance violation, δhad,0 = 1e-21 (Sec. 3.2). Without it, the GZK energy-loss mean free path for protons at this energy is O(~10 Mpc), as the paper itself states in Sec. 1, so the association collapses. This coefficient is not derived from the model; it is chosen within the Auger upper limit δhad,0 < 1e-19. A one-sided bound is not positive evidence, and the value sits at the low end of what is needed: the paper's own consistency argument in Sec. 3.2 uses the single-event TA normalization and the SED-inferred Lp = 1.3e45 erg/s to infer log10(Rmax/V) ≈ 19.0–19.2, while the Auger combined-fit cutoff at δhad,0 = 1e-21 is quoted as ~18.6 and increases with δ, implying the internally preferred δ may be somewhat larger than 1e-21. This is not a fatal internal contradiction because the paper labels the scenario as plausible rather than established, but it means the entire propagation leg rests on an unverified BSM parameter. The strong EGMF (B_EG ≈ 1.5 nG) and proton-primary assumption are additional load-bearing choices, but they are secondary: even if both were correct, without the LIV coefficient the 244 EeV proton cannot reach Earth from z = 0.137.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript investigates whether the Telescope Array 'Amaterasu' event at 244 EeV could originate from the blazar PKS 1717+177 at z=0.137, assuming a proton primary and a 2.5-degree angular offset. The authors fit the multiwavelength SED with a one-zone leptohadronic jet model using GAMERA, infer a proton luminosity L_p = 1.3e45 erg/s, and predict a subdominant muon-neutrino flux of about 0.1 events per 10 years at IceCube. They then show that, for a chosen Lorentz invariance violation coefficient δ_had,0 = 1e-21 and an extragalactic magnetic field strength B_EG ≈ 1.5 nG with coherence length 500 kpc, a 244 EeV proton can propagate over the 590 Mpc comoving distance and be deflected by approximately 2.5 degrees, making the association plausible under these assumptions. The paper concludes that PKS 1717+177 is a viable candidate source and motivates multimessenger follow-up.","tokens_in":16514,"tokens_out":4664,"duration_ms":46326,"significance":"If the underlying assumptions are accepted, the manuscript provides a self-consistent multimessenger framework connecting an UHECR event, a blazar SED, and a predicted neutrino flux, with concrete falsifiable consequences for CTA, LHAASO, IceCube-Gen2, and KM3NeT. The use of public codes (GAMERA, CRPropa) and archival multiwavelength data makes the analysis reproducible, and the paper is commendably explicit about the key assumptions (proton primary, LIV, strong EGMF). However, the central claim is conditional: the propagation of the 244 EeV proton from z=0.137 is enabled entirely by an unverified BSM parameter, and the SED fit is underconstrained. These limitations substantially bound the significance of the result, but the paper may still be valuable as a scenario study if the load-bearing assumptions are clearly framed and quantitatively tested.","major_comments":[{"comment":"The entire propagation leg rests on an assumed LIV coefficient δ_had,0 = 10^-21, which is not independently constrained. As the manuscript itself states in Sec. 1, the GZK energy-loss mean free path for protons at this energy is ~10 Mpc, so without LIV a 244 EeV proton cannot travel from z=0.137 (d_c ≈ 590 Mpc) to Earth. The paper notes that the Auger upper limit is δ_had,0 < 10^-19, but a one-sided bound is not positive evidence, and the chosen value sits at the low end of what is needed. Moreover, the manuscript's own consistency argument in Sec. 3.2 infers log10(Rmax/V) ≈ 19.0–19.2 from the TA single-event normalization and the SED-inferred L_p, while the Auger combined-fit cutoff at δ_had,0 = 10^-21 is quoted as ~18.6 and increases with δ, implying that the internally preferred δ may be larger than 10^-21. Please provide a quantitative propagation calculation for this specific source and energy, including the LIV threshold and GZK losses, and discuss the prior or external evidence that justifies the chosen δ; without this, the association collapses.","section":null},{"comment":"The SED fit is heavily underconstrained. Table 1 lists more than a dozen free parameters with no uncertainties and no goodness-of-fit statistic, and the text says only that 'parameter degeneracy was assessed by scanning over a broad range of values' without reporting the scan or its results. The claim that the hadronic component is required to reproduce the VHE gamma-ray and soft X-ray data is not quantitatively supported; a purely leptonic EC model may also fit the data given the sparse VHE points. Please quantify the improvement of the leptohadronic fit over a leptonic-only fit (e.g., a chi-square or likelihood comparison), report confidence intervals on the key derived quantities (L_p, neutrino flux), and show the degeneracy scan explicitly.","section":null},{"comment":"The multimessenger predictions are partly circular. The neutrino flux is computed from the same proton spectrum whose normalization is fixed by the SED fit, and the UHECR event rate is used as an input to fix Rmax in Fig. 2; therefore the neutrino flux and the Rmax–LIV consistency check are not independent tests of the model. Please clarify which quantities are inputs and which are predictions, and ideally compute the expected TA event rate from the SED-normalized proton spectrum and compare it with the single observed Amaterasu event, rather than using that event to set the cutoff.","section":null},{"comment":"The directional association depends on several tuned assumptions: a proton primary, B_EG ≈ 1.5 nG, λ_c = 500 kpc, and the small-angle scattering formula of Eq. (7). The required EGMF is an order of magnitude higher than typical cosmic-void fields (~0.1 nG), and the manuscript does not justify that the line of sight to PKS 1717+177 indeed passes through magnetized filaments with these properties. Please also quantify the chance-coincidence probability of finding a blazar within 2.5 degrees (or within the backtracked angular uncertainty) and test the sensitivity of the inferred B_EG to λ_c and to the scattering model.","section":null}],"minor_comments":[{"comment":"The caption states that 'black and gray error bars correspond to the neutrino flux from TXS 0506+056', but the main text describes black data points and gray upper limits for the electromagnetic SED; please clarify which colors and line styles refer to the neutrino reference flux.","section":null},{"comment":"The equation for the proton injection spectrum is missing a closing brace and the exponential factor is ambiguous; please write it as exp[1 - E_p/(Z R_max)] with proper bracketing.","section":null},{"comment":"The sentence 'We constrain the value of Rmax and hence the required LIV coefficient from the proton luminosity required to explain the blazar SED' is misleading: the LIV coefficient is not derived but chosen, and the argument only shows consistency for a specific δ. Please rephrase to distinguish a consistency check from a measurement.","section":null},{"comment":"The discussion of the cosmogenic gamma-ray spectrum would benefit from a more explicit statement of how the strong EGMF affects the flux and why most interactions must occur near the source, since this is central to the proposed CTA/LHAASO test.","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper. It's a scenario, not a discovery, but it's an honest one and probably the most detailed attempt yet to tie this specific blazar to Amaterasu. What's new: a leptohadronic one-zone SED model for PKS 1717+177, the first for this source, and a propagation scenario with LIV that gets a 244 EeV proton to Earth from z=0.137. The authors are upfront that this requires a proton primary, δhad,0=1e-21, and B_EG≈1.5 nG. They use public codes and clearly state the tensions: large Rext, sparse X-ray data, low neutrino rate.\n\nThe main soft spot is the LIV coefficient. It is chosen, not derived. The paper quotes the Auger upper bound, but a one-sided bound is not positive evidence. There is also an internal consistency issue: from their SED-inferred Lp and alpha=2, they infer log10(Rmax/V)≈19.0–19.2, while the Auger combined-fit cutoff at δhad,0=1e-21 is ≈18.6 and increases with δ. That suggests their preferred scenario may actually need a larger LIV coefficient than the one they adopt. Not fatal, but it means the propagation leg rests on a parameter that is both unverified and possibly not optimally chosen.\n\nThe SED fit itself is plausible but underconstrained. Table 1 has many parameters with no uncertainties or goodness-of-fit, and the hadronic component is not uniquely required; the paper admits a leptonic IC model remains feasible. The neutrino flux is a consistency output, and the UHECR rate is normalized to one event. So the multimessenger 'predictions' are really consequences of the assumptions.\n\nThat said, the paper does what a good scenario paper should: it makes concrete, testable statements. The SED model predicts soft X-ray and sub-TeV gamma-ray signatures that CTA and LHAASO can check, and the next-generation neutrino telescopes could see an event during flaring. The GMF deflection analysis is modern (UF23, CRPropa) and the conclusion that a strong EGMF is needed is clear.\n\nBottom line: worth a serious referee. The core argument is coherent given the assumptions, and the authors are transparent about what is assumed. A referee should push for fit statistics, a sensitivity study of the LIV/EGMF parameters, and a clearer statement that the association is conditional on unverified BSM physics.","headline":"A transparent scenario paper: it plausibly ties PKS 1717+177 to the Amaterasu event, but only under an assumed LIV coefficient and strong EGMF, and the SED fit is underconstrained.","tokens_in":17094,"tokens_out":2628,"would_cite":false,"duration_ms":24786,"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":"The 244 EeV Amaterasu cosmic ray can be a proton from blazar PKS 1717+177 if Lorentz invariance violation lets it cross 590 Mpc.","keywords":["ultra-high-energy cosmic rays","blazars","active galactic nuclei","gamma-ray astronomy","neutrino astronomy","Lorentz invariance violation","extragalactic magnetic fields","leptohadronic jet models"],"falsifier":"Measure the extragalactic magnetic field toward PKS 1717+177; if it is much weaker than 1.5 nG, or if the hadronic LIV bound is pushed below $10^{-21}$, the proposed origin cannot hold.","tokens_in":1964,"feed_emoji":"🌌","tokens_out":6271,"duration_ms":153750,"temperature":0.7,"pith_summary":"The Amaterasu event, a 244 EeV cosmic ray observed by the Telescope Array, arrives from a direction within 2.5 degrees of the blazar PKS 1717+177 at redshift $z=0.137$. This paper argues that this blazar is a plausible source, provided the primary is a proton and hadronic Lorentz invariance violation (LIV) with coefficient $\\delta_{\\mathrm{had},0}=10^{-21}$ lengthens the proton's energy-loss mean free path. A one-zone leptohadronic jet model reproduces the source's radio-to-TeV spectrum, attributes the very-high-energy gamma rays to hadronic cascades, and predicts a neutrino flux an order of magnitude lower than from TXS 0506+056. An extragalactic magnetic field of about 1.5 nG would bend the proton by the observed 2.5 degrees. If correct, the same proton population produces measurable hadronic signatures in gamma rays, neutrinos, and the highest-energy cosmic-ray flux.","feed_headline":"One blazar may explain the 244 EeV Amaterasu cosmic ray","feed_subtitle":"A jet model links the 244 EeV proton, gamma rays, and a neutrino signal—if Lorentz invariance violation opens the way.","key_machinery":"The argument rests on a one-zone leptohadronic jet model that reproduces the SED from radio to TeV energies; it includes synchrotron and external Compton emission from primary electrons plus photohadronic cascades from protons, with the proton normalization fixed by the very-high-energy gamma-ray data. The second pillar is the leading-order LIV dispersion relation $E^2-p^2=m^2+\\delta_{\\mathrm{had},0}E^2$ with $\\delta_{\\mathrm{had},0}=10^{-21}$, which increases the proton energy-loss length so that a 244 EeV proton can traverse 590 Mpc. A small-angle scattering formula then converts the 2.5-degree offset into the required extragalactic field strength of about 1.5 nG.","core_discovery":"The central claim is that PKS 1717+177, a blazar at $z=0.137$, can accelerate and release a proton that is observed as the 244 EeV Amaterasu event. The evidence is a simultaneous fit to the source's multiwavelength spectrum using a one-zone jet model with electron and proton injection; the same proton population that escapes the jet produces the sub-TeV gamma-ray excess through pion-decay cascades inside the jet. Propagation from the source to Earth is enabled by hadronic Lorentz invariance violation with $\\delta_{\\mathrm{had},0}=10^{-21}$, which suppresses photopion losses on the cosmic microwave background, and by a $\\sim1.5$ nG extragalactic magnetic field that deflects the proton by the observed $2.5^\\circ$. The model predicts a subdominant neutrino flux, about an order of magnitude lower than TXS 0506+056, and implies negligible Galactic magnetic deflection.","pith_inferences":["If the association holds, a single cosmic ray becomes a probe of Lorentz invariance violation; a confirmed source would imply that all $\\gtrsim100$ EeV protons can travel across cosmological distances, dramatically widening the search for UHECR sources.","The required 1.5 nG field is characteristic of cosmic filaments rather than voids, suggesting that lines of sight through filaments should show systematically larger UHECR deflections; a population study could test this prediction.","The same modeling template could be applied to future extreme-energy events with candidate blazars, turning single-event associations into a statistical test of LIV parameters and magnetic field models."],"forward_implications":["The very-high-energy gamma-ray emission of PKS 1717+177 can be produced by hadronic cascades inside the jet, not only by leptonic inverse-Compton scattering.","The predicted muon-neutrino flux is about 0.1 events per decade at current IceCube effective area, roughly an order of magnitude below the TXS 0506+056 level.","The same proton population that produces the hadronic signature escapes the jet above about 0.1 EeV, and the single-event detection implies a proton kinetic power of about $2\\times10^{45}$ erg s$^{-1}$ with a rigidity cutoff above $10^{19}$ V.","Galactic magnetic fields deflect the proton by less than the angular uncertainty, so the observed 2.5-degree offset must be imprinted by an extragalactic magnetic field of about 1.5 nG.","If the source's X-ray or gamma-ray activity increases by an order of magnitude, next-generation neutrino telescopes could detect a neutrino within a few years, providing a direct multimessenger confirmation."],"supporting_citations":[{"why":"Supplies the Amaterasu event's energy, arrival direction, exposure, and photon-primary exclusion that define the problem.","marker":"Telescope Array Collaboration 2023"},{"why":"Provides the hadronic LIV bound and the energy-loss-length curves used to justify propagation from $z=0.137$.","marker":"Pierre Auger Collaboration 2022"},{"why":"Introduces the modified dispersion relation that underlies the LIV-based escape from GZK opacity.","marker":"Coleman & Glashow 1997"},{"why":"Gives the short proton energy-loss mean free path at these energies and the small-angle deflection formula linking offset to EGMF strength.","marker":"Dermer et al. 2009"},{"why":"Supplies the parameterizations for pion-decay and Bethe-Heitler secondary spectra used to compute the hadronic cascade and neutrino flux.","marker":"Kelner & Aharonian 2008"},{"why":"Provides the numerical method for modeling hadronic emission in the one-zone jet that the SED fit follows.","marker":"Das et al. 2022a"},{"why":"Extends that hadronic modeling method and is cited alongside Das et al. for the cascade calculation.","marker":"Prince et al. 2024"},{"why":"The CRPropa backtracking simulations that show negligible Galactic deflection for the event directions.","marker":"Alves Batista et al. 2016, 2022"},{"why":"The UF23 Galactic magnetic field model used in the backtracking to establish that the 2.5-degree offset must be extragalactic.","marker":"Unger & Farrar 2024b"},{"why":"The TXS 0506+056 multimessenger association that serves as the comparison point for the predicted neutrino flux.","marker":"IceCube Collaboration et al. 2018"}],"fun_headline_variants":["One blazar may explain 244 EeV Amaterasu cosmic ray","Blazar PKS 1717+177: a source for the Amaterasu event?","Does this blazar solve the 244 EeV cosmic ray puzzle?","Amaterasu's origin: nearby blazar with exotic physics","Blazar jet model links Amaterasu to gamma rays and neutrinos"],"cache_read_input_tokens":18944,"weakest_assumption_plain":"The chain depends on the existence of hadronic Lorentz invariance violation with coefficient near $\\delta_{\\mathrm{had},0}=10^{-21}$; without it, a proton from the blazar cannot survive the 590 Mpc journey.","fun_headline_variants_meta":{"raw":{"variants":["One blazar may explain 244 EeV Amaterasu cosmic ray","Blazar PKS 1717+177: a source for the Amaterasu event?","Does this blazar solve the 244 EeV cosmic ray puzzle?","Amaterasu's origin: nearby blazar with exotic physics","Blazar jet model links Amaterasu to gamma rays and neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00044,"raw_usage":{"total_tokens":2290,"prompt_tokens":1061,"completion_tokens":1229,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":677,"completion_tokens_details":{"reasoning_tokens":1130}},"tokens_in":677,"tokens_out":1229,"duration_ms":8179,"temperature":1.0,"reasoning_tokens":1130,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:13:22.590367+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the extragalactic magnetic field toward PKS 1717+177; if it is much weaker than 1.5 nG, or if the hadronic LIV bound is pushed below $10^{-21}$, the proposed origin cannot hold.","supporting_citations":[],"review_version":1}