We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations.

We investigate superfluid phase transitions of asymmetric nuclear matter at finite temperature ($T$) and density ($\rho$) with a low proton fraction ($Y_{\rm p} \le 0.2$) which is relevant to the inner crust and outer core of neutron stars. A strong-coupling theory developed for two-component atomic Fermi gases is generalized to the four-component case and is applied to the system of spin-$1/2$ neutrons and protons. The empirical phase shifts of neutron-neutron (nn), proton-proton (pp) and neutron-proton (np) interactions up to $k = 2$ ${\rm fm}^{-1}$ are described by multi-rank separable potentials. We show that (i) the critical temperature of the neutron superfluidity $T_{\rm c}^{\rm nn}$ at $Y_{\rm p}=0$ agrees well with Monte Carlo data at low densities and takes a maximum value $T_{\rm c}^{\rm nn}=1.68$ MeV at $\rho/\rho_0 = 0.14$ with $\rho_0=0.17$ fm$^{-3}$, (ii) the critical temperature of the proton superconductivity $T_{\rm c}^{\rm pp}$ for $Y_{\rm p} \le 0.2$ is substantially suppressed at low densities due to np-pairing fluctuations and starts to dominate over $T_{\rm c}^{\rm nn}$ only above $\rho/\rho_0 = 0.70$ $(0.77)$ for $Y_p =0.1$ $(0.2)$, and (iii) the deuteron condensation temperature $T_{\rm c}^{\rm d}$ is suppressed at $Y_{\rm p}\le 0.2$ due to the large mismatch of the two Fermi surfaces.

A moderately strong vector repulsion between quarks in dense quark matter is needed to explain how a quark core can support neutron stars heavier than two solar masses. We study this repulsion, parametrized by a four-fermion interaction with coupling g_V, in terms of non-perturbative gluon exchange in QCD in the Landau gauge. Matching the energy of quark matter, g_V n_q^2 (where n_q is the number density of quarks) with the quark exchange energy calculated in QCD with a gluon propagator parametrized by a finite gluon mass m_g and a frozen coupling alpha_s, at moderate quark densities, we find that gluon masses m_g in the range 200 - 600 MeV and alpha_s = 2 - 4 lead to a g_V consistent with neutron star phenomenology. Estimating the effects of quark masses and a color-flavor-locked (CFL) pairing gap, we find that g_V can be well approximated by a flavor-symmetric, decreasing function of density. We briefly discuss similar matchings for the isovector repulsion and for the pairing attraction.

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.

If the gravitino mass is in the region from a few GeV to a few 10's GeV, the scalar lepton X such as stau is most likely the next lightest supersymmetry particle. The negatively charged and long-lived X^- may form a Coulomb bound state (A X) with a nucleus A and may affect the big-bang nucleosynthesis through catalyzed fusion process. We calculate a production cross section of Li6 from the catalyzed fusion (He4 X^-) + d \to Li6 + X^- by solving the Schr\"{o}dinger equation exactly for three-body system of He4, d, and X. We utilize the state-of-the-art coupled-channel method, which is known to be very accurate to describe other three-body systems in nuclear and atomic reactions. The importance of the use of appropriate nuclear potential and the exact treatment of the quantum tunneling in the fusion process are emphasized. We find that the astrophysical S-factor at the Gamow peak corresponding to T=10 keV is 0.038 MeV barn. This leads to the Li6 abundance from the catalyzed process as Li6|_{CBBN}\simeq 4.3\times 10^{-11} (D/2.8\times 10^{-5}) ([n_{X^-}/s]/10^{-16}) in the limit of long lifetime of X. Particle physics implication of this result is also discussed.