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

Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions 1 first-author 1 2 3 corresp Corresponding author, gxpeng@ucas.ac.cn Corresponding author, gxpeng@ucas.ac.cn gxpeng@ucas.ac.cn 4 1 1 School of Physics, University of Chinese Academy of Sciences, 380 Huaibeizhen, Beijing 101408, China School of Physics, University of Chinese Academy of Sciences, 380 Huaibeizhen, Beijing 101408, China Theoretical Physics Center for Science Facilities, Institute of High Energy Physics, Beijing 100049, China Theoretical Physics Center for Science Facilities, Institute of High Energy Physics, Beijing 100049, China Synergetic Innovation Center for Quantum Effects and Application, Hunan Normal University, Changsha 410081, China Synergetic Innovation Center for Quantum Effects and Application, Hunan Normal University, Changsha 410081, China Key Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, Beijing 100190, China Key Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, Beijing 100190, China 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 scaling Magnetic field Quark matter Compact star 1 Introduction 1 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 1 ]. Since the early works of many authors [ 2 2 - 4 4 2-4 ], especially Witten’s conjecture that SQM might be absolutely stable [ 5 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 6 ]. Although many meaningful works have been done in the past decades, there are still many aspects left open [ 7 7 ]. For example, the stability and equation of state has not yet been definitively fixed and are still under active investigation [ 8 8 ]. A magnetic field has strong effects on the stability and properties of SQM [ 9 9 , 10 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 12 G [ 11 11 ]. The magnetic field in a magnetar could be as high as 10 13 -10 15 G or even higher [ 12 12 , 13 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 13 ]. Another model demonstrates magnetic flux conservation: the relatively small magnetic fields were amplified during the star collapse [ 14 14 , 15 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 1 , 16 16 ], conventional perturbation models [ 17 17 - 19 19 17-19 ], the Richard potential model [ 20 20 , 21 21 ], the Nambu and Jona-Lasinio model [ 22 22 ], the field correlator method [ 23 23 ], the quark-cluster model [ 24 24 ], the model with magnetic-field dependent coupling [ 25 25 ], the unification way to describe both slets and strange stars [ 26 26 , 27 27 ], and many other models [ 28 28 - 39 39 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 40 - 50 50 40-50 ], the other is the density and/or temperature dependent [ 51 51 - 54 54 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 8 , 55 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 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 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 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 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 51 ] for light quarks, and later extended to strange quarks [ 52 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 35 ]. This cubic scaling was extended to a finite temperature [ 36 36 ]. It has also been extensively applied to the investigations of SQM-related physics, e.g. the QCD phase diagram [ 55 55 ], properties of SQM and slets [ 56 56 ] at zero and finite temperatures [ 36 36 , 57 57 ], the damping time scale due to the coupling of the viscosity and r mode [ 58 58 ], and the quark-diquark equation of state and compact star structure [ 59 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 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 1 Let’s start with the energy of a relativistically moving free particle, i , with mass, m i , i.e. <math display="block" id="M1"><mrow><msub><mi>ε</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><msubsup><mi>p</mi><mi>z</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>p</mi><mo>⊥</mo><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt><mn>,</mn></mrow></math> where m i is the particle mass, and p z is the momentum along the z direction while p ⊥ is the momentum component in the x-y plane or perpendicular to the z axis: <math id="M2"><mrow><msubsup><mi>p</mi><mo>⊥</mo><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>p</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>p</mi><mi>y</mi><mn>2</mn></msubsup></mrow></math> . The corresponding thermodynamic potential density is <math display="block" id="M3"><mrow><msub><mi>Ω</mi><mtext>f</mtext></msub><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mfrac><mrow><mn>2</mn><msub><mi>g</mi><mi>i</mi></msub></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mtext>​</mtext><mtext>​</mtext><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mtext>​</mtext><mtext>​</mtext><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mrow><mo>(</mo><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>ε</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mtext> </mtext><mi>d</mi><msub><mi>p</mi><mi>x</mi></msub><mi>d</mi><msub><mi>p</mi><mi>y</mi></msub><mi>d</mi><msub><mi>p</mi><mi>z</mi></msub><mn>,</mn></mrow></math> where g i 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 q i ( 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 p z can be still taken as a continuous variable, the component p ⊥ perpendicular to the z direction should be quantized according to Sommerfield’s rule: <math display="block" id="M4"><mrow><msubsup><mi>p</mi><mo>⊥</mo><mn>2</mn></msubsup><mo>=</mo><mn>2</mn><mi>ν</mi><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><msub><mi>B</mi><mtext>m</mtext></msub><mn>,</mn></mrow></math> where the integer quantum number, ν , is <math display="block" id="M5"><mrow><mi>ν</mi><mo>≡</mo><mi>l</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn><mo>−</mo><mtext> </mtext><mi>s</mi><mi>g</mi><mi>n</mi><mtext> </mtext><mo stretchy="false">(</mo><msub><mi>q</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mi>S</mi><mn>,</mn></mrow></math> 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 <math display="block" id="M6"><mrow><msub><mi>ε</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msqrt><mrow><msubsup><mi>p</mi><mi>z</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><mi>ν</mi><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><msub><mi>B</mi><mtext>m</mtext></msub></mrow></msqrt><mo>≡</mo><msqrt><mrow><msubsup><mi>p</mi><mi>z</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><mn>,</mn></mrow></math> where <math display="block" id="M7"><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><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><msub><mi>B</mi><mtext>m</mtext></msub></mrow></msqrt><mn>.</mn></mrow></math> For a given ν , one has two sets of ( l , S ) pair satisfying Eq. (4): <math display="block" id="M8"><mrow><mi>l</mi><mo>=</mo><mi>ν</mi><mn>,</mn><mi>S</mi><mo>=</mo><mtext>sign</mtext><mo stretchy="false">(</mo><msub><mi>q</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow></math> <math display="block" id="M9"><mrow><mi>l</mi><mo>=</mo><mi>ν</mi><mo>−</mo><mn>1,</mn><mi>S</mi><mo>=</mo><mo>−</mo><mtext> </mtext><mtext>sign</mtext><mtext> </mtext><mo stretchy="false">(</mo><msub><mi>q</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mn>2.</mn></mrow></math> 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 ⊥ d p ⊥ =| q i | B m d ν . Accordingly, we have ∫∫d p x d p y =∫∫ p ⊥d p ⊥d θ =2 π | q i | B m ∫ d ν , where θ is the direction angle of <math id="M10"><mrow><msub><mover accent="true"><mi>p</mi><mo>→</mo></mover><mo>⊥</mo></msub></mrow></math> . For the quantized case, the following substitution rule is thus obvious: <math display="block" id="M11"><mrow><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>​</mtext><mtext>​</mtext><mtext>​</mtext><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext> </mtext><mtext>d</mtext><msub><mi>p</mi><mi>x</mi></msub><mtext> </mtext><mtext>d</mtext><msub><mi>p</mi><mi>y</mi></msub><mo>=</mo><mn>2</mn><mi>π</mi><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><msub><mi>B</mi><mtext>m</mtext></msub><mstyle displaystyle="true"><munder><mo>∑</mo><mi>ν</mi></munder></mstyle><mn>.</mn></mrow></math> 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: <math display="block" id="M12"><mtable columnalign="left"><mtr><mtd><msub><mi>Ω</mi><mtext>M</mtext></msub><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><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"><munder><mo>∑</mo><mi>ν</mi></munder></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><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mo stretchy="false">(</mo><msub><mi>μ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>ε</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">)</mo><mtext>d</mtext><msub><mi>p</mi><mi>z</mi></msub></mtd></mtr><mtr><mtd><mo>≡</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mstyle displaystyle="true"><munder><mo>∑</mo><mi>ν</mi></munder></mstyle><mfrac><mrow><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mrow><mo>(</mo><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>−</mo><msqrt><mrow><msubsup><mi>p</mi><mi>z</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><mo>)</mo></mrow><mtext> </mtext><mtext>d</mtext><msub><mi>p</mi><mi>z</mi></msub><mn>,</mn></mtd></mtr></mtable></math> where we have used f i , ν g i | q i | B m (2- δ ν ,0 )/(2 π 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 p z at a given ν . The normal way is to give the up limit of the respective single particle energy, i.e. the Fermi energy ϵ i , <math display="block" id="M13"><mrow><msqrt><mrow><msubsup><mi>p</mi><mi>z</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><mo>≤</mo><msub><mi>ϵ</mi><mi>i</mi></msub><mn>.</mn></mrow></math> In this case, p z is in the range <math id="M14"><mrow><mo>|</mo><msub><mi>p</mi><mi>z</mi></msub><mo>|</mo><mo>≤</mo><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></math> . Therefore Ω M becomes <math display="block" id="M15"><mrow><msub><mi>Ω</mi><mtext>M</mtext></msub><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><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><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><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></msubsup></mrow></mstyle><mo stretchy="false">[</mo><msub><mi>μ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>ε</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mtext> </mtext><mi>d</mi><mtext> </mtext><msub><mi>p</mi><mi>z</mi></msub><mn>.</mn></mrow></math> The Fermi energy should be linked to the chemical potential by minimizing the thermodynamic potential density with ∂Ω f /∂є i =0 which requires <math display="block" id="M16"><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><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>μ</mi><mi>i</mi></msub><mo>−</mo><msub><mi>ϵ</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mtext> </mtext><mfrac><mrow><msub><mi>ϵ</mi><mi>i</mi></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><mo>=</mo><mn>0.</mn></mrow></math> This equation immediately gives є i = μ i , or the Fermi momentum is <math id="M17"><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mn>2</mn></msubsup><mo>−</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></math> . 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, m i , is a constant. Now to consider medium effects, we let the quark masses be density dependent, i.e. <math display="block" id="M18"><mrow><msub><mi>m</mi><mi>i</mi></msub><mo>=</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><mn>.</mn></mrow></math> We know from a recent study [ 8 8 ] that the original chemical potentials in this case should be replaced by the corresponding effective ones, i.e. <math id="M19"><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>→</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> . By this substitution process, we immediately have <math display="block" id="M20"><mrow><msub><mi>Ω</mi><mn>0</mn></msub><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><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><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></msubsup></mrow></mstyle><mo stretchy="false">[</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>−</mo><msub><mi>ε</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>p</mi><mi>z</mi></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mtext> </mtext><mi>d</mi><msub><mi>p</mi><mi>z</mi></msub><mn>.</mn></mrow></math> The maximum ν value is limited to make all <math id="M21"><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></math> meaningful. For this we should take <math display="block" id="M22"><mrow><msub><mi>ν</mi><mrow><mtext>max</mtext></mrow></msub><mo>=</mo><mtext> </mtext><mi>i</mi><mi>n</mi><mi>t</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><mo>−</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mrow><mn>2</mn><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><msub><mi>B</mi><mtext>m</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where the ’int’ means taking the integer part. Explicitly carrying out the integration gives <math display="block" id="M23"><mtable columnalign="left"><mtr><mtd><msub><mi>Ω</mi><mn>0</mn></msub><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><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><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow><mn>2</mn></mfrac><mo stretchy="false">[</mo><mo>|</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>|</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><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><mrow><mrow><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup><mi>ln</mi><mfrac><mrow><mo>|</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>|</mo><mo>+</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><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></mrow><mo>]</mo></mrow><mn>.</mn></mtd></mtr></mtable></math> Then the particle number densities are obtained from <math id="M24"><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> , giving <math display="block" id="M25"><mrow><msub><mi>n</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><mo>|</mo><msub><mi>q</mi><mi>i</mi></msub><mo>|</mo><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><mi>ν</mi><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>i</mi><mn>,0</mn></mrow></msub><mo stretchy="false">)</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><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> The real chemical potential, μ i , is connected to the effective one, <math id="M26"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> , by <math display="block" id="M27"><mrow><msub><mi>μ</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></mfrac><mo>≡</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>−</mo><msub><mi>μ</mi><mtext>I</mtext></msub><mn>,</mn></mrow></math> 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 <math display="block" id="M28"><mrow><msub><mi>μ</mi><mtext>I</mtext></msub><mo>=</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>−</mo><msub><mi>μ</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfrac><mrow><mtext>d</mtext><msub><mi>m</mi><mtext>I</mtext></msub></mrow><mrow><mtext>d</mtext><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow></mfrac><mn>,</mn></mrow></math> where the summed derivative with respect to the quark masses is <math display="block" id="M29"><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow></mfrac><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><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><msub><mi>m</mi><mi>i</mi></msub><mtext> </mtext><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub><mtext> </mtext><mtext> </mtext><mi>a</mi><mi>r</mi><mi>c</mi><mi>c</mi><mi>o</mi><mi>s</mi><mi>h</mi><mtext> </mtext><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></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> 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 <math id="M30"><mrow><mi>E</mi><mo>=</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><msub><mo>∑</mo><mi>i</mi></msub><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msub><mi>n</mi><mi>i</mi></msub></mrow></math> : <math display="block" id="M31"><mrow><mi>E</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>ν</mi><mo>=</mo><mn>0</mn></mrow></munder></mstyle><mfrac><mrow><msub><mi>f</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow></msub></mrow><mn>2</mn></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><mo>−</mo><msubsup><mi>M</mi><mrow><mi>i</mi><mn>,</mn><mi>ν</mi></mrow><mn>2</mn></msubsup></mrow></msqrt></mrow></msubsup></mrow></mstyle><msqrt><mrow><msubsup><mi>p</mi><mi>z</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><mtext>d</mtext><msub><mi>p</mi><mi>z</mi></msub><mn>.</mn></mrow></math> After completing the integration, we can get <math display="block" id="M32"><mtable columnalign="left"><mtr><mtd><mi>E</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><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><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><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><mo>*</mo></msubsup><msup><mrow/><mn>2</mn></msup><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> 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 <math display="block" id="M33"><mtable columnalign="left"><mtr><mtd><mi>P</mi><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mn>,</mn><mi>j</mi></mrow></munder></mstyle><msub><mi>n</mi><mi>i</mi></msub><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><msub><mi>n</mi><mtext>b</mtext></msub><mfrac><mrow><mtext>d</mtext><msub><mi>m</mi><mtext>I</mtext></msub></mrow><mrow><mtext>d</mtext><msub><mi>n</mi><mtext>b</mtext></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow></mfrac><mn>.</mn></mtd></mtr></mtable></math> Comparing Eq. (24) with the constant-mass case where P =-Ω 0 , 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 1 In the equivparticle model or the previous quark-mass-density-dependent model, the quark mass can be divided into two parts, i.e. <math display="block" id="M34"><mrow><msub><mi>m</mi><mi>i</mi></msub><mo>=</mo><msub><mi>m</mi><mrow><mn>0,</mn><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>m</mi><mtext>I</mtext></msub><mn>,</mn></mrow></math> 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 51 ]. Later, it was suggested to be inversely proportional to the cubic root of the density based on in-medium chiral condensates [ 35 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 8 ] as <math display="block" id="M35"><mrow><msub><mi>m</mi><mtext>I</mtext></msub><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><mo>+</mo><mi>C</mi><msubsup><mi>n</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mn>.</mn></mrow></math> 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 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 60 ]. This way, the quark mass formula becomes <math display="block" id="M36"><mrow><msub><mi>m</mi><mtext>I</mtext></msub><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><mo>+</mo><mi>C</mi><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>n</mi><mtext>b</mtext></msub><mo>−</mo><mi>w</mi></mrow><mrow><msub><mi>n</mi><mtext>a</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><msubsup><mi>n</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mn>,</mn></mrow></math> 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) <math display="block" id="M37"><mrow><msub><mi>μ</mi><mtext>u</mtext></msub><mo>+</mo><msub><mi>μ</mi><mtext>e</mtext></msub><mo>=</mo><msub><mi>μ</mi><mtext>d</mtext></msub><mo>=</mo><msub><mi>μ</mi><mtext>s</mtext></msub><mn>.</mn></mrow></math> Because the difference between the real and effective chemical potentials is the same for all quark favors, the effective chemical potentials <math id="M38"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> ( i = u,d,s) satisfy similar equalities as in Eq. (28): <math display="block" id="M39"><mrow><msubsup><mi>μ</mi><mtext>u</mtext><mo>*</mo></msubsup><mo>+</mo><msub><mi>μ</mi><mtext>e</mtext></msub><mo>=</mo><msubsup><mi>μ</mi><mtext>d</mtext><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>μ</mi><mtext>s</mtext><mo>*</mo></msubsup><mn>.</mn></mrow></math> 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 <math display="block" id="M40"><mrow><msub><mi>n</mi><mtext>b</mtext></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>(</mo><mrow><msub><mi>n</mi><mtext>u</mtext></msub><mo>+</mo><msub><mi>n</mi><mtext>d</mtext></msub><mo>+</mo><msub><mi>n</mi><mtext>s</mtext></msub></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> At the same time, according to the charge of each flavor of quarks, the total electric charge density is <math display="block" id="M41"><mrow><mi>Q</mi><mo>=</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><msub><mi>n</mi><mtext>u</mtext></msub><mo>−</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><msub><mi>n</mi><mtext>d</mtext></msub><mo>−</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><msub><mi>n</mi><mtext>s</mtext></msub><mo>−</mo><msub><mi>n</mi><mtext>e</mtext></msub><mn>.</mn></mrow></math> 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, <math id="M42"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> ( 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 Fig. 1 , we show the particle distribution, n i /(3 n 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, n i , ( 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. The corresponding quark chemical potentials and their effective ones, μ i and <math id="M43"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> ( i = u, d, s), are given in Fig. 2 Fig. 2 . In Fig. 3 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 18 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. 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 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 Fig. 3 with D 1/2 to be, respectively, 140 MeV, 133 MeV, and 129 MeV. 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 31 , 39 39 , 61 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., <math id="M44"><mrow><msubsup><mi>B</mi><mtext>m</mtext><mn>2</mn></msubsup><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><mi>E</mi><mo stretchy="false">)</mo></mrow></math> . In Fig. 5 Fig. 5 , we give the maximum value of this ratio, i.e., <math id="M45"><mrow><msubsup><mi>B</mi><mtext>m</mtext><mn>2</mn></msubsup><mo>/</mo><mo stretchy="false">(</mo><mn>2</mn><msub><mi>E</mi><mrow><mtext>min</mtext></mrow></msub><mo stretchy="false">)</mo></mrow></math> , where E min is the energy density corresponding to the minimum-energy state (the solids dots in Fig. 3 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 18 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 18 Gauss. 4 MASS-RADIUS RELATION OF MAGNETIZED STRANGE QUARK STAR 1 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 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 <math display="block" id="M47"><mrow><mfrac><mrow><mi>d</mi><mi>P</mi></mrow><mrow><mi>d</mi><mi>r</mi></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mi>G</mi><mi>m</mi><mi>E</mi></mrow><mrow><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac><mfrac><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>P</mi><mo>/</mo><mi>E</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup><msup><mi>r</mi><mn>3</mn></msup><mi>P</mi><mo>/</mo><mi>m</mi><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>−</mo><mn>2</mn><mi>G</mi><mi>m</mi><mo>/</mo><mi>r</mi></mrow></mfrac><mn>,</mn></mrow></math> with the subsidiary condition <math display="block" id="M48"><mrow><mi>m</mi><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>r</mi></msubsup></mrow></mstyle><mn>4</mn><mi>π</mi><mi>E</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><msup><mi>r</mi><mn>2</mn></msup><mtext> </mtext><mi>d</mi><mtext> </mtext><mi>r</mi><mn>.</mn></mrow></math> 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 62 ]. In Fig. 6 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. 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 63 , 64 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 1 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 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 http://doi.org/10.1103/PhysRevD.92.025025 Strange matter and strange stars in a thermodynamically self-consistent perturbation model with running coupling and running strange quark mass N. Itoh . Hydrostatic Equilibrium of Hypothetical Quark Stars . Prog. Theor. Phys . 44 , 291 ( 1970 ). doi: 10.1143/PTP.44.291 http://doi.org/10.1143/PTP.44.291 Hydrostatic Equilibrium of Hypothetical Quark Stars A.R. Bodmer . Collapsed nuclei . Phys. Rev. D . 4 , 1601 ( 1971 ). doi: 10.1103/PhysRevD.4.1601 http://doi.org/10.1103/PhysRevD.4.1601 Collapsed nuclei H. Terazawa . INS-Report 336 , Univ. of Tokyo ( 1979 ). INS-Report 336 E. Witten . Cosmic separation of phases . Phys. Rev. D . 30 , 272 ( 1984 ). doi: 10.1103/PhysRevD.30.272 http://doi.org/10.1103/PhysRevD.30.272 Cosmic separation of phases M. Angeles Perez-Garcia , J. Silk and J.R. Stone . Dark Matter, Neutron Stars, and Strange Quark Matter . Phys. Rev. Lett ., 105 , 141101 ( 2010 ). doi: 10.1103/PhysRevLett.105.141101 http://doi.org/10.1103/PhysRevLett.105.141101 Dark Matter, Neutron Stars, and Strange Quark Matter V. Dexheimer , J.R. Torres and D. Menezes . Stability windows for proto-quark stars . Eur. Phys. J. C , 73 , 2569 ( 2013 ). doi: 10.1140/epjc/s10052-013-2569-5 http://doi.org/10.1140/epjc/s10052-013-2569-5 Stability windows for proto-quark stars C.J. Xia , G.X. Peng , S.W. Chen , Z.Y. Lu and J.F. Xu . 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 . 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 K. Sato and H. Suzuki . Analysis of neutrino burst from the supernova 1987A in the Large Magellanic Cloud . Phys. Rev. Lett ., 58 , 2722 ( 1987 ). doi: 10.1103/PhysRevLett.58.2722 http://doi.org/10.1103/PhysRevLett.58.2722 Analysis of neutrino burst from the supernova 1987A in the Large Magellanic Cloud L. Dong and S.L. Shapiro . Cold equation of state in a strong magnetic field - Effects of inverse beta-decay . Astrophys. J ., 383 , 745 ( 1991 ). doi: 10.1086/170831 http://doi.org/10.1086/170831 Cold equation of state in a strong magnetic field - Effects of inverse beta-decay R. Duncan and C. Thompson . Formation of very strongly magnetized neutron stars-implications for gamma-ray bursts . Astron. J ., 392 , L9 ( 1992 ). doi: 10.1086/186413 http://doi.org/10.1086/186413 Formation of very strongly magnetized neutron stars-implications for gamma-ray bursts C. Kouveliotou et al ., An X - ray pulsar with a superstrong magnetic field in the soft gamma-ray repeater SGR 1806-20 . Nature . 393 , 235 ( 1988 ). doi: 10.1038/30410 http://doi.org/10.1038/30410 An X - ray pulsar with a superstrong magnetic field in the soft gamma-ray repeater SGR 1806-20 T. Tatsumi , T. Maruyama , E. Nakano and K. Nawa . Ferromagnetism in quark matter and origin of the magnetic field in compact stars . Nucl. Phys. A . 774 , 827 ( 2006 ). doi: 10.1016/j.nuclphysa.2006.06.145 http://doi.org/10.1016/j.nuclphysa.2006.06.145 Ferromagnetism in quark matter and origin of the magnetic field in compact stars B. Feng , R. Huang and Y. Jia . Gauge Amplitude Identities by On-shell Recursion Relation in S-matrix Program . Phys. Lett. B . 695 , 350 ( 2011 ). doi: 10.1016/j.physletb.2010.11.011 http://doi.org/10.1016/j.physletb.2010.11.011 Gauge Amplitude Identities by On-shell Recursion Relation in S-matrix Program G.X. Peng . Thermodynamic correction to the strong interaction in the perturbative regime . Europhys. Lett ., 72 ( 1 ), 69 ( 2005 ). doi: 10.1209/epl/i2005-10189-8 http://doi.org/10.1209/epl/i2005-10189-8 Thermodynamic correction to the strong interaction in the perturbative regime E.S. Fraga , R.D. Pisarski and J. Schaffner-Bielich . Small, dense quark stars from perturbative QCD . Phys. Rev. D . 63 , 121702(R) ( 2001 ). doi: 10.1103/PhysRevD.63.121702 http://doi.org/10.1103/PhysRevD.63.121702 Small, dense quark stars from perturbative QCD B.A. Freedman and L.M. McLerran . Fermions and Gauge Vector Mesons at Finite Temperature and Density. 1. Formal Techniques . Phys. Rev. D . 16 , 1130 ( 1977 ). doi: 10.1103/PhysRevD.16.1130 http://doi.org/10.1103/PhysRevD.16.1130 Fermions and Gauge Vector Mesons at Finite Temperature and Density. 1. Formal Techniques V. Baluni . Nonabelian Gauge Theories of Fermi Systems: Chromotheory of Highly Condensed Matter . Phys. Rev. D . 17 , 2092 ( 1978 ). doi: 10.1103/PhysRevD.17.2092 http://doi.org/10.1103/PhysRevD.17.2092 Nonabelian Gauge Theories of Fermi Systems: Chromotheory of Highly Condensed Matter M. Bagchi , S. Ray , M. Dey and J. Dey . Compact strange stars with a medium dependence in gluons at finite temperature . Astron. Astrophys . 450 , 431 ( 2006 ). doi: 10.1051/0004-6361:20053732 http://doi.org/10.1051/0004-6361:20053732 Compact strange stars with a medium dependence in gluons at finite temperature M. Sinha , X.G. Huang and A. Sedrakian . Strange quark matter in strong magnetic fields within a confining model . Phys. Rev. D . 88 , 025008 ( 2013 ). doi: 10.1103/PhysRevD.88.025008 http://doi.org/10.1103/PhysRevD.88.025008 Strange quark matter in strong magnetic fields within a confining model Y. Nambu and G. Jona-Lasino . Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity . Phys. Rev . 122 , 345 ( 1961 ). doi: 10.1103/PhysRev.122.345 http://doi.org/10.1103/PhysRev.122.345 Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity S. Plumari , G.F. Bugio and V. Greco . Quark matter in Neutron Stars within the Field Correlator Method . Phys. Rev. D . 88 , 083005 ( 2013 ). doi: 10.1103/PhysRevD.88.083005 http://doi.org/10.1103/PhysRevD.88.083005 Quark matter in Neutron Stars within the Field Correlator Method R.X. Xu . Solid quark matter . Astrophys. J . 596 , L59 ( 2003 ); doi: 10.1086/379209 http://doi.org/10.1086/379209 ; Solid quark matter R.X. Xu , Can cold quark matter be solid? Int. J. Mod. Phys. D 19 , 1437 ( 2003 ). doi: 10.1142/S0218271810017767 http://doi.org/10.1142/S0218271810017767 Can cold quark matter be solid? C.F. Li , L. Yang , X.J. Wen and G.X. Peng . Magnetized quark matter with a magnetic-field dependent coupling . Phys. Rev. D . 93 , 054005 ( 2016 ). doi: 10.1103/PhysRevD.93.054005 http://doi.org/10.1103/PhysRevD.93.054005 Magnetized quark matter with a magnetic-field dependent coupling C.J. Xia , G.X. Peng , E.G. Zhao and S.G. Zhou . To appear in Phys . Rev. D . 94 ( 2016 ). To appear in Phys C.J. Xia , G.X. Peng , E.G. Zhao and S.G. Zhou . From strangelets to strange stars: A unified description . Sci. Bull ., 61 ( 2 ), 172 - 177 ( 2016 ); doi: 10.1007/s11434-015-0982-x http://doi.org/10.1007/s11434-015-0982-x From strangelets to strange stars: A unified description G.S. Khadekar and R. Wanjari . Geometry of quark and strange quark matter in higher dimensional general relativity . Int. J. Theor. Phys . 51 , 1408 ( 2012 ). doi: 10.1007/s10773-011-1016-3 http://doi.org/10.1007/s10773-011-1016-3 Geometry of quark and strange quark matter in higher dimensional general relativity P.K. Sahoo and B. Mishra . Axially Symmetric Space-time With Strange Quark Matter Attached To String Cloud In Bimetric Theory . Int. J. Pure Appl. Math . 82 , 87 ( 2013 ). Axially Symmetric Space-time With Strange Quark Matter Attached To String Cloud In Bimetric Theory X.J. Wen . Color-flavor locked strange quark matter in a strong magnetic field . Phys. Rev. D . 88 , 034031 ( 2013 ). doi: 10.1103/PhysRevD.88.034031 http://doi.org/10.1103/PhysRevD.88.034031 Color-flavor locked strange quark matter in a strong magnetic field A.A. Isayev and J. Yang . Anisotropic pressure in strange quark matter under the presence of a strong magnetic field . J. Phys. G . 40 , 035105 ( 2013 ). doi: 10.1088/0954-3899/40/3/035105 http://doi.org/10.1088/0954-3899/40/3/035105 Anisotropic pressure in strange quark matter under the presence of a strong magnetic field X.Y. Wang and I.A. Shovkovy . Bulk viscosity of spin-one color superconducting strange quark matter . Phys. Rev. D . 82 , 085007 ( 2010 ). doi: 10.1103/PhysRevD.82.085007 http://doi.org/10.1103/PhysRevD.82.085007 ; Bulk viscosity of spin-one color superconducting strange quark matter M. Huang and I.A. Shovkovy . Chromomagnetic instability in dense quark matter . Phys. Rev. D . 70 , 051501(R) ( 2004 ). doi: 10.1103/PhysRevD.70.051501 http://doi.org/10.1103/PhysRevD.70.051501 Chromomagnetic instability in dense quark matter T. Bao , G.Z. Liu , E.G. Zhao and M.F. Zhu . Self-consistently thermodynamic treatment for strange quark matter in the effective mass bag model . Eur. Phys. J. A . 38 , 287 ( 2008 ). doi: 10.1140/epja/i2008-10682-6 http://doi.org/10.1140/epja/i2008-10682-6 Self-consistently thermodynamic treatment for strange quark matter in the effective mass bag model G.X. Peng , H.C. Chiang , P.Z. Ning and B.S. Zhou . Charge and critical density of strange quark matter . Phys. Rev. C . 59 , 3452 ( 1999 ). doi: 10.1103/PhysRevC.59.3452 http://doi.org/10.1103/PhysRevC.59.3452 Charge and critical density of strange quark matter G.X. Peng , H.C. Chiang , J.J. Yang , L. Li and B. Liu . Mass formulas and thermodynamic treatment in the mass density dependent model of strange quark matter . Phys. Rev. C . 61 , 015201 ( 2000 ). doi: 10.1103/PhysRevC.61.015201 http://doi.org/10.1103/PhysRevC.61.015201 Mass formulas and thermodynamic treatment in the mass density dependent model of strange quark matter X.J. Wen , Z.H. Zhong , G.X. Peng , P.N. Shen and P.Z. Ning . Thermodynamics with density and temperature dependent particle masses and properties of bulk strange quark matter and strangelets . Phys. Rev. C . 72 , 015204 ( 2005 ). doi: 10.1103/PhysRevC.72.015204 http://doi.org/10.1103/PhysRevC.72.015204 Thermodynamics with density and temperature dependent particle masses and properties of bulk strange quark matter and strangelets J.X. Hou , G.X. Peng , C.J. Xia and J.F. Xu . Magnetized strange quark matter in a mass-density-dependent model . China Phys. C . 39 ( 1 ), 015101 ( 2015 ). doi: 10.1088/1674-1137/39/1/015101 http://doi.org/10.1088/1674-1137/39/1/015101 Magnetized strange quark matter in a mass-density-dependent model 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 A.A. Isayev . Stability of magnetized strange quark matter in the MIT bag model with a density dependent bag pressure . Phys. Rev. C . 91 , 015208 ( 2015 ). doi: 10.1103/PhysRevC.91.015208 http://doi.org/10.1103/PhysRevC.91.015208 Stability of magnetized strange quark matter in the MIT bag model with a density dependent bag pressure V. Goloviznin and H. Satz . The Refractive properties of the gluon plasma in SU(2) theory . Z. Phys. C . 57 , 671 ( 1993 ). doi: 10.1007/BF01561487 http://doi.org/10.1007/BF01561487 The Refractive properties of the gluon plasma in SU(2) theory A. Peshier , B. Kämpfer , O.P. Pavlenko and G. Soff . An Effective model of the quark-gluon plasma with thermal parton masses . Phys. Lett. B . 337 , 235 ( 1994 ). doi: 10.1016/0370-2693(94)90969-5 http://doi.org/10.1016/0370-2693(94)90969-5 An Effective model of the quark-gluon plasma with thermal parton masses M.I. Gorenstein and S.N. Yang . Gluon plasma with a medium dependent dispersion relation . Phys. Rev. D . 52 , 5206 ( 1995 ). doi: 10.1103/PhysRevD.52.5206 http://doi.org/10.1103/PhysRevD.52.5206 Gluon plasma with a medium dependent dispersion relation M. Bluhm , B. Kämpfer and G. Soff . The QCD equation of state near T(c) within a quasi-particle model . Phys. Lett. B . 620 , 131 ( 2005 ). doi: 10.1016/j.physletb.2005.05.083 http://doi.org/10.1016/j.physletb.2005.05.083 The QCD equation of state near T(c) within a quasi-particle model V.M. Bannur . Self-consistent quasiparticle model for quark-gluon plasma . Phys. Rev. C . 75 , 044905 ( 2007 ). doi: 10.1103/PhysRevC.75.044905 http://doi.org/10.1103/PhysRevC.75.044905 Self-consistent quasiparticle model for quark-gluon plasma F.G. Gardim and F.M. Steffens . Thermodynamics of Quasi-Particles at Finite Chemical Potential . Nucl. Phys. A . 825 , 222 ( 2009 ). doi: 10.1016/j.nuclphysa.2009.05.001 http://doi.org/10.1016/j.nuclphysa.2009.05.001 Thermodynamics of Quasi-Particles at Finite Chemical Potential H. Li , X.L. Luo and H.S. Zong . Bag model and quark star . Phys. Rev. D . 82 , 065017 ( 2010 ); doi: 10.1103/PhysRevD.82.065017 http://doi.org/10.1103/PhysRevD.82.065017 ; Bag model and quark star L.J. Luo , J.C. Yan , W.M. Sun and H.S. Zong . A thermodynamically consistent quasi-particle model without density-dependent infinity of the vacuum zero point energy . Eur. Phys. J. C . 73 , 2626 ( 2013 ). doi: 10.1140/epjc/s10052-013-2626-0 http://doi.org/10.1140/epjc/s10052-013-2626-0 A thermodynamically consistent quasi-particle model without density-dependent infinity of the vacuum zero point energy M. Ruggieri , P. Alba , P. Castorina et al ., Polyakov Loop and Gluon Quasiparticles in Yang-Mills Thermodynamics . Phys. Rev. D . 86 , 054007 ( 2012 ). doi: 10.1103/PhysRevD.86.054007 http://doi.org/10.1103/PhysRevD.86.054007 Polyakov Loop and Gluon Quasiparticles in Yang-Mills Thermodynamics P. Alba , W. Alberico , M. Bluhm et al ., Polyakov loop and gluon quasiparticles: A self-consistent approach to Yang CMills thermodynamics . Nucl. Phys. A . 934 , 41 ( 2015 ). doi: 10.1016/j.nuclphysa.2014.11.011 http://doi.org/10.1016/j.nuclphysa.2014.11.011 Polyakov loop and gluon quasiparticles: A self-consistent approach to Yang CMills thermodynamics C. Wu and R.L. Xu . Strange quark matter and strangelets in the improved quasiparticle model . Eur. Phys. J. A . 51 ( 10 ), 124 ( 2015 ). doi: 10.1140/epja/i2015-15124-x http://doi.org/10.1140/epja/i2015-15124-x Strange quark matter and strangelets in the improved quasiparticle model Z.Y. Lu , G.X. Peng , J.F. Xu and S.P. Zhang . Properties of quark matter in a new quasiparticle model with QCD running coupling . Sci. China Phys. Mech. Astron ., 59 ( 6 ), 662001 ( 2016 ). doi: 10.1007/s11433-015-0524-2 http://doi.org/10.1007/s11433-015-0524-2 Properties of quark matter in a new quasiparticle model with QCD running coupling G.N. Fowler , S. Raha and R.M. Weiner . Confinement and Phase Transitions . Z. Phys. C . 9 , 271 ( 1981 ). doi: 10.1007/BF01410668 http://doi.org/10.1007/BF01410668 Confinement and Phase Transitions S. Chakrabarty , S. Raha and B. Sinha . Strange Quark Matter and the Mechanism of Confinement . Phys. Lett. B . 229 , 112 ( 1989 ). doi: 10.1016/0370-2693(89)90166-4 http://doi.org/10.1016/0370-2693(89)90166-4 Strange Quark Matter and the Mechanism of Confinement O.G. Benvenuto and G. Lugones . Strange matter equation of state in the quark mass density dependent model . Phys. Rev. D . 51 , 1989 ( 1995 ); doi: 10.1103/PhysRevD.51.1989 http://doi.org/10.1103/PhysRevD.51.1989 ; Strange matter equation of state in the quark mass density dependent model G. Lugones and O.G. Benvenuto . Strange matter equation of state and the combustion of nuclear matter into strange matter in the quark mass density dependent model at T > 0 . Phys. Rev. D . 52 , 1276 ( 1995 ). doi: 10.1103/PhysRevD.52.1276 http://doi.org/10.1103/PhysRevD.52.1276 Strange matter equation of state and the combustion of nuclear matter into strange matter in the quark mass density dependent model at T > 0 G.X. Peng , H.C. Chiang , B.S. Zhou , P.Z. Ning and S.J. Luo . Thermodynamics, strange quark matter, and strange stars . Phys. Rev. C . 62 , 025801 ( 2000 ). doi: 10.1103/PhysRevC.62.025801 http://doi.org/10.1103/PhysRevC.62.025801 Thermodynamics, strange quark matter, and strange stars G.X. Peng , A. Li and U. Lombardo . Deconfinement phase transition in hybrid neutron stars from the Brueckner theory with three-body forces and a quark model with chiral mass scaling . Phys. Rev. C . 77 , 065807 ( 2008 ). doi: 10.1103/PhysRevC.77.065807 http://doi.org/10.1103/PhysRevC.77.065807 Deconfinement phase transition in hybrid neutron stars from the Brueckner theory with three-body forces and a quark model with chiral mass scaling G.X. Peng , X.J. Wen and Y.D. Chen . New solutions for the color-flavor locked strangelets . Phys. Lett. B . 633 , 314 ( 2006 ). doi: 10.1016/j.physletb.2005.11.081 http://doi.org/10.1016/j.physletb.2005.11.081 New solutions for the color-flavor locked strangelets G.X. Peng , X.J. Wen and Y.D. Chen . Charge, strangeness and radius of strangelets . J. Phys. G: Nucl. Part. Phys ., 34 , 1697 ( 2007 ). doi: 10.1088/0954-3899/34/7/010 http://doi.org/10.1088/0954-3899/34/7/010 Charge, strangeness and radius of strangelets X.P. Zheng , X.W. Liu , M. Kang and S.H. Yang . 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 G. Lugones and J.E. Horvath . Quark-diquark equation of state and compact star structure . Int. J. Mod. Phys. D . 12 , 495 ( 2003 ). doi: 10.1142/S0218271803002755 http://doi.org/10.1142/S0218271803002755 Quark-diquark equation of state and compact star structure R.D. Woods and D.S. Saxon . Diffuse Surface Optical Model for Nucleon-Nuclei Scattering . Phys. Rev , 95 , 577 ( 2003 ). doi: 10.1103/PhysRev.95.577 http://doi.org/10.1103/PhysRev.95.577 Diffuse Surface Optical Model for Nucleon-Nuclei Scattering P.C. Chu , L.W. Chen and 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 G.X. Peng , H.Q. Chiang , B.S. Zou , P.Z. Ning and S.J. Luo . Thermodynamics, strange quark matter, and strange stars . Phys. Rev. C . 62 , 025801 ( 2000 ). doi: 10.1103/PhysRevC.62.025801 http://doi.org/10.1103/PhysRevC.62.025801 Thermodynamics, strange quark matter, and strange stars P. Demorest , T. Pennucci , S. Ransom , M. Roberts and J. Hessels . Shapiro Delay Measurement of A Two Solar Mass Neutron Star . Nature (London) . 467 , 1081 ( 2010 ). doi: 10.1038/nature09466 http://doi.org/10.1038/nature09466 Shapiro Delay Measurement of A Two Solar Mass Neutron Star J. Antoniadis et al ., A Massive Pulsar in a Compact Relativistic Binary . Science . 340 , 1233232 ( 2013 ). doi: 10.1126/science.1233232 http://doi.org/10.1126/science.1233232 A Massive Pulsar in a Compact Relativistic Binary 10.1007/s41365-016-0095-5 This work was supported by the National Natural Science Foundation of China (Nos. 11135011, 11475110, and 11575190) . 2016-03-01 2016-05-13 2016-05-20 2016-07-15 2016-08-20 2021-04-11 © 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. 2016 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-016-0095-5 http://dx.doi.org/10.1007/s41365-016-0095-5 . 2016 Nuclear Science and Techniques 04 27