Pith. sign in

REVIEW 1 major objections 4 minor 1 cited by

A sector-wise energy reduction finds collisional neutrino flavor instabilities without multi-energy integrals or spurious singularities.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review

2026-07-13 14:37 UTC pith:6LH4V6NI

load-bearing objection Solid, usable CFI reduction that fixes the known singularities of methods A/B and works past k=0; near-mode errors are real but already diagnosed and not fatal to the identification claim. the 1 major comments →

arxiv 2604.01096 v1 pith:6LH4V6NI submitted 2026-04-01 astro-ph.HE hep-ph

Approximate Energy-Integration Method for Identifying Collisional Neutrino Flavor Instabilities

Jiabao Liu , Hiroki Nagakura This is my paper
classification astro-ph.HE hep-ph
keywords collisional flavor instabilityneutrino oscillationsdispersion relationenergy integrationcore-collapse supernovaebinary neutron star mergersquantum kinetic equations
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

In dense environments such as core-collapse supernovae and neutron-star mergers, neutrinos can undergo collisional flavor instabilities whose growth rates are set by the energy-dependent collision rates. Finding those modes from the exact dispersion relation requires costly multi-dimensional phase-space integrals, so earlier surveys relied on crude energy averages that are limited to homogeneous modes and can produce unphysical rates or large errors. This paper constructs a new reduction, method C, that splits the signed neutrino and antineutrino spectra into positive and negative sectors, defines well-behaved effective densities and collision rates for each sector, and collapses the energy integrals to a simple rational dispersion relation. Direct comparison with multi-energy solutions shows that the reduced relation recovers both real frequencies and growth rates for isotropic and anisotropic distributions, and for both homogeneous and inhomogeneous wave numbers. The method therefore supplies a practical, accurate tool for scanning large simulation data sets for collisional instabilities.

Core claim

Method C, obtained by decomposing the energy spectra into positive and negative sectors and expanding the collision denominators to first order, yields a reduced dispersion relation that remains free of singularities and reproduces the exact real frequencies and growth rates of collisional flavor instabilities across isotropic, anisotropic, homogeneous and inhomogeneous regimes.

What carries the argument

The positive/negative sector decomposition of Δf(E) and Δf̄(E), which produces four non-negative moments N± and G± and two effective collision rates Γ± = G±/N±; these replace the multi-energy integrals by the compact rational form N+/(ω+iΓ+) − N−/(ω+iΓ−).

Load-bearing premise

The first-order expansion of the energy-dependent denominators stays accurate only when the characteristic frequency scale of a mode is larger than the effective collision rate; modes sitting near the origin systematically violate that condition.

What would settle it

Apply method C and the exact multi-energy dispersion relation to a neutrino distribution whose near-mode growth rate is comparable to the microscopic collision rate; if the reduced growth rate deviates by more than a few tens of percent while far or resonant modes remain accurate, the expansion premise fails for that class of modes.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 4 minor

Summary. The paper develops Method C, an approximate energy-integration scheme for collisional neutrino flavor instabilities (CFI). Starting from the linearized quantum kinetic equations and the polarization tensor (Eqs. 13–15), the authors decompose signed neutrino and antineutrino spectra into positive and negative sectors (Eqs. 30–34), expand the energy-dependent denominators to first order, and re-sum into sector-wise effective densities and collision rates (Eqs. 36–39). The resulting reduced dispersion relations (Eqs. 40–45 for isotropic k=0; schematic form Eq. 59 and angle-dependent rates Eq. 62 for anisotropic and k eq0 cases) are algebraically simple and free of the signed-integral singularities that afflict Method A and the uncontrolled averaging of Method B. Extensive numerical comparisons against multi-energy exact roots (Secs. V A–C, Figs. 5–11, Appendix A) show that Method C reproduces real frequencies and growth rates for isotropic and axisymmetric distributions, resonance and non-resonance regimes, and both homogeneous and inhomogeneous modes, with the largest residual errors confined to near modes.

Significance. If the accuracy claims hold, Method C supplies a practical, scalable tool for systematic CFI surveys in CCSN and BNSM simulations, where full multi-energy multi-angle root searches remain expensive. The construction is derived rather than purely empirical, eliminates the pathologies of Methods A and B that recent literature has already flagged, and is shown to work beyond the isotropic k=0 limit that restricted earlier schemes. The side-by-side exact-versus-reduced comparisons across a controlled suite of models, together with an explicit geometric diagnosis of the near-mode limitation (Sec. IV C), constitute a solid and falsifiable contribution that the community can adopt or further refine.

major comments (1)
  1. [Sec. IV C; Abstract; Figs. 5, 9, 10] Sec. IV C and the empirical criterion after Fig. 2 correctly identify that the first-order re-summation is controlled only when |ω−kv| ≳ Γ_eff. Near modes (Re ω ≈ 0) systematically violate this condition; Figs. 5, 9 (left) and 10 (left) show fractional errors of tens of percent in Im ω precisely there. The Abstract and Sec. VI claim “good performance … across a wide range of regimes” and “accurate estimates of both the real frequencies and growth rates.” Because mode identification (existence and rough scale of Im ω) is the stated practical goal, the residual accuracy may still be adequate, but the manuscript should either (i) quantify a clear acceptance threshold (e.g., relative error in Im ω for near modes) or (ii) soften the absolute language so that the near-mode caveat is visible in the abstract-level claim rather than only in the body.
minor comments (4)
  1. [Fig. 1] Fig. 1 is schematic but the four labeled areas are never mapped explicitly onto the definitions of N± and G± in Eqs. (36)–(37); a short caption sentence would help readers who skip the algebra.
  2. [Introduction] The phrase “we assess” appears with a lowercase “w” at the start of a sentence in the Introduction (page 2, column 1).
  3. [Sec. V A] In Sec. V A the energy grid is stated as 500 logarithmic points on [10^{-4}, 100] MeV; a one-sentence remark on whether the same grid is used for the anisotropic and k eq0 tests would remove a minor ambiguity.
  4. [Figs. 7–8, 10] The Isotropized-CFI comparison in Figs. 7–8 and 10 is useful; labeling the green/red symbols more consistently with the IP/IB terminology of Sec. III would improve readability.

Circularity Check

0 steps flagged

No circularity: Method C is an explicit first-order reduction of the multi-energy DR, validated by independent exact multi-energy roots on the same models.

full rationale

The derivation chain is self-contained. Sec. IV A–B starts from the exact polarization integrals (Eqs. 14, 48–51), applies a positive/negative sector split of the signed spectra (Eqs. 30–34), expands the energy-dependent denominators to first order, and re-sums to obtain effective N±, G±, Γ± (Eqs. 36–38) and the reduced rational form (Eq. 39). The resulting algebraic DR (Eqs. 40–44) is solved in closed form or by polynomial root-finding. Validation (Sec. V) consists of direct numerical comparison of these reduced roots against independently computed multi-energy exact roots of the same det Π=0 on identical isotropic/anisotropic spectra and collision rates (Figs. 5–11); no free parameters are fitted to the target Im ω or Re ω. Self-citations supply background on CFI phenomenology and prior approximate schemes A/B, but the accuracy claim rests solely on the internal exact-vs-reduced comparisons performed in this work. The acknowledged limitation for near modes (Sec. IV C) is a controlled expansion error, not a circular reduction. No step reduces a claimed prediction to its own input by construction.

Axiom & Free-Parameter Ledger

3 free parameters · 6 axioms · 0 invented entities

The central claim rests on standard QKE linearization plus two modeling choices (relaxation collisions; neglect of vacuum/matter Hamiltonians) and on the paper’s own first-order sector expansion. Numerical experiment parameters are fixed test inputs, not fitted constants that define the claim. No new physical entities are postulated.

free parameters (3)
  • Collision-rate normalizations Γ_να,0 = 1.8e-5, 0.8e-5, 0.2e-5 cm^{-1}
    Fixed by hand to 1.8, 0.8, 0.2 × 10^{-5} cm^{-1} for the numerical experiments; they set the absolute scale of growth rates in the test suite but are not fitted to external data.
  • Spectral normalizations g_να and Legendre scalings g_να,ℓ = varied in [0,1]; c_ν̄e ≃ 0.938 for resonance
    Hand-chosen to control energy crossings and angular moments in Models I/II and NF/F series; used only to generate controlled test spectra.
  • Fermi–Dirac temperatures and chemical potentials = T=3.4–5.6 MeV, μ=0–3 MeV
    Fixed thermodynamic parameters of the test spectra (Table I); not inferred from observation.
axioms (6)
  • domain assumption Linear stability analysis of the quantum kinetic equations with plane-wave ansatz yields the polarization-tensor dispersion relation det Π=0.
    Standard framework for CFI/FFI (Sec. II, Eqs. 6–15); assumed throughout.
  • domain assumption Collision term is a flavor-diagonal relaxation toward the initial flavor eigenstate with rates Γ_να(E).
    Adopted explicitly in Eq. (5); full collision kernels are not used.
  • domain assumption Vacuum and matter Hamiltonians may be neglected for CFI identification.
    Stated in Sec. II; they do not drive CFI in this framework.
  • ad hoc to paper First-order expansion of energy-dependent denominators followed by re-summation into sector-wise effective rates is a controlled approximation when |ω−kv| ≳ Γ_eff.
    Core construction of method C (Eqs. 35–39, 61–62; Sec. IV C).
  • ad hoc to paper Positive and negative spectral sectors may be grouped as ([Δf]+,[Δf̄]−) and ([Δf]−,[Δf̄]+) without loss of the leading CFI physics.
    Defining organizational step of method C (Sec. IV A, Fig. 1).
  • domain assumption Neutrinos are massless and ultra-relativistic (|v|=1).
    Stated in Sec. II; standard for this class of analyses.

reviewed 2026-07-13 · how reviews work

0 comments
read the original abstract

We present an approximate energy-integration method for identifying collisional neutrino flavor instabilities. Direct evaluation of the dispersion relation requires multi-dimensional integrals over neutrino phase space, making systematic searches for unstable modes in numerical models of core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs) computationally expensive. In the literature there are some approximate schemes, but they are largely restricted to the homogeneous limit and can exhibit inaccuracies as reported in recent studies. In the current paper, we clarify the origin of the limitations in previous schemes and provide a better approximation method that robustly preserves the key physics of spectral asymmetries and collision rates. It yields a reduced dispersion relation that is inexpensive to evaluate. Comparison with exact solutions demonstrates that our new approximate method shows a good performance in computing both real frequencies and growth rates across a wide range of regimes, including isotropic and anisotropic neutrino distributions for both homogeneous and inhomogeneous modes. This provides a practical, accurate, and scalable framework for identifying collisional flavor instabilities in high-energy astrophysical simulations such as CCSNe and BNSMs.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Collective neutrino-antineutrino pair oscillations

    hep-ph 2026-04 unverdicted novelty 7.0

    In anisotropic neutrino gases, νν-bar pairing instabilities emerge when the excessive pair-occupation number distribution changes sign, producing pair conversions at growth rates comparable to fast flavor instabilities.

Reference graph

Works this paper leans on

61 extracted references · 4 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Capozzi and N

    F. Capozzi and N. Saviano, Universe8, 10.3390/uni- verse8020094 (2022)

    doi:10.3390/uni- 2022

  2. [2]

    Tamborra and S

    I. Tamborra and S. Shalgar, Annual Review of Nuclear and Particle Science71, 165 (2021)

    2021

  3. [3]

    Johns, S

    L. Johns, S. Richers, and M.-R. Wu, Annual Review of Nuclear and Particle Science75, 399 (2025)

    2025

  4. [4]

    M. C. Volpe, Rev. Mod. Phys.96, 025004 (2024)

    2024

  5. [5]

    YAMADA, H

    S. YAMADA, H. NAGAKURA, R. AKAHO, A. HARADA, S. FURUSAWA, W. IWAKAMI, H. OKAWA, H. MATSUFURU, and K. SUMIYOSHI, Proceedings of the Japan Academy, Series B100, 190 (2024)

    2024

  6. [6]

    Fischer, G

    T. Fischer, G. Guo, K. Langanke, G. Mart´ ınez-Pinedo, Y.-Z. Qian, and M.-R. Wu, Progress in Particle and Nu- clear Physics137, 104107 (2024)

    2024

  7. [7]

    G. G. Raffelt, H.-T. Janka, and D. F. G. Fior- illo, Neutrinos from core-collapse supernovae (2026), arXiv:2509.16306 [astro-ph.HE]

    arXiv 2026

  8. [8]

    H. T. Janka, Long-term multidimensional models of core- collapse supernovae: Progress and challenges (2026), 18 arXiv:2502.14836 [astro-ph.HE]

    arXiv 2026

  9. [9]

    Samuel, Phys

    S. Samuel, Phys. Rev. D48, 1462 (1993)

    1993

  10. [10]

    Pantaleone, Physics Letters B287, 128 (1992)

    J. Pantaleone, Physics Letters B287, 128 (1992)

    1992

  11. [11]

    H. Duan, G. M. Fuller, and Y.-Z. Qian, Annual Review of Nuclear and Particle Science60, 569 (2010)

    2010

  12. [12]

    Duan and J

    H. Duan and J. P. Kneller, Journal of Physics G: Nuclear and Particle Physics36, 113201 (2009)

    2009

  13. [13]

    R. F. Sawyer, Phys. Rev. D72, 045003 (2005)

    2005

  14. [14]

    R. F. Sawyer, Phys. Rev. D79, 105003 (2009)

    2009

  15. [15]

    Richers and M

    S. Richers and M. Sen, Fast flavor transformations, in Handbook of Nuclear Physics, edited by I. Tanihata, H. Toki, and T. Kajino (Springer Nature Singapore, Sin- gapore, 2020) pp. 1–17

    2020

  16. [16]

    Wu and I

    M.-R. Wu and I. Tamborra, Phys. Rev. D95, 103007 (2017)

    2017

  17. [17]

    Abbar, H

    S. Abbar, H. Duan, K. Sumiyoshi, T. Takiwaki, and M. C. Volpe, Phys. Rev. D100, 043004 (2019)

    2019

  18. [18]

    Morinaga, H

    T. Morinaga, H. Nagakura, C. Kato, and S. Yamada, Phys. Rev. Res.2, 012046 (2020)

    2020

  19. [19]

    Abbar, H

    S. Abbar, H. Duan, K. Sumiyoshi, T. Takiwaki, and M. C. Volpe, Phys. Rev. D101, 043016 (2020)

    2020

  20. [20]

    Nagakura, A

    H. Nagakura, A. Burrows, L. Johns, and G. M. Fuller, Phys. Rev. D104, 083025 (2021)

    2021

  21. [21]

    Richers, Phys

    S. Richers, Phys. Rev. D106, 083005 (2022)

    2022

  22. [22]

    Grohs, S

    E. Grohs, S. Richers, S. M. Couch, F. Foucart, J. P. Kneller, and G. McLaughlin, Physics Letters B846, 138210 (2023)

    2023

  23. [23]

    Richers, H

    S. Richers, H. Duan, M.-R. Wu, S. Bhattacharyya, M. Za- izen, M. George, C.-Y. Lin, and Z. Xiong, Phys. Rev. D 106, 043011 (2022)

    2022

  24. [24]

    Froustey, S

    J. Froustey, S. Richers, E. Grohs, S. D. Flynn, F. Foucart, J. P. Kneller, and G. C. McLaughlin, Phys. Rev. D109, 043046 (2024)

    2024

  25. [25]

    Mukhopadhyay, J

    P. Mukhopadhyay, J. Miller, and G. C. McLaughlin, The Astrophysical Journal974, 110 (2024)

    2024

  26. [26]

    Xiong, M.-R

    Z. Xiong, M.-R. Wu, M. George, C.-Y. Lin, N. K. Largani, T. Fischer, and G. Mart´ ınez-Pinedo, Phys. Rev. D109, 123008 (2024)

    2024

  27. [27]

    Xiong, M.-R

    Z. Xiong, M.-R. Wu, M. George, and C.-Y. Lin, Phys. Rev. Lett.134, 051003 (2025)

    2025

  28. [28]

    D. F. G. Fiorillo and G. G. Raffelt, J. High Energy Phys. 08, 225, arXiv:2403.12143 [hep-ph]

    Pith/arXiv arXiv

  29. [29]

    D. F. G. Fiorillo and G. G. Raffelt, J. High Energy Phys. 12, 205, arXiv:2406.12345 [hep-ph]

    Pith/arXiv arXiv

  30. [30]

    J. Liu, L. Johns, H. Nagakura, M. Zaizen, and S. Ya- mada, Dynamical equilibria of fast neutrino flavor con- version (2025), arXiv:2509.26418 [astro-ph.HE]

    arXiv 2025

  31. [31]

    Froustey, F

    J. Froustey, F. Foucart, C. Hall, J. P. Kneller, D. Kundu, Z. Lin, G. C. McLaughlin, and S. Richers, Phys. Rev. D 113, 063050 (2026)

    2026

  32. [32]

    Nagakura, T

    H. Nagakura, T. Morinaga, C. Kato, and S. Yamada, The Astrophysical Journal886, 139 (2019)

    2019

  33. [33]

    Nagakura, Phys

    H. Nagakura, Phys. Rev. Lett.130, 211401 (2023)

    2023

  34. [34]

    Nagakura and M

    H. Nagakura and M. Zaizen, Phys. Rev. Lett.129, 261101 (2022)

    2022

  35. [35]

    Akaho, J

    R. Akaho, J. Liu, H. Nagakura, M. Zaizen, and S. Ya- mada, Phys. Rev. D109, 023012 (2024)

    2024

  36. [36]

    Zaizen and H

    M. Zaizen and H. Nagakura, Phys. Rev. D107, 103022 (2023)

    2023

  37. [37]

    D. F. G. Fiorillo and G. G. Raffelt, Phys. Rev. Lett.133, 221004 (2024)

    2024

  38. [38]

    M.-R. Wu, M. George, C.-Y. Lin, and Z. Xiong, Phys. Rev. D104, 103003 (2021)

    2021

  39. [39]

    Richers, D

    S. Richers, D. Willcox, and N. Ford, Phys. Rev. D104, 103023 (2021)

    2021

  40. [40]

    Johns, H

    L. Johns, H. Nagakura, G. M. Fuller, and A. Burrows, Phys. Rev. D101, 043009 (2020)

    2020

  41. [41]

    Bhattacharyya and B

    S. Bhattacharyya and B. Dasgupta, Phys. Rev. Lett.126, 061302 (2021)

    2021

  42. [42]

    Johns, Phys

    L. Johns, Phys. Rev. Lett.130, 191001 (2023)

    2023

  43. [43]

    Xiong, L

    Z. Xiong, L. Johns, M.-R. Wu, and H. Duan, Phys. Rev. D108, 083002 (2023)

    2023

  44. [44]

    J. Liu, M. Zaizen, and S. Yamada, Phys. Rev. D107, 123011 (2023)

    2023

  45. [45]

    Padilla-Gay, I

    I. Padilla-Gay, I. Tamborra, and G. G. Raffelt, Phys. Rev. D106, 103031 (2022)

    2022

  46. [46]

    Johns and Z

    L. Johns and Z. Xiong, Phys. Rev. D106, 103029 (2022)

    2022

  47. [47]

    Lin and H

    Y.-C. Lin and H. Duan, Phys. Rev. D107, 083034 (2023)

    2023

  48. [48]

    D. F. G. Fiorillo, I. Padilla-Gay, and G. G. Raffelt, Phys. Rev. D109, 063021 (2024)

    2024

  49. [49]

    Xiong, M.-R

    Z. Xiong, M.-R. Wu, G. Mart´ ınez-Pinedo, T. Fischer, M. George, C.-Y. Lin, and L. Johns, Phys. Rev. D107, 083016 (2023)

    2023

  50. [50]

    J. Liu, R. Akaho, A. Ito, H. Nagakura, M. Zaizen, and S. Yamada, Phys. Rev. D108, 123024 (2023)

    2023

  51. [51]

    Shalgar and I

    S. Shalgar and I. Tamborra, Phys. Rev. D109, 103011 (2024)

    2024

  52. [52]

    Zaizen, S

    M. Zaizen, S. Richers, H. Nagakura, H. Suzuki, and C. Kato, Inspecting neutrino flavor instabilities during proto-neutron star cooling phase in supernova: I. spher- ically symmetric model (2024), arXiv:2407.20548 [astro- ph.HE]

    Pith/arXiv arXiv 2024

  53. [53]

    T. Wang, H. Nagakura, L. Johns, and A. Burrows, Phys. Rev. D112, 063039 (2025)

    2025

  54. [54]

    Nagakura, K

    H. Nagakura, K. Sumiyoshi, S. Fujibayashi, Y. Sekiguchi, and M. Shibata, Neutrino flavor instabilities in a binary neutron star merger remnant: Roles of a long-lived hy- permassive neutron star (2025), arXiv:2504.20143 [astro- ph.HE]

    Pith/arXiv arXiv 2025

  55. [55]

    Banerjee, A

    A. Banerjee, A. Dighe, and G. Raffelt, Phys. Rev. D84, 053013 (2011)

    2011

  56. [56]

    Airen, F

    S. Airen, F. Capozzi, S. Chakraborty, B. Dasgupta, G. Raffelt, and T. Stirner, Journal of Cosmology and Astroparticle Physics2018(12), 019

  57. [57]

    Sigl and G

    G. Sigl and G. Raffelt, Nuclear Physics B406, 423 (1993)

    1993

  58. [58]

    D. N. Blaschke and V. Cirigliano, Phys. Rev. D94, 033009 (2016)

    2016

  59. [59]

    S. A. Richers, G. C. McLaughlin, J. P. Kneller, and A. Vlasenko, Phys. Rev. D99, 123014 (2019)

    2019

  60. [60]

    Froustey, C

    J. Froustey, C. Pitrou, and M. C. Volpe, Journal of Cos- mology and Astroparticle Physics2020(12), 015

  61. [61]

    Nagakura, Phys

    H. Nagakura, Phys. Rev. D106, 063011 (2022). Appendix A: Coexistence of F ast and Collisional Modes atk= 0 In this appendix, we present two additional axisym- metric models that exhibit both FFIs and CFIs atk= 0. The models are constructed in the same manner as those in Sec. V B, and their parameters are summarized in Ta- ble III. The approximate and exac...

    2022

This paper was first reviewed by grok-4.5 on July 13, 2026.