{"id":"6eaec1ca-44e6-44cb-be7f-21ec4919ce1b","arxiv_id":"2412.01457","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A self-regulated stochastic acceleration model, where fast protons damp the corona's turbulence, reproduces the 30 to 300 TeV proton spectrum inferred from IceCube's NGC 1068 neutrino observations.","lead":"This paper models how protons in the turbulent, magnetized coronae around supermassive black holes get accelerated and produce the TeV neutrinos IceCube sees from nearby Seyfert galaxies like NGC 1068. It finds that when accelerated protons damp the turbulence that feeds them, the process self-regulates and can explain the observed spectrum without hand-tuning the total energy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unvalidated extrapolation of the generalized Fermi kernel from vA=0.4c to 0.2c and from 10 to 150 lc/c is the weakest load-bearing link for the claimed self-regulated spectral fit.","rationale":"The reader's verdict of CONDITIONAL is appropriate. The paper is transparent: it states the difficult extrapolation, shows parameter sensitivity, and does not overclaim. However, the central 'no ad hoc normalization' claim hinges on the damping model and on the acceleration kernel. Between these, the kernel extrapolation is the least secure because it is explicitly unverified at the operating point, while the damping model's energy-balance normalization u_p/u_B ~ c/v_A is a simple consequence of energy conservation and less likely to be off by orders of magnitude. I therefore agree with the reader's weakest_assumption. A dedicated MHD simulation at v_A = 0.2c would directly settle whether the kernel rescaling is valid. One caveat: the self-regulated scenario is also computed with the Fokker-Planck model (a), which does not use the generalized Fermi kernel, so a failure of the kernel would weaken but not necessarily destroy the conclusion. This is why UNCHANGED (CONDITIONAL) is the right verdict rather than REJECT.","tokens_in":21872,"tokens_out":16515,"duration_ms":138604,"concrete_test":"Run the same MHD turbulence simulation that was used to calibrate the kernel at v_A = 0.4c, but with v_A = 0.2c, keeping box size, resolution, and magnetization parameter the same. Directly measure the jump kernel phi(epsilon|epsilon') and the time-dependent Green function G_turb(epsilon; epsilon', tau) out to tau ~ 150 l_c/c, and compare the measured kernel width and late-time Green function to the prediction obtained by rescaling the v_A = 0.4c kernel width by (0.2/0.4)^2 and integrating Eq. (5). If the width differs by more than the statistical uncertainty (say, 20%) or the late-time Green function deviates from the master-equation result, the self-regulated spectrum in Fig. 3 is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central new result is the self-regulated scenario of Sect. 3 (bottom panel of Fig. 3), in which the proton flux normalization is set by balancing particle energy extraction against the turbulent cascade injection rate rather than by a free normalization constant. For that result to hold, the acceleration rate and the shape of the energy gain probability must be correct at the coronal parameters used (v_A ~ 0.2c, l_c = 10 r_g, tau_adv ~ 150 l_c/c). The generalized Fermi kernel phi(epsilon|epsilon') appearing in Eq. (5) was, however, measured in one MHD simulation at v_A = 0.4c over about 10 l_c/c, and the paper rescales its width as (v_A/0.4c)^2 and applies it to v_A = 0.2c and to integration times ~15 times longer, explicitly acknowledging in Sect. 2.2 that this extrapolation is nontrivial. The same kernel underlies model (b) in the self-regulated calculation, where it sets nu_acc, the cutoff location, and the balance condition u_p/u_B ~ c/v_A that fixes the normalization. If the kernel shape or its v_A scaling changes with v_A or with time, the predicted 30-300 TeV slope and cutoff move, and the claimed fit and even the self-regulated flux level are not robust. No benchmark, simulation, or analytic argument is provided to support the extrapolation at coronal conditions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Taking the IceCube detection of TeV neutrinos from NGC 1068 at face value, the paper asks whether stochastic acceleration in a magnetized, turbulent accretion-disk corona can produce the inferred parent proton spectrum (roughly E^{-3} between 30 and 300 TeV, with a pressure ratio of order 0.1). It reformulates the transport problem using a generalized Fermi kernel in addition to the standard Fokker-Planck operator, models advective escape as a hard integration-time limit, and includes turbulent transport in the escape time. It then compares three scenarios: (i) test-particle turbulent acceleration, (ii) turbulent pre-acceleration followed by jet-shear acceleration, and (iii) turbulent acceleration self-regulated by the damping of turbulence by the accelerated protons. The authors conclude that all three can reproduce the inferred spectrum for reasonable parameters, but that only the self-regulated scenario fixes the flux normalization without ad hoc choices, since backreaction forces the suprathermal energy density to saturate near u_p/u_B ~ c/v_A.","tokens_in":22219,"tokens_out":7918,"duration_ms":61359,"significance":"The paper's strength is that it treats the transport equation more carefully than most phenomenological studies: advection is imposed as an acceleration-time limit rather than a leakage term, and the comparison between a Fokker-Planck and a generalized Fermi kernel makes the theoretical uncertainty visible. The authors are also explicit that scenarios (i) and (ii) require an ad hoc proton energy content. If the self-regulated scenario is robust, it is significant because it converts the otherwise surprising neutrino luminosity of NGC 1068 into a consequence of turbulence damping and yields a concrete prediction for how the high-energy cutoff depends on v_A, l_c, and tau_adv. However, the claim is not yet backed by a quantitative fit and rests on an extrapolation of the acceleration kernel that the paper itself flags as non-trivial, so the significance is conditional.","major_comments":[{"comment":"The generalized Fermi kernel phi(epsilon|epsilon') entering Eq. (5) was calibrated in an MHD simulation at v_A = 0.4c over about 10 l_c/c, and is applied here at v_A approximately 0.2c with its width rescaled as (v_A/0.4c)^2, integrated up to tau_adv approximately 150 l_c/c. The authors explicitly acknowledge in Sect. 2.2 that this extrapolation is non-trivial, but it is load-bearing: in the self-regulated scenario of Sect. 3 (bottom panel of Fig. 3) this kernel sets nu_acc, the cutoff energy, and the balance condition u_p/u_B approximately c/v_A that fixes the flux normalization. If the kernel shape or its v_A scaling changes at coronal conditions, the predicted 30-300 TeV slope and cutoff, and hence the claimed fit, change significantly. Please add a robustness test, e.g., vary the kernel width over the uncertainty of the simulation calibration, benchmark the kernel at v_A approximately 0.2c in a PIC/MHD run longer than 10 l_c/c, or provide an analytic argument for the (v_A/0.4c)^2 rescaling.","section":"Section 2.2, Eq. (5)"},{"comment":"The bottom-panel caption of Fig. 3 states beta_p = 0.5, while the text and Appendix A state that beta_p = 1 is used everywhere. This matters because the damping threshold u_p/u_B approximately c/v_A and the quoted energy densities in units of p_gas (3.4 for model (a), 4.1 for model (b)) depend on the magnetic energy density. The inconsistency should be resolved and the reported values recomputed for the intended beta_p. In addition, model (a) uses v_A = 0.21c and model (b) v_A = 0.28c, so the location of the cutoff is partly set by these choices; the claim that the self-regulated scenario reproduces the flux 'without additional fine-tuning' would be strengthened by showing the dependence of the bottom-panel spectrum on v_A and on the injected energy density (one quarter of u_B).","section":"Section 3, Fig. 3 bottom panel; Appendix A"},{"comment":"All three 'satisfactory fits' are established by visual comparison with the butterfly diagram, whose normalization is itself uncertain by a factor of a few (Sect. 2.1). Given the exponential sensitivity of the peak energy to nu_acc tau_adv, the manuscript should provide a quantitative fit statistic over the 30-300 TeV band, or at least a table of trial parameters with the resulting chi-squared or likelihood. It should also specify how the butterfly was represented in the figures and which part of the band drives the comparison.","section":"Section 3 and Fig. 3"},{"comment":"The bottom-panel spectra are obtained with models (a) and (b) 'described by Lemoine et al. (2024)', but the equations for the damped turbulence evolution and for the particle acceleration in the self-regulated regime are not reproduced in this manuscript. Since this is the main new physical ingredient, the reader cannot verify the claimed normalization without going to the earlier paper. Please include the model equations, or a substantial summary in an appendix, so that the calculation is self-contained.","section":"Section 3, self-regulated model"}],"minor_comments":[{"comment":"The first paragraph reads 'recent announcement the arrival directions'; it should be 'recent announcement that the arrival directions'.","section":"Introduction"},{"comment":"The phrase 'we payed particular attention' should be 'we paid particular attention'.","section":"Section 2.2"},{"comment":"The expression 'v_add = alpha_d v_k(r_co)' appears to be a typo for 'v_adv = alpha_d v_k(r_co)'.","section":"Appendix A"},{"comment":"The butterfly diagram's vertical normalization (a pressure ratio of 0.3 corresponding to 0.1, as stated in Sect. 2.1) should be repeated in the caption, since the figure is otherwise hard to interpret.","section":"Fig. 1 caption"},{"comment":"The phrase 'injected particles at gamma_0 about 1 ... with an energy density lower by a factor of four than contained in the magnetized turbulence' should be rephrased as 'a factor of four lower than the magnetic energy density' for clarity.","section":"Section 3"},{"comment":"The reference to Waxman (2024) as 'Talk given at RICAP 24' is not a conventional archival citation; either give a published version or cite it as a private communication.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The reader's stress-test concern about the kernel extrapolation is, in my reading, the correct point to press; the internal beta_p inconsistency in the self-regulated section is an independent problem that is easy to fix. I do not think rejection is warranted, because the transport setup and the self-regulation idea are valuable and the deficits are addressable in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a careful, transparent phenomenological study of stochastic acceleration in Seyfert coronae as the origin of the IceCube NGC 1068 neutrino signal. The genuinely new piece is the third scenario: self-regulated acceleration via turbulence damping, where the proton flux normalization emerges from balancing particle energy extraction against cascade injection, rather than being set by hand. The treatment of advection as a hard time cutoff rather than an escape term is also a real improvement, and the point that turbulent transport dominates escape is worth keeping in mind. The transport equations are set out clearly, and the authors are explicit that the first two scenarios require ad hoc normalization.\n\nThe soft spots are real but mostly acknowledged. The generalized Fermi kernel in Eq. (5) was measured in one MHD simulation at vA=0.4c over about 10 lc/c, then rescaled as (vA/0.4c)^2 and applied at vA~0.2c for integration times up to 150 lc/c. The paper itself flags this extrapolation as nontrivial, but the self-regulated scenario's predicted slope, cutoff, and flux level all sit on top of it. There is no benchmark or independent simulation for the coronal regime, and there is no statistical comparison to the IceCube butterfly, just visual fits. Parameters like vA, lc, vadv, and the cascade rate are chosen by hand; the authors say this, but it means the \"no fine-tuning\" claim is weaker than it sounds. The third scenario still has free parameters; it just removes the normalization constant.\n\nThe paper takes the IceCube-inferred proton spectrum at face value, which is a reasonable choice given the state of the data, and it does not pretend otherwise. The citation pattern looks fair; earlier stochastic acceleration work is credited. The lack of code or numerical details is a minor issue, but it would help a referee to check the damping model.\n\nThis is a paper for people working on AGN neutrino production and particle acceleration. It will get cited for the self-regulation idea and the transport treatment, but the extrapolated kernel means I would treat the central fit as plausible rather than demonstrated. It deserves a serious referee; the physics is sound and the limitations are stated. My recommendation: send it to review, and ask the authors for robustness checks on the kernel extrapolation, or at least a clear statement of how sensitive the self-regulated spectrum is to the kernel shape.","headline":"Honest and useful phenomenology; the self-regulated damping scenario is the genuine new idea, but the extrapolated Fermi kernel makes the central fit plausible rather than proven.","tokens_in":22712,"tokens_out":2057,"would_cite":true,"duration_ms":18820,"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":"A turbulent black hole corona whose accelerated protons damp the turbulence itself can account for the neutrino spectrum IceCube observes from NGC 1068 without fine-tuned proton injection.","keywords":["stochastic particle acceleration","turbulence damping","black hole coronae","Seyfert galaxies","neutrino astrophysics","NGC 1068","Alfvénic turbulence","cosmic ray acceleration"],"falsifier":"Run MHD or particle-in-cell simulations of the energy-jump kernel at an Alfvén velocity near 0.2 of the speed of light for durations up to 150 coherence lengths: if the kernel width does not scale as the assumed velocity-squared law, or if a higher-significance measurement of the NGC 1068 neutrino flux does not track the magnetic energy budget, the central claim is refuted.","tokens_in":21588,"feed_emoji":"🌌","tokens_out":9511,"duration_ms":79241,"temperature":0.7,"pith_summary":"The paper asks whether the proton population that IceCube infers from the neutrino emission of the Seyfert galaxy NGC 1068 can be produced by stochastic acceleration in the magnetically turbulent corona of a supermassive black hole. Its answer is yes, with a caveat: the test-particle versions of turbulent acceleration (with or without a second stage of shear acceleration) reproduce the 30–300 TeV proton spectrum only if the number of injected protons is tuned so that their final energy density approaches the turbulent energy density. The paper's central case is the self-regulated version, in which protons drain enough energy from the turbulence to damp it, so that the acceleration rate stalls and the suprathermal pressure settles naturally near the magnetic pressure. In that version the flux level and shape come out as required without an ad hoc normalization, and the output is controlled by the rate at which turbulent energy is injected on the outer scale (about 10 gravitational radii). This matters because it turns an energy-budget coincidence into a physical prediction and suggests that different Seyfert galaxies will show different neutrino spectra.","feed_headline":"Self-regulated acceleration reproduces the NGC 1068 neutrino spectrum","feed_subtitle":"When accelerated protons damp the turbulence that feeds them, the inferred flux level and shape follow without tuning.","key_machinery":"The load-bearing piece is the generalized Fermi transport kernel, the probability per interaction time that a proton of energy $\\epsilon'$ jumps to energy $\\epsilon$, which was calibrated in magnetohydrodynamic simulations and is assumed to rescale in width with the square of the Alfvén velocity when applied to other conditions. It replaces the purely diffusive Fokker–Planck operator, and the resulting transport equation is integrated with advection implemented as a time limit $\\tau_{\\rm adv} = r_{\\rm co}/v_{\\rm adv}$ on the Green function rather than as an extra escape term. The second key element is the backreaction prescription: once the rate at which particles drain energy from the turbulence exceeds the cascade replenishment rate, the turbulence is damped and the acceleration rate becomes time-dependent, so the system self-regulates; this feedback law, combined with the radiation-loss and escape rates for the fiducial NGC 1068 environment, yields spectra whose slope and normalization emerge from the energy injection rate on the outer coherence scale.","core_discovery":"The paper's central claim is that a self-regulating version of stochastic acceleration is a natural and sufficient mechanism for the 30–300 TeV proton population inferred from IceCube's detection of NGC 1068. In this scenario, accelerated protons extract energy from the turbulent cascade at a rate proportional to the acceleration rate times their energy density; once this exceeds the rate at which the cascade is replenished from the outer scale, the turbulence is damped, the acceleration rate stalls, and the process reaches a stationary state in which the suprathermal proton energy density settles near the magnetic energy density multiplied by the ratio of light speed to the Alfvén velocity. Solving the time-dependent transport equation, with the generalized Fermi kernel as the acceleration law and advection implemented as a strict time cutoff, the authors reproduce the inferred proton spectrum in the 30–300 TeV band with characteristic parameters of an Alfvén velocity near 0.2 of the speed of light, a coherence length of about 10 gravitational radii, a corona of 30 gravitational radii, and an advection speed of 0.02 of the speed of light. In the test-particle and shear-reacceleration alternatives the same fit is possible, but the normalization must be placed by hand at an energy density comparable to that of the turbulence, which the authors argue is itself an indication that the self-regulated regime is the physically relevant one.","pith_inferences":["One can test the extrapolation at the heart of the model by measuring the energy-jump kernel in MHD or PIC simulations at an Alfvén velocity near 0.2 of the speed of light and for runtimes of order 150 coherence lengths; a width that grows more slowly than the assumed velocity-squared scaling would shift the predicted cutoff and break the claimed fit.","The same self-regulation mechanism should apply to other turbulence-dominated cosmic accelerators such as galaxy clusters, radio lobes, and X-ray binary coronae, predicting a near-universal relation between nonthermal particle pressure and magnetic pressure that future multi-messenger observations could test.","The paper's logic implies that the spectral slope seen by IceCube in a single source measures the distribution of acceleration rates across the coronal volume rather than an intrinsic power-law index of the accelerator, so comparing neutrino spectra from several Seyfert galaxies could map coronal inhomogeneity."],"forward_implications":["If the self-regulated scenario is right, the neutrino luminosity of a Seyfert corona is set by the magnetic energy content and the turbulent energy injection rate, not solely by the X-ray luminosity that supplies the photon targets.","Sources with different Alfvén velocities or advection times will show different neutrino spectra, and the paper points to NGC 4151 as a likely example of such variation.","Below the cutoff the self-regulated proton spectra approach a shape close to $d n/d\\epsilon \\propto \\epsilon^{-2}$ per number, so a spectrum steeper than that in the observed band signals that the IceCube window lies near the high-energy cutoff.","The condition for self-regulation, $\\chi_0 \\gtrsim (v_{\\rm A}/c)^{-1}\\beta_p^{-1}(\\epsilon_{\\max}/\\epsilon_{\\rm th})^{-1}$, implies that any source that pushes protons to very high energies should generically operate in the self-regulated regime."],"supporting_citations":[{"why":"It supplies the IceCube detection of NGC 1068 that defines the target neutrino and parent-proton spectrum to reproduce.","marker":"IceCube Collaboration et al. 2022"},{"why":"It provides the inferred 30–300 TeV proton spectrum, its spectral slope, and the pressure normalization used as the benchmark for the fits.","marker":"Murase et al. 2020a"},{"why":"It introduces the generalized Fermi acceleration picture for particles interacting with intermittent magnetic structures.","marker":"Lemoine 2021"},{"why":"It derives the jump kernel used as the acceleration operator in the transport equation.","marker":"Lemoine 2022"},{"why":"It benchmarks the generalized Fermi kernel against particle-in-cell simulations of turbulent particle acceleration.","marker":"Comisso & Sironi 2022"},{"why":"It formulates the turbulence-damping backreaction prescription that self-regulates acceleration in the third scenario.","marker":"Lemoine et al. 2024"}],"fun_headline_variants":["Self-regulated acceleration explains NGC 1068 neutrinos","Black hole turbulence self-damps to match neutrino flux","Neutrino spectrum from black holes: no tuning needed","When protons damp turbulence, neutrinos match IceCube","Self-regulated acceleration: natural NGC 1068 neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire fit rests on the assumption that the energy-jump kernel measured at one Alfvén speed over short simulation times keeps its shape and scales in width as the square of the Alfvén velocity when applied to slower turbulence and to times up to about 150 coherence lengths; the paper flags this extrapolation as nontrivial.","fun_headline_variants_meta":{"raw":{"variants":["Self-regulated acceleration explains NGC 1068 neutrinos","Black hole turbulence self-damps to match neutrino flux","Neutrino spectrum from black holes: no tuning needed","When protons damp turbulence, neutrinos match IceCube","Self-regulated acceleration: natural NGC 1068 neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000325,"raw_usage":{"total_tokens":1904,"prompt_tokens":1108,"completion_tokens":796,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":724,"completion_tokens_details":{"reasoning_tokens":718}},"tokens_in":724,"tokens_out":796,"duration_ms":7003,"temperature":1.0,"reasoning_tokens":718,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:18:27.440842+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run MHD or particle-in-cell simulations of the energy-jump kernel at an Alfvén velocity near 0.2 of the speed of light for durations up to 150 coherence lengths: if the kernel width does not scale as the assumed velocity-squared law, or if a higher-significance measurement of the NGC 1068 neutrino flux does not track the magnetic energy budget, the central claim is refuted.","supporting_citations":[{"cited_title":"2021, Phys","cited_arxiv_id":null,"evidence_quote":"It introduces the generalized Fermi acceleration picture for particles interacting with intermittent magnetic structures."},{"cited_title":"2022, Phys","cited_arxiv_id":null,"evidence_quote":"It derives the jump kernel used as the acceleration operator in the transport equation."},{"cited_title":"& Sironi, L","cited_arxiv_id":null,"evidence_quote":"It benchmarks the generalized Fermi kernel against particle-in-cell simulations of turbulent particle acceleration."},{"cited_title":"2024, Phys","cited_arxiv_id":null,"evidence_quote":"It formulates the turbulence-damping backreaction prescription that self-regulates acceleration in the third scenario."}],"review_version":1}