Cusp in the Symmetry Energy, PREX-II And Pseudo-Conformal Sound Speed in Neutron Stars

Hyun Kyu Lee,1, ∗ Yong-Liang Ma,2, 3, † Won-Gi Paeng,4, ‡ and Mannque Rho5, § 1Department of , , Seoul 133-791, 2School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, UCAS, Hangzhou, 310024, China 3International Center for Theoretical Physics Asia-Pacific (ICTP-AP) (Beijing/Hangzhou), UCAS, Beijing 100190, China 4A I Lab., Clunix, Seoul 07209, Korea 5Institut de Physique Th´eorique,Universit´eParis-Saclay, CNRS, CEA, 91191, Gif-sur-Yvette, France (Dated: July 23, 2021) The recent announcement of the PREX-II measurement of the neutron skin of 208Pb that suggests a stiff symmetry energy near nuclear matter density n0 and its impact on the EoS of massive compact stars raise the issue as to whether the widely accepted lore in nuclear astrophysics that the EoS determined at n0 necessarily gives a stringent “constraint” at high densities relevant to massive compact stars. We present the argument that the “cusp” structure in the symmetry energy at n1/2 >∼ 2n0 predicted by a topology change in dense matter could obstruct the validity of the lore. The topology change, encoding the emergence of QCD degrees of freedom in terms of hidden local and scale symmetries, predicts an EoS that is soft below and stiff above n >∼ n1/2, involving no low-order phase transitions, and yields the macrophysical properties of neutron stars overall consistent with the astrophysical observations including the maximum mass 2.0 <∼ M/M <∼ 2.2 as well as the GW data. Furthermore it describes the interior core of the massive stars constituted of baryon-charge-fractionalized quasi-fermions, that are neither baryonic nor quarkonic, with the 2 2 “pseudo-conformal” sound speed vpcs/c ≈ 1/3 converged from below at n1/2 with a nonzero trace of energy-momentum tensor. In the renormalization-group approach to interacting fermions dubbed GnEFT, the strangeness degrees of freedom play no role in the density regime relevant to the massive stars considered.

I. INTRODUCTION available it has the power to go beyond the standard chi- ral effective field theory (SχEFT), currently heralded as In accessing dense neutron-star matter in terms of a as a possible “first-principles approach” to nuclear the- topology change for the putative -quark continu- ory at low energy and densitiy, and gives extremely sim- ity, it was discovered in 2011 [1] that a cusp is present in ple predictions that have the merit to be unambiguously confronted by experiments in the density regime inacces- the nuclear symmetry energy Esym at a density ∼ (2 − 3) sible by SχEFT. It has thus far accounted with no un- times the normal nuclear matter density n = n0 ' 0.16 fm−3. This cusp structure has been found to play the surmountable tension for all macro-physical observables most important role in the approach to the EoS of dense available in both terrestrial and astrophysical laborato- compact-star matter developed entirely independently of ries. See for the current status, e.g., [2]. other on-going approaches in nuclear astrophysics. 1 For- In this paper, we show that this cusp structure zeroes mulated with the minimum number of degrees of freedom in on the recent issue raised by the PREX-II measure- ment of the neutron skin thickness of 208Pb [3] and the impact on the equation of state (EoS) of massive com- 208 pact stars. An analysis using the new Rskin and certain ∗ [email protected] correlations with the symmetry energy J and its slope L † [email protected] ‡ (to be defined) at n = n0 led to the 1 σ intervals [4] arXiv:2107.01879v2 [nucl-th] 22 Jul 2021 [email protected] § [email protected] J = (38.1 ± 4.7) MeV,L = (106 ± 37) MeV. (1) 1 When the idea of [1] was first submitted to Phys. Rev. Lett., it was summarily rejected by the DA Editor of PRL at the time on the ground that “the approach to was never considered as a viable alternative to more traditional ab initio calculations ... and is now more or less abandoned ...” matter, string theory etc. became much too difficult for many The DAE’s decision followed two concurrent events that took nuclear theorists – except for a few diehard enthusiasts. The sec- place then. First, all the mathematically astute – mostly par- ond event was that at the time the paper [1] was written, chiral ticle – theorists who contributed to resuscitating, in the light EFT was becoming very popular in the nuclear theory commu- of QCD, the Skyrme model which had been left totally ignored nity. It was also a lot easier to access numerically for quick results for two decades in nuclear physics community, quickly moved for publication. What transpired in this paper and is continu- over to string theory which was then beginning to undergo its ing to transpire in current developments, as briefly mentioned first “string revolution for theory of everything.” In the absence in Conclusion regarding the BPS skyrmion and holographic dual of those string theorists, making progress with highly intricate approaches to nuclear physics, aptly illustrates an unfortunately mathematics involved in skyrmion physics as is now being abun- prevalent short-sightedness on truly novel and unconventional dantly recognized in other areas of physics such as condensed ideas in the field too often dismissed as “conjectural.” 2

These values seemingly overshoot greatly the currently where α = (N −P )/(N +P ) is the neutron-proton asym- “accepted” values [5] 2, metry with P (N) standing for the number of protons (neutrons) in A = N + P system. The J and L J = (31.7 ± 1.1) MeV,L = (59.8 ± 4.1) MeV. (2) concerned are This result means the EoS must be a lot stiffer at normal ∂Esym(n) J = Esym(n0),L = 3 |n=n . (4) nuclear matter density than what has been considered up ∂n 0 to date. A similar observation termed as a “dilemma” is The issue associated with the Pb skin thickness puzzle arrived at by Piekarewicz from the electric dipole polar- involves this symmetry energy on which our discussion izability of neutron-rich nuclei [6]. Naively extrapolated will be focused. In standard nuclear physics approaches 208 to the massive compact-star density, the Rskin data could (SNPA) anchored on effective density functionals such as rule out most of, or at least put in serious tension, the the Skyrme potential, relativistic mean field (RMF) and EoS’ currently available in the literature for compact-star varieties thereof as well as SχPT up to manageable chi- physics. ral order, equipped with a certain number of parameters The stiff EoS implied by the dilemma turns out, as fit to available empirical data, the E(n, α) can be more we will discuss later, to have a dramatic effect on the or less reliably determined in the vicinity of the nuclear properties of massive stars such as the composition of matter equilibrium density n0. It has also been extended, the star core and sound speed. We will show that the with albeit significant uncertainty, up to slightly above cusp structure discovered in [1] gives rise to a totally dif- n0 from heavy-ion collision experiments. Thus one can ferent picture. This has a close connection to the lore say that the nuclear symmetry energy Esym is fairly well popularly accepted in certain nuclear astrophysics circles determined up to n0 in SNPAs. It should, however, be that the EoS determined accurately at low density, say, stressed that its slope in density, namely, L and higher as ∼ n0, should make an indispensable constraint to the derivatives remain uncertain, say in SχPT, unless chiral- m EoS at higher densities. Put differently, the lore that order terms up to N LO for m >∼ 4 are fully included. we shall refer to as nuclear-astrophysics lore (nLORE This is closely tied to the fact that the chiral power ex- κ for short) states that what happens in the core of com- pansion (say, in SχEFT) is bound to break down for kF pact stars must be constrained by what happens in nu- for κ >∼ 5 ((e.g., [9]) as the hadron-quark crossover den- clear matter. This of course must be true in a uniguely sity is approached. So the problem is how Esym and given theory, namely, QCD. However given that QCD can its derivatives behave beyond the equilibrium density n0. directly access neither nuclear matter nor compact-star This is where heavy degrees of freedom (HDsF) need to matter, what’s available is effective field theory (EFT) in enter. the sense defined by Weinberg’s Folk Theorem. In EFT, this nLORE cannot be valid if there are phase changes or crossovers. In fact, we will argue the presence of the cusp A. Cusp in Skyrmion Crystal at O(1/Nc) in our approach debunks the nLORE on constraints on EoS. What turns out to importantly figure in our argu- We address this problem exploiting a topological struc- ment is the existence of that cusp at a density >∼ (2−3)n0 ture of dense baryonic matter. This is because in the in the symmetry energy induced by a (what we consider large Nc limit in QCD, the only known non-perturbative to be robust) topology change in dense matter that ef- tool available in strong interaction physics applicable to fectively encodes the putative hadron-quark continuity baryonic matter at large density – in the absence of lat- expected in QCD. This point is signaled also in a dif- tice QCD – is putting (or instantons in holo- ferent context in [7]. It aptly reconciles a soft EoS at graphic QCD) on crystal lattice [10]. Application of the n <∼ n1/2 to a hard EoS at n > n1/2 accounting no- crystal skyrmion lattice method to nuclear matter and tably, among others, for massive >∼ 2M compact-stars dense matter has been around for some time (see for and other macroscopic star properties including the re- an early review [12]) but only recently is the power of cent gravity wave data. skyrmion approach beginning to be recognized in nuclear physics, contrary to condensed matter as well as string theory where the skyrmion structure in various spatial di- II. THE CUSP IN Esym mensions has been having remarkable impacts [11]. This is because of the extreme mathematical subtlety involved The quantity that plays the most important role in the in the skyrmion physics. Furthermore, the condition for EoS for compact-star matter [8] is the symmetry energy the validity of lattice skyrmions in particular, i.e. large Esym in the energy per nucleon given by Nc and large density, is not met at the density where there is a wealth of experimental data, namely finite nu- 2 4 E(n, α) = E(n, α = 0) + Esym(n)α + O(α ) + ··· ,(3) clei3. However the cusp structure in question that takes

3 There is a striking recent development that we will refer to below 2 We will elaborate on these “accepted” values below as the potential power of the skyrmion structure in finite nuclei 3 place at relatively high density – relative to normal nu- n-nucleon system is taken a nearly pure neutron system, clear matter – seems to meet the two conditions as in- Itot ≤ n/2, to the leading order of n for n → ∞, the dicated by the quasi-scale invariance seen in the crystal energy per baryon takes the form simulation in the half-skyrmion phase [13]. 1 2 To illustrate the basic idea, we first take the Skyrme E = Msol + α . (9) model [14] stabilized by the (Skyrme) quartic term for 8λI skyrmions supplemented by a scalar dilaton as first Thus the symmetry energy is shown in [1]. What is crucial is that the Skyrme model encodes the necessary topological structure. But by it- 1 2 Esym = + O(1/Nc ). (10) self, it misses certain nontrivial crucial dynamical char- 8λI acteristics encoded in QCD. We will implement these The moment of inertia λ ∼ O(N ) can be computed in missing ingredients with hidden local symmetry (HLS for I c the leading Nc order as the integral over the single cell short) supplemented with hidden scale symmetry (HSS) and takes the form and incorporate them for quantitative discussions in the generalized chiral effective field theory (EFT) approach f 2 λ = π 4 − 2φ2 + δλ + ··· , (11) that is dubbed GnEFT [2]. We will present the argument I 6 0 I that the HLS and HSS (combined, referred to as sHLS), where the first term comes from the quadratic current the degrees of freedom associated with which are identi- algebra term and the second stands for the contribu- fied with the HDsF involved at high density, are “dual” tion from the Skyrme quartic term which consists of four to QCD (gluons and quarks)4 in the density regime rele- terms involving φ and space derivatives of the chiral field vant to compact stars. It will be argued that the density 0 ξ. Here φ , proportional to the quark condensate hqq¯ i, involved is located far below asymptotic density at which 0 plays a crucial role in the whole development in [2]. hardon-quark continuity presumably does break down (to In the skyrmion crystal formalism, the topology change be specified below). is associated with the behavior of the quark condensate Following [1], we calculate the symmetry energy by at a density labeled n which should, and generically quantizing the crystal as a whole object through a col- 1/2 does, lie above n . The quark condensate Σ ≡ hqq¯ i, lective rotation in iso-space with the rotation angle C(t) 0 nonzero both globally and locally for n < n , goes to acting on the relevant chiral fields U = ξ2 (in unitary 1/2 zero at n when space averaged, φ ≡ Σ¯ → 0. But it gauge) as 1/2 0 is locally non-zero, thus generating chiral density wave † ξc(x) → ξ(x, x) = C(t)ξc(x)C (t), and giving rise to a non-vanishing decay constant, V (x) → V (x, t) = C(t)V (x)C†(t), (5) fπ 6= 0. This transition triggers a skyrmion in the matter µ,c µ µ,c to fractionalize into 2 half-skyrmions. Since the order where the subindex “c” means the static configuration parameter, here the pion decay constant, is non-zero in with the lowest energy for a given crystal size L and C(t) the changeover, there is no low-order phase transition. is a time-dependent unitary SU(2) matrix in isospace. This half-skyrmion “phase”5 resembles what is referred We define the angular velocity through Ω to as “pseudo-gap phase” in condensed matter physics, e.g., in high-T superconductivity. i τ · Ω = C†(t)∂ C(t). (6) An important – and most crucial – property of the sym- 2 0 metry energy in this formulation, ∝ 1/λI , is that it devel- The energy of the n-nucleon system can be written as ops a cusp at O(1/Nc) at the density n1/2 where φ0 → 0. The cusp structure seen in the skyrmion lattice simula- 1 tot 2 Mtot = Mstatic + λI Ω . (7) tion [1] is schematically depicted in Fig. 1 (right panel). 2 The exact location of the cusp depends on the parame- By regarding the angular momentum in isospace, J = ters of the Lagrangian which are `apriori unknown, so it is δMtot/δΩ, as the isospin operator, one can write the total arbitrary. It will be determined later from neutron-star energy of the system as observations to lie within the range 2 <∼ n1/2/n0 < 4. This cusp form comes from an interplay involving the 1 tot tot Mtot = nMsol + I (I + 1), (8) behavior of φ0 between the quadratic derivative current 2nλI algebra term and the countering contribution from the tot Skyrme quartic derivative term. Roughly what happens where Msol, λI and I are, respectively, the mass and moment of inertia of the single skyrmion in the system, is that the increase of λ from the quadratic term as φ0 and the total isospin of the n-nucleon. Given that the goes to zero is stopped by the quartic term at n1/2 and starts dropping, causing the cusp in 1/λ. It will be shown

relative to the nLORE. It will exhibit the indispensable role of HDsF. 5 Lacking a better terminology, we will continue to (mis)use this 4 This notion of hadron-quark duality will be specified below. term. 4

FIG. 1. Panorama of the symmetry energy Esym. Left panel (copied from [15]): Wilderness in both various nuclear models and SχEFTs and bounds given by (up-to-date) observations (solid blue lines). Right panel: Schematic form of the cusp (dotted line) in skyrmion crystal [1]. The solid line caricatures the effect of smearing by heavy degrees of freedom 0 (HDsF). The interval between n >∼ n0 and n1/2 is the density regime that is arguably the most difficult to access by standard χEFT from below and by QCD proper from above as discussed in the text in connection with the sound speed and the tidal deformation. that this picture will be modified in nature by, among We should stress that what’s involved in our approach others, two observations. First, the skyrmion crystal sim- is “hadron-quark duality,” not just hadron-quark conti- ulation cannot be trusted at low density below n/12, and nuity that captures crossover from hadronic degrees of next, the Skyrme quartic term can be taken as what re- freedom to quark/gluon degrees of freedom. In fact the sults from integrating out HDsF from the skyrmion La- notion of hadron-quark duality is a lot more general in grangian. The large Nc consideration gives a remarkably the sense elaborated recently in the Cheshire-Cat Prin- simple Esym. ciple [16]. It represents the necessity at densities ex- To illustrate what is captured in this cusp, we quote ceeding n0 of the “heavy degrees of freedom (HDsF).” in Fig. 1 (left panel) the recent illuminating summary Those HDsF are to encode the quarks/gluons degrees of by B.A. Li et al. [15] of the up-to-date experimental and QCD at some high density without explicit presence of theoretical status of the symmetry energy. It presents a quark/gluons. How to do this precisely is presently un- giant wilderness. All the theoretical models available up known in the density regimes relevant to compact stars. to date, e.g., various energy density functionals, χEFTs This is because the densities involved are too far from the etc,, fit $Esym by fiat to what’s given in nature at ∼ n0. asymptotic regime where perturbative QCD is applicable There are ample parameters available to allow it. The and the only nonperturbative tool known, lattice QCD, swamp sets in beyond n0. Given the absence of trustful is famously inaccessible at high density. So the question models – not to mention theories, there is no guidance is: How does one proceed? how the Esym should move at higher densities. There is nothing to prevent it from going up or down, even plunging below zero. The current experimental observa- B. Cusp Induced by Nuclear Tensor Force tions such as neutron stars (and heavy-ion data limited to only a few times n0) do not fare any better as indicated In [1], the cusp was reproduced by the role played by by the solid (blue) lines in the left panel of Fig. 1. the and the vector mesons in the nuclear tensor What is certain is that he cusp is buried in this jungle. force in standard nuclear structure physics. There the It may not be absurd to think that the cusp structure vector mesons were identified as hidden local fields and could be just an artifact of the lattice simulation. But it the scalar σ as a dilaton. The key idea there was to ex- turns out, we will see, that it is not. When the jungle ploit the vacuum-change-induced density dependence in is cleared up by the symmetries assumed to be involved, the hHLS Lagrangian in the presence of baryonic mat- the cusp yields an extremely simple and portent mech- ter [21]. Here we repeat essentially the same arguments anism needed for the EoS for massive stars. In particu- to bring out certain characteristics of sHLS hidden in lar, we will argue, the cusp represents the hadron-quark the discussions, namely, the “duality” assumed to hold “duality” expressed in topology change [2]. It will lead `ala Seiberg between hidden local gauge fields and QCD to what will be termed “pseudo-conformal sound speed” gluons and a -quarks/gluons duality. The ob- and baryon-charge fractionalized “confined” fermions in jective is to link it to what SχEFT does at low energy the core of neutron stars. (and density) and to extend it to higher densities where 5

SχEFT is to break down. This would make our approach fσ ≈ fπ 6= 0, accommodating massive matter fields, in to compact-star matter in line with the spirit of the Folk particular, light-quark baryons, `ala genuine dilaton sce- Theorem on EFT. At present, the duality assumed is only nario with the dilaton condensate dictating how hadron a conjecture, but there are several compelling indications masses scale in density [21]. The two effects entail the that such duality does most likely hold at high density abrupt changeover at n1/2 in the tensor force in Fig. 2 (and perhaps also at high temperature) [17–20]. In the from the left panel to the right panel. absence of a rigorous proof, we take this as our working To see how this changeover produces the cusp in Esym, assumption. one recalls that the symmetry energy is predominantly Our reasoning relies on two well-known (established) controlled by the tensor force. A quick and simple way facts in nuclear physics in the presence of the HDsF. to estimate the dominant tensor-force contribution to The first is that the symmetry energy is predominantly Esym is to do the closure-sum approximation of the it- controlled by the nuclear tensor force, and the second is erated tensor force terms [26]. This exploits that the that the nuclear tensor force gets principal contributions ground state is strongly coupled by the tensor force (sub- from the exchange of the pseudo-Nambu-Goldstone pion ject, however, to the decreasing strength with density de- π and the ρ meson and coming with an opposite sign, scribed above) to the states of excitation energies ∼ 200 they tend to destructively interfere. It has also been es- MeV, so tablished, within the framework of GnEFT with density- 2 scaling hadron masses [21], that the net tensor force is hVT i Esym ≈ C (12) to decrease with increasing density in the effective range 200 MeV of force in nuclear medium with short-range correlations with C > 0 a known constant. With the NN interactions suitably taken into account. What figures here are the duly screened by short-range correlations (for which the “vector manifestation” of the ρ meson at high density, ω meson enters), it can be seen that hV 2i decreases as the dilaton condensate controlling the hadron masses in T density goes toward n1/2 and then increases afterwards in dense medium and the interplay of the ω-nucleon cou- the precise way as in the skyrmion lattice simulation, thus pling with the nucleon mass [2]. reproducing the cusp Fig. 1 at n1/2. While this argument The resulting tensor force is depicted in Fig. 2. The holds more reliably on the right side of the cusp, namely left panels shows the decreasing tensor force at increasing 6 in the halfs-skrymion phase, it is not expected to to hold density in the absence of topology change . (We note for well in the skyrmion phase away from the cusp. This is later discussion that the net tensor force would vanish because there the effects well described by SχEFT that in the relevant range at n ∼ 3n0). However if there include complicated configurations at high chiral orders intervenes the topology change at n1/2, the tensor force involving other components of the force than the tensor- undergoes a dramatic change as seen Fig. 2 (right panel). force are missing in this treatment. This will become > For n ∼ n1/2, the ρ tensor gets suppressed more or less visible in the GnEFT result shown below. completely so that the net tensor gets abruptly recovered to that of the pionic strength. How this changeover comes about is quite involved re- C. Smoothed Cusp quiring a series of arguments [2], but it is not hard to understand what’a at work in the mechanism with two assumptions. The assumptions are that (A) the vector This calculation for the cusp (with the nuclear ten- mesons introduced as HDsF are hidden local symmet- sor force affected by the topology change) smoothed by ric subject to “composite gauge symmetry” [23, 24] and the HDsF corresponds to the large Nc and quasi-classical (B) the scalar that provides an attractive nuclear force approximation in standard nuclear physics calculations. is the dilaton σ of the “genuine dilaton” structure [25]. In the formulation of GnEFT, this is equivalent to the Now the assumption (A) asserts that at some high den- mean-field approximation with the sHLS Lagrangian [2] sity, the vector mesons become massless, in particular which corresponds to the Landau Fermi-liquid fixed point approximation in the large Nc and large N¯ ≡ kF /(Λ−kF ) with the gauge coupling gρ going to zero [23] associ- ated with the vector manifestation fixed point mentioned limit [27]. As shown in Fig. 3 (right panel) this cusp above, and the assumption (B) admits a (precocious) is made to smoothly cross over in the “Vlowk renormal- emergence of spontaneously broken scale symmetry with ized group (RG) approach” going beyond the Fermi- liquid fixed point approximation in GnEFT employed in [2, 13]. It takes into account 1/N¯ corrections in the “ring-diagram” approximation. It corresponds to a gen- 6 This dropping of the tensor force at increasing density is man- eralized Fermi-liquid theory applicable to the relevant ifested in various nuclear phenomena, the validity of which has range of densities from n0 to the compact-star matter been amply supported. A most spectacular case is the simple and density ∼ (5 − 7)n with the topology change incorpo- elegant explanation of the long lifetime of C-14 beta-decay [22]. 0 There are some ab initio calculations using three-body forces rated at n1/2. It is strictly valid in the large Nc limit that seem to explain more or less equally well, but this should but has been verified to work well for nuclear matter, not be considered belittling the beauty of the simple tensor-force arguably as well as the SχEFT to N3LO. The power of mechanism. Correctly done, both are equally correct in physics. this approach is that while the SχEFT most likely breaks 6

20 20

15 10 10 D D

MeV 5 MeV 0 @ @ T T Ž Ž Þ Þ 0 -10 -5

-10 -20 0 1 2 3 4 0 1 2 3 4 r @fmD r @fmD

FIG. 2. Tensor force vs. density n = (1 − 3)n0: Without topology change (left panel) and with topology change at n1/2 ≈ 2n0 (right panel)

down at n1/2, it becomes more reliable at higher densi- • (A) Up to the topology change density n1/2 >∼ 2n0, ties as the Fermi-liquid fixed point is approached, that is there is only one parameter that is completely as 1/N¯ → 0. The role sHLS plays here in smoothing the determined by how the pion decay constant fπ cusp is analogous to eliminating the cusp singularity in scales with density; it is known up to n ∼ n0. the η0 potential term for the η0 EFT with the HLS fields Within a bit of fine-tuning on this scale parame- becoming topological Chern-Simons fields, giving rise to ter, all the EoS properties come out fully consis- the fractional quantum Hall droplet baryon [18, 20]. Al- tent with the accepted values : They are n0 = −3 though the cusp is smoothed, it makes the symmetry 0.16 fm , BE =16.7 MeV, K0 = 250(240 ± 20) energy that is soft below n1/2 to stiffen above n1/2. This MeV, J = Esym(n0) = 30.2(31.7 ± 3.2) MeV, not only accounts for the observed massive neutron stars L = 67.8(58.7 ± 28.1) MeV. Given in the paren- but as we will show, will render moot the nLORE, hence thesis are quoted – for the illustrative purpose – the possible PREX-II dilemma. What’s even more strik- from the recent compilation by Zhang and Li [28]. ing is that the cusp impacts via the Esym so constructed The same analysis gives the comparison at n = 2n0: the sound speed vpcs as Esym(2n0) = 56.4(50.55 ± 5.99) MeV. This will be relevant for our argument given below. 2 2 Esym(n) → vpcs(n)/c → 1/3 for n >∼ n1/2. (13) • (B) For n > n1/2 ≈ 2n0, there are effectively two This is because Esym(n) = E(n, α = 1) − E(n, α = 0) additional scaling parameters, one for the coupling from Eq. (3) and the trace of the energy-momentum ten- constant gA and the other for the ω-meson gauge sor (TEMT) given by E(n, α) is density-independent for coupling which differs from the ρ gauge coupling α = 0 and α = 1 for n >∼ n1/2 [2]. Hence the crossover that flows to the vector manifestation fixed point in the latter directly impacts the bump in v . These pcs gρ = 0. Both are intricately correlated with the matters will be taken up in Sect. IV B. emergent scale symmetry [2]. This does not affect what follows below, so we won’t go into details here. 7 III. HOW GnEFT FARES The predicted star properties are : Maximum star mass 2.05 >∼ Mmax/M >∼ 2.23 for 2.0 <∼ n1/2/n0 < 4.0, radius R ≈ 12.0 km, Λ < 650, R ≈ In order to give credence to the E obtained in 2.05M 1.4 ∼ 1.4 sym 12.8 km. GnEFT that we will rely on, we summarize what the GnEFT treated in VlowkRG predicts for the EoS for nu- clear matter and how it fares in nature. As stressed in [2], the possible topology change density is constrained < to the range 2 ∼ n1/2/n0 < 4. For simplicity we take 7 There is a possible caveat in what is quoted as our predictions for n1/2 >∼ 2n0 as representing our prediction within a small the relation between the radii R and masses M, particularly for range of uncertainty. GW data. It is argued [29] that to make a quantitatively reliable In what follows, the strangeness flavor degrees of free- analysis on R vs. M, the EoSs of the core and crust should be dom, hyperons as well as kaons, do not enter in the den- treated thermodynamically consistently. This consistency has not been imposed in [13] from which we are quoting the predicted sity regime involved, say, n ∼ (5 − 6)n . The reason for 0 values where the crust-core transition was taken at ncore−crust ≈ this is explained in Conclusion Section. 0.5n0. This caveat might be relevant to the quantities Λ1.4 and We divide the density regime into two: (A) n >∼ n0 and R1.4 but most likely less for other macrophysical quantities of massive stars. (B) n >∼ n1/2 = 2n0. 7

200 0.6

0.5 150 0.4 2 / c

( MeV ) 100 0.3 2 pcm v sym E 0.2 50 0.1

0 0.0 0 1 2 3 4 5 0 1 2 3 4 5

n/n0 n/n0

2 2 FIG. 3. Esym (in MeV) (left panel) and vpcs/c (right panel) for neutron matter (α = 1) calculated in GnEFT for n1/2 ≈ 2.5n0. 0 The bump samples the density range between n and n1/2 in Fig. 1 (right panel). Note that vpcs follows directly and entirely from Esym for n >∼ n1/2 as explained in the text.

IV. HEAVY DEGREES OF FREEDOM AS (right lane). It is the solid line in Fig. 3 (left lane) that DUAL TO GLUONS AND QUARKS: connects the numerically obtained VlowkRG “data” that HADRON-QUARK CONTINUITY 2 2 precisely gives the sound speed vpcs/c that converges to 1/3 at n ∼ 3n0 >∼ n1/2 and stays at 1/3 beyond the The issue of possible resolution to the PREX-II central density ∼ 5n0 of the star. 2 2 dilemma in our approach which is closely linked to also It should be noted that vpcs/c = 1/3 here does not other issues currently in discussion in the literature is en- represent the conformal sound speed associated with capsulated in Esym in Fig. 3 (left panel). It is given by the the vanishing trace of the energy-momentum tensor generalized GnEFT that involves only one Lagrangian (TEMT). It is not the “conformal sound-speed bound” with the HDsF suitably incorporated together with the that is referred to in the literature in addressing the role topology change. It contains no phase transitions in the of “deconfined quarks” in the core of dense neutron stars. sense of Ginzburg-Landau-Wilsonian paradigm, but it is In the system we are dealing with here, the TEMT cannot taken to simulate hadron-quark/gluon continuity. Here go to zero in the density regime involved since it is still we are quoting the result obtained for the crossover den- far from the (putative) IR fixed point [25]. We under- 2 2 sity n1/2 ≈ 2.5n0. For the semi-quantitative aspect we line here that vpcs/c = 1/3 reflects pseudo-conformality are addressing here, the conclusion we arrive at is essen- emergent from strong correlations involving the degrees tially the same for the range 2 <∼ n1/2/n0 < 4. of freedom including the HDsF that give rise to nearly non-interacting quasi-fermions [2]. It embodies hidden scale symmetry. It is far from the state of nearly free “de- A. Heavy Degrees of Freedom and Correlated confined” quarks discussed in the literature where first- Fermions order phase transitions are invoked. It depicts a strongly correlated matter in a way resembling what takes place In addressing the issues involved, there are two impor- in certain condensed matter physics. tant points to note: First, the heavy-degrees of freedom smoothen (or do away with) the cusp “singularity” with the vector mesons B. Bumps of Sound Speed playing the (dual) role of the gluons and induce the crossover from soft-to-hard in the EoS at the transition There are a great deal of discussions currently in the region. As mentioned, the maximum that can be reached literature on the impact on the EoS of dense baryonic in GnEFT is Mmax ≈ 2.23M . At this crossover density, matter in terms of the structure of the sound speed in the however, the maximum of the bump/spike in the sound crossover region in density >∼ 2n0. From the point of view speed exceeds the causality bound at n ≥ 4n0, so may not of nuclear physics, the problem here as mentioned above be physically acceptable although no other star proper- is that from the crossover region indicated between n0 ties seem to go haywire. This implies that our approach and >∼ n1/2 in Fig. 1 for hadrons-to-quarks/gluons to the will get into tension with >∼ 2.5M stars should they be core density of massive stars is the density regime which confirmed to be stable neutron stars. is the hardest to access theoretically, both bottom-up and Second, totally distinctive from all currently available top-down in density. Chiral effective field theory SχEFT ones in the literature, the present EoS unambiguously works reliably for nuclear matter properties, but it is very predicts [2] what we call “pseudo-conformal sound speed” likely to break down at this crossover region. Top-down, 2 2 vpcs/c ≈ 1/3 for density n >∼ n1/2 depicted in Fig. 3 the perturbative QCD must also break down at the star- 8

FIG. 4. Plethora of bumps, spikes, skates and what not in the sound speed cs ≡ vs: (left panel) Wilderness for massive star with M > 2.5M with and without phase transitions assuming it is a stable neutron star instead of a black hole [30]; (right panel) strongly interacting baryonic matter in the core [31]. matter density and will certainly be inapplicable at the the conclusion that the conformal sound speed bound is crossover region. Thus it is not absurd to come up with strongly violated as depicted in Fig. 4(right panel) reach- various wild effects taking place in the region involved. ing the maximum In fact, much discussed are the plethora of bumps, spikes, kinks etc. of various sizes ranging from the crossover to the core density region of star with or with- 2 2 +0.21 out phase transitions. vs /c = 0.79−0.20. (14) Two illustrative cases are given in Fig. 4. The left panel shows the possibilities of massive stars M > 2.5M with bumps of the sound speed all the way from zero to violating the causality bound, typically in- This strong deviation from the conformal sound speed volving phase changes [30]. Our approach, as mentioned is attributed by the authors to “strongly interacting above, cannot access this mass star within the framework hadronic degrees of freedom” that the authors interpret we are working with. It will require a major revamping as “disfavoring” the appearance of “explicit” QCD de- to accommodate such massive stars. grees of freedom in the core of stars. This property is On the contrary, the right panel illustrates the case consistent with the low central density ∼ 3n0 found in without phase changes (or with smooth crossover) that the analysis. The PREX-II dilemma (1) would belong to displays the sound speed largely violating the confor- this class of scenario. mal bound v2/c2 = 1/3 starting from the crossover re- s It should be admitted that given the total paucity of gion [31]. Since the issue of the PREX-II is related to theoretical tools applicable in that density regime, per- what’s treated in [31], this case is highly relevant. The haps one cannot rule out this possibility. analysis of [31] relies on what is called “non-parametric model based on Gaussian processes, un-tied to specific However what GnEFT has predicted is strongly at nuclear models, not subject to systematic errors and pos- odds with the conclusion of [31]. As summarized in sesses wide-range intra-density correlations and targets [2], our pseudo-conformal structure yields the following: wide-range of densities.” While it is not clear to us what Mmax ≈ 2.24M , R1.4 ≈ 12.8 km, ∆R = R2.0 − R1.4 ≈ this model means with respect to our approach, there is a −0.08 km and ncent ≈ 5.1n0. Thus apart from the central striking difference between the two. It is in the structure density which signals the structure of the constituents of of the constituents in the core of massive stars. the core, the overall macrophysical properties predicted For comparison with our prediction, we make a are essentially the same as those of [31]. But the sound list of some of the results reached by the analy- speed of the star (see Fig. 3 (right panel) and Fig. 4 (right sis [31] on the most massive star known so far, panel)) is drastically different. The closeness, but not i.e., J0740+6620 (NICER+XMM-Newton): Mmax = equal, of the pseudo-conformal speed to the conformal +0.34 +1.01 2.24−1.06M , R1.4 = 12.54−1.06 km, ∆R = R2.0 −R1.4 = sound speed bound is consistent with the higher central +0.81 +1.6 −0.04−0.83 km and ncent = 3.0−1.6n0. Based on these density predicted, signaling fractionalized quasi-fermions and other considerations, the authors of [31] arrive at different from baryonic matter. 9

C. Topology Encodes Microscopic Dynamics and/or the dilaton limit fixed point is approached [2], so the scale symmetry should be broken (as indicated by the The principal advantage of our approach on the con- effective dilaton mass which must be substantial coun- trary8 is that it relies on (potentially robust) topological terbalancing the ω mass). However there is an indica- structure which provides a coarse-grained macroscopic tion that scale symmetry can be “emergent,” even if not description of what is presumably taking place in the den- intrinsic, in certain highly correlated nuclear dynamics. sity regime more or less uncontrolled by theory9. In our One prominent evidence for it was seen in the so-called approach the sound speed does produce a simple bump “quenched gA” phenomenon in nuclear beta decay [33]. structure of Fig. 3 (right lane) caused by the interven- The effective gA in nuclear Gamow-Teller transition in ∗ tion of the HDsF dual to QCD degrees of freedom [2] light nuclei gA ≈ 1 arises due to strong nuclear corre- with its characteristics capturing the crossover density lations influenced by hidden scale symmetry reflected in low-energy theorems. The approach g → g = 1 at n1/2. For the case of n1/2 ' 2.0(3.0)n0, it is a bump V A 2 2 high density > 25n0 as the dilaton limit fixed point is reaching vs /c ∼ 0.7(0.8). But for n1/2 ∼ 4n0, as men- ∼ tioned, it is a spike with the maximum of which going approached is closely correlated to how the quenched gA out of the causality bound. Yet despite the different results in finite nuclei [34]. Furthermore that the simple sound-speed structure directly governed by the Esym(n) bump heights in the range 2 < n /n0 < 4 (even vi- ∼ 1/2 ∼ setting in slightly above the crossover density with none olating causality bound for the case n = 4n0), the 1/2 of the compact-star properties significantly disagreeing sound speed v2 /c2 converges in all cases to 1/3 slightly pcs with observations is another indication that the hidden above n with the global star properties not notice- 1/2 scale symmetry is manifested in the density regime of ably different between them. The positive aspect of compact stars. our prediction [2] is that it is an extremely economical description – coarse-graining the microscopic details of what’s found in [30, 31] – that could be capturing the underlying physics. We are arguing that it is precisely D. The PREX-II “Dilemma” and Hadron-Quark the correlated strong interactions leading to the Landau Duality Fermi-liquid quasiparticle structure which becomes more valid with increasing density after the topology change We are now equipped with what’s needed to address as N¯ = kF /(ΛF − kF ) → ∞ [27], manifesting pseudo- the PREX-II dilemma and the issue of whether in EFT conformal symmetry in the sound speed. the EoS at low density near n0 must necessarily con- The question that can be raised is how can the physics strain what happens at higher densities relevant to mas- of the complexity in the sound speed favored by [31] be sive stars. reproduced by the extremely simple structure driven by One can read off from the HDsF-driven Esym ( Fig. 3) the emerging hidden scale symmetry that leads to Fig. 3 that J ≈ 30.2 MeV and L ≈ 67.8 MeV. J is “soft” consis- (right lane)? The possible answer to this could be that tent with (2) but the slope L comes higher than the cen- the macroscopic properties of massive neutron stars are tral value of (2) by ∼ 10 MeV. Reliable SχEFT calcula- in some way insensitive to the microscopic details of the tions to N3LO converge to the central value of ∼ (52−56) bump structure of the sound speed with the emergent MeV [35] which is consistent with (2). symmetries manifesting in the sound speed related to What does this difference of ∼ 10 MeV mean? what’s operative in the “quenching of gA” in baryonic This can be understood as that the “soft” EoS at n <∼ matter mentioned below. If so, the question is: Are there n0 starts to stiffen as the density approaches n1/2 >∼ 2n0 any physically meaningful observable probes for them? reflecting the cusp smoothed by the HDsF. This leads To the best we are aware, there seems to be none at to Esym(2n0) ≈ 56.4 MeV consistent with what is indi- present. cated in nature [28]. This reflects that the SχEFT de- Needless to say, as coarse-grained, there can be fluc- fined with the cutoff Λ <∼ mV starts breaking down at tuations on top of 1/3 coming from corrections to the ∼ 2n0 precisely due to the necessity of the HDsF signal- underlying scale symmetry. That the sound speed con- ing the emergence – vial Seiberg-type duality – of QCD 2 2 verges precociously to vpcs/c ' 1/3 could be an oversim- degrees of freedom. This crossover not only accounts for plification. First of all the density involved < 10n0 is far the massive star masses but also provides a simple mech- from the density at which the vector manifestation limit anism to bring – with additional help with the crust con- sistently treated thermodynamically – Λ1.4 ∼ 650 (which is still consistent with the presently accepted with the upper bound) to a lower value in the vicinity of ∼ 400 8 We believe that what we are discussing here would be little, if which may be favored should the tidal deformability any, affected by the caveat associated with the crust. bound be further tightened in the future measurements. 9 The crucial role of topology played here has an analogy in con- densed matter physics. For instance in the fractional quantum Again there is a logically simple reason for this. The Hall effects Chern-Simons topological field theory captures the heavy-meson-induced smoothing tends to locate the cen- microscopic structure of, say, Kohn-Sham density functional the- tral density of the ∼ 1.4M star for Λ1.4 in the density ory as pointed out (with condensed-matter references) in [32]. regime < n1/2, i.e., the lower side of the cusp – which 10 is soft – that could be in principle accurately calculated argue that the so-called “hyperon problem” is absent in by SχEFTs (and GnEFT) by fine-tuning the crossover the density regime relevant to the core of massive neutron density n1/2 within the range involved. stars. The argument was based on the RG approach to We are then led to suggest that the “strong Rskin-L in interacting protons and neutrons coupled to the HDsF the PREX-II measurement could not constrain the EoS and the kaons on the Fermi surface [42]. Invoking the kF of the core of compact stars. There is, in fact, nothing same large Nc → ∞ and large N¯ = → ∞ limit Λ−kF special about arriving at this sort of conclusion in effec- that underlies the results obtained in this paper (and tive field theories for QCD, the presumed fundamental more generally [2]), it was shown [43] that (1) kaons con- theory of strong interaction physics. An apt example, dense and hyperons appear at the same density and (2) perhaps not widely recognized in nuclear physics com- the kaon condensation threshold density nK satisfies the munity, is the applicability of the skyrmion approach – bound as an EFT – to nuclear physics. Given that the skyrmion theory should be a good low-energy effective field theory nK > Nn¯ 0. (15) of QCD at large Nc limit, it should in principle describe various different low-energy properties of nuclear physics This implies that nK should be considerably higher than valid at large Nc. Indeed in some cases, it works ex- the core density of the stars ∼ (5−6)n0. We note that the tremely well. For instance, the BPS skyrmion is seen bound (15) with nK > 7n0 was arrived at in a different to give an excellent description of nuclear binding ener- but related consideration – short-range correlations – a gies [36] and radii [37] for a wide range of nuclei from light long time ago [44]. A rigorous justification for this could to heavy nuclei A > 200. But the same BPS Lagrangian be given in the Vlowk RG approach to GnEFT ∈ SU(3)f , by itself does not satisfy the soft-pion theorems, the hall- which unfortunately is not feasible at the moment. mark of current algebra and chiral symmetry. This seems at odds with the general belief that nuclear phenomena are governed by chiral symmetry which in the skyrmion V. CONCLUDING REMARKS theory is encoded in the current algebra term in the La- grangian. But it does not necessarily mean that soft the- Starting with the cusp structure found in a skyrmion- orems naively interpreted must give the constraint in the crystal simulation, with the incorporation of heavy de- domain where the BPS structure is more appropriate. grees of freedom considered to be dual to the glu- It is now understood that the infinite tower of vector ons/quarks in the EoS for dense matter, we have ar- mesons, say, HDsF generalized from what we have been rived at the symmetry energy Esym(n) that contraries discussing, subsumed in the BPS Lagrangian is at work the nLORE (standard nuclear astrophysics lore), hence for the particular nuclear properties concerned [38]. the PREX-II “dilemma,” and gives rise to the pseudo- Furthermore one can write down [39] a skyrmion model conformal sound speed for n > n1/2. This result if con- as a sum of BPS submodels, each of which has its own firmed could bring about a potential paradigm change in characteristics applicable to different regions of scales and nuclear physics. dynamics. How to go from one to others is of course an There are several remarks to make to support this pro- open issue that is to be worked out. It is clear however posal. that it is not necessarily constrained by nLORE. The The first is that the topological cusp structure renders recent discovery of Seiberg-type dualities for sHLS [17– moot the necessity to constrain the high-density EoS for 20] indicates that there may intervene more than the massive compact stars by that fixed at normal nuclear skyrmion-half-skyrmion topology change we have been matter density. This of course does not mean that a exploiting in the phase structure of dense hadronic mat- potentially unified theory with no patching cannot be ter in going to high densities in the core of massive stars, 0 connected from n below to above n1/2. In fact in our a notable current example being the phase where the η approach with a single Lagrangian they are connected ring singularity is exposed, say, in the vicinity of the pu- via topological change in a single Lagrangian. tative IR fixed point density [40]. Such a multiple phase The corollary to the first remark, perhaps equally un- structure involving “hadrons” in place of quarks and glu- orthodox, is that the dense star core populated by the ons could persist all the way to the density where the fractionalized fermionic constituents [7] as precisely pre- hadron-quark continuity breaks down [41]. dicted by Esym in the present formulation [2] renders the 2 2 pseudo-conformal sound speed vpcs/c ≈ 1/3 consistent with the picture of “deconfined” quarks with the poly- E. Strangeness Plays No Role tropic index γ approaching 1 [45]. Note however our sys- tem is not at the IR fixed point, hence the constituents of In what’s treated in this paper and elsewhere, the the core are not deconfined. This revamps the common strangeness degrees of freedom played no role. This seems notion of “deconfined quarks” as the only genuine signal at odds with the genuine dilaton scheme [25] where kaons for QCD degrees of freedom. figure on the same footing as the dilaton. This point is The structure of dense matter that emerges from our discussed in [2]. It turns out however to be feasible to work is that what results from highly nonperturbative 11 correlation of QCD degrees of freedom, quarks and glu- asymptotic density way outside of the range of most mas- ons, at the crossover density regime has a dual descrip- sive compact stars stable against gravitational collapse. tion via topology change in terms of hidden local gauge fields and hidden genuine dilaton scalar field in chiral symmetric environment. We suggest that such a dual description stays viable and applicable up to the density Acknowledgments at which the notion of hadron-quark continuity [41] as well as generalized Fermi-liquid structure of weakly in- The work of YLM was supported in part by the Na- teracting fractionalized quasi-fermions, neither baryonic tional Science Foundation of China (NSFC) under Grant nor quarkonic, break down, possibly taking place only at No. 11875147 and 11475071.

[1] H. K. Lee, B. Y. Park and M. Rho, “Half-skyrmions, [15] B. A. Li, B. J. Cai, W. J. Xie and N. B. Zhang, “Progress tensor forces and symmetry energy in cold dense matter,” in constraining nuclear symmetry energy using neutron Phys. Rev. C 83, 025206 (2011) Erratum: [Phys. Rev. C star observables since GW170817,” Universe 7, no.6, 182 84, 059902 (2011)]. (2021). [2] Y. L. Ma and M. Rho, “Towards the hadron–quark con- [16] Y. L. Ma, M. A. Nowak, M. Rho and I. Zahed, “Baryon as tinuity via a topology change in compact stars,” Prog. a quantum Hall droplet and the Cheshire Cat Principle,” Part. Nucl. Phys. 113, 103791 (2020); “Topology change, Phys. Rev. Lett. 123, 172301 (2019). emergent symmetries and compact-star matter,” AAPPS [17] Z. Komargodski, “Vector mesons and an interpretation Bulletin (2021) 31:16; M. Rho and Y. L. Ma, “Manifes- of Seiberg duality,” JHEP 1102, 019 (2011). tation of hidden symmetries in baryonic matter: From [18] A. Karasik, “Skyrmions, quantum Hall droplets, and one finite nuclei to neutron stars,” Mod. Phys. Lett. A 36, current to rule them all,” SciPost Phys. 9, 008 (2020); no.13, 2130012 (2021). “Vector dominance, one flavored baryons, and QCD do- [3] D. Adhikari et al. [PREX Collaboration], “Accurate main walls from the “hidden” Wess-Zumino term,” Sci- determination of the neutron skin thickness of 208Pb Post Phys. 10, 138 (2021), arXiv:2010.10544 [hep-th]. through parity-violation in electron scattering,” Phys. [19] N. Kan, R. Kitano, S. Yankielowicz and R. Yokokura, Rev. Lett. 126, no. 17, 172502 (2021). “From 3d dualities to hadron physics,” Phys. Rev. D 102, [4] B. T. Reed, F. J. Fattoyev, C. J. Horowitz and no. 12, 125034 (2020). J. Piekarewicz, “Implications of PREX-2 on the equa- [20] R. Kitano and R. Matsudo, “Vector mesons on the wall,” tion of state of neutron-rich matter,” Phys. Rev. Lett. JHEP 2103, 023 (2021) 126, no. 17, 172503 (2021). [21] G. E. Brown and M. Rho, “Scaling effective Lagrangians [5] C. Drischler, R. J. Furnstahl, J. A. Melendez and in a dense medium,” Phys. Rev. Lett. 66, 2720 (1991). D. R. Phillips, “How well do we know the neutron-matter [22] J. W. Holt, G. E. Brown, T. T. S. Kuo, J. D. Holt and equation of state at the densities inside neutron stars? A R. Machleidt, “Shell model description of the C-14 dat- Bayesian approach with correlated uncertainties,” Phys. ing beta decay with Brown-Rho-scaled NN interactions,” Rev. Lett. 125, no. 20, 202702 (2020). Phys. Rev. Lett. 100, 062501 (2008). [6] J. Piekarewicz, “Implications of PREX-2 on the [23] M. Bando, T. Kugo, S. Uehara, K. Yamawaki and electric dipole polarizability of neutron rich nuclei,” T. Yanagida, “Is ρ meson a dynamical gauge boson of arXiv:2105.13452 [nucl-th]. hidden local symmetry?,” Phys. Rev. Lett. 54, 1215 [7] M. Rho, “Fractionalized quasiparticles in dense baryonic (1985); M. Harada and K. Yamawaki, “Hidden local sym- matter,” arXiv:2004.09082 [nucl-th]. metry at loop: A New perspective of composite gauge [8] A. W. Steiner, M. Prakash, J. M. Lattimer and P. J. Ellis, boson and chiral phase transition,” Phys. Rept. 381, 1 “Isospin asymmetry in nuclei and neutron stars,” Phys. (2003). Rept. 411, 325 (2005). [24] M. Suzuki, “Inevitable emergence of composite gauge [9] J. W. Holt, M. Rho and W. Weise, “Chiral symmetry and bosons,” Phys. Rev. D 96, no. 6, 065010 (2017). effective field theories for hadronic, nuclear and stellar [25] R. J. Crewther, “Genuine dilatons in gauge theories,” matter,” Phys. Rept. 621, 2 (2016). Universe 6, no. 7, 96 (2020). [10] V. Kaplunovsky, D. Melnikov and J. Sonnenschein, [26] G. E. Brown and R. Machleidt, “Strength of the rho me- “Holographic baryons and instanton crystal,” in [11]. son coupling to ,” Phys. Rev. C 50, 1731 (1994). [11] The Multifaceted Skyrmion (World Scientific, Singa- [27] R. Shankar, “Renormalization group approach to inter- pore,2017) ed. by Mannque Rho and Ismail Zahed. acting fermions,” Rev. Mod. Phys. 66, 129-192 (1994). [12] B-Y. Park and V. Vento, “Skyrmion approach in finite [28] N. B. Zhang and B. A. Li, “Impacts of NICER’s radius density and temperature,” in [11]. measurement of PSR J0740+6620 on nuclear symmetry [13] W. G. Paeng, T. T. S. Kuo, H. K. Lee, Y. L. Ma and energy at suprasaturation densities,” arXiv:2105.11031 M. Rho, “Scale-invariant hidden local symmetry, topol- [nucl-th]. ogy change, and dense baryonic matter. II.,” Phys. Rev. [29] L. Suleiman, M. Fortin, J. L. Zdunik and P. Haensel, D 96, no. 1, 014031 (2017). “Influence of the crust on the neutron star macrophysi- [14] T. H. R. Skyrme, “A unified field theory of mesons and cal quantities and universal relations,” [arXiv:2106.12845 baryons,” Nucl. Phys. 31, 556 (1962). [astro-ph.HE]]. 12

[30] H. Tan, T. Dore, V. Dexheimer, J. Noronha-Hostler Lett. 111, no.23, 232501 (2013). and N. Yunes, “Extreme matter meets extreme gravity: [37] L. A. Ferreira and L. R. Livramento, “A False vacuum Ultra-heavy neutron stars with crossovers and first-order Skyrme model for nuclear matter,” [arXiv:2106.13335 phase transitions,” arXiv:2106.03890 [astro-ph.HE]. [hep-th]]. [31] I. Legred, K. Chatziioannou, R. Essick, S. Han and [38] C. Naya and P. Sutcliffe, “Skyrmions and clustering in P. Landry, “Impact of the PSR J0740+6620 radius light nuclei,” Phys. Rev. Lett. 121, no.23, 232002 (2018). constraint on the properties of high-density matter,” [39] C. Adam, J. Sanchez-Guillen and A. Wereszczynski, arXiv:2106.05313 [astro-ph.HE]. “BPS submodels of the Skyrme model,” Phys. Lett. B [32] Y. L. Ma and M. Rho, “Mapping topology to nuclear 769, 362-367 (2017). dilaton-HLS effective field theory for dense baryonic mat- [40] Y. L. Ma and M. Rho, “Dichotomy of baryons as quan- ter,” arXiv:2103.01860 [nucl-th]. tum Hall droplets and skyrmions In compact-star mat- [33] Y. L. Ma and M. Rho, “Quenched gA in nuclei and emer- ter,” [arXiv:2009.09219 [nucl-th]]. gent scale symmetry in baryonic matter,” Phys. Rev. [41] A. Cherman, T. Jacobson, S. Sen and L. G. Yaffe, “Higgs- Lett. 125, no. 14, 142501 (2020). confinement phase transitions with fundamental repre- [34] M. Rho, “Multifarious roles of hidden chiral-scale sym- sentation matter,” Phys. Rev. D 102, 105021 (2020). metry: Quenching” gA in nuclei,” arXiv:2012.07418 [42] W. G. Paeng and M. Rho, “Kaon condensation in bary- [nucl-th]. onic Fermi liquid at high density,” Phys. Rev. C 91, no.1, [35] R. Essick, P. Landry, A. Schwenk and I. Tews, “A De- 015801 (2015). tailed examination of astrophysical constraints on the [43] M. Rho, “Why explicit strangeness is not relevant in com- symmetry energy and the neutron skin of 208Pb with min- pact stars,” [arXiv:1712.06284 [nucl-th]]. imal modeling assumptions,” [arXiv:2107.05528 [nucl- [44] V. R. Pandharipande, C. J. Pethick and V. Thorsson, th]]. “Kaon energies in dense matter,” Phys. Rev. Lett. 75, [36] C. Adam, C. Naya, J. Sanchez-Guillen and 4567-4570 (1995). A. Wereszczynski, “Bogomolnyi-Prasad-Sommerfield [45] E. Annala, T. Gorda, A. Kurkela, J. N¨attil¨a and Skyrme model and nuclear binding energies,” Phys. Rev. A. Vuorinen, “Evidence for quark-matter cores in mas- sive neutron stars,” Nature Phys. 16, no. 9, 907 (2020).