Rotochemical heating of neutron stars in 2D
Abstract
The matter inside a spinning-down neutron star (NS) undergoes a slow compression that may drive it out of chemical equilibrium and, consequently, trigger non-equilibrium nuclear reactions that can heat up the star. This phenomenon is called rotochemical heating, and it is relatively well-investigated in the literature. Historically, it has been studied almost exclusively via 1D spherically symmetric models. However, rapidly rotating NSs are not spherically symmetric. We present pioneering simulations of rotochemical heating in full General Relativity in 2D axially symmetric models without relying on the slow rotation approximation. We explicitly follow thermal, mechanical, and chemical evolution of rotating NSs as they spin down taking into account both radial and angular distributions of temperature, composition, and mechanical structure. The simplest model is employed: only nucleonic composition, only the non-equilibrium modified Urca process as the main heating/cooling mechanism, and no superfluidity. We analyze the impact of the magnetic field strength as well as the effects of the heat blanketing envelope, initial rotation frequency, NS mass, and core EOS. Most of our conclusions coincide with the results of previous works on this topic, with one important exception. We have found a strong dependence of the cooling curves during the rotochemical heating phases on the core EOS, whereas in the literature it was regarded as very weak. The most likely explanation for this discrepancy is that the authors of the previous works have not considered EOSs with wide enough range of variation of the symmetry energy as a function of density.
Related
- A re-examination of neutron star cooling in transient sources. No shallow heating required: The role of thermonuclear heating leakage
- Structural quasi-universality in highly magnetized, differentially rotating neutron stars
- Torque lags from magnetic unlocking in neutron stars and gravitational-wave smearing
- Thermal evolution of neo-neutron stars. II. Temperature-dependent crusts
- Coexistence of strange quark stars and neutron stars: metastability and nucleation in proto-neutron stars
BibTeX
@misc{beznogov2026rotochemical,
title = {{Rotochemical heating of neutron stars in 2D}},
author = {Mikhail V. Beznogov and Jerome Novak and Dany Page and Adriana R. Raduta},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.20423},
eprint = {2610.20423},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.20423},
}
RIS
TY - UNPB
TI - Rotochemical heating of neutron stars in 2D
AU - Mikhail V. Beznogov
AU - Jerome Novak
AU - Dany Page
AU - Adriana R. Raduta
PY - 2026
DA - 2026/10/14
PB - arXiv
DO - 10.48550/arxiv.2610.20423
UR - https://arxiv.org/abs/2610.20423
AB - The matter inside a spinning-down neutron star (NS) undergoes a slow compression that may drive it out of chemical equilibrium and, consequently, trigger non-equilibrium nuclear reactions that can heat up the star. This phenomenon is called rotochemical heating, and it is relatively well-investigated in the literature. Historically, it has been studied almost exclusively via 1D spherically symmetric models. However, rapidly rotating NSs are not spherically symmetric. We present pioneering simulations of rotochemical heating in full General Relativity in 2D axially symmetric models without relying on the slow rotation approximation. We explicitly follow thermal, mechanical, and chemical evolution of rotating NSs as they spin down taking into account both radial and angular distributions of temperature, composition, and mechanical structure. The simplest model is employed: only nucleonic composition, only the non-equilibrium modified Urca process as the main heating/cooling mechanism, and no superfluidity. We analyze the impact of the magnetic field strength as well as the effects of the heat blanketing envelope, initial rotation frequency, NS mass, and core EOS. Most of our conclusions coincide with the results of previous works on this topic, with one important exception. We have found a strong dependence of the cooling curves during the rotochemical heating phases on the core EOS, whereas in the literature it was regarded as very weak. The most likely explanation for this discrepancy is that the authors of the previous works have not considered EOSs with wide enough range of variation of the symmetry energy as a function of density.
ER -