neutron stars provide a natural laboratory to study extremely dense matter . in the interiors of such stars ,
the density can reach up to several times the nuclear saturation density , @xmath8 . at such high densities quarks
could be squeezed out of nucleons to form quark matter  @xcite .
the true ground state of dense quark matter at high densities and low temperatures remains an open problem due to the difficulty of solving non - perturbative quantum chromodynamics ( qcd ) .
it has been suggested that strange quark matter that consists of comparable numbers of @xmath9 , @xmath10 , and @xmath11 quarks may be the stable ground state of normal quark matter  @xcite .
this led to the conjecture that the family of compact stars may have members consisting entirely of quark matter ( so - called strange stars ) and/or members featuring quark cores surrounded by a hadronic shell ( hybrid stars )  @xcite .
observationally , it is very challenging to distinguish the various types of compact objects , such as the strange stars , hybrid stars , and ordinary neutron stars .
their early cooling behavior is dominated by neutrino emission which is a useful probe of the internal composition of compact stars .
thus , cooling simulations provide an effective test of the nature of compact stars  @xcite .
however , many theoretical uncertainties and the current amount of data on the surface temperatures of neutron stars leave sufficient room for speculations  @xcite .
another useful avenue for testing the internal structure and composition of compact stars is astro - seismology , i.e. , the study of the phenomena related to stellar vibrations  @xcite .
in particular , there are a number of instabilities which are associated with the oscillations of rotating stars . here
we will be concerned with the so - called r - mode instability [ see refs .
@xcite for reviews ] .
this instability is known to limit the angular velocity of rapidly rotating compact stars .
the r - mode and related instabilities in rotating neutron stars are damped by the shear and bulk viscosities of matter , therefore these are important ingredients of theoretically modelling rapidly rotating stars .
such models and their microscopic input can then be constrained via the observations of rapidly rotating pulsars , such as the crab pulsar and the millisecond pulsars .    for quark matter in chemical equilibrium , the shear viscosity is dominated by strong interactions between quarks . the bulk viscosity , however , is dominated by flavor - changing weak processes , whereas strong interactions play a secondary role . for normal ( non - superconducting ) strange quark
matter , the bulk viscosity is dominated by the non - leptonic process  @xcite    @xmath12    @xmath13    since the contributions of the leptonic processes @xmath14 and @xmath15 are suppressed due to much smaller phase spaces .
the bulk viscosity of various phases of quark matter has been studied extensively , see refs .
@xcite .
compact stars are strongly magnetized .
neutron star observations indicate that the magnetic field is of the order of @xmath16 g at the surface of ordinary pulsars .
magnetars - strongly magnetized neutron stars - may feature even stronger magnetic fields of the order of @xmath17 g  @xcite .
an upper limit on the magnetic field can be set through the virial theorem .
gravitational equilibrium of stars is compatible with magnetic fields of the order of @xmath18  g  @xcite .
in such a strong magnetic field , not only the thermodynamical but also the hydrodynamical properties of stellar matter will be significantly affected . in particular , due to the large magnetization of strange quark matter the fluid will be strongly anisotropic in a strong magnetic field ( we note here that the magnetization of ordinary neutron matter is small  @xcite ) . therefore , there is need to develop an anisotropic hydrodynamic theory to describe strongly magnetized matter in compact stars .
as we show below , the matter is completely described in terms of eight viscosity coefficients , which include six shear viscosities and two bulk viscosities .    in this paper
, we will carry out a theoretical study of the anisotropic hydrodynamics of magnetized strange quark matter and will calculate the two bulk viscosities .
we will also discuss the implications of the anisotropic bulk viscosities on the r - mode instability in rotating quark stars .
the paper is organized as follows .
the formalism of anisotropic hydrodynamics for magnetized strange quark matter is developed in sec .  [ 2 ] . in sec .  [ 3 ] we apply the local linear response method to derive explicit expressions for bulk viscosities .
the stability of the fluid under strong magnetic field is analyzed in sec .  [ 4 ] .
section  [ 5 ] contains our numerical results for the bulk viscosities .
the damping of the r - mode instability in rotating quark stars by the bulk viscosity is studied in sec .
section  [ 7 ] contains our summary .
we use natural units @xmath19 .
the metric tensor is @xmath20 .
we will use the si system of units in our equations involving electromagnetism , however we will quote the strength of the magnetic field in cgs units ( gauss ) , as is common in the literature on compact stars .
hydrodynamics arises as an effective theory valid in the long - wavelength , low - frequency limit where the energy - momentum tensor @xmath21 , the conserved baryon current @xmath22 , the conserved electric current @xmath23 , and the entropy density flux @xmath24 , _ etc .
_ , are expanded in terms of gradients of the 4-velocity @xmath25 and the thermodynamic parameters of the system , such as the temperature @xmath26 , baryon chemical potential @xmath27 , _ etc._. the hydrodynamic equations can be expressed as conservation laws for the total energy - momentum tensor @xmath21 , as well as baryon and electric currents , @xmath22 and @xmath23 .
the zeroth order terms in the expansion correspond to an ideal fluid and we shall use the index 0 to label them . in the presence of an electromagnetic field , the zeroth - order terms can be generally written as  @xcite , @xmath28 where @xmath29 , @xmath30 , @xmath31 , @xmath32 , and @xmath11 are the local energy density , thermodynamic pressure , baryon number density , electric charge density , and entropy density , respectively measured in the rest frame of the fluid .
@xmath33 is the projector on the directions orthogonal to @xmath25 .
here @xmath34 is the energy - momentum tensor of the electromagnetic field .
@xmath35 is the field - strength tensor which can be decomposed into components parallel and perpendicular to @xmath25 as @xmath36 where in the second line we have introduced the 4-vectors @xmath37 and @xmath38 with @xmath39 being the totally anti - symmetric levi - civita tensor . in the rest frame of the fluid , @xmath40
, we have @xmath41 , @xmath42 and @xmath43 , which are precisely the electric and magnetic fields in this frame . therefore , @xmath44 and @xmath45 are nothing but the electric and magnetic fields measured in the frame where the fluid moves with a velocity @xmath25 .
the antisymmetric tensor @xmath46 is the polarization tensor which describes the response to the applied field strength @xmath35 .
for example , if @xmath47 is the thermodynamic potential of the system , @xmath48 . for later use
, we also define the in - medium field strength tensor @xmath49
. in analogy to @xmath35 we can decompose @xmath46 and @xmath50 as @xmath51 with @xmath52 , @xmath53 , @xmath54 and @xmath55 .    in the rest frame of the fluid ,
the non - trivial components of these tensors are @xmath56 , @xmath57 , @xmath58 , @xmath59 , @xmath60 , and @xmath61 . here
@xmath62 and @xmath63 are the electric polarization vector and magnetization vector , respectively . in the linear approximation
they are related to the fields @xmath64 and @xmath65 by @xmath66 and @xmath67 , with @xmath68 and @xmath69 being the electric and magnetic susceptibilities .
the 4-vectors @xmath70 are all space - like , @xmath71 , and normalized as @xmath72 , where @xmath73 and @xmath74 .
since the electric field is much weaker than the magnetic field in the interior of a neutron star , we will neglect it in most of the following discussion . upon introducing the 4-vector @xmath75 , which is parallel to @xmath45 and is normalized by the condition @xmath76 , and the antisymmetric tensor @xmath77 , we can write @xmath78 with @xmath79 and @xmath80 .
the maxwell equation @xmath81 takes the form @xmath82 its non - relativistic form , which is known as the induction equation , is given by @xmath83 where @xmath84 is the 3-velocity of the fluid .
contracting eq .
( [ maxwell0 ] ) with @xmath85 gives @xmath86 where @xmath87 and @xmath88 .
the second maxwell equation can be written as @xmath89 whose non - relativistic form is @xmath90 where @xmath91 is the electric charge density and @xmath92 is the corresponding current .
it is useful to re - write the energy - momentum tensor in the following form  @xcite , @xmath93 where @xmath94 is the projection tensor on the direction perpendicular to both @xmath25 and @xmath95 .
we have defined the transverse and longitudinal pressures @xmath96 and @xmath97 relative to @xmath95 ; here @xmath30 is the thermodynamic pressure . in the absence of a magnetic field , the fluid is isotropic and @xmath98 . in the local rest frame of fluid
, we have @xmath99 ( without loss of generality , we choose the @xmath100-axis along the direction of the magnetic field ) , hence the electromagnetic tensor takes the usual form , while @xmath101 .
next we would like to check the consistency of the terms that appear in @xmath102 with the formulae of standard thermodynamics involving electromagnetic fields .
by using the thermodynamic relation @xmath103 and the conservation equations for @xmath104 , @xmath105 , and @xmath106 in ideal hydrodynamics , one can show that the hydrodynamic equation @xmath107 together with the maxwell equation ( [ maxwell ] ) implies @xmath108 which is consistent with the standard thermodynamic relation @xmath109 one should note that the potential energy @xmath110 has already been included in our definition of @xmath29 .
otherwise , new terms @xmath110 , @xmath111 and @xmath112 should be added to the left - hand sides of eq .
( [ thermo ] ) , eq .
( [ de ] ) , and eq .  ( [ de2 ] ) , respectively .
thus , we conclude that our hydrodynamical equations are consistent with well - known thermodynamic relations .      by keeping the first - order terms of the derivative expansion of conserved quantities one obtains the navier - stokes - fourier - ohm theory . in this theory , @xmath21 , @xmath22 , @xmath23 , and @xmath24 can be generally expressed as @xmath113 where @xmath114 , and @xmath115 are the dissipative fluxes .
they all are orthogonal to @xmath25 ; this reflects the fact that the dissipation in the fluid should be spatial .
we shall assume that @xmath115 can be expressed as a linear combination of @xmath116 , and @xmath117  @xcite .
this allows us to incorporate the fact that the entropy flux is determined by the energy - momentum and baryon number diffusion fluxes .
thus , @xmath118 with the coefficients @xmath119 , @xmath120 , and @xmath121 being functions of thermodynamic variables
.    next , the hydrodynamic equations are specified by utilizing the conservation laws of the total energy - momentum @xmath21 , the baryon number density flow @xmath22 , electric current @xmath23 , and the second law of thermodynamics , @xmath122    to discuss the dissipative parts , let us first define the 4-velocity @xmath25 , since it is not unique when energy exchange by thermal conduction is allowed for .
we will use the landau - lifshitz frame in which @xmath25 is chosen to be parallel to the energy density flow , so that @xmath123 . upon projecting the first equation of eq .
( [ hydroeqn ] ) on @xmath124 and after some straightforward manipulations , we find @xmath125 combining eq .
( [ energyeqn ] ) , eq .
( [ thermo ] ) , and the second equation in eq .
( [ hydroeqn ] ) , we arrive at @xmath126 where @xmath127 and @xmath128 . for a thermodynamically and hydrodynamically stable system ,
( [ divs ] ) should be non - negative .
this implies @xmath129 where @xmath130 , @xmath131 is the rank - four tensor of viscosity coefficients , and @xmath132 and @xmath133 are thermal and electrical conductivity tensors with respect to the diffusion fluxes of baryon number density and electric charge density . by definition ,
@xmath131 is symmetric in the pairs of indices @xmath134 and @xmath135 .
it necessarily satisfies the condition @xmath136 , which is onsager s symmetry principle for transport coefficients .
similarly , the tensors @xmath132 and @xmath133 should satisfy the conditions @xmath137 and @xmath138 . furthermore , all the tensors of transport coefficients @xmath131 , @xmath132 , and @xmath133 must be orthogonal to @xmath25 by definition .
as we have seen , the appearance of the magnetic field makes the system anisotropic .
such anisotropy is specified by the vector @xmath95 , so that the tensors @xmath131 , @xmath133 , and @xmath132 should be in general expressed in terms of @xmath25 , @xmath95 , @xmath139 , and @xmath140 .
all independent irreducible tensor combinations having the symmetry of @xmath131 and which are orthogonal to @xmath25 are  @xcite(see also appendix [ proof ] ) @xmath141 all independent irreducible tensor combinations having the symmetry of @xmath132 and @xmath133 and which are orthogonal to @xmath25 are @xmath142 in accordance with the number of tensors ( [ combination ] ) and ( [ combination - k ] ) , a fluid in a magnetic field in general has eight independent viscosity coefficients , three independent thermal conduction coefficients and three independent electrical conductivities .
they may be defined as the coefficients in the following decompositions for the viscous stress tensor , heat flux , and electric charge flux    @xmath143    @xmath144    @xmath145    where @xmath146 , @xmath147 and @xmath148 is constructed so that the @xmath149 s are the coefficients of its traceless parts , _
i.e. _ , they can be regarded as shear viscosities ; @xmath150 s are the coefficients of the parts with non - zero trace and can be considered as bulk viscosities . the @xmath151 s and @xmath152 s are thermal and electrical conductivities , respectively .
now the divergence of entropy density flux ( [ divs ] ) can be explicitly written as    @xmath153    one should note that the terms corresponding to the transport coefficients @xmath154 , @xmath155 , and @xmath156 in eqs .
( [ pmn])-([jbe ] ) do not contribute to the divergence of the entropy density flux .
for stable systems , all the other transport coefficients must be positive definite according to the second law of thermodynamics . in sec.[4
] we will demonstrate explicitly that negative bulk viscosities @xmath0 or / and @xmath1 indeed cause an instability in the hydrodynamic evolution of strange stars .    to conclude this section
, we compare our definition of the viscosity coefficients in eq .
( [ pmn ] ) with the definition given in ref .
@xcite for non - relativistic fluid , which reads    @xmath157    where @xmath158 and the remaining notations are self - explanatory .
our viscosity coefficients in eq .
( [ pmn ] ) are related to the coefficients in eq .
( [ pij ] ) by @xmath159 in ref .
@xcite there is no term that corresponds to our @xmath160 .
the reason is that ref .
@xcite considers the combination of vectors ( viii ) in eq .
( [ combination ] ) as dependent on the others ; in the appendix we will show that ( at least for relativistic fluids ) all the combinations ( i)-(viii ) are linearly independent .
note that the transport coefficients in eq .
( [ divs2 ] ) appear as prefactors of quadratic forms , therefore the second law of thermodynamics requires that these coefficients must be positive definite for stable ensembles .
this is not manifest in eq .
( [ pij ] ) .
the typical oscillation frequency of neutron stars is of the order of magnitude of the rotation frequency , @xmath161s@xmath162s@xmath163 .
the most important microscopic processes which dissipate energy on the corresponding timescales are the weak processes .
the compression and expansion of strange quark matter with nonzero strange quark mass will drive the system out of equilibrium .
the processes ( [ weak1 ] ) and ( [ weak2 ] ) are the most efficient microscopic processes that restore local chemical equilibrium .
therefore , the bulk viscosities are determined mainly by the processes ( [ weak1 ] ) and ( [ weak2 ] ) . in this section
we will derive analytical expressions for the bulk viscosities @xmath0 and @xmath1 . , these processes involve quark - quark scattering via an interchange of a single gluon  @xcite . however , these processes are too fast to contribute to the bulk viscosity at the frequency of the order of the stellar oscillation frequency . ]
let us imagine an isotropic flow @xmath164 which characterizes the stellar oscillation .
if there are no dissipative processes , such an oscillation will drive the system from one instantaneous equilibrium state to another instantaneous equilibrium state .
the appearance of dissipation changes the picture : during the oscillations the thermodynamic quantities will differ from their equilibrium values .
let us explore how the thermodynamic quantities evolve during the flow oscillation .    in general ,
we can write the change of baryon density @xmath165 induced by the oscillation of the fluid as @xmath166 where @xmath167 is the static ( time - independent ) equilibrium value , @xmath168 denotes the equilibrium value shift from @xmath167 due to the volume change and @xmath169 denotes the instantaneous departure from the equilibrium value . because processes ( [ weak1 ] ) and ( [ weak2 ] ) conserve baryon number , @xmath170 can be set to zero , if we neglect other microscopic processes .
then @xmath171 can be determined through the continuity equation of ideal hydrodynamics , @xmath172 since processes ( [ weak1 ] ) and ( [ weak2 ] ) also conserve the sum @xmath173 , a similar argument leads to the relation @xmath174    when the system is driven out of chemical equilibrium , the chemical potential of the @xmath11 quark will be slightly different from that of the @xmath10 quark .
let us denote this difference by @xmath175 , with @xmath176 being the deviation of @xmath177 from its static equilibrium value .
up to linear order in the deviation we find @xmath178 where @xmath179 denotes the number density of quarks of flavor @xmath180 and the subscript 0 indicates that the quantity in the bracket is computed in static equilibrium state .
@xmath181 and @xmath182 are the deviations of @xmath11-quark and @xmath10-quark densities from their static equilibrium value . in the final expressions they should be functions of @xmath183 .
the instantaneous departure from equilibrium is restored by the weak processes ( [ weak1 ] ) and ( [ weak2 ] ) . adopting the linear approximation ,
this can be described by @xmath184 where @xmath185 and @xmath186 are the rates of processes ( [ weak1 ] ) and ( [ weak2 ] ) , respectively .
if the weak processes are turned off , one should have @xmath187 where the dot denotes the time derivative .
after turning on the weak processes , we have @xmath188 this system of coupled linear first - order equations is closed by substituting eqs .
( [ deltanb])-([deltamu ] ) .
it is then easy to obtain the solution , @xmath189 where the coefficient @xmath190 is defined by @xmath191    the parallel and transverse components of the pressure @xmath192 and @xmath193 can be written as @xmath194 and @xmath195 where @xmath196 the small departure of the magnetic field @xmath197 can be calculated by the variation of eq .
( [ maxwell ] ) .
one finds @xmath198 a direct calculation then gives , @xmath199 where we introduced the coefficient @xmath200 as @xmath201 now we obtain @xmath202 with @xmath203 defined as @xmath204    then the deviation of @xmath21 from its equilibrium value can be written as @xmath205 for isotropic flows , we have @xmath206 by comparing the above two expressions , we obtain @xmath207 and @xmath208 expressions ( [ vis1 ] ) and ( [ vis2 ] ) show that the bulk viscosities @xmath0 and @xmath1 are functions of the perturbation frequency @xmath209 , the weak rate @xmath210 and the thermodynamic quantities @xmath211 . from the derivation above we can convince ourselves that these expressions should be valid also in the case of color - superconducting matter .
for zero magnetic field ,
( [ vis2 ] ) reduces to eq .
( [ vis1 ] ) which , with parameters @xmath212 , and @xmath210 taken in the absence of magnetic field gives the expression for the usual bulk viscosity @xmath213 defined in isotropic hydrodynamics .
both @xmath0 and @xmath1 attain their maxima in the limit of zero frequency , @xmath214 ; and the maxima are inversely proportional to the weak interaction rate .
at high frequency , @xmath215 , @xmath0 and @xmath1 fall off as @xmath216 . for practical applications to cold strange stars , where the chemical potential is much larger than the temperature , the quantities @xmath217 , and @xmath190 can be evaluated in the zero - temperature limit .
their dependence on temperature is weak .
contrary to this , the coefficient @xmath210 depends strongly on temperature : for normal quark matter , @xmath210 has a power - law dependence on @xmath26 ( see sec.[5 ] ) ; for a fully paired color - superconducting phase , the weak rate is exponentially suppressed by a boltzmann factor @xmath218 with @xmath219 being the superconducting gap .
consequently , @xmath0 and @xmath1 depend exponentially on @xmath26  @xcite .
before we find the numerical values for the bulk viscosities , we first need to analyze the stability of magnetized strange quark matter .
we observe that according to eq .
( [ vis2 ] ) negative values of @xmath0 are a priori not excluded . in the following
we analyze the consequences and implications of negative @xmath0 on the stability of the system .
stable equilibrium in a self - gravitating fluid , such as in strange quark stars , is attained through the balance of gravity and pressure .
the gravitational equilibrium requires that both components of the pressure @xmath192 and @xmath193 should be positive ( otherwise the star will undergo a gravitational collapse ) . at zero temperature , the one - loop thermodynamic pressure @xmath220 of non - interacting strange quark matter , where @xmath221 is the grand partition function , is given by @xmath222,{\nonumber\\}\end{aligned}\ ] ] where @xmath223 is the absolute value of electric charge , @xmath177 and @xmath224 are the chemical potential and the mass of quark of flavor @xmath180 , @xmath225 labels the landau levels
, @xmath226 is the degree of degeneracy of each landau level and @xmath227 $ ] is the highest landau level for quarks of flavor @xmath180 . by differentiating eq .
( [ pstrong ] ) with respect to @xmath2 one can easily get the magnetization as @xmath228.{\nonumber\\}\end{aligned}\ ] ] in fig .
[ mag_u ] we illustrate the magnetization as a function of magnetic field at zero temperature .
the parameters are chosen as @xmath229 on average , the magnetization increases when @xmath2 grows and eventually becomes constant when @xmath230 , where @xmath231 however , the detailed structure of the magnetization exhibits strong de haas - van alphen oscillations  @xcite .
this oscillatory behavior is of the same origin as the de haas - van alphen oscillations of the magnetization in metals and originates from the quantization of the energy levels associated with the orbital motion of charged particles in a magnetic field .
the irregularity of this oscillation shown in fig .
[ mag_u ] is due to the unequal masses and charges of @xmath232 , and @xmath11 quarks .
when @xmath230 , all quarks are confined to their lowest landau level and their transverse motions are frozen . in this case , the longitudinal pressure @xmath233 , so the magnetization @xmath234 is independent of @xmath2 , and the transverse pressure @xmath96 of the system vanishes .
this behavior is evident in fig .
[ pressure_u ] .
.,width=283 ]     ( dashed , black online ) and transverse @xmath193 ( solid , red online ) pressures of strange quark matter as functions of the magnetic field @xmath2 in units of the pressure @xmath235 for zero magnetic field.,width=283 ]    thus , we conclude that when all the quarks are confined to their lowest landau level , the transverse pressure vanishes .
the system therefore becomes mechanically unstable and would collapse due to the gravity  @xcite .
this phenomenon establishes an upper limit on the magnetic field sustained by a quark star .
given @xmath236 mev , this upper limit is roughly @xmath237 g as shown in fig .
[ pressure_u ] .
the thermodynamic stability requires that the local entropy density should reach its maximum in the equilibrium state  @xcite .
the total energy density of the fluid and the magnetic field is @xmath238 and the corresponding first law of thermodynamics , in variational form , is @xmath239 where @xmath240 is the strength of the magnetic field and @xmath241 stands for a small departure of a given quantity from its equilibrium value .
varying eq .
( [ dtve ] ) on both sides and taking into account that @xmath242 , and @xmath2 are independent variational variables , one obtains @xmath243 where @xmath244 and @xmath245 the thermodynamical stability criteria require that @xmath246 , or , equivalently , @xmath247 is positive definite .
taking into account the relation @xmath248 , it is easy to show that these criteria are equivalent to the requirement @xmath249 from eq .
( [ pstrong ] ) we obtain @xmath250 then it is evident that the condition @xmath251 is always satisfied .
however , @xmath252 is divergent when @xmath2 approaches a @xmath253 landau level for each flavor quark from below , @xmath254 this shows that strange quark matter will be thermodynamically unstable just below each landau level @xmath255 ( @xmath253 ) for every flavor @xmath180 .
the first three thermodynamically unstable windows ( tuw ) associated , respectively , with @xmath256 , and @xmath257 are illustrated in the @xmath258 plane in fig .
[ tiw ] for our parameters ( [ paras ] ) .
the tuw is actually very narrow .
one may conjecture that such an instability may lead to formation of magnetic domains  @xcite .
the presence of possible magnetic domains in neutron star crusts was discussed in ref .
@xcite ; furthermore , such a possibility for color - flavor - locked quark matter was pointed out in ref .
we will not pursue here the study of domain structure and related physics , since among other things , this will require us to specify the geometry of the system
.     plane for @xmath259 and @xmath260 .
, width=283 ]      in this subsection we address the problem of hydrodynamic stability within the theory presented in sec.[2 ] .
the fluid is said to be stable if it returns to its initial state after a transient perturbation . otherwise , _
i.e. _ , when the perturbation grows and takes the fluid into another state , the fluid is unstable .
our particular goal here is to determine whether a small , plane - wave perturbation around a homogeneous equilibrium state grows for nonzero @xmath0 and @xmath1  @xcite .
all the other transport coefficients are set to zero .
for this purpose , it is sufficient to solve the hydrodynamic and maxwell equations that are linearized around the homogeneous equilibrium state @xmath261 where    @xmath262    upon linearizing the normalization conditions @xmath263 , one finds that the perturbed variables need to satisfy the constraints @xmath264 in the equations above the perturbations are assumed to have the form @xmath265 , where @xmath266 is constant , and the unperturbed quantities are independent of space and time .
the hydrodynamic and maxwell equations need to be supplemented by an equation of state ( eos ) in order to close the system .
the linearized eos is given by @xmath267 where @xmath268 is the speed of sound .
in the most general case eqs .
( [ linear])-([eos ] ) constitute 15 independent equations , but the equations associated with the conservation of @xmath31 and @xmath32 are decoupled from the others if the eos is taken in the form ( [ eos ] ) .
therefore , we are left with 13 equations .
we work in the rest frame of the equilibrium fluid , @xmath269 and @xmath99 and choose as independent variables @xmath270 ( n.b .
one can choose other independent variables , but the results do not change ) .
the thirteen linear equations can be collected into the following matrix form @xmath271 the matrix @xmath272 has the following form @xmath273 where @xmath274 is the total enthalpy .
the exponential plane - wave solutions for frequencies @xmath275 and wave - vectors @xmath276 satisfy the dispersion relations given by @xmath277 for the modes propagating parallel ( longitudinal modes ) or perpendicular ( transverse modes ) to the magnetic field , eq .
( [ detg ] ) has simple solutions .
\1 ) transverse modes , @xmath278 .
there are two types of transverse modes .
one is solely determined by the maxwell equations and has the following dispersion relation @xmath279 where @xmath280 , and describes simply an electromagnetic wave .
another solution has the dispersion relation @xmath281 where the second approximate relation is valid in the long - wavelength limit .
this solution represents a sound wave propagating perpendicular to the magnetic field .
the speed of this sonic wave is @xmath282 and is smaller than the speed @xmath283 of a ordinary sound wave .
it is seen that positive @xmath0 implies dissipation of the sonic wave , _
i.e. _ , a decay of the initial disturbance .
we conclude that the fluid flow is stable in the case .
however , we see that for negative @xmath0 , the initial disturbance grows and the fluid is unstable .
thus , we conclude that negative transverse bulk viscosity implies hydrodynamic instability via growth of transverse sound waves .
\2 ) longitudinal modes , @xmath284 .
we find three types of longitudinal modes .
the first one is again the electromagnetic wave with the dispersion relation @xmath285 the second one is a transverse wave oscillating perpendicularly to the magnetic field , but traveling along the magnetic field lines .
it has the dispersion relation @xmath286 this mode is the alfven wave whose speed is equal @xmath287 .
the third longitudinal mode has the following dispersion relation @xmath288 this mode represents an ordinary sound wave with dissipation due to the longitudinal bulk viscosity @xmath1 .
it is obvious that if @xmath289 this mode will not decay , rather grow , thus leading to hydrodynamic instability .    in the next section
, we will show that for certain values of the parameters , the transverse bulk viscosity @xmath0 could be indeed negative .
we emphasize here that this does not imply a violation of the second law of thermodynamics , rather this manifests a hydrodynamic instability of the ground state , _
i.e. _ , small perturbations will take the system via this hydrodynamic instability to a new state .
a candidate state is the one which has inhomogeneous ( domain ) structure .
both the structure of the new state and the transition from the homogeneous to the inhomogeneous state are interesting problems which are beyond the scope of this study .
however , we would like to point out a number analogous cases where a negative transport coefficient indicates instability towards formation of a new state with domain structure .
one such case is the negative resistivity ( also known as the gunn effect ) in certain semiconducting materials  @xcite .
another case is the negative ( effective ) shear viscosity , which is extensively studied in the literature  @xcite .
finally , negative bulk viscosity has been investigated in different contexts in  @xcite .
in order to calculate the bulk viscosities , we need to determine the coefficients @xmath290 and @xmath210 . from eqs .
( [ pstrong ] ) , ( [ m0 ] ) , and ( [ nf0 ] ) we obtain in a straightforward manner @xmath291 where @xmath292    the rate @xmath210 of the weak processes ( [ weak1 ] ) and ( [ weak2 ] ) should also be affected by a strong magnetic field .
the major effect of a magnetic field on @xmath210 is to modify the phase space of weak processes ( [ weak1 ] ) and ( [ weak2 ] )  @xcite . taking this into account ,
one obtains @xmath293 where @xmath294 is the fermi constant .
\1 ) when the magnetic field is much smaller than the typical chemical potential , say , @xmath295 , its effect on the bulk viscosities is negligible . for typical parameters ( [ paras ] )
, this condition holds up to @xmath296 g. in this case , the system is practically isotropic , @xmath1 and @xmath0 are effectively degenerate with the isotropic @xmath213 , the bulk viscosity of unmagnetized matter .
the zero - magnetic field limit results can be obtained easily by replacing @xmath297 and @xmath298 .
the bulk viscosity for zero magnetic field @xmath213 as function of oscillation frequency @xmath209 for various temperature is shown in fig .  [ zeta0 ] .
the  shoulder " structure and the temperature dependence of @xmath213 are easily understood from eq .
( [ vis1 ] ) and have been widely discussed in the literature  @xcite .     at zero magnetic field as function of the oscillation frequency @xmath209 for @xmath299 mev at @xmath300 ( solid , blue online ) @xmath301 ( dashed , red online ) and @xmath302 ( black , dotted online ) mev .
, width=283 ]    \2 ) when the magnetic field is extremely large , say , @xmath303 g , for our choice of parameters ( [ paras ] ) all the quarks are confined in their lowest landau level .
in this case we obtain @xmath304 and therefore @xmath305 we used the parameters ( [ paras ] ) and assumed physically interesting frequencies @xmath306 s@xmath163 .
the bulk viscosity @xmath0 vanishes as a consequence of vanishing @xmath193 when @xmath230 . since @xmath1 is now inversely proportional to @xmath307 it approaches zero for large @xmath2 .
therefore , both @xmath0 and @xmath1 are suppressed for large @xmath2 . in contrary , @xmath192 is enhanced by the extremely large magnetic field , see fig .
[ pressure_u ] .
\3 ) when the magnetic field is strong , but not strong enough to confine all the quarks to their lowest landau level , the situation becomes complicated . for our chosen parameters ( [ paras ] ) , this situation roughly corresponding to the interval @xmath308 g. in this case , a finite number of landau levels is occupied , and the essential observation is that @xmath200 and @xmath203 can be negative .
the behaviors of @xmath200 and @xmath203 are shown in fig .
[ c_u ] as functions of @xmath2 .
let us concentrate on the few levels just above the value @xmath309 g. when @xmath2 grows passing over @xmath310 or @xmath311 for each @xmath225 , both @xmath200 and @xmath203 change their sign .
more importantly , they have always opposite signs .
therefore , in this region , @xmath0 is negative which leads to hydrodynamic instability ( see the analysis in sec.[43 ] ) .
( dashed , blue online ) and @xmath203 ( dotted , red online ) as functions of @xmath2 at zero temperature.,width=283 ]    the numerical values of the bulk viscosities @xmath1 and @xmath0 are shown in fig .  [ zeta_u ] as functions of @xmath2 .
the parameters are those given in eq .
( [ paras ] ) .
we also fix the temperature @xmath312 mev and oscillation frequency @xmath313 .
both @xmath1 and @xmath0 have  quasi - periodic " oscillatory dependence on the magnetic field .
the two boundaries of each  period " correspond to a pair of neighboring @xmath314 and @xmath315 , and hence the period is roughly @xmath316 for large @xmath2 .
therefore , on average , the period increases as @xmath2 grows .
the amplitude of these oscillations also grows with increasing magnetic field until @xmath317 .
thereafter all the quarks are confined to their lowest landau levels and @xmath0 vanishes . from fig .
[ zeta_u ] we see that the magnitudes of @xmath0 and @xmath1 can be 100 to 200 times larger than their zero field value @xmath213 .
due to the unequal masses and charges of @xmath232 and @xmath11 quarks , @xmath1 and @xmath0 behave very irregularly .
we illustrate the zoomed - in curves around @xmath318 and @xmath319 g in the sub - panels , which look more regular .
the quasi - periodic structures are more evident in these sub - panels .     and @xmath0 scaled by the isotropic bulk viscosity @xmath213 as functions of the magnetic field @xmath2 at fixed frequency @xmath320 and temperature @xmath312 mev .
the dashed red curves correspond to viscosities lying in the unstable regions and would be not physically reachable .
sub - panels show the amplifications around @xmath321 and @xmath322 g. our parameters are given in eq .
( [ paras]).,title="fig:",width=340 ]   and @xmath0 scaled by the isotropic bulk viscosity @xmath213 as functions of the magnetic field @xmath2 at fixed frequency @xmath320 and temperature @xmath312 mev .
the dashed red curves correspond to viscosities lying in the unstable regions and would be not physically reachable .
sub - panels show the amplifications around @xmath321 and @xmath322 g. our parameters are given in eq .
( [ paras]).,title="fig:",width=340 ]    the most unusual feature seen in fig .
[ zeta_u ] is that for a wide range of field values , the transverse bulk viscosity @xmath0 is negative .
therefore , strange quark matter in this region is hydrodynamically unstable . besides this hydrodynamical instability , near each @xmath255
, there is a narrow window in which thermodynamical instability arises .
we depict @xmath1 and @xmath0 in these unstable regions by dashed red curves .
the solid blue curves correspond to the stable regime .
the magnetic field in a compact star need not be homogeneous and may have a complicated structure with poloidal and toroidal components .
furthermore , the fields will be functions of position in the star because of the density dependence of the parameters of the theory .
furthermore , the instabilities , described above , may lead to fragmentation of matter and formation of domain structures , where the regions with magnetic fields are separated from those without magnetic field by domain walls .
accordingly , only the averaged viscosities over some range of magnetic field have practical sense for assessing the large - scale behavior of matter . averaging over many oscillation periods in the stable region
, we find that the averaged values of @xmath1 and @xmath0 are much more regular , with their magnitudes restricted from 0 to several @xmath213 , see fig .
[ zetaaver ] . in obtaining the curves in fig .
[ zetaaver ] , we have eliminated the viscosities lying in the unstable regime .
the solid black curves are obtained by averaging over a short period @xmath323 .
the period was chosen such that the most rapid fluctuations are smeared out , but the oscillating structures over larger scale are intact .
the short - dashed red curves correspond to averaging over an even longer period , @xmath324 .
the result of long - period averaging is that @xmath1 first increases slowly and then drops down quickly once @xmath325 g ; similarly , @xmath0 first slowly decreases and then drops down very fast for @xmath326 g. such a dropping behavior reflects the fact that a large number of quarks are beginning to occupy the lowest landau level .     and
@xmath0 scaled by the isotropic bulk viscosity @xmath213 as functions of the magnetic field @xmath2 at fixed frequency @xmath327 and temperature @xmath312 mev .
the black - solid curve corresponds to averaging over a short period ( @xmath323 ) ; while the red short - dashed curve corresponds to averaging over a long period ( @xmath324 ) .
the viscosities lying in the unstable regions have been eliminated in the averaging.,title="fig:",width=283 ]   and @xmath0 scaled by the isotropic bulk viscosity @xmath213 as functions of the magnetic field @xmath2 at fixed frequency @xmath327 and temperature @xmath312 mev .
the black - solid curve corresponds to averaging over a short period ( @xmath323 ) ; while the red short - dashed curve corresponds to averaging over a long period ( @xmath324 ) .
the viscosities lying in the unstable regions have been eliminated in the averaging.,title="fig:",width=283 ]    we note that the appearances of thermodynamic , mechanical , and hydrodynamical instabilities are all induced by the landau quantization of the quark levels , _
i.e. _ , are quantum mechanical in nature .
more precisely , they are all due to the interplay between the landau levels and the fermi momentum ( reflected in the quantity @xmath255 ) .
additionally , the hydrodynamical instability requires that the quark matter is para - magnetized .
although we did our analysis by using the free quark gas approximation , eq .
( [ pstrong ] ) , it should be valid as long as there are sharp fermi surfaces ( low temperature ) , quantized landau levels ( high magnetic field ) and para - magnetization
. the appearances of these instabilities are expected to be robust feature for such systems .
we also checked that if one imposes the neutrality condition , _
i.e. _ , the condition @xmath328 , there is only minor quantitative change , while the qualitative conclusions are almost unchanged .
our choice of chemical potential @xmath329 mev roughly corresponds to the choice of @xmath330 for neutral strange quark matter .
the purpose of this section is to discuss the damping of the r - modes of newtonian models of strange stars by dissipation driven by the bulk viscosities @xmath0 and @xmath1 . as is well known , rotating equilibrium configurations of self - gravitating fluids are susceptible to instabilities at high rotation rates .
starting from the mass - shedding limit and going down with the rotation rate , the first instability point corresponds to the dynamical instability of the @xmath331 and @xmath332 mode .
this bar - mode instability is independent of the dissipative processes inside the star and occurs at values of the kinetic to potential energy ratio @xmath333  @xcite . for smaller rotation rates two secular instabilities with @xmath331 arise , each corresponding to a sign of @xmath334 . for incompressible fluids at constant density
the @xmath335 values for the onset of secular instabilities coincide .
one instability is driven by the viscosity , the other instability is driven by the gravitational radiation .
for realistic stars the @xmath335 values for the onset of these instabilities do not coincide ; relativity and other factors shift the viscosity - driven instability to higher values of @xmath335 . at the same time
the gravitational radiation instability is shifted to lower values of @xmath335 .
the gravitational radiation instability arises for the modes which are retrograde in the co - rotating frame , while prograde in the ( distant ) laboratory frame .
the underlying mechanism is the well established chandrasekhar - friedman - schutz ( cfs ) mechanism  @xcite .
the bulk and shear viscosities can prevent the development of the cfs instabilities , except in a certain window in the rotation and temperature plane .    in the following we shall concentrate on axial modes of newtonian stars , the so - called r - modes , which are known to undergo a cfs - type instability .
our main goal will be to assess the role of strong magnetic fields and bulk viscosity on the stability of these objects .
we shall adopt the formalism of refs .
@xcite for our study of the damping of the r - modes by bulk viscosity . for the sake of simplicity
we shall describe both fluid mechanics and gravity in the newtonian approximation .
the equations that describe the dynamical evolution of the star are    @xmath336    @xmath337    @xmath338    where @xmath339 is defined by the integral @xmath340 the quantity @xmath341 is the mass density of the fluid which is assumed to satisfy a barotropic equation of state , @xmath342 .
@xmath343 is the gravitational potential and @xmath272 is the gravitational constant . the potential @xmath344 is used to determine the velocity field @xmath84 .
the oscillation modes of a uniformly rotating star can be completely described in terms of two perturbation potentials @xmath345 and @xmath346 , where @xmath347 and @xmath348 are the potentials that correspond to the equilibrium configuration of the star .
we assume that the time and azimuthal angular dependence of any perturbed quantity is described by @xmath349 , where @xmath350 is an integer and @xmath351 is the frequency of the mode in laboratory frame .
let @xmath47 denote the rotation frequency of the star and @xmath209 denote the frequency of the perturbed quantity measured in the co - rotating frame [ which corresponds the @xmath209 in eqs .
( [ vis1 ] ) and ( [ vis2 ] ) because we will work in the co - rotating frame ] . for small @xmath47
there is a simple relation between @xmath352 , and @xmath209  @xcite , @xmath353 by linearizing the euler equation ( [ newtonian2 ] ) around the equilibrium configuration , the velocity perturbation @xmath354 is determined by  @xcite @xmath355 the tensor @xmath356 is a function of @xmath351 and the rotation frequency @xmath47 of the star , @xmath357},\;\;\;\end{aligned}\ ] ] where @xmath358 is a unit vector pointing along the rotation axis of the equilibrium star , which we assume to be parallel to the magnetic field , _ i.e. _ , @xmath359 in cartesian coordinate system . here
@xmath360 is the fluid velocity of the equilibrium star .    having the linearized euler equation
, one proceeds to the linearization of the mass continuity equation ( [ newtonian1 ] ) and the equation for the gravitational potential ( [ newtonian3 ] ) ; one finds @xmath361 these equations , together with the appropriate boundary conditions at the surface of the star for @xmath362 and at infinity for @xmath363 , determine the potentials @xmath362 and @xmath363 .    for slowly rotating stars , eq .
( [ newtonian4 ] ) can be solved order by order in @xmath47 , @xmath364 } , { \nonumber\\}{\delta}{\phi}=r^2{\omega}^2{\left[}{\delta}{\phi}_0+{\delta}{\phi}_2\frac{{\omega}^2}{{\pi}g { \rho}}+o({\omega}^4){\right]},\end{aligned}\ ] ] where @xmath365 is the radius of nonrotating star .
since we need only the perturbed velocity , we will focus on @xmath362 in the following discussion .
the zeroth - order contribution to the r - mode is generated by the following form of the potential @xmath366 , @xmath367 @xmath368 where @xmath369 is an arbitrary dimensionless constant and @xmath370 are the associated legendre polynomials .
it has been shown that the most unstable mode is the one with @xmath332  @xcite , therefore we shall consider only this case in the following discussion . substituting @xmath366 into eq .
( [ dv ] ) one obtains the first - order perturbed velocity , @xmath371 where @xmath372 and @xmath373 is the magnetic - type spherical harmonic function , @xmath374    it is straightforward to check that the first - order perturbed velocity satisfies @xmath375 therefore it does not contribute to the dissipation due to the bulk viscosities @xmath0 and @xmath1 . in order to see how the bulk viscosities damp the r - mode instability one must consider next - to - first order , _
i.e. _ , the third - order perturbed velocity which is generated by the potential @xmath376 .
one can not determine analytically @xmath376 from eqs .
( [ newtonian4])-([dudf ] ) , but the angular structure of @xmath376 can be well represented by the following spherical harmonics expansion  @xcite , @xmath377 the functions @xmath378 and @xmath379 have been determined numerically in ref .
. a useful approximation is provided by the following simple expressions    @xmath380    which excellently fit the numerical result .
we will use eqs .
( [ f1 ] ) and ( [ f2 ] ) in the following numerical calculation .
the energy of r - modes comes both from the velocity perturbation and the perturbation of the gravitational potential . for slowly rotating stars
, the main contribution comes from the velocity perturbation  @xcite .
then , the energy of the r - mode measured in the co - rotating frame is @xmath381 assuming spherical symmetry , we have @xmath382 this energy will be dissipated both by gravitational radiation and by the thermodynamic transport in the fluid  @xcite , @xmath383 the dissipation rate due to gravitational radiation is given by  @xcite @xmath384},{\nonumber\\}\end{aligned}\ ] ] where @xmath385 ^ 2}.\end{aligned}\ ] ] @xmath386 and @xmath387 are the mass and current multipole moments of the perturbation , @xmath388 taking into account eq .
( [ omega ] ) one obtains @xmath389 which implies that the total sign of @xmath390 is positive : gravitational radiation always increases the energy of the r - modes .    in order to compare the relative strengths of different dissipative processes ,
it is convenient to introduce the dissipative timescales defined by @xmath391 where the index @xmath392 labels the dissipative process .
the lowest - order contribution to @xmath393 comes from the current multipole moment @xmath394 . for the most important case @xmath395
, this leads to the following timescale ( derived for a simple polytropic equation of state @xmath396 )  @xcite @xmath397    the bulk viscosities @xmath0 and @xmath1 dissipate the energy of the r - mode according to @xmath398 accordingly , the time scales @xmath399 and @xmath400 are given by @xmath401    currently , the shear viscosities @xmath402 of strange quark matter are not known . in order to determine the damping of the r - mode by shear viscosity
, we take as a crude estimate the value of @xmath403 in the absence of a magnetic field  @xcite @xmath404 where @xmath405 is the coupling constant of strong interaction
. we will choose the value @xmath406 and apply eq .
( [ eta0 ] ) to highly degenerate 3-flavor quark matter with equal chemical potentials of all flavors , ( @xmath407 ) .
the contribution to the energy dissipation rate @xmath408 due to shear viscosity @xmath149 now becomes @xmath409 assuming a uniform mass density star , the time scale @xmath410 can be simply expressed as  @xcite @xmath411    the total time scale @xmath412 is given by the following sum @xmath413 which characterizes how fast the r - mode decays .
most importantly , if the sign of @xmath414 is negative the amplitude of the r - mode will not decay , rather it will increase with time .
thus , it is important to determine the critical angular velocity @xmath415 for the onset of instability @xmath416 at a given temperature , stars with @xmath417 will be unstable due to gravitational radiation .     of a strange quark star as a function of magnetic field @xmath2 at temperature @xmath418 mev .
the red dashed lines correspond to the unstable , while the blue solid lines correspond to the stable regime.,width=283 ]    figure.[rb ] shows the critical angular velocity @xmath415 of a strange quark star with mass @xmath419 and radius @xmath420 km as a function of the magnetic field @xmath2 near @xmath322 g. the temperature is fixed as @xmath418 mev and other parameters are taken according to eq .
( [ paras ] ) . in obtaining fig .
[ rb ] , we have taken into account the thermodynamical and hydrodynamical stability conditions .
the blue - solid curves correspond to the thermodynamically and hydrodynamically stable region , while the red - dashed curves correspond to unstable regions .
the critical angular velocity is strongly oscillating with increasing @xmath2 .
this behavior is due to the oscillating nature of the bulk viscosities @xmath1 and @xmath0 as shown in fig .  [ zeta_u ]
. thus , this macroscopic behavior originates from a purely quantum mechanical effect , namely the landau quantization of the energy levels of quarks . as discussed for the bulk viscosities @xmath1 and @xmath0 shown in fig .
[ zetaaver ] , averaging is needed to obtain physically relevant quantities . in fig .
[ rbaver ] we show the averaged critical angular velocity at various temperatures.the solid black curves are obtained by averaging over a short period @xmath421 , whereas the short - dashed red curves correspond to averaging over a long period @xmath422 .
it is seen that after short - period averaging , the critical angular velocity ( solid black curves ) shows regular oscillation , the amplitude of which is growing as the @xmath2-field increases .
the critical angular velocity @xmath415 displays a sharp drop for fields @xmath423 g ( short - dashed red curves ) , which is the consequence of the sharp drop of @xmath1 and @xmath0 shown in fig .  [ zetaaver ] .
thus , we conclude that for extremely large magnetic fields , the critical angular velocity at which the r - mode instability sets in could be significantly lower than in the absence of magnetic field .     of a strange quark star as a function of magnetic field @xmath2 at various temperatures @xmath424 , and @xmath425 mev .
the black solid curve corresponds to averaging over a short period ( @xmath421 ) ; while the red short - dashed curve corresponds to averaging over a long period ( @xmath422).,title="fig:",width=245 ]   of a strange quark star as a function of magnetic field @xmath2 at various temperatures @xmath424 , and @xmath425 mev .
the black solid curve corresponds to averaging over a short period ( @xmath421 ) ; while the red short - dashed curve corresponds to averaging over a long period ( @xmath422).,title="fig:",width=245 ]   of a strange quark star as a function of magnetic field @xmath2 at various temperatures @xmath424 , and @xmath425 mev .
the black solid curve corresponds to averaging over a short period ( @xmath421 ) ; while the red short - dashed curve corresponds to averaging over a long period ( @xmath422).,title="fig:",width=245 ]    figure.[rmode ] shows the window of the r - mode instability in the @xmath426 plane for a strange quark star of mass @xmath427 and radius @xmath420 km .
the regions above the respective curves correspond to the parameter space where the r - mode oscillations are unstable , _
i.e. _ , a star in this region will rapidly spin down by emission of gravitational waves .
the dashed green curve corresponds to vanishing bulk viscosities @xmath428 .
the solid black curve represents the ( in)stability window of an un - magnetized strange quark star .
the curves with symbols show the typical instability window for magnetic fields around @xmath321 g ( the red curve marked by triangles ) and @xmath429 g ( the blue curve , marked by circles ) .
the symbolled curves are obtained by using the bulk viscosities averaged over the period @xmath324 . for low temperatures , @xmath430 kev ,
the r - mode instability is suppressed mainly by the shear viscosity ; at these low temperatures the bulk viscosities are an insignificant source of damping , independent of how large the magnetic field is .
however , for larger temperatures the bulk viscosities dominate the damping of r - mode oscillations . for magnetic fields below @xmath431 g ,
the critical rotation frequency is almost independent of the @xmath2-field .
the r - mode instability window increases as the magnetic field grows . for fields @xmath5 g
it is a very sensitive function of the field , as a consequence of the rapid variation of the bulk viscosities with the field .
asymptotically , the instability window can become significantly larger than the window at zero magnetic field ( see also fig .
[ rbaver ] ) . for completeness ,
[ rmode ] also shows the observed distribution of low mass x - ray binaries ( lmxbs ) by the shadowed box , which corresponds to the typical temperatures ( @xmath432 k ) and rotation frequencies ( 300 - 700 hz ) of the majority of observed lmxbs @xcite .
it is seen that even in the case of extremely large magnetic fields , our instability window is consistent with the current lmxb data .    .
the black solid curve represents the window of unmagnetized strange quark matter .
the curves with symbols show the typical behavior of the instability window when the magnetic fields are around @xmath321 g ( red curve ) and @xmath429 g ( blue curve ) .
the shadowed box represents typical temperatures ( @xmath433 k ) and rotation frequencies ( 300 - 700 hz ) of the majority of observed lmxbs @xcite.,width=336 ]
in this paper we have studied anisotropic hydrodynamics of strongly magnetized matter in compact stars .
we find that there are in general eight viscosity coefficients , six of them are identified as shear viscosities , the other two , @xmath0 and @xmath1 , are bulk viscosities [ see eq .
( [ pmn ] ) ] .
we applied our formalism to magnetized strange quark matter and gave explicit expressions for the bulk viscosities [ eq .  ( [ vis1 ] ) and eq .  ( [ vis2 ] ) ] due to the non - leptonic weak reactions ( [ weak1 ] ) and ( [ weak2 ] ) . due to the landau quantization of the energy levels of charged particles in a strong magnetic field ,
the magnetic field dependence of @xmath1 and @xmath0 is very complicated and exhibits  quasi - periodic " de haas - van alphen type oscillations ( see fig .
[ zeta_u ] ) . for a magnetic field
@xmath434  g the effect of the magnetic field on the transport coefficients is small and the bulk viscosities can be well approximated by their zero - field values . for large fields
@xmath435  g the viscosities are substantially modified , @xmath0 may even become negative for some values of the @xmath2-field .
we showed that negative @xmath0 render the fluid hydrodynamically unstable .    for a number of reasons ( density dependence of parameters along the star profile , formation of domains , intrinsic multicomponent nature of the magnetic field ) the dependence of the transport coefficients on the magnetic fields
are needed at different resolutions , _
i.e. _ , they require some suitable averaging over a range of magnetic field . we have provided such averages over an increasingly larger scale .
we find that if the averaging period is small the bulk viscosities show regular oscillations , the amplitudes of which increase with magnetic field ( see the black solid curves in fig .
[ zetaaver ] ) .
these oscillations are smoothed out if we further increase the averaging scale . at this larger scale
the most interesting feature is the rapid drop in the bulk viscosity of the matter due to the confinement of quarks to the lowest landau level ; this occurs for magnetic fields in excess of @xmath436 g ( see the short - dashed red curves in fig .  [ zetaaver ] ) .    as an application
, we utilized our computed anisotropic bulk viscosities to study the problem of damping of r - mode oscillations in rotating newtonian stars .
we find that the instability window increases as the magnetic field is increased above the value @xmath437 g. by increasing the field one covers the entire range of parameter space which lies between the two extremes : the case when bulk viscosity vanishes ( dashed green curve in fig .
[ rmode ] , which corresponds to extremely large magnetic fields @xmath438 g for which the bulk viscosity drops to zero ) , and the case when the magnetic field is absent ( solid black curve in fig .  [ rmode ] ) .
the found novel dependence of the r - mode instability window on the magnetic field may help to distinguish quark stars from ordinary neutron stars with strong magnetic fields , since the latter are much more difficult to magnetize .
it would be interesting to see whether the objects that lie in between these extremes , _ e.g. _ , hybrid configurations featuring quark cores and hadronic envelopes  ( see ref .
@xcite and references therein ) , may interpolate smoothly between the physics of ordinary and strange compact objects .
we thank t.  brauner , t.  koide , b.  sad , a.  schmitt , and i.  shovkovy for helpful discussions .
this work is supported , in part , by the helmholtz alliance program of the helmholtz association , contract ha216/emmi `` extremes of density and temperature : cosmic matter in the laboratory '' and the helmholtz international center for fair within the framework of the loewe ( landesoffensive zur entwicklung wissenschaftlich - konomischer exzellenz ) program launched by the state of hesse .
m.  h. is supported by cas program
outstanding young scientists abroad brought - in " , cas key project kjcx3-syw - n2 , nsfc10735040 , nsfc10875134 , and by k.  c.  wong education foundation , hong kong .
the purpose of this appendix is to show explicitly that the eight different decompositions in eq .
( [ combination ] ) are linearly independent . to this end , let us write down a general linear combination of the eight different decompositions , @xmath439 if @xmath440 are linearly independent , the coefficients @xmath441 should all vanish for any values of the vectors @xmath25 and @xmath95 .
firstly , it is easy to see that the components @xmath442 and @xmath443 are independent of the other components , because they have odd parity under reflection @xmath444 , whereas the other six components have even parity under this transformation . besides that it is obvious that @xmath445 and @xmath443 are independent of each other .
therefore , we only need to treat the linear equation @xmath446 by contracting the indices @xmath447 and @xmath448 , we obtain the following three conditions @xmath449 contracting eq .
( [ linear0 ] ) with @xmath95 and @xmath450 we find the following two additional conditions @xmath451 contracting the indices @xmath448 and @xmath369 in eq .
( [ linear ] ) , we obtain one further condition @xmath452 the non - trivial solution of the set of eqs .
( [ linear2])-([linear5 ] ) is @xmath453 then , we have the following condition , @xmath454 the only possible solution is @xmath455 , which thus proves the independence of @xmath440 .
20 d.  d.  ivanenko and d.  f.  kurdgelaidze , astrophysics * 1 * , 251 ( 1965 ) [ astrofiz .   * 1 * , 479 ( 1965 ) ] .
n.  itoh , prog .
phys .   *
44 * , 291 ( 1970 ) .
j.  c.  collins and m.  j.  perry , phys .
lett .   * 34 * , 1353 ( 1975 ) .
e.  witten , phys .
d * 30 * , 272 ( 1984 ) .
f.  weber , prog .
phys .   * 54 * , 193 ( 2005 ) [ arxiv : astro - ph/0407155 ] .
n.  iwamoto , phys .
* 44 * , 1637 ( 1980 ) .
d.  blaschke , t.  klahn and d.  n.  voskresensky , astrophys .
j.   * 533 * , 406 ( 2000 ) [ arxiv : astro - ph/9908334 ] .
d.  blaschke , h.  grigorian and d.  n.  voskresensky , astron .
astrophys .
* 368 * , 561 ( 2001 ) [ arxiv : astro - ph/0009120 ] .
t.  schafer and k.  schwenzer , phys .
d * 70 * , 114037 ( 2004 ) [ arxiv : astro - ph/0410395 ] .
m.  alford , p.  jotwani , c.  kouvaris , j.  kundu and k.  rajagopal , phys .
d * 71 * , 114011 ( 2005 ) [ arxiv : astro - ph/0411560 ] .
p.  jaikumar , c.  d.  roberts and a.  sedrakian , phys .
c * 73 * , 042801 ( 2006 ) [ arxiv : nucl - th/0509093 ] .
a.  schmitt , i.  a.  shovkovy and q.  wang , phys .
d * 73 * , 034012 ( 2006 ) [ arxiv : hep - ph/0510347 ] .
r.  anglani , g.  nardulli , m.  ruggieri and m.  mannarelli , phys .
d * 74 * , 074005 ( 2006 ) [ arxiv : hep - ph/0607341 ] .
x.  g.  huang , q.  wang and p.  f.  zhuang , phys .
d * 76 * , 094008 ( 2007 ) [ arxiv:0706.4381 [ hep - ph ] ] .
x.  g.  huang , q.  wang and p.  f.  zhuang , int .
j.  mod .
e * 17 * , 1906 ( 2008 ) [ arxiv:0710.3634 [ hep - ph ] ] . c.  kouvaris , phys .
d * 79 * , 123008 ( 2009 ) [ arxiv:0901.1028 [ astro-ph.sr ] ] .
d.  g.  yakovlev and c.  j.  pethick , ann .
astrophys .   * 42 * , 169 ( 2004 ) [ arxiv : astro - ph/0402143 ] .
a.  sedrakian , prog .
* 58 * , 168 ( 2007 ) [ arxiv : nucl - th/0601086 ] .
d.  blaschke and h.  grigorian , prog .
phys .   * 59 * , 139 ( 2007 ) [ arxiv : astro - ph/0612092 ] .
n.  andersson , astrophys .
j.   * 502 * , 708 ( 1998 ) [ arxiv : gr - qc/9706075 ] .
j.  l.  friedman and s.  m.  morsink , astrophys .
j.   * 502 * , 714 ( 1998 ) [ arxiv : gr - qc/9706073 ] .
n.  andersson and k.  d.  kokkotas , int .
j.  mod .
d * 10 * , 381 ( 2001 ) [ arxiv : gr - qc/0010102 ] .
l.  lindblom , arxiv : astro - ph/0101136 .
p.  jaikumar , g.  rupak and a.  w.  steiner , phys .
d * 78 * , 123007 ( 2008 ) [ arxiv:0806.1005 [ nucl - th ] ] .
b.  knippel and a.  sedrakian , phys .
d * 79 * , 083007 ( 2009 ) [ arxiv:0901.4637 [ astro-ph.sr ] ] .
q.  d.  wang and t.  lu , phys .
b * 148 * , 211 ( 1984 ) .
r.  f.  sawyer , phys .
b * 233 * , 412 ( 1989 ) [ erratum - ibid .
b * 237 * , 605 ( 1990 ) ] .
j.  madsen , phys .
d * 46 * , 3290 ( 1992 ) .
x.  p.  zheng , s.  h.  yang and j.  r.  li , phys .
b * 548 * , 29 ( 2002 ) [ arxiv : hep - ph/0206187 ] .
z.  xiaoping , k.  miao , l.  xuewen and y.  shuhua , phys .
c * 72 * , 025809 ( 2005 ) [ arxiv : hep - ph/0507216 ] .
j.  d.  anand , n.  chandrika devi , v.  k.  gupta and s.  singh , pramana * 54 * , 737 ( 2000 ) .
b.  a.  sad , i.  a.  shovkovy and d.  h.  rischke , phys .
d * 75 * , 065016 ( 2007 ) [ arxiv : astro - ph/0607643 ] . m.  g.  alford and a.  schmitt , j.  phys .
g * 34 * , 67 ( 2007 ) [ arxiv : nucl - th/0608019 ] .
b.  a.  sad , i.  a.  shovkovy and d.  h.  rischke , phys .
d * 75 * , 125004 ( 2007 ) [ arxiv : astro - ph/0703016 ] . h.  dong , n.  su and q.  wang , phys .
d * 75 * , 074016 ( 2007 ) [ arxiv : astro - ph/0702104 ] . m.  g.  alford , m.  braby , s.  reddy and t.  schafer , phys .  rev .
c * 75 * , 055209 ( 2007 ) [ arxiv : nucl - th/0701067 ] . c.  manuel and f.  j.  llanes - estrada , jcap * 0708 * , 001 ( 2007 ) [ arxiv:0705.3909 [ hep - ph ] ] .
m.  mannarelli , c.  manuel and b.  a.  sad , phys .
* 101 * , 241101 ( 2008 ) [ arxiv:0807.3264 [ hep - ph ] ] .
m.  mannarelli and c.  manuel , arxiv:0909.4486 [ hep - ph ] .
h.  dong , n.  su and q.  wang , j.  phys .
g * 34 * , s643 ( 2007 ) [ arxiv : astro - ph/0702181 ] . m.  g.  alford , m.  braby and a.  schmitt , j.  phys .
g * 35 * , 115007 ( 2008 ) [ arxiv:0806.0285 [ nucl - th ] ] .
r.  c.  duncan and c.  thompson , astrophys .  j.   * 392 * , l9 ( 1992 ) . c.  thompson and r.  c.  duncan , astrophys .  j.   * 408 * , 194 ( 1993 ) . c.  thompson and r.  c.  duncan , mon .  not .
soc .   * 275 * , 255 ( 1995 ) . c.  thompson and r.  c.  duncan , astrophys .  j.   * 473 * , 322 ( 1996 ) . c.  kouveliotou
_ et al . _ ,
astrophys .  j.   * 510 * , l115 ( 1999 ) [ arxiv : astro - ph/9809140 ] . c.  y.  cardall , m.  prakash and j.  m.  lattimer , astrophys .
j.   * 554 * , 322 ( 2001 ) [ arxiv : astro - ph/0011148 ] .
a.  e.  broderick , m.  prakash and j.  m.  lattimer , phys .
b * 531 * , 167 ( 2002 ) [ arxiv : astro - ph/0111516 ] .
d.  lai and s.  l.  shapiro , astrophys .
j.   * 383 * , 745 ( 1991 ) .
s.  chakrabarty , d.  bandyopadhyay and s.  pal , phys .
* 78 * , 2898 ( 1997 ) [ arxiv : astro - ph/9703034 ] .
d.  bandyopadhyay , s.  chakrabarty and s.  pal , phys .
lett .   * 79 * , 2176 ( 1997 ) [ arxiv : astro - ph/9703066 ] .
a.  broderick , m.  prakash and j.  m.  lattimer , arxiv : astro - ph/0001537 .
s.  r.  de groot , _ the maxwell equations .
non - relativistic and relativistic derivations from electron theory _ , amsterdam , north - holland pub .
co. , 1969 .
n.  sadooghi , arxiv:0905.2097 [ hep - ph ] .
w.  a.  hiscock and l.  lindblom , phys .
d * 31 * , 725 ( 1985 ) .
w.  israel and j.  m.  stewart , annals phys .
* 118 * , 341 ( 1979 ) .
e.  m.  lifshitz and l.  p.  pitaevskii , _ physcial kinetics , course of theoretical physics , vol 10 _ , pergamon press ltd . , 1981 .
m.  chaichian , s.  s.  masood , c.  montonen , a.  perez martinez and h.  perez rojas , phys .
lett .   * 84 * , 5261 ( 2000 ) [ arxiv : hep - ph/9911218 ] .
d. kondepudi and i. prigogine , _ modern thermodynamics _ , john wiley and sons , 1998 .
s.  chakrabarty , d.  bandyopadhyay and s.  pal , int .
j.  mod .
a * 13 * , 295 ( 1998 ) .
z.  xiao - ping , z.  xia and l.  xue - wen , arxiv : astro - ph/0507390 .
j. r. ipser and l. lindblom , astrophys .
j. * 355 * , 226 ( 1990 ) .
l.  lindblom , g.  mendell and b.  j.  owen , phys .
d * 60 * , 064006 ( 1999 ) [ arxiv : gr - qc/9902052 ] .
l.  lindblom , b.  j.  owen and s.  m.  morsink , phys .
lett .   * 80 * , 4843 ( 1998 ) [ arxiv : gr - qc/9803053 ]
. b.  j.  owen , l.  lindblom , c.  cutler , b.  f.  schutz , a.  vecchio and n.  andersson , phys .
d * 58 * , 084020 ( 1998 ) [ arxiv : gr - qc/9804044 ] .
b.  a.  sad , arxiv:0806.3359 [ astro - ph ] .
e. f. brown , l. bildsten and p. chang , astrophys .
j. * 574 * , 920 ( 2002 ) .