Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions

Cheng Peng，

Guang-Xiong Peng，

Cheng-Jun Xia，

Jian-Feng Xu，

Shi-Peng Zhang

Nuclear Science and Techniques

Vol.27, No.4

Article number 98

Published in print 20 Aug 2016

Available online 15 Jul 2016

DOI：10.1007/s41365-016-0095-5

48600

We investigated the properties of strange quark matter in an external strong magnetic field with both confinement and leading-order perturbative interactions considered. It was found that the leading-order perturbative interaction can stiffen the equation of state of magnetized quark matter, while the magnetic field lowers the minimum energy per baryon. By solving the Tolman-Oppenheimer-Volkoff equations, we obtain the internal structure of strange stars. The maximum mass of strange stars can be as large as 2 times the solar mass.

New mass scalingMagnetic fieldQuark matterCompact star

1 Introduction

Strange quark matter (SQM) is a new form of matter that contains deconfined up (u), down (d), and strange (s) quarks in β-equilibrium, with electric and color charge neutrality [1]. Since the early works of many authors [2-4], especially Witten’s conjecture that SQM might be absolutely stable [5], SQM has been an important topic in nuclear, astrophysics, and cosmology due to its far-reaching theoretical significance.

SQM may be produced in energetic heavy-ion collisions, or exist in cosmic rays and/or in the inner part of neutron stars. A neutron star could be converted to a quark star or hybrid star due to leptonic weak interactions, or seeded with slets by the self-annihilating weakly interacting massive particles [6]. Although many meaningful works have been done in the past decades, there are still many aspects left open [7]. For example, the stability and equation of state has not yet been definitively fixed and are still under active investigation [8].

A magnetic field has strong effects on the stability and properties of SQM [9, 10]. It is generally known that there exists a strong magnetic field in many compact objects. For example, the typical strength on the surface of a pulsar could be in the order of ∼10¹² G [11]. The magnetic field in a magnetar could be as high as 10¹³-10¹⁵ G or even higher [12, 13]. The origin of such strong magnetic fields is presently not well understood. A widely accepted theory is the magneto hydrodynamic dynamo mechanism, where a large magnetic field is generated by the rotating plasma of a photoneutron star [13]. Another model demonstrates magnetic flux conservation: the relatively small magnetic fields were amplified during the star collapse [14, 15].

In principle, the fundamental theory of strong interaction, i.e. quantum chromodynamics (QCD), can solve all problems in investigating SQM. Due to the non-perturbative difficulty at relatively lower densities, however, the motion equations are unable to be solved exactly from the first principle theory. Therefore, people usually resort to various models, e.g., the thermodynamically enhanced pQCD model [1, 16], conventional perturbation models [17-19], the Richard potential model [20, 21], the Nambu and Jona-Lasinio model [22], the field correlator method [23], the quark-cluster model [24], the model with magnetic-field dependent coupling [25], the unification way to describe both slets and strange stars [26, 27], and many other models [28-39].

In order to include interactions between quarks, it is useful to adopt medium-dependent quark masses. There are two kinds of such masses. One is the chemical-potential and/or temperature dependent, which is normally called the quasiparticle model [40-50], the other is the density and/or temperature dependent [51-54]. In both cases, there are a lot of discussions on consistent thermodynamic treatments. It can be generally shown that the pressure at the minimum energy per baryon should be exactly zero [8, 55]. With this standard, many phenomenological models have to be ruled out.

In the case of chemical-potential dependent masses, one can have thermodynamic consistency by adding an additional term to the thermodynamic density [41]. In the density-dependent case, it is now clear that the original chemical potentials should be replaced with effective ones when the quark masses become density dependent [55]. These kind of models were previously called the density-dependent quark mass models. But its thermodynamic treatment is inconsistent. The most recent version is called the equivparticle model [8]. It is fully thermodynamically self-consistent due to the introduction of effective chemical potentials. The concept of an effective chemical potential is also shown to be useful in the extended bag model with a density-dependent bag constant [39].

Another important issue in the equivparticle model is the quark mass scaling, i.e. how to parameterize the density dependence of the quark masses. In the early stage, people mainly emphasize quark confinement, and the interaction part in the quark mass is parameterized to be inversely proportional to the baryon number density [51] for light quarks, and later extended to strange quarks [52]. According to the in-medium chiral condensate, it was shown that the interaction quark mass should be inversely proportional to the cubic root of the density [35]. This cubic scaling was extended to a finite temperature [36]. It has also been extensively applied to the investigations of SQM-related physics, e.g. the QCD phase diagram [55], properties of SQM and slets [56] at zero and finite temperatures [36, 57], the damping time scale due to the coupling of the viscosity and r mode [58], and the quark-diquark equation of state and compact star structure [59] etc.

It is true that confinement interactions dominate at lower densities. With increasing densities, however, the perturbative interaction sets in, and becomes more and more important, and thus should be included. Recently, a new quark-mass scaling considering both the confinement and perturbative interactions was derived, and applied to the investigation of SQM and strange stars [8]. But the effect of a strong magnetic field was not considered there.

The purpose of the present paper is to study the properties of SQM in a strong magnetic field with the newly developed equivparticle model, which considers both the confinement and perturbative interactions. In the Sect. 2, we describe the thermodynamic formulas with emphasis on the thermodynamic consistency, while in Sect. 3, we study properties of magnetized SQM and present the results and discussions. Then, in Sect. 4, we study the mass-radius relation of magnetized strange stars. Section 5 is a short summary.

2 Themodynamic Formulae

Let’s start with the energy of a relativistically moving free particle, i, with mass, mi, i.e.

ε_{i} = \sqrt{p_{z}^{2} + p_{⊥}^{2} + m_{i}^{2}},

(1)

where mi is the particle mass, and pz is the momentum along the z direction while p⊥ is the momentum component in the x-y plane or perpendicular to the z axis: $p_{⊥}^{2} = p_{x}^{2} + p_{y}^{2}$ . The corresponding thermodynamic potential density is

Ω_{f} = - \sum_{i} \frac{2 g_{i}}{{(2 π)}^{3}} \int ​ ​ ​ \int ​ ​ ​ \int (μ_{i} - ε_{i}) d p_{x} d p_{y} d p_{z},

(2)

where gi is the degeneracy factor with g_u=g_d=g_s=3 and g_e=1, while the extra factor 2 takes care of the spin degeneracy.

If a constant magnetic field with strength B_m is present along the z direction, classically the movement of a fermion with a charge qi (q_u=2/3, q_d=q_s=-1/3, and q_e=-1) is a regular helix whose axis is along the z axis and the projection on the x-y plane is a circle. Although the pz can be still taken as a continuous variable, the component p⊥ perpendicular to the z direction should be quantized according to Sommerfield’s rule:

p_{⊥}^{2} = 2 ν | q_{i} | B_{m},

(3)

where the integer quantum number, ν, is

ν \equiv l + 1 / 2 - s g n (q_{i}) S,

(4)

with l being the principle quantum number, S=1/2 for spin-up particles, and S=-1/2 for spin-down particles. The sign function sgn equals to 1 with a positive argument and to -1 with a negative argument. So the energy level of the particle becomes

ε_{i, ν} (p_{z}) = \sqrt{p_{z}^{2} + m_{i}^{2} + 2 ν | q_{i} | B_{m}} \equiv \sqrt{p_{z}^{2} + M_{i, ν}^{2}},

(5)

where

M_{i, ν} \equiv \sqrt{m_{i}^{2} + 2 ν | q_{i} | B_{m}} .

(6)

For a given ν, one has two sets of (l, S) pair satisfying Eq. (4):

l = ν, S = sign (q_{i}) / 2

(7)

l = ν - 1, S = - sign (q_{i}) / 2.

(8)

But for ν=0, one has only the first pair. Therefore, the lowest energy level is singly degenerate, while others are all doubly degenerate.

Eq. (3) means p_⊥dp_⊥=|qi|B_mdν. Accordingly, we have ∫∫dpxdpy=∫∫p⊥dp⊥dθ=2π|qi|B_m ∫ dν, where θ is the direction angle of ${\vec{p}}_{⊥}$ . For the quantized case, the following substitution rule is thus obvious:

\int ​ ​ ​ \int d p_{x} d p_{y} = 2 π | q_{i} | B_{m} \sum_{ν} .

(9)

Now applying this substitution to Eq. (2), and simultaneously replacing the degeneracy fact 2 for a spin with 2-δ_ν,0, we immediately get the thermodynamic potential density for a system of free Fermions in a strong magnetic field:

\begin{array}{l} Ω_{M} = - \sum_{i} \frac{g_{i} | q_{i} | B_{m}}{4 π^{2}} \sum_{ν} (2 - δ_{ν,0}) \int (μ_{i} - ε_{i, ν}) d p_{z} \\ \equiv - \sum_{i} \sum_{ν} \frac{f_{i, ν}}{2} \int (μ_{i} - \sqrt{p_{z}^{2} + M_{i, ν}^{2}}) d p_{z}, \end{array}

(10)

where we have used f_i,ν gi|qi|B_m(2-δ_ν,0)/(2π²) to simplify notations.

For the system to be determined with a given set of chemical potentials, we need to determine the limits of the integration on pz at a given ν. The normal way is to give the up limit of the respective single particle energy, i.e. the Fermi energy ϵi,

\sqrt{p_{z}^{2} + M_{i, ν}^{2}} \leq ϵ_{i} .

(11)

In this case, pz is in the range $| p_{z} | \leq \sqrt{ϵ_{i}^{2} - M_{i, ν}^{2}}$ . Therefore Ω_M becomes

Ω_{M} = - \sum_{i} \sum_{ν = 0}^{ν_{max}} f_{i, ν} \int_{0}^{\sqrt{ϵ_{i}^{2} - M_{i, ν}^{2}}} [μ_{i} - ε_{i, ν} (p_{z})] d p_{z} .

(12)

The Fermi energy should be linked to the chemical potential by minimizing the thermodynamic potential density with ∂Ω_f/∂є_i=0 which requires

\sum_{ν = 0}^{ν_{max}} f_{i, ν} (μ_{i} - ϵ_{i}) \frac{ϵ_{i}}{\sqrt{ϵ_{i}^{2} - M_{i, ν}^{2}}} = 0.

(13)

This equation immediately gives єi=μi, or the Fermi momentum is $p_{i} = \sqrt{μ_{i}^{2} - m_{i}^{2}}$ . Then the εi in Eq. (12) should be replaced by μi.

In the above derivation, we have, in fact, implicitly assumed that the quark mass, mi, is a constant. Now to consider medium effects, we let the quark masses be density dependent, i.e.

m_{i} = m_{i} (n_{b}) .

(14)

We know from a recent study [8] that the original chemical potentials in this case should be replaced by the corresponding effective ones, i.e. $μ_{i} \to μ_{i}^{*}$ . By this substitution process, we immediately have

Ω_{0} = - \sum_{i} \sum_{ν = 0}^{ν_{max}} f_{i, ν} \int_{0}^{\sqrt{μ_{i}^{*}^{2} - M_{i, ν}^{2}}} [μ_{i}^{*} - ε_{i, ν} (p_{z})] d p_{z} .

(15)

The maximum ν value is limited to make all $\sqrt{μ_{i}^{2} - M_{i, ν}^{2}}$ meaningful. For this we should take

ν_{max} = i n t (\frac{μ_{i}^{*}^{2} - m_{i}^{2}}{2 | q_{i} | B_{m}}),

(16)

where the ’int’ means taking the integer part.

Explicitly carrying out the integration gives

\begin{array}{l} Ω_{0} = - \sum_{i} \sum_{ν = 0}^{ν_{max}} \frac{f_{i, ν}}{2} [| μ_{i}^{*} | \sqrt{μ_{i}^{*}^{2} - M_{i, ν}^{2}} \\ - M_{i, ν}^{2} \ln \frac{| μ_{i}^{*} | + \sqrt{μ_{i}^{*}^{2} - M_{i, ν}^{2}}}{M_{i, ν}}] . \end{array}

(17)

Then the particle number densities are obtained from $n_{i} = - \partial Ω_{0} / \partial μ_{i}^{*}$ , giving

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

(18)

The real chemical potential, μi, is connected to the effective one, $μ_{i}^{*}$ , by

μ_{i} = μ_{i}^{*} + \sum_{j} \frac{\partial Ω_{0}}{\partial m_{j}} \frac{\partial m_{j}}{\partial n_{i}} \equiv μ_{i}^{*} - μ_{I},

(19)

where the second term is due to the density dependence of quark masses. Because n_b=(n_u+n_d+n_s)/3, it is easy to get

μ_{I} = μ_{i}^{*} - μ_{i} = - \frac{1}{3} \frac{d m_{I}}{d n_{b}} \sum_{j} \frac{\partial Ω_{0}}{\partial m_{j}},

(20)

where the summed derivative with respect to the quark masses is

\sum_{j} \frac{\partial Ω_{0}}{\partial m_{j}} = \sum_{i} \sum_{ν = 0}^{ν_{max}} m_{i} f_{i, ν} a r c c o s h (\frac{μ_{i}^{*}}{M_{i, ν}}) .

(21)

Therefore, the interaction chemical potential is the same for all three flavor quarks. This is understandable due to the fact that the strong interaction is a color interaction.

The energy density is obtained by $E = Ω_{0} + \sum_{i} μ_{i}^{*} n_{i}$ :

E = \sum_{i} \sum_{ν = 0} \frac{f_{i, ν}}{2} \int_{0}^{\sqrt{μ_{i}^{*}^{2} - M_{i, ν}^{2}}} \sqrt{p_{z}^{2} + M_{i, ν}^{2}} d p_{z} .

(22)

After completing the integration, we can get

\begin{array}{l} E = \sum_{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}

(23)

In principle, the longitudinal pressure and transverse pressure become different due to the magnetic field. But this difference only becomes obvious at an extremely high value of the magnetic field strength. Here we give the thermodynamically self-consistent pressure in the framework of the equivparticle model

\begin{array}{l} P = - Ω_{0} + \sum_{i, j} n_{i} \frac{\partial Ω_{0}}{\partial m_{j}} \frac{\partial m_{j}}{\partial n_{i}} \\ = - Ω_{0} + n_{b} \frac{d m_{I}}{d n_{b}} \sum_{j} \frac{\partial Ω_{0}}{\partial m_{j}} . \end{array}

(24)

Comparing Eq. (24) with the constant-mass case where P=-Ω₀, we find an additional term with the density derivative. It is because of this term that we can guarantee thermodynamic consistency.

3 THE PROPERTIES OF MSQM

In the equivparticle model or the previous quark-mass-density-dependent model, the quark mass can be divided into two parts, i.e.

m_{i} = m_{0, i} + m_{I},

(25)

where m_0,i is the current mass, and the interacting part, m_I, reflects the medium effects. The quark mass plays a key role in describing quark confinement. Originally, m_I was parameterized as inversely proportional to the density [51]. Later, it was suggested to be inversely proportional to the cubic root of the density based on in-medium chiral condensates [35]. There are two forms of parameterizations related to the confinement interaction that dominates at lower densities. As is well known, perturbative interactions become more and more important with increasing density. To consider the perturbative effect, a new parametrization was recently formulated in Ref. [8] as

m_{I} = \frac{D}{n_{b}^{1 / 3}} + C n_{b}^{1 / 3} .

(26)

The new parametrization considers both confinement and perturbative interactions. But it does not show explicit asymptotic freedom. One can easily find that, with increasing density, there exist a special density value (D/C)^3/2: when the density, n_b, is lower than this value, the interaction mass, m_I, as a whole decreases while it begins to increase unlimitedly when the density exceeds this special value. Therefore, the constant, C, in the second term of Eq. (26) should be a decaying function of the density. In fact, as has already been shown in Ref. [8], the quantity C is proportional to the square root of the QCD running coupling at the first-order level in perturbative QCD. However, the higher-order effects on the C variation is presently not available.

To eliminate the catastrophic increase at an extremely high density, we phenomenologically replace the constant, C, in Eq. (26) with a decaying function of the Woods-Saxon type [60]. This way, the quark mass formula becomes

m_{I} = \frac{D}{n_{b}^{1 / 3}} + C {[1 + \exp (\frac{n_{b} - w}{n_{a}})]}^{- 1} n_{b}^{1 / 3},

(27)

where we have two new parameters: the asymptotic size, w, and diffuseness parameter, n_a. They are obviously not available from QCD presently, and have to be determined phenomenologically. In the present calculations, we take w=4 fm^-3 and n_a=0.6 fm^-3.

Because of the weak interactions, the chemical potentials, μi (i= u, d, s, e), satisfy the weak equilibrium conditions (neutrinos can hardly interplay with other particles, so they can enter and leave the system freely)

μ_{u} + μ_{e} = μ_{d} = μ_{s} .

(28)

Because the difference between the real and effective chemical potentials is the same for all quark favors, the effective chemical potentials $μ_{i}^{*}$ (i= u,d,s) satisfy similar equalities as in Eq. (28):

μ_{u}^{*} + μ_{e} = μ_{d}^{*} = μ_{s}^{*} .

(29)

Because electrons do not participate in weak interactions, the corresponding real chemical potential is equal to the effective one. It is, therefore, not necessary to distinguish them for electrons.

To be in accordance with conventional nuclear physics, the baryon number density of a SQM system is defined to be

n_{b} = \frac{1}{3} (n_{u} + n_{d} + n_{s}) .

(30)

At the same time, according to the charge of each flavor of quarks, the total electric charge density is

Q = \frac{2}{3} n_{u} - \frac{1}{3} n_{d} - \frac{1}{3} n_{s} - n_{e} .

(31)

The charge neutrality condition requires Q=0.

At a given density, n_b, we can numerically solve the equations in Eq. (29)–(31) with the help of Eq. (18) for the four effective chemical potentials, $μ_{i}^{*}$ (i = u, d, s, e). Then, all other thermodynamic quantities can be calculated according to the formulas presented in the preceding section.

In Fig. 1, we show the particle distribution, ni/(3n_b), the respective particle number density in the unit of the total particle number density, as a function of the baryon number density for the parameters C=0.6 and D^1/2=133 MeV. Because the electron content is much smaller than the quark’s, it is therefore separately plotted on the right axis. At extremely high density, the electron density tends to zero. Accordingly, all three quark fractions, ni, (i= u, d, s), approach to one third of the total baryon number density. With decreasing density, however, the difference between the quark fractions becomes gradually obvious: the d quark fraction increases, the s quark fraction decreases, while the u quark fraction is nearly constant.

Fig. 1.

Particle distributions in magnetized SQM. The quark fractions are plotted on the left axis. The electron number to the total quark number is very small, and is thus given on the right axis.

The corresponding quark chemical potentials and their effective ones, μi and $μ_{i}^{*}$ (i = u, d, s), are given in Fig. 2.

Fig. 2.

The chemical potentials of u and d quarks and the corresponding effective ones as functions of the density. The chemical potential for s quarks equal to that of d quarks, while the chemical potential of electron is the difference between u and u quarks.

In Fig. 3, the energy per baryon of the MSQM system is plotted with different values of the parameters C and D, as indicated in the legend, for the fixed magnetic field strength B_m=2 × 10¹⁸ G. On each curve, there is a minimum where the pressure calculated by Eq. (24) is exactly zero, showing the requirement of thermodynamic consistency. The position of the minimum is influenced by the parameters chosen. Generally speaking, the minimum energy per baryon will decease with increasing values of C because, in the stable regime, when there is a larger value for C, the corresponding D value can be smaller.

Fig. 3.

The energy per baryon as functions of the density for different parameters. The perturbative strength, parameter, the confinement parameter, and the magnetic filed strength are indicated in the legend, while other parameters are fixed to be m_0,u=5 MeV, m_0,d=10 MeV, m_0,s=80 MeV, w=4 fm^-3, n_a=0.6 fm^-3.

The minimum energy per baryon is strongly influenced by the external strong magnetic field strength. If there is no magnetic field, i.e. B_m=0, one will find that the lowest energy becomes higher. With increasing magnetic field strength, the minimum energy per baryon decreases. For the parameters where SQM is stable, the MSQM becomes more stable. To explicitly show the effects of the external magnetic field to the minimum energy, we give, in Fig. 4, the variation of the minimum energy per baryon with the external magnetic field strength for the three sets of parameters in Fig. 3 with D^1/2 to be, respectively, 140 MeV, 133 MeV, and 129 MeV.

Fig. 4.

The minimum energy per baryon varies with the magnetic field strength.

It should be pointed out that the pure magnetic field contribution is not included in the present calculations. This contribution should, in principle, be accounted for in a consistent study. In fact, many interesting investigations have considered both the matter and field contributions [31, 39, 61]. Because the origin of the strong magnetic fields in compact stars is still under discussion, the present paper is mainly interested in the effects on the mass-radius relation of MSQM stars when the asymptotic freedom is respected in the new quark mass scaling with perturbative interactions. To have an impression on the validity range of the present calculation, we use the relative importance ratio of the energy contribution from the pure magnetic field to that of the magnetized SQM, i.e., $B_{m}^{2} / (2 E)$ . In Fig. 5, we give the maximum value of this ratio, i.e., $B_{m}^{2} / (2 E_{min})$ , where E_min is the energy density corresponding to the minimum-energy state (the solids dots in Fig. 3) for each set of parameters. It is obvious that the ratio increases by increasing the magnetic field strength. At the present range of interest for B_m < 2× 10¹⁸ Gauss, however, the importance ratio is smaller than the ten percent at the minimum energy per baryon where the relative importance of the magnetic field contribution reaches its maximum. With increasing density, such as in a compact star, this ratio becomes much smaller. Therefore, the present calculations should be valid when the magnetic field strength does not obviously exceed the order of 10¹⁸ Gauss.

Fig. 5.

The energy density of the magnetic field to that of the magnetized SQM at the minimum energy per baryon. The ratio

B_{m}^{2} / (2 E_{min})

is obviously smaller than 0.1 if the field strength is less than 2 × 10¹⁸ Gauss.

4 MASS-RADIUS RELATION OF MAGNETIZED STRANGE QUARK STAR

Neutron stars are interesting astronomical compact objects consisting of neutrons. Because the inner density is high, while SQM is possible if absolutely stable, it has long been speculated that there exist objects consisting completely of quarks, the so-called strange stars. In fact, many pulsars might be candidates for strange stars [26]. The structure of strange stars plays an important role in astrophysics, cosmology, and relevant nuclear physics.

As has been done in many investigations, we need to solve the Tolman-Oppenheimer-Volkov equation

\frac{d P}{d r} = - \frac{G m E}{r^{2}} \frac{(1 + P / E) (1 + 4 π^{2} r^{3} P / m)}{1 - 2 G m / r},

(32)

with the subsidiary condition

m = \int_{0}^{r} 4 π E (r) r^{2} d r .

(33)

We now apply the equations of state obtained in this paper to solve for the mass-radius relation of magnetized strange stars. For a concrete description of the solving process, one may refer to Ref. [62].

In Fig. 6, we plot the mass-radius relation of magnetized strange stars for different sets of parameters. On each curve, there is a maximum mass indicated by a full dot. This dot separates the full cure into two parts. The quark stars on the right are mechanically stable, while those on the left are not mechanically stable, and thus surely do not exist.

Fig. 6.

The mass-radius relation of magnetized strange stars for the equations of state obtained in the last section (Fig. 3). The full dot on each cures is the maximum mass for that set of parameters.

The maximum star mass depends on the model parameters. For the parameter set with C=0.6 and D^1/2=133 MeV, the maximum star mass is close to 2 times the solar mass, while for the parameters with C=0.7 and D^1/2=129 MeV, the maximum mass can be even bigger [63, 64]. Only for the parameters C=0.4 and D^1/2=156 MeV, when MSQM is not absolutely stable, does the maximum mass become smaller.

5 Summary

We have investigated the properties of SQM in an external strong magnetic field with a new mass scaling. The new mass scaling not only considers both confinement and perturbative interactions, but also shows explicit asymptotic freedom.

The thermodynamic treatment follows the general formulas of the equivparticle model by the inclusion of effective chemical potentials. In this treatment, all the fundamental thermodynamic relations are still valid and is, thus, fully self-consistent.

With these new mass formulas and thermodynamic treatments, we studied the equation of state of magnetized quark matter and the structure of magnetized quark stars. It is found that a strong magnetic field can lower the minimum energy per baryon of MSQM. For parameters where SQM is not stable, MSQM can become absolutely stable. Because of perturbative interactions, the equations of state of MSQM becomes harder, which make the maximum star mass as large as more than two times the solar mass.

References

[1]

J.F. Xu, G.X. Peng, F. Liu, D.F. Hou and L.W. Chen.

Strange matter and strange stars in a thermodynamically self-consistent perturbation model with running coupling and running strange quark mass

. Phys. Rev. D. 92, 025025 (2015). doi: 10.1103/PhysRevD.92.025025