Whether fast cooling processes occur or not is crucial for the thermal evolution of neutron stars. In particular, the threshold of the direct Urca process, which is one of the fast cooling processes, is determined by the interior proton fraction $Y_p$, or the nuclear symmetry energy. Since recent observations indicate the small radius of neutron stars, a low value is preferred for the symmetry energy. In this study, simulations of neutron star cooling are performed adopting three models for equation of state (EoS): Togashi, Shen, and LS220 EoSs. The Togashi EoS has been recently constructed with realistic nuclear potentials under finite temperature, and found to account for the small radius of neutron stars. As a result, we find that, since the direct Urca process is forbidden, the neutron star cooling is slow with use of the Togashi EoS. This is because symmetry energy of Togashi EoS is lower than those of other EoSs. Hence, in order to account for observed age and surface temperature of isolated neutron stars (INS) with use of the Togashi EoS, other fast cooling processes are needed regardless of the surface composition.

MAXI J0556-332 is the hottest transient accreting neutron star at the beginning of its quiescence. A theoretical model with crustal heating indicates that an additional shallow heat source of Qshallow > 6 MeV per accreted nucleon is required in the shallow outer crust with respect to the deeper star crust by considering the observed decline in accretion rate at the end of outburst. However, the physical source of this shallow heating is still unclear. In the present investigation, we performed stellar evolutionary calculations, adopting the effects of outburst behavior of the accretion rate. As a consequence, we find that the quiescent light curve of MAXI J0556-332 can be well explained by the nuclear energy generation due to the hot CNO cycle.

Quiescent X-ray luminosities are presented in low mass X-ray binaries with use of evolutionary calculations. The calculated luminosities are compared with observed ones in terms of timeaveraged mass accretion rate. It is shown that neutrino emission by strong pion condensation can explain quiescent X-ray luminosity of SAX J1808.4-3658 and we do not need direct Urca processes concerning nucleons and/or hyperons.

We show a scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star with high effective temperature, and it is consistent with the cooling without exotic phases. The Cas A observation also gives the mass range of $M \geq 1.5 M_\odot$. It may conflict with the current cooling scenarios of compact stars that heavy stars show rapid cooling. We include the effect of the color superconducting (CSC) quark matter phase on the thermal evolution of compact stars. We assume the gap energy of CSC quark phase is large ($\Delta \gtrsim \mathrm{10 MeV}$), and we simulate the cooling of compact stars. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency.

We present a new scenario for the cooling of compact stars considering the central source of Cassiopeia A (Cas A). The Cas A observation shows that the central source is a compact star that has high effective temperature, and it is consistent with the cooling without exotic phases. The observation also gives the mass range of $M \geqslant 1.5 M_\odot$, which may conflict with the current plausible cooling scenario of compact stars. There are some cooled compact stars such as Vela or 3C58, which can be barely explained by the minimal cooling scenario, which includes the neutrino emission by nucleon superfluidity (PBF). Therefore, we invoke the exotic cooling processes, where a heavier star cools faster than lighter one. However, the scenario seems to be inconsistent with the observation of Cas A. Therefore, we present a new cooling scenario to explain the observation of Cas A by constructing models that include a quark color superconducting (CSC) phase with a large energy gap; this phase appears at ultrahigh density region and reduces neutrino emissivity. In our model, a compact star has CSC quark core with a low neutrino emissivity surrounded by high emissivity region made by normal quarks. We present cooling curves obtained from the evolutionary calculations of compact stars: while heavier stars cool slowly, and lighter ones indicate the opposite tendency without considering nucleon superfluidity. Furthermore, we show that our scenario is consistent with the recent observations of the effective temperature of Cas A during the last 10 years, including nucleon superfluidity.

One of the most promising possibilities may be the appearance of quark matter in astrophysical phenomena in the light of recent progress in observations. The mechanism of deconfinement is not well understood, but the thermodynamical aspects of the hadron-quark (HQ) phase transition have been extensively studied in recent years. Then the mixed phase of hadron and quark matter becomes important; the proper treatment is needed to describe the HQ phase transition and derive the equation of state (EOS) for the HQ matter, based on the Gibbs conditions for phase equilibrium. We here adopt a EOS based on the baryon-baryon interactions including hyperons for the hadron phase, while we use rather simple EOS within the MIT bag model in the quark phase. For quark matter we further try to improve the previous EOS by considering other effective models of QCD. One of the interesting consequences may be the appearance of the inhomogeneous structures called "pasta", which are brought about by the surface and the Coulomb interaction effects. We present here a comprehensive review of our recent works about the HQ phase transition in various astrophysical situations: cold catalyzed matter, hot matter and neutrino-trapped matter. We show how the pasta structure becomes unstable by the charge screening of the Coulomb interaction, thermal effect or the neutrino trapping effect. Such inhomogeneous structure may affect astrophysical phenomena through its elasticity or thermal properties. Here we also discuss some implications on supernova explosion, gravitational wave and cooling of compact stars.

The effects of a new triple-$\alpha$ reaction rate (OKK rate) on the helium flash of a helium accreting neutron star in a binary system have been investigated. Since the ignition points determine the properties of a thermonuclear flash of type I X-ray bursts, we examine the cases of different accretion rates, $dM/dt (\dot{M})$, of helium from $3\times10^{-10} M_{\odot} \rm yr^{-1}$ to $3\times10^{-8} M_{\odot} \rm yr^{-1}$, which could cover the observed accretion rates. We find that for the cases of low accretion rates, nuclear burnings are ignited at the helium layers of rather low densities. As a consequence, helium deflagration would be triggered for all cases of lower accretion rate than $\dot{M}\simeq 3\times10^{-8} M_{\odot} \rm yr^{-1}$. We find that OKK rate could be barely consistent with the available observations of the X-ray bursts on the helium accreting neutron star. However this coincidence is found to depend on the properties of crustal heating and the neutron star model.We suggest that OKK rate would be reduced by a factor of $10^{2-3}$ for $10^8$ K in the range of the observational errors.