{"id":"75e4de76-a143-44fc-8db8-ed4e25aeff68","arxiv_id":"2506.14659","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A jet-plus-corona model of NGC 1068, constrained by IceCube neutrinos, gives ALP-photon coupling limits g_aγ ≲ 7×10^-11 GeV^-1 for masses below 10^-9 eV, weaker than current bounds but offering a new multi-messenger route.","lead":"A team of physicists modeled the jets of the Seyfert galaxy NGC 1068 to explain its gamma-ray and neutrino emissions, then used that model to search for axion-like particles (ALPs). The resulting limits on the ALP-photon coupling are weaker than existing bounds, but the approach shows how multi-messenger data can sharpen new-physics searches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central limit is not reproducible: Sec. III never specifies how B0, n0, and Rem in the GammaALPs Jet model are derived from the fitted HadJet parameters, so the g_aγ bound carries an unquantified order-of-magnitude uncertainty from the propagation environment.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the unspecified mapping from the HadJet fitted parameters to the GammaALPs jet environment. This is a genuine reproducibility gap. Without the mapping, the central limit is not defined, and any shift in B0 or Rem changes the derived coupling bound approximately linearly (or quadratically in the perturbative regime). I considered whether the post-hoc bin selection or the hard-cut neutrino constraint was more fundamental, but those affect the intrinsic SED fit rather than the ALP propagation physics, and the paper's parameter ranges partially cover them. The internal inconsistency of the stated best-fit point (units swapped and value above the claimed limit) is suspicious but likely a typographical error; the missing mapping is a structural omission, not a typo. The paper uses a public, reproducible code (GammaALPs) and provides a large simulation set, which is a credit, but the missing input specification prevents independent reproduction of the headline limit. A CONDITIONAL verdict is therefore appropriate: the analysis may be sound once the mapping is stated and its sensitivity quantified, but as written the central claim is not verifiable.","tokens_in":13751,"tokens_out":14357,"duration_ms":135766,"concrete_test":"Ask the authors to supply the exact mapping from HadJet parameters (r0, σf, pe, fsc) to the GammaALPs inputs B0, n0, Rem used in Sec. III. Then recompute the 95% limit using the GammaALPs Jet model with B0 and n0 set to the best-fit values from Table I and Rem = 5 Rg, and compare to a run with B0 and n0 set to the extreme values of their derived ranges (B0 = 4709 G, n0 = 2.22e9 cm^-3 and B0 = 157 G, n0 = 4e7 cm^-3) at the same r0 and σf. If the limit shifts by more than 20% between runs, the missing mapping is decisive and the paper must either quantify this systematic or weaken its claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The limit g_aγ ≲ 7e-11 GeV^-1 is computed from ALP-photon attenuation in the jet magnetic field using the GammaALPs 'Jet' model (Sec. III, Eqs. 9-10), with B(r)=B0(r/Rem)^-1 and n(r)=n0(r/Rem)^-2. The paper does not state how B0, n0, and Rem are fixed from the HadJet parameters (r0, σf, pe, fsc) that were fitted to the Fermi-LAT and IceCube data. Table I reports B and n at the dissipation region for the best fit (B=196 G, n=1.013e8 cm^-3) but never connects these to the GammaALPs inputs, and Rem is not defined anywhere. Since the conversion probability scales as (g B L)^2 in the perturbative regime, an order-of-magnitude uncertainty in B or the propagation length translates directly into an order-of-magnitude shift in the derived coupling limit. The marginalization over θ_astro in Sec. III.A varies only the intrinsic SED, not the propagation environment; if B0, n0, Rem are held fixed at best-fit values, the limit ignores astrophysical uncertainty in the jet magnetic field. If they are instead varied with r0 and σf, the mapping is absent. This is the least secure condition for the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models the lepto-hadronic jet emission of the Seyfert galaxy NGC 1068, fitting the Fermi-LAT gamma-ray spectrum subject to an IceCube neutrino-band selection, and then uses the GammaALPs \"Jet\" model to compute ALP-photon conversion in the jet magnetic field. Marginalizing over four astrophysical jet parameters, it derives a 95% CL upper limit g_aγ <~ 7e-11 GeV^-1 for m_a <~ 1e-9 eV, and also presents projected limits for improved Fermi-LAT precision and for a population of ten NGC 1068-like sources. The paper is explicit that the best-fit ALP point is only 1.02 sigma and that the derived limits are weaker than existing astrophysical bounds, so the main claim is the demonstration that multi-messenger modeling can provide a new, complementary ALP probe.","tokens_in":14127,"tokens_out":4044,"duration_ms":42713,"significance":"If the result is robust, the paper makes a useful methodological contribution: it uses a concrete, physically motivated jet+corona model for NGC 1068, scans 10^6 parameter combinations, and marginalizes over four astrophysical nuisance parameters in the ALP limit, which is good practice. The paper also honestly reports that the best-fit ALP point is not significant and that the limits are weaker than existing constraints. However, the numerical value of the central limit is not reproducible from the manuscript as written, because the GammaALPs propagation environment (B0, n0, Rem) is never connected to the fitted HadJet parameters. Since the conversion probability scales roughly as (g_aγ B L)^2 in the perturbative regime, this missing mapping leaves an order-of-magnitude uncertainty in the derived coupling limit. The central claim is therefore defensible only after the propagation-environment mapping is specified and its uncertainty is propagated.","major_comments":[{"comment":"The GammaALPs Jet model requires the magnetic field profile B(r) = B0 (r/Rem)^-1 and electron density profile n(r) = n0 (r/Rem)^-2, but the paper never states how B0, n0, and Rem are fixed from the HadJet parameters (r0, sigma_f, pe, fsc) that were fitted to Fermi-LAT and IceCube data. Table I reports B = 196 G and n = 1.013e8 cm^-3 \"at the dissipation region\" for the best fit, but it does not say that these are the GammaALPs B0 and n0, and Rem is not given a numerical value anywhere. Because the ALP-photon conversion probability depends on (g_aγ B L)^2 in the perturbative regime, an unspecified normalization or propagation length translates directly into an unquantified shift of the derived limit, potentially by an order of magnitude. Please specify the exact mapping, state the value of Rem used, and either marginalize over the propagation-environment parameters or demonstrate explicitly that the limit is insensitive to them.","section":"Section III, Eqs. (9)-(10) and Table I"},{"comment":"The IceCube neutrino data are used only as a hard selection: models are kept if their predicted neutrino flux falls inside the observed 79+22/-20 event band, but the neutrino information does not enter the chi-square in Eq. (3) or the profile likelihood in Sec. III.A. This discards information and can bias the selected region of astrophysical parameter space. The authors should either include a Poisson likelihood term for the IceCube event count in the profile likelihood or demonstrate that the final ALP limit is robust to the chosen band-cut criterion.","section":"Section II.B and II.C"},{"comment":"The Fermi-LAT fit uses only the first 8 of the 14 energy bins, with the stated motivation that higher-energy bins receive a starburst contribution that is not modeled. The best-fit SED and therefore the derived ALP limit depend on this post-hoc bin selection, yet no test of this dependence is shown. Please show how the limit changes with the number of included bins or include a starburst template, since this choice directly affects the baseline spectrum against which ALP attenuation is constrained.","section":"Section II.B"}],"minor_comments":[{"comment":"The reported best-fit point has its units interchanged: it should read (m_a, g_aγ) = (2.2e-9 eV, 4.6e-10 GeV^-1), not (2.2e-9 GeV^-1, 4.6e-10 eV).","section":"Section III.B"},{"comment":"The text attributes the Perseus Galaxy bound to reference [28], which is actually the MAGIC paper on NGC 1068; the Perseus bound appears to be reference [19]. Please correct the citation.","section":"Section III.B and Fig. 3"},{"comment":"Equation (13) writes delta_i ~ N(0,1), while the surrounding text states that the perturbations have standard deviation sigma = 0.1; please use N(0, 0.1) or otherwise clarify the notation.","section":"Section III.B, Eq. (13)"},{"comment":"The symbol Rem is introduced as the \"emitting region\" but is never related to zdiss = 5 Rg or to the corona radius Rcor = 50 Rg; a definition and numerical value would make the propagation calculation reproducible.","section":"Section III, Eqs. (9)-(10)"},{"comment":"The row for pe contains a formatting artifact (\"log(dn/dE)\" over \"log(dE)\"); please define the power-law index cleanly, e.g., dn/dE ∝ E^-pe.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the missing mapping between the fitted HadJet parameters and the GammaALPs Jet-model environment; this is fixable but is load-bearing for the quoted limit. I would encourage the editor to request, as part of the revision, a precise statement of the mapping and ideally the relevant parameter files or code, since the current manuscript does not allow the central numerical result to be reproduced. I do not see evidence of circularity: the ALP coupling is scanned rather than fitted, and the best-fit point is only 1.02 sigma."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThe paper does something genuinely new: it takes the IceCube neutrino data on NGC 1068 as a hard anchor for the astrophysical baseline and then asks what ALP-photon conversion in the jet would do to the gamma-ray spectrum. The HadJet extension with corona targets is reasonable, and the authors are honest that their limit g_aγ ≲ 7e-11 is weaker than Mrk 421, NGC 1275, or Perseus. The main value is the demonstration that a neutrino-bright Seyfert can serve as a controlled environment for ALP searches.\n\nWhat the paper does well: it fits the full SED, marginalizes over four astrophysical parameters, uses IceCube to filter jet models rather than fitting them, and checks that the best-fit model gives 92 muon events against the observed 79. The profile-likelihood treatment and the projection to future sources with 10% precision are sensible. The limit itself is not circular in any obvious way, since the ALP coupling is scanned and the best-fit point is only 1.02σ.\n\nThe soft spots are real but not fatal. The biggest one is the missing mapping between the HadJet parameters and the GammaALPs propagation environment. The paper never states how B0, n0, and Rem in Eqs. (9–10) are fixed from (r0, σf, pe, fsc). Table I gives B and n at the dissipation region; if the intended mapping is B0=B(zdiss), n0=n(zdiss), Rem=zdiss, that should be written down explicitly. And the propagation environment should be varied along with the nuisance parameters, otherwise the marginalization misses the magnetic-field uncertainty, which in turn shifts the limit by potentially an order of magnitude. The post-hoc choice of the first 8 Fermi-LAT bins is a bit ad hoc, and the IceCube constraint is applied as a hard band cut rather than a full likelihood term. Neither of those changes the qualitative conclusion, but they should be documented as limitations.\n\nThe stress-test note about the missing mapping is fair, though I would not call the paper irreproducible in the strong sense. An expert can guess the intended mapping, but the paper should state it. The central limit is plausible, and the authors themselves downplay it; that honesty earns credit.\n\nWho is this for? Someone working on ALP constraints from gamma-ray sources, or on multi-messenger modeling of Seyferts. It deserves a serious referee: the methodology is new, the data handling is mostly careful, and the implications for future sources are worth airing.\n\nRecommendation: send to peer review; ask the authors to clarify the GammaALPs mapping, extend the marginalization to the propagation environment, and discuss the bin selection. After that, it could be a solid contribution.\n\nBest,\n\n[Your name]","headline":"A methodologically interesting but weak ALP limit from NGC 1068; the central number is credible but the propagation-environment mapping needs to be made explicit.","tokens_in":14591,"tokens_out":3203,"would_cite":true,"duration_ms":31518,"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":"Joint Fermi-LAT and IceCube data from NGC 1068 cap the axion–photon coupling at 7×10⁻¹¹ GeV⁻¹ for axion masses below 10⁻⁹ eV.","keywords":["axion-like particles","ALP-photon coupling","NGC 1068","multi-messenger astrophysics","high-energy neutrinos","gamma-ray attenuation","Seyfert galaxy","jet model"],"falsifier":"Measure NGC 1068's gamma-ray spectrum between 100 GeV and 1 TeV with percent-level precision: the predicted ALP survival probability for $g_{a\\gamma}=7\\times10^{-11}$ GeV$^{-1}$ would produce a specific energy-dependent attenuation, and its absence would falsify the claimed limit. A direct determination of the jet magnetic-field normalization from Faraday rotation or synchrotron self-absorption would test the input that sets the bound.","tokens_in":13568,"feed_emoji":"🔭","tokens_out":10163,"duration_ms":91604,"temperature":0.7,"pith_summary":"This paper uses the joint gamma-ray and neutrino observations of the Seyfert galaxy NGC 1068 to build a new probe of axion-like particles. It models the jet near the supermassive black hole with a lepto-hadronic emission model that reproduces both the Fermi-LAT spectrum and the IceCube neutrino excess, then lets gamma rays oscillate into ALPs in the jet's magnetic field. The result is a 95% upper limit of $g_{a\\gamma} \\lesssim 7\\times10^{-11}$ GeV$^{-1}$ for $m_a \\lesssim 10^{-9}$ eV, after marginalizing over the astrophysical jet parameters. The limit is somewhat weaker than existing bounds, but it demonstrates that multi-messenger sources can disentangle intrinsic emission from new-physics attenuation, and projects that ten such sources at 10% precision would match the best current constraints.","feed_headline":"NGC 1068 caps axion coupling at 7×10⁻¹¹ GeV⁻¹","feed_subtitle":"Gamma-ray and neutrino data together pin down the ALP coupling; ten similar galaxies would rival the best limits.","key_machinery":"The load-bearing object is the ALP–photon mixing matrix that couples the two photon polarizations with the ALP field, integrated over a jet magnetic field that decreases as $B(r)=B_0(r/R_{\\rm em})^{-1}$ and an electron density $n(r)=n_0(r/R_{\\rm em})^{-2}$ across the emission region. The intrinsic spectrum entering the mixing calculation comes from a lepto-hadronic jet model whose four parameters—jet base radius, magnetization at the dissipation region, particle power-law index, and acceleration efficiency—are fit to Fermi-LAT data subject to the IceCube neutrino constraint, then marginalized over when deriving the ALP limit. The mixing matrix yields the photon survival probability as a function of energy, axion mass, and coupling, which is what converts the observed spectrum into a bound.","core_discovery":"The central claim is that ALP–photon oscillations in the magnetic field of the inner jet of NGC 1068 would over-attenuate its gamma-ray spectrum for couplings above $g_{a\\gamma}\\simeq7\\times10^{-11}$ GeV$^{-1}$ when $m_a\\lesssim10^{-9}$ eV, at 95% confidence. To reach this, the paper constructs a source model in which protons accelerated in the jet collide with corona particles, producing the gamma rays observed by Fermi-LAT and the neutrinos observed by IceCube; the model predicts 92 muon neutrino events, consistent with the reported $79^{+22}_{-20}$. This source model fixes the intrinsic spectrum, so the ALP search looks for extra energy-dependent attenuation of that spectrum. No significant ALP preference is found (best-fit at $1.02\\sigma$), so the paper reports an upper limit rather than a detection.","pith_inferences":["The same marginalization procedure could be applied to other neutrino-bright Seyferts such as NGC 7469 once their gamma-ray spectra are available; the paper only sketches this possibility.","Because the conversion probability scales with $B^2$ and propagation length, an independent handle on the jet magnetic field from radio interferometry or spectral breaks would sharpen the limit considerably.","Applying the same pipeline to the starburst component of NGC 1068 at higher energies, which the paper intentionally excludes, could probe heavier ALP masses and a different magnetic environment.","If future Cherenkov telescopes see no spectral irregularities in NGC 1068, the combined dataset would push the bound below $10^{-11}$ GeV$^{-1}$, into the region where ALPs could account for dark matter."],"forward_implications":["ALPs with $m_a\\lesssim10^{-9}$ eV and $g_{a\\gamma}\\gtrsim7\\times10^{-11}$ GeV$^{-1}$ would erase the NGC 1068 gamma-ray flux, so the joint Fermi-LAT and IceCube data exclude them.","The IceCube neutrino measurement is what breaks the degeneracy: the source model must simultaneously account for gamma rays and neutrinos, leaving no room to hide ALP attenuation in a reshaped intrinsic spectrum.","A population of ten NGC 1068-like sources measured at 10% precision would push ALP limits from gamma-ray sources to the level of current best constraints, making multi-messenger galaxies a competitive probe.","The absence of a significant ALP preference (best fit at $1.02\\sigma$) means the result is a clean upper limit, not evidence for new physics."],"supporting_citations":[{"why":"Supplies the base jet model with the magnetic and density profiles used as the ALP propagation environment.","marker":"[51]"},{"why":"Extends the jet model with hadronic processes that produce the IceCube neutrinos and the pp/pγ gamma rays.","marker":"[52]"},{"why":"Provides the IceCube neutrino measurement that the source model must reproduce, anchoring the multi-messenger constraint.","marker":"[27]"},{"why":"Supplies the Fermi-LAT gamma-ray data used in the fit.","marker":"[29]"},{"why":"Provides the dedicated Fermi-LAT analysis of NGC 1068 whose energy bins are reused in the fit.","marker":"[34]"},{"why":"Computes the photon survival probability in the jet environment, the core of the ALP limit.","marker":"[63]"},{"why":"Establishes the ALP-photon mixing formalism that the paper applies.","marker":"[6]"},{"why":"Sets the corona density, optical depth, and extent that determine the target for hadronic collisions.","marker":"[33]"},{"why":"Provides the likelihood-ratio criterion used for the 95% CL limit.","marker":"[68]"}],"fun_headline_variants":["NGC 1068 multi-messenger data cap axion-photon coupling","Neutrino and gamma-ray data together bound axion coupling","Axion coupling limited by NGC 1068's neutrinos and gamma rays","No axion signal from NGC 1068, limit set at 7e-11 GeV^-1","Multi-messenger NGC 1068 gives axion coupling ceiling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The limit rests on the assumed power-law radial profile of the jet magnetic field and its unstated normalization; if the true field inside the emission region differs, the coupling bound would shift.","fun_headline_variants_meta":{"raw":{"variants":["NGC 1068 multi-messenger data cap axion-photon coupling","Neutrino and gamma-ray data together bound axion coupling","Axion coupling limited by NGC 1068's neutrinos and gamma rays","No axion signal from NGC 1068, limit set at 7e-11 GeV^-1","Multi-messenger NGC 1068 gives axion coupling ceiling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000684,"raw_usage":{"total_tokens":3123,"prompt_tokens":985,"completion_tokens":2138,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":601,"completion_tokens_details":{"reasoning_tokens":2036}},"tokens_in":601,"tokens_out":2138,"duration_ms":14681,"temperature":1.0,"reasoning_tokens":2036,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:49:10.229185+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure NGC 1068's gamma-ray spectrum between 100 GeV and 1 TeV with percent-level precision: the predicted ALP survival probability for $g_{a\\gamma}=7\\times10^{-11}$ GeV$^{-1}$ would produce a specific energy-dependent attenuation, and its absence would falsify the claimed limit. A direct determination of the jet magnetic-field normalization from Faraday rotation or synchrotron self-absorption would test the input that sets the bound.","supporting_citations":[{"cited_title":"Lucchini, C","cited_arxiv_id":null,"evidence_quote":"Supplies the base jet model with the magnetic and density profiles used as the ALP propagation environment."},{"cited_title":"Kantzas, S","cited_arxiv_id":null,"evidence_quote":"Extends the jet model with hadronic processes that produce the IceCube neutrinos and the pp/pγ gamma rays."}],"review_version":2}