Universal relations for compact stars with exotic degrees of freedom

The nature of the highly dense matter inside the supernova remnant compact star is not constrained by terrestrial experiments and hence modeled phenomenologically to accommodate the astrophysical observations from compact stars. The observable properties of the compact stars are highly sensitive to the microscopic model of highly dense matter. However, some universal relations exist between some macroscopic properties of compact stars independent of the matter model. We study the universal relation including the stars containing exotic degrees of freedom such as heavier strange and non-strange baryons, strange quark matter in normal and superconducting phases, etc. We examine the universal relations for quantities moment of inertia - tidal love number - quadrupole moment. We also study the correlation of non-radial f-mode and p-mode frequencies with stellar properties. We find the f-mode frequency observes the universal relation with dimensionless tidal deformability but the p-mode frequency does not show a good correlation with stellar properties. The p-mode frequency is sensitive to the composition of the matter. We find that universal relation is also applicable for stars with exotic matter in the core of the star with several models of exotic matter.


Introduction
The massive stars end their lives by core-collapse supernova explosion leaving highly compact stars (CSs) as central objects.The CSs have an average density ∼ 10 14 gm/cm 3 which is several times the normal nuclear saturation density (n 0 ).These compact objects accommodate the extremely dense matter in the universe and thus allow insights into various new aspects of dense matter physics.Such high densities are not attainable in any laboratory or terrestrial experiments and hence the exact composition of matter as well as the inter-particle interaction inside the compact objects is unknown.The equation of state (EOS), which describes the relationship between energy density (ε) and pressure (p), at extremely high densities, is among the most uncertain aspects of nuclear physics.Consequently, many phenomenological models are discussed with different compositions and different inter-particle interactions.The most studied and discussed possibility of matter at such high density is pure nucleonic matter -the matter composed of mostly neutrons with some admixture of protons and electrons.However, the appearance of exotic degrees of freedom in the interior of massive NSs where the density is a few times n 0 tends to be feasible, although it remains an open question.For example, there are possibilities of the appearance of strange baryons [1][2][3][4][5][6], non-strange heavier baryons [7][8][9][10][11][12][13], Boson condensates [14][15][16][17][18][19] etc.Another possibility is that at much higher density, the baryons, may get decomposed into their constituent quarks and generate a region of deconfined strange quark matter (SQM) [20,21].The star may be entirely composed of SQM up to the surface due to some stability conditions [22][23][24].Another possibility is a hybrid star (HS) composed of SQM at the core surrounded by baryonic matter up to the surface [25][26][27][28][29].
The CSs have an interior structure that is highly sensitive to their EOS.Consequently, these factors dictate the external characteristics of the objects in ques-tion, including their mass and radius; their deformability, as indicated by their quadrupole moment Q and tidal Love number Λ; and their rotation rate, which is defined by their moment of inertia I.The macroscopic properties of CSs are highly sensitive to the microscopic properties of highly dense matter which are still unknown and model-dependent as of now.Hence, the only way to minimize the theoretical uncertainties regarding highly dense matter is to fit the proposed models with the astrophysical observations coming from these CSs.
However, Yagi and Yunes in 2013 discovered that some combinations of physical parameters do not depend on the EOS and follow universal relations [30,31].They showed that the relations between the moment of inertia (I), tidal Love number (Λ), and quadrupole moment (Q) of CSs are independent of EOS of highly dense matter irrespective of whether the matter is pure nucleonic or deconfined strange quark matter.The combinations of the above-mentioned variables obey a universal relation.Recently, universal relations have been getting some attention .Moreover, later on some empirical relations of non-radial oscillation frequencies with stellar parameters and universal relation have been reported [55][56][57][58].
Universal relations are crucial because they enable us to calculate the others if we know one parameter.When it comes to CSs, it is possible to compute the love number Λ and Q without physically measuring them if I can be measured in some way.This is extremely useful, as both are challenging to quantify for binaries separated by a significant distance.By accurately measuring the love number in coalescing binaries, it is possible to infer the Q and, consequently, gain insight into the spins of the coalescing binary CSs.Most of the studies regarding universal relations have been carried out with neutron stars composed of pure nucleonic matter.The I − Q relation for rapidly rotating NSs was first studied by Doneva [59].They found that the I −Q relation is broken and becomes more EOSdependent for NSs with a fixed frequency.However, it was soon found by Pappas & Apostolatos [60], and Chakrabarti et al. [43] that the I − Q relation can remain approximately EOS-insensitive if one chooses suitable dimensionless parameters instead of dimensional quantities.Bandyopadhyay et al. [54] studied the universal relations with Λ-hyperons and anti-kaon condensates.Eemeli Annala et al. [61] reported a maximum deviation of 15% to 20% in universal relations due to strong first-order phase transition at low density.Raduta et al. [47] discussed the universal relation considering hot stars with heavier strange and non-strange baryons and quarks in the core.The universal relation for hot and cold rapidly rotating HSs has been considered by Largani et al. [51].The implication of GW170817 with universal relations for HS has been studied by Paschalidis et al. [33].At high enough density (µ >> T ), there may exist the Color-Flavor-Locked (CFL phase), which is a superconducting phase of SQM [62][63][64].In recent works, we extend the study of universalities in the case of CS with strange and heavier nonstrange baryons as well as SQM in the Color-Flavor-Locked (CFL) phase inside the interior of the star along with the pure NS.
There are several quasi-normal modes like the fundamental f-mode, pressure p-modes, gravity g-modes, spacetime w-modes, etc. [65][66][67], each classified based on the restoring forces that act to bring it back to equilibrium.For example, the f-and p-modes, which are acoustic waves in the star, are restored by fluid pressure while g-modes which arise due to discontinuities in density or temperature are restored by gravity (buoyancy).Various previous works have explored non-radial modes in cold and finite temperature neutron stars [68][69][70][71].Several years ago, Andersson and Kokkotas [55] put forward an empirical relation for the f-mode oscillation frequency ω with the stellar parameters mass and radius, based on the Newtonian theory of stellar perturbations.They observed that in full GR, ω depends almost linearly on the square root of the average density.Another relation, based on estimates using the quadrupole formula, was established for the damping time due to gravitational wave emission, τ.Later, Benhar et al. [56] presented further results that included more and newer equations of state, updating the fits from [55].The average frequencies ω were systematically lower than the one for the old EOS sample, which they attributed to the fact that the new sample included stiffer EOSs.Later [57] thoroughly studied the universality.They considered a wide range of masses and the EOSs to minimize the uncertainty.In this work, we took it further and computed the universality for non-radial oscillation frequencies for all the possible families of CSs.In this work, we adopt the Cowling classification, where the various modes are separated by the number of radial nodes [72,73].The f-mode, whose frequency lies between those of the pand g-modes, has no radial node number.Quadrupolar oscillations (l=2) of all modes lead to the emission of GWs.With the advent of enhanced next-generation telescopes like the Cosmic Explorer and the Einstein telescope which carry about 10 times the sensitivity of Advanced LIGO, the possibility of detection of these modes increase [74].Although the high-frequency p-modes aren't expected to be detected by the next-generation GW detectors, we include their study for the sake of completeness.
In the Sect.2, we describe the models of highly dense matter considered in the work and structure of stars composed of matter within these models.In Sect.3, we discuss the non-radial modes and determine their frequencies for CS with heavier baryons and SQM in the CFL phase inside the core.Next, in Sect. 4 we study the universal relations between different quantities.

Matter model and consequent stellar structure
In this work, we study the universal relation for different observable quantities related to CS composed of matter with different possible compositions.As mentioned in the Sect. 1, the most discussed stellar model is with nucleonic matter.With this, we consider the appearance of heavier strange and non-strange baryons and SQM in the core of the star.In this context, first, we briefly discuss the star with baryonic matter only in the subsequent subsection.Then, in the next subsection, we discuss the HS with SQM in normal and CFL phases inside the core.

Baryonic star
Inside CS highly dense matter is thought to be composed of mostly neutrons with a small admixture of protons and electrons.With the advent of astrophysical observation, now it is established that the CS can be as massive as 2 M ⊙ or even more than that which leads to the assumption of the existence of matter at several times nuclear saturation density near the core of the CS.In that case, near the core, there is a possibility of the appearance of heavier baryons.We consider the Relativistic Mean Field (RMF) approach for baryonic matter at high density.Within the RMF model, we assume the interaction between nucleons is mediated by the exchange of isoscalar-scalar σ , isoscalar-vector ω and φ , and isovector-vector ρ mesons.The recent observations of massive CS candidates indicate that the density inside the core is much higher than n 0 with which the appearance of heavier strange and non-strange baryons is highly probable.On the other hand, the appearance of heavier baryons softens the matter reducing the maximum attainable mass of CS which contradicts the recent astrophysical observations.To come out of this problem, one should consider the baryon interaction density-dependent.Hence, in our current work, we stick to that formalism within the RMF approach.We consider the matter composed of only nucleons as well as the possibility of the appearance of strange baryons belonging to baryon octet and ∆-resonances [13].
In the RMF approach, the strength of baryon coupling to mesons is taken as a parameter.The values of the parameters are fixed from observed nuclear matter properties at n 0 , such as values of n 0 , binding energy per nucleon E/A at n 0 , symmetry energy E sym at n 0 and its derivative with respect to n 0 , incompressibility K, etc. Different sets of parameter values can reproduce these observed nuclear mat- Table 1: The range of nuclear parameters at n 0 .The lower row represents the values (along with references) of the saturation properties obtained from various terrestrial experiments and methods.
ter properties.Depending on this several models of parametrizations are there for baryonic matter models.In this work, we employ the parameterizations namely DD2 [75], DD-ME2 [76] and DD-MEX [77] for baryonic matter which covers the nuclear saturation properties as shown in Table 1.The readers may refer to Ref. [80] for further details of the comparison of experimental data from finite nuclei and heavy-ion collisions with different microscopic model calculations.
The equation of state (EOS) for only nucleon matter is shown by the curves with circles in Fig. 1.When we consider the possibility of the appearance of strange baryons (hyperons) then the EOS softens as shown in the curves with the thin dotted line.The appearance of non-strange heavier baryons ∆ further softens the EOS at low density and affects the maximum attainable mass as evident from curves with the thick dotted line in the same figure.

Hybrid star
With the increase in density towards the center, the deconfinement of quarks is another possibility that leads to the formation of SQM at the core of the star, making an HS.Naturally, inside HS, the deconfined SQM at the core is surrounded by pure hadronic matter.In our recent work we consider the outer baryonic part of the HS to be purely nucleonic.Inside the core for quark matter, we consider two possibilities: one is normal SQM matter and another possibility is superconducting quark matter in the CFL phase.Various models are available for the description of SQM.However, the MIT bag model and NJL model are well established.The original MIT bag model considers free quarks confined within a bag.Subsequently, some interaction between the quarks has been incorporated [21,81].In our work we consider the MIT bag model with interaction and denote it as MIT bag model.In another model, repulsive vector interaction terms are included within free quark MIT bag model which is known as vector bag model (vBAG) [24,[82][83][84].For the superconducting CFL phase, we use the MIT bag model including the gap term as mentioned in Ref. [85][86][87].For the universal relations we consider the inner core SQM inside the HS considering these models mentioned above, which are briefly discussed in the subsequent sections.

Normal quark matter
The normal SQM inside the core of the HS is composed of up(u), down(d), strange(s) quarks, and electrons(e).
MIT bag model: Here we briefly describe the MIT bag model with interaction as introduced by Farhi and Jaffe [81].In this model the thermodynamic potential of each quark is given by with µ q and m q the chemical potential and mass of each quark q respectively, α s the coupling parameter for strong interaction and σ the renormalization scale.The thermodynamic potential of the electron is with µ e the electron chemical potential.Then the EOS can be obtained through the relations Here, B is the bag pressure that separates the quarks confined within the bag from the vacuum.where and the number density is We consider two sets of parameters with this model one is α s = 0.8 fm 2 with B 1/4 = 160 MeV and another is for softer EOS α s = 0. vBAG model: The Lagrangian density of SQM with the vBAG model is given by the following equation as where ψ q is the field of quark q considering it as a fermion and m q is its mass.Similarly, V µ and m V are the field of the mediator and its mass respectively.Θ is the heavy side function that is unity inside and zero outside the bag.After solving equations of motions for the Lagrangian density, the chemical potential gets shifted as The energy density of SQM matter within the vBAG model is here V 0 is the mean value of V µ in ground state.From energy density, pressure can be obtained by the following equation where n q signifies number density of quark q.For simplicity, we define the parameter G V = (g qv /m V ) 2 .In this work, we take two sets of parameters, one is a higher value of this parameter G V = 0.25 fm 2 with B 1/4 = 180 MeV and another is a lower value G V = 0.17 fm 2 with B 1/4 = 172 MeV for softer EOS.With the vBAG model, we consider masses quarks the same as the MIT bag model.The inclusion of vector interaction stiffens the EOS compared to the MIT bag model.For HS with SQM core considering this model, the stiffer EOS is shown by the dash double dotted line in Fig. 1 with DDME2 nuclear EOS.
NJL model: The Lagrangian density describing the dynamics of SQM within the NJL model with vector interactions and the 't Hooft six-quark interaction term is given by [88][89][90][91]: where λ a are Gell-Mann matrices.ψ = (u, d, s) T is fermionic field of quarks and m0 is diagonal matrix with elements (m u0 , m d0 , m s0 ).The scalar interaction and 't Hooft interaction strengths are denoted by the parameters G S and K respectively.In Fig. 1, the HS EOS with NJL model for quark matter is represented by the dashed line.
Here, the constituent quark mass M f and the constituent quark mass for vacuum (n u = n d = n s = 0) M f 0 are given by and respectively with the chiral condensates We consider the mass of quarks as m u0 = m d0 = 5.5 and m s0 = 140.7 MeV.Value of ultraviolet cutoff parameter in momentum space (Λ) is considered 602.3 MeV.Parameters G S and K are related to Λ as G S = 1.835/Λ 2 and K = 12.36/Λ 5 respectively.The vector interaction parameter is G V = 0.2G S .

Superconducting quark matter
In this phase, quarks are paired up in such a way that they form color-neutral and flavor-neutral condensates, and this leads to a unique and stable ground state of matter.In such a scenario, these paired quarks would move together with correlated momentum, somewhat analogous to Cooper pairs in superconductors.
Studying the CFL phase in neutron stars is intimately linked to quantum chromodynamics (QCD), the theory of strong interactions.Here, each quark has a common Fermi momentum (ν) and number density (n) to form Cooper pairs Here, ∆ is the CFL energy gap and quark number chemical potential is µ = (µ u + µ d + µ s )/3.The thermodynamic potential of this phase is given by with The pressure (p), and energy density (ε) of CFL phase are The quark matter EOS with the CFL phase is represented by the solid line.With these models of SQM we construct the core of the HS surrounded by nucleonic matter with the models discussed in the Subsect.2.1.To construct the HS we consider Gibb's construction in which equilibrium between the baryonic phase and quark phase is established in a mixed phase during the phase transition.In this phase, global charge neutrality is maintained rather than local charge neutrality as where χ is the volume fraction of quark matter to baryonic matter and ρ c is the charge density.Superscripts B and Q denote the baryonic phase and quark phase respectively.The energy density of mixed-phase is given by the following equation here ε Q and ε H are the energy densities of quark matter and baryonic matter respectively.The baryonic density of matter in the mixed phase is Baryonic matter density for quark matter is ρ Q = (n u + n d + n s )/3.The pressure of mixed-phase is equal to the pressure of quark matter and baryonic matter.With this, we construct the HS structure with different combinations of the baryonic matter model and SQM models as discussed above.

Comparison of the models
It is well known that nuclear matter is the stiffest among all possible compositions of matter.The appearance of exotic matter at higher density softens the matter undoubtedly.From the comparison of the appearance of different exotic degrees of freedom, we notice that the appearance of the quark softens more compared to the Here we show the nuclear matter with three parametrizations as discussed in sec.2.1.We see that DDMEX is the stiffest and DD2 is the softest among them for nuclear matter.To understand the effect of exotic matter, we consider the nucleonic part of the matter with DDME2 parametrization.When considering the appearance of hyperonic matter it makes the EOS softer.In the case of hyperonic matter with ∆ baryons, ∆ starts appearing very early and makes matter soft at a low-density regime.Without ∆ other strange baryons appear a little later.It also softens the matter more at high density.SQM appears at comparatively higher density depending on the model and model parametrizations considered.The stiffness of the matter decreases progressively from the vBAG model to the MITBAG model, then further to the NJL model, culminating in the CFL phase.In very high-density regions matter with SQM is softer than the matter with hyperons.The parametrization for this plot is detailed in Table 2.These models incorporate a broad range of soft to stiff EOSs, making them ideal for testing the robustness of our universal relations.In Fig. 1, every EOS has central energy density corresponding to maximum mass less than 1250 MeV/fm 3

Non-radial modes
Now we describe the methods of finding the frequencies of non-radial modes relevant for the CSs.The fully general relativistic calculation is somewhat cumbersome, the most used approximation is the Cowling approximation.We have used that as well as the full general relativistic calculation to find the non-radial frequencies.

Relativistic Cowling approximation
First, we describe the Cowling approximation.In this approximation, the metric is assumed to be unperturbed.Since the non-radial oscillation formalism is developed in linearised gravity, the line element of a non-rotating neutron star is taken as: The fluid Lagrangian displacement vector is assumed to be: [92] where Y lm (θ , φ ) are the spherical harmonics.W and V are perturbative fluid variables, which are functions of r with a harmonic time dependence W (t, r) = W (r)e iωt and V (t, r) = V (r)e iωt .The mode frequencies ω are thus found by solving the following system of ordinary differential equations [92]: where ε is the energy density and p the pressure.The dash ( ′ ) represents a derivative with respect to the radius.The boundary conditions at the center of the star are W (r → 0) = Ar l+1 and V (r → 0) = −Ar l /l, where A is an arbitrary constant, usually taken to be 1.
The surface boundary condition corresponds to the fluid pressure vanishing at the surface of the star (r = R): The lowest frequency solution to the oscillation mode equation has no radial node and thus corresponds to the f-mode frequency.The next highest solution has one radial node and is thus the first p-mode (p1-mode).Although [73] suggests that this classification fails before 0.4 seconds of the bounce, for the models in consideration in this work, the Cowling classification holds good.The number of radial nodes was found by counting the number of times the perturbative variables W and V become 0 within r < R since this condition ensures that the three-velocity of the fluid (which is the time derivative of the fluid Lagrangian displacement vector), becomes 0 at those points.The correctness of our code was ensured when our results matched with those reported in [68,71].

Fully general relativistic calculations
In fully general relativistic calculations, we consider the perturbation in metric too.The equations that need to be solved and the various techniques required to solve them have been examined in great detail in several previous works [66,67,69,[93][94][95] and we simply state them for completeness.
In this work, we follow the method of direct numerical integration, first proposed by [66,67] and refined by [96].The perturbed metric is taken as: while the fluid displacement vector in terms of the perturbation functions W (r) and V (r) is: Lindblom and Detweiller [66,67] introduced a new function X(r) to replace V (r) and the relations between all these functions are as follows: The system of differential and algebraic equations, 31 completely describes the perturbations inside the star.The differential equations to be solved can be stored in an array Y = {H 1 , K,W, X}.This system is clearly singular at r = 0 and numerically, it will blow up at values of r close to 0. Thus near the center, Y (r) is approximated as ) and the various terms of this approximation are given in [96].At the surface of the star, the pressure perturbations, and thus X must be 0. To solve equation 31, we follow the method outlined in [66].We start off with 3 linearly independent solutions at the surface, and 2 linearly independent solutions at the center and integrate them to some point inside the star where they get matched.A linear combination of these solutions, with the coefficients obtained after matching provides the true values of H 1 and K at the star's surface.These variables and H 0 are the only variables defined outside the star, where the perturbation equations reduce to the Zerilli equation [69,96]: where V (r * ) is the Zerilli potential, r * is the tortoise coordinate, r * = r + 2M ln (r/2M − 1), n = (l − 1)(l + 2)/2 and b = M/r, with M being the total mass of the star.
In case of a first-order phase transition inside the HS, we impose additional conditions that ensure the continuity of H 1 , K, W , and X across the radius of discontinuity [95].
The perturbed metric outside the star describes a combination of outgoing and incoming gravitational waves, which is the general solution to the Zerilli equation.We are interested in the case of purely outgoing waves, representing the quasinormal modes (QNM) of the star.At the surface of the star, where r = R, the fluid variables can be converted to the Zerilli ones using [69]: Here After continuing the integration of the Zerilli Eq. 32 to sufficiently far away from the star (∼ 50ω −1 ), the solution can be approximated as a linear combination of incoming and outgoing waves as Z(r * ) = A − (ω)Z − (r * ) + A + (ω)Z + (r * ) where Z − represents the outgoing wave, Z + the incoming wave and A − and A + their amplitudes.At a large enough radius, Here Z + is the complex conjugate of Z − (and hence A + the complex conjugate of A − ) and, [69] β 0 can be any complex number that represents an overall phase.By matching the solution of Z(r * ) and dZ(r * ) dr * obtained from 34 with the above equation, we can find the amplitude A + with a simple matrix inversion [69].The frequency of the QNM corresponds to that ω which gives A + = 0.
To find the QNM frequency and its damping time we first find A + , which in general will be a complex number, for several real values of ω close to the original guess.We then perform a complex polynomial fitting to approximate a parabola passing through the A + points corresponding to the ω values.The root of this parabola which has a positive imaginary part is our QNM.We then take the real   part of this ω and repeat the entire procedure several more times till the desired tolerance is reached.The real part of this final ω is the frequency of the QNM.The inverse of the imaginary part is the corresponding damping time.
As with the Cowling case, the validity of our code was checked thoroughly by comparing our results with those in [68,95].Moreover, we also calculated the oscillation modes and damping times using Thorne [97] Ferrari's [94] Breit-Wigner resonance fitting approach and although the damping times could not be obtained with good accuracy for all cases, the frequency values matched exactly to those obtained by the Lindblom [66,67] approach.

Universal relations
With the CS models discussed above, we study the universal relations for the CSs.First, we study the I-Love-Q relations.For that, as usual, we compute the dimensionless moment of inertia Ī = I/M 3 , dimensionless quadrupole moment Q = −Q/M 3 χ 2 L , and dimensionless tidal love number or tidal deformability Λ = λ = λ /M 5 .Here, I is the moment of inertia, Q the spin-induced quadrupole moment [98,99] and χ L = J/M 2 (J being total angular moment) the dimensionless angular momentum.We compute these quantities for CSs with pure nucleonic matter, hyperonic matter and for HS with normal and superconducting SQM at the core with the help of RNS code [100,101], which is built on solutions of Einstein's field equations for axially symmetric and stationary space-time in spherical coordinates.In this work, we choose the minimum possible frequency of rotation frequency for all the EOSs to find the solutions Ω ∼ 480 Hz which is much below the mass shedding frequency.
We show the universality in Ī − Q, Q − Λ and Ī − Λ plot along with the error due to polynomial fitting various possible composition and parametrization of matter EOS as discussed above in different panels of Fig. 2. The dashed line represents the fitted line by Yagi and Yunes (YY) [30].The effectiveness of curve fitting can be evaluated using the coefficient of determination For I − Q relation, shown in the upper left panel of Fig. 2, we get the value of R 2 = 0.9998 very close to 1, indicating a good fitting.The maximum deviation, being less than 5%, indicates the universality of the I − Q relationship.All the Ī values depend on Q rather than properties of matter as expected in a universal relation.The YY fitted line shows some deviation from our fitted line at low density due to a difference in rotation speed [102].The coefficients for the polynomial fit are given in Table 3.The variation of Q with Λ is shown in the upper right panel of Fig. 2. We achieve a good fit with R 2 = 0.9997, indicating a high degree correlation.Despite a maximum deviation of approximately 4.65% with the matter composition of nucleonic matter with hyperon in DDME2 parametrization, though it remains below the 5%, these findings suggest a good universal relationship between these parameters.Here, the YY fitted line is very close to our fitted line.Similarly, we plot the relation between Ī and Λ in the lower right panel of Fig. 2. We find R 2 = 0.9993 and the maximum deviation is 1.56% at low density,  which means these parameters are strongly correlated.Our fitted line is almost parallel and close to the YY fitted line.
The universal I − Λ − Q relations have been observed in previous literature [4,102] for CSs with nuclear matter, nuclear hyperonic matter, nuclear hyperon-∆ matter and SQM with NJL model.We observe here that the same universality is applicable for HSs with SQM in the MIT bag model with ad-hoc and vector interaction as well as the CFL phase at the core surrounded by the nucleonic matter.
Next, we focus on correlations in QNMs of non-radial oscillation.We show the relationship of the dimensionless angular frequency ω = 2πGMν/c 3 of nonradial oscillations in f-and p-modes with Λ.Here, M is the star's mass, and ν is the frequency of the respective non-radial modes in Hz.
We plot the variation of f-mode frequencies ω f c using Cowling approximation  and ω f with Λ using full GR in the upper left and right panels respectively along with the analytic fit and fractional error in the angular frequencies ω f c and ω f in Fig. 3.The value of R 2 with Cowling approximation is 0.9993 and the maximum deviation is 2.65% with DD2 MITBAGs near the maximum mass.If we calculate f-mode frequency in a fully general relativistic framework, this becomes a robust universal relation.Now, R 2 increases to the value 0.9999, which is very close to 1, and the maximum deviation reduces to 0.07%.The strength of this correlation has been previously investigated in the references with only pure nucleonic matter and SQM [69,103].
Next, we plot the variation of p-mode frequencies ωpc using Cowling approximation and ωp with Λ using full GR in the lower left and right panels respectively along with the analytic fit and fractional error in the angular frequencies ωpc and ωp in the lower panels of the same figure.Here, the maximum deviation is about 13.2% from the fitted line with R 2 = 0.9730.The low value of R 2 and the high value of maximum deviation suggest that it's a poor fit.Further, we plot p-mode frequency with full GR calculations in the lower panel.Now, the value of R 2 further drops to 0.9597, and the maximum deviation increases to approximately 15%.This means the inclusion of perturbation due to gravitational potential or full GR calculation makes it a poorer fit.However, the variation of p-mode frequencies ωpc using Cowling approximation and ω using full GR with the stellar mass shows better correlations.This is evident from the left and right panels of Fig. 4. In this case, we obtain the value of R 2 = 0.9976 with Cowling approximation and R 2 = 0.9917 with full GR calculations.These values of R  a better fit as compared to variation with Λ.The maximum deviation also drops to 6.68% with Cowling approximation and 7.07% with full GR.It seems that the p-mode frequency of oscillation is sensitive to matter composition rather than the structure of the star.
If we examine the correlation between the frequencies of non-radial oscillation from Cowling approximation and full general relativistic calculations, we observe, that the effect of metric perturbation increases with stellar mass as expected.Moreover, for f-mode frequency, the effect of metric perturbation is almost independent of matter composition while for p-mode frequency, the metric perturbation effect does depend on the matter composition.This is evident from the left and right panels of Fig. 5.This is because the p-modes are more sensitive to the matter distribution inside the star [55], possibly because of a radial node.The position of the node would be affected by the distribution of matter, which is in turn affected by the back-reaction of the perturbing spacetime on the star.It is worth mentioning that the p-mode frequency is sensitive to pressure at low density as the p-mode frequency correlates with pressure and other matter parameters [68].We plot the relationship between frequency by Cowling approximation and by applying full GR calculations for both f-and p-modes in Fig. 5.This explains the poorer correlation for the p-mode oscillation frequency in the case of full general relativistic calculation.For relations we discussed above, we use the following polynomial for curve fitting y = a 5 x 5 + a 4 x 4 + +a 3 x 3 + a 2 x 2 + a 1 x + a 0 (39) The coefficients for the curve fitting are given in Table 3.The components of 4th-order are found to be consistent with the ones reported in ref. [104].For a bit higher accuracy than the 4th-order, we approach the 5th-order polynomial component fit in cases of non-radial oscillations.The contribution from the 5thorder component can be seen in Table 3.The higher degree of polynomials can also be taken for a little bit more accuracy like Zhao et al. [69] but the coefficient of the 6th-order becomes very small.

Conclusions and summary
The CS born after the supernova explosion contains matter at a density a few times n 0 .At that high density, the matter properties are not well constrained from experimental as well as theoretical points of view.Consequently, the matter at high density is modeled for different compositions and interactions between constituent particles.The macroscopic properties of CS are highly sensitive to highly dense matter models.However, a few relations between some macroscopic quantities are independent of microscopic models of highly dense matter.Yagi and Yunes first pointed this out in 2013 [30,31] and further supported by many studies [32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][49][50][51][52].Most of the studies of universal relations have been done with stars composed of pure nucleonic matter -the neutron star.Other possible configurations of compact stars are baryonic stars with heavier baryons in the inner part of the star and the HS.The universal relation with HS configuration has also been studied in refs.[33,51] with Maxwell construction of phase transition to SQM in normal phase.In this current work, we examine the universal relation for HS with Gibbs construction to SQM in the normal phase with MIT bag model with ad-hoc or vector interactions and in the CFL phase and we find that this kind of HS also maintains the same universality with the CS made up of pure baryonic matter.
Then we find the universal relations between the frequencies of non-radial f-mode oscillation with dimensionless tidal deformability for the HS configuration which is compatible with the pure baryonic star configuration.However, the frequency of non-radial p-mode oscillation does not correlate well with dimensionless tidal deformability.The correlation improves with mass but is not up to the mark.We find that p-mode oscillation frequency is sensitive to the composition of matter in both Cowling approximation and full general relativistic calculations.
3 fm 2 with B 1/4 = 174 MeV.With this quark model, we consider the up quark mass m u = 4 MeV, down quark mass m d = 7 MeV, and strange quark mass m s = 100 MeV.The value of renormalization is 300 MeV here.The appearance of SQM naturally softens the EOS.With a small strength of interaction parameters and a high value of bag parameter, the EOS is comparatively softer.The stiffer EOS of HS is shown by the dash-dotted line in Fig. 1 with this model.

Figure 1 :
Figure 1: The variation of pressure with energy density

Figure 2 :
Figure 2: I − Λ − Q relations along with the analytic fit and the fractional error from the fitting function.Upper left panel: for I − Q relation, upper right panel: for Q − Λ relation, and lower right panel: for I − Λ relation.The solid line represents the fitted curve for all data points.The sub-figures display error bars, with the dashed line indicating a 5% error.The nomenclature for the shapes is depicted in the lower left panel.

Figure 3 :
Figure 3: The relation of dimensionless frequency with tidal deformability.Upper left panel: for f-mode frequency calculated by Cowling approximation and upper right panel: with full GR framework.Similarly, lower panels are for p-mode.The solid line represents the fitted curve for all data points.The sub-figures display error bars, with the dashed line indicating a 5% error.The nomenclature for the shapes remains the same as the Fig. 2.

Figure 4 :
Figure 4: The relation of dimensionless p-mode frequency with star's mass.Left panel: for frequency calculated by Cowling approximation and Right panel: with full GR framework.The solid line represents the fitted curve for all data points.The sub-figures display error bars, with the dashed line indicating a 5% error.The nomenclature for the shapes remains the same as the Fig. 2.

Figure 5 :
Figure 5: Variation of frequency calculated in full GR framework with frequency calculated by Cowling approximation.Left panel: for f-mode and right panel: for p-mode.The solid line depicts the fitted curve for all data points.The nomenclature for the shapes remains the same as the Fig. 2.

Table 2 :
and pure SQM does not appear inside the core of HS.So, for HS we include softer EOSs with vBAG and MITBAG models.The different combinations of parameters are given in Table2.In this table, we denote the softer parameter sets within a model with the notation (s) beside the name of the models.Set of parameters we used for HSs EOS, B.M. denotes baryonic matter and Q.M. denotes quark matter