Topology change, emergent symmetries and compact star matter

Topology effects have being extensively studied and confirmed in strongly correlated condensed matter physics. In the large color number limit of QCD, baryons can be regarded as topological objects -- skyrmions -- and the baryonic matter can be regarded as a skyrmion matter. We review in this paper the generalized effective field theory for dense compact-star matter constructed with the robust inputs obtained from the skyrmion approach to dense nuclear matter, relying to possible ``emergent"scale and local flavor symmetries at high density. All nuclear matter properties from the saturation density $n_0$ up to several times $n_0$ can be fairly well described. A uniquely novel -- and unorthdox -- feature of this theory is the precocious appearance of the pseudo-conformal sound velocity $v^2_{s}/c^2 \approx 1/3$, with the non-vanishing trace of the energy momentum tensor of the system. The topology change encoded in the density scaling of low energy constants is interpreted as the quark-hadron continuity in the sense of Cheshire Cat Principle (CCP) at density $\gsim 2n_0$ in accessing massive compact stars. We confront the approach with the data from GW170817 and GW190425.


I. INTRODUCTION
The structure of dense nuclear matter relevant to compact stars has been investigated for several decades but still remains largely uncharted. Unlike at high temperature, so far, the physics at high density can be accessed by neither terrestrial experiments nor lattice simulation. Recently, the observation of massive neutron stars with mass ∼ > 2.0M and detection of gravitational waves from neutron star mergers provide indirect information of nuclear matter at low temperature and high density, say, up to ∼ 10 times the normal nuclear matter density n 0 0.16 fm −3 [1][2][3][4][5][6]. These new developments offer the powerful means to explore the nuclear matter in the interior of the compact stars, for example, the patterns of the symmetries involved therein, what is in the core of the stars, say, baryons and/or quarks and a combination thereof. For recent discussions on these aspects, we suggest, e.g., Refs. [7][8][9][10][11][12] and some relevant references therein.
The study of nuclear matter in the literature has largely relied on either phenomenological approaches anchored on density functionals or effective field theoretical models implemented with assumed QCD symmetries and degrees of freedom appropriate for the cutoff to which the theory is applicable. For finite nuclei as well as the infinite nuclear matter up to ∼ n 0 , the physics can be described very well by using the nuclear effective theory with or without pion, in addition to the nucleon [7,13] (denoted as sχEFT). However, in the dense system relevant to the compact stars at ∼ 10n 0 , the sχEFT is symmetries are not explicit in the matter free space but it seems reasonable to think that they get (partially) restored in the dense system. At least there is nothing glaring at odds with the presently available observations. Prior to QCD, Skyrme suggested that baryons can be described by the topology solution of a mesonic theory, skyrmion [20]. After the arrival of QCD, it was argued that when the number of color N c is infinitely large, baryons in the constituent quark model share the same N c scaling properties as skyrmions [21,22]. Since then, the Skyrme(-type) model 1 has become one of the models in the study of nucleon, nuclei as well as nuclear matter [23][24][25].
In the skyrmion approach to dense nuclear matter obtained by putting skyrmions on crystal lattice, a robust observation independent of the model and crystal structure-at least what has been checked so far-is the topology change where the skyrmions with the integer winding number transit to half-skyrmions with the halfinteger winding number. The density at which this takes place is denoted as n 1/2 . This model-independent topology change gives rise to several interesting density dependences of hadron properties that have not been found in other approaches.
Although the Skyrme model approach can describe the nucleon, nuclei as well as nuclear matter in a unified way, it is a daunting task to put this approach into practice since the calculation depends on the efficiency of the computer and the results are valid in the large N c limit. Therefore, in practice, one resorts to chiral effective models that incorporate baryons as explicit degrees of freedom. In our GnEFT, we incorporate the robust characteristics of topology in the low energy constants of the model. The effect of the change of the degrees of freedom is formulated in terms of the possible topology change at a density n 1/2 encoded in the behavior of the parameters of the GnEFT Lagrangian as one moves from below to above the changeover density n 1/2 . After making the V lowk RG approach implementing the strategy of Wilsonian renormalization group flow [26], we construct the pseudo-conformal model (PCM) of dense nuclear matter [27,28](see Ref. [10] for a review).
The PCM that satisfy all the constraints from astrophysics turns out to have a peculiar feature that has not been found in any other approaches: The sound velocity approaches the conformal limit v 2 s /c 2 ≈ 1/3 at the density relevant to compact stars although the trace of the energy-momentum tensor does not vanish. This is in stark contrast to the standard scenario favored in the field [29]. This conceptually novel approach predicts that the core of massive compact stars is populated by confined quasi-fermions of fractional baryon charge [30], not "deconfined quarks" expected in perturbative QCD [31]. We suggest that this phenomenon, together with the 1 Hereafter, for convenience, we use Skyrme model with the pion field only to represent the Skyrme model and its extensions.
"quenched g A problem" in nuclei, shows that hidden symmetries hidden in medium-free vacuum of QCD emerge in nuclear dynamics [32].

II. TOPOLOGY CHANGE AND HADRON-QUARK CONTINUITY
It has long been discussed that in the large number of color N c limit, baryons can be regarded as topological objects-solitons, namely skyrmions. In the skyrmion approach, the dense nuclear matter can be accessed by putting the skyrmions onto the crystal lattice [25,33,34].
Here we exploit the Skyrme model with the Lagrangian connected to QCD in the sense of Weinberg "folk theorem" on effective field theories [35]. For a development quite different in spirit from ours, we refer to, e.g., review [36] and the references therein.

A. Topology change
Topology change is a novel phenomenon that has not been observed in any approach other than the skyrmion crystal approach to dense nuclear matter.
To have an intuitive idea, let us look at the distribution of the baryon number density in a specific lattice, say, face-centered cubic crystal. The distribution of the baryon number density looked along an axis is illustrated in the left panel of Fig. 1. The winding number is 1 if one integrates out the blue volume. Now, squeeze the system. One finds that, after a critical density n 1/2 (or equally, the crystal size L 1/2 ), the distribution of the baryon number density changes to the right panel of Fig. 1. What happens is that when the increasing matter density surpasses n 1/2 (or the crystal size drops below L 1/2 ), the constituents of the matter given in the blue square transit from winding number-1 objects (left panel) to winding number-1/2 objects, half-skyrmions (right panel). (How this happens in the numerical simulation can be seen in [34] .) In the half-skyrmion configuration, as a consequence of symmetry, the space average of where U 0 is the static configuration of the chiral field U = exp(2iπ a T a /f π ) with T a = σ a /2. This means that the quark condensate qq vanishes when the space is averaged. Therefore, one can use this quantity as a signal of the skyrmion-half-skyrmion transition.
It should be noted that the location of n 1/2 cannot be pined down theoretically because it is model-dependent. Since nuclear dynamics at low density can be well described by sχEFT, we set n 1/2 ∼ > 2.0n 0 . Later, we will see that astrophysical observations observations indicate 2.0n 0 ∼ < n 1/2 ∼ < 4.0n 0 .
B. Implications of topology change Chiral symmetry breaking.-In the skyrmion crystal approach to dense nuclear matter, the pion decay constant can be calculated through the axial-vector currelator [37] iG ab At the leading order of fluctuations, we can express the medium modified pion decay constant as In the skyrmion phase, since φ 2 0 decreases with density, f * 2 π decreases with density. After passing n 1/2 from below, since φ 0 = 0 in the chiral limit, φ 2 0 0. Thus a nonzero constant although φ 0 = 0. This argument is supported by explicit numerical calculation. In terms of current algebra, the generalized Gell-Mann-Oakes-Renner relation tells us [38] where for convenience, we have kept the current quark mass. F n stands for the contribution from multiquark condensation. Since the pion mass scales little with density, when going to the half-skyrmion matter, which is in qualitative agreement with the result from skyrmion crystal calculation. Equation (4) means that the chiral symmetry is only partially restored in the half-skyrmion matter and we are still in the Nambu-Goldstone phase. This means that the skyrmion-half-skyrmion transition is not a Landau-Ginzburg-type phase transition. Although it is not a paradigmatic phase change, in what follows, we will use the term "half-skyrmion phase" for simplicity.
Chiral doublet structure.-It is found that when the system goes to the half-skyrmion medium, the nucleon mass becomes a density-independent constant [39]. Therefore, one can decompose the nucleon mass as where ∆( qq ) is the sector of the nucleon mass coming from the quark condensate which becomes zero in the half-skyrmion medium. m 0 is the sector of the nucleon mass independent of qq and has a magnitude about (50−70)% of the nucleon mass in vacuum. The existence of m 0 = 0 implies that there is a part of the nucleon mass that is chiral invariant. It should be noted that, the decomposition (7) can also be inferred from other approaches. The lattice calculation found that, when the chiral symmetry is unbroken, baryons are still massive and one should not expect a drop of the mass in dense medium [40]. The same behavior was found in Ref. [41] in a renormalization group (RG) analysis of hidden local symmetric Lagrangian with baryons. Moreover, in Ref. [42], by using a chiral effective model with parity doubler, it was found that, to reproduce the nuclear matter around saturation density, the nucleon mass should has a sizable chiral invariant component. So far, it is not clear to us whether m 0 reflects a fundamental feature of QCD or an emergent symmetry via correlations in medium as in condensed matter as indicated in this crystal calculation.
Symmetry energy.-The symmetry energy of nuclear matter E sym (n) which plays the most important role in the equation of state (EoS) for compact stars is not under control at the density relevant to compact stars [43,44]. It is given by the term proportional to α 2 in the energy per nucleon E(n, α) E(n, α) = E(n, α = 0) + E sym (n)α 2 + O(α 4 ), (8) where α = (N −P )/(N +P ) with P (N ) being the number of protons (neutrons).
Since the symmetry energy arises from the protonneutron asymmetry, to calculate it from the skyrmion crystal approach, the crystal lattice should be rotated through a single set of collective coordinates [45]. A tedious but straightforward calculation yields where λ I is the isospin moment of inertia. The density dependence of the symmetry energy obtained from the skyrmion crystal approach is schematically plotted in dotted curve in Fig. 2. What is interesting is the appearance of the cusp structure locked at n 1/2 , i.e., the symmetry energy first decreases with density and then increases when the density passes n 1/2 . To understand the density dependence of symmetry energy, we consider the expression of λ I [46] λ I = f 2 where · · · stands for the contribution from the Skyrme term and · · · indicates the space average of the quantity inside. As discussed above, with the increasing of the density, φ 2 0 decreases to zero. So, 1/λ I , or equivalently E sym , decreases going toward n 1/2 . After n 1/2 , the tendency of E sym is highly involved. Since at n ∼ > n 1/2 , φ 2 0 ≈ 0, the density dependence from the quartic term in the Skyrme model which represents massive excitations -such as the vector mesons in the HLS modelsintervene. It gives the cusp structure.
It should be stressed that the crystal description of baryonic matter at low density cannot be reliable, so the density dependence of symmetry energy obtained at density n n 0 cannot be taken seriously. The cusp structure at n 1/2 is present in nuclear correlations as is shown below in terms of nuclear tensor forces. What is important in the skyrmion crystal calculation is that the symmetry energy decreases toward the cusp density after which it increases. We will see later that this cusp sheds light on the medium modified-hadron properties.
Nuclear tensor force.-We have shown that, the robust characteristic in the skyrmion crystal approach is the existence of the cusp structure in the symmetry energy. A natural question is what is the implication of this cusp in GnEFT including nucleon as an explicit degree of freedom or equivalently, how to reproduce this cusp in GnEFT. To address this question, we consider the tensor force between nucleons that is mediated by one boson exchange.
The symmetry energy is dominated by the nuclear tensor force V T and can be written in the closure approximation as [47] Therefore, the behavior of the symmetry energy is controlled by the absolute value of the tensor force between nucleons carried by the exchanged mesons. For the present purpose, it suffices to consider the onepion and one-ρ contributions to two-body tensor forces. The scalar meson, here dilaton, does not contribute directly at the tree level but affects indirectly on the scaling relations of the masses and coupling constants in the Lagrangian. In the non-relativistic limit, the tensor forces are given by where M = π, ρ, S ρ(π) = +1(−1) and with the Pauli matrices τ i and σ i for the isospin and spin of the nucleons with i = 1, 2, 3. The density dependence enters through the scaling parameters in the in-medium quantities marked with asterisk [48]. The strength f * N M scales as where g M N N are the effective meson-nucleon couplings. What is significant in Eq. (12) is that given the same radial dependence, the two forces (through the pion and ρ meson exchanges) come with an opposite sign and therefore cancels each other.
As discussed in Ref. [49], if the hadron scales from low to high densities with no topology change, the tensor force will decrease monotonically with density. There will then be no cusp in the symmetry energy. This feature will be in conflict with what happens in Nature.

Now let us see what happens
if there is the topology change at n 1/2 . For illustration we take R * ρ ≈ Φ 2 at n > n 1/2 but with all others the same as in the case without topology change. The results are plotted in Fig. 3. It shows that the topology change effect is dramatic. Due to the cancellation between these two tensor forces, in the range of nuclear forces relevant for the nuclear interaction, r ∼ > 1 fm, the magnitude of the net force first decreases and then, after passing n 1/2 , increases and the force from the rho meson is nearly totally suppressed. Then, from Eq. (11), one concludes that going toward to n 1/2 from below the symmetry energy is to drop and more or less abruptly turn over at n 1/2 and then increase beyond n 1/2 . This reproduces precisely the cusp predicted in the crystal calculation. As a result, the cusp structure in E sym -a consequence of topology change with the onset of the half-skyrmion phase-is signaling the different density scaling property of the gauge coupling from n ≤ n 1/2 to n > n 1/2 .
In summary, the topology change found in the skyrmion crystal approach to density nuclear matter indicates that the hadron properties, such as nucleon mass, meson masses, pion decay constant and hidden gauge coupling and so on, have different density scaling in the skyrmion and half-skyrmion phases. We will see later that, this observation has a drastic effect on the dense nuclear matter for compact stars.
We should mention here that higher correlation corrections brought in the V lowk renormalization flow calculation "smoothen" the cusp in the form represented in solid line in Fig. 2.

C. Quark-hadron continuity
We have argued that the topology change is a robust feature in the skyrmion crystal approach to dense nuclear matter. The question is whether or how the topology change represents the "quark deconfinement" process in QCD. There is no clear answer at present, so we can only offer a conjecture on how one can establish the connection in the sense of Cheshire Cat Principle (CCP) based on the chiral bag model of nucleon.
For the number of flavors N f ≥ 2, baryons can be described by chiral bags [14,15]. Inside the bag, the degrees of freedom are quarks and gluons, and the baryon number is carried by the quarks. Outside of the bag, mesons are the relevant degrees of freedom, and the baryon number is carried by topology in the winding number. When the bag is shrunk, all the quarks drop into the inifinite hotel and turn into skyrmions with only the Cheshire Cat smile remaining. That physics should not depend on the bag size is the CCP.
In the case of single flavor N f = 1, the situation is quite different because there is no N f = 1 skyrmion. It turns out that the baryon should be a soliton resembling a pancake [50] or pita [51] having a fractional quantum Hall (FQH) topology structure. There is a Cheshire Cat description for this in terms of an anomaly flow [52]. But what is puzzling is that there are two Cheshire Cats, one involving 3D ball and the other 2D sheets. It seems very plausible that at low density baryonic matter is in skyrmions in 3D with the metastable 2D FQH pan-cakes/pitas suppressed. However it seems indispensable at high density that the FQH topology structure be taken into account. This is because at high density where chiral transition takes place, the vector mesons in hidden local symmetry become the Chern-Simons fields (via Seibergtype duality). This part of the high density story is not yet understood, so we can only say that we really do not understand what happens at high density. In what we have done, we are simply assuming that the Chern-Simons fields do not figure importantly in the range of compact-star densities. We will simply ignore this "dichotomy problem." This aspect of the problem is discussed in [53,54].
The topology change in the skyrmion crystal approach appears at the density at which the profiles of solitons overlap and the valence quarks inside the baryons rearrange to form different clusters, here configurations with baryon number-1/2. This picture resembles the quarkyonic matter proposed in [8,55] and the hard-core realization of the deconfinement from nuclear to quark matter phrased in Ref. [56].
As mentioned above and will be discussed later, owing to the topology change implemented in the parameters in the Lagrangian of GnEFT, the symmetry energy E sym , as it approaches n 1/2 from slightly below, softens and after passing n 1/2 , hardens. This generates a spike in the density dependence of the sound velocity. In Ref. [57], this spike was attributed to the enhancement and then suppression of the ω 0 condensate in the low and high density region. We suppose that this behavior of ω 0 condensate can be naturally explained using the scale-chiral effective theory beyond the leading order scale symmetry in which not only the ω meson mass but also the ω-N -N coupling scales with density [58].

III. EMERGENT SYMMETRIES
After the discussion on the topology change which serves as one of the key ingredient of the PCM, let us now turn to two other essential ingredients, the hidden local gauge symmetry and hidden scale symmetry which are invisible in the vacuum of QCD. Our approach is to exploit the possible emergence of these symmetries as density increases to the regime relevant to compact stars, say, ∼ < 10n 0 . We use these symmetries to include the higher-energy degrees of freedom-the lowest-lying vector mesons V = (ρ, ω) and the scalar meson f 0 (500). Here, we focus on the points directly relevant to the PCM construction, leaving the details to [10,59] A. Emergent hidden local symmetry To bring in the lowest-lying vector mesons ρ and ω into the chiral effective theory, we adopt the strategy of hidden local symmetry (HLS) [16][17][18] which at low density is gauge equivalent to nonlinear sigma model, the basis of sχEFT.
By decomposing the chiral field U (x) as U (x) = ξ † L ξ R , one can introduce a redundant local symmetry h(x) under which ξ L,R transforms as but keeps the chiral properties of U (x) intact. When a chiral effective theory is expressed in terms of ξ L,R , the gauge fields of local symmetry h(x)-V (x)-enter the theory. After higssing the gauge symmetry, the gauge fields V (x) obtain masses. In HLS, the field content depends on the symmetry h(x). If one chooses (2) and ω ∈ U (1). It is assumed that the kinetic terms of V µ (x) can be generated by underlying dynamics of QCD or quantum corrections, thus V µ (x) become dynamical gauge bosons [16]. Compared to other approaches of vector mesons, with HLS, one can establish a systematic power counting by treating the vector mesons on the same footing as the Nambu-Goldstone boson, pions [18]. Now, come back to the nuclear matter. At low density where the nucleons are far from each other, the vector mesons are massive objects and can be exchanged between them. Using the equations of motion of the vector mesons, their effects are accounted for as a two-pion exchange effect, i.e., one-loop contribution in sχEFT. The question is in whst sense the vector mesons can be regarded as hidden local gauge fields. The Suzuki's theorem [60] states that "when a gauge-invariant local field theory is written in terms of matter fields alone, a composite gauge boson or bosons must be inevitably formed dynamically." If we assume the "vector manifestation (VM)" [18,61] that m 2 ρ ∝ f 2 π g 2 ρ → 0 since g ρ → 0 at certain scale valid at some theoretically unknown high density n vm , the hidden local gauge symmetry emerges in dense system. We will see below that n vm ∼ > 25n 0 is indicated for the emergence of the pseudo-conformal sound velocity in stars.
Moreover, it was argued that the HLS fields could be (Seiberg-)dual to the gluons [62][63][64]-the intrinsic quantity in QCD. At this moment, we do not know how could this happen. But, if this is right, we believe it means that the HLS gets un-hidden at high density. Since this duality indicates a Higgs phase-to-topological phase transition coinciding with the quark deconfinement at asymptotic density, it is most likely irrelevant to the compact stars we are concerned with [53] .

B. Emergent scale symmetry
It is well known that the scalar meson f 0 (500) is essential for providing the attractive force between nucleons. In our approach, it figures in as the Nambu-Goldstone boson of scale symmetry, the dilaton χ. Actually, in Ref. [65], the trace anomaly has been applied as a source of the scalar meson to construct the effective model of scalar meson by using anomaly matching. Here, we adopt the "genuine dilaton" (GD) structure proposed in [19].
The key premise of the GD idea is the existence of an infrared fixed point (IRFP) with the beta function β(α IR ) = 0 for flavor number N f = 3 and we are living slightly away from this IRFP. Both the distance from the IRFP and current quark masses account for the dilaton mass. Explicitly, the dilaton mass is expressed as where · · · stands for the contribution from quark mass and higher order of ∆ IR = α IR − α s with α s . This is in analogy to the Gell-Mann-Oakes-Renner relation in the pseudoscalar meson sector. Since, unlike the unflavored hadrons, the effective masses of strange hadrons do not drop so much in dense medium, we will not consider the strangeness here.
Whether the proposed IRFP exists in QCD is still under debate. In Ref. [66], it was argued that in the IR region there is a nonperturbative scale invariance different from that in the UV region. This is argued to lead to the possibility of massless glueballs in the fluid. What may be significant is the possible zero-mass glueball excitation. If we simply assume this picture works in dense system, this can be regarded as an indirect support of our theme. Anyway, we did not find any contradiction with nature in using this GD idea.

IV. PSEUDO-CONFORMAL MODEL OF COMPACT STAR MATTER
By using the GnEFT discussed above we are now in the position to calculate the nuclear matter properties. We shall focus on the EOS of the baryonic matter, leaving out such basic issues as corrections to gravity, dark matters etc. Unless otherwise stated the role of leptons-electrons, muons, neutrinos, etc.,-is included in the EOS. Hereafter, we mainly focus on the effect of topology change. For other aspects, we refer to [10].

A. Density scaling
In the construction of the PCM, we incorporate the medium modified hadron properties (dubbed as "intrinsic density dependence (IDD)" to the GnEFT we constructed above by using the Brown-Rho scaling [48] for n ≤ n 1/2 (R-I) and the topological inputs for n > n 1/2 (R-II).
Density scaling in R-I.-In R-I, only one parameter Φ in Eq. (18) fixes all the IDDs. To the leading order in the chiral-scale counting [67], the density scaling in R-I can be written as [10] where V = (ρ, ω). Since there is no first-principle information on this quantity, for convenience, we fix it by taking the form Φ I = 1 1 + c I n n0 (19) with c I a constant. The range of c I that gives a good fit to nuclear matter properties as shown in Table. I is found to be [27,49] c I ≈ 0.13 − 0.20 (20) with the upper value giving the measured pion decay constant [68]. Of course, it is expected as would be agreed by all nuclear physicists that certain fine-tuning in the parameters be required for ground-state properties of nuclear matter. Density scaling in R-II.-Due to the topology change at n 1/2 > n 0 , the density dependence of some parameters are drastically different from that in R-I.
Since the hidden local gauge coupling g ρ and the ρ meson mass are related to each other through the KSRF relation, we take the simplest form where n VM is the putative VM fixed-point density. How to join the Φ ρ from Φ I for n ≤ n 1/2 is discussed in Ref. [27]. To have a result consistent with that from skyrmion crystal approach discussed above and mean field approach based on the leading order scale symmetry (LOSS) [72], we take n VM ∼ > 25n 0 . The density scaling of the ω meson is more involved and different from that of ρ meson which flows to the VM fixed point [41,49]. It should be fine-tuned to match to the well constrained nuclear matter properties around the saturation density. Here, we take where g ω is the U (1) gauge coupling and with d ≈ 0.05. This reflects the predicted break-down in R-II of the flavor U (2) symmetry for the vector mesons which holds well in R-I. As for other parameters, we simply adopt the inputs from the skyrmion crystal approach, that is The dilaton mass also goes proportional to the dilaton condensate. This follows from the partially conserved dilatation current (PCDC) [19] m * σ m σ ≈ κ.
It follows also from low-energy theorems that The dilaton coupling to nucleon and other fields is unscaling to the leading order in scale-chiral symmetry, so it is a constant in R-II as in R-I. V lowk renormalization group approach.-Equipped with the IDD, we are ready to calculate the EoS of nuclear matter. Here, to take into account the hadron fluctuation effects, we apply the V lowk renormalization group technique [26] which accounts for higher-order corrections to the Landau Fermi-liquid approximations [10]. In this procedure, in addition to the IDD implemented in the density scaling of the parameters, the induced density dependence from the nucleon correlation denoted as DD induced 2 is also included. Therefore, the density dependence in the obtained EoS includes both IDD and DD induced . We denote the sum of IDD and DD induced as IDD.
We would like to point out that, owing to the IDD of the two-nucleon potentials, our calculation amounts to doing roughly N 3 LO SχEFT including chiral 3-body potentials which are essential for the nuclear matter stabilized at the proper equilibrium density [73]. The same mechanism has been found to work for the C-14 dating Gamow-Teller matrix element where the three-body potential effect in SχEFT is reproduced by the IDD.

B. The pseudo-conformal model of dense nuclear matter
Using the density scaling discussed above, we can calculate the nuclear matter properties now. First we see from Table I that the empirical values of the normal nuclear matter properties can be well reproduced. Now, go to a higher density. Due to the topology change at n 1/2 , there is a drastic change in the scaling of the parameters of GnEFT leading to a qualitative impact on the structure of the EoS. So far, there is no theoretical argument to pin down n 1/2 . Phenomenologically, we can estimate its range as 2.0n 0 < n 1/2 < 4.0n 0 by using various astrophysical observations available, such as the maximum mass, the gravity-wave data and specially the star's sound speed, and so on. Sound velocity.-One of the most striking predictions that is in stark contrast to the conventional picture is the precocious appearance of the conformal sound velocity of the compact star matter. From Figure 4, one can see that while the sound speed increases steadily and overshoots the conformal velocity at presumed ∼ n 1/2 = 2n 0 , it comes down and converges to v 2 s ≈ 1/3. It should be noted that the appearance of the conformal sound velocity at some high density is not so peculiar. Some reasonable sχEFT results resemble more or less this picture. But they show much broader and bigger bumps not exceeding the causality bound v s = 1 before converging to the conformal speed v 2 s = 1/3 but at an asymptotic density ∼ > 50n 0 [29]. After all, the convergence to the conformal speed at asymptotically high density is expected in perturbative QCD. What is striking and in a way unorthodox is the precocious onset of, and the convergence to, v 2 s ≈ 1/3 before reaching to an asymptotically high density despite that the trace of the energy-momentum tensor is nonzero. See below. It is somewhat like the "quenched g A " going to 1 in light nuclei [32], reflecting the pervasive imprint of hidden scale symmetry.
In our approach, the conformal sound speed follows as a logical outcome of the propositions [10], different from the parameter scanning done in [74]. These propositions yield that, going toward the DLFP [41], the trace of the energy-momentum tensor θ µ µ is a function of only the dilaton condensate χ * . Now if the condensate goes to a constant ∼ m 0 due to the emergence of parity-doubling as we learned after the topology change, the θ µ µ will become (more or less) independent of density. In this case, we will have ∂ ∂n θ µ µ = 0. (27) This would imply that where v 2 s = ∂P (n) ∂n / ∂ ∂n and and P are, respectively, the energy density and the pressure. If we assume ∂ (n) ∂n = 0, i.e., no Lee-Wick-type states in the range of densities involved, we can then conclude This means that the dilaton condensate χ * goes to the density-independent constant m 0 due to the parity for n vm ∼ > 25n 0 . This suggests the parity doubling at high density is linked to the ρ decoupling from the nucleon together with the vector manifestation [41].
The above chain of reasoning is confirmed in the full V lowk RG formalism specifically for the case of n 1/2 = 2n 0 . In Fig. 5 is shown the trace of the energy momentum tensor (left panel) that gives the conformal velocity for n ∼ > 3n 0 (right panel). This feature of both the TEMT and the sound velocity are expected to hold for any n 1/2 at which the topology change sets in, i.e., within the range 2 ∼ < n 1/2 /n 0 ∼ < 4. Equation of state.-We now focus on the EoS of compact stars. It turns out that at density n ≥ n 1/2 , the conformality of the sound velocity can be captured by a simple two-parameter formula for the energy per-particle where X, Y are parameters to be fixed.
What we refer to as the pseudo-conformal model (PCM for short) for the EoS is then E/A given by the union of that given by V lowk in R-I (n < n 1/2 ) and that given by Eq. (30) in R-II ( n ≥ n 1/2 ) with the parameters X α and Y α fixed by the continuity at n = n 1/2 of the chemical potential and pressure µ I = µ II , P I = P II at n = n 1/2 .
This formulation is found to work very well for both α = 0 and 1 in the entire range of densities appropriate for massive compact stars, say up to n ∼ (6 − 7)n 0 , for the case n 1/2 = 2n 0 where the full V lowk RG calculation is available [27]. We apply this PCM formalism for the cases where n 1/2 > 2n 0 . Since a neutron star with mass 1.4M for which the tidal deformability Λ obtained for n 1/2 = 2.0n 0 is Λ 1.4 790 [27,28] that corresponds to the upper bound set by the gravity-wave data, we take the lower bound for the topology change density FIG. 5. θ µ µ (upper panel) and vs (lower panel) vs. density for α = 0 (nuclear matter) and α = 1 (neutron matter) in V lowk RG for n 1/2 = 2n0.
Next, let's see how the sound velocity comes out for n 1/2 /n 0 = 3 and 4 [28,75]. (The case for n 1/2 = 2n 0 was given in Fig. 5.) The results for neutron matter are summarized in Fig. 6. It is clear from Fig. 6 that, when n 1/2 = 4n 0 , the sound velocity violates the causality bound v 2 s < 1. The spike structure could very well be an artifact of the sharp connection made at the boundary. It may also be the different behaviour of the ω 0 condensation at the low and high densities [57]. What is however physical is the rapid increase of the sound speed at the transition point signaling the changeover of the degrees of freedom. Significantly, together with the lower bound (32), this allows us to pinpoint the region of the topology change Later, we will explore whether or how the waveforms of the gravitational waves emitted from the binary neutron star mergers respond to the location of n 1/2 which in our formulation corresponds to the point of hadron-quark continuity in QCD. At this moment, we cannot obtain a more precise constraint than (33). The important point is that it is an order of magnitude lower than the asymptotic density ∼ > 50n 0 that perturbative QCD predicts and signals the precocious emergence of pseudo-conformality in compact stars. However the robustness of the topological inputs figuring in the formulation convinces us that the precocious onset of the pseudo-conformal structure can be trusted at least qualitatively. In this connection, a recent detailed analysis of currently available data in the quarkyonic model is consistent with the possible onset density of v 2 c ≈ 1/3 at ∼ 4n 0 [76]. Predicted pressure for neutron matter (α = 1) vs density compared with the available experimental bound (shaded) given by Ref. [77] and the bound at 6n0(blue band). Fig. 7 is the predicted pressure P vs. density for n 1/2 /n 0 = 3, 4 compared with the presently available heavy-ion data [77]. The case of n 1/2 = 4n 0 , while consistent with the bound at n ∼ 6n 0 , goes outside of the presently available experimental bound at n ∼ 4n 0 . This may again be an artifact of the sharp matching, but that it violates the causality bound seems to put it in tension with Nature. Nonetheless, without a better understanding of the cusp singularity present in the symmetry energy mentioned above it would be too hasty to rule out the threshold density n 1/2 = 4n 0 .

V. STAR PROPERTIES AND GRAVITATIONAL WAVES
Star mass.-The solution of the TOV equation with the pressures of leptons in beta equilibrium duly taken into account as in Ref. [27] yields the results for the star mass M vs. the radius R and the central density n cent as given in Fig. 8. The maximum mass comes out to be roughly 2.04M ∼ 2.23M for 2.0 n 1/2 /n 0 4.0, the higher the n 1/2 , the greater the maximum mass. This bound is consistent with the observation of the massive neutron stars M = 1.908 ± 0.016M for PSR J1614 − 2230 [1], = 2.01 ± 0.04M for PSR J0348 + 0432 [2], = 2.14 +0.10 −0.09 M for PSR J0740 + 6620 [3]. Note that this is not at odds with the conclusion of Ref. [78] since in our model, the sound velocity exceeds the conformal limit in the intermediate density. Fig. 8 shows that, when n 1/2 ≥ 3.0n 0 , changing the position of n 1/2 affects only the compact stars with mass 2.0M although the mass-radius relation is affected by the topology change when 2.0n 0 ≤ n 1/2 ≤ 3.0n 0 .
Tidal deformability.-Next, we confront our theory with what came out of the LIGO/Virgo gravitational observations-the dimensionless tidal deformability Λ. We will consider the dimensionless tidal deformability Λ i for the star M i andΛ defined bỹ Λ = 16 13 for M 1 and M 2 constrained to the well-measured "chirp mass" GW190425 [6]. (35) We plot our predictions forΛ in Fig. 9 and for Λ 1 vs. Λ 2 in Fig. 10 and compare our predictions with the results obtained with the parametrization of the EoS from the sound velocity constraints [74]. As it stands, our prediction with n 1/2 ∼ > 2n 0 is compatible with the LIGO/Virgo constraint. Although there seems to be some tension with the pressure, the result for n 1/2 = 4n 0 is of quality comparable to that of n 1/2 = 2n 0 . A detailed analysis of the difference between PCM and [74] will be made later.
Massive star composition.-Recently, combining astrophysical observations and model-independent theoretical ab initio calculations, Annala et al. arrive at the conclusion that the core of the massive stars is populated by "deconfined" quarks [31]. This is based on the observation that, in the core of the maximally massive stars, v s approaches the conformal limit v s /c → 1/ √ 3 and the polytropic index takes the value γ < 1.75 -the value close to the minimal one obtained in hadronic models.
We have seen above that, in the PCM, the predicted pseudo-conformal speed sets in precociously at n ≈ 3n 0 and stays constant in the interior of the star. In addition, it is found that the polytropic index γ drops, again rapidly, below 1.75 at ∼ 3n 0 and approaches 1 at n ∼ > 6n 0 [30]. This can be see from Fig. 11. Microscopic descriptions such as the quarkyonic model [80] typically exhibit more complex structures at the putative hadron-quark transition density than our description, which is not unexpected given our picture is coarse- . The the PCM prediction with n 1/2 = 2.0n0 is plotted in solid line and those by [74] are in dashed and dot-dashed lines (see Ref. [74] for notation). The grey band in the upper panel is the constraint from the low spinΛ = 300 +500 −190 obtained from GW170817 [5] and that in the lower panel is Λ ≤ 600 from GW190425. grained macroscopic description whereas the quarkyonic is a microscopic rendition of what's going on.
To understand the origin of the similarity and difference between [31] and PCM, we compare in Fig. 12 our prediction for P/ with the conformality band obtained by the SV interpolation method [31]. We see that our prediction is close to, and parallel with, the conformality band. There are basic differences between the two. First of all, in our theory, conformality is brokenthough perhaps only slightly at high density-in the system which can be seen from the deviation from the conformal band. Most importantly, the constituents of the matter after topology change in our theory is not (perturbatively) "deconfined" quarks. It is a quasiparticle of fractional baryon charge, neither purely baryonic nor purely quarkonic. In fact it can be anyonic lying on a (2+1) dimensional sheet [53,54]. That the predicted P/ deviates from the conformal band is indicating that the scale symmetry the EoS of our theory is probing is some Tidal deformabilities Λ1 and Λ2 associated with the high-mass M1 and low mass M2 components of the binary neutron star system GW170817 with chirp mass 1.188M (upper panel) and GW190425 with chirp mass 1.44M (lower panel). The constraint from GW170817 at the 90% probability contour is also indicated [79]. distance away from the IR fixed point with non-vanishing dilaton mass.
Gravitational wave.-We finally apply our theory to the description of the waveforms of the gravitational waves [81]. The purpose is to explore whether one can probe the possible continuous crossover from hadrons to quarks represented in terms of the topology change. For this purpose, we consider the typical values n 1/2 = 2n 0 and 3n 0 and the neutron star mass 1.5M .
The dominant mode of GW strain h + 22 multiplied by the distance of the observer to the origin R from BNS mergers is plotted in Fig. 13. The plot shows the location of the topology change affecting the number of the inspiral orbits, i.e., the number of the peaks in the inspi- ral phase which is the number of the peaks before merger, defined as the maximum of the amplitude of the GWs. Explicitly, the larger the n 1/2 , the more the number of peaks. This could be within the detection ability of the on-going and up-coming facilities, especially the groundbased facilities [82]. The effect of the topology change on the waveforms can be understood from the distribution of the matter evolution of the BNS merger shown in Fig. 14. It is found that the matter evolves faster when n 1/2 = 2n 0 (the EoS is softer) than when n 1/2 = 3n 0 (the EoS is stiffer). Therefore the stars merge more easily with a shorter inspiral period. This indicates that the waveforms of the gravitational waves emitted from the merger process as well as the matter evoluation could be sensitive to the EoS of compact stars (see, e.g., [83]). This observation explains the waveforms of Fig. 13. However, there is a caveat: given that no qualitatively striking differences are predicted for all other astrophysical observables so far studied for n 1/2 /n 0 = 2 and 3, it appears unnatural that the waveforms appear so different for only slightly different locations of the topology change. Furthermore since the transition involves no obvious phase change, at least within the framework, the seemingly different impact of the topology change density-which is a coarsegrained description of the phenomenon-seems puzzling. It would be interesting to see whether the "microscopic" models that simulate the quark degrees of freedom for hadron-quark continuity show similar sensitivity on the transition point. If indeed the waveforms were indeed very sensitive to the precise location of the cross-over, it would be extremely interesting.

VI. SUMMARY AND PERSPECTIVE
In this work we reviewed the effect of the topology change representing the putative hadron-quark continuity on dense nuclear matter. The hadron and nuclear matter properties obtained from the skyrmion crystal approach, supplemented with the presumed emergent scale and flavor symmetries, inspired the construction of the pseudo-conformal model of dense nuclear matter relevant to compact stars. Locked to the density dependence of hadron properties effected by the topology change at n 1/2 , the trace of the energy momentum tensor of the model turns out to be a nonzero density-independent quantity and induce the precocious appearance of the pseudo-conformal limit with v 2 s = 1/3, in stark contrast to what's widely accepted in the field [29].
So far, the pseudo-conformal model can describe the nuclear matter properties from low density to high density in a unified way. The nuclear matter properties calculated around the saturation density, the star properties such as the maximum mass, the mass-radius relation, the tidal deformability and so on all satisfy more or less satisfactorily the constraints from terrestrial experiments, astrophysical observations and gravitational wave detection.
Finally we state the possible caveats in and extensions of the model.
One can explicitly see from the above that although the tidal deformability predicted in the approach satisfies the currently cited constraint from the gravitational wave detection, it lies at the upper bound. But should the bound turn out to go to a substantially lower value than what's given presently, the description of the cusp structure of the symmetry energy would need a serious revamping. In the present framework, the tidal deformability probes the density regime slightly below the topology change density n 1/2 where the EoS is softer, and this is the density regime which is the hardest to control quantitatively in terms of the coarse-grained approach. It would require a more refined V lowk -renormalizationgroup treatment than what has been done so far in [27], including the approximation made for the anomaly effect in the GD (genuine dilaton) scheme and the role of strangeness mentioned below.
One possible way to resolve the above caveat is to include the corrections to the LOSS, applied so far, in such a way that, in addition to the mass parameters, the coupling constants also carry IDDs. This procedure may change the property of the EOS in the vicinity of n 1/2 . As a consequence the sound velocity after the topology change may also expose bumps, i.e., fluctuations from the conformal limit, because of the explicit breaking of the conformal symmetry. However if the corrections from the explicit breaking of the conformal limit are taken as chiral-scale perturbation, the global picture of the compact star discussed would remain more or less intact.
Another point is the density at which the hidden scale and local flavor symmetries emerge. This is encoded in the IDDs of the hadron parameters such as pion decay constant, dilaton decay constant, ρ-N-N coupling and meson masses. By checking the effect of the location of the emergent symmetries on the star properties, one can also extract the information on the emergent symmetries and the phase structure of QCD at low temperature.
Lastly we have left out the strangeness in the present discussion. It seems to have worked well without it in our approach up to now. But there is of course no strong reason to ignore it. It could very well be that strangeness does play a crucial role but indirectly, buried in the coarse-graining in the approach. Or it could also be that strangeness does not play a significant role up to the density involved in compact stars. The chiral-scale effective theory that the pseudo-conformal description relies on is based on three flavor QCD with the scalar f 0 (500) taken on the same footing as the pseudoscalar mesons pion and kaon. There is however a good reason to believe that in nuclear dynamics, the dilaton scalar is strongly affected by medium whereas the kaon is not. Implementing the strangeness in our approach would require doing the V lowk RG for 3-flavor systems with the hyperons treated on the Fermi sea together with the nucleons as Fermi-liquid theory. This would then involve the kaon condensation as well as the hyperons as bound states of skyrmions and kaons. As argued in [84], it could postpone the role of strangeness to a much higher density than relevant to the most massive compact stars stable against gravitational collapse. What happens beyond, such as color-flavor-locking, would be irrelevant to the problem.
One excuse for ignoring the strangeness could be that the whole thing works without it, so why not adopt the spirit "Damn the torpedoes! Full speed ahead!" 3 and proceed until hit by a torpedo?