As neutron stars spin-down and contract, the deconfinement phase transition can continue to occur, resulting in energy release(so-called deconfinement heating) in case of the first-order phase transition. The thermal evolution of neutron stars is investigated to combine phase transition and the related energy release self-consistently. We find that the appearance of deconfinement heating during spin-down result in not only the cooling delay but also the increase of surface temperature of stars. For stars characterized by intermediate and weak magnetic field strength, a period of increasing surface temperature could exist. Especially, a sharp jump in surface temperature can be produced as soon as quark matter appears in the core of stars with a weak magnetic field. We think that this may serve as evidence for the existence of deconfinement quark matter. The results show that deconfinement heating facilitates the emergence of such characteristic signature during the thermal evolution process of neutron stars.

We study the cooling of hybrid stars coupling with spin-down. Due to the spin-down of hybrid stars, the interior density continuously increases, different neutrino reactions may be triggered(from the modified Urca process to the quark and nucleon direct Urca process) at different stages of evolution. We calculate the rate of neutrino emissivity of different reactions and simulate the cooling curves of the rotational hybrid stars. The results show the cooling curves of hybrid stars clearly depend on magnetic field if the direct urca reactions occur during the spin-down. Comparing the results of the rotational star model with the transitional static model, we find the cooling behavior of rotational model is more complicated, the temperature of star is higher, especially when direct urca reactions appear in process of rotation. And then we find that the predicted temperatures of some rotating hybrid stars are compatible with the pulsar's data which are contradiction with the results of transitional method.