We study thermal evolution of isolated neutron stars in scalar-tensor theories for the first time. Whether the rapid cooling due to the direct Urca process occurs or not is an interesting question in the viewpoint of the temperature observation of isolated neutron stars. The cooling effect is typically influenced by the proton fraction and the central density. If a fifth force is mediated due to modification of gravity, the relation between the central density and mass of neutron stars differs from one in general relativity, and the cooling curve is also naively expected to be varied. We find that an unscreened fifth force near the surface of neutron stars changes mass-central density relation, and the direct Urca process can be triggered even for neutron stars with smaller mass. This might enable us to test gravitational theories with neutron star cooling.

We investigate X-ray bursts during the thermal evolution of an accreting neutron star which corresponds to the X-ray burster GS\ 1826-24. Physical quantities of the neutron star are included using an equation of state below and above the nuclear matter density. We adopt an equation of state and construct an approximate network that saves the computational time and calculates nuclear energy generation rates accompanying the abundance evolutions. The mass and radius of the neutron star are got by solving the stellar evolution equations from the center to the surface which involve necessary information such as the nuclear energy generation in accreting layers, heating from the crust, and neutrino emissions inside the stellar core. We reproduce the light curve and recurrence time of the X-ray burst from GS 1826-24 within the standard deviation of 1$\sigma$ for the assumed accretion rate, metallicity, and equation of state. It is concluded that the observed recurrence time is consistent with the theoretical model having metallicity of the initial CNO elements $Z_{\rm CNO}$ = 0.01. We suggest that the nuclear reaction rates responsible for the $rp$-process should be examined in detail, because the rates may change the shape of the light curve and our conclusion.

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.