{"id":"7630a7a7-a1ef-45bb-a762-c31fc27745d1","arxiv_id":"2501.06677","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Origin likelihood functions for extreme-energy cosmic rays can be computed analytically from photodisintegration cascade rates, giving distance and source rigidity estimates for events like Amaterasu without Monte Carlo sampling.","lead":"Astrophysicists usually cannot trace the most energetic cosmic rays back to their birthplaces, because the particles randomly fragment and bend along the way. This paper proposes an analytic way to compute the likely distance and rigidity of origin for an extreme-energy cosmic ray such as the Amaterasu event, without running Monte Carlo simulations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation 2 as written gives f(L)=0 because each row of Λ sums to zero; the exit vector is undefined, so the central formula is internally inconsistent until e is defined or absorbing states are removed.","rationale":"I read the paper in good faith: it proposes an elegant phase-type formalism for photodisintegration cascades and applies it to the Amaterasu event. The reader identified boost conservation as the weakest assumption. My stress-test found a more elementary issue: the central equation, as printed, is algebraically degenerate. Because λtot_Sn is defined as the sum of all transition rates in row n, and final states are absorbing, every row of Λ sums to zero. With the usual all-ones vector e, Λe = 0 and Eq. (2) returns zero density for all distances. The standard phase-type construction avoids this by using a transient-only subintensity matrix with exit vector −Λ1, or by including absorbing states and selecting them with an indicator vector. The manuscript does neither, so the central formula cannot be evaluated as written. This is a fixable write-up error, likely inherited from prior work, so I would not reject the underlying method; but the paper must define the exit vector precisely and re-derive the figures. I also note the Amaterasu energy inconsistency (1.64e20 vs 244 EeV) and the qualitative CRPropa comparison, but those are secondary to the Eq. (2) issue. The proposed concrete test—evaluating Eq. (2) on a one-step chain—exposes the missing definition immediately. Because the error is correctable without changing the conceptual approach, the reader's CONDITIONAL verdict remains appropriate; my analysis does not move the verdict.","tokens_in":5947,"tokens_out":11558,"duration_ms":127469,"concrete_test":"Apply Eq. (2) to a minimal one-step chain, A→B with rate λ, using the displayed matrix Λ = [[−λ, λ], [0, 0]] and α0 = [1, 0]. With e = (1,1)^T, Λe = (0,0)^T and f(L) = 0 for all L, whereas the correct first-passage density is λ exp(−λL). If the author intended e to be the indicator of final species, the paper must state this and the test should be repeated with e_f = (0,1)^T, which gives the correct density. This single check settles whether Eq. (2) is a notational omission or a substantive error in the central formula.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2, Eq. (2) is the central result, but as printed it is degenerate. In the displayed matrix (3), every row sums to zero by construction: λtot_Sn = Σ_j λ_{S_n→S_j}, and the final species have no outgoing transitions (their diagonal is −λtot_Sf with λtot_Sf = 0). If e is the standard vector of ones, then Λe = 0 for every row and f(L) = α0 exp(ΛL) Λe = 0 identically. The correct phase-type density for first passage to a final set requires either (i) restricting Λ to transient states and using the exit vector t = −Λ1, or (ii) keeping absorbing states in Λ and multiplying by an indicator vector e_f that selects the final species. The manuscript does not define e, so the key equation cannot be evaluated as written. This is not a matter of physical approximation or external consensus; it is an internal algebraic inconsistency. The boost-conservation approximation is a separate, additional concern, but Eq. (2) fails even under the paper's own assumptions.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an analytic method for computing origin likelihood functions of extreme-energy cosmic rays, replacing Monte Carlo propagation with matrix-exponential distributions for photodisintegration cascades. Using the Telescope Array Amaterasu event as a case study, it presents two-dimensional likelihood surfaces over source distance and rigidity for a range of final nuclear species, and argues that for intermediate charge groups (2 < Z ≤ 7) the origin distribution is nearly independent of the injected composition. The method is compared qualitatively with the CRPropa-based analysis of Unger and Farrar (2024b), and the analytic approach is claimed to offer higher resolution and much lower computational cost.","tokens_in":6115,"tokens_out":4920,"duration_ms":51273,"significance":"If the framework is correct, it provides a fast, deterministic alternative to Monte Carlo for a problem where event statistics are extremely sparse, and it yields likelihood surfaces with arbitrary numerical precision. The observation that some final charge groups give composition-independent origin distributions is practically valuable for single-event source localization. The paper also benefits from building on explicit interaction-rate inputs and from a clear application to a concrete event. However, the central formula as printed is internally inconsistent because the exit vector is undefined, and the main physical approximation—exact rigidity boost conservation across the cascade—is asserted without quantitative validation. These issues must be resolved before the claimed advantages can be accepted.","major_comments":[{"comment":"The vector e in f(L) = α0 exp(ΛL) Λe is never defined. With the standard phase-type convention e = (1,...,1)^T, the expression vanishes identically because each row of the displayed matrix Λ sums to zero (including rows for the final species, whose total rate is zero). The correct density requires either restricting Λ to the transient states and using the exit vector t = −Λ1, or keeping absorbing states and multiplying by an indicator vector selecting the final species. As written, Eq. (2) cannot be evaluated and is the central result of the paper; this must be corrected and the notation made explicit.","section":"Section 2, Eq. (2)-(3)"},{"comment":"The paper treats the entire photodisintegration cascade as a Markov chain at a single fixed boost γ = ηR, justified by 'the boost conservation characteristic in photodisintegration interactions.' However, η = Z/A changes between fragments (e.g., 0.464 for 56Fe and 0.5 for 4He), and the interaction rates in Figure 1 vary strongly with γ, especially in the GDR-dominated range. The claim that η changes only slightly is not quantified, and no validation is offered against a propagation model that tracks each fragment's own boost after each disintegration. Since all likelihood surfaces in Figures 3-5 depend on this single-boost assumption, please provide a quantitative error estimate or a comparison with a per-species boost cascade.","section":"Section 2, boost-conservation approximation"},{"comment":"The text states that the observed nuclei energies are distributed as a Gaussian with mean E_low = 1.64 × 10^20 eV, but the Amaterasu event is reported at 244 EeV = 2.44 × 10^20 eV. This discrepancy directly affects the rigidity constraints applied in Figures 3-5: a wrong mean shifts all reported likelihood surfaces. If E_low is not a typographical error, its role must be explained; if it is a typo, it must be corrected to 2.44 × 10^20 eV (or the appropriate value from Unger and Farrar 2024b).","section":"Section 2, Gaussian energy constraint"},{"comment":"For final-state protons and deuterons, the paper explicitly neglects photopion energy losses ('this effect is expected to be reduced and was not included in the results shown in here'), yet the likelihood functions in Figures 3-5 include p and d panels. Because the observed-energy Gaussian is applied at Earth, energy losses during propagation change which source rigidities and distances can contribute to the observed energy window. Please either include the losses within the rewards-transform framework or demonstrate quantitatively that their effect on the displayed likelihoods is negligible.","section":"Section 2, proton energy losses"}],"minor_comments":[{"comment":"The phrase 'makes difficult the task' is ungrammatical; use 'makes the task difficult' or 'complicates the task.'","section":"Abstract"},{"comment":"The notation λ_{S_n→S_f} = Σ_{j=1}^{k} λ_{S_n→S_f_j} uses S_f for both the set of final species and an individual final species; please disambiguate, for example by writing S_f^set and S_f^j.","section":"Section 2, after Eq. (3)"},{"comment":"The phrase 'a maximal rigidity of Rmax ∈ [3, 5] EV (?)' contains a placeholder '?' and no citation; this must be completed.","section":"Section 1"},{"comment":"'corner stone' should be spelled 'cornerstone.'","section":"Section 3"},{"comment":"The label 'starting species (blue)' in the text is ambiguous; specify in each caption that blue denotes the likelihood for the injected nucleus itself.","section":"Figures 3-5"},{"comment":"The statement that the origin is 'independent of the original composition' is stronger than the evidence presented, since only three starting species (28Si, 40Ca, 56Fe) are tested and the independence is shown only for some final charge groups; please soften the wording accordingly.","section":"Section 1 / Section 3"},{"comment":"Since Eq. (2) depends on the rate matrix Λ, and the rates are taken from Morejon (2023) without a self-contained definition, adding a short appendix or explicit references to the equations in that work would improve reproducibility.","section":"Overall"}],"recommendation":"major_revision","confidential_remarks":"This is a short proceedings contribution with a promising analytic framework. The undefined exit vector in Eq. (2) is a serious internal inconsistency in the central result, and the boost-conservation approximation is a load-bearing assumption that needs quantitative support. The corrections appear feasible within the scope of a revision, so I recommend major revision rather than rejection. The paper should also be checked for consistency with the stated Amaterasu energy in the Gaussian constraint before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the core formula, Eq. (2), is broken as written. In the displayed matrix (3), every row sums to zero because the diagonal is the negative sum of the off-diagonal rates, including the final species whose total rate is zero. With e as the usual vector of ones, Λ e = 0 and f(L) = 0 identically. The paper never defines e, so the key equation cannot be evaluated. This is not a physics disagreement; it is an internal algebraic inconsistency. The standard phase-type fix (restrict Λ to transient states, use t = −Λ1 as the exit vector, or select final species by an indicator vector) is well known, but as submitted the central result collapses.\n\nNow the credit. The idea of using matrix-exponential distributions to get analytic origin likelihoods for individual extreme-energy cosmic rays is a genuine step beyond CRPropa brute force. The composition-independence observation for the 2 < Z ≤ 7 group (Figures 3–5) is interesting and is a real, checkable claim. The resolution advantage over Monte Carlo is also real: the likelihood is computed with arbitrary precision, not limited by sampling.\n\nThe soft spots are substantial and mostly in the presentation. The rates and the full matrix-exponential machinery come from Morejon (2023) without derivation. For a proceedings paper that is tolerable, but it makes this text unverifiable in isolation. The Amaterasu energy used in the Gaussian constraint (E_low = 1.64e20 eV) does not match the 244 EeV quoted everywhere else; that shifts the rigidity axes and could change the figures. No code or tables are released, and the comparison to Unger & Farrar is qualitative. Proton photopion losses are neglected, but the text admits this.\n\nThe biggest problem remains Eq. (2). Because it is the foundation for every figure, the paper's numerical results cannot be trusted until the equation is repaired and the code or tables are made available. The energy typo alone would be minor; the undefined exit vector is not.\n\nNet: the paper is not ready as is. It shows a promising method, but the central formula needs a straightforward fix and the numerics need a pass. I would send a revised version to a referee. As is, I would not cite it or base any analysis on it.","headline":"Central Eq. (2) is degenerate as printed, but the underlying matrix-exponential method is promising enough to deserve a referee once the author fixes it.","tokens_in":6655,"tokens_out":3264,"would_cite":false,"duration_ms":33780,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Extreme-energy cosmic-ray origins can be located analytically: the photodisintegration cascade is a matrix-exponential distribution, yielding arbitrarily precise likelihoods over source distance and rigidity, with some observed charge…","keywords":["ultra-high-energy cosmic rays","extreme-energy cosmic rays","photodisintegration","matrix-exponential distributions","origin likelihood","Amaterasu event","cosmic-ray composition","stochastic interactions"],"falsifier":"Run a Monte Carlo propagation that allows fragments to lose energy between disintegrations (relaxing exact boost conservation) and compare the resulting two-dimensional origin-likelihood surface for an Amaterasu-like event; if the peak distances shift by more than a few megaparsecs relative to the matrix-exponential prediction, the method's resolution claim fails.","tokens_in":47,"feed_emoji":"🔭","tokens_out":7224,"duration_ms":204792,"temperature":0.7,"pith_summary":"Ultra-high-energy cosmic rays above $10^{20}$ eV are so strongly absorbed by the cosmic microwave background that any detected event must originate within a few tens of megaparsecs, but stochastic photodisintegration of the nuclear species makes pinpointing the source difficult. This paper claims that the full probability distribution for the distance and rigidity at the origin can be computed analytically as a matrix-exponential expression built from photodisintegration transition rates, requiring no Monte Carlo sampling. Applied to the Amaterasu event, the method yields arbitrarily precise normalized likelihood surfaces, resolving features that event-by-event Monte Carlo codes cannot constrain. A striking consequence is that for some observed charge groups, such as $2 < Z \\le 7$, the origin likelihood is nearly independent of whether the injected nucleus was silicon, calcium, or iron, suggesting that source localization may be possible without knowing the initial composition.","feed_headline":"Analytic method pinpoints extreme cosmic-ray origins","feed_subtitle":"Matrix-exponential likelihoods give arbitrarily precise source maps for events like Amaterasu, without Monte Carlo sampling.","key_machinery":"The load-bearing object is the matrix-exponential distribution $f(L) = \\alpha_0 \\exp(\\Lambda L) \\Lambda e$ for the distance to complete a photodisintegration cascade. $\\Lambda(\\gamma)$ is an upper-triangular matrix whose off-diagonal entries are the rates $\\lambda_{S_i \\to S_j}$ for one nucleus to disintegrate into another and whose diagonal entries are the total removal rates; the boost $\\gamma$ (energy per nucleon) is treated as constant because photodisintegration conserves it to good approximation, and the rigidity at the source follows from the observed energy and charge. For light products such as protons that are produced at multiple cascade steps, the same distribution is transformed by a 'rewards' construction that counts the yield along the chain. The likelihood surfaces in Figures 3–5 are these densities, normalized and restricted by a Gaussian on the observed energy, evaluated over source distance and rigidity.","core_discovery":"The central discovery is that the origin likelihood of an extreme-energy cosmic ray is a matrix-exponential distribution. Writing the photodisintegration cascade as a continuous-time Markov chain in which each species has an interaction rate $\\lambda(\\gamma)$ depending on the Lorentz boost $\\gamma = E/A$, the probability density for the distance $L$ needed to go from an initial species $S_0$ to a set of final species is $f(L) = \\alpha_0 \\exp(\\Lambda L) \\Lambda e$, where $\\Lambda$ is the matrix of transition rates and $\\alpha_0$ is the initial-state vector. Because photodisintegration conserves the boost almost exactly, the same $\\Lambda$ applies along the whole path, and the observed particle energy fixes the rigidity at the source through $R = E_{\\rm obs}/(A_{\\rm obs}\\eta)$. The paper evaluates these functions for three injected species ($^{28}$Si, $^{40}$Ca, $^{56}$Fe) and compares the resulting likelihood maps with the Monte Carlo–based analysis of the Amaterasu event, showing higher resolution and the ability to assign relative probabilities in regions where Monte Carlo has few or no events.","pith_inferences":["The composition-independence claim for $2 < Z \\le 7$ is stated for single injected species; a natural testable extension is to mix the injected species with weights from realistic source spectral fits and check whether the likelihood maps still coincide.","The single-boost Markov chain approximation ignores energy loss of protons via photopion production after they are produced mid-cascade; the author notes this is expected to be small, but including a distance-dependent boost for light fragments would show how much the likelihood peaks shift.","The matrix-exponential machinery may allow analytic marginalization over unknown source parameters such as the maximum rigidity $R_{\\rm max}$ or the injection spectrum, offering a direct route from event likelihoods to Bayesian source inference.","If the method's resolution survives the inclusion of energy losses and galactic magnetic fields, the same formalism could be applied to smaller statistical ensembles of sub-$10^{20}$ eV events for which the extragalactic magnetic field cannot be neglected."],"forward_implications":["Origin likelihoods for extreme-energy cosmic rays can be computed with arbitrary numerical precision, eliminating the sparse-event phase-space artifacts of Monte Carlo methods.","For observed charge groups such as $2 < Z \\le 7$, the origin likelihood is nearly independent of the injected composition, so source localization may not require knowing the initial nucleus.","Combined with an event-by-event composition measurement at the highest energies, the method can narrow the origin to within a few megaparsecs, potentially identifying a specific host galaxy.","The method reproduces the qualitative features of the Monte Carlo–based Amaterasu analysis while providing continuous, normalized likelihood values across the full distance–rigidity plane."],"supporting_citations":[{"why":"Supplies the matrix-exponential distributions for UHECR propagation that form the basis of Equation (2).","marker":"(Morejon, 2023)"},{"why":"Provides the rewards transform used to derive likelihoods for light products like protons.","marker":"(Bladt and Nielsen, 2017)"},{"why":"Reports the Amaterasu event and its energy, the target of the likelihood analysis.","marker":"(TA Collaboration, 2023)"},{"why":"Gives the Monte Carlo–based origin analysis and the Gaussian energy constraint used to restrict the likelihoods.","marker":"(Unger and Farrar, 2024b)"},{"why":"Describes the standard Monte Carlo propagation framework that the method is compared against.","marker":"(Batista et al., 2022)"},{"why":"Establishes the photodisintegration horizon for nuclei, motivating the local-volume restriction.","marker":"(Puget et al., 1976)"}],"fun_headline_variants":["Analytic method locates extreme cosmic-ray origins instantly","Matrix-exponential likelihoods yield exact cosmic-ray origins","No Monte Carlo: high-res origin maps for extreme cosmic rays","Exact source likelihoods for extreme cosmic rays"],"cache_read_input_tokens":8832,"weakest_assumption_plain":"The calculation assumes that, as a cosmic-ray nucleus breaks apart, every fragment continues with the same energy per nucleon as the original, so that one fixed set of interaction rates governs the whole journey from source to Earth.","fun_headline_variants_meta":{"raw":{"variants":["Analytic method locates extreme cosmic-ray origins instantly","Matrix-exponential likelihoods yield exact cosmic-ray origins","No Monte Carlo: high-res origin maps for extreme cosmic rays","Exact source likelihoods for extreme cosmic rays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000677,"raw_usage":{"total_tokens":3143,"prompt_tokens":1075,"completion_tokens":2068,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":691,"completion_tokens_details":{"reasoning_tokens":2004}},"tokens_in":691,"tokens_out":2068,"duration_ms":15968,"temperature":1.0,"reasoning_tokens":2004,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:55:29.427575+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a Monte Carlo propagation that allows fragments to lose energy between disintegrations (relaxing exact boost conservation) and compare the resulting two-dimensional origin-likelihood surface for an Amaterasu-like event; if the peak distances shift by more than a few megaparsecs relative to the matrix-exponential prediction, the method's resolution claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the matrix-exponential distributions for UHECR propagation that form the basis of Equation (2)."},{"cited_title":"and Nielsen, B","cited_arxiv_id":null,"evidence_quote":"Provides the rewards transform used to derive likelihoods for light products like protons."},{"cited_title":"An extremely energetic cosmic ray observed by a surface detector array","cited_arxiv_id":null,"evidence_quote":"Reports the Amaterasu event and its energy, the target of the likelihood analysis."},{"cited_title":"A., Tjus, J","cited_arxiv_id":null,"evidence_quote":"Describes the standard Monte Carlo propagation framework that the method is compared against."},{"cited_title":"L., Stecker , F","cited_arxiv_id":null,"evidence_quote":"Establishes the photodisintegration horizon for nuclei, motivating the local-volume restriction."}],"review_version":1}