Symmetry energy of strange quark matter and tidal deformability of strange quark stars

Symmetry energy of strange quark matter and tidal deformability of strange quark stars 1 2 corresp first-author † xujf@squ.edu.cn † xujf@squ.edu.cn xujf@squ.edu.cn https://orcid.org/0000-0003-3821-025X 3 corresp ‡ cjxia@yzu.edu.cn ‡ cjxia@yzu.edu.cn cjxia@yzu.edu.cn corresp § luzhenyan@hnust.edu.cn § luzhenyan@hnust.edu.cn luzhenyan@hnust.edu.cn https://orcid.org/0000-0002-2629-3329 5 6 corresp ¶ gxpeng@ucas.ac.cn ¶ gxpeng@ucas.ac.cn gxpeng@ucas.ac.cn 7 College of Information Engineering, Suqian University, Suqian 223800, China College of Information Engineering, Suqian University, Suqian 223800, China School of Physics and Electrical Engineering, AnYang Normal University, Anyang 455000, China School of Physics and Electrical Engineering, AnYang Normal University, Anyang 455000, China Center for Gravitation and Cosmology, College of Physical Science and Technology, Yangzhou University, Yangzhou 225009, China Center for Gravitation and Cosmology, College of Physical Science and Technology, Yangzhou University, Yangzhou 225009, China Hunan Provincial Key Laboratory of Intelligent Sensors and Advanced Sensor Materials, School of Physics and Electronics, Hunan University of Science and Technology, Xiangtan 411201, China Hunan Provincial Key Laboratory of Intelligent Sensors and Advanced Sensor Materials, School of Physics and Electronics, Hunan University of Science and Technology, Xiangtan 411201, China School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 100049, China School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 100049, China Institute of High-Energy Physics, Chinese Academy of Sciences, Beijing 100049, China Institute of High-Energy Physics, Chinese Academy of Sciences, Beijing 100049, China School of Mathematics and Physics, Henan Urban Construction University, Pingdingshan 467036, China School of Mathematics and Physics, Henan Urban Construction University, Pingdingshan 467036, China Research performed during the past decade revealed an important role of symmetry energy in the equation of state (EOS) of strange quark matter (SQM). By introducing an isospin-dependent term into the quark mass scaling, the SQM stability window in the equivparticle model was studied. The results show that a sufficiently strong isospin dependence C I can significantly widen the SQM region of absolute stability, yielding results that simultaneously satisfy the constraints of the astrophysical observations of PSR J1614-2230 with 1.928 ± 0.017 M ⊙ and tidal deformability 70 ≤ Λ 1.4 ≤ 580 measured in the event GW170817. With increasing C I , the difference between the u , d , and s quark fractions for the SQM in β -equilibrium becomes inconspicuous for C > 0, leading to small isospin asymmetry δ , and further resulting in similar EOS and structures of strange quark stars (SQSs). Moreover, unlike the behavior of the maximum mass of u – d QSs, which varies with C I depending on the sign of the parameter C , the maximum mass of the SQSs decreases monotonously with increasing C I . Symmetry energy Strange quark matter Strange stars Tidal deformability Introduction 1 At sufficiently high baryon number densities, strange quark matter (SQM) [ 1 1 - 4 4 1-4 ] with deconfined u , d , and s quarks is expected to be the true ground state of quantum chromodynamics (QCD). Owing to its possible implications in many fields, e.g., nuclear physics, astrophysics, and cosmology, it has garnered significant attention in the last few decades. Theoretically, SQM is a system dominated mainly by strong interactions, which can be studied in detail using QCD. Unfortunately, owing to the well-known quark confinement phenomenon at relatively low Densities, which could be relevant to neutron-star interiors and quark stars, perturbative QCD calculations of SQM are not reliable [ 5 5 ]. In addition, although the ab initio lattice QCD calculations were very successful in the finite-temperature regime with a zero chemical potential (baryon number density), they still suffer from the famous intractable Fermion sign problem [ 6 6 ]. Therefore, researchers typically resort to phenomenological models, such as the density-dependent model [ 7 7 - 11 11 7-11 ], quark cluster model [ 12 12 , 13 13 ], quasi-particle model [ 14 14 - 18 18 14-18 ], Nambu-Jona-Lasinio model [ 19 19 ], and so on. In Refs. [ 20 20 , 21 21 ], strong magnetic field effect on the properties of SQM had been investigated, and numerous interesting results have been obtained in other studies [ 22 22 - 25 25 22-25 ]. To study the SQM properties using phenomenological models, one of the most important issues is determining reasonable values of the model parameters. Given the conjecture that SQM is absolutely stable, the most commonly used constraint is that, for SQM, the energy per baryon should be smaller than 930 MeV, while for the bulk u - d quark matter ( u - d QM) it should be larger than 930 MeV. Moreover, significant progress has been made in the observations of global properties of compact objects, such as mass, radius, and tidal deformability. For example, the measurements of PSR J1614-2230 and PSR J0348+0432 several years ago precisely determined masses of 1.97 ± 0.04 M ⊙ [ 26 26 ] (a more precise mass measurement yielded 1.928 ± 0.017 M ⊙ [ 27 27 ] with a 68% credible level) and 2.01 ± 0.04 M ⊙ [ 28 28 ], respectively. In addition, the more recently measured gravitational mass of PSR J0740 + 6620 reached 2.14 <math id="M1"><mrow><msubsup><mrow/><mrow><mo>−</mo><mn>0.18</mn></mrow><mrow><mo>+</mo><mn>0.20</mn></mrow></msubsup></mrow></math> M ⊙ with 95.4% credibility, likely making it the the largest compact star observed up to date [ 29 29 ] (the updated result for PSR J0740 + 6620 is 2.08 <math id="M2"><mrow><msubsup><mrow/><mrow><mo>−</mo><mn>0.07</mn></mrow><mrow><mo>+</mo><mn>0.07</mn></mrow></msubsup></mrow></math> M ⊙ [ 30 30 ] with 68.3% credibility). Therefore, these observations of pulsar masses have been widely used to constrain the EOS of compact stars and the onset of the phase transition from hadronic matter to quark matter in hybrid stars [ 31 31 , 32 32 ]. In addition, with the coming of multi-messenger observations of neutron stars, owing to the detection of gravitational waves for events GW170817 [ 33 33 ] and GW190814 [ 34 34 ], additional stringent constraints on the EOS of compact star matter have been appeared [ 35 35 - 37 37 35-37 ]. Specifically, it was determined that for GW190814, the newly discovered compact binary merger possesses a secondary component with mass as large as 2.50-2.67 M ⊙ at the 90% credibility level. If it is confirmed to be a compact star, rather than a light black hole, it will be possible to rule out many EOSs. Moreover, this detection of gravitational waves for GW170817 offers a novel opportunity to directly probe the properties of matter under extreme conditions found in the interior of compact stars, yielding improved estimates of tidal deformability <math id="M3"><mrow><msub><mi>Λ</mi><mrow><mn>1.4</mn></mrow></msub><msubsup><mrow><mo>=</mo><mn>190</mn></mrow><mrow><mo>−</mo><mn>120</mn></mrow><mrow><mo>+</mo><mn>390</mn></mrow></msubsup></mrow></math> [ 38 38 ]. This measurement of tidal deformability has been widely used for constraining the properties of compact star matter [ 39 39 - 46 46 39-46 ]. Recently, an interesting proposition was made, suggesting that, instead of the SQM, the u - d QM may be the true ground state of QCD [ 47 47 ]. According to this proposition, it has been found that the properties of the u - d QS are in accordance with various astrophysical observations [ 48 48 ], where quark symmetry energy effects are included in the EOS of the u - d QM. In practice, it has been found that symmetry energy plays an important role in determining the properties of isospin asymmetry matter [ 49 49 - 53 53 49-53 ]. For nuclear matter, although the symmetry energy behaviors have been well determined below and around the saturation density, it remains a big challenge to determine the density behaviors of the symmetry energy at suprasaturation densities, where phase transitions of strong interaction matter may occur, and various new forms of matter may appear. In fact, due to the different symmetry terms in different phases, isospin effects are rather significant [ 54 54 ]. In this study, considering the isospin effects of the SQM, we investigated the stability window and symmetry energy of SQM, and employed some of the constraints mentioned above in the equivparticle model [ 23 23 ] to examine whether the SQS properties are consistent with the reported astrophysical observations. It is found that inter-quark isospin effects are necessary for reconciling the SQSs obtained in the equivparticle model with the aforementioned constraints. The remainder of this paper is organized as follows. Section 2.1 briefly describes the equivparticle model, where quark mass scaling with isospin effects is introduced. Then in Sect. 2.2, the equations for calculating the tidal deformability of compact stars are presented. In Sect. 3, the numerical results of the properties of SQM and the structures of SQSs are given and discussed. Finally, we conclude this paper with a brief summary in Sect. 4. Equations for symmetry energy and tidal deformability 1 Equivparticle model with isospin asymmetry 2 In the equivparticle model, by introducing an effective quark chemical potential <math id="M4"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> , the thermodynamic potential density of the quark matter system can be written as [ 23 23 ] <math display="block" id="M5"><mrow><msub><mi>Ω</mi><mn>0</mn></msub><mo stretchy="false">(</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mn>,</mn><msub><mi>m</mi><mi>i</mi></msub><mo stretchy="false">)</mo><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><msubsup><mi>m</mi><mi>i</mi><mn>4</mn></msubsup></mrow><mrow><mn>24</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mi>f</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where <math id="M6"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mi>x</mi><msqrt><mrow><msup><mi>x</mi><mn>2</mn></msup><mo>−</mo><mn>1</mn></mrow></msqrt><mo stretchy="false">(</mo><msup><mi>x</mi><mn>2</mn></msup><mo>−</mo><mfrac><mn>5</mn><mn>2</mn></mfrac><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>3</mn><mn>2</mn></mfrac><mtext>arch</mtext><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></math> , g i (=6 for quarks and 2 for electrons) is the particle degeneracy, m i is the isospin density-dependent equivalent quark mass, and <math id="M7"><mrow><msub><mi>ν</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt></mrow></math> is the Fermi momentum of particle species i . Using the basic thermodynamic relation <math id="M8"><mrow><msub><mi>ρ</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mtext>d</mtext><msub><mi>Ω</mi><mn>0</mn></msub><mo>/</mo><mtext>d</mtext><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> , the particle number density ρ i is readily obtained <math display="block" id="M9"><mtable><mtr><mtd><msub><mi>ρ</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><mfrac><mrow><mtext>d</mtext><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mtext>d</mtext><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac><mo>−</mo><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>ν</mi><mi>i</mi></msub></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>ν</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msubsup><mi>ν</mi><mi>i</mi><mn>3</mn></msubsup></mrow><mrow><mn>6</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub></mrow><mrow><mn>6</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><msup><mrow><mo>(</mo><mrow><msubsup><mi>μ</mi><mi>i</mi><mrow><mo>*</mo><mn>2</mn></mrow></msubsup><mo>−</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><mn>.</mn></mtd></mtr></mtable></math> Regarding the equivalent quark mass m i , different types of quark mass scaling were constructed to account for strong inter-quark interactions [ 7 7 , 10 10 , 55 55 , 23 23 ]. Conventionally, in the equivparticle model, m i is a function of the baryon number density ρ b . To include the isospin effect on the EOS of quark matter, the isospin asymmetry term <math id="M10"><mrow><mi>δ</mi><mo>≡</mo><mn>3</mn><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>−</mo><msub><mi>ρ</mi><mi>u</mi></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>u</mi></msub><mo stretchy="false">)</mo></mrow></math> is introduced into the mass scaling framework. In Ref. [ 56 56 ], m i is parameterized as <math id="M11"><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><mrow><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><mi>δ</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mi>D</mi><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>−</mo><msub><mi>τ</mi><mi>i</mi></msub><mi>δ</mi><msub><mi>D</mi><mtext>I</mtext></msub><msubsup><mi>ρ</mi><mtext>b</mtext><mi>α</mi></msubsup><msup><mtext>e</mtext><mrow><mo>−</mo><mi>β</mi><msub><mi>n</mi><mtext>b</mtext></msub></mrow></msup></mrow></math> with τ i being the third component of the isospin for quark flavor i . In this study, we used the parameterization that was recently proposed in Ref. [ 48 48 ], i.e., <math display="block" id="M12"><mtable><mtr><mtd><msub><mi>m</mi><mi>i</mi></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>m</mi><mrow><mtext>I</mtext><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mi>D</mi><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><mi>C</mi><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msub><mi>C</mi><mrow><mtext>I</mtext><mi>i</mi></mrow></msub><msup><mi>δ</mi><mn>2</mn></msup><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn></mtd></mtr></mtable></math> where m i 0 is the quark current mass, D is the confinement parameter indicating the linear confinement effect, C represents the one-gluon exchange interaction for C < 0 and the leading-order perturbative interaction for C > 0, and parameter C I i ( <math id="M13"><mrow><msub><mi>C</mi><mrow><mtext>Iu</mtext></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mtext>Id</mtext></mrow></msub><mn>,</mn><msub><mi>C</mi><mrow><mtext>Is</mtext></mrow></msub><mo>=</mo><mn>0</mn></mrow></math> ) is used for adjusting the strength of the isospin-dependent equivalent quark mass and thus the magnitude of the quark-matter symmetry energy. In addition, this equivalent quark mass satisfies the exchange symmetry between u and d quarks in quark matter in the absence of inter-quark electromagnetic interactions. The energy density and pressure of the SQM are, respectively, given as <math display="block" id="M14"><mrow><mi>E</mi><mo>=</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></mfrac><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msubsup><mi>m</mi><mi>i</mi><mn>4</mn></msubsup></mrow><mrow><mn>16</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mi>g</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow></math> where <math id="M15"><mrow><mi>g</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mi>x</mi><msqrt><mrow><msup><mi>x</mi><mn>2</mn></msup><mo>−</mo><mn>1</mn></mrow></msqrt><mo stretchy="false">(</mo><mn>2</mn><msup><mi>x</mi><mn>2</mn></msup><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mo>−</mo><mtext>arch</mtext><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow></math> and <math display="block" id="M16"><mrow><mi>P</mi><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><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><mtext>j</mtext></msub></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mn>,</mn><mi>j</mi></mrow></munder></mstyle><msub><mi>ρ</mi><mi>i</mi></msub><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>i</mi></msub></mrow></mfrac><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><mtext>j</mtext></msub></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mn>.</mn></mrow></math> Based on the definition of the isospin asymmetry δ , the derivatives of δ with respect to ρ u and ρ d are, respectively <math display="block" id="M17"><mrow><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>u</mi></msub></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mn>6</mn><msub><mi>ρ</mi><mi>d</mi></msub></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>u</mi></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mtext>and</mtext><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>d</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>6</mn><msub><mi>ρ</mi><mi>u</mi></msub></mrow><mrow><msup><mrow><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>u</mi></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mn>.</mn></mrow></math> Accordingly, the summation of the index i in the third term on the right-hand side of Eq. (4) is zero, which leads to the following simple pressure form: <math display="block" id="M18"><mrow><mi>P</mi><mo>=</mo><mo>−</mo><msub><mi>Ω</mi><mn>0</mn></msub><mo>+</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder></mstyle><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mn>,</mn></mrow></math> with <math display="block" id="M19"><mrow><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mi>D</mi><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub><msup><mrow/><mrow><mn>4</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mi>C</mi><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub><msup><mrow/><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></mfrac><mo>+</mo><msub><mi>C</mi><mrow><mtext>I</mtext><mi>i</mi></mrow></msub><msup><mi>δ</mi><mn>2</mn></msup><mn>,</mn></mrow></math> and <math display="block" id="M20"><mrow><mfrac><mrow><mo>∂</mo><msub><mi>Ω</mi><mn>0</mn></msub></mrow><mrow><mo>∂</mo><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>g</mi><mi>i</mi></msub><msub><mi>m</mi><mi>i</mi></msub></mrow><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><msub><mi>ν</mi><mi>i</mi></msub><mo>−</mo><msubsup><mi>m</mi><mi>i</mi><mn>2</mn></msubsup><mi>ln</mi><mrow><mo>(</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msub><mi>ν</mi><mi>i</mi></msub></mrow><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mn>.</mn></mrow></math> To obtain the EOS of the SQM, one has to solve the following equations, which are respectively the electric-charge neutrality conditions <math display="block" id="M21"><mrow><mn>2</mn><msub><mi>ρ</mi><mi>u</mi></msub><mo>−</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>−</mo><msub><mi>ρ</mi><mi>s</mi></msub><mo>−</mo><mn>3</mn><msub><mi>ρ</mi><mi>e</mi></msub><mo>=</mo><mn>0,</mn></mrow></math> the conservation of the baryon number density <math display="block" id="M22"><mrow><msub><mi>ρ</mi><mtext>b</mtext></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>u</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>d</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mn>,</mn></mrow></math> and the β -equilibrium condition <math display="block" id="M23"><mrow><msub><mi>μ</mi><mi>u</mi></msub><mo>+</mo><msub><mi>μ</mi><mi>e</mi></msub><mo>=</mo><msub><mi>μ</mi><mi>d</mi></msub><mo>=</mo><msub><mi>μ</mi><mi>s</mi></msub><mn>.</mn></mrow></math> It should be noted that the chemical potentials in Eq. (11) are real ones, instead of the effective ones. The relationship between the real chemical potential μ i and effective chemical potential <math id="M24"><mrow><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup></mrow></math> in the equivparticle model is <math display="block" id="M25"><mtable><mtr><mtd><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>ρ</mi><mi>i</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><msubsup><mi>μ</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></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><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>i</mi></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><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> with <math display="block" id="M26"><mrow><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mo>=</mo><mn>2</mn><msub><mi>C</mi><mrow><mi>I</mi><mi>j</mi></mrow></msub><mi>δ</mi><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>.</mn></mrow></math> Using Eq. (12, 13), Eq. (11) can be rewritten as <math display="block" id="M27"><mrow><msubsup><mi>μ</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>μ</mi><mi>u</mi><mo>*</mo></msubsup><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>u</mi></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><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mo>+</mo><msub><mi>μ</mi><mi>e</mi></msub><mn>,</mn></mrow></math> and <math display="block" id="M28"><mrow><msubsup><mi>μ</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>μ</mi><mi>d</mi><mo>*</mo></msubsup><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>δ</mi></mrow><mrow><mo>∂</mo><msub><mi>ρ</mi><mi>d</mi></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><mfrac><mrow><mo>∂</mo><msub><mi>m</mi><mi>j</mi></msub></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mn>.</mn></mrow></math> Using Eqs. (9) and (10), Eqns. (14) and (15) can be considered as a set of equations for ρ u and ρ d . By solving these equations for a given ρ b , we obtain ρ u and ρ d , and then, accordingly, the EOS of the SQM. There are two methods for obtaining the symmetry energy. The first approach uses the parabolic law, approximately given by [ 56 56 , 57 57 ] <math display="block" id="M29"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>≃</mo><mfrac><mn>1</mn><mn>9</mn></mfrac><mo stretchy="false">[</mo><mi>E</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><mi>δ</mi><mo>=</mo><mn>3,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>−</mo><mi>E</mi><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><mi>δ</mi><mo>=</mo><mn>0,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mn>,</mn></mrow></math> which implies that the symmetry energy is an estimate of the energy cost required for converting all u quarks in the symmetric u - d QM to d quarks, for given ρ b and ρ s . The other one can be analytically deduced from the following definition, i.e., <math display="block" id="M30"><mrow><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><msub><mi>E</mi><mrow><mi>u</mi><mi>d</mi><mi>s</mi></mrow></msub><mo>/</mo><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow><mrow><mo>∂</mo><msup><mi>δ</mi><mn>2</mn></msup></mrow></mfrac><msub><mo>|</mo><mrow><mi>δ</mi><mo>=</mo><mn>0</mn></mrow></msub><mn>,</mn></mrow></math> where <math display="block" id="M31"><mtable columnalign="left"><mtr><mtd><msub><mi>E</mi><mrow><mi>u</mi><mi>d</mi><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mi>g</mi><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>δ</mi><mo>/</mo><mn>3</mn><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow></msup><mi>ν</mi></mrow></msubsup></mrow></mstyle><msqrt><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>u</mi><mn>2</mn></msubsup></mrow></msqrt><msup><mi>k</mi><mn>2</mn></msup><mtext>d</mtext><mi>k</mi></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mi>g</mi><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>δ</mi><mo>/</mo><mn>3</mn><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow></msup><mi>ν</mi></mrow></msubsup></mrow></mstyle><msqrt><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>d</mi><mn>2</mn></msubsup></mrow></msqrt><msup><mi>k</mi><mn>2</mn></msup><mtext>d</mtext><mi>k</mi></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mi>g</mi><mrow><mn>2</mn><msup><mi>π</mi><mn>2</mn></msup></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mrow><msub><mi>ν</mi><mi>s</mi></msub></mrow></msubsup></mrow></mstyle><msqrt><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><mi>s</mi><mn>2</mn></msubsup></mrow></msqrt><msup><mi>k</mi><mn>2</mn></msup><mtext>d</mtext><mi>k</mi><mn>.</mn></mtd></mtr></mtable></math> Here, the Fermi momenta of u and d quarks are expressed, respectively, as <math display="block" id="M32"><mtable columnalign="left"><mtr><mtd><msub><mi>ν</mi><mi>u</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mi>δ</mi><mo>/</mo><mn>3</mn></mrow><mo>)</mo></mrow><mrow><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow></msup><mi>ν</mi><mn>,</mn></mtd></mtr><mtr><mtd><msub><mi>ν</mi><mi>d</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mi>δ</mi><mo>/</mo><mn>3</mn></mrow><mo>)</mo></mrow><mrow><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow></msup><mi>ν</mi><mn>,</mn></mtd></mtr></mtable></math> where ν is the quark Fermi momentum of the symmetric u - d QM at the quark number density. ρ =2 ρ u =2 ρ d . According to Eq. (18), Eq. (17) can be explicitly written as <math display="block" id="M33"><mtable columnalign="left"><mtr><mtd><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>9</mn><msub><mi>C</mi><mi>I</mi></msub><msub><mi>ρ</mi><mtext>b</mtext></msub><mi>m</mi></mrow><mrow><mn>2</mn><msup><mi>ν</mi><mn>3</mn></msup></mrow></mfrac><mi>A</mi><mo>+</mo><mfrac><mn>1</mn><mrow><mn>18</mn></mrow></mfrac><mfrac><mrow><msup><mi>ν</mi><mn>2</mn></msup></mrow><mrow><msqrt><mrow><msup><mi>ν</mi><mn>2</mn></msup><mo>+</mo><msup><mi>m</mi><mn>2</mn></msup></mrow></msqrt></mrow></mfrac></mrow><mo>)</mo></mrow><mfrac><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mi>s</mi></msub></mrow><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> with <math display="block" id="M34"><mrow><mi>A</mi><mo>=</mo><mi>ν</mi><msqrt><mrow><msup><mi>ν</mi><mn>2</mn></msup><mo>+</mo><msup><mi>m</mi><mn>2</mn></msup></mrow></msqrt><mo>−</mo><msup><mi>m</mi><mn>2</mn></msup><mi>a</mi><mi>r</mi><mi>c</mi><mi>s</mi><mi>h</mi><mrow><mo>(</mo><mrow><mfrac><mi>ν</mi><mi>m</mi></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where <math id="M35"><mrow><mi>m</mi><mo>=</mo><msub><mi>m</mi><mn>0</mn></msub><mo>+</mo><mi>D</mi><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msub><mi>C</mi><mrow><mi>I</mi><mi>i</mi></mrow></msub><msubsup><mi>ρ</mi><mi>b</mi><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup></mrow></math> with <math id="M36"><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>≡</mo><msub><mi>m</mi><mrow><mi>u</mi><mn>0</mn></mrow></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>d</mi><mn>0</mn></mrow></msub></mrow></math> . Here, we have assumed m u 0 = m d 0 . If m u 0 ≠ m d 0 , then <math display="block" id="M37"><mtable columnalign="left"><mtr><mtd><msub><mi>E</mi><mrow><mtext>sym</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>b</mtext></msub><mn>,</mn><msub><mi>ρ</mi><mi>s</mi></msub><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>3</mn><msub><mi>C</mi><mtext>I</mtext></msub><msub><mi>ρ</mi><mtext>b</mtext></msub><msub><mi>m</mi><msup><mi>u</mi><mo>′</mo></msup></msub></mrow><mrow><mn>2</mn><msup><mi>ν</mi><mn>3</mn></msup></mrow></mfrac><msub><mi>A</mi><mi>u</mi></msub><mo>+</mo><mfrac><mn>1</mn><mrow><mn>54</mn></mrow></mfrac><mfrac><mrow><msup><mi>ν</mi><mn>2</mn></msup></mrow><mrow><msqrt><mrow><msup><mi>ν</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><msup><mi>u</mi><mo>′</mo></msup><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac></mrow><mo>)</mo></mrow><mfrac><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mo>+</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mrow><mfrac><mrow><mn>3</mn><msub><mi>C</mi><mtext>I</mtext></msub><msub><mi>ρ</mi><mtext>b</mtext></msub><msub><mi>m</mi><msup><mi>d</mi><mo>′</mo></msup></msub></mrow><mrow><mn>2</mn><msup><mi>ν</mi><mn>3</mn></msup></mrow></mfrac><msub><mi>A</mi><mi>d</mi></msub><mo>+</mo><mfrac><mn>1</mn><mrow><mn>54</mn></mrow></mfrac><mfrac><mrow><msup><mi>ν</mi><mn>2</mn></msup></mrow><mrow><msqrt><mrow><msup><mi>ν</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><msup><mi>d</mi><mo>′</mo></msup><mn>2</mn></msubsup></mrow></msqrt></mrow></mfrac></mrow><mo>)</mo></mrow><mfrac><mrow><mn>3</mn><msub><mi>ρ</mi><mtext>b</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mi>s</mi></msub></mrow><mrow><mn>2</mn><msub><mi>ρ</mi><mtext>b</mtext></msub></mrow></mfrac><mn>,</mn></mtd></mtr></mtable></math> with <math display="block" id="M38"><mrow><msub><mi>m</mi><msup><mi>i</mi><mo>′</mo></msup></msub><mo>=</mo><msub><mi>m</mi><mrow><mi>i</mi><mn>0</mn></mrow></msub><mo>+</mo><mi>D</mi><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msub><mi>C</mi><mrow><mtext>I</mtext><mi>i</mi></mrow></msub><msubsup><mi>ρ</mi><mtext>b</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup></mrow></math> and <math display="block" id="M39"><mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><mi>ν</mi><msqrt><mrow><msup><mi>ν</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>m</mi><msup><mi>i</mi><mo>′</mo></msup><mn>2</mn></msubsup></mrow></msqrt><mo>−</mo><msubsup><mi>m</mi><msup><mi>i</mi><mo>′</mo></msup><mn>2</mn></msubsup><mtext>arcsinh</mtext><mrow><mo>(</mo><mrow><mfrac><mi>ν</mi><mrow><msub><mi>m</mi><msup><mi>i</mi><mo>′</mo></msup></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>.</mn></mrow></math> We found that the discrepancy between the symmetry energies given by Eq. (20) and Eq. (22) is essentially negligible. Therefore, in the following calculations, we adopted Eq. (20) to calculate the SQM symmetry energy. Equations for calculating structures and tidal deformability of SQSs 2 To calculate the structures and tidal deformability of SQSs, we briefly introduce relevant equations. To obtain the SQS structures, one should numerically solve the Tolman-Oppenheimer-Volkov (TOV) equation: <math display="block" id="M40"><mrow><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext> d</mtext><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><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mn>4</mn><mi>π</mi><msup><mi>r</mi><mn>3</mn></msup><mi>P</mi><mo>/</mo><mi>m</mi></mrow><mo>)</mo></mrow></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 subsidiary conditions: <math display="block" id="M41"><mrow><mfrac><mrow><mtext>d</mtext><mi>m</mi></mrow><mrow><mtext> d</mtext><mi>r</mi></mrow></mfrac><mo>=</mo><mn>4</mn><mi>π</mi><msup><mi>r</mi><mn>2</mn></msup><mi>E</mi><mn>,</mn></mrow></math> where G =6.707 × 10 -45 MeV -2 denotes the gravitational constant. Dimensionless tidal deformability is defined as [ 58 58 - 61 61 58-61 ] <math display="block" id="M42"><mrow><mi>Λ</mi><mo>≡</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><msub><mi>k</mi><mn>2</mn></msub><msup><mi>C</mi><mrow><mo>−</mo><mn>5</mn></mrow></msup><mn>,</mn></mrow></math> where <math id="M43"><mrow><mi>C</mi><mo>≡</mo><mfrac><mi>M</mi><mi>R</mi></mfrac><mi>G</mi></mrow></math> is the compactness and k 2 is the gravitational Love number with its value usually in the range of 0.2–0.3. To calculate Λ as a function of the star mass M or radius R , we need the value of k 2 , which is defined as <math display="block" id="M44"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>k</mi><mn>2</mn></msub></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><mn>8</mn><msup><mi>C</mi><mn>5</mn></msup></mrow><mn>5</mn></mfrac><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mn>2</mn><mi>C</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mrow><mo>[</mo><mrow><mn>2</mn><mo>+</mo><mn>2</mn><mi>C</mi><mrow><mo>(</mo><mrow><msub><mi>y</mi><mtext>R</mtext></msub><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>−</mo><msub><mi>y</mi><mtext>R</mtext></msub></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>×</mo><mrow><mo>{</mo><mrow><mn>2</mn><mi>C</mi><mrow><mo>[</mo><mrow><mn>6</mn><mo>−</mo><mn>3</mn><msub><mi>y</mi><mtext>R</mtext></msub><mo>+</mo><mn>3</mn><mi>C</mi><mrow><mo>(</mo><mrow><mn>5</mn><msub><mi>y</mi><mtext>R</mtext></msub><mo>−</mo><mn>8</mn></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mo>+</mo><mn>4</mn><msup><mi>C</mi><mn>3</mn></msup></mrow></mrow></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>×</mo><mrow><mo>[</mo><mrow><mn>13</mn><mo>−</mo><mn>11</mn><msub><mi>y</mi><mtext>R</mtext></msub><mo>+</mo><mi>C</mi><mrow><mo>(</mo><mrow><mn>3</mn><msub><mi>y</mi><mtext>R</mtext></msub><mo>−</mo><mn>2</mn></mrow><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msup><mi>C</mi><mn>2</mn></msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><msub><mi>y</mi><mtext>R</mtext></msub></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><mn>3</mn><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mn>2</mn><mi>C</mi><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mrow><mo>[</mo><mrow><mn>2</mn><mo>−</mo><msub><mi>y</mi><mtext>R</mtext></msub><mo>+</mo><mn>2</mn><mi>C</mi><mrow><mo>(</mo><mrow><msub><mi>y</mi><mtext>R</mtext></msub><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msup><mrow><mrow><mrow><mi>ln</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mn>2</mn><mi>C</mi><mo stretchy="false">)</mo></mrow><mo>}</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where <math id="M45"><mrow><msub><mi>y</mi><mtext>R</mtext></msub><mo>≡</mo><mi>y</mi><mo stretchy="false">(</mo><mi>R</mi><mo stretchy="false">)</mo></mrow></math> and y ( r ) is obtained as the solution to the following differential equation: <math display="block" id="M46"><mrow><mi>r</mi><mfrac><mrow><mtext>d</mtext><mi>y</mi></mrow><mrow><mtext>d</mtext><mi>r</mi></mrow></mfrac><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup><mo>+</mo><mi>y</mi><mi>F</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><msup><mi>r</mi><mn>2</mn></msup><mi>Q</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mn>0.</mn></mrow></math> Here, F ( r ) and Q ( r ) are, respectively, defined as <math display="block" id="M47"><mrow><mi>F</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>≡</mo><mfrac><mrow><mn>1</mn><mo>−</mo><mn>4</mn><mi>π</mi><msup><mi>r</mi><mn>2</mn></msup><mo stretchy="false">[</mo><mi>ϵ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>−</mo><mi>P</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mi>G</mi></mrow><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac><mn>,</mn></mrow></math> and <math display="block" id="M48"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>Q</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>≡</mo></mrow></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mn>4</mn><mi>π</mi></mrow><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac><mrow><mo>[</mo><mrow><mn>5</mn><mi>ϵ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mi>G</mi><mo>+</mo><mn>9</mn><mi>P</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mi>G</mi><mo>+</mo><mfrac><mrow><mi>ϵ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><mi>P</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><msubsup><mi>v</mi><mtext>s</mtext><mn>2</mn></msubsup><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac><mi>G</mi></mrow></mrow></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mrow><mrow><mo>−</mo><mfrac><mn>6</mn><mrow><mn>4</mn><mi>π</mi><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>]</mo></mrow><mo>−</mo><mn>4</mn><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><mi>m</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><mn>4</mn><mi>π</mi><msup><mi>r</mi><mn>3</mn></msup><mi>P</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><msup><mi>r</mi><mn>2</mn></msup><mi>f</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></mfrac><mi>G</mi></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where <math id="M49"><mrow><msubsup><mi>v</mi><mtext>s</mtext><mn>2</mn></msubsup><mo>=</mo><mo>∂</mo><mi>P</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>/</mo><mo>∂</mo><mi>ϵ</mi><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></math> denotes the squared sound velocity. In the case of a bare SQS, the energy density on the star surface is finite, which requires a correction to be added to the calculation of y R to account for the energy discontinuity, reading <math display="block" id="M50"><mrow><msub><mi>y</mi><mtext>R</mtext></msub><mo>→</mo><msub><mi>y</mi><mtext>R</mtext></msub><mo>−</mo><mfrac><mrow><mn>4</mn><mi>π</mi><msup><mi>R</mi><mn>3</mn></msup><msub><mi>ϵ</mi><mtext>s</mtext></msub></mrow><mi>M</mi></mfrac><mn>,</mn></mrow></math> where ϵ s is the energy density difference between internal and external regions. Numerical results and discussion 1 To study the isospin effects on the EOS of the SQM and the SQS structures, the values of the model parameters should be designated. First, the current masses of quarks are given as m u 0 =5 MeV, m d 0 =10 MeV, and m s 0 =100 MeV. If Eq. (20) is used for computing the symmetry energy, the values of the up and down quark current masses are m u 0 = m d 0 =7.5 MeV. The values of the parameters C , D , and C I are determined by studying the SQM stability window, combined with the astrophysical observation data of both massive compact star PSR J1614-2230 with 1.97 ± 0.04 M ⊙ and tidal deformability 70 ≤ Λ 1.4 ≤ 580 measured for GW170817. Shown in Fig. 1 Fig. 1 are the SQM stability windows for different C I , with the value labeled in each panel. The black solid line in each panel indicates 930 MeV per baryon for the SQM; the SQM becomes unstable if the energy per baryon exceeds 930 MeV. Similarly, the dashed line in each panel corresponds to the 930 MeV (per baryon) level for the u - d QM; below this level, the energy per baryon of the bulk u - d QM is smaller than 930 MeV, and nuclei are likely to decay to the u - d QM. This case contradicts the modern nuclear physics and should be forbidden. Therefore, only the area between the black solid and dashed lines is the SQM’s absolute-stability region. In addition, the blue line in each panel represents the configuration that can yield the maximum SQS mass with 2 M ⊙ . However, when C I and/or C are negative, and baryon number density is sufficiently large, further decreasing C I will lead to a negative quark mass and the requirement of maximum mass with 2 M ⊙ will be no longer satisfied. Therefore, the blue lines in Fig. 1 Fig. 1 will end at a certain point that is indicated by a star. The red line indicates the upper limit of Λ 1.4 =580 measured for GW170817. Moreover, it is worth mentioning that, if we decrease the values of the maximum mass M max or tidal deformability Λ 1.4 , the lines representing them will shift upward, as indicated in panel (a) with up arrows. Comparing the three panels in Fig. 1 Fig. 1 , it is evident that, by increasing the isospin dependence parameter, C I , the dashed line will move significantly downward, whereas the solid line will be barely affected. Consequently, this will widen the SQS stable region. However, comparing panels (b) and (c), it can be observed that, for sufficiently high isospin parameter C I values (e.g., for C I =400 MeV ⋅ fm 3 ), further increasing it does not affect the absolute-stability region. In particular, we verified that beyond <math id="M51"><mrow><msub><mi>C</mi><mtext>I</mtext></msub><mo>≳</mo><mn>10000</mn></mrow></math> MeV ⋅ fm 3 , the absolute-stability region barely changed as the value of the parameter was further increased. By contrast, the blue lines for massive SQSs with 2 M ⊙ and red lines for the SQS tidal deformability with Λ 1.4 =580 tend slightly downward as C I increases, but neither is affected significantly by varying C I . To satisfy these two astrophysical observations (i.e., M max =2 M ⊙ and Λ 1.4 =580), the cross-points of the blue and red lines should be located in the absolute-stability region. Unfortunately, from these three panels, we can easily observe that, although the cross-point approaches the stable region with increasing C I , it does not end up in the stable region, even for C I values as high as 10000 MeV ⋅ fm 3 . Given that the absolute-stability region does not change significantly by increasing C I when it is already sufficiently large, we conclude that in the equivparticle model with mass scaling as given by Eq. (2), the two observables M max =2 M ⊙ and Λ 1.4 =580 can not be simultaneously satisfied by one set of model parameters. However, if we loosen the constraint on tidal deformability to its lower limit, namely decreasing Λ 1.4 to 70, the red line moves upward to the position of the red line with circles, as shown in panel (c). Therefore, it can be inferred that as Λ 1.4 decreases, the intersection point of the blue and red lines moves upward along the blue line. Therefore, we can choose the point ( C,D 1/2 /MeV)=(-0.6,164.5) indicated by the black dot in panel (c) as one of the typical parameter sets, which is located in the SQM absolute-stability region. For comparison, we select parameter sets ( C,D 1/2 /MeV)=(-0.7,168) and (0.08,139) as typical parameter sets for investigating the SQM properties and the SQS structures, which are respectively indicated with a red square and blue diamond in panel (c). To study the effect of isospins on the SQM properties and SQS structures, we also chose points ( C,D 1/2 /MeV)=(-0.6,164.5) in panels (a) and (b) as typical model parameter sets indicated by black dots, which are in the SQM unstable region. In Fig. 2 Fig. 2 , the density behavior of the energy per baryon is shown for typical parameters. The figure reveals that in the equivparticle model, the minimum energy (circles) for each curve is at exactly the same density as zero pressure (asterisks), satisfying the requirement of thermodynamic consistency. Moreover, the one-gluon exchange interaction ( C < 0) can significantly lower the minimum energy per baryon, which makes SQM much more stable than normal nuclear matter. In addition, comparing the three black lines, we find that with increasing C I , the minimum energy per baryon increases. However, for high values of C I , the difference between the solid black line ( C I =10000 MeV ⋅ fm 3 ) and dashed black line ( C I =400 MeV ⋅ fm 3 ) becomes small. This is because, with increasing C I , the difference between u and d quark fractions becomes small in β -equilibrated SQM, which results in vanishing isospin asymmetry δ . We can examine this point in Fig. 3 Fig. 3 . With the parameter set ( C,D 1/2 /MeV)=(-0.6,164.5), the quark fractions are expressed as different C I values indicated in each panel in Fig. 3 Fig. 3 . Evidently, the differences between u , d and s quark fractions become smaller as C I increases. In particular, for C I =10000 MeV ⋅ fm -3 , the difference between the u and d quark fractions is small, which leads to vanishing δ . Additionally, in this case, it is interesting to note that when <math id="M52"><mrow><msub><mi>ρ</mi><mtext>b</mtext></msub><mo>≳</mo><mn>0.6</mn></mrow></math> fm -3 , the u , d , and s quark fractions are essentially equal and approach 0.33, similar to the case of the color-flavor-locked (CFL) phase. The results for the density behavior of quark fractions in this study resemble those reported in Ref. [ 56 56 ]. Figure 4 Figure 4 shows the density behavior of the symmetry energy for typical parameters. From the left panel of Fig. 4 Fig. 4 , it can be seen that the value of C plays an important role in determining the symmetry energy. Generally, for the nuclear matter symmetry energy there are essentially two groups of density behavior, where one increases monotonously with density, and the other first increases and then decreases with increasing density [ 57 57 , 62 62 ]. The latter is typically regarded as soft. Specifically, if the symmetry energy of nuclear matter becomes negative at suprasaturation densities, neutron stars exhibit pure neutron matter cores [ 63 63 ]. Similar to nuclear matter symmetry energy, the behavior of the SQM symmetry energy in the present work can be divided into two groups, depending on the value of the model parameter C . For nonnegative C , the symmetry energy increases with increasing density. Otherwise, it first increases and then decreases with increasing density, which means that it becomes soft at high densities. However, unlike the nuclear matter symmetry energy, the quark matter symmetry energy in our model does not become negative, because the quark masses given by Eq. (2) become negative at sufficiently high densities with negative C . However, the soft quark matter symmetry energy will affect the chemical composition of the SQS core, and may have important implications for the structures and dynamical evolution of SQSs; this issue requires further studies. Moreover, as will be discussed later, to satisfy the requirement of the two-solar-mass constraint of the SQS maximum mass, the symmetry energies shown in the left panel can be as high as 3500 MeV, as indicated by the solid black line. In addition, in the left panel of Fig. 4 Fig. 4 we also show the results from Fig. 1 Fig. 1 in Ref. [ 56 56 ]; these are shown here by the red dashed line indicated with DI-300×30. The reason for choosing the parameter set DI-300 for comparison was that the SQM symmetry energy with DI-300 is of the same order of magnitude as the nuclear matter symmetry given by the RMF model with interaction NLρδ. Therefore, from this panel, we can see that the SQM symmetry energy is approximately 30 times higher than the nuclear matter symmetry energy, which is consistent with the findings in Ref. [ 56 56 ]. In addition, it is interesting to note that at relatively low densities, the symmetry energies shown by the blac k solid and dashed black lines are indistinguishable. The right panel of Fig. 4 Fig. 4 shows that, on one hand, when C I =0, the symmetry energy increases with increasing density, while when C I =400 or 10000 MeV ⋅ fm 3 , the symmetry energy first increases and then decreases for densities above 0.4 fm -3 ; on the other hand, the symmetry energy increases drastically with increasing strength of the isospin dependence parameter C I , e.g., the symmetry energy for C I =10000 MeV ⋅ fm 3 is approximately 20 times that for C I =400 MeV ⋅ fm 3 , or even more than 150 times that for vanishing C I . Left and right panels of Fig. 5 Fig. 5 show, respectively, the SQS structures and tidal deformability. The left panel of Fig. 5 Fig. 5 shows that, as the isospin dependence parameter C I increases, the maximum mass and the corresponding radius of the SQS decrease, which is different from some previously reported results [ 56 56 ]. Moreover, although all maximum mass values in the left panel are above 1.93 M ⊙ , namely, the lower limit of the observed mass of PSR J1614-2230, only the black lines with parametric sets are in the SQM absolute-stability region, which is illustrated in Fig. 1 Fig. 1 . This finding reflects the crucial role of the symmetry energy on the SQM EOS in the equivparticle model. Typically, the black solid line with ( C,D 1/2 /MeV)=(-0.6,164.5) and C I =10000 MeV ⋅ fm 3 yields the maximum mass M max =1.97 M ⊙ , corresponding to the measured mass of PSR J1614-2230. In addition, owing to the fact that with increasing SQS central density the quark equivalent mass given by Eq. (2) becomes negative, the dashed black line where ( C,D 1/2 /MeV)=(-0.7,168) does not reach the maximum mass. However, the mass of the most massive star given by the black dashed line can exceed 2.01 M ⊙ , which is in accordance with the measured mass of PSR J0348 + 0432 [ 28 28 ]. If we constrain the parameters to the SQM absolute-stability region in panel (c) of Fig. 1 Fig. 1 , we obtain that the SQS maximum mass can be as large as 2.02 M ⊙ in the present model, which is still within the limits of the measured mass of PSR J0740 + 6620 with values of 2.08 ± 0.07 M ⊙ [ 29 29 ]. Nevertheless, for GW190814 [ 34 34 ] it was shown that the binary merger possesses a secondary compact object with the mass as high as 2.50-2.67 M ⊙ at the 90% credibility level, which is far beyond the maximum mass of 2.02 M ⊙ in the present model. If the secondary compact object is a neutron star or a quark star, this may imply that other inter-quark interactions should be taken into consideration, or a novel quark mass scaling should be proposed for the equivparticle model. Specifically, Table 1 Table 1 lists the typical quantities and symmetry energy for the SQS core. From Table 1 Table 1 , all the maximum masses shown are approximately 2 M ⊙ , and with a positively valued model parameter C , the SQS maximum mass yields a lower central density. However, it is worth mentioning that the quark matter symmetry energy at the SQS center for positive C is much higher than that for negative C , although the central density is lower. 1 1 ( C,D 1/2 /MeV) 1 1 (0.08,139) 1 1 (-0.6,164.5) 1 1 (-0.7,168) 1 1 M max / M ⊙ 1 1 1.93 1 1 1.97 1 1 2.01 1 1 R (km) 1 1 11.39 1 1 10.44 1 1 10.61 1 1 n 0 (fm -3 ) 1 1 0.98 1 1 1.20 1 1 1.11 1 1 E sym (MeV) 1 1 7493.97 1 1 874.38 1 1 215.93 The right panel in Fig. 5 Fig. 5 shows the relation between the dimensionless tidal deformability and the SQS mass, for the equivparticle model with different parameter sets in Fig. 5 Fig. 5 . The horizontal blue dashed line in the right panel of Fig. 5 Fig. 5 corresponds to the SQS with 1.4 M ⊙ , whereas the two blue solid lines indicate the range of the tidal deformability, that is, 70 ≤ Λ ≤ 580, based on the improved analysis of GW170817 using LIGO and Virgo collaboration [ 38 38 ]. Note that the tidal deformability increases rapidly as the SQS mass decreases, which means that lighter SQSs can be deformed easier than heavier SQSs. Meanwhile, it is evident that all intersections of the typical lines and the blue dashed line are located between the two solid blue lines. Given the maximum masses in the left panel, the typical configurations with the parameter sets we have chosen in the present model can satisfy both the astrophysical observations of the massive compact star PSR J1614-2230 with 1.97 ± 0.04 M ⊙ and the tidal deformability 70 ≤ Λ 1.4 ≤ 580 measured for GW170817. Shown in Fig. 6 Fig. 6 is the dependence of the maximum mass on C I , for the u - d QS (left panel) and SQS (right panel). Comparing the two panels, it is evident that, for the same model parameters, the maximum u - d QS mass is higher than that of the SQS. In addition, unlike the maximum SQS mass, which decreases monotonously with increasing C I , the variation in the maximum mass of the u - d QS is complex. In particular, for nonnegative C the u - d QS maximum mass first increases with C I at C I < 100 MeV ⋅ fm 3 and then decreases with C I ; by contrast, for negative C the maximum u - d QS mass first decreases with increasing C I to approximately 400 MeV ⋅ fm 3 , and then increases with increasing C I . From these two panels, it is evident that for sufficiently large C I (e.g., for C I =10000 MeV ⋅ fm 3 ), further increasing C I will not significantly affect the maximum mass. Furthermore, Fig. 6 Fig. 6 shows that for <math id="M53"><mrow><mo stretchy="false">(</mo><mi>C</mi><mn>,</mn><msqrt><mi>D</mi></msqrt><mo>/</mo><mtext>MeV</mtext><mo stretchy="false">)</mo></mrow></math> =(-0.6,164.5), the maximum u - d QS mass can be as large as 2.3 M ⊙ , exceeding the observed 2.08 M ⊙ for J0740 + 6602. Nevertheless, in this case, the tidal deformability of the u - d QS yields Λ 1.4 ≈ 720, not in the 70 ≤ Λ 1.4 ≤ 580 range. This is illustrated in Table 2 Table 2 , where the typical parameters ( C,D 1/2 /MeV), the maximum mass M max , the tidal deformability Λ 1.4 , and the energy per baryon E/n b are shown for the u - d QS. Moreover, the data in the table suggest that, if the value of <math id="M54"><mrow><msqrt><mi>D</mi></msqrt></mrow></math> is increased to 170, 175, or 180, the tidal deformability will be located in a reasonable range, but the maximum u - d QS mass will decrease as well. Specifically, for <math id="M55"><mrow><mo stretchy="false">(</mo><mi>C</mi><mn>,</mn><msqrt><mi>D</mi></msqrt><mo>/</mo><mtext>MeV</mtext><mo stretchy="false">)</mo></mrow></math> =(-0.6,170), the energy per baryon for the u - d QM will exceed 930 MeV, implying an unstable state for the u - d QM, although the parameter set will yield proper values of <math id="M56"><mrow><msub><mi>M</mi><mrow><mtext>max</mtext></mrow></msub><mo>≈</mo><mn>2.16</mn><msub><mtext>M</mtext><mo>⊙</mo></msub></mrow></math> and Λ 1.4 ≈ 520, satisfying the constraints of astrophysical observations. Therefore, for determining whether the u - d QS can simultaneously satisfy all of the above-mentioned constraints, further detailed investigations are needed. 1 1 ( C,D 1/2 /MeV) 1 1 (-0.6,164.5) 1 1 (-0.6,170) 1 1 (-0.6,175) 1 1 (-0.6,180) 1 1 M max / M ⊙ 1 1 ~2.30 1 1 ~2.16 1 1 ~2.04 1 1 ~1.93 1 1 Λ 1.4 1 1 ~720 1 1 ~520 1 1 ~380 1 1 ~280 1 1 E/n b (MeV) 1 1 ~934 1 1 ~965 1 1 ~993 1 1 ~1021 summary 1 In this study, we investigated the isospin effects on the SQM stability window by introducing an isospin-dependent term into the quark mass scaling to account for the quark-matter symmetry energy, and further studied the EOS of the isospin-asymmetric SQM. We found that by increasing the isospin dependence parameter C I , the SQM absolute-stability region becomes significantly wider, therefore enabling support of the massive compact star PSR J1614-2230 with 1.97 ± 0.04 M ⊙ and the tidal deformability 70 ≤ Λ 1.4 ≤ 580 measured for GW170817. As for quark fractions, we found that with increasing C I the difference between u and d quark fractions became negligible, yielding vanishing isospin asymmetry δ ; in addition, as the baryon number density increases, the u , d , and s quark fractions approach approximately 0.33, similar to the CFL phase results. The symmetry energy properties were then examined by considering different values of C I . It is interesting to note that for C < 0, with increasing C I the SQM symmetry energy first increases and then decreases, implying soft symmetry energy. In addition, the symmetry energy increases drastically with increasing isospin dependence parameter C I , especially for C > 0. Owing to the widening of the SQM absolute-stability region by incorporating the isospin effects into the SQM EOS, the measured mass of PSR J1614-2230 with 1.97 M ⊙ and the dimensionless tidal deformability 70 ≤ Λ 1.4 ≤ 580 measured for GW170817 could be well explained as SQSs by choosing proper model parameter values in the SQM absolute-stability region. However, considering that the measured mass of PSR J0740 + 6620 can be as large as 2.08 M ⊙ , the present equivparticle model still cannot explain more massive SQSs with masses exceeding 2.02 M ⊙ . This may imply that other SQM effects should be considered, e.g., the presence of u - d quark Cooper pairs (2SC phase) [ 51 51 ], or a novel quark mass scaling framework should be developed [ 56 56 ]. Finally, unlike the SQS maximum mass, which decreases monotonously with increasing C I , the variation in the u - d QS maximum mass is uncertain, depending on the sign of the model parameter C . References A.R. Bodmer , Collaped nuclei . Phys. Rev. D 4 , 1601 ( 1971 ). doi: 10.1103/PhysRevD.4.1601 http://doi.org/10.1103/PhysRevD.4.1601 Collaped nuclei 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 H. Terazawa , Super-Hypernuclei in the Quark-Shell Model . J. Phys. Soc. Jpn. 58 , 3555 ( 1989 ). doi: 10.1143/JPSJ.58.3555 http://doi.org/10.1143/JPSJ.58.3555 Super-Hypernuclei in the Quark-Shell Model E. Farhi , R.L. Jaffe , Strange matter . Phys. Rev. D 30 , 2379 ( 1984 ). doi: 10.1103/PhysRevD.30.2379 http://doi.org/10.1103/PhysRevD.30.2379 Strange matter K. Fukushima , T. Hatsuda , The phase diagram of dense QCD . Rept.Prog.Phys. 74 , 014001 ( 2011 ). doi: 10.1088/0034-4885/74/1/014001 http://doi.org/10.1088/0034-4885/74/1/014001 The phase diagram of dense QCD I. Barbour , N.-E. Behilil , E. Dagotto , et al. , Problems with finite density simulations of lattice QCD . Nucl. Phys. B 275 , 296 ( 1986 ). doi: 10.1016/0550-3213(86)90601-2 http://doi.org/10.1016/0550-3213(86)90601-2 Problems with finite density simulations of lattice QCD G. Fowler , S. Raha , R. 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 , 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. Benvenuto , 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.X. Peng , H.C. Chiang , J.J. Yang , et al. , Mass formulas and thermodynamic treatment in the mass-density-dependent model of strange quark matter . Phys. Rev. C 61 , 015201 ( 1999 ). 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 G.X. Peng , H.C. Chiang , B.S. Zou , et al. , 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 R.X. Xu , Solid Quark Stars? . Astrophys. J. Lett. 596 , L59 ( 2003 ). doi: 10.1086/379209 http://doi.org/10.1086/379209 Solid Quark Stars? R.X. Xu , Can cold quark matter be solid? Int . J. Mod. Phys. D 19 , 1437 ( 2010 ). doi: 10.1142/S0218271810017767 http://doi.org/10.1142/S0218271810017767 Can cold quark matter be solid? Int A. Peshier , B. Kmpfer , O. Pavlenko , et al. , 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. Bluhm , B. Kmpfer , 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. Gardim , F. 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 X.J. Wen , Z.Q. Feng , N. Li , et al. , Strange quark matter and strangelets in the quasiparticle model . J. Phys. G 36 , 025011 ( 2009 ). doi: 10.1088/0954-3899/36/2/025011 http://doi.org/10.1088/0954-3899/36/2/025011 Strange quark matter and strangelets in the quasiparticle model T. Xia , L. He , P. Zhuang , Three-flavor Nambu–Jona-Lasinio model at finite isospin chemical potential . Phys. Rev. D 88 , 056013 ( 2013 ). doi: 10.1103/PhysRevD.88.056013 http://doi.org/10.1103/PhysRevD.88.056013 Three-flavor Nambu–Jona-Lasinio model at finite isospin chemical potential J.-F. Xu , Bulk viscosity of interacting magnetized strange quark matter . Nucl. Sci. Tech. 32 , 111 ( 2021 ). doi: 10.1007/s41365-021-00954-3 http://doi.org/10.1007/s41365-021-00954-3 Bulk viscosity of interacting magnetized strange quark matter C. Peng , G.-X. Peng , C.-J. Xia , et al. , Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions . Nucl. Sci. Tech. 27 , 98 ( 2016 ). doi: 10.1007/s41365-016-0095-5 http://doi.org/10.1007/s41365-016-0095-5 Magnetized strange quark matter in the equivparticle model with both confinement and perturbative interactions L. Luo , J. Cao , Y. Yan , et al. , A thermodynamically consistent quasi-particle model without density-dependent infinity of the vacuum zero-point energy . Euro. 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 C.J. Xia , G.X. Peng , S.W. Chen , et al. , Thermodynamic consistency, quark mass scaling, and properties of strange matter . Phys. Rev. D 89 , 105027 ( 2014 ). doi: 10.1103/PhysRevD.89.105027 http://doi.org/10.1103/PhysRevD.89.105027 Thermodynamic consistency, quark mass scaling, and properties of strange matter J.F. Xu , G.X. Peng , F. Liu , et al. , 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 Z.-Y. Lu , G.-X. Peng , S.-P. Zhang , et al. , Quark mass scaling and properties of light-quark matter . Nucl. Sci. Tech. 27 , 148 ( 2016 ). doi: 10.1007/s41365-016-0148-9 http://doi.org/10.1007/s41365-016-0148-9 Quark mass scaling and properties of light-quark matter P. Demorest , T. Pennucci , S. Ransom , et al. , Shapiro delay measurement of a two solar mass neutron star . Nature 467 , 1081 ( 2010 ). doi: 10.1038/nature09466 http://doi.org/10.1038/nature09466 Shapiro delay measurement of a two solar mass neutron star E. Fonseca , T.T. Pennucci , J.A. Ellis , et al. , The NANOGrav nine-year data set: mass and geometric measurements of binary millisecond pulsars . Astrophys. J. 832 , 167 ( 2016 ). doi: 10.3847/0004-637X/832/2/167 http://doi.org/10.3847/0004-637X/832/2/167 The NANOGrav nine-year data set: mass and geometric measurements of binary millisecond pulsars J. Antoniadis , P.C.C. Freire , N. Wex , 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 H.T. Cromartie , E. Fonseca , S.M. Ransom , et al. , Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar . Nature Astronomy 4 , 72 ( 2020 ). doi: 10.1038/s41550-019-0880-2 http://doi.org/10.1038/s41550-019-0880-2 Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar E. Fonseca , H.T. Cromartie , T.T. Pennucci , et al. , Refined mass and geometric measurements of the high-mass PSR J0740+ 6620 . Astrophys. J. Lett. 915 , L12 ( 2021 ). doi: 10.3847/2041-8213/ac03b8 http://doi.org/10.3847/2041-8213/ac03b8 Refined mass and geometric measurements of the high-mass PSR J0740+ 6620 D. Blaschke , M. Cierniak , Studying the onset of deconfinement with multi-messenger astronomy of neutron stars . Astron. Nachr. 342 , 227 ( 2021 ). doi: 10.1002/asna.202113909 http://doi.org/10.1002/asna.202113909 Studying the onset of deconfinement with multi-messenger astronomy of neutron stars H. Liu , J. Xu , P.-C. Chu , Symmetry energy effects on the properties of hybrid stars . Phys. Rev. D 105 , 043015 ( 2022 ). doi: 10.1103/PhysRevD.105.043015 http://doi.org/10.1103/PhysRevD.105.043015 Symmetry energy effects on the properties of hybrid stars B.P. Abbott , R. Abbott , T.D. Abbott , et al. , GW170817: observation of gravitational waves from a binary neutron star inspiral . Phys. Rev. Lett. 119 , 161101 ( 2017 ). doi: 10.1103/PhysRevLett.119.161101 http://doi.org/10.1103/PhysRevLett.119.161101 GW170817: observation of gravitational waves from a binary neutron star inspiral The LIGO Scientific Collaboration and the Virgo Collaboration , GW190814: gravitational waves from the coalescence of a 23 solar mass black hole with a 2.6 solar mass compact object . Astrophys. J. Lett. 896 , L44 ( 2020 ). doi: 10.3847/2041-8213/ab960f http://doi.org/10.3847/2041-8213/ab960f GW190814: gravitational waves from the coalescence of a 23 solar mass black hole with a 2.6 solar mass compact object C. Drischler , S. Han , J.M. Lattimer , et al. , Limiting masses and radii of neutron stars and their implications . Phys. Rev. C 103 , 045808 ( 2021 ). doi: 10.1103/PhysRevC.103.045808 http://doi.org/10.1103/PhysRevC.103.045808 Limiting masses and radii of neutron stars and their implications G.A. Contrera , D. Blaschke , J.P. Carlomagno , et al. , Quark-nuclear hybrid equation of state for neutron stars under modern observational constraints . Phys. Rev. C 105 , 045808 ( 2022 ). doi: 10.1103/PhysRevC.105.045808 http://doi.org/10.1103/PhysRevC.105.045808 Quark-nuclear hybrid equation of state for neutron stars under modern observational constraints M. Cierniak , D. Blaschke , Hybrid neutron stars in the mass-radius diagram . Astrono. Nachr. 342 , 819 ( 2021 ). doi: 10.1002/asna.202114000 http://doi.org/10.1002/asna.202114000 Hybrid neutron stars in the mass-radius diagram B.P. Abbott , R. Abbott , T.D. Abbott , et al. (The LIGO Scientific Collaboration and the Virgo Collaboration), GW170817: Measurements of Neutron Star Radii and Equation of State . Phys. Rev. Lett. 121 , 161101 ( 2018 ). doi: 10.1103/PhysRevLett.121.161101 http://doi.org/10.1103/PhysRevLett.121.161101 (The LIGO Scientific Collaboration and the Virgo Collaboration), GW170817: Measurements of Neutron Star Radii and Equation of State P.G. Krastev , B.-A. Li , Imprints of the nuclear symmetry energy on the tidal deformability of neutron stars . J. Phys. G 46 , 074001 ( 2019 ). doi: 10.1088/1361-6471/ab1a7a http://doi.org/10.1088/1361-6471/ab1a7a Imprints of the nuclear symmetry energy on the tidal deformability of neutron stars Z. Carson , A.W. Steiner , K. Yagi , Constraining nuclear matter parameters with GW170817 . Phys. Rev. D 99 , 043010 ( 2019 ). doi: 10.1103/PhysRevD.99.043010 http://doi.org/10.1103/PhysRevD.99.043010 Constraining nuclear matter parameters with GW170817 N.-B. Zhang , B.-A. Li , Extracting nuclear symmetry energies at high densities from observations of neutron stars and gravitational waves . Eur. Phys. J. A 55 , 39 ( 2019 ). doi: 10.1140/epja/i2019-12700-0 http://doi.org/10.1140/epja/i2019-12700-0 Extracting nuclear symmetry energies at high densities from observations of neutron stars and gravitational waves S. Thakur , V. Thakur , R. Kumar , et al. , Structural properties of rotating hybrid compact stars with color-flavor-locked quark matter core and their tidal deformability . Eur. Phys. J. A 58 , 93 ( 2022 ). doi: 10.1140/epja/s10050-022-00744-4 http://doi.org/10.1140/epja/s10050-022-00744-4 Structural properties of rotating hybrid compact stars with color-flavor-locked quark matter core and their tidal deformability L. Liu , M. Kang , The properties of hybrid stars in a low symmetry energy model . New Astron. 95 , 101817 ( 2022 ). doi: 10.1016/j.newast.2022.101817 http://doi.org/10.1016/j.newast.2022.101817 The properties of hybrid stars in a low symmetry energy model H.-M. Jin , C.-J. Xia , T.-T. Sun , et al. , Deconfinement phase transition and quark condensate in compact stars . ( 2022 ). arXiv:2205.04801 . Deconfinement phase transition and quark condensate in compact stars W.-L. Yuan , A. Li , Z. Miao , et al. , Interacting ud and uds quark matter at finite densities and quark stars . Phys. Rev. D 105 , 123004 ( 2022 ). doi: 10.1103/PhysRevD.105.123004 http://doi.org/10.1103/PhysRevD.105.123004 Interacting ud and uds quark matter at finite densities and quark stars Z.-Q. Miao , C.-J. Xia , X.-Y. Lai , et al. , A bag model of matter condensed by the strong interaction . Int. J. Mod. Phys. E 31 , 2250037 ( 2022 ). doi: 10.1142/S0218301322500379 http://doi.org/10.1142/S0218301322500379 A bag model of matter condensed by the strong interaction B. Holdom , J. Ren , C. Zhang , Quark matter may not be strange . Phys. Rev. Lett. 120 , 222001 ( 2018 ). doi: 10.1103/PhysRevLett.120.222001 http://doi.org/10.1103/PhysRevLett.120.222001 Quark matter may not be strange L. Wang , J. Hu , C.-J. Xia , et al. , Stable up-down quark matter nuggets, quark star crusts, and a new family of white dwarfs . Galaxies 9 , 70 ( 2021 ). doi: 10.3390/galaxies9040070 http://doi.org/10.3390/galaxies9040070 Stable up-down quark matter nuggets, quark star crusts, and a new family of white dwarfs G. Pagliara , J. Schaffner-Bielich , Hadron-quark phase transition at nonzero isospin density: The effect of quark pairing . Phys. Rev. D 81 , 094024 ( 2010 ). doi: 10.1103/PhysRevD.81.094024 http://doi.org/10.1103/PhysRevD.81.094024 Hadron-quark phase transition at nonzero isospin density: The effect of quark pairing G.Y. Shao , M. Colonna , M. Di Toro , et al. , Influence of vector interactions on the hadron-quark/gluon phase transition . Phys. Rev. D 85 , 114017 ( 2012 ). doi: 10.1103/PhysRevD.85.114017 http://doi.org/10.1103/PhysRevD.85.114017 Influence of vector interactions on the hadron-quark/gluon phase transition K.S. Jeong , S.H. Lee , Symmetry energy in cold dense matter . Nucl. Phys. A 945 , 21 ( 2016 ). doi: 10.1016/j.nuclphysa.2015.09.010 http://doi.org/10.1016/j.nuclphysa.2015.09.010 Symmetry energy in cold dense matter L.-W. Chen , Nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula . Phys. Rev. C 83 , 044308 ( 2011 ). doi: 10.1103/PhysRevC.83.044308 http://doi.org/10.1103/PhysRevC.83.044308 Nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula A.W. Steiner , M. Prakash , J.M. Lattimer , et al. , Isospin asymmetry in nuclei and neutron stars . Phys. Rept. 411 , 325 ( 2005 ). doi: 10.1016/j.physrep.2005.02.004 http://doi.org/10.1016/j.physrep.2005.02.004 Isospin asymmetry in nuclei and neutron stars M.Di Toro , B. Liu , V. Greco , et al. , Symmetry energy effects on the mixed hadron-quark phase at high baryon density . Phys. Rev. C 83 , 014911 ( 2011 ). doi: 10.1103/PhysRevC.83.014911 http://doi.org/10.1103/PhysRevC.83.014911 Symmetry energy effects on the mixed hadron-quark phase at high baryon density S.-W. Chen , L. Gao , G.-X. Peng , One-gluon-exchange effect on the properties of quark matter and strange stars . Chin. Phys. C 36 , 947 ( 2012 ). doi: 10.1088/1674-1137/36/10/005 http://doi.org/10.1088/1674-1137/36/10/005 One-gluon-exchange effect on the properties of quark matter and strange stars P.-C. Chu , L.-W. Chen , Quark matter symmetry energy and quark stars . Astrophys. J. 780 , 135 ( 2014 ). doi: 10.1088/0004-637X/780/2/135 http://doi.org/10.1088/0004-637X/780/2/135 Quark matter symmetry energy and quark stars L.-W. Chen , Symmetry Energy in Nucleon and Quark Matter . Nucl. Phys. Rev. 34 , 20 ( 2017 ). doi: 10.11804/NuclPhysRev.34.01.020 http://doi.org/10.11804/NuclPhysRev.34.01.020 Symmetry Energy in Nucleon and Quark Matter T. Damour , A. Nagar , Relativistic tidal properties of neutron stars . Phys. Rev. D 80 , 084035 ( 2009 ). doi: 10.1103/PhysRevD.80.084035 http://doi.org/10.1103/PhysRevD.80.084035 Relativistic tidal properties of neutron stars T. Hinderer , B.D. Lackey , R.N. Lang , et al. , Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral . Phys. Rev. D 81 , 123016 ( 2010 ). doi: 10.1103/PhysRevD.81.123016 http://doi.org/10.1103/PhysRevD.81.123016 Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral S. Postnikov , M. Prakash , J.M. Lattimer , Tidal Love numbers of neutron and self-bound quark stars . Phys. Rev. D 82 , 024016 ( 2010 ). doi: 10.1103/PhysRevD.82.024016 http://doi.org/10.1103/PhysRevD.82.024016 Tidal Love numbers of neutron and self-bound quark stars E.-P. Zhou , X. Zhou , A. Li , Constraints on interquark interaction parameters with GW170817 in a binary strange star scenario . Phys. Rev. D 97 , 083015 ( 2018 ). doi: 10.1103/PhysRevD.97.083015 http://doi.org/10.1103/PhysRevD.97.083015 Constraints on interquark interaction parameters with GW170817 in a binary strange star scenario L.-W. Chen , C.M. Ko , B.-A. Li , Determination of the stiffness of the nuclear symmetry energy from isospin diffusion . Phys. Rev. Lett. 94 , 032701 ( 2005 ). doi: 10.1103/PhysRevLett.94.032701 http://doi.org/10.1103/PhysRevLett.94.032701 Determination of the stiffness of the nuclear symmetry energy from isospin diffusion Y. Zhou , L.-W. Chen , Ruling out the supersoft high-density symmetry energy from the discovery of PSR J0740+6620 with mass <math id="M57"><mrow><msubsup><mrow><mn>2.14</mn></mrow><mrow><mo>−</mo><mn>0.09</mn></mrow><mrow><mo>+</mo><mn>0.10</mn></mrow></msubsup><msub><mi>M</mi><mo>⊙</mo></msub></mrow></math> . Astrophys. J. 886 , 52 ( 2019 ). doi: 10.3847/1538-4357/ab4adf http://doi.org/10.3847/1538-4357/ab4adf Ruling out the supersoft high-density symmetry energy from the discovery of PSR J0740+6620 with mass 2.14−0.09+0.10M⊙ All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Jian-Feng Xu, Cheng-Jun Xia,Zhen-Yan Lu,Guang-Xiong Peng and Ya-Peng Zhao. The first draft of the manuscript was written by Jian-Feng Xu and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. 10.1007/s41365-022-01130-x This work was supported by the National Natural Science Foundation of China (Nos. 12005005 and 11875052), the National SKA Program of China (No. 2020SKA0120300), the Hunan Provincial Nature Science Foundation of China (No. 2021JJ40188), and the Scientific Research Start-up Fund of Talent Introduction of Suqian University (No. Xiao2022XRC061). 2022-07-28 2022-09-23 2022-09-26 2022-11-11 2022-12-16 © 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. 2022 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-022-01130-x http://dx.doi.org/10.1007/s41365-022-01130-x . Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law. 2022 Nuclear Science and Techniques 11 33