We review the equation of state of matter in neutron stars from the solid crust through the liquid nuclear matter interior to the quark regime at higher densities. We focus in detail on the question of how quark matter appears in neutron stars, and how it affects the equation of state. After discussing the crust and liquid nuclear matter in the core we briefly review aspects of microscopic quark physics relevant to neutron stars, and quark models of dense matter based on the Nambu--Jona-Lasinio framework, in which gluonic processes are replaced by effective quark interactions. We turn then to describing equations of state useful for interpretation of both electromagnetic and gravitational observations, reviewing the emerging picture of hadron-quark continuity in which hadronic matter turns relatively smoothly, with at most only a weak first order transition, into quark matter with increasing density. We review construction of unified equations of state that interpolate between the reasonably well understood nuclear matter regime at low densities and the quark matter regime at higher densities. The utility of such interpolations is driven by the present inability to calculate the dense matter equation of state in QCD from first principles. As we review, the parameters of effective quark models -- which have direct relevance to the more general structure of the QCD phase diagram of dense and hot matter -- are constrained by neutron star mass and radii measurements, in particular favoring large repulsive density-density and attractive diquark pairing interactions. We describe the structure of neutron stars constructed from the unified equations of states with crossover. Lastly we present the current equations of state -- called "QHC18" for quark-hadron crossover -- in a parametrized form practical for neutron star modeling.

Bulk properties of cold and hot neutron stars (NSs) are studied on the basis of the hadron-quark crossover picture where a smooth transition from the hadronic phase to the quark phase takes place at finite baryon density. By using a phenomenological equation of state (EOS) "CRover" which interpolates the two phases at around 3 times the nuclear matter density, it is found that the cold NSs with the gravitational mass larger than 2-solarmass can be sustained. This is in sharp contrast to the case of the first-order hadron-quark transition. The radii of the cold NSs with the CRover EOS are in the narrow range which is insensitive to the NS masses. Due to the stiffening of the EOS induced by the hadron-quark crossover, the central density of the NSs is at most 4 times the nuclear matter density and the hyperon-mixing barely occurs inside the NS core. This constitutes a solution of the long-standing hyperon puzzle. The effect of color superconductivity (CSC) on the NS structures is also examined with the hadron-quark crossover. For the typical strength of the diquark attraction, a slight softening of the EOS due to two-flavor CSC (2SC) takes place and the maximum mass is reduced by about 0.2-solarmass. The CRover EOS is generalized to the supernova matter at finite temperature to describe the hot NSs at birth. The hadron-quark crossover is found to decrease the central temperature of the hot NSs under isentropic condition. The gravitational energy release and the spin-up rate during the contraction from the hot NS to the cold NS are also estimated.

We construct a new isentropic equation of state (EOS) at finite temperature "CRover" on the basis of the hadron-quark crossover at high density. By using the new EOS, we study the structure of hot neutron stars at birth with the typical lepton fraction ($Y_l=0.3-0.4$) and the typical entropy per baryon ($S=1-2$). Due to the gradual appearance of quark degrees of freedom at high density, the temperature T and the baryon density at the center of the hot neutron stars with the hadron-quark crossover are found to be smaller than those without the crossover by a factor of 2 or more. Typical energy release due to the contraction of a hot neutron star to a cold neutron star with 1.4 solarmass is shown to be about 0.04 solarmass with the spin-up rate about 14%.

On the basis of the percolation picture from the hadronic phase with hyperons to the quark phase with strangeness, we construct a new equation of state (EOS) with the pressure interpolated as a function of the baryon density. The maximum mass of neutron stars can exceed $2M_{\odot}$ if the following two conditions are satisfied; (i) the crossover from the hadronic matter to the quark matter takes place at around three times the normal nuclear matter density, and (ii) the quark matter is strongly interacting in the crossover region and has stiff equation of state. This is in contrast to the conventional approach assuming the first order phase transition in which the EOS becomes always soft due to the presence of the quark matter at high density. Although the choice of the hadronic EOS does not affect the above conclusion on the maximum mass, the three-body force among nucleons and hyperons plays an essential role for the onset of the hyperon mixing and the cooling of neutron stars.

Using the idea of a smooth crossover from the hadronic matter with hyperons to quark matter with strangeness, we show that the maximum mass of neutron stars with quark matter core can be larger than those without quark matter core. This is in contrast to the conventional softening of equation of state due to exotic components at high density. Essential conditions to reach our conclusion are (i) the crossover takes place at relatively low densities, around 3 times the normal nuclear density, and (ii) the quark matter is strongly interacting in the crossover region. By these, the pressure of the system can be greater than that of purely hadronic matter in the crossover region and leads to the maximum mass greater than 2 solar mass. This conclusion is insensitive to the different choice of the hadronic equation of state with hyperons. Several implications of this result to the nuclear incompressibility, the hyperon mixing, and the neutrino cooling are also remarked.

Neutrino processes caused by the neutral current are studied in octet baryon matter. Previous confusion about the baryonic matrix elements of the neutral current interaction is excluded, and a correct table for them improved by consideration of the proton spin problem is presented instead.

With the successful launch of Chandra and XMM/Newton X-ray space missions combined with the lower-energy band observations, we are in the position where careful comparison of neutron star cooling theories with observations will make it possible to distinguish among various competing theories. For instance, the latest theoretical and observational developments already exclude both nucleon and kaon direct URCA cooling. In this way we can now have realistic hope for determining various important properties, such as the composition, degree of superfluidity, the equation of state and steller radius. These developments should help us obtain better insight into the properties of dense matter.