Effects of Hadron-Quark Phase Transitions in Hybrid Stars within the NJL Model

We study local and non-local Polyakov Nambu-Jona-Lasinio models and analyze their respective phase transition diagram. We construct hybrid stars using the zero temperature limit of the local and non-local versions of Nambu-Jona-Lasinio model for quark matter and the modern GM1(L) parametrization of the non-linear relativistic mean field model for hadronic matter. We compare our models with data from PSR J1614-2230 and PSR J0343+0432 and also from GW170817 and its electromagnetic counterpart GRB170817A and AT2017gfo. We study observational signatures of the appearance of a mixed phase as a result of modeling a phase transition that mimics the Gibbs formalism and compare the results with the sharp first-order phase transition obtained using the Maxwell construction. We also study in detail the $g$-mode associated with discontinuities in the equation of state, and calculate non-radial oscillation modes using relativistic Cowling approximation.


Introduction
One of the purposes of Quantum Chromodynamics (QCD) is to describe how quarks and gluons interact to form the bound states composing hadronic matter, which can basically be classified into baryons and mesons. Experience indicates that under ordinary conditions of temperature and density, the phenomenon known as confinement occurs, in which quarks and gluons can only form such bound states. However, at sufficiently high densities and/or temperatures matter of unbound quarks and gluons, known by the generic name of quark-gluon plasma (QGP), can be formed. This has motivated the theoretical study of the possible phases of QCD and has revealed a complex phase structure. Thus, the study of the QCD phase diagram represents a theoretical and experimental challenge to date. The extreme conditions in which QGP is formed are difficult to produce under laboratory conditions, thus leaving few physical situations for experimental study, for example through heavy ion-collisions [1]. One of the natural scenarios where deconfined phases of QCD could appear is the interior of the so-called compact objects like neutron stars (NSs), which correspond to the region of high densities and low temperatures in the QCD phase diagram.
NSs are the densest known objects in the Universe with a mass, M∼1.4 M and a radius, R∼10 km. Astronomical observations of these objects can be used to test several high-energy physical theories [2]. In the interior of NSs, matter is compressed several times the nuclear saturation density. Due such extreme conditions, it is expected the appearance of new degrees of freedom as hyperons and deconfined quark matter. However, the equation of state (EoS) describing the dense matter composing the cores of NSs is still unknown.
The observation of the 2 M binary pulsars (PSR J1614-2230 [3,4] and PSR J0343+0432 [5,6]) has imposed a lower bound to the maximum mass of NSs. Therefore, a stiff enough EoS describing the dense matter in such objects should be able to fulfill this mass constraint. These mass measurements challenge high-density EoS and turn attention to the strange quark matter hypothesis, which suggests that the true ground state of matter is quark matter [7,8].
On the other hand, the event of gravitational waves (GW) GW170817 emitted during a binary NS merger and its subsequent electromagnetic counterpart (GRB170817A and AT2017gfo) impose an upper limit to the effective tidal deformability,Λ, of the binary system. This parameter depends on the microphysical properties of the matter composition of the star and it is a measure of the influence of a star's internal structure on the waveform. Effects on this quantity related to the appearance of quark matter in the inner cores of compact stars have been studied in Refs. [9,10]. Analysis of the data also serve to constraint the chirp mass of system, which determines (to leading-order) the evolution of binary NSs orbit and also the evolution of the GWs frequency signal during the inspiral phase. According theΛ < 800 constraint, several works conclude, to a 90% level, that the radius of a 1.4 M NS, R 1.4 , can not exceed ∼13.6 km (see, for example, Refs. [11][12][13][14][15][16][17]). This is a clear example of the powerful tool that represents the detection of GWs for NS astrophysics, setting additional constraints on the EOS of dense matter composing such stars. It is expected that during LIGO's third observing run (planned to begin in April 2019) more NS-NS mergers would be detected from which new updates of the high-density EoS constraints can be established.
The possible existence of quark matter inside of NSs remains an open issue because QCD cannot be solved for dense matter in astrophysical environments with the tools usually employed by relativistic quantum field theories. Due the severe technical difficulties related to the non-perturbative nature of QCD at low and intermediate densities, several effective models have been proposed (MIT Bag model [18], Nambu-Jona-Lasinio (NJL) model [19,20] non-linear sigma model [21], Field Correlator Method [22,23]). These models incorporate some (but not all) of the most relevant features of QCD.
The possibility of a hadron-quark matter phase transition in the inner core of NSs has been studied in several works (see, for example, Refs. [24][25][26]) and some of the astrophysical implications of the existence of such dense matter have also been studied. In this paper we will work with both, the standard NJL and a non-local extension of the NJL model for the description of quark matter. We analyze the possibility of a hadron-quark phase transition in the interior of cold NSs using a modified GM1 relativistic mean field model, the modern GM1(L) parametrization, to describe the hadronic phase. In Refs. [27,28] the authors performed similar calculations constructing hybrid stars with quark matter cores and focusing on the finite-size effects, i.e, a geometric structure of the mixed phase or pasta phase structures with non-vanishing surface tension, which might be relevant in the hadron-quark interface. Also the authors of Ref. [29] investigate the effect that a crystalline quark-hadron mixed phase can have on the neutrino emissivity from the cores of neutron stars considering the formation of spherical blob, rod, and slab rare phase geometries in the mixed phase. In this work we simulate the mixed phase through a continuous interpolation between the hadronic and quark matter phases that mimics the mixing or percolation. This interpolation function simulates the Gibbs construction and it will allow us to compare the results with the obtained considering a sharp phase transition through Maxwell construction.
The paper is organized as follows. In Section 2 we present the NJL models at finite temperature, with the inclusion of the Polyakov loop, and the zero limit equations to model the quark phase of cold NSs. We also present the phase diagram for each model. In Section 3 we analyze some astrophysical applications of the models focusing on the hybrid star EoS. We also describe the hadronic model used to study the hadron-quark phase transition, and we discuss briefly generalities for the occurrence of such transition. In addition, we present the formalism to study isolated NSs oscillation modes and calculate some frequencies, paying especial attention to the g-modes. Moreover we re-examine the universal relationship linked to the Constant Speed of Sound (CSS) parametrization obtained in Ref. [30]. Finally, Section 4 is devoted to summarize the obtained results and discuss their astrophysical relevance.

Quark Matter Description within NJL Models
In its simplest form, the Lagrangian of the NJL model contains local terms of scalar-isoscalar and pseudoscalar-isovectorial interaction between quarks, which are those that reproduce the chiral symmetry breaking dynamics of QCD. The development of this model was done without taking into account the phenomenon of confinement, since it was originally constructed to describe properties of nucleons, which do not exhibit this characteristic. Therefore, the absence of confinement is an important limitation of the theory. However, the lack of confinement can be partially remedied by including a dynamic variable known as Polyakov loop, which is treated as a background field that accounts for gluonic degrees of freedom that can reproduce the deconfinement transition. There are several anzats in the literature for the effective Polyakov-loop potential, which accounts for gauge field self-interactions [31]. In this work we use [32] U where the coefficients a 1 (T, T 0 ) and a 2 (T, T 0 ) are fixed taking into account Lattice QCD simulations of gluon dynamics and the traced Polyakov loop Φ = [2 cos(φ 3 /T) + 1]/3 [33,34]. The quantity φ 3 is related with the color background field, as described in Section 2.1. The effective potential is key to describe the phase transition from the color confined state (T < T 0 , with the minimum of the effective potential being at Φ = 1) to the color deconfined state (T > T 0 , with the minima of the effective potential at Φ = 0), where T 0 is the critical temperature of the deconfinement phase transition. This parameter is the only free parameter of the Polyakov loop once the effective potential is fixed. Therefore, NJL-like models supplemented with the Polyakov loop accounts comfortably for chiral symmetry as well as for the deconfinement-confinement transition in QCD.
In addition to the confinement problem, the NJL model is non-renormalizable. Thus, it is necessary to apply some form of ultraviolet regularization with the introduction of a finite cutoff parameter. Many regularization procedures have been discussed in the literature, such as three and four-momentum cutoffs, Pauli-Villars regularization, proper-time regularization [35]. However, the introduction of a finite cutoff can be problematic to regularize quark loops integrals preserving some physics amplitudes in the meson sector. The way solve the drawbacks of the local NJL is to use non-local rather local interactions. Non-locality arises from several successful approaches to low-energy quark dynamics, such as one-gluon exchange descriptions [36][37][38], the instanton liquid model [39,40], and the Schwinger-Dyson resummation technique [41,42]. The non-local extensions of the NJL model are designed to remove the deficiencies of the local theory.
The general expression for the mean-field thermodynamic potential in NJL-like models can be written as where Ω 0 is defined by the condition that Ω vanishes at T = µ = 0. Once the thermodynamical grand potential is determined, the system's pressure P = −Ω, quark number density n q = ∑ f n f , where f runs over all quark flavors. With these quantities, the EoS for the quark matter (P) = −P + TS + ∑ f µ f n f (where S = ∂P ∂T and n f = ∂P ∂µ f ), can be also calculated.

Local SU(3) PNJL Model
The regularized and free thermodynamic potentials for the SU(3) local NJL model at finite temperature are given by where, α = ψ u ψ u , β = ψ d ψ d and γ = ψ s ψ s are the condensates corresponding to each quark flavor. Since we are consider three quark colors, N c = 3 and E f = p 2 + M 2 f . The color background fields due the coupling to the Polyakov loop are φ c = cφ 3 = n 3 φ 3 , i.e., φ r = −φ g = φ 3 and φ b = 0. The sums over flavor and color indices run over f = (u, d, s) and c = (r, g, b), respectively. The constituent quark masses M f are given by The scalar coupling constant, G s , the 't Hooft coupling constant, H, quark masses and the three-momentum ultraviolet cutoff, Λ, are model parameters. Their values, m u = m d = 5.5 MeV, m s = 140.7 MeV, Λ = 602.3 MeV, G s Λ 2 = 3.67 and HΛ 5 = −12.36, are taken from Ref. [43].
The minimization of the thermodynamic potential with respect to the quark condensates and the Polyakov loop color field φ 3 leads to a system of coupled non-linear equations that can be solved numerically.
The solution of the system equations allow us to construct the phase diagram of the model and to get the thermodynamical quantities to determine the quark matter EoS.

Zero Temperature Limit Including Vector Interaction
At zero temperature, the Polyakov loop vanishes. In this limit, we include a repulsive vector interaction among quarks. The repulsive character of the vector coupling in NJL-like models shifts chiral restoration to larger values of the quark chemical potential in the phase diagram [44,45]. This fact affects the quark-hadron phase transition. Therefore, if the quark deconfined transition in the cores of NSs is modeled through NJL-like models, it is expected that the vector coupling contribution modifies the hybrid EoS and hence the mass-radius relationship of NSs. Vector interactions are crucial to stiffen the quark EoS for the astrophysical application of the NJL model. The regularized thermodynamic grand potential is the same of Equation (2). The free thermodynamic potential is modified at zero temperature where The contribution to the grand potential due vector interaction is given by where the quark number density of flavor f in the mean field approximation is given by The vector coupling constant G v is treated as a free parameter and it is usually described in terms of the strong coupling constant G s . From here on, we will use ξ L v = G v /G s . Vector interaction shifts the quark chemical potential by The system of non-linear equations at zero temperature are solved together with the condition of charge neutrality and baryonic number conservation to obtain the EoS which can be applied to describe quark matter in the interior of NSs.

Non-Local SU(3) Model
The quantities Ω reg and Ω free are given by where S f are auxiliary fields obtained by minimizing the thermodynamic potential, Ω, with respect to the mean field valuesσ f [38,46]. Similar to the local case, these equations plus minimizing Ω with respect to the Polyakov loop color field φ 3 leads to a system of coupled non-linear equations that can be solved numerically to obtain the mean field valuesσ f and the traced Polyakov loop. The momentum dependent constituent quark masses are given by To regulate the non-local interactions we use a Gaussian form factor R(w 2 f nc ) = e

Zero Temperature Limit Including Vector Interaction
As in the local NJL model, at zero temperature we include vector interaction to get a stiffer non-local quark matter EoS for the construction of hybrid stars. The mean-field regularized and free thermodynamic potentials for the SU(3) non-local NJL model can be expressed as The constituent quark masses, M f , are treated as momentum-dependent quantities expressed as where R(ω 2 f ) denotes the Gaussian form factor we use, R(ω 2 f ) = exp −ω 2 f /Λ 2 (for results using different form factors, see for example Ref. [48]) and the quantity ω is affected by the vector interaction contribution [26], which shifts the quark chemical potential whereθ f denotes the vector mean fields. The contribution to the grand potential due vector interaction is given by where ξ NL v = G v /G s . Please note that the quark chemical potential shift in Equation (10) does not affect the non-local form factor, avoiding a recursive problem (see Ref. [26] and references therein).
As in the local model the charge neutrality and baryonic number conservation conditions, together with the minimization of the grand canonical potential with respect to theσ f fields, determine the astrophysical quark matter EoS.

QCD Phase Diagram
According the QCD phase diagram, it is known that the chiral symmetry is spontaneously broken at low temperature, but it is recovered above a certain value. The local order parameter, the quark condensate, is different from zero at low temperature, where the chiral symmetry is broken, and is zero above the chiral phase transition. Lattice QCD simulations suggest that for fermions without mass, the chiral restoration and deconfinement take place at the same temperature, at least in the case of zero chemical potential [49]. In this case, the chiral symmetry restoration temperature is T c 155 − 205 MeV, depending on the number of flavors. As a way to compare both the local and non-local models we have chosen parametrizations in which the values of the light quark condensates at zero chemical potential are comparable (α 1/3 = β 1/3 = −241.9 MeV for the local model and MeV for the non-local model). It is important to note that for moderate quark mass values, the chiral transition does not have a well-defined defined order parameter, and a pure phase transition does not occur but only a rapid change, called crossover, present in the gray line and the labeled "Local" lines of Figure 1. From ab-initio lattice QCD simulations follow that in the QCD phase diagram the transition on the temperature axis is a crossover at zero chemical potential [49]. The QCD phase transition occurring at high temperature and low densities is extremely relevant as the spectrum of primordial GWs changes dramatically if such phase transition is sharp or a smooth crossover. The main reason of these differences are based of the fact that if the primordial hadron-quark phase transition is sharp, bubble nucleation process produce larger perturbations which excite more efficiently the production of GWs [50]. The local NJL, show a crossover phase transition compatible with previous results obtained using different parametrizations of the NJL-like models. The results obtained with the non-local version of the NJL model is a first order phase transition which is qualitatively similar with the result obtained with the standard MIT bag model [50]. In this way, data from primordial GWs could serve to probe the nature of the phase transition and also the existence or not of the QCD critical end point.

Astrophysical Applications of NJL EoS
High-energy astrophysics, especially X-ray and Radio astronomy, and more recently GW astronomy, can be used to shed some light into the nature of matter inside a NS (for details, see Ref. [17] and references therein). There are several works where the NJL model has been employed to describe a possible quark matter phase in NSs (see for example Refs. [51][52][53]. Local and non-local extensions of the NJL model at zero temperature can be used to study hybrid stars with and/or without a quark-hadron mixed phases in their inner cores to explore the role of quark deconfinement in NSs.

Hybrid Stars
Matter in the inner core of NSs is compressed to densities several times higher than the saturation density, ρ 0 . Under these conditions, a phase transition in which nuclei dissolve to form a QGP may occur. When pure QGP appears in the interior of a compact object we will call it a quark-hybrid star. Objects in which the QGP appears forming a mixed phase which coexist with hadronic matter will be denoted as hybrid stars.
In this work, we will construct hybrid and quark hybrid stars, using the NJL models described in the previous section to model the QGP phase in their inner core.
To describe the hadronic phase (outer core) of the hybrid stars we will use the GM1(L) extension of the classical GM1 parametrization in the Relativistic mean field theory (RMF), including both hyperons and Delta isobars [54].

Hadronic EoS
To describe matter in the outer core of our hybrid star models, we will use the relativistic mean field theory taking into account the GM1(L) parametrization. The main change from the classical GM1 model to this improved version is the inclusion of a density-dependent isovector coupling constant for the meson-baryon interaction [54]. The exponential coupling constant introduces an extra parameter that is fixed to fit the experimentally determined slope of the asymmetry energy, L, at saturation density.
The mediators of the interaction are the mesons σ (scalar field), ω (vector field) and ρ (isovector field), with bare masses m σ , m ω and m r ho, respectively. The Lagrangian of this model read as where g σB , g ωB and g ρB (n) are the meson-baryon coupling constants, being the last one density dependent, where n = ∑ B n B is the (total) baryon number density. In this work we consider that the baryons B include the spin 1/2 baryon octet compounded by the nucleons N = (n, p) and hyperons Y = (Λ, Σ − , Σ 0 , Σ + , Ξ − , Ξ 0 ), as well as the spin 3/2 delta isobar quartet (∆ − , ∆ 0 , ∆ + , ∆ ++ ). In addition to the interaction terms meson-baryon, two non-linear terms were included taking into account the meson-meson interaction for the sigma meson. The density dependence of the ρ meson coupling constants is given by This choice of parametrization accounts for nuclear medium effects by making the meson-baryon coupling constants dependent on the local baryon number density [55]. The parameter a ρ is adjusted to satisfy the constraints on the asymmetry energy slope at saturation density L 0 without affecting the other saturation properties, leaving the same values for the rest of the GM1 parameters. The baryon and meson field equations are obtained by evaluating the Euler-Lagrange equations for the fields in Equation (12). Then, applying the RMF approximation to the system, the meson mean-field equations are given by whereσ ≡ σ ,ω ≡ ω ,ρ ≡ ρ and I 3B is the 3-component of isospin. n s B and n B are the scalar and particle number densities for each baryon B, respectively given by where γ B = 2J B + 1 is the corresponding spin degeneration factor, k B is the Fermi momentum, m * B = m B − g σB (n)σ is the effective mass and E B = k 2 B + m * 2 B (σ) is the Fermi energy. Now, to reproduce the NS matter conditions, after including the leptons in the model as a Fermi free gas, electron charge and baryon number must be conserved, according respectively to where the subscripts q B and q l corresponds to baryons and leptons electric charges, respectively, in units of the elementary charge. This constitutes a system of five coupled non-linear equations that are solved simultaneously to determine the meson mean-fields (σ,ω andρ) and the neutron and electron Fermi momenta (k n and k e ). The Fermi momenta of the rest of the baryons are governed by the condition that NS matter be in chemical equilibrium, where µ B is the baryon chemical potential given by R is a rearrangement term, due the density dependence of the coupling constants. The inclusion of Equation (19) is mandatory for thermodynamical consistency [56]. The rearrangement term also affects the baryonic pressure, but it is not present in the energy density, The saturation properties including the nuclear saturation density, n 0 , energy per nucleon, E 0 , nuclear incompressibility, K 0 , effective nucleon mass, m * /m N , asymmetry energy, J, asymmetry energy slope, L 0 , and nucleon potential, U N , and the parameters of the GM1(L) model used in this work are listed in Tables 1 and 2, respectively. Besides the RMF, there are other approaches to construct the EoS that describes hadronic matter. One option is to calculate the hadronic EoS performing microscopic calculations through the Brueckner-Bethe-Goldstone many-body theory with two and three body nuclear interactions [57], which satisfies several requirements of a "realistic" EoS (adequate saturation point for symmetric nuclear matter, appropriates symmetry energy and incompressibility). Effective field theory with and without pions considering effective potentials to describe two-body low-energy observables [58] has also been used to describe nuclear matter and NSs. Alternatively, the new microscopic EoS of dense asymmetric and β−stable nuclear matter at zero temperature presented in Ref. [59] is another option, where two-body and three-body nuclear interactions derived in the context of chiral perturbation theory and Delta isobars were considered. Table 1. Properties of nuclear matter at saturation density for the hadronic GM1(L) parametrization [29,54].

Saturation Property GM1L
65.5 Table 2. Parameters of GM1(L) that produce the properties of symmetric nuclear matter at saturation density given in Table 1.

Phase Transition Formalism
Depending on the value of the hadron-quark surface tension, σ HQ , the phase transition between hadronic and quark matter might be sharp (described via the Maxwell construction, which assumes that σ HQ is infinitely large) or soft (that could be described via the Gibbs construction). On the other hand, "Bulk" Gibbs construction assumes that σ HQ is zero while for intermediate values of σ HQ , the "full" Gibbs formalism must be used [67]. A critical value for σ HQ has been estimated to be ∼5-40 MeV/fm 2 [60,61]. If the hadron-quark surface tension were grater than this value, a sharp phase transition would be favored. On the contrary, a mixed phase in which hadrons and quarks coexist may occur.
In this work we present results using the Maxwell construction. In addition, we mimic mixing and percolation of the Gibbs construction using an interpolation function, I ± (p), between the hadronic, H (p) and quark phases Q (p) where + (−) is used to characterize the function used for pressures greater (smaller) than the transition pressure, p trans , and b can be interpreted as a mixing length [30,[62][63][64][65][66]. The mixing length serves to quantify the extension of the mixed phase inside the compact object. Results obtained using this approach are qualitatively similar to those presented in Ref. [67] in which the authors include finite-size effects on the hadron-quark mixed phase in the context of compact objects. The extension of the mixed phase increases when both σ HQ and b −1 decrease. This fact, although it is not clear (as there are several model-dependent phenomena), seems to indicate the existence of a relationship between these two quantities. Finding and analyzing this possible relationship is well beyond the scope of this work. Considering the mimicking of the mixed phase, the hybrid EoS reads The inclusion of the interpolating function presented in Equation (22) allows us to qualitatively study and analyze how the nature of the phase transition alter the structure and oscillation frequencies of hybrid stars. To study potential observable effects on such quantities, we use two different mixing lengths: b = 0 (Maxwell construction) and b = 2. To calculate the onset of a mixed phase region, we assume that it appears at the pressure where the energy density obtained using Equation (23) differs from the pure hadronic EOS by a 5%. Using a similar approach, we also estimate the appearance of a pure quark core. Alternative methods to describe the hadron-quark phase transition using different interpolation prescriptions have been presented and discussed in Refs. [68,69]. Our results are in general agreement with those obtained by such approximations. In left panel of Figure 2 we present the hybrid EoS constructed using the local NJL quark matter model. With solid lines we show hybrid EoS where the phase transition is constructed using Maxwell criteria for two different values of the vector coupling constant, ξ L v . The dashed lines correspond to the hybrid EoS with a mixed phase mimicked using the interpolating function (22). In the left panel of Figure 3 we present the same but for hybrid EoS constructed using non-local NJL quark model with vector interaction, ξ NL v . The structure of NSs strongly depends on the EoS used to describe matter inside the star. Given a particular EoS, the TOV equations allow to obtain the corresponding family of stationary stellar models. These theoretical curves can be used to discard models that are unable to reproduce astronomical observations [70].
In right panel of Figure 2 we present the M-R relationships obtained using the local NJL quark matter model. In right panels of Figures 2 and 3 we present the M-R relationship obtained using local and non-local NJL models, respectively. The lines colors and styles are the same as used in the corresponding left panels. In both figures the horizontal bars are the measured masses and error bars of pulsars PSR J0348+0432 and PSR J1614-2230. The horizontal arrow represents R 1.4M < 13.76 km the constraint calculated using data of GW170817 [13]. For the b = 2 cases, we have computed the appearance of a mixed phase obtaining an extended branch of hybrid stellar configurations. In all the considered cases these branches spans over a range ∆R∼0.3 km and ∆M∼0.04 M . For these cases, none of the stable stellar configurations present a pure quark matter core. This result is similar to the one obtained using Gibbs formalism to treat the hadron-quark phase transition [25,26].

Stellar Oscillations
There are two distinct families of fluid perturbation modes: odd (or axial), which produce toroidal deformations and even (or polar) which produce spheroidal ones. From an observational point of view, the most promising gravitational waves (GWs) are those related to the fundamental ( f ), the pressure (p) and the gravity (g) fluid modes. The g-modes are particularly interesting as they appear in cold, non-rotating NSs only if a discontinuous EoS is needed to describe matter in their interior [71]. For this reason, a detection of such modes would be a clear indicator of the existence of a sharp hadron-quark phase transition in the inner core of such compact object.
When the metric perturbations are negligible, the equations that govern linear perturbations are greatly simplified. This idea, first applied in Newtonian theory [72], was then extended to study relativistic scenarios [73]. This simplification is called relativistic Cowling's approximation.
Isolated NSs are (in almost all situations) very accurately described using spherically symmetric metrics. In such cases, the equations needed to solve in order to find these frequencies can be written as The functions Φ(r) and Λ(r) characterize the background spacetime and the functions V(r) and W(r), along with the frequency ω, characterize the Lagrangian perturbation vector of the fluid, where Y m (θ, φ) is the m-spherical harmonic [74]. For details related to the physically correct boundary conditions and the numerical procedure used to solve Equations (24), see Ref. [30]. The frequencies obtained, for quadrupolar perturbations, = 2, using the relativistic Cowling's approximation are in general in agreement with those obtained using linearized general relativity showing both qualitatively and quantitative good results for the p and g-modes and only a qualitatively good results for the f -modes with errors ∼15-20% (for details see, [30,75,76], and references therein). We have calculated the frequencies of the fundamental, f , and first pressure, p 1 -modes for each of the EoSs presented in Section 3.1 using both the sharp and mixed phase formalism to construct the phase transition. For the discontinuous EoSs we have also calculated the gravity, g, modes associated with the energy gap, ∆ , between the hadronic and quark phases.
The results for the f and g-modes are presented in the left (right) panel of Figure 4 for the hybrid EoSs constructed using local (non-local) NJL model. For the local (non-local) NJL model, the results of the first pressure are presented in the left (right) panel of Figure 5. The general result is that a detection of f and p 1 modes is not enough to determine nor the quark matter model used nor the nature of the phase transition. Only a simultaneous detection of a g and f -mode would indicate the occurrence of a sharp hadron-quark phase transition in the inner core of NSs. Similar results where obtained in Ref. [30].    In addition, in Figure 6, we present a revisited universal relationship obtained after the inclusion of the additional frequencies associated with the EoSs used in this work. The universal relationship proposed in Ref. [30] between x CSS ≡ ∆ / trans , one of parameters of the Constant Speed of Sound (CSS) parametrization for quark matter [53,77], and the frequency of the g-modes. The original fitting function is y old fit = c old 1 log(x CSS ) + c old 2 with c old 1 = 0.454 ± 0.031 and c old 2 = 0.235 ± 0.023. The root mean square of the residuals of the fit is R old rms = 0.048 showing that the proposed correlation is strong. According the additional frequencies, the relationship suffers only small quantitative changes:

Results and Discussion
In this work we have constructed hybrid EoSs using the GM1(L) parametrization to describe the hadronic phase and local and non-local NJL models including vector interactions to model the quark matter phase. The phase transition is constructed both using the Maxwell construction considering a sharp phase transition and an interpolation function that mimics mixed phase effects for the Gibbs construction. The main effect of the mixed phase is to reduce the maximum mass of the stationary family of stellar configurations. Models in which we consider sharp phase transition show stellar configurations with pure quark matter in their cores, but its appearance inmediately destabilizes the star as it has been shown in several previous works (see, for example, Refs. [53,78,79]). It is worth to mention that in this work we have not included color superconductivity effects that are expected to be relevant at large densities (see, for example, Refs. [80,81] and references therein). Color superconductivity is known to lower the hadron-quark transition pressure giving rise to longer branches of quark hybrid stars (see for example, Ref. [79] and references therein).
In addition, we have studied the QCD phase diagram by compute the phase transition curve for local and non-local NJL models including Polyakov loop. The non-local NJL parametrization used in this work produces a first order phase transition and has no critical end point. This result do not agree with LQCD results. On the contrary, the local model do present a crossover as it is expected, the CEP is located at T CEP = 103.2 MeV and µ CEP = 953.4 MeV. We have considered T 0 = 195 MeV as the critical temperature of the Polyakov loop effective potential. A different anzat of the Polyakov loop potential [44,45] as well as other choices of T 0 could improve the obtained results regarding the critical temperature at zero chemical potential (T c = 164.9 MeV for the non-local model and T c = 198.5 MeV for the local model).
The models developed in this work have been compared to Lattice QCD data calculations and experimental chemical freeze-out results, showing better agreement when non-locality is taken into account. In general, several studies using non-local NJL models have shown a more realistic and consistent description of the strongly interacting matter at extreme conditions of temperature and/or density than the local NJL models (see for example Refs. [33,82,83] and references therein). However, the parametrization of the non-local model used in this work leads to results that are not compatible with the predictions of Lattice QCD [49] regarding the QCD phase diagram. The new parametrization used in Ref. [17] solves completely these discrepancies, as it is shown in Figure 1.
We have also presented, for quadrupolar perturbations, = 2, the frequency modes ( f , p 1 and g, if they exists) for the given EoSs. The detection of such modes is not a powerful tool to discriminate between different NJL models since the differences between the frequencies of the f and p 1 -modes for a given stellar mass is less that ∼1% for the most massive stars. Moreover, unless a g-mode is detected, gravitational waves would not be enough to determine the nature of the phase transition since the effects of the appearance of a mixed phase produce no observable effects on such frequencies. It is expected, as presented in Ref. [30], the frequencies of the calculated g-modes become higher for stars with a color superconducting quark matter phase in their cores. Additionally, we have obtained slight modifications to the fit parameters of the universal relationship between the frequency of the g-modes and the x CSS parameter obtained in [30]. This relationship, and other existing ones (see, for example, [84,85]) could shed some light into the problem of determining NSs parameters (such as mass, radius and internal composition) from GWs observations. The analysis of neutron stars with the presence of deconfined quark mater in their cores, shows that for the NJL models considered in this work a quark matter core is only possible for the last stable star in the M-R diagram. The presence of diquarks can give longer stable hybrid branches [79].