{"id":"5d45328f-87fa-40c3-9816-70cf9014d4ba","arxiv_id":"2411.17632","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"A super-Eddington accretion flow around the black hole in TXS 0506+056 can reproduce both its 2014-2015 neutrino flare and its ten-year steady neutrino emission.","lead":"The paper tests whether neutrinos from the known source TXS 0506+056 could be made in the hot gas falling onto its black hole, not in the jet that points toward Earth. If right, it would link this blazar's neutrinos to the same accretion-powered mechanism proposed for other galaxies, and predict new signals for future telescopes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model underproduces the adopted flare luminosity at its stated mdot=10: Eqs. (3)-(4) require mdot ~20-50, so the paper's central 'mdot~10' is not self-consistent.","rationale":"The paper is a clearly written feasibility study with honest caveats: it parameterizes the acceleration timescale and the cosmic-ray loading, varies eta over two decades, and acknowledges that eta < 30000 is needed. I do not find a fatal flaw in the qualitative claim that accretion-flow proton-proton interactions could produce the observed neutrinos, and the gamma-ray absorption argument for hidden cascades is a plausible route around GeV-TeV constraints. The load-bearing weakness I identify is quantitative and internal: the normalization algebra in Eqs. (3)-(4) demands a significantly higher accretion rate than the value stated in the Summary and used in the figures. This is not a disagreement with external consensus; it is a self-consistency issue in the paper's central number. The reader's weakest assumption about eps_CR and eta is related but secondary, because eps_CR enters the same normalization and eta mainly controls the cutoff; a direct luminosity integration would settle whether the factor-of-two discrepancy is absorbed by observational uncertainties or is a genuine underproduction. Since the scenario remains qualitatively viable at somewhat higher mdot, I keep the CONDITIONAL verdict rather than moving to REJECT, but the paper should correct or justify the 'mdot~10' statement and report the required accretion rate as a derived range rather than a fixed input.","tokens_in":15134,"tokens_out":17842,"duration_ms":168682,"concrete_test":"Recompute the model at mdot=10 using Eqs. (3)-(4) and the f_pp(E) curves shown in Figs. 1-4, integrating the resulting per-flavor neutrino spectrum over 32 TeV-4 PeV. If the integrated luminosity is below 1e47 erg/s by more than a factor of 2, repeat the calculation with mdot=22, and as a sensitivity check with eps_CR=0.03 and 0.3, reporting the required mdot for each scenario. This directly settles whether 'mdot~10' is a derived result or an input chosen so the curves pass near the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing problem is internal arithmetic. Eq. (3) gives L_nu_mu = (1/8) f_pp L_CR = 5e45 f_pp (mdot/1) erg/s for MBH=3e8 Msun and eps_CR=0.1. Eq. (4) therefore requires mdot ~ 20 f_pp^-1 to reach the quoted 2014-2015 flare luminosity of 1e47 erg/s. With f_pp ~ 0.9 at low energies and smaller at the multi-PeV proton energies that produce the highest-energy observed neutrinos, the implied mdot is roughly 20-50, not 10. Yet the figures adopt mdot=10 and the Summary states Mdot ~ 10 Mdot_Edd. At mdot=10, the predicted per-flavor neutrino luminosity is ~4.5e46 erg/s, a factor ~2.2 below the adopted 1e47 erg/s. The agreement in Figs. 1-4 therefore relies on the large IceCube error bars, while the quoted accretion-rate requirement is not actually derived from the stated luminosity. This matters because the viability of the scenario hinges on how extreme the super-Eddington accretion must be: if eps_CR is moderately below 0.1 or the effective f_pp in the IceCube band is ~0.3, the required mdot rises to ~50-150, a qualitatively different and much less plausible regime.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the neutrino flare of TXS 0506+056 in 2014-2015 and its time-integrated neutrino emission could originate from the accretion flow around the supermassive black hole, rather than from the relativistic jet. The authors parameterize cosmic-ray acceleration and cooling in magnetically arrested (MAD) and standard/normal (SANE) accretion regimes, compute the steady-state proton spectrum, and derive the pp/pγ neutrino spectra. They report that a super-Eddington accretion rate of approximately 10 times Eddington can explain the flare, while a sub-Eddington rate can sustain the steady emission, and they present the resulting spectra against IceCube data.","tokens_in":15461,"tokens_out":6347,"duration_ms":57722,"significance":"If the model is correct, it offers an alternative to jet-dominated neutrino production in blazars and connects the TXS 0506+056 signal to the accretion-flow neutrino production proposed for Seyfert galaxies such as NGC 1068. The paper is methodically useful because it provides explicit scaling relations for the required accretion rate, compares MAD and SANE regimes, and makes falsifiable predictions for hidden MeV gamma-ray emission and for neutrino emission from misaligned AGNs. The main caveat is that the quantitative claim about the required accretion rate is undermined by an internal arithmetic inconsistency, and several key parameters are assumed rather than derived.","major_comments":[{"comment":"The paper's central claim that a super-Eddington accretion with mdot ~ 10 mdot_Edd explains the 2014-2015 flare is inconsistent with its own Eq. (4). With f_pp ~ 0.9, as shown in Fig. 1b, Eq. (4) gives mdot ~ 22 to reach L_nu_mu = 1e47 erg/s; for f_pp ~ 0.3-0.5 (the values relevant for the multi-PeV protons that produce the highest-energy neutrinos) the required rate is mdot ~ 40-70. At the adopted mdot = 10, Eq. (3) gives L_nu_mu ~ 4.5e46 erg/s for f_pp = 0.9, which is a factor of roughly 2.2 below the quoted flare luminosity. Because the viability of the scenario depends on how extreme the super-Eddington accretion must be, the figures and summary should adopt the higher mdot implied by Eq. (4) and recompute the spectra, or the target luminosity should be explicitly revised downward with a discussion of the consequences.","section":null},{"comment":"The model's ability to explain the observed neutrino flux and cutoff rests on two parameterized, unverified assumptions: the proton acceleration efficiency eta = 300 in Eq. (5) and the cosmic-ray loading fraction eps_CR = 0.1 in Eq. (2). Equation (4) shows that the required accretion rate scales as eps_CR^{-1}, so a modest reduction to eps_CR = 0.03 would push mdot to roughly 60-150 in the flare state. The paper varies eta (Figs. 5-6) and finds that eta < 30000 is needed, but it does not vary eps_CR or justify the fiducial value with simulation-based estimates of the CR loading fraction in MAD/SANE flows. A sensitivity discussion for eps_CR, or at least a citation to relevant particle-in-cell or MHD simulation results, is needed to support the central parameter claim.","section":null},{"comment":"The agreement between the model spectra and the IceCube flux points is partly by construction: the accretion rate is normalized to the neutrino luminosity (Eq. 4) and the dissipation radius R is varied to reproduce the spectral cutoff (e.g., R = 30 vs. 60 R_g in Fig. 1c). This means the fit does not independently test the model's flux normalization or cutoff energy. The non-trivial predictions are the hard, pp-dominated spectral shape, the relative pp vs. pγ contribution, and the hidden MeV cascade. The paper should state this distinction explicitly and frame the results as constraints on mdot and R rather than as an independent reproduction of the observed flux. A concrete test of the model would be to fix mdot from an independent accretion-rate estimate (e.g., from the X-ray Eddington ratio) and check whether the predicted neutrino flux then matches IceCube.","section":null}],"minor_comments":[{"comment":"The phrase 'no comic rays' appears in both the abstract and the introduction; it should read 'no cosmic rays'.","section":null},{"comment":"The text says 'protron' instead of 'proton' in the description of the injection spectrum.","section":null},{"comment":"The reference to Zathul et al. (2024) gives the arXiv number as 'arXiv:24.14598'; the correct identifier is arXiv:2411.14598, and the author list and title should be checked.","section":null},{"comment":"The caption lists 'eta = 300, beta = 10 and mdot = 10' but also repeats 'beta = 10' in the sentence 'We adopt the parameters ... beta = 10, eta = 300, beta = 10 and mdot = 10.' The duplication should be removed.","section":null}],"recommendation":"major_revision","confidential_remarks":"The arithmetic inconsistency between Eq. (4) and the adopted mdot = 10 is straightforward to fix but changes the central quantitative parameter of the paper. The qualitative conclusion that a super-Eddington accretion flow is required would likely survive revision, so I view this as a major-revision issue rather than a rejection. The authors should also address the sensitivity to eps_CR and clearly distinguish construction from prediction in their fits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What should you know? This paper is a serious, clearly written feasibility study of an accretion-flow origin for TXS 0506+056's neutrinos, with a quantitative two-regime (MAD and SANE) treatment. The central claim is not that the jet is ruled out, but that a super-Eddington accretion flow can produce the 2014-2015 flare and a sub-Eddington flow can explain the steady excess. That is worth taking seriously.\n\nWhat's genuinely new: the specific pp-dominated super-Eddington MAD scenario with the resulting hard neutrino spectrum, the demonstration that the dissipation radius must be outside ~30 Rg in MAD to avoid proton synchrotron losses, and the predictions for hidden MeV cascades and neutrino emission from radio galaxies. The authors also cite Zathul et al. (2024) for the same broad idea and are clear about what their own model adds.\n\nThe soft spots. The most concrete problem is an internal arithmetic inconsistency. Eq. (4) correctly states that mdot ~ 20 f_pp^{-1} is needed to reach the adopted flare luminosity of 1e47 erg/s. Yet the figures and the summary use mdot=10, and the summary says 'mdot ~ 10' is needed. With f_pp ~0.9 at low energies and lower at the multi-PeV proton energies that dominate the highest-energy neutrinos, the required mdot is more like 20-50, not 10. The plotted spectra still sit within the large IceCube error bars, so the model is not dead, but the quoted accretion-rate requirement is understated by about a factor of two to five. That should be reconciled in revision.\n\nBeyond that, the acceleration efficiency (eta=300, eps_CR=0.1) is parameterized, not derived; the paper acknowledges this and explores eta, but a skeptic will want a physical argument that reconnection or turbulence in a MAD actually delivers that. The dissipation radius is chosen partly to match the spectral cutoff, so the spectral agreement is partly by construction; this is common in the field but worth stating more explicitly. No code or data files are provided, which limits reproducibility of the spectra.\n\nWho is this for? People working on AGN neutrino production and the NGC 1068-TXS connection. It deserves a serious referee; it's a testable model with concrete predictions (MeV cascade, radio-galaxy neutrinos) and the internal inconsistency is fixable. I'd send it to review with the request that the authors address the mdot discrepancy explicitly.","headline":"A serious feasibility study with a real internal arithmetic inconsistency that understates the required accretion rate; still worth refereeing.","tokens_in":16039,"tokens_out":3679,"would_cite":true,"duration_ms":30047,"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":"This paper argues that the neutrino flare and steady emission from TXS 0506+056 can be produced in the accretion flow around its supermassive black hole, requiring a super-Eddington accretion rate during the 2014-2015 outburst.","keywords":["neutrino astronomy","TXS 0506+056","accretion flow","supermassive black hole","IceCube","proton-proton collisions","magnetically arrested accretion","blazar"],"falsifier":"Measure the bolometric luminosity or Eddington ratio of TXS 0506+056 during the 2014-2015 neutrino flare: if direct accretion diagnostics show an Eddington ratio well below ten, the required super-Eddington rate would be excluded. Alternatively, detect the predicted MeV cascade bump accompanying the neutrino flare with a future MeV telescope; its absence would contradict the hidden-source expectation.","tokens_in":14911,"feed_emoji":"🔭","tokens_out":9153,"duration_ms":76971,"temperature":0.7,"pith_summary":"TXS 0506+056 is the first blazar associated with high-energy neutrinos, and those neutrinos are usually attributed to its relativistic jet. This paper argues they could instead come from the accretion flow around the central supermassive black hole, which matters because AGN jets may be Poynting-flux-dominated and carry few cosmic rays. The authors show that a super-Eddington accretion rate of about ten times Eddington during 2014-2015 can power the observed neutrino flare through proton-proton collisions, while a sub-Eddington flow can account for the decade-long steady neutrino flux. Because the emitting region is compact, the gamma rays produced alongside neutrinos would be absorbed and re-emerge at MeV energies, making TXS 0506+056 a hidden gamma-ray source that future MeV telescopes can test.","feed_headline":"Neutrinos may come from the black hole's accretion flow, not its jet","feed_subtitle":"A super-Eddington disk around TXS 0506+056 can explain IceCube's 2014-2015 flare and steady flux.","key_machinery":"The machinery is a one-zone accretion-flow model built on a parametrized proton acceleration timescale $t_{\\rm acc} = \\eta\\, r_L/c$ with fiducial $\\eta = 300$, and a cosmic-ray luminosity $L_{\\rm CR} = \\epsilon_{\\rm CR} \\dot{M} c^2$ with $\\epsilon_{\\rm CR}=0.1$. The proton population is evolved through a transport equation that balances injection, cooling (proton-proton, photohadronic, Bethe-Heitler, and proton synchrotron losses), and escape (diffusion plus infall), and the resulting neutrino spectra are computed for two accretion regimes: MAD, a highly magnetized flow with slow radial infall and strong magnetic fields, and SANE, a lower-field regime. The argument proceeds by comparing timescales: in MAD, proton synchrotron cooling is so rapid that the dissipation radius must be pushed to roughly $30$-$60\\,R_g$, while in SANE the weaker field allows smaller radii and makes photohadronic production competitive at high energies.","core_discovery":"The paper's central claim is that the neutrino emission from TXS 0506+056 need not originate in the relativistic jet at all. Modelling the inner region as a magnetically arrested accretion flow (MAD) or a standard and normal evolution flow (SANE), with protons accelerated by magnetic reconnection or turbulence and with ten percent of the accretion power going into cosmic rays, the authors find that an accretion rate $\\dot{M} \\sim 10\\,\\dot{M}_{\\rm Edd}$ reproduces the 2014-2015 IceCube neutrino outburst with per-flavor luminosity $\\sim 10^{47}\\,{\\rm erg\\,s^{-1}}$, while $\\dot{M} \\sim 0.1\\,\\dot{M}_{\\rm Edd}$ reproduces the ten-year time-integrated flux. In the flare state, proton-proton collisions dominate and give a hard neutrino spectrum with a cutoff set by pion cooling; in the sub-Eddington steady state, photohadronic interactions can dominate, especially in the SANE regime. In the MAD case the strong magnetic field forces the dissipation site out to tens of gravitational radii so that proton and pion synchrotron cooling does not suppress the neutrino flux.","pith_inferences":["A natural extension is to apply the same accretion-flow picture to other IceCube blazar candidates; if neutrino flares trace super-Eddington accretion episodes, they should correlate with X-ray and optical disk-state changes rather than with GeV gamma-ray flares.","A multi-zone treatment that follows the radial dependence of density and magnetic field inside a MAD could sharpen the predicted cutoff energy and the MeV cascade luminosity, making the model easier to confirm or rule out.","If the accretion-flow origin is correct, the neutrino luminosity per unit Eddington luminosity may be similar across radio-loud and radio-quiet AGN, consistent with the reported X-ray-neutrino correlation; stacking neutrino data on AGN with measured Eddington ratios could test this.","The pp-dominated flare spectrum predicts a relatively flat neutrino energy flux up to a cutoff around the pion-cooling energy; comparing the detailed spectral shape with the jet-based photohadronic models should be possible with the next generation of neutrino telescopes."],"forward_implications":["If the accretion flow, not the jet, powers the neutrinos of TXS 0506+056, then a Poynting-flux-dominated or pair-dominated jet remains consistent with IceCube's detection.","The model predicts a hard neutrino spectrum from proton-proton collisions during the 2014-2015 flare, distinct from the softer photohadronic spectrum usually invoked for sources like NGC 1068.","Gamma rays produced alongside the neutrinos are absorbed by the dense disk and corona photon field, so the hadronic core should appear as a hidden MeV source rather than a bright GeV-TeV emitter.","Misaligned radio galaxies, whose jets do not point at Earth, would also produce neutrinos through the accretion flow if this mechanism is generic, enlarging the population of detectable neutrino sources.","The acceleration efficiency cannot be arbitrarily low: with $\\eta = 30000$ the neutrino cutoff softens below the observed spectrum, so the model requires $\\eta < 30000$ and works with the fiducial $\\eta = 300$."],"supporting_citations":[{"why":"Provides the 2014-2015 neutrino flare detection with about 13 high-energy events and the per-flavor luminosity the flare model must reproduce.","marker":"IceCube Collaboration et al. 2018b"},{"why":"Provides the ten-year time-integrated neutrino flux that the steady-state sub-Eddington model must reproduce.","marker":"IceCube Collaboration et al. 2022"},{"why":"Provides the original high-energy neutrino event coincident with flaring TXS 0506+056, establishing the source as a neutrino emitter.","marker":"IceCube Collaboration et al. 2018a"},{"why":"Supplies the average hard X-ray luminosity used to normalize the coronal photon field and motivates a common neutrino production mechanism across AGN types.","marker":"Kun et al. 2024"},{"why":"Supplies the gamma-gamma absorption treatment for compact coronae and disks and the effective Bethe-Heitler cross section used in the cooling calculation.","marker":"Murase et al. 2020"},{"why":"Defines the magnetically arrested accretion regime and the slow infall velocity used to set the gas density and magnetic field in the MAD scenario.","marker":"Narayan et al. 2003"},{"why":"Supplies the magnetic-field estimate for super-Eddington MAD flows and the precedent of neutrino production in that regime.","marker":"Hayasaki & Yamazaki 2019"},{"why":"Provides the proton-proton cross section and inelasticity used to compute hadronic interaction efficiencies and neutrino yields.","marker":"Kelner et al. 2006"}],"fun_headline_variants":["Accretion flow, not jet, may power TXS 0506+056 neutrinos","Accretion disk, not jet, may explain blazar's neutrino flare","Black hole's feeding disk, not jet, may produce neutrinos","Neutrino emission may trace black hole accretion, not jet","Super-Eddington disk may feed TXS 0506+056 neutrino flare"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model holds only if protons in the inner accretion flow can actually be accelerated to tens-to-hundreds of PeV with about ten percent of the accretion power going into cosmic rays; if acceleration is much slower or less efficient, the required accretion rate would become unrealistically high.","fun_headline_variants_meta":{"raw":{"variants":["Accretion flow, not jet, may power TXS 0506+056 neutrinos","Accretion disk, not jet, may explain blazar's neutrino flare","Black hole's feeding disk, not jet, may produce neutrinos","Neutrino emission may trace black hole accretion, not jet","Super-Eddington disk may feed TXS 0506+056 neutrino flare"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000974,"raw_usage":{"total_tokens":4187,"prompt_tokens":1040,"completion_tokens":3147,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":3046}},"tokens_in":656,"tokens_out":3147,"duration_ms":59826,"temperature":1.0,"reasoning_tokens":3046,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:54:09.089230+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the bolometric luminosity or Eddington ratio of TXS 0506+056 during the 2014-2015 neutrino flare: if direct accretion diagnostics show an Eddington ratio well below ten, the required super-Eddington rate would be excluded. Alternatively, detect the predicted MeV cascade bump accompanying the neutrino flare with a future MeV telescope; its absence would contradict the hidden-source expectation.","supporting_citations":[],"review_version":1}