Data and code from: Critical slowing down and bulk viscosity in binary neutron star mergers
Data files
Jun 09, 2026 version files 474.34 KB
-
EVHRG.nb
29.26 KB
-
EVHRG.pdf
33.47 KB
-
PDG2021Plus_ThFIST_JMK_anti.dat
40.91 KB
-
README.md
1.25 KB
-
zeta_estimation.nb
201.58 KB
-
zeta_estimation.pdf
167.87 KB
Abstract
Hydrodynamic simulations of neutron star mergers rely on the clear separation between the strong-interaction, weak-interaction, and hydrodynamic timescales. In this effective framework, weak Urca interactions are typically the slowest microscopic processes, and therefore the Urca rate determines the bulk-viscous dissipation. This assumed hierarchy of dissipative mechanisms can be decisively altered, without invalidating hydrodynamics, if the trajectory of the matter in a neutron star merger passes through the vicinity of a possible low-temperature QCD critical point. The enhanced density fluctuations lead to critical slowing down and rapid growth of transport coefficients, including bulk viscosity. While this growth is regulated by finite-time effects, finite-size effects, and the breakdown of hydrodynamic scale separation, which bound the correlation length, we demonstrate that the QCD contribution to bulk viscosity can rival the electroweak contribution in realistic conditions. Thus, critical dynamics could leave observable imprints on the hydrodynamic evolution of neutron star mergers.
Dataset DOI: 10.5061/dryad.qv9s4mwvw
Description of the data and file structure
Run the notebook zeta_estimation.nb in Mathematica to reproduce the results of arXiv:2603.13474.
Files and variables
File: EVHRG.nb
Description: calculates EoS-dependent quantities using and excluded-volume HRG model.
File: EVHRG.pdf
Description: PDF version of EVHRG.nb.
File: zeta_estimation.nb
Description: main notebook for estimating the critical contribution to bulk viscosity, as described in arXiv:2603.13474.
File: zeta_estimation.pdf
Description: PDF version of zeta_estimation.nb.
File: PDG2021Plus_ThFIST_JMK_anti.dat
Description: list of resonances for HRG calculation from the PDG2021+ resonance list (see arXiv:2309.01737). Each column contains:
#1 ID
#2 Name
#3 Stable = 1/Unstable = 0
#4 Mass (in GeV)
#5 Degeneracy
#6 Statistics : Baryon = 1/Meson = -1
#7 Baryon number
#8 Electric charge
#9 Strangeness
#10 Charmness
#11 | S |
#12 | C |
#13 Width (in GeV)
#14 Threshold
Code/software
Notebooks were written in Wolfram 14.2 Mathematica.
