arxiv: 2604.28062 · v1 · submitted 2026-04-30 · ✦ hep-ph · astro-ph.HE

Recognition: unknown

Semi-analytic bounds on axion-like-particle supernovae emission

Ana Luisa Foguel , Eduardo S. Fraga

Authors on Pith no claims yet

Pith reviewed 2026-05-07 06:04 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HE

keywords axion-like particlessupernova coolingSN1987Aproto-neutron starsemi-analytic boundsALP mass and couplingnuclear uncertainties

0 comments

The pith

A semi-analytic model using six proto-neutron star parameters yields axion-like particle bounds from supernovae that match numerical simulations.

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

The paper develops a semi-analytic method to constrain axion-like particles produced in core-collapse supernovae by extending a framework that parameterizes the proto-neutron star with six global parameters. It incorporates finite ALP mass into the calculation of production and escape, then derives limits in the coupling-mass plane from the SN1987A cooling bound. This approach is faster than numerical simulations and the resulting bounds agree well with prior simulation results, showing that the method captures the essential physics despite astrophysical uncertainties. A sympathetic reader would care because it offers a simpler way to explore how different modeling choices affect new physics constraints from supernovae.

Core claim

By expressing proto-neutron star observables in terms of six global parameters and including finite mass effects for axion-like particles, the authors obtain bounds on the axion-nucleon coupling versus mass from the requirement that ALP emission does not over-cool the supernova. These bounds are shown to be in good agreement with those from full numerical simulations, demonstrating the robustness of the reduced description.

What carries the argument

The semi-analytic framework that expresses proto-neutron star observables through six global parameters and calculates finite-mass ALP emission rates to derive cooling constraints.

Load-bearing premise

The six global proto-neutron star parameters together with the selected nuclear effects and cooling exclusion criteria are sufficient to describe the dominant processes of ALP production and escape during the supernova cooling phase.

What would settle it

A full numerical supernova simulation that incorporates the exact same six-parameter calibration but yields ALP bounds differing by more than the reported agreement margin would falsify the robustness claim.

Figures

Figures reproduced from arXiv: 2604.28062 by Ana Luisa Foguel, Eduardo S. Fraga.

**Figure 1.** Figure 1: FIG. 1. Upper panel: ALP nucleon-nucleon Bremsstrahlung view at source ↗

**Figure 2.** Figure 2: FIG. 2. ALP (solid) and neutrino (dashed) luminosities at view at source ↗

**Figure 3.** Figure 3: FIG. 3. Supernova cooling bounds in the axion-nucleon cou view at source ↗

**Figure 5.** Figure 5: FIG. 5. Supernova-cooling bounds in the ( view at source ↗

read the original abstract

Core-collapse supernovae provide natural laboratories for the production of new light particles. In particular, axion-like particles (ALPs) can be constrained via SN1987A cooling arguments. However, significant astrophysical and nuclear uncertainties imply that such bounds may vary strongly depending on modeling choices, even when expensive simulations are employed. In this context, semi-analytic methods offer a simple and fast alternative for deriving new-physics constraints. Building on a previous semi-analytic framework, in which proto-neutron star (PNS) observables are expressed in terms of six global PNS parameters, we include a finite ALP mass in the calculation and derive bounds in the axion-nucleon coupling versus mass plane. The obtained bounds are in good agreement with previous results from numerical simulations, demonstrating the robustness of the method. We also illustrate the sensitivity of the bounds to different PNS parameter calibrations, nuclear effects and cooling exclusion criteria.

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.

Desk Editor's Note private letter to a colleague

This extends semi-analytic supernova bounds to finite-mass ALPs with results close to simulations, but the six global PNS parameters remain the key assumption to watch.

read the letter

This paper gives a workable semi-analytic route to bounds on finite-mass axion-like particles from supernovae, and the results track the numerical simulations reasonably well. They take an existing six-parameter description of the proto-neutron star and extend it to include ALP mass in the production and escape calculations. That produces new limits in the g_aN - m_a plane. The abstract notes good agreement with prior simulation-based bounds, and they check how the limits change when you vary the PNS parameters, nuclear effects, and the cooling cutoff. The strength is the speed and the explicit sensitivity study. If you need to scan a lot of ALP models without running full supernova codes every time, this looks practical. The fact that they recover similar bounds to the heavy simulations is the main selling point. The potential issue is whether six global parameters really capture enough of the radial structure. For m_a > 0 the optical depth and the bremsstrahlung rate vary with local T and rho, so any mismatch in the assumed profiles could affect the integrated luminosity at the level of the quoted uncertainties. They do vary the global params and show the effect, but a direct test against a full numerical density-temperature profile for the massive case would strengthen the robustness claim. Without seeing the full equations I can't tell how they handled the mass dependence in the phase space integrals, but the abstract suggests they did it. This work is aimed at the particle astrophysics community that sets limits on light new physics. It is incremental but cleanly executed, so it should go to peer review. Ask the referees to look at the profile modeling for the finite-mass case and whether the agreement holds when the PNS structure is taken from a real simulation rather than the parameterized one.

Referee Report

2 major / 2 minor

Summary. The manuscript extends a prior semi-analytic framework for proto-neutron star (PNS) observables, now incorporating finite axion-like particle (ALP) mass into nucleon bremsstrahlung production and escape probability calculations. Using a PNS model reduced to six global parameters (mass, radius, core and surface temperatures, density scale, etc.), it derives bounds on the ALP-nucleon coupling g_aN versus ALP mass m_a from SN1987A cooling arguments. The central result is that these bounds agree with those obtained from full numerical simulations, with additional exploration of sensitivity to PNS parameter choices, nuclear effects, and cooling exclusion criteria.

Significance. If the agreement with numerical simulations holds under detailed scrutiny, the work provides a computationally efficient and transparent alternative to expensive simulations for placing ALP constraints from supernovae. This enables rapid sensitivity studies to astrophysical and nuclear uncertainties, which is valuable in a field where such modeling choices can significantly affect bounds. The approach builds on an established six-parameter reduction and demonstrates its utility for finite-mass cases, potentially aiding future explorations of new-physics scenarios in core-collapse supernovae.

major comments (2)

[Results section] Results section: The central claim that the semi-analytic bounds are in good agreement with previous numerical simulations (and thus demonstrate robustness) requires stronger substantiation for m_a > 0. The manuscript reports sensitivity to the six global PNS parameters but does not present a direct side-by-side comparison of ALP luminosity or derived g_aN bounds computed from the analytic profiles versus integration over full numerical radial T(r) and ρ(r) profiles for any finite-mass point. Since both the production rate and the optical depth τ(r) for massive ALPs depend on local conditions in the decoupling region, profile deviations at the 10-20% level could shift the bounds comparably to the quoted agreement.
[Method section] Method section: The reduction of the PNS to six global parameters (with assumed analytic profiles such as isothermal core plus exponential atmosphere) is load-bearing for the finite-mass calculation. While the paper illustrates sensitivity to parameter calibrations, it should quantify the fidelity of these profiles against specific simulation outputs in the regions controlling ALP trapping and emission, e.g., by reporting the fractional difference in integrated luminosity for benchmark m_a values.

minor comments (2)

[Abstract] The abstract refers to 'six global PNS parameters' without enumerating them; adding a short list (e.g., mass, radius, core temperature, surface temperature, density scale, and one additional) would improve immediate clarity for readers.
[Methods] Ensure that the cooling exclusion criteria are defined with an explicit equation or threshold value in the methods, as their impact on the bounds is highlighted in the sensitivity analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which help clarify the presentation of our semi-analytic approach. We address each major comment below and describe the revisions we will make to strengthen the substantiation of our results.

read point-by-point responses

Referee: [Results section] Results section: The central claim that the semi-analytic bounds are in good agreement with previous numerical simulations (and thus demonstrate robustness) requires stronger substantiation for m_a > 0. The manuscript reports sensitivity to the six global PNS parameters but does not present a direct side-by-side comparison of ALP luminosity or derived g_aN bounds computed from the analytic profiles versus integration over full numerical radial T(r) and ρ(r) profiles for any finite-mass point. Since both the production rate and the optical depth τ(r) for massive ALPs depend on local conditions in the decoupling region, profile deviations at the 10-20% level could shift the bounds comparably to the quoted agreement.

Authors: We agree that a direct side-by-side comparison of the ALP luminosity (and resulting optical depth) for finite m_a would provide stronger substantiation of the agreement with numerical simulations. The current manuscript demonstrates agreement at the level of the final derived bounds on g_aN versus m_a (see Figure 5 and associated discussion), which already incorporate the integrated production and escape probabilities. However, to address the referee's concern explicitly, we will add a new panel or appendix figure in the revised manuscript that shows the fractional difference in computed ALP luminosity for benchmark finite-mass values (e.g., m_a = 10 MeV and m_a = 100 MeV) when evaluated with the six-parameter analytic profiles versus direct integration over the radial T(r) and ρ(r) profiles from the reference numerical simulation. This will quantify any profile-induced discrepancies in the decoupling region and confirm that they remain smaller than the overall uncertainty band on the bounds. revision: yes
Referee: [Method section] Method section: The reduction of the PNS to six global parameters (with assumed analytic profiles such as isothermal core plus exponential atmosphere) is load-bearing for the finite-mass calculation. While the paper illustrates sensitivity to parameter calibrations, it should quantify the fidelity of these profiles against specific simulation outputs in the regions controlling ALP trapping and emission, e.g., by reporting the fractional difference in integrated luminosity for benchmark m_a values.

Authors: We acknowledge that the six-parameter reduction with analytic profiles is central to the finite-mass extension, and its fidelity in the trapping and emission regions must be quantified for full transparency. The parameters are calibrated to reproduce global PNS observables (total energy loss, neutrino luminosities, and average temperatures) from the reference numerical simulations, as established in our prior work. To directly address the request, we will include in the revised Methods section (or a dedicated subsection) the fractional difference in integrated ALP luminosity for the same benchmark m_a values, computed using the analytic profiles versus the specific simulation radial profiles in the density range 10^13–10^14 g/cm^3 and temperature range 10–30 MeV that dominate the decoupling. This will be presented alongside the existing sensitivity studies to PNS parameters. revision: yes

Circularity Check

0 steps flagged

Finite-mass ALP bounds derived from prior PNS parameterization but validated externally against simulations

full rationale

The derivation extends an existing semi-analytic PNS model (six global parameters) by incorporating finite ALP mass into production and escape integrals, then computes new g_aN bounds. These bounds are compared to independent numerical simulations for validation rather than being forced by the inputs. The prior framework supplies the base profiles and calibration, but the mass-dependent optical depth and luminosity integrals introduce new content not reducible to the calibration by construction. No self-citation chain or fitted-input renaming carries the central result; the agreement with simulations serves as an external check on the approximation.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the six global PNS parameters whose values are calibrated to observations or simulations, plus standard assumptions about ALP production mechanisms in hot dense matter and the validity of the cooling argument for SN1987A.

free parameters (1)

six global PNS parameters
Used to express all relevant proto-neutron star observables; their specific values are chosen or fitted to match known supernova properties.

axioms (2)

domain assumption ALP emission can be treated as an additional cooling channel whose rate depends on the axion-nucleon coupling and ALP mass.
Invoked when extending the semi-analytic cooling model to include finite-mass ALPs.
domain assumption The six-parameter description adequately captures the essential thermodynamics and neutrino/ALP transport in the PNS.
Core modeling choice stated in the abstract as the basis for the semi-analytic framework.

pith-pipeline@v0.9.0 · 5455 in / 1517 out tokens · 58484 ms · 2026-05-07T06:04:11.124291+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

43 extracted references · 25 canonical work pages

[1]

G. G. Raffelt,Stars as laboratories for fundamental physics: The astrophysics of neutrinos, axions, and other weakly interacting particles(1996)

1996
[2]

Hirataet al.(Kamiokande-II), Phys

K. Hirataet al.(Kamiokande-II), Phys. Rev. Lett.58, 1490 (1987)

1987
[3]

R. M. Biontaet al., Phys. Rev. Lett.58, 1494 (1987)

1987
[4]

R. D. Peccei and H. R. Quinn, Phys. Rev. Lett.38, 1440 (1977)

1977
[5]

Weinberg, Phys

S. Weinberg, Phys. Rev. Lett.40, 223 (1978)

1978
[6]

Wilczek, Phys

F. Wilczek, Phys. Rev. Lett.40, 279 (1978)

1978
[7]

Chikashige, R

Y. Chikashige, R. N. Mohapatra, and R. D. Peccei, Phys. Lett. B98, 265 (1981)

1981
[8]

G. B. Gelmini and M. Roncadelli, Phys. Lett. B99, 411 (1981)

1981
[9]

R. T. Co and K. Harigaya, Phys. Rev. Lett.124, 111602 (2020), arXiv:1910.02080 [hep-ph]

work page arXiv 2020
[10]

Bittar, G

P. Bittar, G. Burdman, and G. M. Salla, Phys. Rev. D 112, 055039 (2025), arXiv:2410.00964 [hep-ph]

work page arXiv 2025
[11]

Wilczek, Phys

F. Wilczek, Phys. Rev. Lett.49, 1549 (1982)

1982
[12]

Sikivie, Phys

P. Sikivie, Phys. Rev. Lett.51, 1415 (1983)

1983
[13]

Di Luzio, M

L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli, Phys. Rept.870, 1 (2020), arXiv:2003.01100 [hep-ph]

work page arXiv 2020
[14]

C. A. J. O’Hare, PoSCOSMICWISPers, 040 (2024), arXiv:2403.17697 [hep-ph]

work page arXiv 2024
[15]

Bharucha, F

A. Bharucha, F. Br¨ ummer, N. Desai, and S. Mutzel, JHEP02, 141 (2023), arXiv:2209.03932 [hep-ph]

work page arXiv 2023
[16]

Armando, P

G. Armando, P. Panci, J. Weiss, and R. Ziegler, Phys. Rev. D109, 055029 (2024), arXiv:2310.05827 [hep-ph]

work page arXiv 2024
[17]

Beenakker, D

W. Beenakker, D. Mikkers, A. V. Phan, and S. Westhoff, (2025), arXiv:2511.19619 [hep-ph]

work page arXiv 2025
[18]

Brandt, F

B. Brandt, F. Cuteri, G. Endr˝ odi, J. J. H. Hern´ andez, and G. Mark´ o, PoSLA TTICE2022, 174 (2023), arXiv:2212.03385 [hep-lat]

work page arXiv 2023
[19]

B. B. Brandt, G. Endr˝ odi, J. J. Hern´ andez Hern´ andez, G. Mark´ o, and L. Pannullo, (2026), arXiv:2603.29153 [hep-lat]

work page arXiv 2026
[20]

M. S. Turner, Phys. Rev. Lett.60, 1797 (1988)

1988
[21]

Carenza, T

P. Carenza, T. Fischer, M. Giannotti, G. Guo, G. Mart´ ınez-Pinedo, and A. Mirizzi, JCAP10, 016 (2019), [Erratum: JCAP 05, E01 (2020)], arXiv:1906.11844 [hep-ph]

work page arXiv 2019
[22]

Lella, P

A. Lella, P. Carenza, G. Lucente, M. Giannotti, and A. Mirizzi, Phys. Rev. D107, 103017 (2023), arXiv:2211.13760 [hep-ph]

work page arXiv 2023
[23]

Lella, P

A. Lella, P. Carenza, G. Co’, G. Lucente, M. Giannotti, A. Mirizzi, and T. Rauscher, Phys. Rev. D109, 023001 (2024), arXiv:2306.01048 [hep-ph]

work page arXiv 2024
[24]

Payez, C

A. Payez, C. Evoli, T. Fischer, M. Giannotti, A. Mirizzi, and A. Ringwald, JCAP02, 006 (2015), arXiv:1410.3747 [astro-ph.HE]

work page arXiv 2015
[25]

J. S. Lee, (2018), arXiv:1808.10136 [hep-ph]

work page arXiv 2018
[26]

Lucente, P

G. Lucente, P. Carenza, T. Fischer, M. Giannotti, and A. Mirizzi, JCAP12, 008 (2020), arXiv:2008.04918 [hep- ph]

work page arXiv 2020
[27]

A. L. Foguel and E. S. Fraga, Astropart. Phys.151, 102855 (2023), arXiv:2209.14318 [hep-ph]

work page arXiv 2023
[28]

Y. Suwa, A. Harada, K. Nakazato, and K. Sumiyoshi, PTEP2021, 013E01 (2021), arXiv:2008.07070 [astro- ph.HE]

work page arXiv 2021
[29]

Weinberg,Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity(John Wiley and Sons, New York, 1972)

S. Weinberg,Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity(John Wiley and Sons, New York, 1972)

1972
[30]

Raffelt and L

G. Raffelt and L. Stodolsky, Phys. Rev. D37, 1237 (1988)

1988
[31]

Chakraborty, A

S. Chakraborty, A. Gupta, and M. Vanvlasselaer, Phys. Rev. D110, 063032 (2024), arXiv:2403.12169 [hep-ph]

work page arXiv 2024
[32]

Carena and R

M. Carena and R. D. Peccei, Phys. Rev. D40, 652 (1989)

1989
[33]

Hadronic axion window and the big-bang nucleosynthesis

S. Chang and K. Choi, Phys. Lett. B316, 51 (1993), arXiv:hep-ph/9306216

work page Pith review arXiv 1993
[34]

R. P. Brinkmann and M. S. Turner, Phys. Rev. D38, 2338 (1988)

1988
[35]

M. S. Turner, H. S. Kang, and G. Steigman,Axions, SN 1987A, and one pion exchange, Tech. Rep. (Chicago Univ., IL (USA), 1988)

1988
[36]

Suwa, Publ

Y. Suwa, Publ. Astron. Soc. Jap.66, L1 (2014), arXiv:1311.7249 [astro-ph.HE]

work page arXiv 2014
[37]

Pengyu, Data for ”achieving optimal-distance atom-loss correction via pauli envelope”, 10.5281/zen- odo.19339056 (2026)

H. Dembinski and P. O. et al., (2020), 10.5281/zen- odo.3949207. 9

work page doi:10.5281/zen- 2020
[38]

Fischer, S

T. Fischer, S. Chakraborty, M. Giannotti, A. Mirizzi, A. Payez, and A. Ringwald, Phys. Rev. D94, 085012 (2016), arXiv:1605.08780 [astro-ph.HE]

work page arXiv 2016
[39]

Burrows, M

A. Burrows, M. T. Ressell, and M. S. Turner, Phys. Rev. D42, 3297 (1990)

1990
[40]

Raffelt and D

G. Raffelt and D. Seckel, Phys. Rev. Lett.60, 1793 (1988)

1988
[41]

G. G. Raffelt, Lect. Notes Phys.741, 51 (2008), arXiv:hep-ph/0611350

work page arXiv 2008
[42]

Betranhandy and E

A. Betranhandy and E. O’Connor, Phys. Rev. D106, 063019 (2022), arXiv:2204.00503 [astro-ph.HE]

work page arXiv 2022
[43]

Alonso-Gonz´ alez, D

D. Alonso-Gonz´ alez, D. Cerde˜ no, M. Cerme˜ no, and A. D. Perez, Phys. Rev. D111, 083029 (2025), arXiv:2412.19890 [hep-ph]

work page arXiv 2025