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.

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.