pith. machine review for the scientific record. sign in

arxiv: 2605.04814 · v1 · submitted 2026-05-06 · ⚛️ nucl-th · astro-ph.HE

Recognition: unknown

Charged current neutrino processes in hot nuclear matter with a recent Skyrme parametrization constrained by microscopic calculations

Michael Urban, Mingya Duan

Authors on Pith no claims yet

Pith reviewed 2026-05-08 16:41 UTC · model grok-4.3

classification ⚛️ nucl-th astro-ph.HE
keywords charged current neutrino processesSkyrme parametrizationhot nuclear matterproto-neutron starsrandom phase approximationdirect Urca processbeta equilibrium
0
0 comments X

The pith

Neutrino rates in hot nuclear matter with the Sky3s Skyrme parametrization differ by up to an order of magnitude from SLy4 results.

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

The paper computes charged current neutrino production, absorption, and direct Urca emission in hot nuclear matter using the random phase approximation. It compares the recent Sky3s Skyrme force, whose effective masses and spin terms were adjusted to microscopic calculations, against the SLy4 force from earlier studies. Rates obtained with Sky3s differ from SLy4 by up to one order of magnitude for some processes and energies. For densities above the direct Urca threshold, Sky3s yields stationary compositions where the standard beta equilibrium condition is less severely violated at finite temperature than literature predictions based on SLy4.

Core claim

Employing the Sky3s parametrization in RPA calculations of neutrino processes yields rates that differ from SLy4 by up to one order of magnitude. For densities above the direct Urca threshold, the electron, muon, and proton fractions leading to stationary composition show that the standard beta equilibrium condition is not as badly violated at finite temperature as previously found in the literature.

What carries the argument

The Sky3s extended Skyrme parametrization, with effective masses and spin-dependent terms constrained by microscopic calculations, used inside the random phase approximation to evaluate neutrino rates and stationary matter compositions.

Load-bearing premise

The Sky3s Skyrme parametrization, constrained only at zero temperature by microscopic calculations, remains accurate for the finite-temperature and finite-density regime of proto-neutron-star matter.

What would settle it

An independent finite-temperature calculation of the same neutrino rates and Urca thresholds using an ab initio many-body method, followed by direct numerical comparison to both the Sky3s and SLy4 results.

Figures

Figures reproduced from arXiv: 2605.04814 by Michael Urban, Mingya Duan.

Figure 1
Figure 1. Figure 1: FIG. 1. Response functions in HF, Landau approximation, and full RPA for Sky3s (top panels) at view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Similar to Fig. 1 but for view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Differential rates of (anti)neutrino production and absorption as a function of the (anti)neutrino energy view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Similar to Fig. 3 but for view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Response functions for Sky3s at view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Differential production rates of (anti)neutrinos per energy and volume as a function of the energy view at source ↗
read the original abstract

Neutrino processes are important in the modeling of supernova explosions, proto-neutron star evolution, and binary neutron star mergers. We study neutrino production and absorption in proto-neutron star and supernova matter and direct Urca neutrino emission of neutron star matter in the framework of the random phase approximation (RPA). As interactions, we employ the recent extended Skyrme parametrization Sky3s whose effective masses and spin-dependent terms were adjusted to microscopic calculations, and the SLy4 parametrization that was used in previous calculations of neutrino rates. The rates obtained for Sky3s differ from those for SLy4 by up to one order of magnitude for some processes and energy regions. We also determine the electron, muon, and proton fractions that lead to a stationary composition of matter for a density above the direct Urca threshold, and find that with Sky3s the standard $\beta$ equilibrium condition is not as badly violated at finite temperature as predicted in the literature. There are also minor differences between the full RPA and the common Landau approximation, but they are probably not significant for astrophysical simulations. We conclude that it would be worthwhile to repeat the calculation of neutrino rates for the use in astrophysical simulations, and the corresponding simulations, with several and better constrained interactions than SLy4, such as Sky3s.

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

3 major / 2 minor

Summary. The manuscript computes charged-current neutrino production, absorption, and direct Urca emission rates in hot, dense nuclear matter within the random-phase approximation (RPA), employing the Sky3s Skyrme parametrization (whose effective-mass and spin-dependent terms were fitted to zero-temperature microscopic calculations) and comparing results to the SLy4 parametrization used in prior work. It reports that Sky3s yields neutrino rates differing by up to an order of magnitude from SLy4 in selected processes and energy ranges, determines the electron/muon/proton fractions that produce stationary composition above the direct Urca threshold, and finds that the standard beta-equilibrium condition is less severely violated at finite temperature with Sky3s than claimed in the literature. Minor differences between full RPA and the Landau approximation are noted but deemed insignificant for simulations.

Significance. If the central results hold, the work demonstrates that neutrino rates in proto-neutron-star and supernova matter are sensitive to the choice of effective interaction, implying that astrophysical simulations relying on SLy4 may carry substantial systematic uncertainty. The explicit determination of stationary composition fractions above the direct Urca threshold provides concrete input for modeling. The use of a parametrization constrained by microscopic calculations at T=0 is a methodological strength that improves upon purely phenomenological fits, though the finite-temperature transferability remains untested within the manuscript.

major comments (3)
  1. [Methods / Skyrme parametrization description] The Sky3s effective masses and spin-orbit parameters are constrained exclusively by zero-temperature microscopic calculations (as stated in the introduction and methods). The RPA response functions and neutrino rates at finite temperature are computed by direct insertion of these parameters with no additional finite-T renormalization, validation against microscopic finite-T benchmarks, or sensitivity analysis; this assumption is load-bearing for the claimed order-of-magnitude rate differences and the reduced beta-equilibrium violation.
  2. [Results section] The abstract and results section assert rate differences 'up to one order of magnitude' for some processes and energy regions, yet the manuscript provides neither tabulated numerical values with uncertainties, explicit comparison plots with error bands, nor a quantitative breakdown by density, temperature, and composition; without these, the robustness of the central claim cannot be assessed.
  3. [Section on beta equilibrium and composition] The stationary-composition fractions (electron, muon, proton) above the direct Urca threshold are obtained by solving the same RPA framework; any systematic bias arising from the unverified finite-T transferability of Sky3s therefore propagates directly into the reported beta-equilibrium conclusions.
minor comments (2)
  1. [Figure captions] Figure captions and text should explicitly state the density, temperature, and composition ranges over which the order-of-magnitude differences are observed, rather than leaving the statement at the level of the abstract.
  2. [Results / RPA vs Landau comparison] The manuscript notes 'minor differences' between full RPA and the Landau approximation but does not quantify them (e.g., percentage deviation in rates); a short table or inset would clarify whether these differences are negligible for all relevant astrophysical conditions.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment point by point below, indicating planned revisions where appropriate.

read point-by-point responses

  1. Referee: The Sky3s effective masses and spin-orbit parameters are constrained exclusively by zero-temperature microscopic calculations (as stated in the introduction and methods). The RPA response functions and neutrino rates at finite temperature are computed by direct insertion of these parameters with no additional finite-T renormalization, validation against microscopic finite-T benchmarks, or sensitivity analysis; this assumption is load-bearing for the claimed order-of-magnitude rate differences and the reduced beta-equilibrium violation.

    Authors: We agree that Sky3s was fitted exclusively to zero-temperature microscopic calculations and that its finite-temperature application assumes transferability without additional renormalization or benchmarks in the present work. This is a standard practice for Skyrme parametrizations in the literature, and the comparison with SLy4 is intended to illustrate sensitivity to the choice of interaction. To address the concern, we will expand the Methods section with an explicit discussion of this assumption, its potential limitations for finite-T results, and the implications for the reported differences. We will also note the desirability of future finite-T microscopic validations. revision: partial

  2. Referee: The abstract and results section assert rate differences 'up to one order of magnitude' for some processes and energy regions, yet the manuscript provides neither tabulated numerical values with uncertainties, explicit comparison plots with error bands, nor a quantitative breakdown by density, temperature, and composition; without these, the robustness of the central claim cannot be assessed.

    Authors: The order-of-magnitude differences are visible in the comparison plots of the Results section for selected processes and kinematic regions. We acknowledge that a more quantitative presentation would improve assessability. In the revision we will add a summary table listing the largest relative differences (Sky3s vs. SLy4) for neutrino production, absorption, and direct Urca processes at representative densities, temperatures, and compositions, together with additional textual breakdown by conditions. Model uncertainties are systematic (arising from interaction choice) rather than statistical, and we will clarify this in the text. revision: yes

  3. Referee: The stationary-composition fractions (electron, muon, proton) above the direct Urca threshold are obtained by solving the same RPA framework; any systematic bias arising from the unverified finite-T transferability of Sky3s therefore propagates directly into the reported beta-equilibrium conclusions.

    Authors: The stationary compositions are obtained self-consistently within the same RPA framework used for the rates, so the finite-T transferability assumption applies equally. We will revise the relevant section to state explicitly that these fractions are specific to Sky3s, to cross-reference the expanded discussion of model assumptions added in response to the first comment, and to emphasize the model dependence of the reduced beta-equilibrium violation. revision: partial

Circularity Check

0 steps flagged

No circularity: rates and fractions derived forward from external parametrization

full rationale

The derivation computes neutrino rates and stationary-composition fractions directly from the RPA response functions using the Sky3s Skyrme force whose parameters were fixed by external zero-temperature microscopic calculations. No step renames a fit as a prediction, defines a quantity in terms of itself, or reduces the central claims to a self-citation chain. The comparison with SLy4 and with literature β-equilibrium results is an independent output of the same framework applied to two different interactions. The finite-temperature extrapolation of Sky3s is an assumption about model validity, not a circular reduction of the reported differences.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The results rest on the validity of the Skyrme effective interaction and the random-phase approximation for charged-current processes at finite temperature and density; Sky3s parameters are taken from prior microscopic fits rather than derived here.

free parameters (1)
  • Sky3s effective-mass and spin-orbit parameters
    Adjusted to match microscopic calculations; used directly in the RPA response functions.
axioms (2)
  • domain assumption Random phase approximation accurately captures the collective response of nuclear matter to weak probes
    Invoked for all rate calculations without further justification in the abstract.
  • domain assumption Skyrme force remains a good effective description at the temperatures and densities of proto-neutron-star matter
    Central modeling choice for both Sky3s and SLy4.

pith-pipeline@v0.9.0 · 5538 in / 1562 out tokens · 77709 ms · 2026-05-08T16:41:39.674654+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

53 extracted references

  1. [1]

    Notice that, in the charged-current case, we define τ± = 1 2(τ1 ±iτ 2),(26) whereτ i (i= 1,2,3) are the isospin Pauli matrices

    but that are needed in the case of generalized Skyrme functionals. Notice that, in the charged-current case, we define τ± = 1 2(τ1 ±iτ 2),(26) whereτ i (i= 1,2,3) are the isospin Pauli matrices. Then, besides theτ 3 contribution to⃗ ρ 2 1 as in Eq. (20) of Ref. [21], there are alsoτ + andτ − contributions, i.e., ⟨ ˆψ†τ+ ˆψ⟩⟨ ˆψ†τ− ˆψ⟩=⟨ ˆψ† n ˆψp⟩⟨ ˆψ† p ...

  2. [2]

    +v 0 3p1 ·p 2 + [v0 4(q) +v 0 5(p2 1 +p 2

  3. [3]

    However, as we will see below, the complete procedure is more complicated than in the neutral-current case because of the general- ized Lindhard functions

    +v 0 6p1 ·p 2]σ1 ·σ 2 +v 0 8iq·(p 1 −p 2)×(σ 1 +σ 2),(30) where the coefficientsv 0 i are now given by v0 1(q) = 4C ρ 1 + (−Cτ 1 −4C ∆ρ 1 + 4C∇ρ 1 )q2, v0 2 = 2C τ 1 , v0 3 =−4C τ 1 , v0 4(q) = 4C s 1 + (−CsT 1 −4C ∆s 1 + 4C∇⊗s 1 )q2, v0 5 = 2C sT 1 , v0 6 =−4C sT 1 , v0 8 = 2C ∇J 1 .(31) We can see that the expressions ofv 0 i here are simpler than in th...

  4. [4]

    S. W. Bruenn, Astrophs. J. Suppl58, 771 (1985)

    1985

  5. [5]

    Rampp and H.-T

    M. Rampp and H.-T. Janka, Astrophys. J. Lett.539, L33 (2000)

    2000

  6. [6]

    Mezzacappa, M

    A. Mezzacappa, M. Liebend¨ orfer, O. E. Messer, W. R. Hix, F.-K. Thielemann, and S. W. Bruenn, Phys. Rev. Lett.86, 1935 (2001)

    1935

  7. [7]

    Rampp and H.-T

    M. Rampp and H.-T. Janka, Astron. Astrophys.396, 361 (2002)

    2002

  8. [8]

    T. A. Thompson, A. Burrows, and P. A. Pinto, Astro- phys. J.592, 434 (2003)

    2003

  9. [9]

    Buras, M

    R. Buras, M. Rampp, H.-T. Janka, and K. Kifonidis, Astron. Astrophys.447, 1049 (2006)

    2006

  10. [10]

    C. J. Horowitz and A. Schwenk, Phys. Lett. B642, 326 (2006)

    2006

  11. [11]

    L. F. Roberts and S. Reddy, Phys. Rev. C95, 045807 (2017)

    2017

  12. [13]

    Rrapaj, J

    E. Rrapaj, J. W. Holt, A. Bartl, S. Reddy, and A. Schwenk, Phys. Rev. C91, 035806 (2015)

    2015

  13. [14]

    Vida˜ na, D

    I. Vida˜ na, D. Logoteta, and I. Bombaci, Phys. Rev. C 106, 035804 (2022)

    2022

  14. [15]

    E. Shin, E. Rrapaj, J. W. Holt, and S. K. Reddy, Phys. Rev. C109, 015804 (2024)

    2024

  15. [16]

    Reddy, M

    S. Reddy, M. Prakash, and J. M. Lattimer, Phys. Rev. D58, 013009 (1998)

    1998

  16. [17]

    Reddy, M

    S. Reddy, M. Prakash, J. M. Lattimer, and J. A. Pons, Phys. Rev. C59, 2888 (1999)

    1999

  17. [18]

    A. A. Dzhioev and G. Mart´ ınez-Pinedo, inEuropean Physical Journal Web of Conferences, European Phys- ical Journal Web of Conferences, Vol. 194 (EDP, 2018) p. 02006

    2018

  18. [19]

    Oertel, A

    M. Oertel, A. Pascal, M. Mancini, and J. Novak, Phys. Rev. C102, 035802 (2020)

    2020

  19. [20]

    Pascal, J

    A. Pascal, J. Novak, and M. Oertel, Mon. Not. R. Astron. Soc.511, 356 (2022)

    2022

  20. [21]

    Chabanat, P

    E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A635, 231 (1998)

    1998

  21. [22]

    Duan and M

    M. Duan and M. Urban, Phys. Rev. C108, 025813 (2023)

    2023

  22. [23]

    Duan and M

    M. Duan and M. Urban, Phys. Rev. C112, 045808 (2025)

    2025

  23. [24]

    Duan and M

    M. Duan and M. Urban, Phys. Rev. C110, 065806 (2024)

    2024

  24. [25]

    Baldo, G

    M. Baldo, G. F. Burgio, H. J. Schulze, and G. Taranto, Phys. Rev. C89, 048801 (2014)

    2014

  25. [26]

    W. Zuo, C. Shen, and U. Lombardo, Phys. Rev. C67, 037301 (2003)

    2003

  26. [27]

    W. Zuo, U. Lombardo, and C.-W. Shen, Communica- tions in Theoretical Physics39, 219 (2003)

    2003

  27. [28]

    W. Zuo, U. Lombardo, C.-w. Shen, J.-Y. Liu, and J.- Q. Li, High Energy Physics and Nuclear Physics27, 999 (2003)

    2003

  28. [29]

    Vida˜ na, A

    I. Vida˜ na, A. Polls, and V. Durant, Phys. Rev. C94, 054006 (2016)

    2016

  29. [30]

    Bigdeli, G

    M. Bigdeli, G. H. Bordbar, and Z. Rezaei, Phys. Rev. C 80, 034310 (2009)

    2009

  30. [31]

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

    2024

  31. [32]

    Janka, inphysics of Neutron Stars of Supernovae, edited by A

    H.-T. Janka, inphysics of Neutron Stars of Supernovae, edited by A. W. Alsabti and P. Murdin (Springer, Cham,

  32. [33]

    Fischer, G

    T. Fischer, G. Guo, A. A. Dzhioev, G. Mart´ ınez-Pinedo, M.-R. Wu, A. Lohs, and Y.-Z. Qian, Phys. Rev. C101, 025804 (2020)

    2020

  33. [34]

    Iwamoto and C

    N. Iwamoto and C. J. Pethick, Phys. Rev. D25, 313 (1982)

    1982

  34. [35]

    J. D. Bjorken and S. D. Drell,Relativistic quantum me- chanics(McGraw-Hill, New York, 1964)

    1964

  35. [36]

    J. J. Sakurai and J. Napolitano,Modern Quantum Me- chanics(Cambridge University Press, Cambridge, 2017)

    2017

  36. [37]

    Vautherin and D

    D. Vautherin and D. M. Brink, Phys. Rev. C5, 626 (1972)

    1972

  37. [38]

    Beiner, H

    M. Beiner, H. Flocard, N. Van Giai, and P. Quentin, Nucl. Phys. A238, 29 (1975)

    1975

  38. [39]

    H. S. K¨ ohler, Nucl. Phys. A258, 301 (1976)

    1976

  39. [40]

    Bender, P.-H

    M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys75, 121 (2003)

    2003

  40. [41]

    Davesne, A

    D. Davesne, A. Pastore, and J. Navarro, Phys. Rev. C 100, 064301 (2019)

    2019

  41. [42]

    Urban and S

    M. Urban and S. Ramanan, Phys. Rev. C101, 035803 (2020)

    2020

  42. [43]

    E. S. Hern´ andez, J. Navarro, and A. Polls, Nucl. Phys. A658, 327 (1999)

    1999

  43. [44]

    Pastore, D

    A. Pastore, D. Davesne, and J. Navarro, Phys. Rep.563, 1 (2015)

    2015

  44. [45]

    Margueron,Effet du milieu sur la propagation des neu- trinos dans la mati` ere nucl´ eaire, Thesis, Universit´ e Paris Sud - Paris XI (2001)

    J. Margueron,Effet du milieu sur la propagation des neu- trinos dans la mati` ere nucl´ eaire, Thesis, Universit´ e Paris Sud - Paris XI (2001)

    2001

  45. [46]

    Burrows and R

    A. Burrows and R. F. Sawyer, Phys. Rev. C59, 510 (1999). 17

    1999

  46. [47]

    Haensel, A

    P. Haensel, A. Y. Potekhin, and D. G. Yakovlev,Neutron Stars 1: Equation of State and Structure, Astrophysics and Space Science Library, Vol. 326 (Springer, New York, 2007)

    2007

  47. [48]

    Degenaar and V

    N. Degenaar and V. F. Suleimanov, inThe Physics and Astrophysics of Neutron Stars, Astrophysics and Space Science Library, Vol. 457, edited by L. Rezzolla, P. Piz- zochero, D. I. Jones, N. Rea, and I. Vida˜ na (Springer, New York, 2018) p. 185

    2018

  48. [49]

    D. G. Yakovlev and C. J. Pethick, Annu. Rev. Astron. Astrophys.42, 169 (2004)

    2004

  49. [50]

    J. M. Lattimer and D. F. Swesty, Nucl. Phys. A535, 331 (1991)

    1991

  50. [51]

    Chabanat, P

    E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A627, 710 (1997)

    1997

  51. [52]

    M. G. Alford and S. P. Harris, Phys. Rev. C98, 065806 (2018)

    2018

  52. [53]

    Sedrakian, Phys

    A. Sedrakian, Phys. Rev. Lett.133, 171401 (2024)

    2024

  53. [54]

    Flowers, M

    E. Flowers, M. Ruderman, and P. Sutherland, Astro- phys. J.205, 541 (1976)

    1976