Next Article in Journal
Photon Acceleration by Superluminal Ionization Fronts
Next Article in Special Issue
Nuclear Matter and Finite Nuclei: Recent Studies Based on Parity Doublet Model
Previous Article in Journal
Numerical Analysis of the Transfer Dynamics of Heavy Metals from Soil to Plant and Application to Contamination of Honey
Previous Article in Special Issue
On the Anomalous Dimension in QCD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constraints on Phase Transitions in Neutron Star Matter

Physics Department, TUM School of Natural Sciences, Technical University of Munich, 85747 Garching, Germany
*
Authors to whom correspondence should be addressed.
Symmetry 2024, 16(1), 111; https://doi.org/10.3390/sym16010111
Submission received: 21 December 2023 / Revised: 11 January 2024 / Accepted: 15 January 2024 / Published: 18 January 2024

Abstract

:
Recent inference results of the sound velocity in the cores of neutron stars are summarized. Implications for the equation of state and the phase structure of highly compressed baryonic matter are discussed. In view of the strong constraints imposed by the heaviest known pulsars, the equation of state must be very stiff in order to ensure the stability of these extreme objects. This required stiffness limits the possible appearance of phase transitions in neutron star cores. For example, a Bayes factor analysis quantifies strong evidence for squared sound velocities c s 2 > 0.1 in the cores of 2.1 solar-mass and lighter neutron stars. Only weak first-order phase transitions with a small phase coexistence density range Δ ρ / ρ < 0.2 (at the 68% level) in a Maxwell construction still turn out to be possible within neutron stars. The central baryon densities in even the heaviest neutron stars do not exceed five times the density of normal nuclear matter. In view of these data-based constraints, much discussed issues such as the quest for a phase transition towards restored chiral symmetry and the active degrees of freedom in cold and dense baryonic matter, are reexamined.

1. Introduction

Within recent decades, the observational data base for neutron stars has expanded substantially. The detailed analysis of this data has sharpened the empirical constraints on the equation of state (EoS) of matter deep inside the cores of these extreme objects [1]. In particular, the existence of heavy neutron stars with masses around and above two solar masses implies that the EoS, i.e., pressure P ( ε ) as a function of energy density ε , must be sufficiently stiff in order to support such massive compact stars against gravitational collapse. Some previously discussed simple forms of exotic matter could, therefore, be excluded when their corresponding EoSs turned out to be too soft.
The quest for a possible phase transition in cold and dense baryonic matter—from hadronic to quark degrees of freedom—has been a topic of prime interest for a long time [2,3]. In this context, a summary will be given of the present state of knowledge about the neutron star EoS as inferred directly from the existing observational data. Of particular interest is the speed of sound in the cores of neutron stars. Its behavior as a function of energy density is a sensitive indicator for phase transitions or continuous changes (crossovers) of degrees of freedom in a star’s composition.
Primary sources of information for the inference of the EoS are neutron star masses deduced from Shapiro delay measurements [4,5,6,7,8,9], combined determinations of masses and radii inferred from X-ray data detected with the NICER telescope [10,11,12,13,14], and gravitational wave signals from binary neutron star mergers observed by the LIGO and Virgo Collaborations [15,16,17]. These data have served as basic input for a large variety of Bayesian and other inference analyses [18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33] in search of a range of data-compatible equations of state with controlled uncertainties. Low-density constraints from nuclear physics around the equilibrium density of nuclear matter, ρ 0 = 0.16 fm 3 , are commonly introduced by referring to results of chiral effective field theory (ChEFT) [34]. At asymptotically high densities, perturbative QCD (pQCD) is applicable and provides a constraint to be matched in extrapolations far beyond the conditions realized in neutron star cores [35,36]. Bayesian inference studies have recently been further extended [31] by the supplementary observational data reported in [37] (the black widow pulsar PSR J0952–0607 with a mass of M = 2.35 ± 0.17 M in solar mass units).
This report presents a survey of results obtained in the detailed inference analysis of Ref. [31] and discusses possible interpretations. A central topic is the evaluation of constraints and possible evidence for a first-order phase transition within the posterior uncertainties of the EoS deduced from the existing empirical neutron star data. Later parts of this manuscript refer to ongoing considerations about the structure and composition of dense matter in neutron star cores. In view of the constraints emerging from the analysis of the observational data, especially of the heaviest known neutron stars, we discuss the quest for a chiral phase transition, quark–hadron continuity and its realization as a soft crossover, a Fermi liquid approach to neutron star matter, and other related issues.

2. Equation of State of Neutron Star Matter

2.1. Observational Constraints

2.1.1. Neutron Star Masses and Radii

Information about the masses of the heaviest neutron stars derives primarily from precise Shapiro time delay measurements of pulsars orbiting in binary systems with white dwarf companions. Three such massive objects, PSR J1614–2230 [4,6,7], PSR J0348+0432 [5], and PSR J0740+6620 [8,9], have been established in the past:
PSR J1614–2230 M = 1.908 ± 0.016 M ,
PSR J0348+0432 M = 2.01 ± 0.04 M ,
PSR J0740+6620 M = 2.08 ± 0.07 M .
The heaviest neutron star observed so far was recently reported [37]—a black widow pulsar,
PSR J0952–0607 M = 2.35 ± 0.17 M .
This result was obtained based on the detection of the black widow’s companion star. PSR J0952–0607 is also one of the fastest rotating pulsars. With a spin period of 1.4 ms, it requires radius-dependent corrections for rotational effects, as described in [31]. For example, the rotating mass of a 2.35 M star with an assumed radius of R = 12 km decreases by 3% to an equivalent non-rotating mass of about 2.28 M .
Together with their masses the radii of neutron stars can be inferred from X-ray profiles of rotating hot-spot patterns measured with the NICER telescope. Two neutron stars have been investigated in this way [10,12]:
PSR J0030+0451 M = 1.34 0.16 + 0.15 M , R = 12.17 1.19 + 1.14 km ,
PSR J0740+6620 M = 2.072 0.066 + 0.067 M , R = 12.39 0.98 + 1.30 km .
All of the mass and radius data are included as inputs when setting constraints for the neutron star EoS. In the case of PSR J0740+6620, the analysis of the NICER measurement includes the Shapiro delay result. Therefore, only the data from NICER have been used in the Bayesian inference dataset. An alternative analysis of the NICER data was also performed by a second independent group [11,13]. Both results are mutually compatible within their uncertainties.

2.1.2. Binary Neutron Star Mergers and Tidal Deformabilities

Gravitational wave signals produced by a merger of two neutron stars in a binary have been detected by the LIGO and Virgo collaborations [15,16,17]. These signals are interpreted using theoretical waveform models that depend on the mass ratio of the two neutron stars, M 2 / M 1 , and a mass-weighted combination of their tidal deformabilities. Based on information from the GW170817 event, the (dimensionless) tidal deformability for a 1.4 solar mass neutron star was deduced [16]:
Λ 1.4 = 190 120 + 390 .
The GW170817 event was further investigated together with electromagnetic signals. In one such analysis [38], the following masses and tidal deformabilities were reported for the individual neutron stars in the binary:
M 1 = 1.46 0.09 + 0.13 M , Λ 1 = 255 171 + 416 , M 2 = 1.26 0.12 + 0.09 M , Λ 2 = 661 375 + 858 .
The GW190425 event, which could have also been a binary neutron star merger, was also included in the inference analysis [17].

2.2. Inference of Sound Velocity and the EoS in Neutron Stars

Based on the data in the previous subsection, the squared sound velocity and the equation of state
c s 2 ( ε ) = P ( ε ) ε and P ( ε ) = 0 ε d ε c s 2 ( ε ) ,
of neutron star matter can be inferred using Bayesian methods. In our recent work [29,31], the ChEFT constraint at low baryon density was used in the Bayesian inference procedure as a likelihood (not as a prior) within a conservative window of baryon densities [39], ρ 1.3 ρ 0 . The equations of state were interpolated to pQCD results at asymptotically high densities. Unlike the conclusions drawn in [36], and as explained in detail in [31], we did, however, not find a strongly constraining influence of pQCD extrapolations down to neutron star densities. Possible reasons for this discrepancy were also discussed in [40].
The result of the inferred sound velocity in Figure 1 shows, within the 68 % highest posterior density credible bands, a rapid increase in c s 2 beyond the conformal limit of 1/3 at energy densities relevant for a wide range of neutron stars with masses M 1.4–2.3 M . A quantification in terms of a corresponding Bayes factor [21,29,31] demonstrates extreme evidence for c s 2 exceeding the conformal bound inside neutron stars.
Figure 2 presents the inferred posterior bands for the EoS. Notably, the additional inclusion of the heaviest pulsar PSR J0952–0607 in the dataset leads to a more rapidly rising pressure as compared to previous inference results, which do not incorporate these data. Within the range of energy densities 0.5–0.8 GeV/fm 3 reached in the core of a 2.3 solar mass star and corresponding to baryon densities between 3.0 ρ 0 and 4.6 ρ 0 , the pressure even exceeds that of the time-honored APR EoS [41].
A further quantity of interest is the chemical potential associated with the conserved baryon number,
μ ( ρ ) = ε ( ρ ) ρ = P + ε ρ ,
where the latter equality is the Gibbs–Duhem relation at a temperature of zero. By solving the equation ρ = ρ l exp ε l ε ( ρ ) d ε [ ε + P ( ε ) ] 1 for ε ( ρ ) with a boundary condition ε l = ε ( ρ l ) chosen at a suitably low density ρ l , the baryon chemical potential can be computed for any given equation of state P ( ε ) within the uncertainty bands in Figure 2. The resulting μ ( ρ ) bands in Figure 3 show a steep rise at baryon densities of ρ > 2 ρ 0 , indicating that strongly repulsive forces are at work. The agnostic approach for the inference of the sound speed on which these results are based does not permit distinction between different species of constituents inside neutron stars, so μ represents the total chemical potential from all active degrees of freedom carrying baryon number: μ = i x i μ i , where x i = ρ i / ρ is the fraction corresponding to each baryonic species i. However, irrespective of whether these species are nucleons, other baryonic composites, or quarks, the empirically inferred behavior of μ demonstrates that these degrees of freedom must be correlated by strongly repulsive forces as the density increases in order to build up the necessary high pressure in the cores of neutron stars to support two solar masses and beyond. We return to this discussion again at a later stage.
The approach to conformality in strongly interacting matter at high baryon density has been an important theme in several recent investigations [42,43,44,45]. A key quantity is the trace anomaly measure,
Δ = g μ ν T μ ν 3 ε = 1 3 P ( ε ) ε .
Conformal matter is characterized by an equation of state P = ε / 3 , so the signature for its occurrence is Δ 0 over an extended energy density range. The posterior credible bands for the trace anomaly measure deduced from the data-based EoS are shown in Figure 4. Starting from Δ = 1 / 3 at zero energy density, the median of Δ turns negative at ε 0.7 GeV / fm 3 , entering a high-pressure domain with P > ε / 3 . This crossing appears within the range of energy densities possibly realized in the cores of the heaviest neutron stars. At much higher energy densities beyond those displayed in Figure 3, pQCD then implies a switch back to positive Δ before the asymptotic limit Δ 0 is approached.

2.3. Selected Neutron Star Properties

Numerically solving the Tolman–Oppenheimer–Volkov (TOV) equations leads to posterior credible bands for the mass–radius relation of neutron stars, as displayed in Figure 5. Notably, the median of R ( M ) suggests an almost constant radius R, which is independent of M.
The density profiles ρ ( r ) computed for neutron stars in the mass range 1.4–2.3 M reach central baryon densities, ρ c = ρ ( r = 0 ) , that are systematically below 5 ρ 0 in the cores of even the heaviest stars. For example,
ρ c ( 1.4 M ) = ( 2.6 ± 0.4 ) ρ 0 and ρ c ( 2.3 M ) = ( 3.8 ± 0.8 ) ρ 0 ,
at the 68% level [31]. In a baryonic picture of neutron stars, this implies that the average distance between baryons even in the highly compressed star center exceeds 1 fm, more than twice the characteristic short-range hard-core distance of 1/2 fm in nucleon–nucleon interactions. We return to this point at a later stage.
For each given EoS, the TOV equations can be solved in combination with equations for the tidal deformability Λ . This leads to posterior probability bands, Λ ( M ) , as a function of neutron star mass, as shown in Figure 6. The large uncertainties of the empirical values deduced from the analysis of the gravitational wave events do not allow one to constrain the tidal deformability bands to a stronger degree so far.

3. Constraints on Phase Transitions in Neutron Stars

The nature and location of the transition from dense nuclear or neutron matter to quark and gluon degrees of freedom is still largely unknown. Lattice computations have been successful in exploring the QCD phase structure at high temperatures and vanishing baryon chemical potential, establishing that the transition is a continuous crossover at a temperature of T c 155 MeV. However, at nonzero baryon densities, the notorious sign problem in the Euclidean action prevents a similarly systematic approach. As a consequence, studies of the phase diagram at low temperatures and high baryon densities are primarily based on models with varying assumptions about the relevant active degrees of freedom [47].
A key element of this discussion is the question about a possible phase transition from spontaneously broken to restored chiral symmetry in QCD. At high baryon densities beyond ρ > 2 ρ 0 , various hypotheses have been explored, ranging from a chiral first-order phase transition to continuous hadron–quark crossovers [2,48,49,50,51,52]. The behavior of the sound speed c s is a prime indicator for phase transitions or crossovers. This is schematically illustrated in Figure 7, which shows typical patterns of such transitions as they would show up in the sound speed. For example, a first-order phase transition with Maxwell construction is characterized by a region of constant pressure over an interval of density (or energy density) in which two phases coexist. Between the lower and upper endpoints of this phase coexistence interval, the squared sound velocity c s 2 = P / ε jumps to zero and back. In a first-order phase transition with Gibbs construction, the pressure in the mixed phase is not constant, but P / ε still changes discontinuously [53]. A crossover has no such discontinuities but features a pronounced maximum in c s 2 . Purely nucleonic scenarios lead instead to a continuously rising sound speed.

3.1. Evidence against a Very Low Squared Sound Speed in Neutron Stars

One indication of the possible appearance of a first-order phase transition would be a very low squared sound velocity ( c s 2 0.1 ) within the range of energy densities relevant for neutron stars. Let us introduce c s , m i n as the minimum speed of sound following a maximum at lower densities. Given the inferred bands of probability distributions for c s 2 in Section 2.2 and Figure 1, one can quantify the evidence for two competing scenarios: minimum c s , m i n 2 > 0.1 versus c s , m i n 2 0.1 ; in terms of corresponding Bayes factors, B c s , m i n 2 0.1 c s , m i n 2 > 0.1 . A detailed analysis of such Bayes factors was carried out in [31] for the minimum speed of sound up to different maximum masses of neutron stars, M m a x . Bayes factors exceeding 10 or 100 indicate strong or extreme evidence, respectively, for the scenario where the minimum c s 2 remains larger than 0.1. This is the case for all masses M m a x 2.1 M , as demonstrated in Table 1. Only for even heavier neutron stars does this evidence turn out to be lower. It is worth noting that the inclusion of the heavy black widow pulsar PSR J0952–0607 plays an important role that strengthens the evidence against c s 2 0.1 further in comparison with previous analyses, which did not include this object [29]. Squared sound speeds c s 2 > 0.1 in the cores of neutron stars with M m a x 2.0 M were also found in [24,26] without the inclusion of the black widow pulsar.
These limiting bounds on very low speeds of sound suggest that the appearance of a strong first-order phase transition with Maxwell construction is unlikely for most of the known neutron stars. An additional feature of a Maxwell-constructed first-order transition is an extended phase coexistence region. This will be examined more closely in the following subsection.

3.2. Evidence against a Strong First-Order Phase Transition in the Cores of Neutron Stars

The relatively moderate baryon densities inferred in the cores of neutron stars together with the evidence against very small sound speeds can be complemented by an additional study that is specific to first-order phase transitions with Maxwell construction. Such phase transitions are characterized by a domain of phase coexistence that extends over a certain range of baryon densities, Δ ρ . The width of this domain, Δ ρ / ρ (with ρ the density at which the coexistence interval starts), is a measure of the strength of the phase transition, i.e., the magnitude of surface tension between the two coexisting phases. We refer to a `strong’ first-order phase transition if the width of the mixed phase is Δ ρ / ρ > 1 . One example is the liquid–gas phase transition in symmetric nuclear matter at temperatures of T < 15 MeV [54,55]. In contrast, a `weak’ first-order phase transition has a value of Δ ρ / ρ that is small compared to unity.
The posterior credible bands inferred from neutron star data, as displayed in Figure 2, permit a systematic study of the maximum possible phase coexistence widths over the range of relevant energy densities. Starting from P ( ε ) , the Gibbs–Duhem relation is used to re-express pressure as a function of baryon density, P ( ρ ) . In this way, posterior credible bands for P ( ρ ) can be derived. These constrain the maximum possible widths of mixed-phase domains, ( Δ ρ / ρ ) max . An example is shown in Figure 8. In the analysis in [31], it turned out that these maximum possible phase coexistence regions were narrow: ( Δ ρ / ρ ) max 0.2 at the 68% level. Even within the 95% credible bands, this width did not exceed 0.3. In fact, one finds that ( Δ ρ / ρ ) max stays nearly constant as a function of baryon density ρ (taken at the starting point of the mixed-phase interval) over the whole region of densities ρ 2 5 ρ 0 that are relevant for neutron stars.
We note that the study performed in [56] also discussed restrictive constraints on strong first-order phase transitions but did not include the NICER and J0952-0607 black widow pulsar data. Even with this much more limited database, the authors found possible phase coexistence intervals of only Δ ρ / ρ 0.5 for first-order phase transitions within neutron stars. As can be expected, the additional radius constraints imposed by the NICER data constrain the maximum extension of the phase coexistence region even further. Moreover, we reiterate the additional restrictive impact of J0952-0607 in constraining the possible occurrence of phase transitions. This is at variance with conclusions drawn about a strong deconfinement transition in neutron star matter [18,26] on the basis of a dataset that did not yet include the heaviest pulsar observed so far.

3.3. Intermediate Summary

The conclusion drawn in our reported analysis is that only a weak first-order phase transition with a mixed-phase width Δ ρ / ρ 0.2 0.3 in a Maxwell construction can be realized inside neutron stars within the posterior credible bands of the inferred equation of state. A strong first-order phase transition is, thus, very unlikely to occur [31]. It is worth pointing out that this conclusion is entirely based on the analysis of observational data, independent of any specific EoS model. As will be discussed later, this rules out some relativistic mean-field models that suggested the occurrence of a strong (chiral) first-order phase transition already at densities below those encountered in neutron star cores. On the other hand, the empirical uncertainty bands still leave room for a continuous transition, such as a hadron-to-quark crossover [49,50,51,52]. The possibility of a purely baryonic equation of state with a monotonically rising speed of sound is also not excluded [29].

4. Phenomenology and Models

The previous sections have shown how present and forthcoming neutron star data (masses, radii, and tidal deformabilities) provide constraints for the equation of state and for possible phase transitions in cold, dense baryonic matter. The enhanced stiffness of the EoS at high baryon densities is a necessary feature; sufficiently high pressures are required to support neutron stars with masses around and beyond 2 M at relatively large radii ( R 11–13 km). The central core densities are not as extreme as was previously imagined. Even in the core of a 2.3 M neutron star, the baryon density at its center stays below five times the equilibrium density of nuclear matter (at the 68% credible level; see Equation (12)). The upper border of the 95% credible interval of the inferred central density ρ c ( 2.3 M ) / ρ 0 = 3.8 1.3 + 1.6 [31] lies slightly above 5 ρ 0 , but with a very low probability. This observation has consequences for the interpretation of the possible structure and composition of neutron star matter, which we discuss in the following.

4.1. Reminder of Low-Energy Nucleon Structure and a Two-Scale Scenario

Spontaneously broken chiral symmetry, as the long-wavelength manifestation of QCD, governs the low-energy structure and dynamics of hadrons, including nucleons and pions. As chiral Nambu–Goldstone bosons, pions play a distinguished role in this context. Models based on chiral symmetry often view the nucleon as a complex system of two scales [57]: a compact hard core that hosts the three valence quarks and, thus, encloses the baryon number and a surrounding quark–antiquark cloud in which pions figure prominently as the `soft’ degrees of freedom.
Such a two-scale scenario is manifested in empirical form factors and sizes of nucleons. Consider, for example, the proton and neutron electromagnetic form factors and their slopes at zero momentum transfer, which determine the corresponding mean squared radii. The empirical r.m.s. proton charge radius, r p 2 1 / 2 = 0.840 ± 0.003 fm, combined with six times the slope of the neutron electric form factor, r n 2 = 0.105 ± 0.006 fm 2 , gives the isoscalar and isovector mean square radii of the nucleon, r S , V 2 = r p 2 ± r n 2 , with the following resulting values [58]:
r S 2 0.78 fm , r V 2 0.90 fm .
Each of the nucleon form factors G i ( q 2 ) related to a current with index i has a representation in terms of a dispersion relation,
G i ( q 2 ) = G i ( 0 ) + q 2 π t 0 d t Im G i ( t ) t ( t q 2 i ϵ ) ,
with the squared four-momentum transfer q 2 = q 0 2 q 2 . The mean square radii are given as
r i 2 = 6 G i ( 0 ) d G i ( q 2 ) d q 2 | q 2 = 0 = 6 π t 0 d t t 2 S i ( t ) ,
where the distribution S i ( t ) = Im G i ( t ) / G i ( 0 ) represents the spectrum of intermediate hadronic states through which the external probing field couples to the respective nucleon current. For example, the isovector charge radius reflects the interacting two-pion cloud of the nucleon governed by the ρ meson and a low-mass tail extending down to the π π threshold, t 0 = 4 m π 2 . The isoscalar charge radius is related to the three-pion spectrum, which is strongly dominated by the narrow ω meson [59] and starts at t 0 = 9 m π 2 . The isoscalar charge form factor of the nucleon, G S E ( q 2 ) (with G S E ( 0 ) = 1 ), is particularly suitable for discussing a delineation between the `core’ and `cloud’ parts of the nucleon [60,61]. The vector meson dominance principle implies, in its simplest version, a representation of the form
G S E ( q 2 ) = F B ( q 2 ) 1 + | q 2 | / m ω 2 .
The form factor F B ( q 2 ) of the baryon number distribution in the nucleon core acts as a source for the ω field that propagates with its mass m ω . Introducing the mean square radius of the baryon core, r B 2 = 6 d F B ( q 2 ) d q 2 | q 2 = 0 , the mean square isoscalar charge radius becomes
r S 2 = r B 2 + 6 m ω 2 .
Using m ω = 783 MeV, the estimated baryonic core radius is
r B 2 0.47 fm .
A nucleon core size of about 1/2 fm is characteristic of chiral `core + cloud’ models. It also holds up in more detailed treatments of the spectral distributions governing the nucleon form factors [58]. The inclusion of additional ϕ meson and π ρ continuum contributions in the spectral function of G S E ( q 2 ) moves the core radius to just slightly larger values.
Consider as another example the form factor associated with the axial vector current of the nucleon. The corresponding mean-square axial radius deduced from neutrino–deuteron scattering data is reported as [62]
r A 2 = ( 0.46 ± 0.22 ) fm 2 .
A schematic axial vector dominance picture would assign a dominant part of the `cloud’ contribution to this form factor through the spectrum of the a 1 meson with its large width. If an approximate scale of this `cloud’ part is identified with δ r A 2 6 / m a 2 using the physical a 1 mass, m a 1.23 GeV, one finds
r A 2 core = [ r A 2 δ r A 2 ] 1 / 2 0.55 fm
with an estimated uncertainty of about 25%.
Yet another interesting piece of information is the mass radius of the proton deduced from J / Ψ photoproduction data [63]. It involves the form factor of the trace, T μ μ , of the nucleon’s energy–momentum tensor and is supposed to be dominated by gluon dynamics at the center of the nucleon. Low-mass (e.g., two-pion) components are suppressed. The result quoted in [63],
r 2 mass 1 / 2 = ( 0.55 ± 0.03 ) fm ,
is, once again, remarkably close to an assumed `core’ size scale of 1 / 2 fm.
These empirical considerations motivate a picture of the nucleon as containing a compact `hard’ core in which the valence quarks and their baryon number are confined, and a `soft’ surface of quark–antiquark pairs forming a mesonic cloud. This structure has implications for the behavior of nucleons in dense baryonic matter. With a core size of R core 1 / 2 fm and a cloud range given, e.g., by the proton charge radius, R cloud 0.84 fm, there is a significant separation of volume scales in vacuum: ( R cloud / R core ) 3 5 .
This scale separation is expected to increase further in dense baryonic matter for the following reasons. The properties of the soft multi-pion cloud are closely tied to spontaneously broken chiral symmetry and the approximate Nambu–Goldstone boson nature of the pion. The size of this cloud is expected to increase with density, along with the decreasing in-medium pion decay constant, f π * ( ρ ) , which acts as a chiral order parameter. The baryonic core, on the other hand, is governed by gluon dynamics, without a leading connection to chiral symmetry in QCD. This core is, therefore, expected to be quite stable against changes from increasing density up until the compact hard cores begin to touch and overlap.
What arises in this way is a two-scale scenario for dense baryonic matter, as described in [51] and sketched in Figure 9. At ρ ρ 0 , the tails of the meson clouds of nucleon pairs overlap, resulting in two-body exchange forces. As the average distance between nucleons decreases with increasing density, around ρ 2 3 ρ 0 , the soft clouds of q q ¯ pairs start to be delocalized. Their mobility expands over larger distances in a percolation process involving larger numbers of nucleons. In the terms of conventional nuclear physics, this corresponds to the emergence of many-body forces, the strength of which grows with increasing density. At that stage, the baryonic cores are still separated but subject to increasingly repulsive Pauli principle effects. The cores begin to touch and overlap at average nucleon–nucleon distances of d 1 fm, corresponding to densities ρ 5 ρ 0 . Further compression of baryonic matter would still have to overcome the strongly repulsive hard core in the nucleon–nucleon interaction. Recalling the inferred credible intervals of the central densities in heavy neutron stars (Equation (12)), one concludes that a phase transition to valence quark deconfinement does not seem likely in a two-scale picture and under the conditions in the cores of neutron stars unless they are extremely heavy.
An interesting and closely related result emerges from a detailed analysis of y scaling in electron–nucleus scattering at large momentum transfers ( | q | 1 GeV) and low energy transfers [64]. The persistently observed y scaling under these kinematical conditions implies that the electrons still scatter from strongly correlated pairs of nucleons, rather than quarks, at short distances. The conclusion drawn in [64] is that at local densities as large as five times ρ 0 , nuclear matter still appears to behave as a collection of nucleons.

4.2. Quark–Hadron Continuity and Crossover

While a strong first-order phase transition in neutron star matter seems unlikely based on the current empirical observations, a continuous crossover from hadrons to quarks is still possible within the present data-driven constraints. Such a scenario is realized, for example, in the QHC21 equation of state [52]. It features a smooth interpolation between low and high densities from nuclear to quark matter regimes. The quark matter phase is described by a three-flavor Nambu–Jona-Lasinio (NJL) model that includes pairing (diquark) degrees of freedom and a strongly repulsive vector coupling between quarks. In order to build up the necessary pressure at high densities, this vector coupling must be comparable to or larger than the strength of the standard scalar–pseudoscalar interaction in the NJL model. The density at which the interpolated turnover to quark matter takes place in the QHC21 EoS is chosen as ρ 3.5 ρ 0 , within the range of central densities that can be reached in M 2 M neutron stars. A crossover from hadrons to quarks could also involve an intermediate phase of quarkyonic matter [49,65].

4.3. Restoration of Chiral Symmetry in Dense Matter: From First-Order Phase Transition to Crossover

The quest for chiral symmetry restoration at high baryon densities—a transition from the spontaneously broken Nambu–Goldstone realization to the unbroken Wigner–Weyl mode—has been a persistent issue for a long time. A possible first-order chiral phase transition and the existence of a corresponding critical endpoint in the QCD phase diagram have always been topics of prime interest [66,67,68]. Early hypotheses concerning the occurrence of a first-order phase transition were frequently based on Nambu–Jona-Lasinio (NJL)-type models in mean-field approximation [69,70,71], which were later extended by incorporating confinement aspects through added Polyakov-loop degrees of freedom [72,73,74].
The empirical constraints on strong first-order phase transitions in dense neutron star matter, as described in Section 3.2, are strikingly at variance with previous mean-field (e.g., NJL model) predictions. These suggested that a first-order chiral phase transition should already appear at relatively low baryon densities around ρ 2 3 ρ 0 . A possible explanation for this discrepancy can be found in [55,75], where a chiral nucleon–meson (ChNM) field theory was used to explore the effects of fluctuations beyond mean-field (MF) approximation. The starting point was a relativistic chiral Lagrangian, L ( Ψ ; π , σ ; v μ ) , shaped around a linear sigma model with nucleons ( Ψ ) , pions ( π ) , and a scalar ( σ ) field. Short-range dynamics were parametrized in terms of heavy isoscalar and isovector vector fields ( v μ ) . The expectation value σ of the scalar field acted as a chiral order parameter normalized in the vacuum to the pion decay constant, f π 92 MeV. Two classes of fluctuations beyond MF were then systematically studied: first, vacuum fluctuations that introduced an additional term proportional to σ 4 ln ( σ / f π ) in the MF partition function; secondly, fluctuations involving pion loops and nucleon particle–hole excitations. The vacuum fluctuations can be included in an extended mean-field (EMF) approximation [76]. Fluctuations involving pion and nucleon loops are computed using non-perturbative functional renormalization group (FRG) methods. The parameters of the ChNM model—in particular, those related to short-distance dynamics—are fixed to reproduce empirical nuclear physics data [55,75].
Figure 10 demonstrates the important role of fluctuations beyond the mean-field approximation for the chiral order parameter σ . In symmetric nuclear matter, the mean-field approximation of the ChNM model correctly reproduces the first-order liquid–gas phase transition at low density. However, the chiral order parameter also displays a strong first-order chiral phase transition with a Maxwell-constructed phase coexistence region starting already below 2 ρ 0 . For neutron matter, which has no liquid–gas phase transition, the MF approximation nevertheless predicts a first-order chiral phase transition at densities around ρ 3 ρ 0 , which is well within the range of densities realized in neutron stars. However, in both nuclear and neutron matter, the inclusion of fermionic vacuum fluctuations (i.e., the effect of the ground state zero-point energy) in the extended mean-field (EMF) approximation converts the first-order chiral phase transition into a smooth crossover and shifts it to densities of ρ > 5 ρ 0 . This effect is further enhanced by the additional fluctuations included in the full FRG calculation, as demonstrated in Figure 10. As a result, the restoration of chiral symmetry is relegated to very high baryon densities beyond the inferred core densities in even the heaviest neutron stars (see Equation (12)). A comparable impact of fluctuations on the phase structure is seen in alternative chiral models [77,78].
Another approach based on chiral symmetry is the parity-doublet model. In this model, the active coupled baryonic degrees of freedom are the nucleon with spin-parity 1 / 2 + and its chiral partner with spin-parity 1 / 2 , where the latter is identified with the N * ( 1535 ) resonance. Spontaneous chiral symmetry breaking in vacuum manifests itself in the mass splitting of these two states, while the in-medium restoration pattern of this symmetry is signaled by the N and N* masses becoming degenerate. A recent detailed analysis [79] of the chiral order parameter in this model using extended mean-field approximation found a chiral phase transition in nuclear matter, but at extremely high densities ( ρ > 10 ρ 0 ) that are far beyond the density scales reached in neutron stars.

4.4. Dense Baryonic Matter: A Fermi Liquid Picture

The constraints on the equation of state inferred from the current empirical data still permit an interpretation of neutron star core matter in terms of baryonic degrees of freedom, such as a system dominated by neutrons [41] with small fractions of protons and perhaps hyperons [80,81,82]. The inferred baryon chemical potential μ ( ρ ) = ε ( ρ ) ρ shown in Figure 3 does not distinguish between different species of baryons. Its behavior nonetheless displays increasingly strong correlations at high densities. It is instructive to analyze the gross properties of this state of matter using the Landau theory of relativistic Fermi liquids [83]. Here, we perform a schematic study assuming “neutron-like” quasiparticles [84] with a baryon number of B = 1 and a density-dependent mass of m ( ρ ) while ignoring other small admixtures in the composition of the dense medium. These quasiparticles are characterized by their (relativistic) Landau effective mass m L * at a Fermi momentum of p F = ( 3 π 2 ρ ) 1 / 3 ,
m L * ( ρ ) = p F 2 + m 2 ( ρ ) ,
together with an effective potential, U ( ρ ) , so that the baryon chemical potential can be written as [84]
μ ( ρ ) = m L * ( ρ ) + U ( ρ ) .
The median of the μ ( ρ ) posterior credible bands in Figure 3 is now taken as a guiding starting point to extract the baryonic quasiparticle properties. With an educated ansatz for m ( ρ ) , the density dependence of the potential U ( ρ ) can then be deduced and further discussed. One possible choice is to take m ( ρ ) (with m ( 0 ) = 939 MeV) from the non-perturbative FRG calculation employing the chiral nucleon–meson field-theoretical model [75] discussed in Section 4.1. The resulting Landau effective mass m L * ( ρ ) is shown in Figure 11 together with the potential U ( ρ ) = μ ( ρ ) m L * ( ρ ) . The 95% credible band of μ ( ρ ) in Figure 3 leads to an uncertainty of about 15% for U at high densities. It is instructive to fit the resulting quasiparticle potential by a series in powers of baryon density for ρ 5 ρ 0 (with ρ 0 = 0.16 fm 3 , as before):
U ( ρ ) = n u n ρ ρ 0 n .
The coefficients fitted to the median of μ ( ρ ) are
u 1 = 90.9 MeV , u 2 = 15.3 MeV , u 3 = 3.2 MeV , u 4 = 0.4 MeV .
This pattern reflects a hierarchy of many-body correlations, recalling that the term linear in density represents two-body interactions, the term of order ρ 2 corresponds to short-range three-body forces, and so forth. The role of the repulsive N-body terms with N > 2 is quite significant; at ρ 4 ρ 0 , corresponding to an average distance of about 1 fm between the baryonic quasiparticles, these terms contribute as much as the two-body forces to the potential U and generate the strong pressure to support heavy neutron stars. Of course, these statements rely on the ansatz for the density-dependent mass m ( ρ ) , which is guided by the FRG calculations of pure neutron matter. A small fraction of protons in beta-equilibrated matter will not substantially change this picture. However, neutron star core compositions that qualitatively deviate from this simplified picture may lead to different conclusions.
Finally, consider the dimensionless Fermi liquid parameters, F 0 and F 1 , of the spin-averaged quasiparticle interaction. In terms of the quasiparticle mass and potential, they are given as [84]
F 0 ( ρ ) = p F π 2 m ( ρ ) m ρ + m L * ( ρ ) U ρ , F 1 ( ρ ) = 3 U ( ρ ) μ ( ρ ) .
Further useful relations are 1 + F 0 = N ( 0 ) ( μ / ρ ) with the density of quasiparticle states at the Fermi surface, N ( 0 ) = m L * p F / π 2 , and 1 + F 1 / 3 = 1 U / μ = m L * / μ .
The result in Figure 12 shows a strongly increasing F 0 at high densities. This reflects once again the growing importance of many-body correlations as matter becomes more and more compact. Such repulsive correlations are responsible for the increase in the sound velocity beyond its canonical conformal limit, c s > 1 / 3 , as seen in Figure 1. The Fermi liquid parameter F 1 is smaller in magnitude and has a negative slope, indicating the decreasing effective mass at higher densities.
The Landau parameter F 0 displays the typical behavior of a strongly correlated Fermi system. However, it is interesting to observe that, in comparison with the leading Landau parameters in liquid 3 He, the correlations in neutron star core matter are not extraordinarily strong. Values of F 0 9.3 were reported for 3 He at zero pressure, and F 0 68.2 at a pressure of 27 bar [85]. Accordingly, the Fermi liquid parameters for matter in the density range that is reached inside even the heaviest neutron stars, ρ 5 ρ 0 , are significantly smaller. From this perspective, the dense baryonic medium encountered in the center of neutron stars is perhaps not as extreme as often imagined.

5. Concluding Remarks

The neutron star database, which was recently extended by adding the heaviest pulsar (PSR J0952–0607) observed so far, requires a much stiffer equation of state for neutron star matter compared to those of previous analyses. As a consequence, the baryon densities reached in the cores of even very heavy neutron stars do not exceed about five times the density of normal nuclear matter. In a baryon-dominated system, such densities correspond to average distances between any two baryons of 1 fm or larger. Chiral model descriptions of the nucleon suggest a scale separation between a compact valence quark core and a mesonic cloud, so the nucleonic cores with radii of about 1/2 fm would begin to touch and overlap only in the deep interior of extremely heavy neutron stars. The appearance of a strong first-order phase transition becomes unlikely under such conditions. This is consistent with the empirical results, which suggest that extended phase coexistence regions and small minimum sound speeds are not favored by the current data. This interpretation is further endorsed by the fact that fluctuations beyond the mean-field approximation convert a possible first-order chiral phase transition into a smooth crossover at high densities. However, this is at variance with conclusions about a strong deconfinement phase transition that have been drawn from datasets that do not yet include the 2.35 solar mass black widow pulsar. On the other hand, a continuous hadron-to-quark crossover scenario—or alternatively, baryonic matter as a strongly correlated relativistic Fermi liquid—remains possible within the inferred data-driven constraints.

Author Contributions

Both authors contributed equally to this research article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by Deutsche Forschungsgemeinschaft (DFG grant TRR110) and National Natural Science Foundation of China (NSFC grant no. 11621131001) through the Sino-German CRC110 “Symmetries and the Emergence of Structure in QCD” and by the DFG Excellence Cluster ORIGINS.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Stimulating discussions and communications with Jean-Paul Blaizot, Kenji Fukushima, Norbert Kaiser, and Mannque Rho are gratefully acknowledged. One of the authors (L.B.) thanks the Japan Society for the Promotion of Science and Tte University of Tokyo for their support and hospitality during his International Fellowship for Research in Japan in the summer of 2023.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lovato, A.; Dore, T.; Pisarski, R.D.; Schenke, B.; Chatziioannou, K.; Read, J.S.; Landry, P.; Danielewicz, P.; Lee, D.; Pratt, S.; et al. Long range plan: Dense matter theory for heavy-ion collisions and neutron stars. arXiv 2022, arXiv:2211.02224. [Google Scholar]
  2. Fukushima, K.; Hatsuda, T. The phase diagram of dense QCD. Rep. Prog. Phys. 2011, 74, 014001. [Google Scholar] [CrossRef]
  3. Aarts, G.; Aichelin, J.; Allton, C.; Athenodorou, A.; Bachtis, D.; Bonanno, C.; Brambilla, N.; Bratkovskaya, E.; Bruno, M.; Caselle, M.; et al. Phase transitions in particle physics—Results and perspectives from lattice quantum chromo-dynamics. Prog. Part. Nucl. Phys. 2023, 133, 104070. [Google Scholar]
  4. Demorest, P.B.; Pennucci, T.; Ransome, S.M.; Roberts, M.S.E. A two-solar mass neutron star measured using Shapiro delay. Nature 2010, 467, 1081–1083. [Google Scholar] [CrossRef]
  5. Antoniadis, J.; Freire, P.C.C.; Wex, N.; Tauris, T.M.; Lynch, R.S.; van Kerkwijk, M.H.; Kramer, M.; Bassa, C.; Dhillon, V.S.; Driebe, T.; et al. A massive pulsar in a compact relativistic binary. Science 2013, 340, 448–502. [Google Scholar] [CrossRef]
  6. Fonseca, E.; Pennucci, T.; Ellis, J.A.; Stairs, I.H.; Nice, D.J.; Ransom, S.M.; Demorest, P.B.; Arzoumanian, Z.; Crowther, K.; Dolch, T.; et al. The NANOGrav nine-year data set: Mass and geometric measurements of binary millisecond pulsars. Astrophys. J. 2016, 832, 167–189. [Google Scholar]
  7. Arzoumanian, Z.; Brazier, A.; Burke-Spolaor, S.; Chamberlin, S.; Chatterjee, S.; Christy, B.; Cordes, J.M.; Cornish, N.J.; Crawford, F.; Cromartie, H.T.; et al. The NANOGrav 11-year data set: High-precision timing of 45 millisecond pulsars. Astrophys. J. Suppl. 2018, 235, 37–78. [Google Scholar] [CrossRef]
  8. Cromartie, H.T.; Fonseca, E.; Ransom, S.M.; Demorest, P.B.; Arzoumanian, Z.; Blumer, P.R.; Brook, P.R.; DeCesar, M.E.; Dolch, T.; Ellis, J.A.; et al. Relativsitic Shapiro delay measurements of an extremely massive millisecond pulsar. Nat. Astron. 2020, 4, 72–76. [Google Scholar] [CrossRef]
  9. Fonseca, E.; Cromartie, H.T.; Pennucci, T.T.; Ray, P.S.; Kirichenko, A.Y.; Ransom, S.M.; Demorest, P.B.; Stairs, I.H.; Arzoumanian, Z.; Guillemot, L.; et al. Refined mass and geometric measurements of the high-mass PSR J0740+6620. Astrophys. J. Lett. 2021, 915, L12. [Google Scholar] [CrossRef]
  10. Riley, T.E.; Watts, A.L.; Bogdanov, S.; Ray, P.S.; Ludlam, R.M.; Guillot, S.; Arzoumanian, Z.; Baker, C.L.; Bilous, A.V.; Chakrabarty, D.; et al. A NICER view of PSR J0030+0451: Millisecond pulsar parameter estimation. Astrophys. J. Lett. 2019, 887, L21. [Google Scholar]
  11. Miller, M.C.; Lamb, F.K.; Dittmann, A.J.; Bogdanov, S.; Arzumanian, Z.; Gendreau, K.C.; Guillot, S.; Harding, A.K.; Ho, W.C.G.; Lattimer, J.M.; et al. PSR J0030+0451 mass and radius from NICER data and implications for the properties of neutron stars. Astrophys. J. Lett. 2019, 887, L24. [Google Scholar] [CrossRef]
  12. Riley, T.E.; Watts, A.L.; Ray, P.S.; Bogdanov, S.; Guillot, S.; Morsink, S.M.; Bilous, A.V.; Arzoumanian, Z.; Choudury, D.; Daneva, J.S.; et al. A NICER view of the massive pulsar PSR J0740+6620 informed by radio timing and XMM-Newton spectroscopy. Astrophys. J. Lett. 2021, 918, L27. [Google Scholar]
  13. Miller, M.C.; Lamb, F.K.; Dittmann, A.J.; Bogdanov, S.; Arzumanian, Z.; Gendreau, K.C.; Guillot, S.; Ho, W.C.G.; Lattimer, J.M.; Loewenstein, M.; et al. The radius of PSR J0740+6620 from NICER and XMM-Newton data. Astrophys. J. Lett. 2021, 918, L28. [Google Scholar] [CrossRef]
  14. Salmi, T.; Vinciguerra, S.; Choudury, D.; Riley, T.E.; Watts, A.L.; Remillard, R.A.; Ray, P.S.; Bogdanov, S.; Guillot, S.; Arzumanian, Z.; et al. The radius of PSR J0740+6620 from NICER with NICER background estimates. Astrophys. J. 2022, 941, 150–173. [Google Scholar]
  15. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. (LIGO Scientific Collaboration and Virgo Collaboration). Properties of the binary neutron star merger GW170817. Phys. Rev. 2019, X9, 011001. [Google Scholar] [CrossRef]
  16. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. (LIGO Scientific Collaboration and Virgo Collaboration). GW170817: Measurement of neutron star radii and equation of state. Phys. Rev. Lett. 2018, 121, 161101. [Google Scholar] [CrossRef] [PubMed]
  17. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Abraham, S.; Acernese, F.; Ackley, K.; Adams, C.; Adhikari, R.X.; Adya, V.B.; Affeldt, C.; et al. (LIGO Scientific Collaboration and Virgo Collaboration). GW190425: Observation of a compact binary coalescence with total mass ∼ 3.4 M. Astrophys. J. Lett. 2020, 892, L3. [Google Scholar] [CrossRef]
  18. Annala, E.; Gorda, T.; Kurkela, A.; Nättilä, J.; Vuorinen, A. Evidence for quark-matter cores in massive neutron stars. Nature Phys. 2020, 16, 907. [Google Scholar] [CrossRef]
  19. Raaijmakers, G.; Greif, S.K.; Hebeler, K.; Hinderer, T.; Nissanke, S.; Schwenk, A.; Riley, T.E.; Watts, A.; Lattimer, J.M.; Ho, W.C.G. Constraints on the dense matter equation of state and neutron star properties from NICER’s mass-radius estimate of PSR J0740+6620 and multimessenger observations. Astrophys. J. Lett. 2021, 918, L29. [Google Scholar]
  20. Pang, P.T.H.; Tews, I.; Coughlin, M.W.; Bulla, M.; Van Den Broeck, C.; Dietrich, T. Nuclear physics multimessenger constraints on the neutron star equation of state: Adding NICER’s PSR J0740+6620 measurement. Astrophys. J. 2021, 922, 14. [Google Scholar] [CrossRef]
  21. Legred, I.; Chatziiannou, K.; Essick, R.; Han, S.; Landry, P. Impact of the PSR J0740+6620 radius constraint on the properties of high-density matter. Phys. Rev. 2021, D104, 063003. [Google Scholar] [CrossRef]
  22. Biswas, B.; Datta, S. Constraining neutron star properties with a new equation of state insensitive approach. Phys. Rev. 2022, D106, 043012. [Google Scholar] [CrossRef]
  23. Ecker, C.; Rezzolla, L. A general, scale-independent description of the sound speed in neutron stars. Astrophys. J. Lett. 2022, 939, L35. [Google Scholar] [CrossRef]
  24. Altiparmak, S.; Ecker, C.; Rezzolla, L. On the sound speed in neutron stars. Astrophys. J. Lett. 2022, 939, L34. [Google Scholar] [CrossRef]
  25. Huth, S.; Pang, P.T.H.; Tews, I.; Dietrich, T.; Le Fevre, A.; Schwenk, A.; Trautmann, W.; Agarwal, K.; Bulla, M.; Coughlin, M.W.; et al. Constraining neutron star matter with microscopic and macroscopic collisions. Nature 2022, 606, 276–280. [Google Scholar] [CrossRef] [PubMed]
  26. Annala, E.; Gorda, T.; Hirvonen, J.; Komoltsev, O.; Kurkela, A.; Nättilä, J.; Vuorinen, A. Strongly interacting matter exhibits deconfined behavior in massive neutron stars. arXiv 2023, arXiv:2303.11356. [Google Scholar] [CrossRef] [PubMed]
  27. Somasundran, R.; Tews, I.; Margueron, J. Investigating signatures of phase transitions in neutron-star cores. Phys. Rev. 2023, C107, 025801. [Google Scholar]
  28. Essick, R.; Legred, I.; Chatziioannou, K.; Han, S.; Landry, P. Phase transition phenomenology with nonparametric representations of the neutron star equation of state. Phys. Rev. 2023, D108, 043013. [Google Scholar] [CrossRef]
  29. Brandes, L.; Weise, W.; Kaiser, N. Inference of the sound speed and related properties of neutron stars. Phys. Rev. 2023, D107, 014011. [Google Scholar] [CrossRef]
  30. Lim, Y.; Holt, J.W. Neutron star radii, deformabilities, and moments of inertia from experimental and ab-initio theory constraints of the 208Pb neutron skin thickness. Galaxies 2022, 10, 99. [Google Scholar] [CrossRef]
  31. Brandes, L.; Weise, W.; Kaiser, N. Evidence against a strong first-order phase transition in neutron star cores: Impact of new data. Phys. Rev. 2023, D108, 094014. [Google Scholar]
  32. Mroczek, D.; Miller, M.C.; Noronha-Hostler, J.; Yunes, N. Nontrivial features in the speed of sound inside neutron stars. arXiv 2023, arXiv:2309.02345. [Google Scholar]
  33. Han, M.-Z.; Huang, Y.-J.; Tang, S.-P.; Fan, Y.-Z. Plausible presence of new state in neutron stars with masses above 0.98 MTOV. Sci. Bull. 2023, 68, 913–919. [Google Scholar] [CrossRef] [PubMed]
  34. Drischler, C.; Han, S.; Reddy, S. Large and massive neutronn stars: Implications for the sound speed within QCD of dense matter. Phys. Rev. 2022, C105, 035808. [Google Scholar]
  35. Gorda, T.; Kurkela, A.; Paatelainen, R.; Säppi, S.; Vuorinen, A. Soft interactions in cold quark matter. Phys. Rev. Lett. 2021, 127, 162003. [Google Scholar] [PubMed]
  36. Komoltsev, O.; Kurkela, A. How perturbative QCD constrains the equation of state at neutron-star densities. Phys. Rev. Lett. 2022, 128, 202701. [Google Scholar]
  37. Romani, R.W.; Kandel, D.; Filippenko, A.W.; Brink, T.G.; Zheng, W.K. PSR J0952–0607: The fastest and heaviest known galactic neutron star. Astrophys. J. Lett. 2022, 934, L17. [Google Scholar]
  38. Fasano, M.; Abdelsalhin, T.; Maselli, A.; Ferrari, V. Constraining the neutron star equation of state using multiband independent measurements of radii and tidal deformabilities. Phys. Rev. Lett. 2019, 123, 141101. [Google Scholar] [CrossRef]
  39. Essick, R.; Tews, I.; Landry, P.; Reddy, S.; Holz, D.E. Direct astrophysical tests of chiral effective field theory at supranuclear densities. Phys. Rev. 2020, C102, 055803. [Google Scholar] [CrossRef]
  40. Komoltsev, O.; Somasundaram, R.; Gorda, T.; Kurkela, A.; Margueron, J.; Tews, I. Equation of state at neutron-star densities and beyond from perturbative QCD. arXiv 2023, arXiv:2312.14127. [Google Scholar]
  41. Akmal, A.; Pandharipande, V.R.; Ravenhall, D.G. The equation of state of nucleon matter and neutron star structure. Phys. Rev. 1998, C58, 1804–1828. [Google Scholar]
  42. Fujimoto, Y.; Fukushima, K.; McLerran, L.D.; Praszalowicz, M. Trace anomaly as signature of conformality in neutron stars. Phys. Rev. Lett. 2022, 129, 252702. [Google Scholar] [CrossRef]
  43. Marczenko, M.; McLerran, L.D.; Redlich, K.; Sasaki, C. Reaching percolation and conformal limits in neutron stars. Phys. Rev. 2023, C107, 025802. [Google Scholar] [CrossRef]
  44. Rho, M. Dense baryonic matter predicted in “pseudo-conformal model”. Symmetry 2023, 15, 1271. [Google Scholar] [CrossRef]
  45. Ma, Y.-L.; Yang, W.-C. Topology and emergent symmetries in dense compact star matter. Symmetry 2023, 15, 776. [Google Scholar] [CrossRef]
  46. Nättilä, J.; Miller, M.C.; Steiner, A.W.; Kajava, J.J.E.; Suleimanov, V.F.; Poutanen, J. Neutron star mass and radius measurements from atmospheric model fits to x-ray burst cooling tail spectra. Astron. Astrophys. 2017, 608, A31. [Google Scholar]
  47. Sumiyoshi, K.; Kojo, T.; Furusawa, S. Equation of state in neutron stars and supernovae. In Handbook of Nuclear Physics; Tanihata, I., Toki, H., Kajino, T., Eds.; Springer: Singapore, 2023; Volume 5, pp. 3127–3177. [Google Scholar]
  48. Holt, J.W.; Rho, M.; Weise, W. Chiral symmetry and effective field theories for hadronic, nuclear and stellar matter. Phys. Reports 2016, 621, 2–75. [Google Scholar]
  49. McLerran, L.; Reddy, S. Quarkyonic matter and neutron stars. Phys. Rev. Lett. 2019, 122, 122701. [Google Scholar] [CrossRef]
  50. Baym, G.; Furusawa, S.; Hatsuda, T.; Kojo, T.; Togashi, H. New neutron star equation of state with quark-hadron crossover. Astrophys. J. 2019, 885, 42–49. [Google Scholar]
  51. Fukushima, K.; Kojo, T.; Weise, W. Hard-core deconfinement and soft-surface delocalization from nuclear to quark matter. Phys. Rev. 2020, D102, 096017. [Google Scholar] [CrossRef]
  52. Kojo, T.; Baym, G.; Hatsuda, T. Implications of NICER for neutron star matter: The QHC21 Equation of State. Astrophys. J. 2022, 934, 46–58. [Google Scholar] [CrossRef]
  53. Han, S.; Mamun, M.A.A.; Lalit, S.; Constantinou, C.; Prakash, M. Treating quarks within neutron stars. Phys. Rev. 2019, D100, 103022. [Google Scholar] [CrossRef]
  54. Wellenhofer, C.; Holt, J.W.; Kaiser, N.; Weise, W. Nuclear thermodynamics from chiral low-momentum interactions. Phys. Rev. 2014, C89, 064009. [Google Scholar]
  55. Brandes, L.; Kaiser, N.; Weise, W. Fluctuations and phases in baryonic matter. Eur. Phys. J. 2021, A57, 243. [Google Scholar] [CrossRef]
  56. Gorda, T.; Hebeler, K.; Kurkela, A.; Schwenk, A.; Vuorinen, A. Constraints on strong phase transitions in neutron stars. Astrophys. J. 2023, 955, 100. [Google Scholar] [CrossRef]
  57. Thomas, A.W.; Weise, W. The Structure of the Nucleon; Wiley-VCH: Berlin, Germany, 2001; pp. 232–246. [Google Scholar]
  58. Lin, Y.H.; Hammer, H.W.; Meißner, U.-G. New insights into the nucleon’s electromagnetic structure. Phys. Rev. Lett. 2022, 128, 052002. [Google Scholar]
  59. Kaiser, N.; Passemar, E. Spectral functions of nucleon form factors: Three-pion continua at low energies. Eur. Phys. J. 2019, A55, 16. [Google Scholar] [CrossRef]
  60. Brown, G.E.; Rho, M.; Weise, W. Phenomenological delineation of the quark-gluon structure from nucleon electromagnetic form factors. Nucl. Phys. 1986, A454, 669–690. [Google Scholar] [CrossRef]
  61. Meissner, U.-G.; Kaiser, N.; Weise, W. Nucleons as Skyrme solitons with vector mesons. Nucl. Phys. 1987, A466, 685–723. [Google Scholar] [CrossRef]
  62. Hill, R.J.; Kammel, P.; Marciano, W.J.; Sirlin, A. Nucleon axial radius and muonic hydrogen—A new analysis and review. Rep. Prog. Phys. 2018, 81, 096301. [Google Scholar] [CrossRef]
  63. Kharzeev, D.E. Mass radius of the proton. Phys. Rev. 2021, D104, 054015. [Google Scholar]
  64. Benhar, O. Testing the paradigm of nuclear many-body theory. Particles 2023, 6, 611–621. [Google Scholar] [CrossRef]
  65. Fujimoto, Y.; Kojo, T.; McLerran, L.D. Momentum shell in quarkyonic matter from explicit duality: A solvable model analysis. arXiv 2023, arXiv:2306.04304. [Google Scholar]
  66. Stephanov, M.; Rajagopal, K.; Shuryak, E. Signatures of the tricritical point in QCD. Phys. Rev. Lett. 1998, 81, 4816–4819. [Google Scholar]
  67. Fukushima, K. Spectral functions in the σ channel near the critical end point. Phys. Rev. 2003, C67, 025003. [Google Scholar]
  68. Fischer, C.S.; Luecker, J.; Welzbacher, C.A. Phase structure of three and four flavor QCD. Phys. Rev. 2014, D90, 034022. [Google Scholar]
  69. Asakawa, M.; Yazaki, K. Chiral restoration at finite density and temperature. Nucl. Phys. 1989, A504, 668–684. [Google Scholar]
  70. Klimt, S.; Lutz, M.F.M.; Weise, W. Chiral phase transition in the SU(3) Nambu and Jona-Lasinio model. Phys. Lett. 1990, B249, 386–390. [Google Scholar]
  71. Scavenius, O.; Mocsy, A.; Mishustin, I.N.; Rischke, D.H. Chiral phase transition within effective models with constituent quarks. Phys. Rev. 2001, C64, 045202. [Google Scholar] [CrossRef]
  72. Rößner, S.; Ratti, C.; Weise, W. Polyakov loop, diquarks and the two-flavour phase diagram. Phys. Rev. 2007, D75, 034007. [Google Scholar]
  73. Rößner, S.; Hell, T.; Ratti, C.; Weise, W. The chiral and deconfinement crossover transitions: PNJL model beyond mean field. Nucl. Phys. 2008, A814, 118–143. [Google Scholar]
  74. Hell, T.; Rößner, S.; Christoforetti, M.; Weise, W. Thermodynamics of the three-flavor nonlocal Polyakov–Nambu–Jona-Lasinio model. Phys. Rev. 2010, D81, 074023. [Google Scholar]
  75. Drews, M.; Weise, W. Functional renormalization group studies of nuclear and neutron matter. Prog. Part. Nucl. Phys. 2017, 93, 69–107. [Google Scholar]
  76. Skokov, V.; Friman, B.; Nakano, E.; Redlich, K.; Schaefer, B.-J. Vacuum fluctuations and thermodynamics od chiral models. Phys. Rev. 2010, D82, 034029. [Google Scholar]
  77. Zacci, A.; Schaffner-Bielich, J. Effects of renormalizing the chiral SU(2) quark-meson model. Phys. Rev. 2018, D97, 074011. [Google Scholar] [CrossRef]
  78. Gupta, U.S.; Tiwari, V.K. Revisiting the phase structure of the Polyakov-quark-meson model in the presence of vacuum fermion fluctuation. Phys. Rev. 2012, D85, 014010. [Google Scholar]
  79. Eser, J.; Blaizot, J.-P. Thermodynamics of the parity-doublet model: Symmetric nuclear matter and the chiral transition. arXiv 2023, arXiv:2309.06566. [Google Scholar]
  80. Lonardoni, D.; Lovato, A.; Gandolfi, S.; Pederiva, F. Hyperon puzzle: Hints from Quantum Monte Carlo calculations. Phys. Rev. Lett. 2015, 114, 092301. [Google Scholar] [CrossRef]
  81. Gerstung, D.; Kaiser, N.; Weise, W. Hyperon-nucleon three-body forces and strangeness in neutron stars. Eur. Phys. J. 2020, A56, 175. [Google Scholar]
  82. Leong, J.; Motta, T.F.; Thomas, A.W.; Guichon, P.A.M. Dense nuclear matter with phenomenological short distance repulsion. Phys. Rev. 2023, C108, 015804. [Google Scholar]
  83. Baym, G.; Chin, S.A. Landau theory of relativistic Fermi liquids. Nucl. Phys. 1976, A262, 527–538. [Google Scholar]
  84. Friman, B.; Weise, W. Neutron star matter as a relativistic Fermi liquid. Phys. Rev. 2023, C100, 065807. [Google Scholar]
  85. Baym, G.; Pethick, C. Landau Fermi-Liquid Theory: Concepts and Applications; Wiley: New York, NY, USA, 1991; p. 117. [Google Scholar]
Figure 1. Squared speed of sound [31] as a function of the energy density inferred from the empirical dataset listed in Section 2.1. Median (solid curve) and highest posterior density credible bands at levels of 68% (gray band) and 95% (dashed curves) are displayed. Intervals are indicated for the 68% ranges of central energy densities in the cores of 1.4 and 2.3 M neutron stars.
Figure 1. Squared speed of sound [31] as a function of the energy density inferred from the empirical dataset listed in Section 2.1. Median (solid curve) and highest posterior density credible bands at levels of 68% (gray band) and 95% (dashed curves) are displayed. Intervals are indicated for the 68% ranges of central energy densities in the cores of 1.4 and 2.3 M neutron stars.
Symmetry 16 00111 g001
Figure 2. Equation of state P ( ε ) [31] deduced from the inferred sound velocity based on the empirical dataset listed in Section 2.1. The median, 68%, and 95% posterior credible bands are displayed, as in Figure 1. The APR EoS [41] (dotted line) is shown for comparison. Also shown for orientation is the baryon density scale, ρ / ρ 0 (in units of the equilibrium density of nuclear matter ρ 0 = 0.16 fm 3 ), which was computed for the median of P ( ε ) .
Figure 2. Equation of state P ( ε ) [31] deduced from the inferred sound velocity based on the empirical dataset listed in Section 2.1. The median, 68%, and 95% posterior credible bands are displayed, as in Figure 1. The APR EoS [41] (dotted line) is shown for comparison. Also shown for orientation is the baryon density scale, ρ / ρ 0 (in units of the equilibrium density of nuclear matter ρ 0 = 0.16 fm 3 ), which was computed for the median of P ( ε ) .
Symmetry 16 00111 g002
Figure 3. Baryon chemical potential μ ( ρ ) in neutron star matter, normalized at μ ( ρ = 0 ) = 939 MeV. Posterior credible bands [31] at the 68% level (gray band) and 95% level (dashed lines) are shown.
Figure 3. Baryon chemical potential μ ( ρ ) in neutron star matter, normalized at μ ( ρ = 0 ) = 939 MeV. Posterior credible bands [31] at the 68% level (gray band) and 95% level (dashed lines) are shown.
Symmetry 16 00111 g003
Figure 4. Median and credible bands [31] of the trace anomaly measure, Δ = 1 3 P ε , at the 68% level (gray band) and 95% level (dashed lines).
Figure 4. Median and credible bands [31] of the trace anomaly measure, Δ = 1 3 P ε , at the 68% level (gray band) and 95% level (dashed lines).
Symmetry 16 00111 g004
Figure 5. Posterior credible bands [31] of the radius R as a function of the neutron star mass M at the 68% level (gray band) and 95% level (dashed line) compared to the analysis of NICER data by Riley et al. for PSR J0030+0451 and PSR J0740+6620 [10,12]. In addition, the mass–radius credible interval at the 68% level of the thermonuclear burster 4U 1702-429 (blue) is displayed [46] (which was not included in the Bayesian analysis).
Figure 5. Posterior credible bands [31] of the radius R as a function of the neutron star mass M at the 68% level (gray band) and 95% level (dashed line) compared to the analysis of NICER data by Riley et al. for PSR J0030+0451 and PSR J0740+6620 [10,12]. In addition, the mass–radius credible interval at the 68% level of the thermonuclear burster 4U 1702-429 (blue) is displayed [46] (which was not included in the Bayesian analysis).
Symmetry 16 00111 g005
Figure 6. Posterior credible bands [31] of the tidal deformability as a function of neutron star mass at the 68% level (gray band) and 95% level (dashed lines) compared to values deduced from the merger event GW170817 [15,38].
Figure 6. Posterior credible bands [31] of the tidal deformability as a function of neutron star mass at the 68% level (gray band) and 95% level (dashed lines) compared to values deduced from the merger event GW170817 [15,38].
Symmetry 16 00111 g006
Figure 7. Characteristic behaviors of the squared sound velocity in the presence of a phase transition or a crossover.
Figure 7. Characteristic behaviors of the squared sound velocity in the presence of a phase transition or a crossover.
Symmetry 16 00111 g007
Figure 8. Illustration of the constraint on the maximum width of a Maxwell-constructed coexistence region for a first-order phase transition within the 68% credible band of P ( ε ) . The upper (baryon density) scale refers to the median of the P ( ε ) distribution, as in Figure 2.
Figure 8. Illustration of the constraint on the maximum width of a Maxwell-constructed coexistence region for a first-order phase transition within the 68% credible band of P ( ε ) . The upper (baryon density) scale refers to the median of the P ( ε ) distribution, as in Figure 2.
Symmetry 16 00111 g008
Figure 9. Sketch of low- and high-density baryonic matter. Baryons (e.g., nucleons) are viewed as valence quark cores surrounded by clouds of quark–antiquark pairs (e.g., chiral meson clouds). At densities of ρ 2 3 ρ 0 , the percolation of quark–antiquark pairs over larger distances starts, as indicated. Valence quark cores begin to touch and overlap at baryon densities ρ 5 ρ 0 .
Figure 9. Sketch of low- and high-density baryonic matter. Baryons (e.g., nucleons) are viewed as valence quark cores surrounded by clouds of quark–antiquark pairs (e.g., chiral meson clouds). At densities of ρ 2 3 ρ 0 , the percolation of quark–antiquark pairs over larger distances starts, as indicated. Valence quark cores begin to touch and overlap at baryon densities ρ 5 ρ 0 .
Symmetry 16 00111 g009
Figure 10. Chiral order parameters in symmetric nuclear matter and neutron matter at a temperature of T = 0 as a function of baryon density in units of nuclear ground state equilibrium density, ρ 0 = 0.16 fm 3 . Dotted lines: liquid–gas phase transition (L-G) in symmetric nuclear matter. Dashed lines: first-order chiral phase transitions emerging from the mean-field (MF) approximation of a relativistic chiral nucleon–meson (ChNM) field-theoretical model. Solid lines show the results of extended mean-field (EMF) calculations (with the inclusion of fermionic vacuum fluctuations) and functional renormalization group (FRG) calculations based on the same ChNM model. The figures were adapted from Refs. [55,75].
Figure 10. Chiral order parameters in symmetric nuclear matter and neutron matter at a temperature of T = 0 as a function of baryon density in units of nuclear ground state equilibrium density, ρ 0 = 0.16 fm 3 . Dotted lines: liquid–gas phase transition (L-G) in symmetric nuclear matter. Dashed lines: first-order chiral phase transitions emerging from the mean-field (MF) approximation of a relativistic chiral nucleon–meson (ChNM) field-theoretical model. Solid lines show the results of extended mean-field (EMF) calculations (with the inclusion of fermionic vacuum fluctuations) and functional renormalization group (FRG) calculations based on the same ChNM model. The figures were adapted from Refs. [55,75].
Symmetry 16 00111 g010
Figure 11. The Landau effective mass m L * ( ρ ) = p F 2 + m 2 ( ρ ) and potential U ( ρ ) of quasiparticles representing the median of the posterior distribution of the baryon chemical potential μ ( ρ ) (see Figure 3).
Figure 11. The Landau effective mass m L * ( ρ ) = p F 2 + m 2 ( ρ ) and potential U ( ρ ) of quasiparticles representing the median of the posterior distribution of the baryon chemical potential μ ( ρ ) (see Figure 3).
Symmetry 16 00111 g011
Figure 12. The Landau Fermi liquid parameters derived from quasiparticle properties that are based on the median of the data-inferred baryon chemical potential μ (see Figure 3).
Figure 12. The Landau Fermi liquid parameters derived from quasiparticle properties that are based on the median of the data-inferred baryon chemical potential μ (see Figure 3).
Symmetry 16 00111 g012
Table 1. Bayes factors B c s , m i n 2 0.1 c s , m i n 2 > 0.1 comparing EoS samples with competing scenarios: (a) minimum squared speed of sound (following a maximum) with c s , m i n 2 > 0.1 versus (b) equations of state with c s , m i n 2 0.1 . The minimum speeds of sound are computed up to the maximum neutron star masses, as indicated (taken from [31]).
Table 1. Bayes factors B c s , m i n 2 0.1 c s , m i n 2 > 0.1 comparing EoS samples with competing scenarios: (a) minimum squared speed of sound (following a maximum) with c s , m i n 2 > 0.1 versus (b) equations of state with c s , m i n 2 0.1 . The minimum speeds of sound are computed up to the maximum neutron star masses, as indicated (taken from [31]).
M m a x / M 1.92.02.12.22.3
B c s , m i n 2 0.1 c s , m i n 2 > 0.1 500.9229.815.03.62.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Brandes, L.; Weise, W. Constraints on Phase Transitions in Neutron Star Matter. Symmetry 2024, 16, 111. https://doi.org/10.3390/sym16010111

AMA Style

Brandes L, Weise W. Constraints on Phase Transitions in Neutron Star Matter. Symmetry. 2024; 16(1):111. https://doi.org/10.3390/sym16010111

Chicago/Turabian Style

Brandes, Len, and Wolfram Weise. 2024. "Constraints on Phase Transitions in Neutron Star Matter" Symmetry 16, no. 1: 111. https://doi.org/10.3390/sym16010111

APA Style

Brandes, L., & Weise, W. (2024). Constraints on Phase Transitions in Neutron Star Matter. Symmetry, 16(1), 111. https://doi.org/10.3390/sym16010111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop