{"id":"12d92ca9-6bac-4c02-990f-2591400a9589","arxiv_id":"2601.01999","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"PIC simulations show collisionless shocks in black-hole coronae convert ~10% of energy into non-thermal protons even at low Mach numbers, supporting a corona origin for IceCube neutrinos.","lead":"Simulations of collisionless shocks modeled on supermassive black hole coronae show protons are efficiently accelerated, with about 10% of shock energy going into energetic ions even for weak shocks. This supports the idea that the high-energy neutrinos IceCube sees from Seyfert galaxies like NGC 1068 are produced by shock-accelerated protons in the corona.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The universal 10% ion efficiency and C^-1~0.1 rest on 1D3V PIC with m_R=100 and θ=30°, with no 3D or full-mass convergence test; these values could be geometry- or mass-ratio-dependent.","rationale":"The paper convincingly demonstrates that DSA operates in its 1D3V reduced-mass-ratio simulations: energy grows linearly in time, the NRI dominates, and the non-thermal ion fraction is consistently around 10% in the surveyed parameter space. However, the central astrophysical conclusion—that coronal shocks can satisfy both the spectral and energetic requirements of the IceCube neutrinos—depends on transferring these quantitative values to a real 3D, full-mass, variable-obliquity corona. The reader's weakest assumption identifies exactly this gap: the absence of any convergence test with respect to dimensionality, mass ratio, or obliquity. My re-reading confirms that Section 3 introduces the reduced mass ratio and 1D setup without a validation against 2D/3D or m_R=1836, and Section 5.1 uses the resulting C^-1≈0.1 and E_nt/E_tot≈10% as firm inputs. The uncertainty is load-bearing because the energy-budget margin is only a factor of ~3 (required 3% vs. simulated 10%). A modest change in efficiency or acceleration rate could invalidate the conclusion, even though the order-of-magnitude nature of the application is acknowledged. I therefore support the CONDITIONAL verdict and propose a concrete simulation-based test to settle whether the numbers survive the more realistic setup.","tokens_in":18978,"tokens_out":22342,"duration_ms":222548,"concrete_test":"Run a 2D3V PIC simulation of the reference case T10 (M_s=8, M_A=26, T_i/T_e=10, v_pt=0.25c) and, if feasible, the low-Mach case S2T1, keeping the same dimensionless parameters and box length; separately run a 1D simulation of T10 with m_R=400 or 1836. Compare the downstream non-thermal ion energy fraction (Eq. 8) and the normalized acceleration rate C^-1 (Section 4.2.3). If either changes by more than ~50% from the 1D m_R=100 values, the claimed universalities in Figures 7–8 are not converged and the coronal application in §5.1 must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Central claim (Sections 4.2.3, 4.3, 5.1): non-thermal ion energy fraction ~10% and normalized acceleration rate C^-1≈0.1 are universal in coronal shocks, and this is sufficient for the NGC 1068 neutrino energy budget. These numbers are measured from 1D3V PIC simulations with reduced mass ratio m_R=100 and fixed obliquity θ=30° (Section 3). No convergence study against 3D geometry, larger mass ratio, or varying obliquity is presented; the only numerical check mentioned is spatial/temporal resolution for run A4. In 1D, the magnetic turbulence is restricted to modes along the shock normal, and particles cannot diffuse in transverse directions; in 3D, oblique modes and cross-field transport can alter the saturated turbulence level and the injection efficiency. With m_R=100, the ion/electron scale separation is d_i/d_e=10 instead of 43, which changes the growth and saturation of the NRI/RI responsible for upstream scattering. If the true 3D full-mass efficiency were a factor ≳3 lower, the required non-thermal proton luminosity L_nt_p≈0.03 L_sh (Eq. 13) would exceed the available efficiency and the energetic conclusion in §5.1.2 would fail. The paper itself calls the application an order-of-magnitude assessment (§5.1.2), but the efficiency value is used as a firm quantitative constraint.","agreement_with_reader":"agree"},"referee_report":null,"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuine parameter survey of quasi-parallel, trans-relativistic shocks in the SMBH corona regime, and it gives the field something useful: a reported ion acceleration efficiency around 10% and a normalized acceleration rate C^-1 around 0.1 across a range of sonic Mach, Alfvenic Mach, temperature ratio, and shock speed. The DSA interpretation is coherent: maximum ion energy grows roughly linearly in time, spectra are hard with s~2, and the NRI/RI distinction in the upstream turbulence matches prior work. The NGC 1068 application is an honest order-of-magnitude consistency check: the required non-thermal proton luminosity is about 3% of the shock power, below the 10% efficiency, and the timescale argument for reaching ~100 TeV works for B~100 G. What is new here is the parameter coverage itself: low sonic Mach down to ~2, trans-relativistic speeds, and two-temperature Ti/Te. Earlier PIC DSA studies mostly did not target this corner of parameter space, so the survey fills a real gap. The soft spots, in proportion. First, the headline numbers come from 1D3V runs with mass ratio mR=100 and a single obliquity of 30 degrees. There are no convergence tests against 3D geometry or larger mass ratio. In 1D, magnetic turbulence is restricted to the shock normal and cross-field transport is absent, so the saturation level of the streaming instabilities could change. This might shift the efficiency or acceleration rate by a factor of a few. That does not sink the paper, but it means the word 'universal' is doing more work than the evidence supports. Second, the thermal/non-thermal break is set by fiat: gamma_br,i = 10 E_sh and gamma_br,e = 10 T_de. The 10% number depends on where you cut, and there is no sensitivity test or reported uncertainty. An error bar on the efficiency would be a real improvement. Third, the paper is candid in Section 5.1.2 that the application is an order-of-magnitude assessment, but the abstract and conclusions state the efficiency as a firm result. The energy-budget margin is only about a factor of three; if 3D effects cut the efficiency below ~3%, the NGC 1068 conclusion starts to weaken. That is worth saying explicitly in the discussion. The simulation setup is documented well enough to reproduce: SMILEI, resolution, box sizes, particle numbers. The citation pattern is appropriate; the self-citations are to the authors' own coronal work and to the Caprioli/Spitkovsky/Gupta lineage, which is normal in this subfield. For anyone working on AGN neutrino hadronic models or PIC shock acceleration, this is worth reading carefully. I would send it to peer review. The main requests should be convergence tests or at least an explicit estimate of dimensionality and mass-ratio systematics, plus a robustness check on the break definitions. I would not desk-reject it.","headline":"A serious PIC survey makes a plausible case for ~10% ion acceleration efficiency in coronal-like quasi-parallel shocks, but the 'universal' rate rests on 1D reduced-mass runs and should be treated as provisional until convergence is shown.","tokens_in":816,"tokens_out":875,"would_cite":true,"duration_ms":31623,"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":"Collisionless shocks in black hole coronae efficiently accelerate protons to the energies that produce IceCube neutrinos.","keywords":["diffusive shock acceleration","particle-in-cell simulation","black hole corona","active galactic nuclei","high-energy neutrinos","NGC 1068","cosmic ray acceleration","collisionless shock"],"falsifier":"A single 3D or full-mass-ratio PIC simulation of a quasi-parallel coronal shock in which the downstream non-thermal ion energy fraction deviates strongly from ~10% (e.g., below a few percent) would falsify the claimed universality; alternatively, a future measurement showing that the neutrino spectrum of NGC 1068 requires a proton luminosity exceeding 10% of the shock power, or a proton cutoff below ~30 TeV contradicting the derived acceleration timescale, would falsify the application to NGC 1068.","tokens_in":18860,"feed_emoji":"🌌","tokens_out":8497,"duration_ms":76584,"temperature":0.7,"pith_summary":"The paper aims to establish that diffusive shock acceleration in hot, collisionless coronae around supermassive black holes can accelerate protons to the ~100 TeV energies required to produce the high-energy neutrinos IceCube detects from Seyfert galaxies like NGC 1068. Using first-principles particle-in-cell simulations over a wide parameter range, it finds that about 10% of the shock's kinetic energy consistently goes into non-thermal ions, even at sonic Mach numbers as low as ~2, and that the normalized proton acceleration rate varies by less than a factor of two. Because the required non-thermal proton luminosity for NGC 1068 is only a few percent of the estimated shock power, the authors conclude that coronal shocks can satisfy both the spectral and energetic requirements of the neutrino observations while the accompanying gamma-rays are absorbed. This would resolve the puzzle of why neutrinos are seen without a corresponding gamma-ray flux.","feed_headline":"Coronal shocks give protons 10% energy, powering IceCube neutrinos","feed_subtitle":"Kinetic simulations show ~10% of shock energy becomes protons, enough to feed IceCube's neutrinos from Seyferts.","key_machinery":"The load-bearing mechanism is diffusive shock acceleration (DSA) operating in quasi-parallel shocks: back-streaming protons excite upstream magnetic turbulence (predominantly the non-resonant Bell instability), which scatters particles so they repeatedly cross the shock and gain energy. The quantitative anchors are the normalized acceleration coefficient C^{-1} ≈ 0.1 and the ~10% non-thermal ion energy fraction, both extracted from the downstream spectra and used to scale the acceleration timescale t_acc ≈ 11.1 s (E_p/100 TeV)(C^{-1}/0.1)^{-1}(β_sh/0.32)^{-2}(B_0/100 G)^{-1}.","core_discovery":"The central discovery is that quasi-parallel collisionless shocks in black-hole-corona conditions consistently channel about 10% of the shock kinetic energy into non-thermal ions, with a normalized acceleration rate C^{-1} ≈ 0.1 that varies by less than a factor of two across all surveyed parameters. This persists even at sonic Mach number M_s ≈ 2, previously considered too weak for efficient acceleration. Electron acceleration is subdominant and variable. From these values, the authors derive an acceleration timescale for NGC 1068 shorter than all cooling timescales up to ~1 PeV and an energy budget requiring only ~3% of the shock power to match the observed neutrino luminosity.","pith_inferences":["A direct multi-dimensional, full-mass-ratio PIC simulation would test whether the 10% ion efficiency and C^{-1} ≈ 0.1 survive outside 1D geometry and reduced mass ratio; if they do, the same scaling could be applied to other Seyferts with coronal shocks, predicting their neutrino fluxes from measured coronal parameters.","The predicted proton spectrum with slope ~2 and a Hillas-limited cutoff near 0.1–1 PeV implies a spectral softening in the neutrino flux above ~10 TeV, which future high-energy neutrino detectors could observe.","If the electron acceleration suppression at high T_i/T_e holds, radio/millimeter synchrotron from primary electrons in these coronae should be subdominant; upcoming millimeter observations could distinguish this DSA scenario from reconnection or stochastic-acceleration models that predict stronger leptonic emission.","The acceleration timescale scaling could be applied to coronal shocks in other AGN classes (e.g., radio-loud or changing-look AGN) to assess whether they also contribute to the diffuse neutrino background."],"forward_implications":["Even weak shocks with sonic Mach number ~2 convert about 10% of shock energy into non-thermal protons, so low-Mach coronal shocks remain viable hadronic accelerators.","The normalized acceleration rate C^{-1} ≈ 0.1 gives a concrete, simulation-calibrated timescale; at B0 = 100 G and β_sh = 0.32, protons reach 100 TeV in about 11 seconds, shorter than cooling and free-fall times.","The required non-thermal proton luminosity for NGC 1068's neutrino signal is ~2×10^43 erg/s, only ~3% of the shock power, within the measured 10% efficiency.","Because electrons receive <1% of the energy and even less for high T_i/T_e, primary lepton emission is weak; the observed gamma-ray deficit can be explained by pair-production opacity rather than by a lack of hadronic acceleration.","The nearly parameter-independent C^{-1} allows hadronic models to replace arbitrary acceleration-efficiency assumptions with a physically derived value."],"fun_headline_variants":["10% of shock energy goes to protons even in weak coronal shocks","Black hole corona shocks: 10% energy to protons, feeding IceCube","Coronal shocks convert 10% energy to protons, explaining IceCube neutrinos","Weak shocks in black hole coronae still make protons for neutrinos","Even Mach-2 shocks in black hole coronae channel 10% energy to protons"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative transfer of the 10% efficiency and C^{-1} ≈ 0.1 to real coronae rests on the unverified premise that 1D3V particle-in-cell simulations with a reduced ion-to-electron mass ratio of 100 and fixed 30° magnetic obliquity reproduce the true diffusive shock acceleration efficiency and rate of a 3D, full-mass coronal shock.","fun_headline_variants_meta":{"raw":{"variants":["10% of shock energy goes to protons even in weak coronal shocks","Black hole corona shocks: 10% energy to protons, feeding IceCube","Coronal shocks convert 10% energy to protons, explaining IceCube neutrinos","Weak shocks in black hole coronae still make protons for neutrinos","Even Mach-2 shocks in black hole coronae channel 10% energy to protons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000787,"raw_usage":{"total_tokens":3325,"prompt_tokens":775,"completion_tokens":2550,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":519,"completion_tokens_details":{"reasoning_tokens":2449}},"tokens_in":519,"tokens_out":2550,"duration_ms":18846,"temperature":1.0,"reasoning_tokens":2449,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T12:37:57.966843+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single 3D or full-mass-ratio PIC simulation of a quasi-parallel coronal shock in which the downstream non-thermal ion energy fraction deviates strongly from ~10% (e.g., below a few percent) would falsify the claimed universality; alternatively, a future measurement showing that the neutrino spectrum of NGC 1068 requires a proton luminosity exceeding 10% of the shock power, or a proton cutoff below ~30 TeV contradicting the derived acceleration timescale, would falsify the application to NGC 1068.","supporting_citations":[],"review_version":1}