{"id":"1f72c1de-a42f-449b-b5f4-8c1a3d762935","arxiv_id":"2603.27749","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A new hybrid simulation code models proton acceleration in turbulent black hole coronae and reproduces the IceCube neutrino signal of NGC 1068 and other Seyferts with standard coronal parameters.","lead":"This paper introduces Turb-AM3, a code coupling stochastic proton acceleration in turbulent plasma with full radiative transfer, and applies it to model neutrinos from the black hole corona in NGC 1068. The model reproduces the IceCube neutrino signal with standard coronal parameters while keeping gamma-ray emission hidden, supporting the corona as the neutrino source.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"IceCube reproduction depends on the imposed electron/X-ray field, not on solved coronal physics","rationale":"The reader's weakest-assumption statement explicitly identifies the imposed Maxwellian electron distribution and the resulting X-ray target field as a fragile premise. My stress-test pass converges on that same point as the single most load-bearing concern: the dominant neutrino channel depends on a target photon field that is an input assumption, not a solved output of the self-consistent calculation. This is not a claim of internal inconsistency; the timescale separation argument in §2.4 motivates holding the electron distribution fixed, and the paper is transparent about deferring detailed electron modelling. It is, however, a direct qualification of the central 'self-consistent reproduction' claim, because a modest change in the assumed electron temperature or spectral shape can shift the pγ neutrino flux and peak energy. The proposed test isolates this dependence without requiring a full rerun of the coupled turbulence calculation. Since the reader already recommended CONDITIONAL for related reasons, my assessment does not move the verdict; it reinforces the condition. I also credit the paper for its honest presentation of parameter degeneracies and the clear statement in the Conclusions that electron energy-loss/energization processes are left to future work.","tokens_in":26944,"tokens_out":6997,"duration_ms":85736,"concrete_test":"Take the fiducial proton spectrum of Fig. 2 and recompute the neutrino spectrum with AM3 while replacing the internally produced X-ray field by: (a) the same thermal Comptonization spectrum with Θ_e = 0.1 and 0.4; and (b) the observed 3–79 keV NuSTAR spectrum of NGC 1068, normalized to the same 2–10 keV luminosity L_X = 6.8×10^43 erg/s, keeping all other parameters fixed. If E^2Φ_ν at 10 TeV, or the 1–30 TeV spectral slope, changes by more than ~30–50%, the IceCube reproduction is controlled by the imposed electron/X-ray field rather than by the turbulent acceleration physics. A complementary run with a few-percent nonthermal electron tail would test sensitivity to that as well.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that a standard set of coronal parameters reproduces the IceCube neutrino flux. But the dominant 1–30 TeV neutrino emission is produced through pγ interactions (Fig. 2, §3.1), and the pγ rate is proportional to the coronal X-ray photon density. That photon density is not a solved quantity: it is generated in AM3 by inverse-Compton scattering of disk photons off an electron distribution that is imposed as a fixed Maxwellian with Θ_e = 0.2 (§2.1, §2.4, Table 1). The paper explicitly postpones detailed electron energization and notes that the electron distribution around the thermal peak is not finely resolved (§3.1, Conclusions). Thus the neutrino normalization and spectral peak are conditional on an assumed electron temperature, on the absence of a nonthermal electron tail, and on the specific Comptonization implementation. If the real corona has a harder X-ray tail, a different cutoff energy, or a nonthermal electron component, the pγ optical depth and the proton energy at which losses balance acceleration both change, moving both the neutrino peak and the flux normalization. The claim that a 'standard parameter set' reproduces IceCube is therefore load-bearing on an input assumption rather than on a self-consistently derived target radiation field. This is an internal limitation acknowledged by the authors, but it directly qualifies the headline claim of self-consistency for the neutrino prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces Turb-AM3, a hybrid code coupling the AM3 lepto-hadronic radiative solver with a stochastic proton acceleration module that includes turbulent damping by accelerated particles. After presenting the code and physical setup in §2, the authors apply it to a one-zone model of the NGC 1068 corona in §3.1. They claim that a standard parameter set (R_cor ≈ 15 rg, v_A ≈ 0.25c, B ≈ 10^3 G, L_d ≈ 5e44 erg/s, ξ_p = 0.1) reproduces the IceCube neutrino flux while remaining compatible with gamma-ray upper limits, and that for ξ_p ≳ 10^-2 the proton–turbulence system self-regulates, producing a universal neutrino spectral template. The paper also extends the model to NGC 4151 and NGC 7469 by varying only the disk luminosity, and explores sensitivities to ξ_p, t_esc/t_acc, and t_adv/t_acc.","tokens_in":27214,"tokens_out":3147,"duration_ms":39980,"significance":"If the central claim holds, the paper would provide a concrete physical mechanism—turbulent acceleration with self-consistent back-reaction—that explains both the normalization and shape of the IceCube neutrino signal from AGN coronae, while keeping the gamma-ray flux suppressed via electromagnetic cascades. The prediction of a universal neutrino template and the explicit energy-budget consistency checks are useful and falsifiable. However, the significance is conditional: the 'self-consistent' result explicitly relies on an imposed, fixed Maxwellian electron distribution and on an X-ray photon field that is therefore an input assumption rather than a solved outcome. The paper also does not ship code or data, and it does not provide convergence tests, which weakens the reproducibility of the numerical claims. These issues are correctable within the paper's scope, so I regard the work as promising but not yet ready for acceptance as is.","major_comments":[{"comment":"The neutrino flux normalization and peak are driven by pγ interactions whose target photons are the coronal X-ray field. That field is computed from an assumed, fixed Maxwellian electron distribution at Θ_e = 0.2, not solved self-consistently. The paper itself states in §3.1 that 'the electron distribution around the thermal peak is not finely resolved' and postpones electron energization to future work. The headline 'self-consistent' reproduction of the IceCube flux is therefore conditional on an input assumption. The authors should quantify the sensitivity: e.g., vary Θ_e around 0.2 or add a modest nonthermal electron tail and recompute the neutrino spectrum, to show how much the normalization and peak move.","section":"§2.1, §2.4, §3.1"},{"comment":"The comparison procedure in Fig. 4 tunes v_A to a fixed pivot point (Eν = 10 TeV, E²Φν = 1.4e-11 erg/cm²/s). Since v_A controls both the acceleration rate and, through Eq. (3), B and the coronal magnetic energy, this tuning strongly influences the predicted neutrino flux. The fiducial set in §3.1 is not tested against a quantitative goodness-of-fit statistic; the claim of 'satisfactorily reproduced' rests on visual overlap. The authors should present untuned predictions and a statistical measure (e.g., chi-square or likelihood) for the fiducial and tuned cases, given the degeneracies displayed in Figs. 4–6.","section":"§3.2, Fig. 4"},{"comment":"No convergence tests are presented for the numerical method: no resolution study in momentum space, no timestep convergence, no comparison of the coupled acceleration–cascade solver against an analytic limit. The universal spectral shape and saturation behavior claimed in §3.2–3.3 could be numerical artifacts if the grid does not resolve the acceleration–loss balance near the TeV cutoff. The authors should add convergence tests and ideally release the code or a reproducible input/output package, as this is central to validating the new Turb-AM3 framework.","section":"§2.7, §3.1, §3.3"},{"comment":"The 'universal' spectral shape is asserted on the basis of visual inspection of Figs. 4, 5, and A.1. No quantitative definition or metric is given for convergence to this template (e.g., a fractional deviation of the peak energy or spectral slope across ξ_p and t_esc/t_acc). Since the universality is a central new result, the authors should define and compute a shape-difference measure for the neutrino spectrum in the self-regulated regime.","section":"§3.2, §3.3"}],"minor_comments":[{"comment":"Typographical errors: 'excape' in §2.1, 'respectivly' in the Introduction, and 'articially' in §3.2. The manuscript would benefit from a careful proofread.","section":"Abstract, passim"},{"comment":"The figure contains many curves and labels; the color/line assignments are complex. A table of components (or a simplified version separating photon and neutrino panels) would improve readability.","section":"Fig. 2"},{"comment":"The footnote about Eq. (10) is helpful, but the limiting case t_adv ≫ t_esc is stated only verbally; it would be clearer to write the explicit limiting form.","section":"§2.5, Eq. (10)"},{"comment":"The kernel φ(k,p) is said to be 'normalized to unity' but the precise normalization variable in the Gaussian is ambiguous (ln k vs k). Clarify the argument of the Gaussian and the normalization measure.","section":"§2.6, Eq. (11)"},{"comment":"Several 2026 references are cited as 'arXiv e-prints' with no journal page numbers; if any have since been published, the bibliographic entries should be updated.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid code paper, but the central 'self-consistent reproduction' claim is weakened by the imposed electron distribution and target photon field, and the numerical verification standard is below what I would expect for a new code paper. The concerns are fixable with additional calculations and a more cautious framing; this is not a rejection. I would encourage the editor to ask for the code/data or at least a detailed convergence study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: the new thing here is Turb-AM3, a coupling of the AM3 lepto-hadronic solver to a stochastic acceleration module that includes the back-reaction of accelerated protons on the turbulent cascade. That self-regulation is physical and is the best part of the paper. For ξ_p above about 1e-2, the proton population saturates near the turbulent energy density and the neutrino spectrum converges to a template shape that is insensitive to injection details. The Appendix showing that Fokker-Planck, generalized Fermi, and advective prescriptions give similar spectra strengthens that claim. This is a real step beyond the semi-analytic and steady-state treatments in the prior literature.\n\nThe soft spot is the X-ray field. The 1–30 TeV neutrinos are pγ-dominated, and the target photon density is not solved for: it is produced by inverse Compton scattering of disk photons off an imposed Maxwellian electron population at Θ_e ≈ 0.2, held fixed and admittedly not finely resolved at the thermal peak (Sec. 3.1). So the neutrino normalization and peak position are conditional on an input assumption. The paper is honest about this, but it means \"self-consistent\" applies to the proton–turbulence loop, not to the full radiation field. If the real corona has a harder X-ray tail or a nonthermal electron component, the pγ optical depth changes and the predicted flux shifts. The match to IceCube is therefore a demonstration that the framework can accommodate the data within a plausible parameter set, not a unique prediction of the model. That is still worth having.\n\nSecond soft spot: no code, no data, no convergence tests, and the quality of the IceCube match is judged visually rather than with a fit statistic. For a paper whose main deliverable is a new numerical tool, that is a real gap. A referee should ask for a public version with parameter files, at least a grid-convergence check, and a quantitative comparison against the IceCube model fit. Also note that the parameter exploration tunes v_A to hit a pivot point (Sec. 3.2); that is fine for extracting templates, but it weakens the \"standard parameters reproduce the flux\" phrasing.\n\nWho is this for? Anyone modeling neutrino production in compact, radiatively dense environments—coronae, hidden sources, and possibly jets. The template spectra are a concrete prediction for next-generation telescopes, and the paper maps which parameters control the spectral shape. I would send it to peer review. Not because the central fit is watertight—it is not—but because the framework is new, the self-regulation mechanism is well motivated, and the template predictions deserve scrutiny. A good referee will push for the code release and a proper treatment of the electron population.","headline":"A genuinely new modeling framework with a physical self-regulation loop, but the IceCube 'reproduction' leans on an imposed X-ray field and no code is out yet — worth refereeing, not worth taking as a hard prediction.","tokens_in":27777,"tokens_out":2434,"would_cite":true,"duration_ms":28225,"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":"Self-consistent turbulent corona model reproduces the IceCube neutrino flux of NGC 1068","keywords":["active galactic nuclei","NGC 1068","neutrino astrophysics","stochastic particle acceleration","turbulent corona","IceCube","multi-messenger astrophysics","turbulence damping"],"falsifier":"A firm IceCube detection of neutrinos from NGC 1068 above ~100 TeV at a level exceeding the model's sharp cutoff, or a precise measurement of the <1 TeV neutrino flux that does not match the predicted flat pp component, would show that the mechanism or the assumed coronal parameters are wrong.","tokens_in":26711,"feed_emoji":"🧊","tokens_out":4033,"duration_ms":37485,"temperature":0.7,"pith_summary":"This paper tries to show that the TeV neutrinos IceCube sees from NGC 1068 can be produced by protons accelerated in the turbulent corona around the supermassive black hole, without invoking exotic parameters or violating gamma-ray upper limits. The authors claim that when the back-reaction of accelerated protons on the turbulent cascade is included, the proton population self-regulates: energetic protons damp the turbulence that accelerates them, capping their own energy density near the turbulent magnetic pressure. Under these conditions the neutrino spectrum takes a universal template shape, and the observed flux is reproduced for standard coronal parameters. The same template, scaled only by X-ray luminosity, also accounts for the IceCube signals from NGC 4151 and NGC 7469. The significance is that neutrino observations become a direct probe of coronal plasma conditions.","feed_headline":"Turbulence feedback explains NGC 1068's IceCube neutrinos","feed_subtitle":"Self-consistent proton acceleration in the black hole corona yields the observed TeV neutrinos without exceeding gamma-ray limits.","key_machinery":"Turb-AM3, a hybrid code coupling the time-dependent lepto-hadronic radiative solver AM3 with a stochastic acceleration module. The proton transport equation uses either a Fokker-Planck diffusion operator or a generalized Fermi (master) operator, and is solved together with a turbulent cascade equation that includes a damping kernel phi(k,p) transferring turbulent energy from wavenumber k to protons of momentum p. This coupling makes the acceleration rate time-dependent: as protons extract energy from the cascade, the turbulence is quenched at small scales, which automatically limits the non-thermal proton energy density.","core_discovery":"The central claim is that stochastic proton acceleration in the magnetized corona of NGC 1068, treated self-consistently with the damping of turbulence by the accelerated protons, reproduces the IceCube neutrino flux for R_cor ~ 15 r_g, v_A = 0.25c, B ~ 1e3 G, and L_d = 5e44 erg/s, while remaining consistent with gamma-ray upper limits. When the non-thermal proton fraction xi_p is at least about 1e-2, the system self-regulates and the proton spectrum converges to a near-universal shape with equal energy per decade, cut off at a few tens of TeV by photohadronic losses; the resulting neutrino spectrum has a broad peak around 1-10 TeV with a flat low-energy extension from pp interactions. The s","pith_inferences":["If the universal template is correct, the neutrino spectral shape alone cannot pin down the proton injection fraction xi_p; the main observable discriminators are the low-energy pp component and the exact cutoff energy.","The model implies a tight relation between X-ray luminosity and neutrino peak energy: brighter coronae produce more neutrinos but at lower energies; this can be tested with IceCube's growing sample of Seyfert galaxies.","Extended to other compact environments (X-ray binaries, tidal disruption events), the same code would predict neutrino spectra whose shape is set by the same self-regulation mechanism, giving a way to test the generality of the picture.","Because the electron distribution is imposed rather than solved, a fuller test of the scenario requires coupling to a thermal Comptonization model, which the authors flag as future work."],"forward_implications":["If the claim holds, standard corona parameters (R_cor ~ 15 r_g, v_A ~ 0.25c, B ~ 1e3 G, L_d ~ 5e44 erg/s) are sufficient to explain the IceCube neutrino signal from NGC 1068.","For xi_p greater than roughly 1e-2, the neutrino flux normalization becomes insensitive to the exact injection fraction, because turbulent damping self-regulates the proton energy density.","The neutrino spectrum follows a universal template: broad peak near 1-10 TeV, flat low-energy extension from pp interactions, and a sharp cutoff around tens of TeV set by photohadronic losses.","Detections near ~100 TeV (as hinted for NGC 7469) are accommodated by varying X-ray luminosity, while emission above ~100 TeV is suppressed.","The sub-TeV pp component is a diagnostic of proton density and hence of the pair-loading factor in the corona."],"fun_headline_variants":["Turbulence self-regulation yields NGC 1068's IceCube neutrinos","Self-consistent corona model fits IceCube neutrino signal","Proton acceleration in corona reproduces neutrino flux","Corona turbulence damping shapes neutrino spectrum","How black hole corona turbulence yields TeV neutrinos"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The corona is treated as a single, uniform, one-zone region with a constant size, constant advection speed, and a fixed thermal electron population; if the real corona is strongly stratified or has a non-thermal electron component, the target photon field and hence the p-gamma neutrino spectrum could change substantially.","fun_headline_variants_meta":{"raw":{"variants":["Turbulence self-regulation yields NGC 1068's IceCube neutrinos","Self-consistent corona model fits IceCube neutrino signal","Proton acceleration in corona reproduces neutrino flux","Corona turbulence damping shapes neutrino spectrum","How black hole corona turbulence yields TeV neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1137,"prompt_tokens":836,"completion_tokens":301,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":224}},"tokens_in":580,"tokens_out":301,"duration_ms":3506,"temperature":1.0,"reasoning_tokens":224,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T17:08:56.281578+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A firm IceCube detection of neutrinos from NGC 1068 above ~100 TeV at a level exceeding the model's sharp cutoff, or a precise measurement of the <1 TeV neutrino flux that does not match the predicted flat pp component, would show that the mechanism or the assumed coronal parameters are wrong.","supporting_citations":[],"review_version":1}