{"id":"80c2f676-23a4-428a-b1f5-428e15669f10","arxiv_id":"2504.13639","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Starburst-galaxy black hole mergers with precessing jets may dominate the ultra-high-energy cosmic ray flux, and a flux-density-plus-luminosity ranking can identify the likely active sources.","lead":"The paper argues that merging stellar-mass black holes in starburst galaxies, with precessing jets, produce ultra-high-energy cosmic rays, neutrinos, and gravitational waves. It proposes a two-step way to find these sources: pick the brightest infrared galaxy in a given flux bin, because it is most likely to be the active one.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative basis for the claimed UHECR dominance and for the sample-size formula rests on a single candidate source, 41.9+58 in M82, interpreted as a precessing stellar-mass BBH jet; if that identification is wrong or does not scale with FIR luminosity, the central claim and the two-step…","rationale":"I read the paper as proposing a speculative but coherent scenario: stellar-mass BBH mergers with spin-flip precessing jets in star-forming galaxies produce episodic UHECR and neutrino emission, and because the event rate is assumed proportional to FIR luminosity, the highest-luminosity galaxy in a flux-density bin is the most probable host. The reader's weakest-assumption diagnosis matches my own: the entire rate and duty-cycle normalization reduces to one candidate source in M82, whose physical interpretation as a precessing-jet BBH merger is not established beyond the authors' own prior work. The Poisson-probability issue in Section 3.8 is real but secondary: it changes the sample size by a factor of a few, whereas a misidentification of 41.9+58 would remove the empirical anchor completely. I agree with the CONDITIONAL verdict because the concern is concrete, testable, and does not require rejecting the paper's logical structure; it requires better empirical support. No change to the reader's verdict is needed.","tokens_in":16780,"tokens_out":5470,"duration_ms":54886,"concrete_test":"Take the available multi-epoch VLBI images of 41.9+58 in M82 and quantify proper motion, spectral-index changes, and polarization structure against the predictions of (a) a precessing relativistic jet from a coalescing stellar-mass black-hole binary and (b) a supernova remnant expanding in a density gradient. If the observed morphology and time evolution are consistent with a shell-like SNR rather than a collimated precessing jet, the rate normalization of one per 2,500 years (Section 3.1) is unsupported and the Section 3.8 sample-size formula must be re-derived or abandoned. As a secondary check within the same analysis, verify whether the source's radio luminosity and size are compatible with the jet-power consistency argument in Section 3.4.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim that starburst galaxies may dominate the UHECR flux, and the quantitative sample-size rule in Section 3.8, are anchored in Section 3.1 by identifying one compact radio source, 41.9+58 in M82, as a second-generation stellar-mass black-hole merger with precessing jets, yielding one merger per 2,500 years. The authors themselves state this rate has 'an error range of probably at least an order of magnitude' (Section 3.2), and the identification is not independently supported: the double-cone morphology could in principle be an expanding supernova remnant in a stratified medium or a jet from a single accreting black hole, rather than a coalescing binary. If the identification is wrong or atypical, the derived rate, the duty-cycle probability proportional to L_FIR, and the 'approximately 100 percent expectation' sample size (tau_3.4/tau_ep = 2,500/5) all lose their quantitative basis. The claim also requires that the event rate scales linearly with FIR luminosity across galaxies and cosmic epochs; this is assumed, not tested, and no comparison is made with the independent LIGO/Virgo merger-rate density or with the absence of analogous double-cone radio sources in other starbursts. This is a load-bearing empirical anchor, not a minor parameter uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that stellar-mass black-hole mergers in starburst and normal galaxies, when preceded by spin-flip jet precession, act as episodic accelerators of ultra-high-energy cosmic rays and neutrinos, and that the same mergers eventually produce gravitational waves. The authors generalize an earlier proposal based on the compact radio source 41.9+58 in M82 to the cosmic population of starburst galaxies, and argue that these sources may dominate over AGN in the observed UHECR flux. They then outline a two-step observational strategy: within a FIR flux-density interval, rank candidate host galaxies by FIR luminosity, under the assumption that the probability of being currently active is proportional to FIR luminosity. The paper presents analytic estimates for binary orbital angular-momentum loss via magnetic winds, a rate normalization of one merger per 2,500 years in M82, an episode duration of about 5 years, and a sample-size threshold in Section 3.8 based on the ratio of these timescales.","tokens_in":17134,"tokens_out":9133,"duration_ms":85310,"significance":"If the central assumptions hold, the paper offers a simple, falsifiable prioritization rule for multi-messenger follow-up: at fixed FIR flux density, select the highest-FIR-luminosity galaxies because they have the largest duty-cycle probability of harboring an active episodic accelerator. This is a testable prediction that could be confronted with future neutrino and gravitational-wave events, and the authors are transparent about some of the order-of-magnitude uncertainties. The manuscript does not supply machine-checked proofs or code, but it does make its quantitative skeleton explicit enough to be challenged. The significance is currently limited by the fact that the rate normalization, the FIR-luminosity scaling, and the near-maximal-spin fraction are all anchored to one morphological candidate and a series of undocumented assumptions, so the claimed UHECR dominance over AGN is best read as a conditional scenario rather than an established result.","major_comments":[{"comment":"The entire quantitative framework, including the 2,500-year repetition time used in the Section 3.8 sample-size inequality, is normalized to a single compact radio source, 41.9+58 in M82, interpreted as a second-generation stellar-mass black-hole merger with precessing jets. The paper itself states that the resulting rate is 'a very uncertain estimate' and has 'an error range of probably at least an order of magnitude.' The only discriminator offered is a qualitative morphological argument against a stratified-atmosphere explosion, which does not exclude other jet-producing configurations such as an accreting black hole in a binary or an unusual supernova remnant. Because this one source carries the rate normalization for the claimed UHECR dominance and for the strategy's sample-size rule, the issue is load-bearing. The authors should either provide a quantitative precessing-jet model that reproduces the radio morphology and spectrum of 41.9+58 and rules out alternatives, or demonstrate that analogous double-cone sources appear in a FIR-luminosity-selected sample of nearby starbursts with the expected rate, or cross-check the implied merger rate per unit stellar mass against the LIGO/Virgo merger-rate density.","section":"Sections 3.1, 3.2, and 3.8"},{"comment":"The inequality Σ LFIR,j,i > L_FIR,M82 τ3.4/τep is accompanied by the statement that 'there is a ≃ 100 percent expectation, that some galaxy is in an active phase of an episode.' This conflates an expected number with a probability. If each galaxy's active probability is p_i = (τep/τ3.4) (LFIR,i/L_FIR,M82) and the sum of p_i equals unity, the Poisson probability of at least one active galaxy is only 1 - e^{-1} ≈ 63 percent, not approximately 100 percent. The sample-size rule therefore does not provide the claimed confidence unless the duty-cycle interpretation is revised or additional assumptions are stated explicitly.","section":"Section 3.8"},{"comment":"The conclusion that stellar-mass black-hole mergers in starburst and normal galaxies 'may dominate over AGN' in the UHECR flux assumes that a large fraction of massive-star black holes are born near maximally spinning, that magnetic winds efficiently shrink the orbit, and that relativistic jets form with an efficiency of about one-third. The manuscript itself states that the fraction of massive stars that produce such black holes 'is unknown at present' (Section 2.1), and the 1/3 efficiency is an adopted value with no propagated uncertainty. Because this unknown fraction multiplies the entire energy-budget comparison with AGN, the dominance claim is an upper-limit scenario rather than a demonstrated inference. The authors should either derive a plausible range for the fraction from stellar-evolution and binary-population calculations, or explicitly rewrite the claim as a conditional statement and propagate the fraction through the comparison with the AGN contribution.","section":"Sections 2.1 and 3.4"}],"minor_comments":[{"comment":"The statement 'The Gravitational Wave (GW) data support such a scenario' is not backed by any quantitative comparison in the text; either add the specific GW observables being referenced or rephrase as 'are not inconsistent with such a scenario.'","section":"Abstract"},{"comment":"The displayed formula for Jorb, Jorb = π^{-1/3} M^{5/3} G^{2/3} P^{1/3}, differs from the standard equal-mass circular orbital angular momentum by a numerical factor of 4^{1/3} when r is half the binary separation; the subsequent conclusion for rdot/r is insensitive to this constant, but the formula should be corrected for internal consistency.","section":"Section 2.1"},{"comment":"The example '10^12/10^10.6 × 1/2,500 per year, so about 1 in about 100 years' should explicitly state that this is a present-day, local-universe estimate that excludes efficiency factors and redshift evolution, so that readers do not mistake it for a cosmic rate.","section":"Section 3.3"},{"comment":"The sentence 'at any given flux density the sources with the highest luminosity, so highest redshift, have the highest probability to contribute' assumes that FIR luminosity is the correct proxy for the currently-active probability and that no strong redshift-dependent evolution of the near-maximal-spin fraction is present; these are assumptions that should be listed as caveats.","section":"Section 3.6"},{"comment":"The sentence 'This is fully consistent with new Auger results [Auger-Coll. (2024)]' is unsupported as written; no quantitative confrontation of the model with the Auger spectrum or composition constraints appears in the paper, so this claim should be substantiated or removed.","section":"Section 3.4"},{"comment":"Several references and headings contain typographical artifacts, such as 'GA-NIFS: JWST disc overs an offset AGN' and the version-history line at the top of the manuscript, which should be cleaned before publication.","section":"References and front matter"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a speculative but clearly formulated proposal whose quantitative conclusions rest on a single morphological candidate and a series of unvalidated scaling assumptions. I would encourage the editor to require the authors to either add an external cross-check of the M82-based rate or explicitly reframe the paper as a conditional strategy whose success probability is not yet quantified. Correcting the Section 3.8 expectation-value error is essential, as is softening the abstract's claim that GW data support the scenario without a specific citation or calculation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the two-step selection strategy. Take a flux-density-limited sample of galaxies, then rank by FIR luminosity, because at fixed flux density the higher-luminosity objects are at higher redshift and have a higher duty-cycle probability of being active right now. That is clean, implementable, and independent of the details of the emission mechanism. It is a useful tool for anyone doing neutrino or gravitational-wave follow-up.\n\nThe paper earns credit for candor. The authors explicitly flag the M82 rate as one per 2,500 years with 'probably at least an order of magnitude' error, and they do not pretend the identification of 41.9+58 as a second-generation stellar-mass BH merger is airtight. They cite the alternative interpretations and explain why they find them unlikely. The step from one galaxy to a cosmic population is natural, and the proposal to use FIR luminosity as a proxy for episode rate is testable.\n\nThe soft spots are in the quantitative superstructure. The load-bearing anchor is one source in M82. If 41.9+58 is an expanding supernova remnant in a stratified medium or a single accreting BH jet, the whole rate normalization collapses. The scaling of event rate with FIR luminosity is assumed, not demonstrated; there is no comparison with the LIGO/Virgo merger-rate density or with the fact that no other starburst shows a similar double-cone source. Section 3.8 claims a '~100 percent expectation' that a sample contains an active source when the expectation value is one; the probability of at least one is about 63 percent, not 100. That is a fixable error but it inflates the confidence. And the paper never compares its predicted neutrino flux to IceCube stacking limits, so the claim that starbursts may dominate the UHECR flux is a possibility, not a result.\n\nThe citation pattern is heavy on the group's own prior work, but that is because the model is developed in that series; the strategy itself is generic and does not require accepting the M82 interpretation.\n\nFor whom: multi-messenger astronomers who need a cheap way to rank candidate host galaxies. The two-step method is worth having even if the specific merger scenario is wrong. It deserves a serious referee. I would ask for a revision that fixes the Poisson statement, adds the LIGO/Virgo and IceCube comparisons, and scales back the dominance claim to match the evidence.","headline":"A testable two-step selection strategy for neutrino/GW host candidates, built on a speculative single-source rate anchor; the strategy itself is worth engaging, but the quantitative claims need revision.","tokens_in":17737,"tokens_out":3305,"would_cite":true,"duration_ms":29058,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Starburst galaxies may dominate the highest-energy cosmic rays","keywords":["neutrinos","starburst galaxies","black hole mergers","gravitational waves","particle acceleration","ultra-high-energy cosmic rays","multi-messenger astronomy"],"falsifier":"Repeated high-resolution radio imaging of 41.9+58 should show its double-cone structure changing or rotating on timescales of years if it is an active precessing jet; a static structure, or a second such source in M82 appearing at a rate incompatible with one per 2,500 years, would break the rate anchor. On the population side, if a sample of identified neutrino or gravitational-wave hosts shows no preference for the highest far-infrared luminosity galaxy within each flux-density interval, the probability ranking is refuted.","tokens_in":16580,"feed_emoji":"🌌","tokens_out":11560,"duration_ms":87919,"temperature":0.7,"pith_summary":"The paper proposes that mergers of stellar-mass black holes in starburst and normal galaxies, whose precessing jets sweep through gas, are episodic factories for ultra-high-energy cosmic rays, high-energy neutrinos, and gravitational waves. It argues that these galaxies may outshine active galactic nuclei in the flux of ultra-high-energy particles reaching Earth, and that they feed both the cosmic neutrino background and the gravitational-wave background. The practical centerpiece is a two-step identification strategy: for any detected neutrino or gravitational-wave event, rank the candidate host galaxies by far-infrared flux density, then within that sample rank by far-infrared luminosity, because at a given flux density the most luminous galaxy has the highest probability of being active right now. If the claim holds, multi-messenger searches gain a concrete way to assign otherwise orphaned events to a specific galaxy.","feed_headline":"Starburst galaxies may dominate the highest-energy cosmic rays","feed_subtitle":"A two-step infrared ranking could say which galaxy produced any neutrino or gravitational-wave event.","key_machinery":"The central object is a precessing relativistic jet pair produced by two coalescing stellar-mass black holes with misaligned spins: as the spins align before merger, the jet direction sweeps through a cone and continuously encounters fresh interstellar gas to accelerate. The load-bearing identity is the rate budget that connects the M82 anchor to the cosmic population: the inter-event time scales inversely with far-infrared luminosity, so the probability that a galaxy is active in the observer frame is proportional to its far-infrared luminosity, while the probability of detecting it is proportional to its flux density. This ratio, expressed as the repetition time of 2,500 years divided by the episode duration of 5 years, sets the sample size required to have near-certainty that at least one galaxy in the flux-density interval is currently active. That machinery turns the concept of episodic acceleration into a concrete two-step ranking procedure.","core_discovery":"The discovery claim is that stellar-mass black hole mergers with spin-flips, where two black holes with misaligned spins spiral together and their jets precess, accelerate particles to ultra-high energies, and that such events recur throughout cosmic star formation. A single compact radio source, 41.9+58 in M82, is interpreted as one such merger, yielding an estimated rate of one per 2,500 years in that galaxy; scaling this rate with far-infrared luminosity gives one active episode roughly every 100 years in the most luminous starbursts, with each acceleration episode lasting about 5 years. Because the probability that a galaxy hosts an active episode right now scales with its luminosity while detection probability scales with flux density, the highest-luminosity galaxies at any given flux density are the most likely sources. The paper concludes that starburst and normal galaxies may dominate over active galactic nuclei in the ultra-high-energy cosmic-ray flux observed at Earth, and that ranking candidates first by far-infrared flux density and then by far-infrared luminosity is the correct search strategy for identifying the hosts of neutrino and gravitational-wave events.","pith_inferences":["A direct test of the strategy's rate anchor is to survey a volume-limited sample of starburst galaxies for double-cone radio structures; the fraction with one active source should track far-infrared luminosity if the model is right.","The same flux-density-then-luminosity ranking could be applied to other episodic transients whose rate scales with a galaxy luminosity proxy, such as core-collapse supernovae, a case the paper notes would give the same conclusion.","If stellar-mass mergers dominate ultra-high-energy cosmic rays, the arrival directions and composition of these particles should correlate with the cosmic star-formation history rather than with the distribution of active galactic nuclei; next-generation observatories could test this.","The sample-size formula could be inverted: recording how often the highest-luminosity galaxy is confirmed as the host would measure the true inter-event time, turning the search strategy into a population-rate measurement."],"forward_implications":["For a detected high-energy neutrino or gravitational-wave event, ranking candidate galaxies first by far-infrared flux density and then by far-infrared luminosity gives a concrete, ordered list of host candidates, with the top-ranked galaxy the most probable source.","Starburst and normal galaxies would join active galactic nuclei, and may dominate them, as contributors to the ultra-high-energy cosmic-ray flux observed at Earth, which would reshape models of cosmic-ray origin and composition.","The cosmic neutrino background and the gravitational-wave background would receive a previously unquantified contribution from stellar-mass black hole mergers, with episodes lasting about 5 years and repeating on timescales of roughly 100 years in the most luminous starbursts.","Because neutrino emission is relativistically boosted along the jet while gravitational waves are emitted in all directions, a single merger episode will rarely be detectable in both messengers at once; searches should treat neutrino and gravitational-wave samples as complementary rather than requiring coincidence."],"supporting_citations":[{"why":"Discovers the 43 compact radio sources in M82, including the double-cone source 41.9+58 that anchors the rate estimate.","marker":"Kronberg et al. (1985)"},{"why":"Provides radio spectra of the M82 compact sources, supporting the interpretation of 41.9+58 as a jet-driven structure rather than an ordinary supernova.","marker":"Allen & Kronberg (1998)"},{"why":"Develops the earlier proposal that 41.9+58 is a stellar-mass black hole merger with precessing jets and supplies the one-per-2,500-year rate and power estimates.","marker":"Biermann et al. (2018)"},{"why":"Supplies the precession and inspiral timescales used to set the roughly 5-year duration of the acceleration episode.","marker":"Gergely & Biermann (2009)"},{"why":"Documents the high multiplicity of massive stars, the population basis for second-generation mergers with misaligned spins.","marker":"Chini et al. (2012)"},{"why":"Provides the magnetic-wind angular momentum loss formalism used to show how binary stars can be driven together into close orbits.","marker":"Weber & Davis (1967)"},{"why":"Gives the minimum jet power of radio-loud quasars, used as a consistency check for the jet power derived from the M82 source.","marker":"Punsly & Zhang (2011)"},{"why":"Presents the ultra-high-energy cosmic-ray arrival-direction, spectrum, and composition data that the paper says are consistent with a starburst-dominated contribution.","marker":"Auger-Coll. (2024)"}],"fun_headline_variants":["Starbursts may beat AGN for ultra-high-energy cosmic rays","Two-step galaxy ranking finds neutrino and GW sources","Spin-flip black hole mergers in starbursts power cosmic rays","Luminosity and flux density rank likely neutrino hosts","Episodic starburst sources may dominate cosmic-ray flux"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire rate normalization rests on reading one compact radio source in M82 as a precessing-jet merger of two black holes and assuming that such events become more frequent in proportion to a galaxy's far-infrared luminosity; if that reading is wrong or atypical, the claimed dominance and the search strategy lose their anchor.","fun_headline_variants_meta":{"raw":{"variants":["Starbursts may beat AGN for ultra-high-energy cosmic rays","Two-step galaxy ranking finds neutrino and GW sources","Spin-flip black hole mergers in starbursts power cosmic rays","Luminosity and flux density rank likely neutrino hosts","Episodic starburst sources may dominate cosmic-ray flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001033,"raw_usage":{"total_tokens":4358,"prompt_tokens":960,"completion_tokens":3398,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":3313}},"tokens_in":576,"tokens_out":3398,"duration_ms":24156,"temperature":1.0,"reasoning_tokens":3313,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:03:12.170695+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeated high-resolution radio imaging of 41.9+58 should show its double-cone structure changing or rotating on timescales of years if it is an active precessing jet; a static structure, or a second such source in M82 appearing at a rate incompatible with one per 2,500 years, would break the rate anchor. On the population side, if a sample of identified neutrino or gravitational-wave hosts shows no preference for the highest far-infrared luminosity galaxy within each flux-density interval, the probability ranking is refuted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Discovers the 43 compact radio sources in M82, including the double-cone source 41.9+58 that anchors the rate estimate."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides radio spectra of the M82 compact sources, supporting the interpretation of 41.9+58 as a jet-driven structure rather than an ordinary supernova."},{"cited_title":"Supernova explosions of massive stars and cosmic rays","cited_arxiv_id":"1803.10752","evidence_quote":"Develops the earlier proposal that 41.9+58 is a stellar-mass black hole merger with precessing jets and supplies the one-per-2,500-year rate and power estimates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the high multiplicity of massive stars, the population basis for second-generation mergers with misaligned spins."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the magnetic-wind angular momentum loss formalism used to show how binary stars can be driven together into close orbits."},{"cited_title":"& Zhang, S","cited_arxiv_id":null,"evidence_quote":"Gives the minimum jet power of radio-loud quasars, used as a consistency check for the jet power derived from the M82 source."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents the ultra-high-energy cosmic-ray arrival-direction, spectrum, and composition data that the paper says are consistent with a starburst-dominated contribution."}],"review_version":1}