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.

We consider the energy release associated with first-order transition by Gibbs construction and present such energy release as an accumulation of a series of tiny binding energy differences between over-compressed states and stable ones. Universal formulae for the energy release from one homogeneous phase to the other is given. We find the energy release per converted particle varies with number density. As an example, the deconfinement phase transition at supranuclear densities is discussed in detail. The mean energy release per converted baryon is of order 0.1MeV in RMF theory and MIT bag descriptions for hadronic matter and strange quark matter for a wider parameter region.

The deconfinement phase transition from hadronic matter to quark matter can continuously occur during spins down of neutron stars. It will lead to the release of latent heat if the transition is the first-order one. We have investigated the energy release of such deconfinement phase transition for rotating hybrid stars model which include mixed phase of hadronic matter and quark matter. The release of latent heat per baryon is calculated through studying a randomly process of infinitesimal compressing. Finally, we can self-consistently get the heating luminosity of deconfinement phase transition by imputing the EOS of mixed phase, and based on the equation of rotation structure of stars.

As a neutron star spins down, the nuclear matter continuously is converted into quark matter due to the core density increase and then latent heat is released. We have investigated the thermal evolution of neutron stars undergoing such deconfinement phase transition. We have taken into account the conversion in the frame of the general theory of relativity. The released energy has been estimated as a function of change rate of deconfinement baryon number. Numerical solutions to cooling equation are obtained to be very different from the without heating effect. The results show that neutron stars may be heated to higher temperature which is well-matched with pulsar's data despite onset of fast cooling in neutron stars with quark matter core. It is also found that heating effect has magnetic field strength dependence. This feature could be particularly interesting for high temperature of low-field millisecond pulsar at late stage. The high temperature could fit the observed temperature for PSR J0437-4715.

We study the bulk viscosity of strange quark matter (SQM) in a quasiparticle model at finite chemical potential by extrapolating the previous quasiparticle model of finite temperature lattice QCD. The more proper bulk viscosity coefficient can be given in this model where chemical potential $\mu$ and coupling constant $\emph{g}$ are interdependent. We also apply our result to determine the critical rotation of strange stars by r-mode instability window. Our model is compatible to the millisecond pulsar data for a wide range of mass and radius of the stars.

We derive the bulk viscous time scale of neutron stars with quark matter core, i.e. hybrid stars. The r-mode instability windows of the stars show the theoretical result accords with the rapid rotation pulsar data. The fit gives a strong indication for the existence of quark matter in the interior of neutron stars. Hybrid stars instead of neutron or strange stars may result in submillisecond pulsars if they exist.

We study the bulk viscosity of strange quark matter(SQM) in density dependent quark mass model(DDQM) under the background of self-consistent thermodynamics. The correct formulae, with which the viscosity can be evaluated, are derived. We also find that the viscosity in DDQM can be higher by 2 to 3 orders of magnitude than MIT bag model. We calculate the damping time scales due to viscosity coupled to r-modes. The numerical results show the time scale can't be shorter than $10^{-1}$s.