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

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.

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.