pith. sign in

arxiv: 1907.05921 · v1 · pith:BEWVDWBKnew · submitted 2019-07-12 · ✦ hep-ph · astro-ph.HE· nucl-th

Matter And Gravitation In Collisions of heavy ions and neutron stars: equation of state

Pith reviewed 2026-05-24 22:13 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.HEnucl-th
keywords equation of stateneutron star mergersgravitational wavesheavy ion collisionsQCD matterphase structurequark matter
0
0 comments X

The pith

Gravitational wave signals from neutron star mergers can be combined with heavy ion collision data to determine the equation of state of dense QCD matter.

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

The paper shows that the densities and temperatures reached in neutron star mergers match those in heavy ion collisions at existing labs to within 20 percent. A single equation of state for strongly interacting QCD matter can therefore be constructed and applied to both systems despite their different rapidity windows and impact parameters. Gravitational waves from future LIGO and Virgo events then become a direct messenger for the phase structure at high baryon density. This messenger can be cross-checked against multiplicity fluctuations and flow measurements recorded in laboratory detectors. The approach turns two separate experimental domains into complementary constraints on the same underlying physics.

Core claim

The gravitational waves emitted from a binary neutron star merger are sensitive to the appearance of quark matter and the stiffness of the equation of state of QCD matter. These astrophysical extremes match to within 20 percent the densities and temperatures reached in relativistic hydrodynamics of heavy ion collisions. One unified equation of state can therefore be constructed and used for both neutron star physics and hot QCD matter created at laboratory facilities, allowing gravitational wave signals to be combined with heavy ion data to pin down the equation of state and phase structure of dense matter.

What carries the argument

A single unified equation of state for QCD matter that simultaneously describes the inner cores of neutron stars and the matter produced in heavy ion collisions.

If this is right

  • Gravitational wave signals will directly probe the stiffness of the equation of state and the possible appearance of quark matter.
  • High-multiplicity fluctuations and flow measurements in heavy ion detectors will supply independent constraints on the same equation of state.
  • The phase structure of dense QCD matter becomes accessible through two orthogonal experimental routes.
  • Future advanced LIGO and Virgo events can be interpreted with laboratory data already in hand.

Where Pith is reading between the lines

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

  • Discrepancies between the two data sets could indicate that the 20 percent matching window is not sufficient for a single equation of state.
  • The method opens the possibility of using one domain to calibrate uncertainties in the other before new observations arrive.
  • Detector design for both gravitational wave observatories and heavy ion experiments could be guided by the need for matching thermodynamic coverage.

Load-bearing premise

The underlying QCD physics remains similar enough in the two systems that one equation of state applies to both.

What would settle it

A gravitational wave signal from a neutron star merger whose extracted equation of state at a given density is incompatible with the equation of state obtained from heavy ion collision observables at the same density.

Figures

Figures reproduced from arXiv: 1907.05921 by Anton Motornenko, Horst Stoecker, Jan Steinheimer, Stefan Schramm, Volodymyr Vovchenko.

Figure 1
Figure 1. Figure 1: Effective masses of octet baryons (solid) and their respective parity partners (dashed) as functions of temperature T obtained from the CMF model at µB = 0 for isospin symmetric matter, Q/B = 1 2 . With the increase of the temperature the degeneracy in mass between baryons and respective parity partners arises. The Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model [8, 9], CMF, … view at source ↗
Figure 2
Figure 2. Figure 2: The CMF energy per baryon E/B − mN as function of baryon density nB normalized to nuclear saturation density n0 = 0.16 fm−3 , blue – for symmetric nuclear matter ,Q/B = 1 2 , red – for nuclear matter in β-equilibrium, total electric charge is zero Q = 0 that is allowed by presence of leptons. 100 125 150 175 200 225 250 T (MeV) 0.0 0.5 1.0 1.5 2.0 2.5 3.0 P = T 4 CMF model lattice, WB lattice, HotQCD 100 1… view at source ↗
Figure 3
Figure 3. Figure 3: Pressure P I (right) and trace anomaly I (left) at µB = 0 as function of temperature T. Comparison between the CMF model predictions and LQCD results [6, 22]. posite parity. To provide the dynamical mass generation, the octet baryons and their partners are coupled to scalar chiral fields σ and ζ , non-strange and strange chiral condensates, respectively. The σ and ζ fields serve as order parameters for the… view at source ↗
Figure 4
Figure 4. Figure 4: Particle density ratios to the density of baryons ni/nB at T = 0, for quarks a factor of 1/3 is used, presented as functions of baryon density nB. The CMF-results are obtained for isospin symmetric matter (left) and in β-equilibrium (right). of the baryons, where ns is the number of strange quarks in the baryon, and ms = 130 MeV is the mass of the strange quark. The couplings g (j) i are tuned to reproduce… view at source ↗
Figure 5
Figure 5. Figure 5: CMF baryon number kurtosis, ratio of fourth and second order susceptibilities, χ B 4 /χ B 2 in the baryon chemical potential and temperature, µB − T, plane. The three distinctive critical regions and their remnants spread from T = 0 up to T > 200 MeV. Numerous phenomenological features of the CMF model suggest a rather rich phase structure of the model. The critical phenomena associated with the nuclear li… view at source ↗
Figure 6
Figure 6. Figure 6: Colormaps for chiral condensate σ and quark fraction 1 3 nq/nB along the µB −T plane for isospin symmetric matter. Note that significant quark fraction appears only after chiral symmetry is restored, so pure quark matter arises at high µB or T. Mind the transparency of white color, so two quantities are presented on the same plot. MeV, so the change in baryon density is smooth and not step-like. Note, even… view at source ↗
Figure 7
Figure 7. Figure 7: Evolution for different collision energy √ sNN of system excited in heavy-ion collisions along the µB − T (left) and T − nB (right) phase diagram. Initial state as implicit function of √ sNN is calculated by Taub adiabat and is presented by the black line. Colored lines – isentropic lines of constant entropy per baryon S/A at different bombarding energies √ sNN respectively. Heavy ion fixed target experime… view at source ↗
Figure 8
Figure 8. Figure 8: Mass-radius relation (left) and mass - central density relation (right) for NSs calculated using the CMF equation of state. Yellow colorbands present constraint on the maximum mass of NS obtained in [55], red colorband – the mass of NS PSR J0348+0432 [56]. are present as well. The temperatures in the star interiors are negligibly small compared to hot QCD scales, so the T = 0 EoS can be applied. The result… view at source ↗
Figure 9
Figure 9. Figure 9: The CMF tidal deformability Λ (left) and second Love number k2 (right) as functions of NS mass. λ is directly proportional to the second Love number k2: λ = 2 3 k2R 5 . (5.2) For convenience, usually the dimensionless tidal deformability Λ is presented as: Λ = λ M5 = 2 3 k2  R M 5 . (5.3) Here, M and R are the mass and radius of the NS. We present the resulting dimensionless tidal deformability coefficie… view at source ↗
Figure 10
Figure 10. Figure 10: Regions of nB−T phase diagram reachable in lower energy heavy ion collisions and in NS merg￾ers. Bold line present results from full 3D hydrodynamic simulation with CMF EoS used. Thin lines indicate isentropes where entropy per baryon, S/A was calculated by RRHT adiabat. Red and blue ellipses indicate regions of maximal temperature and density respectively, reachable in NS mergers, extracted from [50]. on… view at source ↗
read the original abstract

The gravitational waves emitted from a binary neutron star merger, as predicted from general relativistic magneto-hydrodynamics calculations, are sensitive to the appearance of quark matter and the stiffness of the equation of state of QCD matter present in the inner cores of the stars. This is a new messenger observable from outer space, which does provide direct signals for the phase structure of strongly interacting QCD matter at high baryon density and high temperature. These astrophysically created extremes of thermodynamics do match, to within 20\%, the values of densities and temperatures which we find in relativistic hydrodynamics and transport theory of heavy ion collisions at the existing laboratories, if though at quite different rapidity windows, impact parameters and bombarding energies of the heavy nuclear systems. We demonstrate how one unified equation of state can be constructed and used for both neutron star physics and hot QCD matter excited at laboratory facilities. The similarity in underlying QCD physics allows the gravitational wave signals from future advanced LIGO and Virgo events to be combined with the analysis of high multiplicity fluctuations and flow measurements in heavy ion detectors in the lab to pin down the EoS and the phase structure of dense matter.

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

1 major / 0 minor

Summary. The manuscript argues that densities and temperatures in binary neutron star mergers and relativistic heavy-ion collisions match to within 20%, permitting construction of a single unified QCD equation of state applicable to both systems; gravitational-wave signals from LIGO/Virgo can then be combined with heavy-ion flow and fluctuation data to constrain the EoS and phase structure.

Significance. If the unification and cross-system applicability were demonstrated, the approach would offer a concrete route to jointly constrain the high-density QCD EoS using complementary observables, which is of clear interest for mapping the phase diagram of dense matter.

major comments (1)
  1. [Abstract] Abstract: the central claim that 'similarity in underlying QCD physics' allows one EoS to be used for both neutron-star mergers and heavy-ion collisions rests on an untested assumption; the text supplies no derivation, no explicit EoS construction, and no quantitative check that differences in isospin asymmetry (beta-equilibrated neutron-rich matter versus near-symmetric HIC) or strangeness content leave the pressure-density relation unchanged within the stated 20% tolerance.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed reading and the constructive comment on the abstract. We respond to the major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that 'similarity in underlying QCD physics' allows one EoS to be used for both neutron-star mergers and heavy-ion collisions rests on an untested assumption; the text supplies no derivation, no explicit EoS construction, and no quantitative check that differences in isospin asymmetry (beta-equilibrated neutron-rich matter versus near-symmetric HIC) or strangeness content leave the pressure-density relation unchanged within the stated 20% tolerance.

    Authors: The manuscript's central point is that the thermodynamic conditions (baryon density and temperature) reached in neutron-star mergers and in heavy-ion collisions overlap to within 20%, as obtained from GRMHD simulations and from relativistic hydrodynamics/transport calculations, respectively. On this basis we advocate constructing and employing a single QCD equation of state for both systems. We agree that the present text contains neither an explicit functional form for such an EoS nor a quantitative assessment of how isospin asymmetry or net strangeness modify the pressure-density relation inside the quoted 20% window. The paper is framed as a perspective advocating the joint use of gravitational-wave and heavy-ion observables rather than a technical derivation of the EoS itself. We will therefore revise the abstract to state more precisely that the unification rests on the overlap of thermodynamic variables and will add a brief paragraph noting that future work must quantify the residual effects of asymmetry and strangeness. revision: partial

Circularity Check

0 steps flagged

No significant circularity; central claim rests on external similarity assumption without self-referential reduction

full rationale

The provided text asserts that thermodynamic conditions match to within 20% and that 'similarity in underlying QCD physics' permits construction of one unified EoS for both neutron-star mergers and heavy-ion collisions. No equations, fitted parameters, or self-citations are exhibited that reduce this unification claim to an input by construction (e.g., no parameter fitted to one domain is then relabeled as a prediction for the other). The similarity statement is an assumption offered as justification rather than a derived result that loops back to itself. The paper therefore contains no load-bearing self-definitional, fitted-input, or self-citation steps of the enumerated kinds; the derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract only, so the ledger is necessarily incomplete; the central proposal rests on an assumed overlap of thermodynamic conditions and on the existence of a single EoS that can be constrained by both data sets.

axioms (2)
  • standard math General relativity and relativistic hydrodynamics accurately describe neutron star mergers and heavy ion collisions
    Invoked implicitly when stating that GW signals are sensitive to the EoS and that hydrodynamics/transport theory applies to HICs.
  • domain assumption QCD matter at high baryon density obeys a single equation of state that is independent of the production mechanism
    The claim that one unified EoS can be used for both systems rests on this assumption.

pith-pipeline@v0.9.0 · 5748 in / 1332 out tokens · 19117 ms · 2026-05-24T22:13:22.077647+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

60 extracted references · 60 canonical work pages · 39 internal anchors

  1. [1]

    LIGO S CIENTIFIC , VIRGO collaboration. Phys. Rev. Lett. 116 (2016) 061102 [1602.03837]

  2. [2]

    E. R. Most, L. J. Papenfort, V . Dexheimer, M. Hanauske, S. Schramm, H. Stoecker et al. Phys. Rev. Lett. 122 (2019) 061101 [1807.03684]

  3. [3]

    Identifying a first-order phase transition in neutron star mergers through gravitational waves

    A. Bauswein, N.-U. F. Bastian, D. B. Blaschke, K. Chatziioannou, J. A. Clark, T. Fischer et al. Phys. Rev. Lett. 122 (2019) 061102 [1809.01116]

  4. [4]

    E. V . Shuryak. Phys. Rept. 61 (1980) 71

  5. [5]

    Phases of Dense Quarks at Large N_c

    L. McLerran and R. D. Pisarski. Nucl. Phys. A796 (2007) 83 [0706.2191]

  6. [6]

    Full result for the QCD equation of state with 2+1 flavors

    S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg and K. K. Szabo. Phys. Lett. B730 (2014) 99 [1309.5258]

  7. [7]

    H OTQCD collaboration. Phys. Rev. D96 (2017) 074510 [1708.04897]

  8. [8]

    An effective chiral Hadron-Quark Equation of State

    J. Steinheimer, S. Schramm and H. Stoecker. J. Phys. G38 (2011) 035001 [1009.5239]

  9. [9]

    Motornenko, J

    A. Motornenko, J. Steinheimer, V . V ovchenko, S. Schramm and H. Stoecker. .1905.00866

  10. [10]

    C. E. Detar and T. Kunihiro. Phys. Rev. D39 (1989) 2805

  11. [11]

    Hatsuda and M

    T. Hatsuda and M. Prakash. Phys. Lett. B224 (1989) 11

  12. [12]

    Phase Transition in the chiral $\sigma$-$\omega$ model with dilatons

    P. Papazoglou, J. Schaffner, S. Schramm, D. Zschiesche, H. Stoecker and W. Greiner. Phys. Rev. C55 (1997) 1499 [nucl-th/9609035]

  13. [13]

    Chiral Lagrangian for strange hadronic matter

    P. Papazoglou, S. Schramm, J. Schaffner-Bielich, H. Stoecker and W. Greiner. Phys. Rev. C57 (1998) 2576 [nucl-th/9706024]

  14. [14]

    Nuclei in a chiral SU(3) model

    P. Papazoglou, D. Zschiesche, S. Schramm, J. Schaffner-Bielich, H. Stoecker and W. Greiner. Phys. Rev. C59 (1999) 411 [nucl-th/9806087]

  15. [15]

    Thermodynamics of dense hadronic matter in a parity doublet model

    C. Sasaki and I. Mishustin. Phys. Rev. C82 (2010) 035204 [1005.4811]

  16. [16]

    The hadronic SU(3) Parity Doublet Model for Dense Matter, its extension to quarks and the strange equation of state

    J. Steinheimer, S. Schramm and H. Stoecker. Phys. Rev. C84 (2011) 045208 [1108.2596]

  17. [17]

    Hybrid Stars in an SU(3) Parity Doublet Model

    V . Dexheimer, J. Steinheimer, R. Negreiros and S. Schramm. Phys. Rev. C87 (2013) 015804 [1206.3086]. 13 Matter And Gravitation In Collisions Anton Motornenko

  18. [18]

    Higher-order baryon number susceptibilities: interplay between the chiral and the nuclear liquid-gas transitions

    A. Mukherjee, J. Steinheimer and S. Schramm. Phys. Rev. C96 (2017) 025205 [1611.10144]

  19. [19]

    QCD at high density: Equation of state for nuclear collisions and neutron stars

    A. Motornenko, V . V ovchenko, J. Steinheimer, S. Schramm and H. Stoecker.Nucl. Phys. A982 (2019) 891 [1809.02000]

  20. [20]

    Light baryons below and above the deconfinement transition: medium effects and parity doubling

    G. Aarts, C. Allton, D. De Boni, S. Hands, B. Jäger, C. Praki et al. JHEP 06 (2017) 034 [1703.09246]

  21. [21]

    Hyperons in thermal QCD: A lattice view

    G. Aarts, C. Allton, D. De Boni and B. Jäger. . 1812.07393

  22. [22]

    H OTQCD collaboration. Phys. Rev. D90 (2014) 094503 [1407.6387]

  23. [23]

    Chiral effective model with the Polyakov loop

    K. Fukushima. Phys. Lett. B591 (2004) 277 [hep-ph/0310121]

  24. [24]

    The problem of repulsive quark interactions - Lattice versus mean field models

    J. Steinheimer and S. Schramm. Phys. Lett. B696 (2011) 257 [1005.1176]

  25. [25]

    Do lattice data constrain the vector interaction strength of QCD?

    J. Steinheimer and S. Schramm. Phys. Lett. B736 (2014) 241 [1401.4051]

  26. [26]

    Phases of QCD: lattice thermodynamics and a field theoretical model

    C. Ratti, M. A. Thaler and W. Weise. Phys. Rev. D73 (2006) 014019 [hep-ph/0506234]

  27. [27]

    M. A. Stephanov. Phys. Rev. Lett. 102 (2009) 032301 [0809.3450]

  28. [28]

    Quarkyonic Matter and Neutron Stars

    L. McLerran and S. Reddy. Phys. Rev. Lett. 122 (2019) 122701 [1811.12503]

  29. [29]

    Stoecker and W

    H. Stoecker and W. Greiner. Phys. Rept. 137 (1986) 277

  30. [30]

    L. P. Csernai, J. Kapusta and L. D. McLerran. Phys. Rev. Lett. 97 (2006) 152303 [nucl-th/0604032]

  31. [31]

    Viscosity Information from Relativistic Nuclear Collisions: How Perfect is the Fluid Observed at RHIC?

    P. Romatschke and U. Romatschke. Phys. Rev. Lett. 99 (2007) 172301 [0706.1522]

  32. [32]

    (3+1)-Dimensional Hydrodynamic Expansion with a Critical Point from Realistic Initial Conditions

    J. Steinheimer, M. Bleicher, H. Petersen, S. Schramm, H. Stoecker and D. Zschiesche. Phys. Rev. C77 (2008) 034901 [0710.0332]

  33. [33]

    A. H. Taub. Phys. Rev. 74 (1948) 328

  34. [34]

    Stoecker, W

    H. Stoecker, W. Greiner and W. Scheid. Z. Phys. A286 (1978) 121

  35. [35]

    K. S. Thorne. Astrophys. J. 179 (1973) 897

  36. [36]

    H. G. Baumgardt, J. U. Schott, Y . Sakamoto, E. Schopper, H. Stoecker, J. Hofmann et al. Z. Phys. A273 (1975) 359

  37. [37]

    Stoecker, G

    H. Stoecker, G. Graebner, J. A. Maruhn and W. Greiner. Phys. Lett. 95B (1980) 192

  38. [38]

    Stoecker, A

    H. Stoecker, A. A. Ogloblin and W. Greiner. Z. Phys. A303 (1981) 259

  39. [39]

    Stoecker, M

    H. Stoecker, M. Gyulassy and J. Boguta. Phys. Lett. 103B (1981) 269

  40. [40]

    Hahn and H

    D. Hahn and H. Stoecker. Nucl. Phys. A476 (1988) 718

  41. [41]

    A. V . Merdeev, L. M. Satarov and I. N. Mishustin. Phys. Rev. C84 (2011) 014907 [1103.3988]

  42. [42]

    Hanauske, J

    M. Hanauske, J. Steinheimer, L. Bovard, A. Mukherjee, S. Schramm, K. Takami et al. J. Phys. Conf. Ser .878 (2017) 012031

  43. [43]

    Hadron resonance gas and mean-field nuclear matter for baryon number fluctuations

    K. Fukushima. Phys. Rev. C91 (2015) 044910 [1409.0698]

  44. [44]

    van der Waals Interactions in Hadron Resonance Gas: From Nuclear Matter to Lattice QCD

    V . V ovchenko, M. I. Gorenstein and H. Stoecker.Phys. Rev. Lett. 118 (2017) 182301 [1609.03975]

  45. [45]

    Critical point of nuclear matter and beam energy dependence of net proton number fluctuations

    V . V ovchenko, L. Jiang, M. I. Gorenstein and H. Stoecker.Phys. Rev. C98 (2018) 024910 [1711.07260]. 14 Matter And Gravitation In Collisions Anton Motornenko

  46. [46]

    Nuclear interactions and net-proton number fluctuations in heavy ion collisions at the SIS18 accelerator

    J. Steinheimer, Y . Wang, A. Mukherjee, Y . Ye, C. Guo, Q. Li et al.Phys. Lett. B785 (2018) 40 [1804.08936]

  47. [47]

    Y . Ye, Y . Wang, J. Steinheimer, Y . Nara, H.-j. Xu, P. Li et al.Phys. Rev. C98 (2018) 054620 [1808.06342]

  48. [48]

    Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

    T. Dietrich, S. Bernuzzi, M. Ujevic and B. Brügmann. Phys. Rev. D91 (2015) 124041 [1504.01266]

  49. [49]

    Probing Extreme-Density Matter with Gravitational Wave Observations of Binary Neutron Star Merger Remnants

    D. Radice, S. Bernuzzi, W. Del Pozzo, L. F. Roberts and C. D. Ott. Astrophys. J. 842 (2017) L10 [1612.06429]

  50. [50]

    Hanauske, J

    M. Hanauske, J. Steinheimer, A. Motornenko, V . V ovchenko, L. Bovard, E. R. Most et al.Particles 2 (2019) 44

  51. [51]

    Hanauske, L

    M. Hanauske, L. Bovard, E. Most, J. Papenfort, J. Steinheimer, A. Motornenko et al. Universe 5 (2019)

  52. [52]

    R. C. Tolman. Phys. Rev. 55 (1939) 364

  53. [53]

    J. R. Oppenheimer and G. M. V olkoff. Phys. Rev. 55 (1939) 374

  54. [54]

    G. Baym, C. Pethick and P. Sutherland. Astrophys. J. 170 (1971) 299

  55. [55]

    Using gravitational-wave observations and quasi-universal relations to constrain the maximum mass of neutron stars

    L. Rezzolla, E. R. Most and L. R. Weih. Astrophys. J. 852 (2018) L25 [1711.00314]

  56. [56]

    R. S. Lynch, J. Boyles, S. M. Ransom, I. H. Stairs, D. R. Lorimer, M. A. McLaughlin et al. Astrophys. J. 763 (2013) 81 [1209.4296]

  57. [57]

    A new quark-hadron hybrid equation of state for astrophysics - I. High-mass twin compact stars

    S. Benic, D. Blaschke, D. E. Alvarez-Castillo, T. Fischer and S. Typel. Astron. Astrophys. 577 (2015) A40 [1411.2856]

  58. [58]

    Y . Song, G. Baym, T. Hatsuda and T. Kojo. . 1905.01005

  59. [59]

    Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral

    T. Hinderer, B. D. Lackey, R. N. Lang and J. S. Read. Phys. Rev. D81 (2010) 123016 [0911.3535]

  60. [60]

    Hanauske and L

    M. Hanauske and L. Bovard. Journal of Astrophysics and Astronomy 39 (2018) 45. 15