Bulk viscosity of interacting magnetized strange quark matter

Bulk viscosity of interacting magnetized strange quark matter 1 corresp first-author jfxu@aynu.edu.cn jfxu@aynu.edu.cn jfxu@aynu.edu.cn School of Physics Electrical Engineering, AnYang Normal University, Anyang 455000, China School of Physics Electrical Engineering, AnYang Normal University, Anyang 455000, China The bulk viscosity of interacting strange quark matter in a strong external magnetic field B m with a real equation of state is investigated. It is found that interquark interactions can significantly increase the bulk viscosity, and the magnetic field B m can cause irregular oscillations in both components of the bulk viscosity, <math id="M1"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> (parallel to B m ) and <math id="M2"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> (perpendicular to B m ). A comparison with non-interacting strange quark matter reveals that when B m is sufficiently large, <math id="M3"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> is more affected by interactions than <math id="M4"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> . Additionally, the quasi-oscillation of the bulk viscosity with changes in density may facilitate the formation of magnetic domains. Moreover, the resulting r -mode instability windows are in good agreement with observational data for compact stars in low-mass X-ray binaries. Specifically, the r -mode instability window for interacting strange quark matter in high magnetic fields has a minimum rotation frequency exceeding 1050 Hz, which may explain the observed very high spin frequency of a pulsar with ν =1122 Hz. Strange quark matter Bulk viscosity Strong magnetic field Strange star R -mode instability window 1　Introduction 1 Since the first direct detection of gravitational waves (GWs)[ 1 1 ] emitted during the coalescence of a binary black hole (BH), dozens of GW events have been observed during the first and second observing runs of the advanced GW detector network[ 2 2 ]. In addition to the observations of binary BH mergers, the first detection of GWs from a binary neutron star (NS) inspiral (the GW170817 event)[ 3 3 ] is extraordinarily significant, as the observation of the GWs emitted in this process, possibly combined with electromagnetic observation of the same source[ 4 4 , 5 5 ], may yield insight into the structure of NSs and the equation of state (EOS) of matter under extreme conditions[ 6 6 - 10 10 6-10 ]. Although transient GWs originate from the coalescence of compact stellar objects, the principal sources of continuous gravitational emission are expected to be spinning NSs and/or quark stars (QSs), which need not be in binary systems. A comprehensive review of the mechanisms of continuous GW emission is given in Ref.[ 11 11 ]. Continuous GWs can typically be generated by various processes that produce asymmetry[ 12 12 ]. A pulsar with a mass quadrupole may emit GWs with a spin frequency equal to or twice that of the pulsar, whereas some NSs may radiate GWs strongly through a current quadrupole via r -modes, which oscillate at approximately four-thirds of the spin frequency. Unstable oscillation modes, in particular r -modes with a sufficiently large saturation amplitude, have attracted considerable attention as potential sources of detectable GWs. Methods of searching for GWs from the r -modes of known pulsars are described in Refs.[ 13 13 - 16 16 13-16 ]. The emission of GWs can generally drive r -mode oscillations of compact stars with a certain spin frequency and temperature via the Chandrasekhar–Friedman–Schutz mechanism[ 17 17 , 18 18 ]. In addition, when r -mode oscillation with a sufficiently large saturation amplitude reaches an unstable state, it can in turn cause strong GW emission, which could carry away the angular momentum of compact stars, resulting in a sharp decrease in the spin frequency. This behavior suggests that r -mode instability is likely to play an important role in the evolution of the post-merger remnant[ 19 19 ]. Moreover, the presence of r -mode instability results in theoretical difficulties in explaining the high spin frequencies of pulsars. To solve this problem, different scenarios have been proposed[ 20 20 - 29 29 20-29 ]. One possible effective method emerged from research in recent decades, which indicated that interactions between quarks can increase the bulk viscosity of strange quark matter (SQM) by 1-2 orders of magnitude[ 30 30 - 34 34 30-34 ]. The large bulk viscosity can reduce the r -mode instability window; consequently, theoretical calculations are consistent with astrophysical observations[ 22 22 ]. In Ref.[ 35 35 ], adopting a quark mass scaling with both linear confinement and perturbative interactions, we investigated the bulk viscosity of SQM in the equivparticle model. When we applied the resulting enhanced bulk viscosity, we found that the r -mode instability window for canonical strange stars with 1.4 M ⊙ is in good agreement with the observational frequencies and temperatures of pulsars in low-mass X-ray binaries (LMXBs). Moreover, it is well known that the EOS of NS matter is still unclear and can be affected by many physical parameters such as the symmetry energy[ 36 36 , 37 37 ] and the strong magnetic fields (on the order of approximately 10 12 -10 13 G)[ 38 38 - 43 43 38-43 ] that may be present on the surface of compact stars. For so-called magnetars, the magnetic field can even be as large as 10 14 -10 15 G[ 44 44 , 45 45 ]. In fact, the largest magnetic field that can be sustained by strange stars is estimated to be approximately 1.5 × 10 20 G[ 46 46 ]. According to a previous study, SQM will be more stable when the magnetic field is included in the EOS[ 42 42 ]. Additionally, a strong magnetic field can strongly suppress the reaction rate of the non-leptonic weak interaction <math id="M5"><mrow><mi>u</mi><mo>+</mo><mi>d</mi><mo>↔</mo><mi>u</mi><mo>+</mo><mi>s</mi></mrow></math> , which is one source of the large bulk viscosity of SQM and is expected to affect the viscosity of SQM[ 47 47 ]. In addition, the bulk viscosity of magnetized NS matter was studied in Ref.[ 48 48 ]. Given the important role of bulk viscosity in the emission of continuous GWs by compact stellar objects, in this study we investigated the bulk viscosity of SQM with both strong interactions and a magnetic field in the equivparticle model[ 42 42 , 43 43 , 49 49 ]. First, in Sect. 2, we illustrate the formulas for calculating the bulk viscosity of the interacting magnetized SQM. Next, in Sect. 3, we report and discuss the numerical results. Finally, a short summary is presented in Sect. 4. 2　Bulk viscosity of magnetized SQM in equivparticle model 1 In the equivparticle model, the quark masses m i , ( i = u,d, and s ) vary with the baryon number density n b , which effectively mimics the strong interactions between quarks. To study the effect of interactions on the bulk viscosity of magnetized SQM, we take the quark mass parameterized as follows[ 49 49 ]: <math display="block" id="M6"><mtable columnalign="left"><mtr><mtd><msub><mi>m</mi><mi>i</mi></msub><mo>≡</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>m</mi><mtext>I</mtext></msub><mo stretchy="false">(</mo><msub><mi>n</mi><mtext>b</mtext></msub><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mi>C</mi><msubsup><mi>n</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><mfrac><mi>D</mi><mrow><msubsup><mi>n</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> where m i 0 is the quark current mass, m I is the interacting part of the quark mass, and the model parameters C and D indicate the strength of perturbative interactions and confinement effects, respectively. Because electrons do not participate in strong interactions, their mass is m e = m e0 =0.511 MeV. It is quite convenient to treat the thermodynamics in the equivparticle model, as the bare chemical potentials of quarks are replaced by effective potentials, that is, <math id="M7"><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>→</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> . Therefore, the thermodynamic potential density for particle species i in a strong external magnetic field B m has the same form as that of the free particles: <math display="block" id="M8"><mtable columnalign="left"><mtr><mtd><msub><mi>Ω</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mo stretchy="false">(</mo><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mtd></mtr><mtr><mtd><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup><mi>ln</mi><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow><mrow><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></mfrac><mo stretchy="false">)</mo><mn>,</mn></mtd></mtr></mtable></math> where g i =6 is the degeneracy factor for quarks (the value is 2 for electrons), e i is the absolute value of the electric charge of the particles, <math id="M9"><mrow><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo>=</mo><msqrt><mrow><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><mi>ν</mi><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow></msqrt></mrow></math> , and the upper bound of the summation of ν max can be determined by ensuring that the argument of the square root is positive. To derive Eq. (2), we have conventionally assumed a constant magnetic field B m = B m, z along the z axis. Consequently, the effective thermodynamic potential density for magnetized SQM can be written as <math display="block" id="M10"><mrow><msub><mi>Ω</mi><mn>0</mn></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mtext>I</mtext></munder></mstyle><msub><mi>Ω</mi><mi>i</mi></msub><mn>,</mn><mtext> </mtext><mo stretchy="false">(</mo><mi>i</mi><mo>=</mo><mi>u</mi><mn>,</mn><mi>d</mi><mn>,</mn><mi>s</mi><mn>,</mn><mtext>and</mtext><mtext> </mtext><mi>e</mi><mo stretchy="false">)</mo><mn>.</mn></mrow></math> To investigate the bulk viscosity of SQM in a magnetic field, it is necessary to calculate the magnetization 𝒨 of SQM. The contribution to magnetization from particle species i can be obtained by combining Eq. (2) with the relation <math id="M11"><mrow><msub><mi mathvariant="script">M</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mo>∂</mo><msub><mi>Ω</mi><mi>i</mi></msub><mo>/</mo><mo>∂</mo><msub><mi>B</mi><mtext>m</mtext></msub></mrow></math> : <math display="block" id="M12"><mtable columnalign="left"><mtr><mtd><msub><mi mathvariant="script">M</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mo stretchy="false">(</mo><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">[</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt><mo>−</mo><mo stretchy="false">(</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><mo>+</mo><mn>4</mn><mi>ν</mi><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub><mo stretchy="false">)</mo><mi>ln</mi><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow><mrow><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></mfrac><mo stretchy="false">]</mo><mn>.</mn></mtd></mtr></mtable></math> Because of the strong magnetic field B m , the longitudinal pressure and transverse pressure become different. They are given by <math display="block" id="M13"><mrow><msub><mi>P</mi><mo>∥</mo></msub><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mi>μ</mi><mtext>I</mtext></msub><msub><mi>n</mi><mi>i</mi></msub></mrow></math> and <math display="block" id="M14"><mrow><msub><mi>P</mi><mo>⊥</mo></msub><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mi>μ</mi><mtext>I</mtext></msub><msub><mi>n</mi><mi>i</mi></msub><mo>−</mo><mi mathvariant="script">M</mi><msub><mi>B</mi><mtext>m</mtext></msub><mn>,</mn></mrow></math> respectively, where μ I results from the requirement of thermodynamic consistency, which is crucial when the quark masses are density-dependent. It can be expressed in the following form: <math display="block" id="M15"><mtable columnalign="left"><mtr><mtd><msub><mi>μ</mi><mtext>I</mtext></msub><mo>=</mo><mo>−</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mtext>I</mtext></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mtext>I</mtext></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mo>≡</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msub><mi>f</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>n</mi><mtext>b</mtext></msub><mn>,</mn><msub><mi>m</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>n</mi><mtext>b</mtext></msub><mo stretchy="false">)</mo><mn>,</mn><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mn>,</mn><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>m</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>n</mi><mtext>b</mtext></msub><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mn>.</mn></mtd></mtr></mtable></math> By definition, the effective chemical potential μ i * and bare chemical potential μ i can be written in terms of μ I as follows: <math display="block" id="M16"><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>−</mo><msub><mi>μ</mi><mtext>I</mtext></msub><mn>.</mn></mrow></math> The particle number density can be given by the conventional thermodynamic relation @ <math id="M17"><mrow><msub><mi>n</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>/</mo><mo>∂</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> , which reads <math display="block" id="M18"><mrow><msub><mi>n</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mo stretchy="false">(</mo><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub><mo stretchy="false">)</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt><mn>.</mn></mrow></math> Further information about the equivparticle model of magnetized SQM is given in Ref.[ 43 43 ]. Like the pressure, the bulk viscosity of SQM becomes anisotropic in a strong magnetic field. By using a local linear response method, Huang et al.[ 28 28 ] presented explicit expressions for <math id="M19"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> and <math id="M20"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> , which are transverse and parallel to the external strong magnetic field B m , respectively. Both components of the bulk viscosity originate from the non-leptonic weak interaction <math id="M21"><mrow><mi>u</mi><mo>+</mo><mi>d</mi><mo>↔</mo><mi>u</mi><mo>+</mo><mi>s</mi></mrow></math> and were studied using a simple bag model that did not include interactions and β -equilibrium between quarks[ 28 28 ]. However, as pointed out above, extensive investigations have shown that interactions between quarks also make an important contribution to the bulk viscosity of SQM[ 30 30 - 33 33 30-33 , 35 35 ]. The bulk viscosities <math id="M22"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M23"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> are given as[ 28 28 ] <math display="block" id="M24"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub><mo>=</mo><mfrac><mrow><mi>λ</mi><msubsup><mi>C</mi><mo>∥</mo><mn>2</mn></msubsup></mrow><mrow><msup><mi>ω</mi><mn>2</mn></msup><mo>+</mo><msup><mi>λ</mi><mn>2</mn></msup><msup><mi>W</mi><mn>2</mn></msup></mrow></mfrac><mtext>and</mtext><msub><mi>ζ</mi><mo>⊥</mo></msub><mo>=</mo><mfrac><mrow><mi>λ</mi><msub><mi>C</mi><mo>⊥</mo></msub><msub><mi>C</mi><mo>∥</mo></msub></mrow><mrow><msup><mi>ω</mi><mn>2</mn></msup><mo>+</mo><msup><mi>λ</mi><mn>2</mn></msup><msup><mi>W</mi><mn>2</mn></msup></mrow></mfrac><mn>,</mn></mrow></math> respectively, where <math display="block" id="M25"><mtable columnalign="left"><mtr><mtd><msub><mi>C</mi><mo>∥</mo></msub><mo>≃</mo><msub><mi>n</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>d</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>d</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>−</mo><msub><mi>n</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>s</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mtd></mtr><mtr><mtd><msub><mi>C</mi><mo>⊥</mo></msub><mo>≃</mo><msub><mi>C</mi><mo>∥</mo></msub><mo>−</mo><mi>X</mi><msub><mi>B</mi><mtext>m</mtext></msub><mn>,</mn><mi>W</mi><mo>=</mo><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>s</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>d</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>d</mi></msub></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> with <math display="block" id="M26"><mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="script">M</mi></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>d</mi></msub></mrow></mfrac><mo>−</mo><mfrac><mrow><mo>∂</mo><mi mathvariant="script">M</mi></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>s</mi></msub></mrow></mfrac><mn>.</mn></mrow></math> In addition, <math display="block" id="M27"><mtable columnalign="left"><mtr><mtd><mi>λ</mi><mo>=</mo><mfrac><mrow><mn>4</mn><msubsup><mi>B</mi><mtext>m</mtext><mn>4</mn></msubsup><msup><mover accent="true"><mi>G</mi><mo>˜</mo></mover><mn>2</mn></msup><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>e</mi><mi>u</mi><mn>2</mn></msubsup><msub><mi>e</mi><mi>d</mi></msub><msub><mi>e</mi><mi>s</mi></msub><msub><mi>μ</mi><mi>d</mi></msub></mrow><mrow><mn>5</mn><msup><mi>π</mi><mn>3</mn></msup></mrow></mfrac><msup><mrow><mo>(</mo><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>u</mi><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>u</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac></mrow><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mo>×</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>d</mi><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>d</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>s</mi><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>s</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> where <math id="M28"><mrow><msup><mover accent="true"><mi>G</mi><mo>˜</mo></mover><mn>2</mn></msup><mo>≡</mo><msubsup><mi>G</mi><mtext>F</mtext><mn>2</mn></msubsup><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>sin</mi></mrow></mstyle><mn>2</mn></msup><msub><mi>θ</mi><mtext>C</mtext></msub><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>cos</mi></mrow></mstyle><mn>2</mn></msup><msub><mi>θ</mi><mtext>C</mtext></msub><msup><mrow><mo>=</mo><mn>6.46</mn><mo>×</mo><mn>10</mn></mrow><mrow><mo>−</mo><mn>24</mn></mrow></msup><msup><mrow><mtext>MeV</mtext></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msup></mrow></math> is the Fermi constant. In addition, ω is the oscillation frequency of compact stars, which is typically on the order of the magnitude of the rotation frequency, <math id="M29"><mrow><mn>1</mn><mtext> </mtext><msup><mtext>s</mtext><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>≲</mo><mi>ω</mi><mo>≲</mo><msup><mrow><mn>10</mn></mrow><mn>3</mn></msup><mtext> </mtext><msup><mtext>s</mtext><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> . Here we should emphasize that when the bulk viscosity is calculated in the equivparticle model, the chemical potentials under the square root in Eq. (13) should be replaced with the effective potentials, that is, <math id="M30"><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mn>2</mn></msubsup><mo>−</mo><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></msqrt><mo>→</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></msqrt></mrow></math> , as a result of the integral in the momentum space. In a spherical coordinate frame, the integral is taken from zero to the Fermi momentum, which can be associated with the effective chemical potentials in the equivparticle model at zero temperature. For more information on the reaction rate of the weak process <math id="M31"><mrow><mi>u</mi><mo>+</mo><mi>d</mi><mo>↔</mo><mi>u</mi><mo>+</mo><mi>s</mi></mrow></math> , readers may refer to Ref.[ 50 50 , 51 51 ]. To calculate the bulk viscosities in the equivparticle model, we must find the derivatives <math id="M32"><mrow><mo>∂</mo><msub><mi>μ</mi><mi>s</mi></msub><mo>/</mo><mo>∂</mo><msub><mi>n</mi><mi>s</mi></msub></mrow></math> and <math id="M33"><mrow><mo>∂</mo><msub><mi>μ</mi><mi>d</mi></msub><mo>/</mo><mo>∂</mo><msub><mi>n</mi><mi>d</mi></msub></mrow></math> . Starting from Eq. (9) and using Eqs. (7) and (8), one can obtain the derivative relations between the bare chemical potential and particle number density after considerable algebra, which gives <math display="block" id="M34"><mrow><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>δ</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><mo stretchy="false">(</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mi>A</mi><mi>B</mi><mo>+</mo><mfrac><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow><mn>3</mn></mfrac><mfrac><mrow><mtext>d</mtext><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mi>A</mi></mrow></mfrac><mn>,</mn></mrow></math> where <math display="block" id="M35"><mrow><mi>A</mi><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mn>,</mn></mrow></math> and <math display="block" id="M36"><mrow><mi>B</mi><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mrow><mo>(</mo><mrow><mfrac><mrow><mo>∂</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac><mfrac><mrow><mtext>d</mtext><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>∂</mo><msub><mi>f</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></mfrac><mfrac><mrow><mtext>d</mtext><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mtext>d</mtext><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mfrac><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><msub><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> As a special case, if the quark mass is independent of the baryon number density, it is not difficult to obtain <math display="block" id="M37"><mrow><mfrac><mrow><mo>∂</mo><msub><mi>μ</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><msub><mi>δ</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>e</mi><mi>i</mi></msub><msub><mi>B</mi><mtext>m</mtext></msub></mrow><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub></mrow></munderover></mstyle><mfrac><mrow><mn>2</mn><mo>−</mo><msub><mi>δ</mi><mrow><mi>ν</mi><mn>0</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac><msub><mi>μ</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> from Eq. (14), which is exactly the same as the results reported in Ref.[ 28 28 ]. 3　Numerical results and discussion 1 To calculate the bulk viscosities, the quark current masses are set to m u 0 =5 MeV, m d 0 =10 MeV, and m s 0 =100 MeV. In addition, the model parameters C and D are set to C =0.7 and <math id="M52"><mrow><msqrt><mi>D</mi></msqrt><mo>=</mo><mn>129</mn></mrow></math> MeV, which can ensure that the SQM is in the absolutely stable state[ 49 49 ]. The bulk viscosity can be understood in some sense as the energy dissipation rate, which naturally is closely related to the pressure. Therefore, in Fig. 1 Fig. 1 , we show the anisotropic pressures <math id="M53"><mrow><msub><mi>P</mi><mo>∥</mo></msub></mrow></math> and <math id="M54"><mrow><msub><mi>P</mi><mo>⊥</mo></msub></mrow></math> as functions of magnetic field B m with baryon number density n b =2 n 0 , where n 0 =0.17 fm -3 is the nuclear saturation density. Regardless of whether interactions are included, when <math id="M55"><mrow><msub><mi>B</mi><mtext>m</mtext></msub><mo>≳</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>18</mn></mrow></msup></mrow></math> G, <math id="M56"><mrow><msub><mi>P</mi><mo>∥</mo></msub></mrow></math> and <math id="M57"><mrow><msub><mi>P</mi><mo>⊥</mo></msub></mrow></math> become distinguishable. Then, with increasing B m , <math id="M58"><mrow><msub><mi>P</mi><mo>∥</mo></msub></mrow></math> first increases and then decreases. However, <math id="M59"><mrow><msub><mi>P</mi><mo>⊥</mo></msub></mrow></math> decreases continuously. The reason is that with increasing B m , increasing numbers of particles are confined to lower Landau levels. Moreover, when interactions between quarks are considered, both <math id="M60"><mrow><msub><mi>P</mi><mo>∥</mo></msub></mrow></math> and <math id="M61"><mrow><msub><mi>P</mi><mo>⊥</mo></msub></mrow></math> are reduced. Figure 2 Figure 2 shows the magnetization <math id="M62"><mi mathvariant="script">M</mi></math> in three cases. The dashed line shows the magnetization of the SQM in the bag model with quarks in β -equilibrium, whereas the dotted line represents magnetized SQM in the equivparticle model with both interactions and β -equilibrium. Like the pressure in Fig. 1 Fig. 1 , the magnetization can be decreased significantly by interactions between quarks. Furthermore, for comparison, we show the results presented in Ref.[ 28 28 ] (solid line), where the chemical potentials of u , d , and s quarks are μ u = μ d = μ s =400 MeV. By contrast, here they are set to <math id="M63"><mrow><msubsup><mi>μ</mi><mi>u</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>μ</mi><mi>d</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>μ</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><mn>300</mn></mrow></math> MeV. From Fig. 2 Fig. 2 , for a sufficiently large magnetic field (e.g., <math id="M64"><mrow><msub><mi>B</mi><mtext>m</mtext></msub><mo>≳</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>19.3</mn></mrow></msup></mrow></math> G), the magnetization remains unchanged without β -equilibrium, whereas it decreases sharply when β -equilibrium is considered. The reason is that as the magnetic field B m increases, an increasing number of particles transition from high Landau levels to low Landau levels, and thus ν max gradually becomes small until ν max =0 is reached. In this process, the oscillation in magnetization gradually becomes distinct and ultimately disappears, leaving the magnetization <math id="M65"><mi mathvariant="script">M</mi></math> as a function of only the chemical potentials, and thus a constant. Figure 3 Figure 3 shows the anisotropic bulk viscosities as functions of baryon number density n b under a constant magnetic field B m =10 18.5 G. The dashed lines correspond to the bulk viscosities without interactions, whereas the solid lines show the bulk viscosities with both perturbative interactions and quark confinement effects. Although the oscillation of the parallel bulk viscosity <math id="M66"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> changes significantly, the magnitude of <math id="M67"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> does not change greatly on average. However, the magnitude of the transverse bulk viscosity <math id="M68"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> increases significantly. Therefore, <math id="M69"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> is likely to be more susceptible to interactions than <math id="M70"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> when B m is sufficiently large. In fact, when interactions are not considered, the values of <math id="M71"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> are even negative, which implies that QSs made of SQM in this state are hydrodynamically unstable[ 28 28 ]. If the interactions are included, this situation can be improved greatly. Moreover, whenever a new Landau level appears, both <math id="M72"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M73"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> suddenly decrease. In addition, these sudden decreases in the bulk viscosities with changes in density may result in the fragmentation of matter and the formation of a magnetic domain in QSs from the deep interior to the surface. Another possible reason for the formation of magnetic domains is the hydrodynamic instability caused by the negative bulk viscosity. After magnetic domain structure is formed, regions with a magnetic field become separated from those without a magnetic field by domain walls. Consequently, in a sense, only the averaged bulk viscosity has practical meaning for the large-scale behavior of matter over some range of magnetic fields. Figure 4 Figure 4 shows <math id="M74"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M75"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> as functions of magnetic field. For fixed n b , at low magnetic field strength B m , <math id="M76"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M77"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> are clearly increased by interactions, in agreement with previous results where the effects of magnetic fields were not taken into account. As the magnetic field becomes stronger, the irregularity of the oscillation, which originates from the decrease in occupied Landau levels and the unequal masses and charges of different types of particles, becomes clear. The most severe problem that appears in Fig. 4 Fig. 4 is the negative values of <math id="M78"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> (dashed line in the lower panel), which can be greatly improved by including interactions between quarks, except at extremely strong magnetic fields. Furthermore, to study the properties of SQM in the stable state, the model parameters C and D should be constrained to the absolutely stable region of the stability window[ 49 49 ], where the approximate relationship between C and D can be roughly fitted as[ 35 35 ] <math display="block" id="M80"><mrow><msup><mi>D</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>/</mo><mtext>MeV</mtext><mo>=</mo><mo>−</mo><mfrac><mrow><mn>27</mn></mrow><mrow><mn>0.7</mn></mrow></mfrac><mi>C</mi><mo>+</mo><mn>156,</mn><mtext> </mtext><mi>C</mi><mo>∈</mo><mo stretchy="false">[</mo><mn>0,0.7</mn><mo stretchy="false">]</mo><mn>.</mn></mrow></math> According to this relationship, when C increases, D 1/2 decreases, which is shown on the upper X axis in Fig. 5 Fig. 5 . In addition, with increasing C and decreasing D , both <math id="M81"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M82"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> decrease simultaneously. However, according to previous studies, the bulk viscosity should increase with increasing interquark interactions, including perturbative interactions and/or quark confinements effects. Therefore, the results shown in Fig. 5 Fig. 5 imply that confinement effects may contribute more to the bulk viscosity than perturbative interactions for the parameters n b =2 n 0 and B m =10 18.5 G. Furthermore, the sudden decreases in both <math id="M83"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M84"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> still originate from the variation of the occupied Landau levels. Next, we discuss the calculations of the r -mode instability window of strange stars using the obtained bulk viscosities of magnetized SQM. To obtain the instability window, the following equation is generally solved: <math display="block" id="M85"><mrow><mfrac><mn>1</mn><mi>τ</mi></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><mtext>gw</mtext></mrow></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><mtext>sv</mtext></mrow></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><mtext>bv</mtext></mrow></msub></mrow></mfrac><mo>+</mo><mo>⋯</mo><mn>,</mn></mrow></math> where τ gw is the characteristic time scale of GW emission; τ sv and τ bv represent the damping time scales of the shear and bulk viscosity, respectively; and the ellipse denotes other dissipation mechanisms, such as surface rubbing[ 16 16 , 52 52 - 54 54 52-54 ]. Here, it should be stressed that the damping time scale of the bulk viscosity arises from both <math id="M86"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> and <math id="M87"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> because of the magnetic field; that is, <math display="block" id="M88"><mrow><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><mtext>bv</mtext></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>τ</mi><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></msub></mrow></mfrac><mn>.</mn></mrow></math> Eq. (19) is solved using the method employed in Ref.[ 35 35 , 55 55 , 56 56 ]. According to the literature[ 57 57 , 58 58 ], the driving time scale due to GW emission is <math display="block" id="M89"><mrow><msub><mi>τ</mi><mrow><mtext>gw</mtext></mrow></msub><mo>=</mo><mo>−</mo><mn>3.26</mn><mtext> s</mtext><msup><mrow><mrow><mo>(</mo><mrow><mi>π</mi><mi>G</mi><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mo>/</mo><msup><mi>Ω</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mn>3</mn></msup><mn>,</mn></mrow></math> where G =6.707 × 10 -45 MeV -2 is the gravitational constant, <math id="M90"><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover></math> is the mean density of a compact star, and Ω is the angular rotation frequency. The damping time scale of the shear viscosity[ 59 59 ] is <math display="block" id="M91"><mrow><msub><mi>τ</mi><mrow><mtext>sv</mtext></mrow></msub><mo>=</mo><mn>5.37</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mn>8</mn></msup><mtext> s</mtext><msup><mrow><mrow><mo>(</mo><mrow><msub><mi>α</mi><mi>s</mi></msub><mo>/</mo><mn>0.1</mn></mrow><mo>)</mo></mrow></mrow><mrow><mn>5</mn><mo>/</mo><mn>3</mn></mrow></msup><msubsup><mi>T</mi><mn>9</mn><mrow><mn>5</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mn>,</mn></mrow></math> where α s is the strong coupling; in the numerical calculations, it is α s =0.1. The time scale of the bulk viscosity can be written in two different forms depending on the temperature. The high- T limit applies for T > 10 9 K, and the low- T limit applies for T < 10 9 K. In the equivparticle model at the high- T limit, <math id="M92"><mrow><msubsup><mi>τ</mi><mrow><mtext>bv</mtext></mrow><mrow><mtext>high</mtext></mrow></msubsup></mrow></math> is given by <math display="block" id="M93"><mrow><msubsup><mi>τ</mi><mrow><mtext>bv</mtext><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow><mrow><mtext>high</mtext></mrow></msubsup><mo>=</mo><msubsup><mover accent="true"><mi>τ</mi><mo>¯</mo></mover><mrow><mtext>bv</mtext><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow><mrow><mtext>high</mtext></mrow></msubsup><msup><mrow><mrow><mo>(</mo><mrow><mi>π</mi><mi>G</mi><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mo>/</mo><msup><mi>Ω</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><msubsup><mi>T</mi><mn>9</mn><mn>2</mn></msubsup><msubsup><mi>m</mi><mrow><mn>100</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msubsup><mn>,</mn></mrow></math> with the prefactor <math id="M94"><mrow><msubsup><mover accent="true"><mi>τ</mi><mo>¯</mo></mover><mrow><mtext>bv</mtext><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow><mrow><mtext>high</mtext></mrow></msubsup><msup><mrow><mo>=</mo><mn>2.03</mn><mo>×</mo><mn>10</mn></mrow><mrow><mo>−</mo><mn>27</mn></mrow></msup><msub><mi>α</mi><mrow><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow></msub><msup><mi>β</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><msubsup><mi>μ</mi><mi>d</mi><mn>3</mn></msubsup><msubsup><mi>m</mi><mi>s</mi><mrow><mo>−</mo><mn>4</mn></mrow></msubsup></mrow></math> . Here, T 9 is the temperature of strange stars in units of 10 9 K, and m 100 is the mass of strange quarks in units of 100 MeV. At the low- T limit, <math id="M95"><mrow><msubsup><mi>τ</mi><mrow><mtext>bv</mtext></mrow><mrow><mtext>low</mtext></mrow></msubsup></mrow></math> is given by <math display="block" id="M96"><mrow><msubsup><mover accent="true"><mi>τ</mi><mo>¯</mo></mover><mrow><mtext>bv</mtext><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow><mrow><mtext>low</mtext></mrow></msubsup><mo>=</mo><msubsup><mover accent="true"><mi>τ</mi><mo>¯</mo></mover><mrow><mtext>bv</mtext><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow><mrow><mtext>low</mtext></mrow></msubsup><msup><mrow><mrow><mo>(</mo><mrow><mi>π</mi><mi>G</mi><mover accent="true"><mi>ρ</mi><mo>¯</mo></mover><mo>/</mo><msup><mi>Ω</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><msubsup><mi>T</mi><mn>9</mn><mrow><mo>−</mo><mn>2</mn></mrow></msubsup><msubsup><mi>m</mi><mrow><mn>100</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msubsup><mn>,</mn></mrow></math> with the prefactor <math id="M97"><mrow><msubsup><mover accent="true"><mi>τ</mi><mo>¯</mo></mover><mrow><mtext>bv</mtext></mrow><mrow><mtext>low</mtext></mrow></msubsup><mo>=</mo><mn>9.44</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>24</mn></mrow></msup><msub><mi>α</mi><mrow><mo stretchy="false">(</mo><mo>∥</mo><mn>,</mn><mo>⊥</mo><mo stretchy="false">)</mo></mrow></msub><msubsup><mi>μ</mi><mi>d</mi><mrow><mo>−</mo><mn>3</mn></mrow></msubsup><msubsup><mi>m</mi><mi>s</mi><mrow><mo>−</mo><mn>4</mn></mrow></msubsup></mrow></math> . The parameters α and β in these prefactors depend on the expression of the bulk viscosity and can be written in cgs units as <math display="block" id="M98"><mrow><msub><mi>α</mi><mo>∥</mo></msub><mo>=</mo><mn>9.39</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>22</mn></mrow></msup><mrow><mo>(</mo><mrow><mfrac><mrow><mn>144</mn></mrow><mn>5</mn></mfrac><msup><mi>π</mi><mn>3</mn></msup><msub><mi>μ</mi><mi>d</mi></msub><mi>Σ</mi><msubsup><mi>C</mi><mo>∥</mo><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mo stretchy="false">(</mo><mtext>g</mtext><mo>⋅</mo><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>⋅</mo><msup><mtext>s</mtext><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M99"><mrow><msub><mi>α</mi><mo>⊥</mo></msub><mo>=</mo><mn>9.39</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>22</mn></mrow></msup><mrow><mo>(</mo><mrow><mfrac><mrow><mn>144</mn></mrow><mn>5</mn></mfrac><msup><mi>π</mi><mn>3</mn></msup><msub><mi>μ</mi><mi>d</mi></msub><mi>Σ</mi><msub><mi>C</mi><mo>∥</mo></msub><msub><mi>C</mi><mo>⊥</mo></msub></mrow><mo>)</mo></mrow><mo stretchy="false">(</mo><mtext>g</mtext><mo>⋅</mo><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>⋅</mo><msup><mtext>s</mtext><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow></math> and <math display="block" id="M100"><mrow><mi>β</mi><mo>=</mo><mn>7.11</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>4</mn></mrow></msup><mrow><mo>(</mo><mrow><mfrac><mrow><mn>576</mn></mrow><mrow><mn>25</mn></mrow></mfrac><msup><mi>π</mi><mn>2</mn></msup><msubsup><mi>μ</mi><mi>d</mi><mn>2</mn></msubsup><msup><mi>Σ</mi><mn>2</mn></msup><msup><mi>W</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo stretchy="false">(</mo><msup><mtext>s</mtext><mrow><mo>−</mo><mn>2</mn></mrow></msup><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where <math id="M101"><mrow><msub><mi>C</mi><mo>∥</mo></msub></mrow></math> , <math id="M102"><mrow><msub><mi>C</mi><mo>⊥</mo></msub></mrow></math> , and W are given in Eq. (11), and ∑ is related to λ in Eq. (13) as <math display="block" id="M103"><mrow><mi>Σ</mi><mo>=</mo><mi>λ</mi><mo>/</mo><mrow><mo>(</mo><mrow><mfrac><mrow><mn>64</mn></mrow><mn>5</mn></mfrac><msup><mi>π</mi><mn>5</mn></msup><msup><mover accent="true"><mi>G</mi><mo>˜</mo></mover><mn>2</mn></msup><msub><mi>μ</mi><mi>d</mi></msub><msup><mi>T</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> Figure 6 Figure 6 shows the r -mode instability window for a typical compact star with mass <math id="M104"><mrow><mi>M</mi><mo>=</mo><mn>1.4</mn><mtext> </mtext><msub><mtext>M</mtext><mo>⊙</mo></msub></mrow></math> and radius R =10 km. The observational data (solid dots with error bars) of the spin frequency ν =Ω/2 π and internal temperature T of compact stars in LMXBs are also given for comparison[ 60 60 ]. The resulting instability window is in very good agreement with the observational data. All the stars appear in the stable region (the region below each curve). Compared with that of non-interacting SQM with low magnetic field strength (dotted lines), the stability window for interacting SQM with high magnetic field strength (solid lines) is much larger and yields a minimum rotation frequency that exceeds 1050 Hz, which may explain the recently observed very high spin frequency of a pulsar with ν =1122 Hz[ 61 61 ]. Moreover, a comparison of the dashed lines (non-interacting SQM with high magnetic field strength) and dotted lines (non-interacting SQM with low magnetic field strength) reveals that although a strong magnetic field can enlarge the instability window, compact stars in LMXBs are still located well within the stable region. 4　summary 1 The bulk viscosity of interacting magnetized SQM was investigated using the equivparticle model. First, it was found that regardless of whether interactions are included, <math id="M105"><mrow><msub><mi>P</mi><mo>∥</mo></msub></mrow></math> and <math id="M106"><mrow><msub><mi>P</mi><mo>⊥</mo></msub></mrow></math> become distinguishable when <math id="M107"><mrow><msub><mi>B</mi><mtext>m</mtext></msub><mo>≳</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>18</mn></mrow></msup></mrow></math> G. Second, compared with that of non-interacting SQM, the magnetization <math id="M108"><mi mathvariant="script">M</mi></math> is significantly decreased by the effects of interquark interactions. In addition, the β -equilibrium condition can modify the behavior of <math id="M109"><mi mathvariant="script">M</mi></math> when B m is extremely high. Then, the anisotropic bulk viscosities were studied at varying baryon number densities n b and magnetic fields B m . The results showed that when B m is sufficiently large, <math id="M110"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> can be more susceptible to interactions than <math id="M111"><mrow><msub><mi>ζ</mi><mo>∥</mo></msub></mrow></math> , and the negative <math id="M112"><mrow><msub><mi>ζ</mi><mo>⊥</mo></msub></mrow></math> can be greatly improved by interquark interactions, which may result in stable QSs with strong magnetic fields. Moreover, the quasi-oscillation of the bulk viscosities with changes in density may facilitate the formation of magnetic domains, which may complicate the magnetic field distribution so that only the averaged bulk viscosity has practical meaning for the large-scale behavior of matter over some range of magnetic fields. Finally, the resulting r -mode instability window for a typical compact star with mass <math id="M113"><mrow><mi>M</mi><mo>=</mo><mn>1.4</mn><mtext> </mtext><msub><mtext>M</mtext><mo>⊙</mo></msub></mrow></math> and radius R =10 km was presented. The r -mode instability windows are in good agreement with the observational data for compact stars in LMXBs. In particular, the instability window for interacting SQM with a high magnetic field has a minimum rotation frequency exceeding 1050 Hz, which may explain the observed very high spin frequency of a pulsar with ν =1122 Hz. References B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , Observation of Gravitational Waves from a Binary Black Hole Merger . Phys. Rev. Lett. 116 , 061102 ( 2016 ). doi: 10.1103/PhysRevLett.116.061102 http://doi.org/10.1103/PhysRevLett.116.061102 Observation of Gravitational Waves from a Binary Black Hole Merger B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs . Phys. Rev. X 9 , 031040 ( 2019 ). doi: 10.1103/PhysRevX.9.031040 http://doi.org/10.1103/PhysRevX.9.031040 GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral . Phys. Rev. Lett. 119 , 161101 ( 2017 ). doi: 10.1103/PhysRevLett.119.161101 http://doi.org/10.1103/PhysRevLett.119.161101 GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , Multi-messenger Observations of a Binary Neutron Star Merger . Astrophys. J. 848 , L12 ( 2017 ). doi: 10.3847/2041-8213/aa91c9 http://doi.org/10.3847/2041-8213/aa91c9 Multi-messenger Observations of a Binary Neutron Star Merger A. Goldstein , P. Veres , E. Burns , M.S. Briggs et al. , An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi -GBM Detection of GRB 170817A . Astrophys. J. 848 , L14 ( 2017 ). doi: 10.3847/2041-8213/aa8f41 http://doi.org/10.3847/2041-8213/aa8f41 An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi -GBM Detection of GRB 170817A B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , GW170817: measurements of neutron star radii and equation of state . Phys. Rev. Lett. 121 , 161101 ( 2018 ). doi: 10.1103/PhysRevLett.121.161101 http://doi.org/10.1103/PhysRevLett.121.161101 GW170817: measurements of neutron star radii and equation of state E. Annala , T. Gorda , A. Kurkela et al. , Gravitational-wave constraints on the neutron-star-matter equation of state . Phys. Rev. Lett. 120 , 172703 ( 2018 ). doi: 10.1103/PhysRevLett.120.172703 http://doi.org/10.1103/PhysRevLett.120.172703 Gravitational-wave constraints on the neutron-star-matter equation of state S. De , D. Finstad , J.M. Lattimer et al. , Tidal deformabilities and radii of neutron stars from the observation of GW170817 . Phys. Rev. Lett. 121 , 091102 ( 2018 ). doi: 10.1103/PhysRevLett.121.091102 http://doi.org/10.1103/PhysRevLett.121.091102 Tidal deformabilities and radii of neutron stars from the observation of GW170817 Z.-Y. Zhu , E.-P. Zhou , A. Li , Neutron star equation of state from the quark level in light of GW170817 . Astrophys. J. 862 , 98 ( 2018 ). doi: 10.3847/1538-4357/aacc28 http://doi.org/10.3847/1538-4357/aacc28 Neutron star equation of state from the quark level in light of GW170817 M. Fasano , T. Abdelsalhin , A. Maselli et al. , Constraining the Neutron Star Equation of State Using Multiband Independent Measurements of Radii and Tidal Deformabilities . Phys. Rev. Lett. 123 , 141101 ( 2019 ). doi: 10.1103/PhysRevLett.123.141101 http://doi.org/10.1103/PhysRevLett.123.141101 Constraining the Neutron Star Equation of State Using Multiband Independent Measurements of Radii and Tidal Deformabilities P.D. Lasky , Gravitational Waves from Neutron Stars: A Review . Pub. Astron. Soc. Aust. 32 , e034 ( 2015 ). doi: 10.1017/pasa.2015.35 http://doi.org/10.1017/pasa.2015.35 Gravitational Waves from Neutron Stars: A Review B.P. Abbott , R. Abbott , T.D. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) , All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO O2 data . Phys. Rev. D 100 , 024004 ( 2019 ). doi: 10.1103/PhysRevD.100.024004 http://doi.org/10.1103/PhysRevD.100.024004 All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO O2 data S. Caride , R. Inta , B.J. Owen et al. , How to search for gravitational waves from r -modes of known pulsars . Phys. Rev. D 100 , 064013 ( 2019 ). doi: 10.1103/PhysRevD.100.064013 http://doi.org/10.1103/PhysRevD.100.064013 How to search for gravitational waves from r-modes of known pulsars B. Haskell , K. Schwenzer , Gravitational waves from isolated neutron stars . arXiv:2104.03137 ( 2021 ). Gravitational waves from isolated neutron stars T. Boztepe , E. Göğüs , T. Güver et al. , Strengthening the bounds on the r -mode amplitude with X-ray observations of millisecond pulsars . Mon. Not. Roy. Astron. Soc. 498 , 2734 ( 2020 ). doi: 10.1093/mnras/staa2503 http://doi.org/10.1093/mnras/staa2503 Strengthening the bounds on the r-mode amplitude with X-ray observations of millisecond pulsars Y.-B. Wang , X. Zhou , N. Wang et al. , The r -mode instability windows of strange stars . Res. Astron. Astrophys. 19 , 030 ( 2019 ). doi: 10.1088/1674-4527/19/2/30 http://doi.org/10.1088/1674-4527/19/2/30 The r-mode instability windows of strange stars S. Chandrasekhar , Solutions of Two Problems in the Theory of Gravitational Radiation . Phys. Rev. Lett. 24 , 611 ( 1970 ). doi: 10.1103/PhysRevLett.24.611 http://doi.org/10.1103/PhysRevLett.24.611 Solutions of Two Problems in the Theory of Gravitational Radiation J.L. Friedman , B.F. Schutz , Lagrangian perturbation theory of nonrelativistic fluids . Astrophys. J. 221 , 937 ( 1978 ). doi: 10.1086/156098 http://doi.org/10.1086/156098 Lagrangian perturbation theory of nonrelativistic fluids J. Lin , R.-J. Lu , Signature of r -mode Gravitational-wave Emission in the X-Ray Afterglow of Short GRB 090510 . Astrophys. J. 871 , 160 ( 2019 ). doi: 10.3847/1538-4357/aaf72f http://doi.org/10.3847/1538-4357/aaf72f Signature of r-mode Gravitational-wave Emission in the X-Ray Afterglow of Short GRB 090510 D.D. Ofengeim , M.E. Gusakov , P. Haensel et al. , Bulk viscosity in neutron stars with hyperon cores . Phys. Rev. D 100 , 103017 ( 2019 ). doi: 10.1103/PhysRevD.100.103017 http://doi.org/10.1103/PhysRevD.100.103017 Bulk viscosity in neutron stars with hyperon cores D.D. Ofengeim , M.E. Gusakov , P. Haensel et al. , R -mode stabilization in neutron stars with hyperon cores . Journal of Physics: Conference Series 1400 , 022029 ( 2019 ). doi: 10.1088/1742-6596/1400/2/022029 http://doi.org/10.1088/1742-6596/1400/2/022029 R-mode stabilization in neutron stars with hyperon cores C.-M. Pi , S.-H. Yang , X.-P. Zheng , R -mode instability of strange stars and observations of neutron stars in LMXBs . Res. Astro. Astrophys. 15 , 871 ( 2015 ). doi: 10.1088/1674-4527/15/6/009 http://doi.org/10.1088/1674-4527/15/6/009 R-mode instability of strange stars and observations of neutron stars in LMXBs B. Haskell , N. Degenaar , W.C.G. Ho , Constraining the physics of the r -mode instability in neutron stars with X-ray and ultraviolet observations . Mon. Not. Roy. Astro. Soc. 424 , 93 ( 2012 ). doi: 10.1111/j.1365-2966.2012.21171.x http://doi.org/10.1111/j.1365-2966.2012.21171.x Constraining the physics of the r-mode instability in neutron stars with X-ray and ultraviolet observations R. Bierkandt , C. Manuel , Bulk viscosity coefficients due to phonons and kaons in superfluid color-flavor locked quark matter . Phys. Rev. D 84 , 023004 ( 2011 ). doi: 10.1103/PhysRevD.84.023004 http://doi.org/10.1103/PhysRevD.84.023004 Bulk viscosity coefficients due to phonons and kaons in superfluid color-flavor locked quark matter M.G. Alford , S. Mahmoodifar , K. Schwenzer , Viscous damping of r -modes: Small amplitude instability . Phys. Rev. D 85 , 024007 ( 2012 ). doi: 10.1103/PhysRevD.85.024007 http://doi.org/10.1103/PhysRevD.85.024007 Viscous damping of r-modes: Small amplitude instability M.G. Alford , K. Schwenzer , What the Timing of Millisecond Pulsars Can Teach us about Their Interior . Phys. Rev. Lett. 113 , 251102 ( 2014 ). doi: 10.1103/PhysRevLett.113.251102 http://doi.org/10.1103/PhysRevLett.113.251102 What the Timing of Millisecond Pulsars Can Teach us about Their Interior M.G. Alford , K. Schwenzer , Probing dense matter in compact star cores with radio pulsar data . Nucl. Phys. A 931 , 740 ( 2014 ). doi: 10.1016/j.nuclphysa.2014.08.077 http://doi.org/10.1016/j.nuclphysa.2014.08.077 Probing dense matter in compact star cores with radio pulsar data X.-G. Huang , M. Huang , D.H. Rischke et al. , Anisotropic hydrodynamics, bulk viscosities, and r -modes of strange quark stars with strong magnetic fields . Phys. Rev. D 81 , 045015 ( 2010 ). doi: 10.1103/PhysRevD.81.045015 http://doi.org/10.1103/PhysRevD.81.045015 Anisotropic hydrodynamics, bulk viscosities, and r-modes of strange quark stars with strong magnetic fields K. Schwenzer , How long-range interactions tune the damping in compact stars . arXiv:1212.5242 ( 2012 ). How long-range interactions tune the damping in compact stars J.-F. Xu , Y.-A. Luo , L. Li et al. , Bulk viscosity of strange quark matter in an enhanced perturbative QCD model . Phys. Rev. D 96 , 063016 ( 2017 ). doi: 10.1103/PhysRevD.96.063016 http://doi.org/10.1103/PhysRevD.96.063016 Bulk viscosity of strange quark matter in an enhanced perturbative QCD model X.-P. Zheng , M. Kang , X.-W. Liu et al. , Running coupling constant from lattice data and bulk viscosity of strange quark matter . Phys. Rev. C 72 , 025809 ( 2005 ). doi: 10.1103/PhysRevC.72.025809 http://doi.org/10.1103/PhysRevC.72.025809 Running coupling constant from lattice data and bulk viscosity of strange quark matter X.-P. Zheng , X.-W. Liu , M. Kang et al. , Bulk viscosity of strange quark matter in a density-dependent quark mass model and dissipation of the r mode in strange stars . Phys. Rev. C 70 , 015803 ( 2004 ). doi: 10.1103/PhysRevC.70.015803 http://doi.org/10.1103/PhysRevC.70.015803 Bulk viscosity of strange quark matter in a density-dependent quark mass model and dissipation of the r mode in strange stars X.-P. Zheng , S.-H. Yang , J.-R. Li et al. , Bulk viscosity of interacting strange quark matter . Phys. Lett. B 548 , 29 ( 2002 ). doi: 10.1016/S0370-2693(02)02840-X http://doi.org/10.1016/S0370-2693(02)02840-X Bulk viscosity of interacting strange quark matter A. Goyal , V. Gupta , Pragya et al. , Bulk viscosity of strange quark matter . Z. Phys. A 349 , 93 ( 1994 ). doi: 10.1007/BF01296338 http://doi.org/10.1007/BF01296338 Bulk viscosity of strange quark matter J.-F. Xu , D.-B. Kang , G.-X. Peng et al. , Bulk viscosity for interacting strange quark matter and r -mode instability windows for strange stars . Chin. Phys. C 45 , 015103 ( 2021 ). doi: 10.1088/1674-1137/abc0cd http://doi.org/10.1088/1674-1137/abc0cd Bulk viscosity for interacting strange quark matter and r-mode instability windows for strange stars P.-C. Chu , L.-W Chen , Quark matter symmetry energy and quark stars . AstroPhys. J. 780 , 135 ( 2014 ). doi: 10.1088/0004-637X/780/2/135 http://doi.org/10.1088/0004-637X/780/2/135 Quark matter symmetry energy and quark stars H. Yu , D.-Q. Fang , Y.-G Ma , Investigation of the symmetry energy of nuclear matter . Nucl. Sci. Tech. 31 , 61 ( 2020 ). doi: 10.1007/s41365-020-00766-x http://doi.org/10.1007/s41365-020-00766-x Investigation of the symmetry energy of nuclear matter S. Chakrabarty , Quark matter in a strong magnetic field . Phys. Rev. D 54 , 1306 ( 1996 ). doi: 10.1103/PhysRevD.54.1306 http://doi.org/10.1103/PhysRevD.54.1306 Quark matter in a strong magnetic field P.-C. Chu , L.-W. Chen , X. Wang , Quark stars in strong magnetic fields . Phys. Rev. D 90 , 063013 ( 2014 ). doi: 10.1103/PhysRevD.90.063013 http://doi.org/10.1103/PhysRevD.90.063013 Quark stars in strong magnetic fields P.-C. Chu , X.-H. Li , H.-Y. Ma et al. , Quark matter and quark stars in strong magnetic fields at finite temperature within the confined-isospin-density-dependent mass model . Phys. Lett. B 778 , 447 ( 2018 ). doi: 10.1016/j.physletb.2018.01.064 http://doi.org/10.1016/j.physletb.2018.01.064 Quark matter and quark stars in strong magnetic fields at finite temperature within the confined-isospin-density-dependent mass model X.-J. Wen , S.-Z. Su , D.-H. Yang et al. , Magnetized strange quark matter in a quasiparticle description . Phys. Rev. D 86 , 034006 ( 2012 ). doi: 10.1103/PhysRevD.86.034006 http://doi.org/10.1103/PhysRevD.86.034006 Magnetized strange quark matter in a quasiparticle description S.-S. Cui , G.-X. Peng , Z.-Y. Lu et al. , Properties of color-flavor locked strange quark matter in an external strong magnetic field . Nucl. Sci. Tech. 26 , 040503 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.040503 http://doi.org/10.13538/j.1001-8042/nst.26.040503 Properties of color-flavor locked strange quark matter in an external strong magnetic field C. Peng , G. X. Peng , C.-J. Xia et al. , Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions . Nucl. Sci. Tech. 27 , 98 ( 2016 ). doi: 10.1007/s41365-016-0095-5 http://doi.org/10.1007/s41365-016-0095-5 Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions C. Kouveliotou , U. D. Desai , T. L. Cline et al. , A 2.2 Second Period in the 1984 August 3 Gamma-Ray Burst . Astrophys. J. 330 , L101 - L105 ( 1988 ). doi: 10.1086/185214 http://doi.org/10.1086/185214 A 2.2 Second Period in the 1984 August 3 Gamma-Ray Burst R. C. Duncan , C. Thompson , Formation of Very Strongly Magnetized Neutron Stars: Implications for Gamma-Ray Bursts . Astrophys. J. 392 , L9 - L13 ( 1992 ). doi: 10.1086/186413 http://doi.org/10.1086/186413 Formation of Very Strongly Magnetized Neutron Stars: Implications for Gamma-Ray Bursts E. J. Ferrer , V. de la Incera , J. P. Keith et al. , Equation of state of a dense and magnetized fermion system . Phys. Rev. C 82 , 065802 ( 2010 ). doi: 10.1103/PhysRevC.82.065802 http://doi.org/10.1103/PhysRevC.82.065802 Equation of state of a dense and magnetized fermion system X. Zhou , X.P. Zheng , Non-leptonic Weak Interaction with Quark Matter in Magnetic Field . High Energy Physics and Nuclear Physics 31 , 1003 - 1009 ( 2007 ). Non-leptonic Weak Interaction with Quark Matter in Magnetic Field J. D. Anand , V. K. Gupta , A. Goyal et al. , Bulk viscosity of magnetized neutron star matter . J. Phys. G 27 , 921 - 931 ( 2001 ). doi: 10.1088/0954-3899/27/4/315 http://doi.org/10.1088/0954-3899/27/4/315 Bulk viscosity of magnetized neutron star matter C. J. Xia , G. X. Peng , S. W. Chen et al. , Thermodynamic consistency, quark mass scaling, and properties of strange matter . Phys. Rev. D 89 , 105027 ( 2014 ). doi: 10.1103/PhysRevD.89.105027 http://doi.org/10.1103/PhysRevD.89.105027 Thermodynamic consistency, quark mass scaling, and properties of strange matter S. Chakrabarty , D. Bandyopadhyay , S. Pal , Weakly interacting quark matter in an ultrastrong magnetic field . Int. J. Mod. Phys. A 13 , 295 - 314 ( 1998 ). doi: 10.1142/S0217751X98000135 http://doi.org/10.1142/S0217751X98000135 Weakly interacting quark matter in an ultrastrong magnetic field X.-P. Zheng , X. Zhou , and X.-W. Liu , Non-leptonic weak interaction in magnetized quark matter . arXiv:astro-ph/0507390 ( 2005 ). Non-leptonic weak interaction in magnetized quark matter J. Madsen , Probing strange stars and color superconductivity by r -mode instabilities in millisecond pulsars . Phys. Rev. Lett. 85 , 10 ( 2000 ). doi: 10.1103/PhysRevLett.85.10 http://doi.org/10.1103/PhysRevLett.85.10 Probing strange stars and color superconductivity by r-mode instabilities in millisecond pulsars L. Bildsten , G. Ushomirsky , Viscous boundary-layer damping of r -modes in neutron stars . Astrophys. J. 529 , L33 - L36 ( 2000 ). doi: 10.1086/312454 http://doi.org/10.1086/312454 Viscous boundary-layer damping of r-modes in neutron stars N. Andersson , D.I. Jones , K.D. Kokkotas et al. , r -Mode Runaway and Rapidly Rotating Neutron Stars . Astrophys. J. 534 , L75 - L78 ( 2000 ) doi: 10.1086/312643 http://doi.org/10.1086/312643 r-Mode Runaway and Rapidly Rotating Neutron Stars J. Madsen , Intermediate mass strangelets are positively charged . Phys. Rev. Lett. 85 , 4687 - 4690 ( 2000 ). doi: 10.1103/PhysRevLett.85.4687 http://doi.org/10.1103/PhysRevLett.85.4687 Intermediate mass strangelets are positively charged X.-P. Zheng , S.-H. Yang , J.-R. Li , The critical rotation of strange stars and rapidly rotating pulsars . Astrophys. J. Lett. 585 , L135 - L138 ( 2003 ). doi: 10.1086/374370 http://doi.org/10.1086/374370 The critical rotation of strange stars and rapidly rotating pulsars L. Lindblom , G. Mendell , B.J. Owen , Second-order rotational effects on the r-modes of neutron stars . Phys. Rev. D 60 , 064006 ( 1999 ). doi: 10.1103/PhysRevD.60.064006 http://doi.org/10.1103/PhysRevD.60.064006 Second-order rotational effects on the r-modes of neutron stars K.D. Kokkotas , N. Stergioulas , Analytic description of the r - mode instability in uniform density stars . Astron. Astrophys. 341 , 110 ( 1999 ). Analytic description of the r - mode instability in uniform density stars H. Heiselberg , C.J. Pethick , Transport and relaxation in degenerate quark plasmas . Phys. Rev. D 48 , 2916 ( 1993 ). doi: 10.1103/PhysRevD.48.2916 http://doi.org/10.1103/PhysRevD.48.2916 Transport and relaxation in degenerate quark plasmas M.E. Gusakov , A.I. Chugunov , E.M. Kantor , Instability windows and evolution of rapidly rotating neutron stars . Phys. Rev. Lett. 112 , 151101 ( 2014 ). doi: 10.1103/PhysRevLett.112.151101 http://doi.org/10.1103/PhysRevLett.112.151101 Instability windows and evolution of rapidly rotating neutron stars P. Kaaret , Z. Prieskorn , J.J.M. In’t Zand et al. , Evidence of 1122 Hz X-Ray Burst Oscillations from the Neutron Star X-Ray Transient XTE J1739-285 . Astrophys. J. 657 , L97 - L100 ( 2007 ). doi: 10.1086/513270 http://doi.org/10.1086/513270 Evidence of 1122 Hz X-Ray Burst Oscillations from the Neutron Star X-Ray Transient XTE J1739-285 10.1007/s41365-021-00954-3 This work was supported by the National Natural Science Foundation of China (Nos. 12005005, 11947098) and Key Research Projects of Universities in Henan Province (No. 20A140003). 2021-05-15 2021-08-05 2021-08-06 2021-10-11 2021-10-01 2021-11-29 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2021 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-021-00954-3 http://dx.doi.org/10.1007/s41365-021-00954-3 . 2021 Nuclear Science and Techniques 10 32