The aim of this work is the study of hybrid stars interiors from the equations of state describing the matter composing them. We used modern hadronic equations of state to describe the matter in the outer core of these stars. The quark phase is modeled by an extended version of the Field Correlator Method, in which vector interactions among quarks and color superconductivity are included. We analyze the effect of these two contributions on the mass-radius relationship and tidal deformability. This is a relevant quantity for neutron stars after the GW170817 event, allowing to impose additional constraints to the equation of state of these compact objetcs.

In the first part of this paper, we investigate the possible existence of a structured hadron-quark mixed phase in the cores of neutron stars. This phase, referred to as the hadron-quark pasta phase, consists of spherical blob, rod, and slab rare phase geometries. Particular emphasis is given to modeling the size othis phase in rotating neutron stars. We use the relativistic mean-field theory to model hadronic matter and the non-local three-flavor Nambu-Jona-Lasinio model to describe quark matter. Based on these models, the hadron-quark pasta phase exists only in very massive neutron stars, whose rotational frequencies are less than around 300 Hz. All other stars are not dense enough to trigger quark deconfinement in their cores. Part two of the paper deals with the quark-hadron composition of hot (proto) neutron star matter. To this end we use a local three-flavor Polyakov-Nambu-Jona-Lasinio model which includes the 't Hooft (quark flavor mixing) term. It is found that this term leads to non-negligible changes in the particle composition of (proto) neutron stars made of hadron-quark matter.

The recent direct observation of gravitational wave event $GW170817$ and its $GRB170817A$ signal has opened up a new window to study neutron stars and heralds a new era of Astronomy referred to as the Multimessenger Astronomy. Both gravitational and electromagnetic waves from a single astrophysical source have been detected for the first time. This combined detection offers an unprecedented opportunity to place constraints on the neutron star matter equation of state. The existence of a possible hadron-quark phase transition in the central regions of neutron stars is associated with the appearance of g-modes, which are extremely important as they could signal the presence of a pure quark matter core in the centers of neutron stars. Observations of g-modes with frequencies between 1 kHz and 1.5 kHz could be interpreted as evidence of a sharp hadron-quark phase transition in the cores of neutron stars. In this article, we shall review the description of the dense matter composing neutron stars, the determination of the equation of state of such matter, and the constraints imposed by astrophysical observations of these fascinating compact objects.

In part one of this paper, we use a non-local extension of the 3-flavor Polyakov-Nambu-Jona-Lasinio model, which takes into account flavor-mixing, momentum dependent quark masses, and vector interactions among quarks, to investigate the possible existence of a spinodal region (determined by the vanishing of the speed of sound) in the QCD phase diagram and determine the temperature and chemical potential of the critical end point. In part two of the paper, we investigate the quark-hadron composition of baryonic matter at zero as well as non-zero temperature. This is of great topical interest for the analysis and interpretation of neutron star merger events such as GW170817. With this in mind, we determine the composition of proto-neutron star matter for entropies and lepton fractions that are typical of such matter. These compositions are used to delineate the evolution of proto-neutron stars to neutron stars in the baryon-mass versus gravitational-mass diagram. The hot stellar models turn out to contain significant fractions of hyperons and $\Delta$-isobars but no deconfined quarks. The latter, are found to exist only in cold neutron stars.

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.

In this work we investigate the effect a crystalline quark-hadron mixed phase can have on the neutrino emissivity from the cores of neutron stars. To this end we use relativistic mean-field equations of state to model hadronic matter and a nonlocal extension of the three-flavor Nambu-Jona-Lasinio model for quark matter. Next we determine the extent of the quark-hadron mixed phase and its crystalline structure using the Glendenning construction, allowing for the formation of spherical blob, rod, and slab rare phase geometries. Finally we calculate the neutrino emissivity due to electron-lattice interactions utilizing the formalism developed for the analogous process in neutron star crusts. We find that the contribution to the neutrino emissivity due to the presence of a crystalline quark-hadron mixed phase is substantial compared to other mechanisms at fairly low temperatures ($\lesssim 10^9$ K) and quark fractions ($\lesssim 30\%$), and that contributions due to lattice vibrations are insignificant compared to static-lattice contributions. There are a number of open issues that need to be addressed in a future study on the neutrino emission rates caused by electron-quark blob bremsstrahlung. Chiefly among them are the role of collective oscillations of matter, electron band structures, and of gaps at the boundaries of the Brillouin zones on bremsstrahlung, as discussed in the summary section of this paper. We hope this paper will stimulate studies addressing these issues.

We study the QCD phase diagram using a non-local SU(3) NJL model with vector interactions among quarks. We analyze several thermodynamic quantities such as entropy and specific heat, and study the influence of vector interactions on the thermodynamic properties of quark matter. Upon imposing electric charge neutrality and baryon number conservation on the field equations, we compute models for the equation of state of the inner cores of proto-neutron stars providing a non-local treatment of quark matter for astrophysics.

In this paper, we discuss the impact of rotation on the particle composition of rotating neutron stars (pulsars). Particular emphasis is put on the formation of quark matter during stellar spin-down, driven by continuous gravitational compression. Our study is based on modern models for the nuclear equation of state whose parameters are tightly constrained by nuclear data, neutron star masses, and the latest estimates of neutron star radii.

In this paper, we use a three flavor non-local Nambu--Jona-Lasinio (NJL) model, an~improved effective model of Quantum Chromodynamics (QCD) at low energies, to investigate the existence of deconfined quarks in the cores of neutron stars. Particular emphasis is put on the possible existence of quark matter in the cores of rotating neutron stars (pulsars). In contrast to non-rotating neutron stars, whose particle compositions do not change with time (are frozen in), the type and structure of the matter in the cores of rotating neutron stars depends on the spin frequencies of these stars, which opens up a possible new window on the nature of matter deep in the cores of neutron stars. Our study shows that, depending on mass and rotational frequency, up to around 8% of the mass of a massive neutron star may be in the mixed quark-hadron phase, if the phase transition is treated as a Gibbs transition. We also find that the gravitational mass at which quark deconfinement occurs in rotating neutron stars varies quadratically with spin frequency, which can be fitted by a simple formula.

We compute models for the equation of state (EoS) of the matter in the cores of hybrid stars. Hadronic matter is treated in the non-linear relativistic mean-field approximation, and quark matter is modeled by three-flavor local and non-local Nambu$-$Jona-Lasinio (NJL) models with repulsive vector interactions. The transition from hadronic to quark matter is constructed by considering either a soft phase transition (Gibbs construction) or a sharp phase transition (Maxwell construction). We find that high-mass neutron stars with masses up to $2.1-2.4 M_\odot$ may contain a mixed phase with hadrons and quarks in their cores, if global charge conservation is imposed via the Gibbs conditions. However, if the Maxwell conditions is considered, the appearance of a pure quark matter core either destabilizes the star immediately (commonly for non-local NJL models) or leads to a very short hybrid star branch in the mass-radius relation (generally for local NJL models).

The discovery of pulsars as heavy as 2 solar masses has led astrophysicists to rethink the core compositions of neutron stars, ruling out many models for the nuclear equations of state (EoS). We explore the hybrid stars that occur when hadronic matter is treated in a relativistic mean-field approximation and quark matter is modeled by three-flavor local and non-local Nambu Jona-Lasinio (NJL) models with repulsive vector interactions. The NJL models typically yield equations of state that feature a first order transition to quark matter. Assuming that the quark-hadron surface tension is high enough to disfavour mixed phases, and restricting to EoSes that allow stars to reach 2 solar masses, we find that the appearance of the quark matter core either destabilizes the star immediately (this is typical for non-local NJL models) or leads to a very short hybrid star branch in the mass-radius relation (this is typical for local NJL models). Using the Constant-Sound-Speed parametrization we can see that the reason for the near-absence of hybrid stars is that the transition pressure is fairly high and the transition is strongly first order.

Numerous theoretical studies using various equation of state models have shown that quark matter may exist at the extreme densities in the cores of high-mass neutron stars. It has also been shown that a phase transition from hadronic matter to quark matter would result in an extended mixed phase region that would segregate phases by net charge to minimize the total energy of the phase, leading to the formation of a crystalline lattice. The existence of quark matter in the core of a neutron star may have significant consequences for its thermal evolution, which for thousands of years is facilitated primarily by neutrino emission. In this work we investigate the effect a crystalline quark-hadron mixed phase can have on the neutrino emissivity from the core. To this end we calculate the equation of state using the relativistic mean-field approximation to model hadronic matter and a nonlocal extension of the three-flavor Nambu-Jona-Lasinio model for quark matter. Next we determine the extent of the quark-hadron mixed phase and its crystalline structure using the Glendenning construction, allowing for the formation of spherical blob, rod, and slab rare phase geometries. Finally we calculate the neutrino emissivity due to electron-lattice interactions utilizing the formalism developed for the analogous process in neutron star crusts. We find that the contribution to the neutrino emissivity due to the presence of a crystalline quark-hadron mixed phase is substantial compared to other mechanisms at fairly low temperatures ($\lesssim 10^9$ K) and quark fractions ($\lesssim 30\%$), and that contributions due to lattice vibrations are insignificant compared to static-lattice contributions.

We study the thermal evolution of neutron stars containing deconfined quark matter in their core. Such objects are generally referred to as quark-hybrid stars. The confined hadronic matter in their core is described in the framework of non-linear relativistic nuclear field theory. For the quark phase we use a non-local extension of the SU(3) Nambu Jona-Lasinio model with vector interactions. The Gibbs condition is used to model phase equilibrium between confined hadronic matter and deconfined quark matter. Our study indicates that high-mass neutron stars may contain between 35 and 40 % deconfined quark-hybrid matter in their cores. Neutron stars with canonical masses of around $1.4\, M_\odot$ would not contain deconfined quark matter. The central proton fractions of the stars are found to be high, enabling them to cool rapidly. Very good agreement with the temperature evolution established for the neutron star in Cassiopeia A (Cas A) is obtained for one of our models (based on the popular NL3 nuclear parametrization), if the protons in the core of our stellar models are strongly paired, the repulsion among the quarks is mildly repulsive, and the mass of Cas A has a canonical value of $1.4\, M_\odot$.

Non-rotating neutron stars are generally treated in theoretical studies as perfect spheres. Such a treatment, however, may not be correct if strong magnetic fields are present (such as for magnetars) and/or the pressure of the matter in the cores of neutron stars is non-isotropic (e.g., color superconducting). In this paper, we investigate the structure of non-spherical neutron stars in the framework of general relativity. Using a parameterized metric to model non-spherical mass distributions, we first derive a stellar structure equation for deformed neutron stars. Numerical investigations of this model equation show that the gravitational masses of deformed neutron stars depend rather strongly on the degree and type (oblate or prolate) of stellar deformation. In particular, we find that the mass of a neutron star increases with increasing oblateness but decreases with increasing prolateness. If this feature carries over to a full two-dimensional treatment of deformed neutron stars, this opens up the possibility that, depending on the type of stellar deformation, there may exist multiple maximum-mass neutron stars for one and for the same model for the nuclear equation of state.

This short review aims at giving a brief overview of the various states of matter that have been suggested to exist in the ultra-dense centers of neutron stars. Particular emphasis is put on the role of quark deconfinement in neutron stars and on the possible existence of compact stars made of absolutely stable strange quark matter (strange stars). Astrophysical phenomena, which distinguish neutron stars from quark stars, are discussed and the question of whether or not quark deconfinement may occur in neutron stars is investigated. Combined with observed astrophysical data, such studies are invaluable to delineate the complex structure of compressed baryonic matter and to put firm constraints on the largely unknown equation of state of such matter.

In this paper, we explore whether or not quark deconfinement may occur in high-mass neutron stars such as J1614-2230 (1.97 \pm 0.04 M_Sun) and J0348+0432 (2.01 \pm 0.04 M_Sun). Our study is based on a non-local extension of the SU(3) Nambu Jona-Lasinio (n3NJL) model with repulsive vector interactions among the quarks. This model goes beyond the frequently used local version of the Nambu Jona-Lasinio (NJL) model by accounting for several key features of QCD which are not part of the local model. Confined hadronic matter is treated in the framework of non-linear relativistic mean field theory. We find that both the local as well as the non-local NJL model predict the existence of extended regions of mixed quark-hadron (quark-hybrid) matter in high-mass neutron stars with masses of 2.1 to 2.4 M_Sun. Pure quark matter in the cores of neutron stars is obtained for certain parametrizations of the hadronic lagrangian and choices of the vector repulsion among quarks. The radii of high-mass neutron stars with quark-hybrid matter and/or pure quark matter cores in their centers are found to lie in the canonical range of 12 to 13 km.

Nonlocal PNJL models allow for a detailed description of chiral quark dynamics with running quark masses and wave function renormalization in accordance with lattice QCD (LQCD) in vacuum. Their generalization to finite temperature T and chemical potential \mu{} allows to reproduce the \mu-dependence of the pseudocritical temperature from LQCD when a nonvanishing vector meson coupling is adjusted. This restricts the region for the critical endpoint in the QCD phase diagram and stiffens the quark matter equation of state (EoS). It is demonstrated that the construction of a hybrid EoS for compact star applications within a two-phase approach employing the nonlocal PNJL EoS and an advanced hadronic EoS leads to the masquerade problem. A density dependence of the vector meson coupling is suggested as a possible solution which can be adjusted in a suitable way to describe hybrid stars with a maximum mass in excess of 2 M_sun with a possible early onset of quark deconfinement even in the cores of typical (M ~ 1.4 M_sun) neutron stars.

Using a nonlocal extension of the SU(3) Nambu-Jona Lasinio model, which reproduces several of the key features of Quantum Chromodynamics, we show that mixed phases of deconfined quarks and confined hadrons (quark-hybrid matter) may exist in the cores of neutron stars as massive as around 2.1 M_Sun. The radii of these objects are found to be in the canonical range of $\sim 12-13$ km. According to our study, the transition to pure quark matter does not occur in stable neutron stars, but is shifted to neutron stars which are unstable against radial oscillations. The implications of our study for the recently discovered, massive neutron star PSR J1614-2230, whose gravitational mass is $1.97 \pm 0.04 M_Sun$, are that this neutron star may contain an extended region of quark-hybrid matter at it center, but no pure quark matter.