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.

We utilize the bulk viscosity of interacting strange quark matter to reevaluate the damping time scale. The presence of medium effect of bulk viscosity leads to a stronger damping of r-modes, which can be over an order of magnitude for realistic parameters. We find that the r-mode instability window is narrowed due to the medium effect, and hence when a pulsar reaches the instability window it will only slow down by gravitational wave emission to a period of 1.78msec instead of 2.5msec given by early estimate. As a theoretical upper rotation limit of pulsars, the period of 1.78msec is very close to the two most rapidly spinning pulsars known, with periods of about 1.6msec.