pith. sign in

arxiv: 2601.18051 · v2 · submitted 2026-01-26 · 🌌 astro-ph.HE · nucl-th

Neutrino opacities in magnetic fields for binary neutron star merger simulations

Mia Kumamoto , Catherine Welch This is my paper

Pith reviewed 2026-05-16 11:30 UTC · model grok-4.3

classification 🌌 astro-ph.HE nucl-th
keywords neutrino opacitiesmagnetic fieldsneutron star mergersLandau quantizationanomalous magnetic momentsneutrino interactionsbinary mergersejecta transport
0
0 comments X

The pith

Approximate neutrino interaction rates in strong magnetic fields are derived for binary neutron star merger simulations, including Landau quantization and anomalous magnetic moments with errors of order sqrt(T/M).

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

The paper develops approximate formulas for neutrino interaction rates when strong magnetic fields are present, as expected in the ejecta of binary neutron star mergers. These formulas incorporate quantum effects from Landau quantization of charged particle orbits and the anomalous magnetic moments of nucleons, while keeping the approximation error at order sqrt(T/M). The work also identifies a production mechanism in which individual neutrons can emit low-energy neutrino-antineutrino pairs even in low-density regions. Such rates matter because magnetic fields up to 10^17 G can alter neutrino transport and flavor evolution, which in turn shape the composition and observable signals of merger ejecta. Without manageable approximations, including these magnetic effects in simulations has been computationally prohibitive.

Core claim

We give approximate interaction rates for neutrinos in the presence of strong magnetic fields, including the effects of Landau quantization and anomalous magnetic moments with errors of order sqrt(T/M). We also comment on a neutrino production channel from individual neutrons that can produce low-energy nu nubar pairs even at low density.

What carries the argument

Approximate neutrino interaction rates that fold in Landau quantization of charged-particle states and anomalous magnetic moments, controlled to errors of order sqrt(T/M).

If this is right

  • Merger simulations can incorporate magnetic-field corrections to neutrino opacities without prohibitive computational cost.
  • Neutrino transport and flavor evolution in the ejecta become sensitive to field strength through the adjusted rates.
  • Low-density regions gain an additional source of low-energy neutrino pairs from single-neutron processes.
  • The approximations apply across the range of conditions typically encountered in post-merger ejecta.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Improved rates could alter predictions for the energy deposition that powers kilonova light curves.
  • The same framework might be adapted to study neutrino behavior in other strong-field environments such as magnetar atmospheres.
  • Embedding the rates in full transport codes would allow quantitative assessment of how magnetic fields shift r-process yields.
  • The low-density pair-production channel may affect early-time neutrino signals observable by future detectors.

Load-bearing premise

The error bound of order sqrt(T/M) remains valid under the density, temperature, and field-strength conditions realized in merger ejecta.

What would settle it

A side-by-side numerical comparison of the approximate rates against exact quantum-field calculations at a representative point in density, temperature, and magnetic field strength drawn from merger ejecta would directly test whether the claimed error scaling holds.

Figures

Figures reproduced from arXiv: 2601.18051 by Catherine Welch, Mia Kumamoto.

Figure 1
Figure 1. Figure 1: FIG. 1: Charged current opacities for a range of temperatures, magnetic field strengths, and neutrino energies. Thin [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Neutral current opacities for a range of temperatures, magnetic field strengths, and neutrino energies. Thin [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Proton scattering kernel at [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Total opacities as a function of incoming neutrino direction as compared to the magnetic field axis. [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Emissivity as a function of density for Urca and synchrotron neutrino emission. Temperature and proton [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Charged current opacities for non-degenerate nucleons comparing Monte Carlo calculation of the full [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Neutral current differential and total opacities for scattering on neutrons comparing Monte Carlo [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Neutral current differential and total opacities for scattering on protons comparing Monte Carlo calculation [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
read the original abstract

Neutrino interactions play a central role in transport and flavor evolution in the ejecta of binary neutron star mergers. Simulations suggest that neutron star mergers may produce magnetic fields as strong as $10^{17}$ G, but computational difficulties have hampered the inclusion of magnetic field effects in neutrino interaction rates. In this paper we give approximate interaction rates for neutrinos in the presence of strong magnetic fields, including the effects of Landau quantization and anomalous magnetic moments with errors of order $\sqrt{T/M}$. We also comment on a neutrino production channel from individual neutrons that can produce low-energy $\nu \bar{\nu}$ pairs even at low density.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper derives approximate neutrino interaction rates and opacities for use in binary neutron star merger simulations under strong magnetic fields up to 10^17 G. The approximations incorporate Landau quantization of charged fermions and nucleon anomalous magnetic moments, with asserted errors of order √(T/M). The work also identifies a low-density neutrino-antineutrino pair production channel from individual neutrons.

Significance. If the error bound holds under merger conditions, the rates would enable direct inclusion of magnetic effects in neutrino transport modules without prohibitive computational cost. This addresses a known gap in current simulations, where magnetic fields are often omitted despite their potential impact on flavor evolution and ejecta dynamics. The controlled expansion in √(T/M) while treating magnetic quantization non-perturbatively is a useful technical advance if validated.

major comments (2)
  1. [§3] §3 (error analysis): The O(√(T/M)) error bound on the phase-space integrals is asserted after expanding in T/M while keeping Landau levels exact, but no explicit numerical or analytic check is provided for the regime eB ≳ T^2 with B ∼ 10^17 G and T ∼ few MeV. This bound is load-bearing for the claim that the rates are ready for merger simulations.
  2. [§2.3] §2.3 (dispersion relations): The treatment of anomalous magnetic moments in the strong-field dispersion relations omits possible O(√(eB)/M) corrections to the density of states that are not parametrically suppressed by √(T/M); if present, they would invalidate the stated accuracy precisely at the field strengths highlighted in the abstract.
minor comments (2)
  1. [Figure 2] Figure 2 caption: the plotted opacity curves lack error bands or shaded regions indicating the claimed √(T/M) uncertainty, making it difficult to assess the approximation visually.
  2. [Introduction] Introduction: the discussion of prior magnetic-field neutrino opacity calculations (e.g., refs. on Landau-level summation) is brief; adding one or two sentences on how the present expansion differs from those works would improve context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below.

read point-by-point responses

  1. Referee: [§3] §3 (error analysis): The O(√(T/M)) error bound on the phase-space integrals is asserted after expanding in T/M while keeping Landau levels exact, but no explicit numerical or analytic check is provided for the regime eB ≳ T^2 with B ∼ 10^17 G and T ∼ few MeV. This bound is load-bearing for the claim that the rates are ready for merger simulations.

    Authors: We agree that an explicit validation would strengthen the paper. In the revised manuscript we will add a numerical check comparing the approximated phase-space integrals to direct quadrature for B = 10^17 G and T = 2–10 MeV, confirming that the relative error remains O(√(T/M)) in this regime. revision: yes

  2. Referee: [§2.3] §2.3 (dispersion relations): The treatment of anomalous magnetic moments in the strong-field dispersion relations omits possible O(√(eB)/M) corrections to the density of states that are not parametrically suppressed by √(T/M); if present, they would invalidate the stated accuracy precisely at the field strengths highlighted in the abstract.

    Authors: We disagree. Section 2.3 solves the Dirac equation exactly with the anomalous magnetic moment term included in the Landau-level energies; the density of states follows directly from the resulting quantization condition with no further approximation. Any O(√(eB)/M) contributions are therefore retained non-perturbatively in B, while the √(T/M) expansion applies only to the thermal integrals over the distribution functions. revision: no

Circularity Check

0 steps flagged

Derivation of approximate neutrino opacities is self-contained

full rationale

The paper presents approximate neutrino interaction rates obtained by incorporating Landau quantization and anomalous magnetic moments into standard phase-space integrals, with an asserted error of order √(T/M). No load-bearing step reduces by construction to a fitted input, self-definition, or self-citation chain; the error bound is stated as an analytical estimate from the expansion rather than a quantity forced by the result itself. The derivation relies on conventional QFT techniques for fermions in magnetic fields and does not rename or smuggle in prior results as new predictions. This is the normal, non-circular outcome for an approximation paper whose central expressions remain independent of their own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard quantum-field-theory treatment of charged particles in uniform magnetic fields plus the assumption that merger ejecta reach the stated field strengths; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Magnetic fields of order 10^17 G are generated in binary neutron star merger ejecta.
    Stated as suggested by prior simulations.

pith-pipeline@v0.9.0 · 5394 in / 1116 out tokens · 28347 ms · 2026-05-16T11:30:23.961119+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · 1 internal anchor

  1. [1]

    Charged current For charged current processes, we drop terms that cancel when the electron momentum is integrated, leaving the following reduced matrix element. M(sn=+,sp=+) ne>0 =1 4(gV +g A)2(1 + cosθν) + 1 4(gV −g A)2(1−cosθ ν) M(sn=+,sp=−) ne>0 =g2 A(1 + cosθν) M(sn=−,sp=+) ne>0 =g2 A(1−cosθ ν) M(sn=−,sp=−) ne>0 =1 4(gV +g A)2(1−cosθ ν) + 1 4(gV −g A)...

    work page

  2. [2]

    (B9) The squared matrix element with equivalent normalization for neutral current scattering with neutrons is analogous

    Neutral current For neutral current interactions with protons, the reduced matrix element is Mss′ N Cp = 1 2 δss′ (1−4 sin 2 θW )2(1 + cosθν cosθ ′ ν + sinθ ν sinθ ′ ν cosϕ) +g 2 A(1 + cosθν cosθ ′ ν −sinθ ν sinθ ′ ν cosϕ)−2g As(1−4 sin 2 θW )(cosθ ν + cosθ ′ ν + (1−δ ss′ )g2 A(2−2 cosθ ν cosθ ′ ν) . (B9) The squared matrix element with equivalent normali...

    work page

  3. [3]

    Shorthand For convenience, here are the notation and conventions used throughout. gp = 5.5858 gn =−3.8263 sp, sn =±1 ˜k= k√ 2M T u= M T eB (1−e −eB/M T) t=e −eB/M T ∆ss′ n =− gn(s−s ′)eB 4M ∆ss′ p =− (gp −2)(s−s ′)eB 4M k0 = M q |q0 −∆ ss′ n | M= 2MnMp Mn +M p −U s δ(n) W M = 1 + 1.1kν M δ(p) W M = 1−7.1 kν M (C1) Us is the scalar potential for nucleons (...

    work page

  4. [4]

    Charged current with non-degenerate neutrons κνn = G2 F cos2 θceBρnδ(n) W M πcosh(g neB/4M T) X sn,sp nFD[µe −E + 0 ]Θ− sn,sp(−kν) ˜Vsn,sp E+2 0 2eB egnsneB/4M T × 1− ρpe(gp−2)speB/4M T cosh[gpeB/4M T] π M T 3/2 M T(1−t) eB+M T(1−t) cosh kzν E+ 0 2M T ×exp − k2 ⊥ν 2 1−t eB+M T(1−t) − k2 zν +E +2 0 4M T (C14) κ¯νp= G2 F cos2 θceBρpeeB/2M T δ(p) W M πcosh[g...

    work page

  5. [5]

    Charged current with degenerate neutrons κνn = G2 F cos2 θceBρnδ(n) W M π X sn,sp nFD[µe −E + 0 ]Θ− sn,sp(−kν) ˜Vsn,sp E+2 0 2eB × 1−(1−t) e(gp−2)speB/4M T 2 cosh[gpeB/4M T] yp 1−y p B r µ T + gnsneB 4M T , M T eB (1−t), egnsneB/8M T cosh(gneB/8M T) (C16) κ¯νp= G2 F cos2 θceBρpeeB/2M T δ(p) W M πcosh[g peB/4M T] X sn,sp nFD[−µe −E − 0 ][1−Θ + sn,sp(kν)] ˜...

    work page

  6. [6]

    Neutral current with non-degenerate neutrons ξn(q0, ⃗ q) = G2 F ρn 4 cosh[gneB/4M T] 1 q r 2πM T X ss′ Mss′ N CnegnseB/4M TΘ(q− |q 0 −∆ ss′ n |) × e−k2 0/2M T − ρnegns′eB/4M T √ 2 cosh[gneB/4M T] π2 M T qe−q2/4M T erf q+ 2k 0 2 √ M T + erf q−2k 0 2 √ M T (C18) κN Cn = X ss′ G2 F (kν + ∆ss′ n )2ρn 4πcosh[g neB/4M T] egnseB/4M T 1 2 Z Mss′ N Cnd cosθνν ′ − ...

    work page

  7. [7]

    Neutral current with degenerate neutrons ξn(q0, ⃗ q) =G2 F M2 8π2q X ss′ Mss′ N CnΘ(q− |q 0 −∆ ss′ n |) q0 +Tlog 1 +e β(ε0−q0−µ) −Tlog 1 +e β(ε0−µ) eβq0 −1 , (C23) 23 where the minimum incoming energy is ε0 = k2 0 2M − gneBs 4M .(C24) κN Cn = G2 F M (2π)3 r M T 2π X ss′ (kν + ∆ss′ n )2 1 2 Z Mss′ N Cnd cosθνν ′ ∆ss′ n eβ∆ss′ n −1 −T e −β(µn+∆ss′ n ) Θ(kν ...

    work page

  8. [8]

    Neutral current with protons ξp(q0, ⃗ q) = G2 F ρp 4 cosh[gpeB/4M T] eeB/2M T s 2πM T q2z X ss′ Mss′ N Cp e(gp−2)seB/4M T r qz q exp − M(q0 −∆ ss′ p )2 2q2z T + (1−cos 2 θq)I¯α(z) exp − q2 ⊥ 2eB 1 +t 1−t − (¯αeB−M|q0 −∆ ss′ p |)2 2q2z M T − ρp sinh(eB/2M T) cosh[gpeB/4M T] 2π3/2 eB √ M T e(gp−2)s′eB/2M T r qz q exp − M(q0 −∆ ss′ p )2 q2z T + (1−cos 2 θq)I...

    work page

  9. [9]

    Neutral current with electrons ξe(q0, ⃗ q) = X h G2 F eBT π e−q2 ⊥/2eBMhh N Ceδ[k ′ ν(h+ cosθ ′ ν)−k ν(h+ cosθ ν)]I(β, q0, µe),(C32) 24 where I(β, q0, µ) =    log 1 +e βµ −log 1 +e β(µ−q0) eβq0 −1 q0 >0 log 1 +e βµ+2q0 −log 1 +e β(µ+q0) eβq0 −1 q0 <0. (C33) Total neutral current opacity for electrons: κN Ce = G2 F (eB)2T 8π3 e−k2 ⊥/2eBX h=± [δh+4 s...

    work page

  10. [10]

    Urca reactions with non-degenerate neutrons ∂2Γ¯νn ∂kν ∂Ω = G2 F eBρnk2 ν 8π4 cosh(gneB/4M T) X sn,sp nFD[µe −E − 0 ]Θ− sn,sp(kν) ˜Vsn,sp E−2 0 2eB egnsneB/4M T × 1− ρpe(gp−2)speB/4M T cosh[gpeB/4M T] π M T 3/2 M T(1−t) eB+M T(1−t) cosh kzν E− 0 2M T ×exp − k2 ⊥ν 2 1−t eB+M T(1−t) − k2 zν +E −2 0 4M T (C38) ∂2Γνp ∂kν ∂Ω = G2 F eBρpk2 ν 8π4 cosh[gpeB/4M T]...

    work page

  11. [11]

    Urca reactions with degenerate neutrons ∂2Γ¯νn ∂kν ∂Ω = G2 F eBρnk2 ν 8π4 X sn,sp nFD[µe −E − 0 ]Θ− sn,sp(kν) ˜Vsn,sp E−2 0 2eB 1− e(gp−2)speB/4M T 2 cosh[gpeB/4M T] (1−t) × yp 1−y p B r µ T + gnsneB 4M T , M T eB (1−t), egnsneB/8M T cosh(gneB/8M T) (C40) 25 ∂2Γνp ∂kν ∂Ω = G2 F eBρpk2 ν 8π4 cosh[gpeB/4M T] X sn,sp nFD[E+ 0 −µ e]Θ− sn,sp(−kν) ˜Vsn,sp E+2 0...

    work page

  12. [12]

    Pianet al., Nature551, 67 (2017)

    E. Pianet al., Nature551, 67 (2017)

    work page 2017

  13. [13]

    J. J. Cowan, C. Sneden, J. E. Lawler, A. Aprahamian, M. Wiescher, K. Langanke, G. Mart´ ınez-Pinedo, and F.-K. Thielemann, Rev. Mod. Phys.93, 015002 (2021)

    work page 2021

  14. [14]

    Burrows, S

    A. Burrows, S. Reddy, and T. A. Thompson, Nuclear Physics A777, 356 (2006), special Isseu on Nuclear Astrophysics

    work page 2006

  15. [15]

    Kiuchi, P

    K. Kiuchi, P. Cerd´ a-Dur´ an, K. Kyutoku, Y. Sekiguchi, and M. Shibata, Phys. Rev. D92, 124034 (2015)

    work page 2015

  16. [16]

    D. J. Price and S. Rosswog, Science312, 719 (2006)

    work page 2006

  17. [17]

    Duan and Y.-Z

    H. Duan and Y.-Z. Qian, Phys. Rev. D69, 123004 (2004)

    work page 2004

  18. [18]

    Duan and Y.-Z

    H. Duan and Y.-Z. Qian, Phys. Rev. D72, 023005 (2005)

    work page 2005

  19. [19]

    Kumamoto and C

    M. Kumamoto and C. Welch, Phys. Rev. D111, 063009 (2025)

    work page 2025

  20. [20]

    L. B. Leinson and A. P´ erez, Journal of High Energy Physics1998, 020 (1998)

    work page 1998

  21. [21]

    D. A. Baiko and D. G. Yakovlev, Astronomy & Astrophysics342, 192 (1999), arXiv:astro-ph/9812071 [astro-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  22. [22]

    Maruyama, A

    T. Maruyama, A. B. Balantekin, M.-K. Cheoun, T. Kajino, M. Kusakabe, and G. J. Mathews, Physics Letters B824, 136813 (2022)

    work page 2022

  23. [23]

    Tambe, D

    P. Tambe, D. Chatterjee, M. Alford, and A. Haber, Phys. Rev. C111, 035809 (2025)

    work page 2025

  24. [24]

    Maruyama, A

    T. Maruyama, A. B. Balantekin, M.-K. Cheoun, T. Kajino, and G. J. Mathews, Physics Letters B805, 135413 (2020)

    work page 2020

  25. [25]

    Vogel, Phys

    P. Vogel, Phys. Rev. D29, 1918 (1984)

    work page 1918

  26. [26]

    C. J. Horowitz, Phys. Rev. D65, 043001 (2002)

    work page 2002

  27. [27]

    Szeg¨ o,Orthogonal Polynomials, 4th ed., American Mathematical Society Colloquium Publications, Vol

    G. Szeg¨ o,Orthogonal Polynomials, 4th ed., American Mathematical Society Colloquium Publications, Vol. 23 (American Mathematical Society, Providence, RI, 1975)

    work page 1975

  28. [28]

    Radice, F

    D. Radice, F. Galeazzi, J. Lippuner, L. F. Roberts, C. D. Ott, and L. Rezzolla, Monthly Notices of the Royal Astronomical Society460, 3255 (2016), https://academic.oup.com/mnras/article-pdf/460/3/3255/8130508/stw1227.pdf

    work page 2016

  29. [29]

    HEA-neutrinos,

    M. Kumamoto and C. Welch, “HEA-neutrinos,” (2025), github repository

    work page 2025

  30. [30]

    NS-landau-quantization,

    M. Kumamoto and C. Welch, “NS-landau-quantization,” (2024), github repository

    work page 2024

  31. [31]

    Canuto and H.-Y

    V. Canuto and H.-Y. Chiu, Phys. Rev.173, 1210 (1968)

    work page 1968

  32. [32]

    I. S. Gradshteyn and I. M. Ryzhik,Table of Integrals, Series, and Products(Academic Press, New York, 1980)

    work page 1980