Bulk viscosity of interacting magnetized strange quark matter

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Bulk viscosity of interacting magnetized strange quark matter

Jian-Feng Xu

Nuclear Science and Techniques

Vol.32, No.10

Article number 111

Published in print 01 Oct 2021

Available online 11 Oct 2021

DOI：10.1007/s41365-021-00954-3

51100

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, $ζ_{∥}$ (parallel to B_m) and $ζ_{⊥}$ (perpendicular to B_m). A comparison with non-interacting strange quark matter reveals that when B_m is sufficiently large, $ζ_{⊥}$ is more affected by interactions than $ζ_{∥}$ . 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 matterBulk viscosityStrong magnetic fieldStrange starR-mode instability window

1　Introduction

Since the first direct detection of gravitational waves (GWs)[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]. 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] is extraordinarily significant, as the observation of the GWs emitted in this process, possibly combined with electromagnetic observation of the same source[4, 5], may yield insight into the structure of NSs and the equation of state (EOS) of matter under extreme conditions[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]. Continuous GWs can typically be generated by various processes that produce asymmetry[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-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, 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]. 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-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-34]. The large bulk viscosity can reduce the r-mode instability window; consequently, theoretical calculations are consistent with astrophysical observations[22]. In Ref.[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, 37] and the strong magnetic fields (on the order of approximately 10¹²-10¹³ G)[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¹⁴-10¹⁵ G[44, 45]. In fact, the largest magnetic field that can be sustained by strange stars is estimated to be approximately 1.5 × 10²⁰ G[46]. According to a previous study, SQM will be more stable when the magnetic field is included in the EOS[42]. Additionally, a strong magnetic field can strongly suppress the reaction rate of the non-leptonic weak interaction $u + d \leftrightarrow u + s$ , which is one source of the large bulk viscosity of SQM and is expected to affect the viscosity of SQM[47]. In addition, the bulk viscosity of magnetized NS matter was studied in Ref.[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, 43, 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

In the equivparticle model, the quark masses mi, (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]:

\begin{array}{l} m_{i} \equiv m_{i 0} + m_{I} (n_{b}) \\ = m_{i 0} + C n_{b}^{1 / 3} + \frac{D}{n_{b}^{1 / 3}}, \end{array}

(1)

where m_i0 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, $μ_{i} \to μ_{i}^{*}$ . 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:

\begin{array}{l} Ω_{i} = - \frac{g_{i} e_{i} B_{m}}{4 π^{2}} \sum_{ν = 0}^{ν_{max}} (2 - δ_{ν 0}) (μ_{i}^{*} \sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}} \\ - M_{i, ν}^{2} \ln \frac{μ_{i}^{*} + \sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}}}{M_{i, ν}}), \end{array}

(2)

where gi=6 is the degeneracy factor for quarks (the value is 2 for electrons), ei is the absolute value of the electric charge of the particles, $M_{i, ν} = \sqrt{m_{i}^{2} + 2 ν e_{i} B_{m}}$ , 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

Ω_{0} = \sum_{I} Ω_{i}, (i = u, d, s, and e) .

(3)

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 $M_{i} = - \partial Ω_{i} / \partial B_{m}$ :

\begin{array}{l} M_{i} = \frac{g_{i} e_{i}}{4 π^{2}} \sum_{ν = 0}^{ν_{max}} (2 - δ_{ν 0}) [μ_{i}^{*} \sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}} - (m_{i}^{2} \\ + 4 ν e_{i} B_{m}) \ln \frac{μ_{i}^{*} + \sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}}}{M_{i, ν}}] . \end{array}

(4)

Because of the strong magnetic field B_m, the longitudinal pressure and transverse pressure become different. They are given by

P_{∥} = - Ω_{0} + \sum_{i} μ_{I} n_{i}

(5)

and

P_{⊥} = - Ω_{0} + \sum_{i} μ_{I} n_{i} - M B_{m},

(6)

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:

\begin{array}{l} μ_{I} = - \frac{1}{3} \frac{\partial m_{I}}{\partial n_{b}} \frac{\partial Ω_{0}}{\partial m_{I}} \\ \equiv \sum_{i} f_{i} (n_{b}, m_{i} (n_{b}), μ_{i}^{*}, M_{i, ν} (m_{i} (n_{b}))) . \end{array}

(7)

By definition, the effective chemical potential μi^* and bare chemical potential μi can be written in terms of μ_I as follows:

μ_{i} = μ_{i}^{*} - μ_{I} .

(8)

The particle number density can be given by the conventional thermodynamic relation @ $n_{i} = - \partial Ω_{0} / \partial μ_{i}^{*}$ , which reads

n_{i} = \frac{g_{i} e_{i} B_{m}}{2 π^{2}} \sum_{ν = 0}^{ν_{max}} (2 - δ_{ν 0}) \sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}} .

(9)

Further information about the equivparticle model of magnetized SQM is given in Ref.[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] presented explicit expressions for $ζ_{⊥}$ and $ζ_{∥}$ , 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 $u + d \leftrightarrow u + s$ and were studied using a simple bag model that did not include interactions and β-equilibrium between quarks[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-33, 35].

The bulk viscosities $ζ_{∥}$ and $ζ_{⊥}$ are given as[28]

ζ_{∥} = \frac{λ C_{∥}^{2}}{ω^{2} + λ^{2} W^{2}} and ζ_{⊥} = \frac{λ C_{⊥} C_{∥}}{ω^{2} + λ^{2} W^{2}},

(10)

respectively, where

\begin{array}{l} C_{∥} ≃ n_{d} (\frac{\partial μ_{d}}{\partial n_{d}}) - n_{s} (\frac{\partial μ_{s}}{\partial n_{s}}), \\ C_{⊥} ≃ C_{∥} - X B_{m}, W = \frac{\partial μ_{s}}{\partial n_{s}} + \frac{\partial μ_{d}}{\partial n_{d}}, \end{array}

(11)

with

X = \frac{\partial M}{\partial n_{d}} - \frac{\partial M}{\partial n_{s}} .

(12)

In addition,

\begin{array}{l} λ = \frac{4 B_{m}^{4} {\tilde{G}}^{2} T^{2} e_{u}^{2} e_{d} e_{s} μ_{d}}{5 π^{3}} {(\sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{u}^{2} - M_{u, ν}^{2}}})}^{2} \\ \times \sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{d}^{2} - M_{d, ν}^{2}}} \sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{s}^{2} - M_{s, ν}^{2}}}, \end{array}

(13)

where ${\tilde{G}}^{2} \equiv G_{F}^{2} \sin^{2} θ_{C} \cos^{2} θ_{C} {= 6.46 \times 10}^{- 24} {MeV}^{- 4}$ 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, $1 s^{- 1} ≲ ω ≲ 10^{3} s^{- 1}$ . 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, $\sqrt{μ_{i}^{2} - M_{i, ν}} \to \sqrt{μ_{i}^{* 2} - M_{i, ν}}$ , 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 $u + d \leftrightarrow u + s$ , readers may refer to Ref.[50, 51].

To calculate the bulk viscosities in the equivparticle model, we must find the derivatives $\partial μ_{s} / \partial n_{s}$ and $\partial μ_{d} / \partial n_{d}$ . 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

\frac{\partial μ_{i}}{\partial n_{j}} = \frac{δ_{i j} + \frac{g_{i} e_{i} B_{m}}{2 π^{2}} \sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}}} (μ_{i}^{*} A B + \frac{m_{i}}{3} \frac{d m_{i}}{d n_{b}})}{\frac{g_{i} e_{i} B_{m}}{2 π^{2}} \sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{i}^{* 2} - M_{i, ν}^{2}}} μ_{i}^{*} A},

(14)

where

A = {(1 + \frac{\partial f_{i}}{\partial μ_{i}^{*}})}^{- 1},

(15)

and

B = \frac{1}{3} \sum_{i} (\frac{\partial f_{i}}{\partial n_{b}} + \frac{\partial f_{i}}{\partial m_{i}} \frac{d m_{i}}{d n_{b}} + \frac{\partial f_{i}}{\partial M_{i, ν}} \frac{d m_{i}}{d n_{b}} \frac{m_{i}}{M_{i, ν}}) .

(16)

As a special case, if the quark mass is independent of the baryon number density, it is not difficult to obtain

\frac{\partial μ_{i}}{\partial n_{j}} = δ_{i j} {(\frac{g_{i} e_{i} B_{m}}{2 π^{2}} \sum_{ν = 0}^{ν_{max}} \frac{2 - δ_{ν 0}}{\sqrt{μ_{i}^{2} - M_{i, ν}^{2}}} μ_{i})}^{- 1}

(17)

from Eq. (14), which is exactly the same as the results reported in Ref.[28].

3　Numerical results and discussion

To calculate the bulk viscosities, the quark current masses are set to m_u0=5 MeV, m_d0=10 MeV, and m_s0=100 MeV. In addition, the model parameters C and D are set to C=0.7 and $\sqrt{D} = 129$ MeV, which can ensure that the SQM is in the absolutely stable state[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, we show the anisotropic pressures $P_{∥}$ and $P_{⊥}$ as functions of magnetic field B_m with baryon number density n_b=2n₀, where n₀=0.17 fm^-3 is the nuclear saturation density. Regardless of whether interactions are included, when $B_{m} ≳ 10^{18}$ G, $P_{∥}$ and $P_{⊥}$ become distinguishable. Then, with increasing B_m, $P_{∥}$ first increases and then decreases. However, $P_{⊥}$ 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 $P_{∥}$ and $P_{⊥}$ are reduced.

Fig. 1

Anisotropic pressures

P_{∥}

and

P_{⊥}

as functions of magnetic field B_m. Regardless of whether interactions are included, when

B_{m} ≳ 10^{18}

P_{∥}

and

P_{⊥}

become distinguishable. Then, with increasing B_m,

P_{⊥}

decreases, and

P_{∥}

first increases and then decreases. Moreover, in contrast to the case for non-interacting SQM (C=D=0), both

P_{∥}

and

P_{⊥}

are reduced when interactions between quarks are considered (C=0.7, D^1/2=129 MeV).

Figure 2 shows the magnetization $M$ 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, the magnetization can be decreased significantly by interactions between quarks. Furthermore, for comparison, we show the results presented in Ref.[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 $μ_{u}^{*} = μ_{d}^{*} = μ_{s}^{*} = 300$ MeV. From Fig. 2, for a sufficiently large magnetic field (e.g., $B_{m} ≳ 10^{19.3}$ 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 $M$ as a function of only the chemical potentials, and thus a constant.

Fig. 2

Magnetization

M

as function of strong magnetic field B_m. Like the pressure, the magnetization can also be decreased by interactions between quarks. For comparison, we also show the results from Ref.[28] (solid line), in which the chemical potentials of u,d, and s quarks are assigned the same value.

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 $ζ_{∥}$ changes significantly, the magnitude of $ζ_{∥}$ does not change greatly on average. However, the magnitude of the transverse bulk viscosity $ζ_{⊥}$ increases significantly. Therefore, $ζ_{⊥}$ is likely to be more susceptible to interactions than $ζ_{∥}$ when B_m is sufficiently large. In fact, when interactions are not considered, the values of $ζ_{⊥}$ are even negative, which implies that QSs made of SQM in this state are hydrodynamically unstable[28]. If the interactions are included, this situation can be improved greatly. Moreover, whenever a new Landau level appears, both $ζ_{∥}$ and $ζ_{⊥}$ 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.

Fig. 3

Bulk viscosities as a function of baryon number density n_b with constant magnetic field B=10^18.5 G. Although strong interactions between quarks has little effect on

ζ_{∥}

, on average, they obviously increase

ζ_{⊥}

. In addition, 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.

Figure 4 shows $ζ_{∥}$ and $ζ_{⊥}$ as functions of magnetic field. For fixed n_b, at low magnetic field strength B_m, $ζ_{∥}$ and $ζ_{⊥}$ 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 is the negative values of $ζ_{⊥}$ (dashed line in the lower panel), which can be greatly improved by including interactions between quarks, except at extremely strong magnetic fields.

Fig. 4

Parallel bulk viscosity

ζ_{∥}

and transverse bulk viscosity

ζ_{⊥}

as functions of magnetic field strength at fixed baryon number density. When interactions are taken into account, both components of the bulk viscosity are increased, especially at low magnetic field strength.

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], where the approximate relationship between C and D can be roughly fitted as[35]

D^{1 / 2} / MeV = - \frac{27}{0.7} C + 156, C \in [0,0.7] .

(18)

According to this relationship, when C increases, D^1/2 decreases, which is shown on the upper X axis in Fig. 5. In addition, with increasing C and decreasing D, both $ζ_{∥}$ and $ζ_{⊥}$ 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 imply that confinement effects may contribute more to the bulk viscosity than perturbative interactions for the parameters n_b=2n₀ and B_m=10^18.5 G. Furthermore, the sudden decreases in both $ζ_{∥}$ and $ζ_{⊥}$ still originate from the variation of the occupied Landau levels.

Fig. 5

Bulk viscosities as functions of model parameters C and D. The relationship between C and D is constrained by the requirement of absolute stability of the SQM in the equivparticle model.

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:

\frac{1}{τ} = \frac{1}{τ_{gw}} + \frac{1}{τ_{sv}} + \frac{1}{τ_{bv}} + \dots,

(19)

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, 52-54]. Here, it should be stressed that the damping time scale of the bulk viscosity arises from both $ζ_{∥}$ and $ζ_{⊥}$ because of the magnetic field; that is,

\frac{1}{τ_{bv}} = \frac{1}{τ_{ζ_{∥}}} + \frac{1}{τ_{ζ_{⊥}}} .

(20)

Eq. (19) is solved using the method employed in Ref.[35, 55, 56]. According to the literature[57, 58], the driving time scale due to GW emission is

τ_{gw} = - 3.26 s {(π G \bar{ρ} / Ω^{2})}^{3},

(21)

where G=6.707 × 10^-45 MeV^-2 is the gravitational constant, $\bar{ρ}$ is the mean density of a compact star, and Ω is the angular rotation frequency. The damping time scale of the shear viscosity[59] is

τ_{sv} = 5.37 \times 10^{8} s {(α_{s} / 0.1)}^{5 / 3} T_{9}^{5 / 3},

(22)

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⁹ K, and the low-T limit applies for T<10⁹ K. In the equivparticle model at the high-T limit, $τ_{bv}^{high}$ is given by

τ_{bv (∥, ⊥)}^{high} = {\bar{τ}}_{bv (∥, ⊥)}^{high} {(π G \bar{ρ} / Ω^{2})}^{2} T_{9}^{2} m_{100}^{- 4},

(23)

with the prefactor ${\bar{τ}}_{bv (∥, ⊥)}^{high} {= 2.03 \times 10}^{- 27} α_{(∥, ⊥)} β^{- 1} μ_{d}^{3} m_{s}^{- 4}$ . Here, T₉ is the temperature of strange stars in units of 10⁹ K, and m₁₀₀ is the mass of strange quarks in units of 100 MeV. At the low-T limit, $τ_{bv}^{low}$ is given by

{\bar{τ}}_{bv (∥, ⊥)}^{low} = {\bar{τ}}_{bv (∥, ⊥)}^{low} {(π G \bar{ρ} / Ω^{2})}^{2} T_{9}^{- 2} m_{100}^{- 4},

(24)

with the prefactor ${\bar{τ}}_{bv}^{low} = 9.44 \times 10^{- 24} α_{(∥, ⊥)} μ_{d}^{- 3} m_{s}^{- 4}$ . The parameters α and β in these prefactors depend on the expression of the bulk viscosity and can be written in cgs units as

α_{∥} = 9.39 \times 10^{22} (\frac{144}{5} π^{3} μ_{d} Σ C_{∥}^{2}) (g \cdot {cm}^{- 1} \cdot s^{- 1}),

(25)

α_{⊥} = 9.39 \times 10^{22} (\frac{144}{5} π^{3} μ_{d} Σ C_{∥} C_{⊥}) (g \cdot {cm}^{- 1} \cdot s^{- 1})

(26)

and

β = 7.11 \times 10^{- 4} (\frac{576}{25} π^{2} μ_{d}^{2} Σ^{2} W^{2}) (s^{- 2}),

(27)

where $C_{∥}$ , $C_{⊥}$ , and W are given in Eq. (11), and ∑ is related to λ in Eq. (13) as

Σ = λ / (\frac{64}{5} π^{5} {\tilde{G}}^{2} μ_{d} T^{2}) .

(28)

Figure 6 shows the r-mode instability window for a typical compact star with mass $M = 1.4 M_{⊙}$ 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]. 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]. 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.

Fig. 6

R-mode instability window (the region below each curve) for a typical compact star with mass

M = 1.4 M_{⊙}

and radius R=10 km. Observational data on spin frequency and internal temperature of compact stars in LMXBs are also presented.

4　summary

The bulk viscosity of interacting magnetized SQM was investigated using the equivparticle model.

First, it was found that regardless of whether interactions are included, $P_{∥}$ and $P_{⊥}$ become distinguishable when $B_{m} ≳ 10^{18}$ G. Second, compared with that of non-interacting SQM, the magnetization $M$ is significantly decreased by the effects of interquark interactions. In addition, the β-equilibrium condition can modify the behavior of $M$ 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, $ζ_{⊥}$ can be more susceptible to interactions than $ζ_{∥}$ , and the negative $ζ_{⊥}$ 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 $M = 1.4 M_{⊙}$ 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

[1]

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

Bulk viscosity of interacting magnetized strange quark matter

1 Introduction

2 Bulk viscosity of magnetized SQM in equivparticle model

3 Numerical results and discussion

4 summary

1　Introduction

2　Bulk viscosity of magnetized SQM in equivparticle model

3　Numerical results and discussion

4　summary