Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Haptic Model for the Quantum Phase of Fermions and Bosons in Hilbert Space Based on Knot Theory
Next Article in Special Issue
Four-Fermion Interaction Model on ℳD−1S1
Previous Article in Journal
An Effective Global Optimization Algorithm for Quadratic Programs with Quadratic Constraints
Previous Article in Special Issue
Nonextensive Quasiparticle Description of QCD Matter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Ignacio Francisco Ranea-Sandoval
1,2,*,
Milva Gabriela Orsaria
1,2,
Germán Malfatti
1,2,
Daniela Curin
1,
Mauro Mariani
1,2,
Gustavo Aníbal Contrera
1,3,† and
Octavio Miguel Guilera
4,5
1
Grupo de Gravitación, Astrofísica y Cosmología, Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, 1900 Buenos Aires, Argentina
2
CONICET, Godoy Cruz 2290, 1425 Buenos Aires, Argentina
3
IFLP, UNLP, CONICET, Facultad de Ciencias Exactas, Diagonal 113 entre 63 y 64, 1900 La Plata, Argentina
4
Instituto de Astrofísica de La Plata, CONICET-UNLP, 1900 Buenos Aires, Argentina
5
Instituto de Astrofísica, Pontificia Universidad Catolica de Chile, 8970117 Santiago, Chile
*
Author to whom correspondence should be addressed.
Current address: Department of Physics, San Diego State University, San Diego, CA 92182, USA.
Symmetry 2019, 11(3), 425; https://doi.org/10.3390/sym11030425
Submission received: 28 February 2019 / Revised: 19 March 2019 / Accepted: 20 March 2019 / Published: 22 March 2019
(This article belongs to the Special Issue Nambu-Jona-Lasinio model and its applications)

Abstract

:
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.

1. 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.

2. 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 ( Φ , T ) = 1 2 a 1 ( T , T 0 ) Φ 2 + a 2 ( T , T 0 ) ln ( 1 6 Φ 2 + 8 Φ 3 3 Φ 4 ) T 4 ,
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
Ω = Ω reg + Ω free + Ω 0 + U ( Φ , T ) ,
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 + T S + f μ f n f (where S = P T and n f = P μ f ), can be also calculated.

2.1. Local SU(3) PNJL Model

The regularized and free thermodynamic potentials for the SU(3) local NJL model at finite temperature are given by
Ω reg = 2 N c f 0 Λ d 3 p 2 π 3 E f + G s ( α 2 + β 2 + γ 2 ) + 4 H α β γ Ω free = 2 T f , c 0 d 3 p 2 π 3 { ln 1 + e E f μ f + i c ϕ 3 T + ln 1 + e E f + μ f i c ϕ 3 T } ,
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 f 2 . 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
M f = m f 2 G s ψ f ¯ ψ f 2 H ψ j ¯ ψ j ψ k ¯ ψ k ,
with f , j , k = u , d , s indicate cyclic permutations.
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]. 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
Ω free ( T 0 ) = f = u , d , s 0 p F f d p p 4 E f ,
where p F f = μ f 2 M f 2 . The contribution to the grand potential due vector interaction is given by
Ω v = G v f n f 2
where the quark number density of flavor f in the mean field approximation is given by
n f = N c 3 π 2 [ ( μ f 2 G v n f ) 2 M f 2 ] 3 / 2 .
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 ξ v L = G v / G s . Vector interaction shifts the quark chemical potential by
μ f μ f 2 ξ v L G s n f .
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.

2.2. Non-Local SU(3) Model

The quantities Ω reg and Ω free are given by
Ω reg = 4 T f , c n = 0 d 3 p ( 2 π ) 3 log w f n c 2 + M f 2 ( w f n c 2 ) w f n c 2 + m f 2 1 2 f σ ¯ f S ¯ f + G s 2 S ¯ f 2 + H 2 S ¯ u S ¯ d S ¯ s Ω free = 2 T f , c d 3 p ( 2 π ) 3 ln 1 + e E f μ f i c ϕ 3 T + ln 1 + e E f + μ f + i c ϕ 3 T ,
where E f = p 2 + m f 2 , w f n c 2 = ( w n i μ f + ϕ c ) 2 + p 2 , w n denote Matsubara frequencies. S ¯ f are auxiliary fields obtained by minimizing the thermodynamic potential, Ω , with respect to the mean field values σ ¯ f [38,45]. 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
M f ( w f n c 2 ) = m f + σ ¯ f R ( w f n c 2 ) .
To regulate the non-local interactions we use a Gaussian form factor R ( w f n c 2 ) = e w f n c 2 / Λ 2 , where Λ is relevant for the stiffness of the chiral transition. The up m u and down m d current quarks masses and the coupling constants G s , H, and Λ are chosen so as to reproduce the phenomenological values of the pion decay constant f π = 92.4 MeV and the meson masses m π = 139.0 MeV, m K = 495.0 MeV, m η = 958.0 MeV [38,46], leading to m s = 127.77 MeV, Λ = 780.63 MeV, G s Λ 2 = 14.48 , and H Λ 5 = 267.24 . Light quark current masses are set to m u = m d = 5.50 MeV.

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
Ω reg ( T 0 ) = N c π 3 f = u , d , s 0 d p 0 0 d p p 2 ln ω ^ f 2 + M f 2 ( ω f 2 ) ω f 2 + m f 2 1 2 f = u , d , s ( σ ¯ f S ¯ f + G S 2 S ¯ f 2 ) + H 2 S ¯ u S ¯ d S ¯ s Ω free ( T 0 ) = N c π 2 f = u , d , s 0 μ f 2 m f 2 d p p 2 ( μ f E f ) Θ ( μ f m f )
where N c = 3 , Θ is the Heaviside function, ω f 2 = ( p 0 + i μ f ) 2 + p 2 and E f = p 2 + m f 2 . The constituent quark masses, M f , are treated as momentum-dependent quantities expressed as
M f ( ω f 2 ) = m f + σ ¯ f R ( ω f 2 ) ,
where R ( ω f 2 ) denotes the Gaussian form factor we use, R ( ω f 2 ) = exp ω f 2 / Λ 2 (for results using different form factors, see for example Ref. [47]) and the quantity ω ^ is affected by the vector interaction contribution [26], which shifts the quark chemical potential
μ f μ ^ f = μ f R ( ω f 2 ) θ ¯ f ,
where θ ¯ f denotes the vector mean fields. The contribution to the grand potential due vector interaction is given by
Ω v = f θ ¯ f 2 4 ξ v N L G s ,
where ξ v N L = 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.

2.3. 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 [48]. 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 ψ u ¯ ψ u 1 / 3 = ψ d ¯ ψ d 1 / 3 = 243.8 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 [48].
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 [49]. 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 [49]. 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.

3. 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. [50,51,52]. 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.

3.1. 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 [53].

3.1.1. 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 [53]. 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 h o , respectively. The Lagrangian of this model read as
L = B ψ ¯ B γ μ [ i μ g ω B ω μ g ρ B ( n ) τ · ρ μ ] [ m B g σ B σ ] ψ B + 1 2 ( μ σ μ σ m σ 2 σ 2 ) 1 3 b σ m N [ g σ N σ ] 3 1 4 c σ [ g σ N σ ] 4 1 4 ω μ ν ω μ ν + 1 2 m ω 2 ω μ ω μ + 1 2 m ρ 2 ρ μ · ρ μ 1 4 ρ μ ν · ρ μ ν
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
g ρ B ( n ) = g ρ B ( n 0 ) exp a ρ n n 0 1 ,
This choice of parametrization accounts for nuclear medium effects by making the meson-baryon coupling constants dependent on the local baryon number density [54]. 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
m σ 2 σ ¯ = B g σ B n B s b σ m N g σ N ( g σ N σ ¯ ) 2 c σ g σ N ( g σ N σ ¯ ) 3 , m ω 2 ω ¯ = B g ω B n B , m ρ 2 ρ ¯ = B g ρ B ( n ) I 3 B n B
where σ ¯ σ , ω ¯ ω , ρ ¯ ρ and I 3 B is the 3-component of isospin. n B s and n B are the scalar and particle number densities for each baryon B, respectively given by
n B s ψ ¯ B ψ B = γ B 2 π 2 0 k B m B * ( σ ¯ ) k 2 + m B * 2 ( σ ¯ ) k 2 d k = B γ B 2 π 2 m B * 2 2 k B E B m B * 2 ln E B + k B m B * , n B ψ B ψ B = γ B 2 π 2 0 k B k 2 d k = γ B 2 π 2 k B 3 ,
where γ B = 2 J 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 B 2 + m B * 2 ( σ ¯ ) 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
B q B n B + l q l n l = 0 and B n B n = 0 ,
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,
μ B = μ n q B μ e
where μ B is the baryon chemical potential given by
μ B = k B 2 + m B * 2 ( σ ¯ ) + g ρ B ( n ) ρ ¯ I 3 B + R ˜ .
R ˜ is a rearrangement term,
R ˜ = B g ρ B ( n ) n I 3 B n B ρ ¯ ,
due the density dependence of the coupling constants. The inclusion of Equation (19) is mandatory for thermodynamical consistency [55]. The rearrangement term also affects the baryonic pressure,
P = B γ B 6 π 2 0 k B k 4 d k k 2 + m B * 2 ( σ ¯ ) + l 1 3 π 2 0 k l k 4 d k k 2 + m l 2 ( σ ¯ ) 1 2 m σ 2 σ ¯ 2 + 1 2 m ω 2 ω ¯ 2 + 1 2 m ρ 2 ρ ¯ 2 1 3 b σ m N ( g σ N σ ¯ ) 3 1 4 c σ ( g σ N σ ¯ ) 4 + n R ˜ ,
but it is not present in the energy density,
ϵ = B γ B 2 π 2 0 k B k 2 + m B * 2 ( σ ¯ ) k 2 d k + l 1 π 2 0 k l k 2 + m l 2 ( σ ¯ ) k 2 d k + 1 2 m σ 2 σ ¯ 2 + 1 2 m ω 2 ω ¯ 2 + 1 2 m ρ 2 ρ ¯ 2 + 1 3 b σ m N ( g σ N σ ¯ ) 3 + 1 4 c σ ( g σ N σ ¯ ) 4 .
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 Table 1 and Table 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 [56], 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 [57] 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. [58] is another option, where two-body and three-body nuclear interactions derived in the context of chiral perturbation theory and Delta isobars were considered.

3.1.2. 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 [59]. 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 )
I ± ( p ) = 1 2 1 ± tanh p p trans 10 b p trans ,
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. [59] 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
ϵ MIX ( p ) = ϵ H ( p ) I ( p ) + ϵ Q ( p ) I + ( p ) .
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. [67,68]. 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, ξ v L . 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, ξ v NL .
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 [69].
In right panel of Figure 2 we present the M-R relationships obtained using the local NJL quark matter model. In right panels of Figure 2 and Figure 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.4 M < 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].

3.2. 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 [70]. 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 [71], was then extended to study relativistic scenarios [72]. 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
d W ( r ) d r = d ϵ d P ω 2 r 2 e Λ ( r ) 2 Φ ( r ) V ( r ) + d Φ ( r ) d r W ( r ) ( + 1 ) e Λ ( r ) V ( r ) , d V ( r ) d r = 2 d Φ ( r ) d r V ( r ) 1 r 2 e Λ ( r ) W ( r ) .
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,
ζ i = e Λ ( r ) W ( r ) , V ( r ) θ , V ( r ) sin 2 θ ϕ r 2 Y m ( θ , ϕ ) ,
where Y m ( θ , ϕ ) is the m -spherical harmonic [73]. 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,74,75], 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 [52,76], and the frequency of the g-modes. The original fitting function is y fit old = c 1 old log ( x CSS ) + c 2 old with c 1 old = 0.454 ± 0.031 and c 2 old = 0.235 ± 0.023 . The root mean square of the residuals of the fit is R rms old = 0.048 showing that the proposed correlation is strong. According the additional frequencies, the relationship suffers only small quantitative changes: y fit new = c 1 new log ( x CSS ) + c 2 new with c 1 new = 0.511 ± 0.036 and c 2 new = 0.2935 ± 0.024 . The root mean square of the residuals of the fit is R rms new = 0.065 . The tight fit is still present, making this universal relationship a potentially powerful tool to perform NS Asteroseismology.

4. 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. [52,77,78]). 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. [79,80] 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. [78] 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 C E P = 103.2 MeV and μ C E P = 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] 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,81,82] 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 [48] 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, [83,84]) 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 [78].

Author Contributions

I.F.R.-S. and M.G.O. supervised the project, wrote most of the article with input from the other authors and performed great part of the analysis of the obtained results. G.M. and D.C. performed calculations related to the phase diagrams of NJL models and helped with the analysis of the results. M.M. calculated the hybrid EOSs at zero temperature and obtained the stationary configurations and produced the corresponding figures. G.A.C. wrote the Subsubsection in which the GM1(L) hadronic model is described and parts of Subsection in which the Non-local SU(3) model is presented. Moreover, he produced some of the figures included in the text. O.M.G. calculated oscillation frequencies within the Cowling approximation and produced the corresponding figures. All the authors discussed the results and commented on the final version of the manuscript.

Funding

This research was funded by CONICET grants number PIP-0714 and PIP-0436 and Universidad Nacional de La Plata under grants number G140, G157, X824 and G144. O.M.G. is also supported by the PICT 2016-0053 from ANPCyT, Argentina.

Acknowledgments

The authors want to thank Tomohiro Inagaki, Emilio Elizalde and Dalia Su for inviting us to publish our work in the special issue of the journal Symmetry “Nambu-Jona-Lasinio model and its applications”. We would also like to thank the anonymous referees for their comments and suggestions that help improve the original manuscript. O.M.G. also acknowledges the hosting as invited researcher from IA-PUC. G.A.C. also acknowledges the hosting at the SDSU, USA, within the CONICET-NSF joint research project titled Structure and properties on neutron star cores.

Conflicts of Interest

There are no conflicts of interest regarding this research work.

References

  1. Busza, W.; Rajagopal, K.; Vander Schee, W. Heavy Ion Collisions: The Big Picture and the Big Questions. Annu. Rev. Nucl. Part Sci. 2018, 68, 339–376. [Google Scholar] [CrossRef]
  2. Weber, F. High Energy Physics. In Pulsars as Astrophysical Laboratories for Nuclear and Particle Physics; Cosmology and Gravitation; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  3. Demorest, P.; Pennucci, T.; Ransom, S.; Roberts, M.; Hessels, J. Shapiro Delay Measurement of A Two Solar Mass Neutron Star. Nature 2010, 467, 1081–1083. [Google Scholar] [CrossRef] [PubMed]
  4. Fonseca, E.; Pennucci, T.T.; Ellis, J.A.; Stairs, I.H.; Nice, D.J.; Ransom, S.M.; Demorest, P.B.; Arzoumanian, Z.; Crowter, K.; Dolch, T.; et al. The nanograv nine-year data set: Mass and geometric measurementsof binary millisecond pulsars. Astrophys. J. 2016, 832, 167. [Google Scholar] [CrossRef]
  5. Antoniadis, J.; Freire, P.C.; Wex, N.; Tauris, T.M.; Lynch, R.S.; van Kerkwijk, M.H.; Kramer, M.; Bassa, C.; Dhillon, V.S.; Driebe, T.; et al. A massive pulsar in a compact relativistic binary. Science 2013, 340, 448. [Google Scholar] [CrossRef] [PubMed]
  6. Lynch, R.S.; Boyles, J.; Ransom, S.M.; Stairs, I.H.; Lorimer, D.R.; McLaughlin, M.A.; Hessels, J.W.T.; Kaspi, V.M.; Kondratiev, V.I.; Archibald, A.M.; et al. The Green Bank Telescope 350 MHz Drift-scan Survey II: Data Analysis and the Timing of 10 New Pulsars, Including a Relativistic Binary. Astrophys. J. 2013, 763, 81. [Google Scholar] [CrossRef]
  7. Bodmer, A.R. Collapsed nuclei. Phys. Rev. 1971, D4, 1601–1606. [Google Scholar] [CrossRef]
  8. Witten, E. Cosmic Separation of Phases. Phys. Rev. 1984, D30, 272–285. [Google Scholar] [CrossRef]
  9. Paschalidis, V.; Yagi, K.; Alvarez-Castillo, D.; Blaschke, D.B.; Sedrakian, A. Implications from GW170817 and I-Love-Q relations for relativistic hybrid stars. Phys. Rev. 2018, 97, 084038. [Google Scholar] [CrossRef] [Green Version]
  10. Christian, J.E.; Zacchi, A.; Schaffner-Bielich, J. Signals in the tidal deformability for phase transitions in compact stars with constraints from GW170817. Phys. Rev. 2019, 99, 023009. [Google Scholar] [CrossRef]
  11. Annala, E.; Gorda, T.; Kurkela, A.; Vuorinen, A. Gravitational-wave constraints on the neutron-star-matter Equation of State. Phys. Rev. Lett. 2018, 120, 172703. [Google Scholar] [CrossRef]
  12. Bauswein, A.; Just, O.; Janka, H.T.; Stergioulas, N. Neutron-star Radius Constraints from GW170817 and Future Detections. Astrophys. J. Lett. 2017, 850, L34. [Google Scholar] [CrossRef] [Green Version]
  13. Fattoyev, F.J.; Piekarewicz, J.; Horowitz, C.J. Neutron Skins and Neutron Stars in the Multimessenger Era. Phys. Rev. Lett. 2018, 120, 172702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Raithel, C.; Özel, F.; Psaltis, D. Tidal deformability from GW170817 as a direct probe of the neutron star radius. Astrophys. J. 2018, 857, L23. [Google Scholar] [CrossRef]
  15. Malik, T.; Alam, N.; Fortin, M.; Providência, C.; Agrawal, B.K.; Jha, T.K.; Kumar, B.; Patra, S.K. GW170817: Constraining the nuclear matter equation of state from the neutron star tidal deformability. Phys. Rev. 2018, C98, 035804. [Google Scholar] [CrossRef]
  16. Tews, I.; Margueron, J.; Reddy, S. Critical examination of constraints on the equation of state of dense matter obtained from GW170817. Phys. Rev. C 2018, 98, 045804. [Google Scholar] [CrossRef]
  17. Orsaria, M.G.; Malfatti, G.; Mariani, M.; Ranea-Sandoval, I.F. Phase transitions in neutron stars and their links to gravitational waves. J. Phys. G 2019. in preparation. [Google Scholar]
  18. Chodos, A.; Jaffe, R.L.; Johnson, K.; Thorn, C.B.; Weisskopf, V.F. New extended model of hadrons. Phys. Rev. D 1974, 9, 3471–3495. [Google Scholar] [CrossRef]
  19. Nambu, Y.; Jona-Lasinio, G. Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity. I. Phys. Rev. 1961, 122, 345–358. [Google Scholar] [CrossRef]
  20. Nambu, Y.; Jona-Lasinio, G. Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity. II. Phys. Rev. 1961, 124, 246–254. [Google Scholar] [CrossRef]
  21. Negreiros, R.; Dexheimer, V.A.; Schramm, S. Modeling hybrid stars with an SU(3) nonlinear σ model. Phys. Rev. C 2010, 82, 035803. [Google Scholar] [CrossRef]
  22. Simonov, Y. Vacuum background fields in QCD as a source of confinement. Nucl. Phys. B 1988, 307, 512–530. [Google Scholar] [CrossRef]
  23. Mariani, M.; Orsaria, M.; Vucetich, H. Constant entropy hybrid stars: A first approximation of cooling evolution. Astron. Astrophys. 2017, 601, A21. [Google Scholar] [CrossRef]
  24. Bonanno, L.; Sedrakian, A. Composition and stability of hybrid stars with hyperons and quark color-superconductivity. Astron. Astrophys. 2012, 539, A16. [Google Scholar] [CrossRef]
  25. Orsaria, M.; Rodrigues, H.; Weber, F.; Contrera, G.A. Quark-hybrid matter in the cores of massive neutron stars. Phys. Rev. 2013, D87, 023001. [Google Scholar] [CrossRef]
  26. Orsaria, M.; Rodrigues, H.; Weber, F.; Contrera, G.A. Quark deconfinement in high-mass neutron stars. Phys. Rev. 2014, C89, 015806. [Google Scholar] [CrossRef]
  27. Logoteta, D.; Bombaci, I.; Providência, C.; Vidaña, I. Chiral model approach to quark matter nucleation in neutron stars. Phys. Rev. 2012, 85, 023003. [Google Scholar] [CrossRef]
  28. Logoteta, D.; Providência, C.; Vidaña, I. Formation of hybrid stars from metastable hadronic stars. Phys. Rev. 2013, 88, 055802. [Google Scholar] [CrossRef]
  29. Spinella, W.M.; Weber, F.; Orsaria, M.G.; Contrera, G.A. Neutrino Emissivity in the Quark-Hadron Mixed Phase. Universe 2018, 4, 64. [Google Scholar] [CrossRef]
  30. Ranea-Sandoval, I.F.; Guilera, O.M.; Mariani, M.; Orsaria, M.G. Oscillation modes of hybrid stars within the relativistic Cowling approximation. JCAP 2018, 1812, 031. [Google Scholar] [CrossRef]
  31. Andersen, J.O.; Naylor, W.R.; Tranberg, A. Chiral and deconfinement transitions in a magnetic background using the functional renormalization group with the Polyakov loop. JHEP 2014, 04, 187. [Google Scholar] [CrossRef]
  32. Rößner, S.; Ratti, C.; Weise, W. Polyakov loop, diquarks, and the two-flavor phase diagram. Phys. Rev. D 2007, 75, 034007. [Google Scholar] [CrossRef]
  33. Contrera, G.A.; Orsaria, M.; Scoccola, N.N. Nonlocal Polyakov-Nambu-Jona-Lasinio model with wave function renormalization at finite temperature and chemical potential. Phys. Rev. D 2010, 82, 054026. [Google Scholar] [CrossRef]
  34. Rößner, S.; Hell, T.; Ratti, C.; Weise, W. The chiral and deconfinement crossover transitions: PNJL model beyond mean field. Nucl. Phys. 2008, 814, 118–143. [Google Scholar] [CrossRef] [Green Version]
  35. Klevansky, S.P. The Nambu—Jona-Lasinio model of quantum chromodynamics. Rev. Mod. Phys. 1992, 64, 649–708. [Google Scholar] [CrossRef]
  36. Ripka, G. Quarks Bound by Chiral Fields: The Quark-Structure of the Vacuum and of Light Mesons and Baryons; Oxford University Press: Oxford, UK, 1997. [Google Scholar]
  37. Dumm, D.G.; Grunfeld, A.G.; Scoccola, N.N. Covariant nonlocal chiral quark models with separable interactions. Phys. Rev. D 2006, 74, 054026. [Google Scholar] [CrossRef] [Green Version]
  38. Contrera, G.A.; Dumm, D.G.; Scoccola, N.N. Nonlocal SU(3) chiral quark models at finite temperature: The role of the Polyakov loop. Phys. Lett. B 2008, 661, 113–117. [Google Scholar] [CrossRef] [Green Version]
  39. Dyakonov, D.; Petrov, V. A theory of light quarks in the instanton vacuum. Nucl. Phys. B 1986, 272, 457–489. [Google Scholar] [CrossRef]
  40. Schaefer, T.; Shuryak, E.V. Instantons in QCD. Rev. Mod. Phys. 1998, 70, 323–426. [Google Scholar] [CrossRef]
  41. Roberts, C.D.; Williams, A.G. Dyson-Schwinger equations and their application to hadronic physics. Prog. Part. Nucl. Phys. 1994, 33, 477–575. [Google Scholar] [CrossRef] [Green Version]
  42. Roberts, C.D.; Schmidt, S.M. Dyson-Schwinger equations: Density, temperature and continuum strong QCD. Prog. Part. Nucl. Phys. 2000, 45, S1–S103. [Google Scholar] [CrossRef] [Green Version]
  43. Rehberg, P.; Klevansky, S.P.; Hüfner, J. Hadronization in the SU(3) Nambu–Jona-Lasinio model. Phys. Rev. C 1996, 53, 410–429. [Google Scholar] [CrossRef] [Green Version]
  44. Fukushima, K. Phase diagrams in the three-flavor Nambu-Jona-Lasinio model with the Polyakov loop. Phys. Rev. D 2008, 77, 114028, Erratum in 2008, 78, 039902. [Google Scholar] [CrossRef]
  45. Scarpettini, A.; Gómez, D.D.; Scoccola, N.N. Light pseudoscalar mesons in a nonlocal SU(3) chiral quark model. Phys. Rev. D 2004, 69, 114018. [Google Scholar] [CrossRef]
  46. Contrera, G.A.; Dumm, D.G.; Scoccola, N.N. Meson properties at finite temperature in a three flavor nonlocal chiral quark model with Polyakov loop. Phys. Rev. D 2010, 81, 054005. [Google Scholar] [CrossRef]
  47. Aguilera, D.N.; Blaschke, D.; Grigorian, H.; Scoccola, N.N. Nonlocality effects on color spin locking condensates. Phys. Rev. D 2006, 74, 114005. [Google Scholar] [CrossRef] [Green Version]
  48. Ratti, C. Lattice QCD and heavy ion collisions: A review of recent progress. Rep. Prog. Phys. 2018, 81, 084301. [Google Scholar] [CrossRef] [PubMed]
  49. Schettler, S.; Boeckel, T.; Schaffner-Bielich, J. Imprints of the QCD phase transition on the spectrum of gravitational waves. Phys. Rev. D 2011, 83, 064030. [Google Scholar] [CrossRef]
  50. Lawley, S.; Bentz, W.; Thomas, A. Neutron star properties from an NJL model modified to simulate confinement. Nucl. Phys. B Proc. Suppl. 2005, 141, 29–33. [Google Scholar] [CrossRef] [Green Version]
  51. Contrera, G.A.; Orsaria, M.; Ranea-Sandoval, I.F.; Weber, F. Hybrid Stars in the Framework of NJL Models. Int. J. Mod. Phys. Conf. Ser. 2017, 45, 1760026. [Google Scholar] [CrossRef] [Green Version]
  52. Ranea-Sandoval, I.F.; Han, S.; Orsaria, M.G.; Contrera, G.A.; Weber, F.; Alford, M.G. Constant-sound-speed parametrization for Nambu-Jona-Lasinio models of quark matter in hybrid stars. Phys. Rev. C 2016, 93, 045812. [Google Scholar] [CrossRef]
  53. Spinella, W.M. A Systematic Investigation of Exotic Matter in Neutron Stars. Ph.D. Thesis, Claremont Graduate University, San Diego, CA, USA, 2017. [Google Scholar]
  54. Fuchs, C.; Lenske, H.; Wolter, H.H. Density dependent hadron field theory. Phys. Rev. 1995, C52, 3043–3060. [Google Scholar] [CrossRef]
  55. Hofmann, F.; Keil, C.M.; Lenske, H. Application of the density dependent hadron field theory to neutron star matter. Phys. Rev. C 2001, 64, 025804. [Google Scholar] [CrossRef]
  56. Baldo, M.; Bombaci, I.; Burgio, G.F. Microscopic nuclear equation of state with three-body forces and neutron star structure. Astron. Astrophys. 1997, 328, 274–282. [Google Scholar]
  57. Bombaci, I.; Logoteta, D. Equation of state of dense nuclear matter and neutron star structure from nuclear chiral interactions. Astron. Astrophys. 2018, 609, A128. [Google Scholar] [CrossRef] [Green Version]
  58. Kievsky, A.; Viviani, M.; Logoteta, D.; Bombaci, I.; Girlanda, L. Correlations imposed by the unitary limit between few-nucleon systems, nuclear matter, and neutron stars. Phys. Rev. Lett. 2018, 121, 072701. [Google Scholar] [CrossRef] [PubMed]
  59. Wu, X.H.; Shen, H. Finite-size effects on the hadron-quark phase transition in neutron stars. Nucl. Theory 2017, 96, 025802. [Google Scholar] [CrossRef]
  60. Alford, M.G.; Rajagopal, K.; Reddy, S.; Wilczek, F. The Minimal CFL nuclear interface. Phys. Rev. D 2001, 64, 074017. [Google Scholar] [CrossRef]
  61. Endo, T. Region of hadron-quark mixed phase in hybrid stars. Phys. Rev. C 2011, 83, 068801. [Google Scholar] [CrossRef]
  62. Macher, J.; Schaffner-Bielich, J. Phase transitions in compact stars. Eur. J. Phys. 2005, 26, 341–360. [Google Scholar] [CrossRef]
  63. Masuda, K.; Hatsuda, T.; Takatsuka, T. Hadron-quark crossover and massive hybrid stars. Prog. Theor. Exp. Phys. 2013, 2013, 073D01. [Google Scholar] [CrossRef]
  64. Masuda, K.; Hatsuda, T.; Takatsuka, T. Hadron-quark Crossover and Massive Hybrid Stars with Strangeness. Astrophys. J. 2013, 764, 12. [Google Scholar] [CrossRef]
  65. Álvarez-Castillo, D.; Blaschke, D. Mixed Phase Effects on High-Mass Twin Stars. Phys. Part. Nucl. 2015, 46, 3. [Google Scholar] [CrossRef]
  66. Alford, M.G.; Harris, S.P.; Sachdeva, P.S. On the stability of strange dwarf hybrid stars. Astrophys. J. 2017, 847, 109. [Google Scholar] [CrossRef]
  67. Abgaryan, V.; Alvarez-Castillo, D.; Ayriyan, A.; Blaschke, D.; Grigorian, H. Two Novel Approaches to the Hadron-Quark Mixed Phase in Compact Stars. Universe 2018, 4, 94. [Google Scholar] [CrossRef]
  68. Maslov, K.; Yasutake, N.; Ayriyan, A.; Blaschke, D.; Grigorian, H.; Maruyama, T.; Tatsumi, T.; Voskresensky, D.N. Hybrid equation of state with pasta phases and third family of compact stars I: Pasta phases and effective mixed phase model. arXiv, 2018; arXiv:1812.11889. [Google Scholar]
  69. Özel, F.; Freire, P. Masses, Radii, and the Equation of State of Neutron Stars. Ann. Rev. Astron. Astrophys. 2016, 54, 401–440. [Google Scholar] [CrossRef] [Green Version]
  70. Finn, L.S. G-modes in zero-temperature neutron stars. Mon. Not. R. Astron. Soc. 1987, 227, 265–293. [Google Scholar] [CrossRef]
  71. Cowling, T.G. The non-radial oscillations of polytropic stars. Mon. Not. R. Astron. Soc. 1941, 101, 367. [Google Scholar] [CrossRef]
  72. McDermott, P.N.; van Horn, H.M.; Scholl, J.F. Nonradial g-mode oscillations of warm neutron stars. Astrophys. J. 1983, 268, 837–848. [Google Scholar] [CrossRef]
  73. Sotani, H.; Yasutake, N.; Maruyama, T.; Tatsumi, T. Signatures of hadron-quark mixed phase in gravitational waves. Phys. Rev. D 2011, 83, 024014. [Google Scholar] [CrossRef]
  74. Yoshida, S.; Kojima, Y. Accuracy of the relativistic Cowling approximation in slowly rotating stars. Mon. Not. R. Astron. Soc. 1997, 289, 117–122. [Google Scholar] [CrossRef] [Green Version]
  75. Chirenti, C.; de Souza, G.H.; Kastaun, W. Fundamental oscillation modes of neutron stars: Validity of universal relations. Phys. Rev. D 2015, 91, 044034. [Google Scholar] [CrossRef]
  76. Alford, M.G.; Han, S.; Prakash, M. Generic conditions for stable hybrid stars. Phys. Rev. D 2013, 88, 083013. [Google Scholar] [CrossRef]
  77. Alford, M.G.; Burgio, G.F.; Han, S.; Taranto, G.; Zappalà, D. Constraining and applying a generic high-density equation of state. Phys. Rev. D 2015, 92, 083002. [Google Scholar] [CrossRef] [Green Version]
  78. Ranea-Sandoval, I.F.; Orsaria, M.G.; Han, S.; Weber, F.; Spinella, W.M. Color superconductivity in compact stellar hybrid configurations. Phys. Rev. C 2017, 96, 065807. [Google Scholar] [CrossRef]
  79. Fujihara, T.; Kimura, D.; Inagaki, T.; Kvinikhidze, A. High density quark matter in the Nambu–Jona-Lasinio model with dimensional versus cutoff regularization. Phys. Rev. D 2009, 79, 096008. [Google Scholar] [CrossRef]
  80. Alford, M.G.; Schmitt, A.; Rajagopal, K.; Schäfer, T. Color superconductivity in dense quark matter. Rev. Mod. Phys. 2008, 80, 1455–1515. [Google Scholar] [CrossRef] [Green Version]
  81. Kashiwa, K.; Hell, T.; Weise, W. Nonlocal Polyakov–Nambu–Jona-Lasinio model and imaginary chemical potential. Phys. Rev. D 2011, 84, 056010. [Google Scholar] [CrossRef]
  82. Gómez Dumm, D.; Izzo Villafañe, M.F.; Noguera, S.; Pagura, V.P.; Scoccola, N.N. Strong magnetic fields in nonlocal chiral quark models. Phys. Rev. D 2017, 96, 114012. [Google Scholar] [CrossRef] [Green Version]
  83. Andersson, N.; Kokkotas, K.D. Towards gravitational wave asteroseismology. Mon. Not. R. Astron. Soc. 1998, 299, 1059–1068. [Google Scholar] [CrossRef] [Green Version]
  84. Andersson, N.; Ferrari, V.; Jones, D.I.; Kokkotas, K.D.; Krishnan, B.; Read, J.S.; Rezzolla, L.; Zink, B. Gravitational waves from neutron stars: Promises and challenges. Gen. Relat. Gravit. 2011, 43, 409–436. [Google Scholar] [CrossRef]
Figure 1. Phase diagram of the NJL models considered in this work. As a comparison, in gray we show the results obtained with a non-local NJL parametrization considering m s = 95 MeV. The solid line corresponds to the non-local parametrization with m s = 127.77 MeV. The critical endpoint, CEP, of each phase equilibrium curve, crossover or first order, is denoted by a round dot. The different gray dots indicate some LQCD results and signals of deconfinement from the determination of the chemical freeze-out (see Ref. [17] and references therein).
Figure 1. Phase diagram of the NJL models considered in this work. As a comparison, in gray we show the results obtained with a non-local NJL parametrization considering m s = 95 MeV. The solid line corresponds to the non-local parametrization with m s = 127.77 MeV. The critical endpoint, CEP, of each phase equilibrium curve, crossover or first order, is denoted by a round dot. The different gray dots indicate some LQCD results and signals of deconfinement from the determination of the chemical freeze-out (see Ref. [17] and references therein).
Symmetry 11 00425 g001
Figure 2. Hybrid EoSs (left) and mass-radius (M-R) relationship (right) for the cases corresponding to the local NJL EoSs. In the EoS curves, the constant pressure regions for the b = 0.0 cases correspond to the transition phase. The b = 2.0 cases does not have an abrupt transition so they have a mixed phase region. In the M-R curves, the rounded dot on the solid line indicates where the quark mater core appears and the triangle on the dashed curve indicates the first stellar configuration in which mixed phase is present in its inner core. After the peaks, towards smaller radii, the stars become unstable. The horizontal bars are the measured masses of the 2 M pulsars with their corresponding errors. The horizontal arrow marks the constraint calculated in Ref. [13] for GW170817, R M = 1.4 M < 13.76 km .
Figure 2. Hybrid EoSs (left) and mass-radius (M-R) relationship (right) for the cases corresponding to the local NJL EoSs. In the EoS curves, the constant pressure regions for the b = 0.0 cases correspond to the transition phase. The b = 2.0 cases does not have an abrupt transition so they have a mixed phase region. In the M-R curves, the rounded dot on the solid line indicates where the quark mater core appears and the triangle on the dashed curve indicates the first stellar configuration in which mixed phase is present in its inner core. After the peaks, towards smaller radii, the stars become unstable. The horizontal bars are the measured masses of the 2 M pulsars with their corresponding errors. The horizontal arrow marks the constraint calculated in Ref. [13] for GW170817, R M = 1.4 M < 13.76 km .
Symmetry 11 00425 g002
Figure 3. Same as Figure 2 but for the non-local NJL.
Figure 3. Same as Figure 2 but for the non-local NJL.
Symmetry 11 00425 g003
Figure 4. Frequencies of the f and g-modes (if present) for the hybrid EoSs constructed with the local model are presented in left panel while those obtained using the non-local model are presented in the right one. With solid lines we present models with sharp phase transition and with dashed lines those with a mixed phase. With circle and triangle the frequency of the quark hybrid stars.
Figure 4. Frequencies of the f and g-modes (if present) for the hybrid EoSs constructed with the local model are presented in left panel while those obtained using the non-local model are presented in the right one. With solid lines we present models with sharp phase transition and with dashed lines those with a mixed phase. With circle and triangle the frequency of the quark hybrid stars.
Symmetry 11 00425 g004
Figure 5. Frequencies of the p 1 -modes as a function of the compact object. In the left panel we show those corresponding to hybrid EoSs constructed with the local model. In the right panel we present those obtained using non-local model. With solid lines we present models with sharp phase transition and with dashed lines those with a mixed phase. With circle and triangle the frequency of the quark hybrid stars.
Figure 5. Frequencies of the p 1 -modes as a function of the compact object. In the left panel we show those corresponding to hybrid EoSs constructed with the local model. In the right panel we present those obtained using non-local model. With solid lines we present models with sharp phase transition and with dashed lines those with a mixed phase. With circle and triangle the frequency of the quark hybrid stars.
Symmetry 11 00425 g005
Figure 6. Represented with dots, decimal logarithm of the frequencies of the g-modes obtained in this work as a function of the logarithm of the Δ ϵ / ϵ trans . With red dashed line we present the revised universal relationship between these quantities. Using gray scale the results presented in [30].
Figure 6. Represented with dots, decimal logarithm of the frequencies of the g-modes obtained in this work as a function of the logarithm of the Δ ϵ / ϵ trans . With red dashed line we present the revised universal relationship between these quantities. Using gray scale the results presented in [30].
Symmetry 11 00425 g006
Table 1. Properties of nuclear matter at saturation density for the hadronic GM1(L) parametrization [29,53].
Table 1. Properties of nuclear matter at saturation density for the hadronic GM1(L) parametrization [29,53].
Saturation PropertyGM1L
n 0 (fm 3 )0.153
E 0 (MeV) 16.30
K 0 (MeV)300.0
m * / m N 0.70
J (MeV)32.5
L 0 (MeV)55.0
U N (MeV)65.5
Table 2. Parameters of GM1(L) that produce the properties of symmetric nuclear matter at saturation density given in Table 1.
Table 2. Parameters of GM1(L) that produce the properties of symmetric nuclear matter at saturation density given in Table 1.
ParameterGM1L
m σ (GeV)0.5500
m ω (GeV)0.7830
m ρ (GeV)0.7700
g σ N 9.5722
g ω N 10.6180
g ρ N 8.9830
b σ 0.0029
c σ 0.0011
a ρ 0.3898

Share and Cite

MDPI and ACS Style

Ranea-Sandoval, I.F.; Orsaria, M.G.; Malfatti, G.; Curin, D.; Mariani, M.; Contrera, G.A.; Guilera, O.M. Effects of Hadron-Quark Phase Transitions in Hybrid Stars within the NJL Model. Symmetry 2019, 11, 425. https://doi.org/10.3390/sym11030425

AMA Style

Ranea-Sandoval IF, Orsaria MG, Malfatti G, Curin D, Mariani M, Contrera GA, Guilera OM. Effects of Hadron-Quark Phase Transitions in Hybrid Stars within the NJL Model. Symmetry. 2019; 11(3):425. https://doi.org/10.3390/sym11030425

Chicago/Turabian Style

Ranea-Sandoval, Ignacio Francisco, Milva Gabriela Orsaria, Germán Malfatti, Daniela Curin, Mauro Mariani, Gustavo Aníbal Contrera, and Octavio Miguel Guilera. 2019. "Effects of Hadron-Quark Phase Transitions in Hybrid Stars within the NJL Model" Symmetry 11, no. 3: 425. https://doi.org/10.3390/sym11030425

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

Article Metrics

Back to TopTop