Your new experience awaits. Try the new design now and help us make it even better

BRIEF RESEARCH REPORT article

Front. Astron. Space Sci., 22 July 2025

Sec. Nuclear Physics​

Volume 12 - 2025 | https://doi.org/10.3389/fspas.2025.1600563

This article is part of the Research TopicStrong and Weak Interactions in Compact StarsView all 7 articles

Confronting recent light compact star observations with color-flavor locked quark matter

  • 1Theoretical Astrophysics, Eberhard-Karls Universität Tübingen, Tübingen, Germany
  • 2Department of Theoretical Physics, Aristotle University of Thessaloniki, Thessaloniki, Greece

Recent analyses on the properties of the central compact object in the HESS J1731-347 remnant and the PSR J1231-1411 pulsar indicated that these two compact objects are characterized by similar (low) masses and possibly different radii. This paper aims at reconciling the aforementioned measurements by utilizing the widely employed color-flavor locked (CFL) MIT bag model. The main objective is related to the examination of the acceptable values for the color superconducting gap Δ and the bag parameter B. Furthermore, our analysis involves two distinct hypotheses for the nature of compact stars. Firstly, we considered the case of absolute stability for strange quark matter and we found that it is possible to explain both measurements, while also respecting the latest astronomical constraints on the masses and radii of compact stars. Secondly, we studied the case of hybrid stellar matter (transition from hadrons to quarks), and concluded that, when early phase transitions are considered, the simultaneous reconciliation of both measurements leads to results that are inconsistent to the existence of massive compact stars. However, we showed that all current constraints may be satisfied under the consideration that the HESS J1731-347 remnant contains a slow stable hybrid star.

1 Introduction

One of the most important unresolved questions in theoretical astrophysics is related to the nature of matter in the cores of compact stars. Notably, compact stars could be composed solely by hadrons (nucleons and hyperons), but the extreme conditions that prevail in their interior may allow for the presence of exotic forms of matter, such as deconfined quarks (Witten, 1984; Annala et al., 2020). The latter opens up intriguing scenarios, such as the existence of strange stars (Weber, 2005), composed purely of quark matter, or hybrid stars (Heiselberg and Hjorth-Jensen, 2000), where a quark core is surrounded by a layer of hadrons. Interestingly, given that different hypotheses for the composition of stellar matter may predict distinct properties for the structure of compact stars (Lattimer and Prakash, 2001; Glendenning and Kettner, 2000), precisely inferred measurements on masses and radii of compact objects may provide important insight into their nature.

In 2022, the analysis of Doroshenko et al. (2022) provided puzzling values for the mass (M) and radius (R) of the central compact object (CCO) in the HESS J1731-347 remnant. More precisely, the authors reported that, at the 1σ level, M=0.770.17+0.20M and R=10.40.78+0.86 km. Notably, the surprisingly low mass and radius values, provided by Doroshenko et al. (2022), led to a wide range of studies attempting their explanation, by considering several different hypotheses for the nature of ultra-dense matter (Di Clemente et al., 2024; Horvath et al., 2023; Oikonomou and Moustakidis, 2023; Das and Lopes, 2023; Rather et al., 2023; Tsaloukidis et al., 2023; Sagun et al., 2023; Laskos-Patkos et al., 2024; Laskos-Patkos et al., 2025; Li et al., 2024; Mariani et al., 2024; Brodie and Haber, 2023; Huang et al., 2024; Li and Sedrakian, 2023; Kubis et al., 2023; Routaray et al., 2024; Char and Biswas, 2024; Tewari et al., 2024). Specifically, it has been shown that the aforementioned CCO could be a strange quark star (Di Clemente et al., 2024; Horvath et al., 2023; Oikonomou and Moustakidis, 2023; Das and Lopes, 2023; Rather et al., 2023; Yang and Pi, 2024; Gholami et al., 2024) (given that there is a mechanism that suppresses rapid cooling, such as color superconductivity (Di Clemente et al., 2024; Horvath et al., 2023)), a hybrid star characterized by an early phase transition (Tsaloukidis et al., 2023; Sagun et al., 2023; Laskos-Patkos et al., 2024; Laskos-Patkos et al., 2025; Li et al., 2024; Mariani et al., 2024; Gao et al., 2024; Pal et al., 2025), or a neutron star described by a soft nuclear EOS (at low densities) (Brodie and Haber, 2023; Huang et al., 2024; Li and Sedrakian, 2023; Kubis et al., 2023). For a relevant critique on the estimation of Doroshenko et al. (2022) and the corresponding required assumptions we refer to the work of Alford and Halpern (2023).

The recent analysis of Salmi et al. (2024a) also provided intriguing values for the properties of the pulsar PSR J1231-1411. Notably, their results appeared to be sensitive to the selection of the radius prior they used. By limiting the radius to be consistent with previous observational constraints and nuclear theory, the authors indicated that M=1.040.03+0.05M and R=12.60.3+0.3 km. In addition, by employing an uninformative prior, and constraining the radius between 10 and 14 km, they reported a similar value for M and a higher value for R. Interestingly, the recently published work of Qi et al. (2025) provided an alternative estimate on the properties of PSR J1231-1411. In particular, by implementing an algorithm based on the spherical star Schwarzschild spacetime and the Doppler approximation, the authors inferred that the mass-radius properties of the aforementioned pulsar are M=1.120.07+0.07M and R=9.910.86+0.88 km. In any case, future studies may be required to obtain a more accurate view of this compact object. In the present work, to enrich and extend our investigation, we will consider both estimates of the properties of PSR J1231-1411, referring to the results of Salmi et al. (2024a) as Case I, and to the results of Qi et al. (2025) as Case II.

Given that color-flavor locked (CFL) quark matter (Alford et al., 1999; Alford, 2001; Alford et al., 2008) has been successfully employed in the reconciliation of HESS J1731-347, not only for its mass and radius but also for its thermal evolution (Di Clemente et al., 2024; Horvath et al., 2023), we aim to consider it in order to examine the possible simultaneous explanation of the PSR J1231-1411 properties. To do so, we consider two distinct hypotheses for strange quark matter (SQM): a) absolutely stable, b) energetically favored at high baryon density. Of utmost importance is to examine if the resulting EOSs satisfy other precisely inferred mass-radius measurements (NICER mission) (Choudhury et al., 2024; Salmi et al., 2024b; Vinciguerra et al., 2024; Riley et al., 2019; Miller et al., 2019) and also allow for the existence of stable massive stars beyond 2M (Antoniadis et al., 2013; Cromartie et al., 2020; Romani et al., 2022). In this way, we aim to provide a general view on the possible existence of CFL quark matter in compact stars by utilizing all available and accurate observational constraints from rotation-powered millisecond pulsars. The main motivation for performing the present work is analyzed in more depth in the two following paragraphs.

When considering Case I (Salmi et al., 2024a), it is rather interesting that, while both the CCO in the HESS J1731-347 remant and the PSR J1231-1411 pulsar have similar masses, their radii do not overlap at the 1σ level. In addition, the fact that PSR J1231-1411 has a slightly higher mass and radius indicates that the M(R) curve may exhibit positive slope (dM/dR>0). Notably, this is a standard feature when considering strange stars and hence we aim to constrain the phenomenological parameters (bag constant and color superconducting gap) of the CFL MIT bag model (Alford et al., 2001; Lugones and Horvath, 2002; Alford et al., 2005a), assuming the absolute stability of strange matter, in light of the aforementioned measurements (Doroshenko et al., 2022; Salmi et al., 2024a; Choudhury et al., 2024; Salmi et al., 2024b; Vinciguerra et al., 2024; Riley et al., 2019; Miller et al., 2019; Antoniadis et al., 2013; Cromartie et al., 2020; Romani et al., 2022). In any case, this characteristic property (dM/dR>0) may not be essential if one considers the possible existence of hybrid stars. In principle, hybrid EOSs, characterized by a large density discontinuity, enable the existence of stars with comparable masses and significantly different radii (Benic et al., 2015; Alford and Sedrakian, 2017; Alvarez-Castillo et al., 2019; Blaschke et al., 2020; Christian and Schaffner-Bielich, 2021; Christian and Schaffner-Bielich, 2022; Li et al., 2020; Li et al., 2023a; Li et al., 2023b; Naseri et al., 2024; Jiménez et al., 2024; Li et al., 2025; Zhang and Li, 2025). Therefore, a detailed analysis on the reconciliation of both HESS J1731-347 and PSR J1231-1411 may provide important insight on the existence of low mass twin star solutions (within the CFL MIT bag model).

In Case II (Qi et al., 2025), the radius of PSR J1231-1411 has an upper bound (at 1σ) that is even lower compared to the one of the CCO in HESS J1731-347. Taking into consideration that the former pulsar is also more massive than the latter, one can deduce that the simultaneous explanation of both measurements may require a sufficiently soft EOS. In this scenario, it is particularly interesting to examine if the CFL model can also fulfill the conservative 2M maximum mass constraint or support the potential existence of even more massive compact stars like PSR J0952-0607 (Romani et al., 2022).

This paper is structured as follows: Section 2 sets the theoretical framework of the present study. Section 3 presents the results for the two distinct hypotheses on the nature of compact stars (quark or hybrid stars) along with a detailed discussion on our findings. Lastly, Section 4 highlights the main insights and conclusions derived from this research.

2 Color-flavor locked equation of state

The equation of state for CFL quark matter can be formulated within the MIT bag model framework. To the order of Δ2 and ms2, where ms represents the mass of the strange quark and μ the quark chemical potential, the pressure and energy density can be expressed as follows (=c=1) (Lugones and Horvath, 2002; Alford et al., 2005a)

P=3μ44π2a4+9αμ22π2B(1)
ε=9μ44π2a4+9αμ22π2+B,(2)

where a4 is a parameter that mimics the impact of perturbative QCD (pQCD) corrections (Alford et al., 2005a), B represents the bag constant and α=ms26+2Δ23. Notably, an analytical expression for ε(P) can be obtained by combining Equations 1, 2:

ε=3P+4B9αμ2π2, where μ2=9α+12π2a4B+P+81α21/23a4(3)

Equation 3 reduces to the one presented in Lugones and Horvath (2002) for a4=1, which indicates no strong interactions. For CFL quark matter to be absolutely stable, its energy per baryon must be lower than the one in the most stable nucleus, 56Fe(mb(Fe)=930MeV), at zero pressure and temperature (Farhi and Jaffe, 1984). This condition establishes an upper limit for the bag constant as a function of ms and Δ, for a fixed value of a4:

B<ms2mb212π2+Δ2mb23π2+mb4108π2a4(4)

Conversely, a lower limit can be determined by requiring that two-flavor quark matter should be less stable than nuclear matter. In the context of the MIT bag model (a4=1), this condition is expressed as (Farhi and Jaffe, 1984):

B57MeV · fm3(5)

This value decreases when considering lower values for a4. In this study, we use the constraint derived with a4=1 (Equation 5), ensuring that all the CFL EOSs used characterize absolutely stable quark matter.

Finally, it is worth noting that in order to ensure that the CFL phase represents the favorable state of matter (compared to 2SC, unpaired or gapless CFL matter) the following condition needs to be met (Alford et al., 2005b):

Δms2/2μ.(6)

Notably, Equation 6 holds for all of the quark EOSs constructed in the present study.

3 Results and discussion

3.1 Strange stars

In this section, we present our results derived under the assumption that SQM represents the true ground state of matter. However, a key question, that requires a proper discussion, immediately arises: if SQM is more stable than nuclear matter why does normal matter persist? One explanation may be that nuclear matter is metastable [see Horvath et al. (1992), Olesen and Madsen (1994), Iida and Sato (1998), Bombaci et al. (2008), Bombaci et al. (2009), Ren and Zhang (2020) and references therein], separated from the favorable SQM state by a significant energy barrier. Thus, at low densities, the conversion to SQM may be suppressed due to quantum tunneling limitations. However, the extreme neutron star environment, characterized by high densities, may enhance the possibility of a SQM droplet appearing via quantum fluctuations. In addition, extreme events such as supernovae and neutron star mergers may also facilitate SQM formation. Given that the appearance of a SQM seed could trigger the conversion of the entire hadronic star into a strange star, it has long been hypothesized that all compact stars may be of SQM nature.

Proceeding with our calculations, we solved the system of the TOV equations (Oppenheimer and Volkoff, 1939), using a wide range of B(57300MeVfm3) and Δ(50300MeV) values to identify the MR curves that are compatible with the light neutron star measurements for HESS J1731-347 and PSR J1231-1411 (at the 1σ level). The latter was examined utilizing the results of both Salmi et al. (2024a) (Case I) and Qi et al. (2025) (Case II). Of course, we only considered B and Δ pairs that lie within the stability window defined by Equation 4 (fixing the strange quark mass at 95 MeV), while to investigate the possible effects of the a4 parameter we employed two distinct values, a4=1 and a4=0.7. The latter value was selected as it results to EOSs similar to the quark matter EOS from pQCD [e.g., see Alford et al. (2005a), Fraga et al. (2001)].

Notably, in the present study we work under the assumption of a universal EOS (a single EOS that accounts for all observations). Therefore, the derived EOSs should be compatible to all state-of-the-art multimessenger constraints on the mass and radius of compact stars. Thus, apart from HESS J1731-347 and PSR J1231-1411, we also utilized the corresponding data related to PSR J0952-0607 (Romani et al., 2022), PSR J0030+0451 (Miller et al., 2019) and PSR J0437-4715 (Choudhury et al., 2024). Lastly, to highlight the peculiar nature of the recent XTE J1814-338 measurement (Kini et al., 2024), we have included it in our analysis, although we did not attempt to interpret it simultaneously with the aforementioned constraints [for some recent works on its reconciliation see Zhou and Huang (2025), Yang et al. (2025), Laskos-Patkos and Moustakidis (2025), Lopes and Issifu (2025), M. Veselský et al. (2025)].

Figure 1 displays the (B,Δ) pairs that are compatible with the above measurements for two different values of a4, accounting for both estimations of PSR J1231-1411’s properties. The red, purple, blue, yellow and gray parameter spaces correspond to XTE J1814-338, PSR J1231-1411 (Case I), PSR J0952-0607, HESS J1731-347 and PSR J0030+0451 measurements, respectively. The last three areas extend all the way to the horizontal axis, though it is not visible in the graphs. The region between the black solid lines corresponds to PSR J0437-4715, while the one between the purple dashed lines to PSR J1231-1411 (Case II). Figure 1a illustrates the parameter spaces of the above objects for a4=1 considering Case I, while Figure 1b shows the same graph for a4=0.7. As is evident, although a smaller a4 value alters the CFL EOS, resulting in a modified stability window, it has only minimal impact on the parameter spaces that are compatible to different observations. Finally, in Figures 1c,d, Case II is examined for a4=1 and a4=0.7 respectively.

Figure 1
www.frontiersin.org

Figure 1. BΔ pairs that generate MR curves that are compatible (at the 1σ level) with the neutron star observations PSR J0952-0607 (Romani et al., 2022) (blue region), HESS J1731-347 (Doroshenko et al., 2022) (yellow region), XTE J1814-338 (Kini et al., 2024) (red region), PSR J1231-1411 (Salmi et al., 2024a) Case I (purple region), PSR J0030+0451 (Miller et al., 2019) (gray region), PSR J0437-4715 (Choudhury et al., 2024) (area between black solid lines) and PSR J1231-1411 (Qi et al., 2025) Case II (area between purple dashed lines). (a) illustrates the parameter spaces when considering Case I and for a4=1, while (b) is identical for a4=0.7. (c,d) examine Case II for a4=1 and a4=0.7 respectively. The dark gray area represents the region where CFL quark matter is not absolutely stable.

Figures 1a,b indicate that HESS J1731-347 and PSR J1231-1411 (Case I), can both be reconciled within the range of parameters that are compatible to the latter, while also respecting the maximum-mass constraints imposed by PSR J0952-0607. This region is characterized by small values for the bag constant, close to the minimum possible one, and large values for the color superconducting gap. As a result, the CFL EOS is extremely stiff, leading to very high maximum masses for all (B,Δ) pairs (above 3M). The latter is demonstrated in Figure 2b, where the maximum masses for all (B,Δ) combinations of the purple region of Figure 1b are calculated (for the a4=0.7 case).

Figure 2
www.frontiersin.org

Figure 2. (a) Mass-radius curves that are compatible with PSR J1231-1411 Case I (solid lines) and Case II (dashed lines), with a4=0.7. The BΔ pairs utilized are representative of the purple region (Case I) and the region between the purple dashed lines (Case II), see Figures 1b,d. The numbers above each curve represent the (B,Δ) pairs in (MeVfm3,MeV). (b) Close-up of PSR J1231-1411’s (Salmi et al., 2024a) parameter space considering Case I with a4=0.7, plotted together with the maximum masses reached in this window. Notably, a minimum maximum mass of 3M is needed. (c) Parameter space that is compatible with PSR J0952-0607 (Romani et al., 2022), PSR J0030+0451 (Miller et al., 2019), HESS J1731-347 (Doroshenko et al., 2022), PSR J0437-4715 (Choudhury et al., 2024) and PSR J1231-1411 (Qi et al., 2025) (Case II) measurements, including the maximum masses reached. This area corresponds to the region between the upper limit of PSR J0030+0451 and the lower one of PSR J1231-1411, see Figure 1d.

Another interesting point from Figures 1a,b is related to the fact that the parameter spaces which are compatible with PSR J0437-4715 and PSR J1231-1411 (Case I) do not overlap at the 1σ level, implying that they cannot be simultaneously explained within this model. However, this tension is mild, as we have checked that a narrow common parameter space exists, when considering the 2σ estimation for PSR J0437-4715 (see Figure 2a, where we have plotted the MR diagrams for representative EOSs which are compatible with PSR J1231-1411).

Considering Case II, Figures 1c,d show that there is a small range of parameters which is compatible with all measurements at the 1σ level (of course except XTE J1814-338). These areas are defined by the lower boundary from PSR J1231-1411 (purple dashed line) and the upper limit from PSR J0030+0451 (gray region). Notably, for PSR J0030+0451 we employed the measurement of Miller et al. (2019). The aforementioned common areas do not change significantly when adopting the measurements of Riley et al. (2019), and they would increase when using the lower (regarding the radius) of the three estimations provided by Vinciguerra et al. (2024).

The parameter space in which all measurements can be explained by pure CFL matter (Case II) is shown in Figure 2c, alongside the corresponding maximum mass predictions (for the a4=0.7 case). The resulting maximum masses span the range of approximately 2.253.25M, which is consistent with current observations of heavy pulsars. This maximum mass range, along with that obtained for the parameter space compatible with PSR J1231-1411 (Case I) for a4=0.7, points toward the possible existence of quark stars with masses exceeding 3M. Notably, such masses reside within the so-called lower mass gap (Shao, 2022; Samsing and Hotokezaka, 2021), in which the nature of compact objects remains undetermined. Interestingly, evidence for compact objects within the mass gap has emerged from multiple recent studies, based on observations of non-interacting binary systems (Thompson et al., 2019; Jayasinghe et al., 2021), radio pulsar surveys (Barr et al., 2024), and gravitational wave detections from binary mergers (Abbott R. et al., 2020; Abbot et al., 2023; Abac et al., 2024; Abbot et al., 2017; Abbot B. P. et al., 2020).

A final remark needs to be made with regards to the possible explanation of the thermal evolution for the CCO in the HESS J1731-347 remnant. Notably, the CCO has a rather high temperature for its estimated age which suggests slow cooling (Di Clemente et al., 2024), analogous to those of purely hadronic stars. Previous works (Di Clemente et al., 2024; Horvath et al., 2023) have qualitatively suggested that superconductivity may suppress the rapid cooling processes that are expected in unpaired SQM and therefore the temperature of the CCO could be explained. Interestingly, according to Horvath et al. (2023) this could occur only under the consideration of vanishingly small gap values Δ<1 MeV. This is related to the fact that for large values of pairing gaps the suppression of standard cooling processes in the stellar core is dramatic. Therefore, there is a potential tension between the necessity for large gaps, dictated by the large maximum mass constraints, and the low gaps suggested by qualitative analyses on the cooling of the CCO in HESS J1731-347. In addition, it remains an open question whether the low central density expected in the low mass CCO is sufficient for CFL matter to appear for vanishingly small pairing gaps [considering the conditions that need to be met for CFL matter to be energetically favorable (Alford et al., 2001; Alford et al., 2005b)]. In any case, a future precise analysis on the cooling of CFL stars incorporating all of the aforementioned issues, including also potential uncertainties on the estimated temperature and age, is essential to safely conclude about whether the CCO in the HESS J1731-347 could be in fact a CFL strange star.

3.2 Hybrid stars

In the previous section, we considered parametrizations that support the absolute stability of CFL quark matter. At this point, we aim to examine the scenario of explaining all current multimessenger constraints by considering a first-order phase transition, from hadronic to CFL quark matter, in the stellar interior.

Notably, the simultaneous reconciliation of both HESS J1731-347 and PSR J1231-1431 (in Case I) measurements would be rather difficult assuming a purely hadronic EOS. In particular, the hadronic model should be rather soft, at low densities, to support the low radius associated with HESS J1731-347 in the sub-solar mass region, and then it should rapidly stiffen to achieve the higher radius of PSR J1231-1431 at slightly larger masses. Interestingly, the low-density domain of the nuclear EOS can be effectively constrained through parity-violating electron scattering experiments, which aim to measure the neutron skin thickness of different nuclei. In particular, the PREX-II collaboration provided a measurement for the neutron skin thickness of lead, which pointed to stiff EOS behavior at low densities (Adhikari et al., 2021; Reed et al., 2021). Then, the subsequent CREX experiment extracted a puzzling value for the neutron skin thickness of calcium (Adhikari et al., 2022), supporting softer models. More precisely, the simultaneous explanation of both CREX and PREX-II was not possible with the use of traditional energy density functionals (Reinhard et al., 2023; Tagami et al., 2022; Miyatsu et al., 2023; Kumar et al., 2023; Burgio et al., 2024). However, recent attempts have proposed sophisticated modifications to the Lagrangians which are used to describe nuclear matter, and they achieved an explanation of both experimental values (Reed et al., 2024; Kumar et al., 2024). Nonetheless, the predicted EOSs are characterized as extremely stiff. Notably, the latter issue was at some extent resolved in Salinas and Piekarewicz (2024).

Considering that current experimental constraints might point to stiff behavior for the hadronic EOS at low density, in the present study, we are going to employ a nuclear model which is sufficiently stiff, so that it crosses the PSR J1231-1411 contour (Case I), and, therefore, potentially incompatible to the HESS J1731-347 constraints. In particular, we are going to use a widely employed Skyrme model, namely Ska (Köhler, 1976; Typel et al., 2022; Gulminelli and Raduta, 2015; Danielewicz and Lee, 2009; Baym et al., 1971). Thus, in our attempt for explaining all current astronomical measurements, by considering a first-order phase transition to CFL quark matter, we will rely on the hybrid branch to account for the existence of the central compact object in the HESS J1731-347 remnant. It is important to comment that, the employed nuclear model is of nucleonic composition. In principle, as density increases, hyperonic degrees of freedom may appear, altering the properties of the EOS. Given the uncertainties related to hyperon-hyperon and hyperon-nucleon interactions, in the present study we work under a simplified framework that neglects their existence.

In Figure 3a, we depicted the mass-radius dependence for hybrid EOSs constructed by varying the values of Δ and B (fixing a4=1). Notably, the pairs of Δ and B are selected so that the phase transition density remains fixed, to allow for the hadronic branch to account for PSR J1231-1411. As one can observe, the hybrid branch crosses the HESS J1731-347 inferred region. However, the EOS softening that is induced due to the phase transition plays a critical role on the determination of the resulting maximum mass. As a consequence, the hybrid branch never reaches the 2M and, therefore, it fails to reproduce one of the most robust constraints derived by compact object observations. In Figure 3b we attempted to resolve this issue by lowering the value of the parameter a4. The motivation to do that relies on the knowledge that reducing a4 reduces the value of pressure for a given baryon chemical potential and, as a consequence, for the same parameter values of B and Δ, shifts the onset of the phase transition to higher densities. Therefore, to recover the desired phase transition onset (i.e., when the hadronic branch crosses the PSR J1231-1411 contour) we need to reduce the value of the bag constant (Baym et al., 2018), leading to EOS stiffening. While the constructed EOSs achieve higher masses (still not the conservative 2M limit), the induced stiffening does not allow for an explanation of the HESS J1731-347 measurement. Finally, we tried to stiffen the quark EOS by adding an ad hoc vector interaction between quarks. We refer to this attempt as ad hoc since the resulting models do not follow from an initial unified description, but are produced by simply adding a vector mean-field term in the CFL MIT bag model expressed in the canonical ensemble [following Kanakis-Pegios et al. (2024), Yang et al. (2021), Yang et al. (2023), the strength of the vector repulsion is quantified via the parameter Gv=(g/mv)2, where g is the coupling constant and mv the mass of the mediating boson]. For more details on the inclusion of this vector term the reader is referred to Kanakis-Pegios et al. (2024). In any case, the results depicted in Figure 3c indicate that such an approach also does not suffice to reconcile all of the considered astronomical constraints. However, while the CFL MIT bag model might fail in reproducing all measurements (when assuming a phase transition), we can not exclude that a different, comprehensive analysis (e.g., via a relativistic mean-field type model), that includes such a stiffening term, may enable the simultaneous explanation of all observations.

Figure 3
www.frontiersin.org

Figure 3. Mass-radius diagrams for hybrid EOSs combining the Ska hadronic EOS and the CFL MIT bag model. The gray contour regions denote mass and radius measurements related to PSR J0740 + 6620 (Salmi et al., 2024b), PSR J0030+0451 (Miller et al., 2019), HESS J1731-347 (Doroshenko et al., 2022), PSR J0437-4715 (Choudhury et al., 2024), PSR J1231-1411 (Salmi et al., 2024a), XTE J1814-338 (Kini et al., 2024) and GW170817 (Abbot et al., 2017). The solid (dashed) contours correspond to 2σ(1σ) confidence. (a) Shows the results for a4=1. (b) Shows the results for a4=0.7. (c) Includes the results for a4=1 including the contribution of a vector interaction among quarks with coupling Gv=0.15fm2.

Up to this moment, we have only considered Case I for PSR J1231-1411. However, it is clear that the consideration of Case II would not alter any conclusion about our inability of reproducing all observational constraints when we combine the CFL model with a stiff hadronic EOS (similar to the employed one). Nevertheless, future refinement of theoretical models or nuclear experiments that point towards a softer nuclear model may alter the current picture.

Interestingly, we can attempt to reconcile all of the considered astronomical constraints by allowing the hadronic branch to reach the two solar masses and then induce an extremely strong phase transition to make the mass-radius diagram drop to cross the HESS J1731-347 contour. Notably, hybrid stars may be stable even at a descending branch of a mass-radius diagram, assuming that the phase conversion is sufficiently slow (Pereira et al., 2018). The characterization slow refers to the magnitude of the phase conversion timescale τconv, compared to the period τosc of radial oscillations that the stellar configurations may undergo. In particular, if τconvτosc, then hybrid stars could be stable even at regions where the well-known turning point criterion would rule them out (Pereira et al., 2018; Rather et al., 2024). The idea of explaining the HESS J1731-347 CCO as a slow stable hybrid star is not new. More precisely, Mariani et al. (2024) established the plausibility of this scenario in their recent work. In the present study, we aim to examine if the CFL MIT bag model can lead to such an explanation. Specifically, the authors of Mariani et al. (2024) considered a constant speed of sound parametrization where there are three free parameters, namely the transition pressure, the speed of sound and the energy density jump. As a consequence, one can construct a sufficiently stiff EOS, by controlling the speed of sound, which would be characterized by a large energy density jump at any chosen transition density. In contrast, in the current framework we have two free parameters (we are keeping the value of a4=0.8 fixed to account for pQCD constraints; see below) and their interplay has to self-consistently provide an appropriate transition density (so that the hadronic branch reaches the 2M), and a large enough energy density discontinuity to trigger an almost vertical drop of the mass-radius curve. Therefore, it is particularly interesting to examine which pairs of the phenomenological parameters B and Δ could lead to a simultaneous explanation of all current astronomical constraints.

In Figure 4a, we depict the mass-radius dependence for hybrid EOSs constructed by varying the color superconducting gap in the range [0,600] MeV. The bag constant was derived by considering that the phase transition occurs when the hadronic mass reaches the 2M. Notably, all of the hybrid branches, in Figure 4a, are plotted until the last slow stable configuration or terminal configuration as named in Mariani et al. (2024). As one can observe, the reconciliation of the HESS J1731-347 constraints is possible, but it would require extremely large values for both B and Δ (see the legend of Figure 4), at least compared to those found usually in the literature. This is, however, somewhat expected as the majority of studies focuses on phase transitions of moderate energy density jumps and low to moderate transition densities. Hence, the extreme conditions required for the explanation of HESS J1731-347 require extreme values for the phenomenological parameters of CFL quark matter. In any case, it is worth noting that recent and robust constraints have been imposed on the color superconduncting gap in Abbott et al. (2025). While a value of around 400 MeV is potentially compatible to their results, larger values may be disfavored [see Figure 2 in Abbott et al. (2025)].

Figure 4
www.frontiersin.org

Figure 4. (a) Mass-radius diagrams for hybrid EOSs combining the Ska and the CFL bag model with a4=1. Parameters B and Δ are chosen so that the phase transition occurs when the hadronic branch reaches 2M. Each curve is drawn until the last slow stable configuration. The observational data are the same as in Figure 3. (b) Comparing the non CFL contribution to the quark EOSs to pQCD results (Kurkela et al., 2010; Fraga et al., 2014). The ‘x’ marks denote the last slow stable configuration.

Notably, the central baryon chemical potential range that is being considered, in order to achieve the reconciliation of the HESS J1731-347 constraints, reaches such high values that it may cross the regime where the results of pQCD, for the EOS of strongly interacting quark matter, are potentially reliable. Typically, pQCD is considered to be credible at densities around 40n0 (where n0 is the nuclear saturation density), which correspond to values of baryon chemical potential that most likely exceed 2 GeV (Abbott et al., 2025; Kurkela et al., 2010; Fraga et al., 2014; Kurkela et al., 2014; Annala et al., 2018; Annala et al., 2020; Kurkela et al., 2024). For that matter, we wanted to test whether the constructed EOSs are compatible to pQCD results (for values beyond 2 GeV), at least for values of baryon chemical potential equal to those found at the center of the last slow stable configurations. In Figure 4b, we plotted the quark EOSs, compared to O(as2) results of pQCD (Kurkela et al., 2010) [employing the fit to the numerical data provided by Fraga et al. (2014)]. Notably, each curve depicts the pressure (divided by the Stefan-Boltzmann pressure PSB=3(μB/3)4/(4π23)) as a function of the chemical potential, neglecting the contribution of color superconductivity [as such a contribution is not encapsulated in the considered pQCD calculations (Fraga et al., 2014)]. Finally, the ’x’ appearing on each curve denotes the terminal configuration for each EOS. As is evident, for values of Δ up to 400 MeV the derived EOSs are compatible to the results of pQCD. However, this is not the case for larger Δ values. The latter is of utmost importance as it places strong constraints on the possible range of parameters that may enable the explanation of the HESS J1731-347 constraints in a slow phase conversion scenario.

A final remark is appropriate regarding the fact that while the existence of slow stable hybrid stars is theoretically intriguing it is not clear how such objects are born. In that sense, future work on possible formation scenarios of such objects would be of utmost importance. Some recent progress on the astrophysical paths that may lead to the existence of twin stars (but for a hybrid branch respecting the turning point criterion) has been made in the work of Naseri et al. (2024). In that direction, we expect that future research will hopefully shed light on the possible existence of slow stable hybrid star branches.

4 Conclusion

In this paper, we have presented a systematic study of the CFL MIT bag model in light of recent observations of low-mass compact stars. We have shown that the intriguing measurements of HESS J1731-347 and PSR J1231-1411 can be simultaneously explained within the framework of pure (absolutely stable) CFL matter, while also satisfying the maximum mass constraint set by PSR J0952-0607 and the latest multimessenger constraints on compact star masses and radii (PSR J0030+0451 and PSR J0437-4715). The parameter space consistent with all these measurements yields MR relations with maximum masses potentially exceeding 3M, a result that opens the discussion of quark stars populating the lower mass gap.

Notably, in the case of absolutely stable SQM, we worked under the assumption that only one state of matter may appear in the stellar interior. However, it is worth mentioning that a new exotic scenario has been reported in the literature suggesting the possible reappearance of hadrons, at large densities, even if SQM represents the true ground state of matter (Zhang and Ren, 2023; Negreiros et al., 2025; Zhang et al., 2024). Thus, an interesting direction for future work would be to investigate how such a consideration might alter the results reported in the present study.

When considering the framework of hybrid stars, we found that although a hybrid branch originating in the PSR J1231-1411 mass-radius region [as provided by Salmi et al. (2024a)] can also accommodate the HESS J1731-347 measurement, the resulting maximum masses remain well below the well-established threshold of 2M. While altering the pQCD related parameter a4 or incorporating a repulsive vector interaction increases the predicted maximum mass, the induced stiffening does not allow for the explanation of the HESS J1731-347 constraints (assuming a stiff hadronic EOS). However, the situation changes if the CCO in the HESS J1731-347 remnant is interpreted as a slow stable hybrid star. In this scenario, we adopted a sufficiently strong phase transition at high central pressures, requiring high values for both Δ and B, to produce a hybrid branch that initiates at 2M and drops thereafter, while maintaining stability. We concluded that such a hybrid EOS can be compatible with all of the aforementioned measurements, while also satisfying the demands of pQCD at extremely high density.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author. The data associated with this study are available from the authors upon reasonable request.

Author contributions

KK: Writing – original draft, Writing – review and editing. PL-P: Writing – original draft, Writing – review and editing. CM: Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. PL-P. acknowledges that the research work was supported by the Hellenic Foundation for Research and Innovation (HFRI) under the fifth Call for HFRI PhD Fellowships (Fellowship Number: 19175).

Acknowledgments

The authors would like to thank Y. Kini for providing the data for the contour regions for XTE J1814-338.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abac, A. G., Abbott, R., Abouelfettouh, I., Acernese, F., Ackley, K., Adhicary, S., et al. (2024). Observation of gravitational Waves from the Coalescence of a 2.5 − −4.5 M compact Object and a neutron star. Astrophys. J. Lett. 970, L34. doi:10.3847/2041-8213/ad5beb

CrossRef Full Text | Google Scholar

Abbot, B. P., Abbott, R., Abbott, T., Acernese, F., Ackley, K., Adams, C., et al. (2017). (LIGO scientific, virgo), GW170817: Observation of gravitational Waves from a binary neutron star inspiral. Phys. Rev. Lett. 119, 161101. doi:10.1103/physrevlett.119.161101

PubMed Abstract | CrossRef Full Text | Google Scholar

Abbott, B. P., Abbott, R., Abbott, T. D., Abraham, S., Acernese, F., Ackley, K., et al. (2020a). (LIGO scientific, Virgo), GW190425: observation of a compact binary coalescence with total mass ∼ 3.4M. Astrophysical J. Lett. 892. doi:10.3847/2041-8213/ab75f5

CrossRef Full Text | Google Scholar

Abbot, R., Abbott, T., Acernese, F., Ackley, K., Adams, C., Adhikari, N., et al. (2023). (KAGRA, VIRGO, LIGO scientific), GWTC-3: compact binary coalescences Observed by LIGO and Virgo during the second Part of the third observing run. Phys. Rev. X 13, 041039. doi:10.1103/physrevx.13.041039

CrossRef Full Text | Google Scholar

Abbott, R., Abbott, T. D., Abraham, S., Acernese, F., Ackley, K., Adams, C., et al. (2020b). (LIGO scientific, virgo), 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. doi:10.3847/2041-8213/ab960f

CrossRef Full Text | Google Scholar

Abbott, R., Detmold, W., Illa, M., Parreño, A., Perry, R. J., Romero-López, F., et al. (2025). NPLQCD collaboration, QCD Constraints on isospin-dense Matter and the nuclear Equation of state. Phys. Rev. Lett. 134, 011903. doi:10.1103/physrevlett.134.011903

PubMed Abstract | CrossRef Full Text | Google Scholar

Adhikari, D., Albataineh, H., Androic, D., Aniol, K., Armstrong, D. S., Averett, T., et al. (2021). Accurate Determination of the neutron skin Thickness of 208Pb through parity-Violation in electron scattering. Phys. Rev. Lett. 126, 172502. doi:10.1103/PhysRevLett.126.172502

PubMed Abstract | CrossRef Full Text | Google Scholar

Adhikari, D., Albataineh, H., Androic, D., Aniol, K. A., Armstrong, D. S., Averett, T., et al. (2022). Precision Determination of the neutral weak form Factor of 48Ca. Phys. Rev. Lett. 129, 042501. doi:10.1103/PhysRevLett.129.042501

PubMed Abstract | CrossRef Full Text | Google Scholar

Alford, J. A. J., and Halpern, J. P. (2023). Do central compact objects have carbon atmospheres? Astrophys. J. 944, 36. doi:10.3847/1538-4357/acaf55

CrossRef Full Text | Google Scholar

Alford, M., Braby, M., Paris, M., and Reddy, S. (2005a). Hybrid stars that masquerade as neutron stars. APJ 629, 969–978. doi:10.1086/430902

CrossRef Full Text | Google Scholar

Alford, M., Kouvaris, C., and Rajagopal, K. (2005b). Evaluating the gapless color-flavor locked phase. Phys. Rev. D. 71, 054009. doi:10.1103/physrevd.71.054009

CrossRef Full Text | Google Scholar

Alford, M., Rajagopal, K., Reddy, S., and Wilczek, F. (2001). Minimal color-flavor-locked–nuclear interface. Phys. Rev. D. 64, 074017. doi:10.1103/physrevd.64.074017

CrossRef Full Text | Google Scholar

Alford, M., Rajagopal, K., and Wilczek, F. (1999). Color-flavor locking and chiral symmetry breaking in high density QCD. Nucl. Phys. B 537, 443–458. doi:10.1016/s0550-3213(98)00668-3

CrossRef Full Text | Google Scholar

Alford, M., and Sedrakian, A. (2017). Compact stars with sequential QCD phase transitions. Phys. Rev. Lett. 119, 161104. doi:10.1103/physrevlett.119.161104

PubMed Abstract | CrossRef Full Text | Google Scholar

Alford, M. G. (2001). Color-superconducting quark matter, ann. Rev. Of nucl. And part. Phys. 51, 131. doi:10.1146/annurev.nucl.51.101701.132449

CrossRef Full Text | Google Scholar

Alford, M. G., Schmitt, A., Rajagopal, K., and Schäfer, T. (2008). Color superconductivity in dense quark matter. Rev. Mod. Phys. 80, 1455–1515. doi:10.1103/revmodphys.80.1455

CrossRef Full Text | Google Scholar

Alvarez-Castillo, D. E., Blaschke, D. B., Grunfeld, A. G., and Pagura, V. P. (2019). Third family of compact stars within a nonlocal chiral quark model equation of state. Phys. Rev. D. 99, 063010. doi:10.1103/physrevd.99.063010

CrossRef Full Text | Google Scholar

Annala, E., Gorda, T., Kurkela, A., Nättilä, J., and Vuorinen, A. (2020). Evidence for quark-matter cores in massive neutron stars. Nat. Phys. 16 (9), 907–910. doi:10.1038/s41567-020-0914-9

CrossRef Full Text | Google Scholar

Annala, E., Gorda, T., Kurkela, A., and Vuorinen, A. (2018). Gravitational-wave constraints on the neutron-star-matter equation of state. Phys. Rev. Lett. 120, 172703. doi:10.1103/physrevlett.120.172703

PubMed Abstract | CrossRef Full Text | Google Scholar

Antoniadis, J., Freire, P., Wex, N., Tauris, T. M., Lynch, R. S., van Kerkwijk, M. H., et al. (2013). A massive pulsar in a compact relativistic binary. Science 340, 448, 1233232. doi:10.1126/science.1233232

PubMed Abstract | CrossRef Full Text | Google Scholar

Barr, E. D., Dutta, A., Freire, P. C. C., Cadelano, M., Gautam, T., Kramer, M., et al. (2024). A pulsar in a binary with a compact object in the mass gap between neutron stars and black holes. Science 383, 275–279. doi:10.1126/science.adg3005

PubMed Abstract | CrossRef Full Text | Google Scholar

Baym, G., Hatsuda, T., Kojo, T., Powell, P. D., Song, Y., and Takatsuka, T. (2018). From hadrons to quarks in neutron stars: a review. Rep. Prog. Phys. 81, 056902. doi:10.1088/1361-6633/aaae14

PubMed Abstract | CrossRef Full Text | Google Scholar

Baym, G., Pethick, P., and Sutherland, P. (1971). The ground state of matter at high densities: equation of state and stellar models. Astrophys. J. 170, 299. doi:10.1086/151216

CrossRef Full Text | Google Scholar

Benic, S., Blaschke, D., Alvarez-Castillo, D. E., Fischer, T., and Typel, S. (2015). A new quark-hadron hybrid equation of state for astrophysics — I. High-mass twin compact stars. Astron. Astrophys. 577, A40. doi:10.1051/0004-6361/201425318

CrossRef Full Text | Google Scholar

Blaschke, D., Ayriyan, A., Alvarez-Castillo, D. E., and Grigorian, H. (2020). Was GW170817 a canonical neutron star merger? Bayesian analysis with a third family of compact stars. Universe 6, 81. doi:10.3390/universe6060081

CrossRef Full Text | Google Scholar

Bombaci, I., Panda, P. K., Providência, C., and Vidaña, I. (2008). Metastability of hadronic compact stars. Phys. Rev. D. 77, 083002. doi:10.1103/physrevd.77.083002

CrossRef Full Text | Google Scholar

Bombaci, I., Logoteta, D., Panda, P. K., Providência, C., and Vidaña, I. (2009). Quark matter nucleation in hot hadronic matter. Phys. Lett. 680B, 448. doi:10.1016/j.physletb.2009.09.039

CrossRef Full Text | Google Scholar

Brodie, L., and Haber, A. (2023). Nuclear and hybrid equations of state in light of the low-mass compact star in HESS J1731-347. Phys. Rev. C 108, 025806. doi:10.1103/physrevc.108.025806

CrossRef Full Text | Google Scholar

Burgio, G. F., Das, H. C., and Vidaña, I. (2024). The nuclear symmetry energy and the neutron skin thickness in nuclei. Front. Astron. Space Sci. 11, 1505560. doi:10.3389/fspas.2024.1505560

CrossRef Full Text | Google Scholar

Char, P., and Biswas, B. (2024). The compact object of HESS J1731-347 and its implication on neutron star matter. arXiv:2408.15220. doi:10.48550/arXiv.2408.15220

CrossRef Full Text | Google Scholar

Choudhury, D., Salmi, T., Vinciguerra, S., Riley, T. E., Kini, Y., Watts, A. L., et al. (2024). A NICER view of the nearest and brightest millisecond pulsar: PSR j0437–4715. Astrophys. J. Lett. 971, L20. doi:10.3847/2041-8213/ad5a6f

CrossRef Full Text | Google Scholar

Christian, J.-E., and Schaffner-Bielich, J. (2021). Supermassive neutron stars rule out twin stars. Phys. Rev. D. 103, 063042. doi:10.1103/physrevd.103.063042

CrossRef Full Text | Google Scholar

Christian, J.-E., and Schaffner-Bielich, J. (2022). Confirming the existence of twin stars in a NICER way. Astrophys. J. 935, 122. doi:10.3847/1538-4357/ac75cf

CrossRef Full Text | Google Scholar

Cromartie, H., Fonseca, E., Ransom, S., Demorest, P. B., Arzoumanian, Z., Blumer, H., et al. (2020). Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar. Nat. Astron. 4, 72–76. doi:10.1038/s41550-019-0880-2

CrossRef Full Text | Google Scholar

Danielewicz, P., and Lee, J. (2009). Symmetry energy I: semi-infinite matter. Nucl. Phys. A 818, 36–96. doi:10.1016/j.nuclphysa.2008.11.007

CrossRef Full Text | Google Scholar

Das, H. C., and Lopes, L. L. (2023). Anisotropic strange stars in the spotlight: unveiling constraints through observational data. Mon. Not. Roy. Astron. Soc. 525 (3), 3571–3575. doi:10.1093/mnras/stad2554

CrossRef Full Text | Google Scholar

Di Clemente, F., Drago, A., and Pagliara, G. (2024). Is the compact object associated with HESS j1731-347 a strange quark star? A possible astrophysical scenario for its formation. Astrophys. J. 967, 159. doi:10.3847/1538-4357/ad445b

CrossRef Full Text | Google Scholar

Doroshenko, V., Suleimanov, V., Pühlhofer, G., and Santangelo, A. (2022). A strangely light neutron star within a supernova remnant. Nat. Astron. 6, 1444–1451. doi:10.1038/s41550-022-01800-1

CrossRef Full Text | Google Scholar

Farhi, E., and Jaffe, R. L. (1984). Strange matter. Phys. Rev. D. 30, 2379–2390. doi:10.1103/physrevd.30.2379

CrossRef Full Text | Google Scholar

Fraga, E. S., Kurkela, A., and Vuorinen, A. (2014). Interacting quark matter equation of state for compact stars. Astrophys. J. Lett. 781, L25. doi:10.1088/2041-8205/781/2/l25

CrossRef Full Text | Google Scholar

Fraga, E. S., Pisarski, R. D., and Schaffner-Bielich, J. (2001). Small, dense quark stars from perturbative QCD. Phys. Rev. D. 63, 121702(R). doi:10.1103/physrevd.63.121702

CrossRef Full Text | Google Scholar

Gao, B., Yan, Y., and Harada, M. (2024). Reconciling constraints from the supernova remnant HESS J1731-347 with the parity doublet model. Phys. Rev. C 109, 065807. doi:10.1103/physrevc.109.065807

CrossRef Full Text | Google Scholar

Gholami, H., Rather, I. A., Hofmann, M., Buballa, M., and Schäffner-Bielich, J. (2024). Astrophysical constraints on color-superconducting phases in compact stars within the RG-consistent NJL model. Phys. Rev. D, 111, 103034. doi:10.1103/PhysRevD.111.103034

CrossRef Full Text | Google Scholar

Glendenning, N. K., and Kettner, C. (2000). Possible third family of compact stars more dense than neutron stars. Astron. Astrophys. 353, L9–L12. Available online at: https://ui.adsabs.harvard.edu/abs/2000A%26A...353L...9G/abstract.

Google Scholar

Gulminelli, F., and Raduta, Ad. R. (2015). Unified treatment of subsaturation stellar matter at zero and finite temperature. Phys. Rev. C 92, 055803. doi:10.1103/physrevc.92.055803

CrossRef Full Text | Google Scholar

Heiselberg, H., and Hjorth-Jensen, M. (2000). Phases of dense matter in neutron stars. Phys. Rep. 328, 237–327. doi:10.1016/s0370-1573(99)00110-6

CrossRef Full Text | Google Scholar

Horvath, J. E., Benvenuto, O. G., and Vucetich, H. (1992). Nucleation of strange matter in dense stellar cores. Phys. Rev. D. 45, 3865–3868. doi:10.1103/physrevd.45.3865

PubMed Abstract | CrossRef Full Text | Google Scholar

Horvath, J. E., Roch, L. S., de Sá, L. M., Moraes, P. H. R. S., Barão, L. G., de Avellar, M. G. B., et al. (2023). A light strange star in the remnant HESS J1731-347: minimal consistency checks. Astron. Astrophys. 672, L11. doi:10.1051/0004-6361/202345885

CrossRef Full Text | Google Scholar

Huang, K., Shen, H., Hu, J., and Zhang, Y. (2024). Hadronic equation of state of low-mass neutron stars from a relativistic mean-field model with tensor couplings. Phys. Rev. D. 109, 043036. doi:10.1103/physrevd.109.043036

CrossRef Full Text | Google Scholar

Iida, K., and Sato, K. (1998). Effects of hyperons on the dynamical deconfinement transition in cold neutron star matter. Phys. Rev. C 58, 2538–2559. doi:10.1103/physrevc.58.2538

CrossRef Full Text | Google Scholar

Jayasinghe, T., Stanek, K. Z., Thompson, T. A., Kochanek, C. S., Rowan, D. M., Vallely, P. J., et al. (2021). A Unicorn in Monoceros: the 3Mdark companion to the bright, nearby red giant V723 Mon is a non-interacting, mass-gap black hole candidate. Mon. Notices R. Astronomical Soc. 504, 2577–2602. doi:10.1093/mnras/stab907

CrossRef Full Text | Google Scholar

Jiménez, J. C., Lazzari, L., and Gonçalves, V. P. (2024). How the QCD trace anomaly behaves at the core of twin stars? Phys. Rev. D. 110, 114014. doi:10.1103/physrevd.110.114014

CrossRef Full Text | Google Scholar

Kanakis-Pegios, A., Petousis, V., Veselský, M., Leja, J., and Moustakidis, Ch.C. (2024). Constraints for the X17 boson from compact objects observations. Phys. Rev. D. 109, 043028. doi:10.1103/physrevd.109.043028

CrossRef Full Text | Google Scholar

Kini, Y., Salmi, T., Vinciguerra, S., Watts, A. L., Bilous, A., Galloway, D. K., et al. (2024). Constraining the properties of the thermonuclear burst oscillation source XTE j1814-338 through pulse profile modelling. Mon. Not. Roy. Astron. Soc. 535, 1507–1525. doi:10.1093/mnras/stae2398

CrossRef Full Text | Google Scholar

Köhler, H. S. (1976). Skyrme force and the mass formula. Nucl. Phys. A 258 (2), 301–316. doi:10.1016/0375-9474(76)90008-7

CrossRef Full Text | Google Scholar

Kubis, S., Wojcik, W., Castillo, D. A., and Zabari, N. (2023). Relativistic mean-field model for the ultracompact low-mass neutron star HESS J1731-347. Phys. Rev. C 108, 045803. doi:10.1103/physrevc.108.045803

CrossRef Full Text | Google Scholar

Kumar, M., Kumar, S., Thakur, V., Kumar, R., Agrawal, B. K., and Dhiman, S. K. (2023). CREX- and PREX-II-motivated relativistic interactions and their implications for the bulk properties of nuclear matter and neutron stars. Phys. Rev. C 107, 055801. doi:10.1103/physrevc.107.055801

CrossRef Full Text | Google Scholar

Kumar, M., Kumar, S., Thakur, V., Kumar, R., Agrawal, B. K., and Dhiman, S. K. (2024). Trends of nuclear matter and neutron-star observables post CREX-PREX-II results within the relativistic mean field model. Phys. Rev. C 110, 055802. doi:10.1103/physrevc.110.055802

CrossRef Full Text | Google Scholar

Kurkela, A., Fraga, E. S., Schaffner-Bielich, J., and Vuorinen, A. (2014). Constraining neutron star matter with quantum chromodynamics. Astrophys. J. 789, 127. doi:10.1088/0004-637x/789/2/127

CrossRef Full Text | Google Scholar

Kurkela, A., Rajagopal, K., and Steinhorst, R. (2024). Astrophysical equation-of-state constraints on the color-superconducting gap. Phys. Rev. Lett. 132, 262701. doi:10.1103/physrevlett.132.262701

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurkela, A., Romatschke, P., and Vuorinen, A. (2010). Cold quark matter. Phys. Rev. D. 81, 105021. doi:10.1103/physrevd.81.105021

CrossRef Full Text | Google Scholar

Laskos-Patkos, P., Koliogiannis, P. S., and Moustakidis, Ch. C. (2024). Hybrid stars in light of the HESS J1731-347 remnant and the PREX-II experiment. Phys. Rev. D. 109, 063017. doi:10.1103/physrevd.109.063017

CrossRef Full Text | Google Scholar

Laskos-Patkos, P., Lalazissis, G. A., Wang, S., Meng, J., Ring, P., and Moustakidis, Ch. C. (2025). Speed of sound bounds and first-order phase transitions in compact stars. Phys. Rev. C 111, 025801. doi:10.1103/physrevc.111.025801

CrossRef Full Text | Google Scholar

Laskos-Patkos, P., and Moustakidis, Ch.C. (2025). XTE J1814-338: a potential hybrid star candidate. Phys. Rev. D. 111, 063058. doi:10.1103/physrevd.111.063058

CrossRef Full Text | Google Scholar

Lattimer, J. M., and Prakash, M. (2001). Neutron star structure and the equation of state. Astrophys. J. 550, 426–442. doi:10.1086/319702

CrossRef Full Text | Google Scholar

Li, J. J., and Sedrakian, A. (2023). Baryonic models of ultra-low-mass compact stars for the central compact object in HESS J1731-347. Phys. Lett. 844B, 138062. doi:10.1016/j.physletb.2023.138062

CrossRef Full Text | Google Scholar

Li, J. J., Sedrakian, A., and Alford, M. (2020). Relativistic hybrid stars with sequential first-order phase transitions and heavy-baryon envelopes. Phys. Rev. D. 101, 063022. doi:10.1103/physrevd.101.063022

CrossRef Full Text | Google Scholar

Li, J. J., Sedrakian, A., and Alford, M. (2023a). Relativistic hybrid stars with sequential first-order phase transitions in light of multimessenger constraints. Astrophys. J. 944, 206. doi:10.3847/1538-4357/acb688

CrossRef Full Text | Google Scholar

Li, J. J., Sedrakian, A., and Alford, M. (2023b). Ultracompact hybrid stars consistent with multimessenger astrophysics. Phys. Rev. D. 107, 023018. doi:10.1103/physrevd.107.023018

CrossRef Full Text | Google Scholar

Li, J. J., Sedrakian, A., and Alford, M. (2024). Hybrid star models in the light of new multimessenger data. Astrophys. J. 967, 116. doi:10.3847/1538-4357/ad4295

CrossRef Full Text | Google Scholar

Li, J. J., Sedrakian, A., and Alford, M. (2025). Confronting new NICER mass-radius measurements with phase transition in dense matter and twin compact stars. JCAP 02, 002. doi:10.1088/1475-7516/2025/02/002

CrossRef Full Text | Google Scholar

Lopes, L. L., and Issifu, A. (2025). XTE J1814-338 as a dark matter admixed neutron star. Phys. Dark. Univ. 48, 101922. doi:10.1016/j.dark.2025.101922

CrossRef Full Text | Google Scholar

Lugones, G., and Horvath, J. E. (2002). Color-flavor locked strange matter. Phys. Rev. D. 66, 074017. doi:10.1103/physrevd.66.074017

CrossRef Full Text | Google Scholar

M. Veselský, , Petousis, V., Koliogiannis, P., Moustakidis, C., and Leja, J. (2025). Simultaneous explanation of XTE J1814-338 and HESS J1731-347 objects using K− and condensates. Phys. Rev. D. 111, L061308. doi:10.1103/physrevd.111.l061308

CrossRef Full Text | Google Scholar

Mariani, M., Ranea-Sandoval, I. F., Lugones, G., and Orsaria, M. G. (2024). Could a slow stable hybrid star explain the central compact object in HESS J1731-347? Phys. Rev. D. 110, 043026. doi:10.1103/physrevd.110.043026

CrossRef Full Text | Google Scholar

Miller, M. C., Lamb, F. K., Dittmann, A. J., Bogdanov, S., Arzoumanian, Z., Gendreau, K. C., et al. (2019). PSR J0030+0451 mass and radius from NICER data and implications for the properties of neutron star matter. Astrophys. J. Lett. 887, L24. doi:10.3847/2041-8213/ab50c5

CrossRef Full Text | Google Scholar

Miyatsu, T., Cheoun, M.-K., Kim, K., and Saito, K. (2023). Can the PREX-2 and CREX results be understood by relativistic mean-field models with the astrophysical constraints? Phys. Lett. B 843, 138013. doi:10.1016/j.physletb.2023.138013

CrossRef Full Text | Google Scholar

Naseri, M., Bozzola, G., and Paschalidis, V. (2024). Exploring pathways to forming twin stars. Phys. Rev. D. 110, 044037. doi:10.1103/physrevd.110.044037

CrossRef Full Text | Google Scholar

Negreiros, R., Zhang, C., and Xu, R. (2025). Rotational properties of inverted hybrid stars. Phys. Rev. D. 111, 063026. doi:10.1103/physrevd.111.063026

CrossRef Full Text | Google Scholar

Oikonomou, P. T., and Moustakidis, Ch.C. (2023). Color-flavor locked quark stars in light of the compact object in the HESS J1731-347 and the GW190814 event. Phys. Rev. D. 108, 063010. doi:10.1103/physrevd.108.063010

CrossRef Full Text | Google Scholar

Olesen, M. L., and Madsen, J. (1994). Nucleation of quark matter bubbles in neutron stars. Phys. Rev. D. 49, 2698–2702. doi:10.1103/physrevd.49.2698

PubMed Abstract | CrossRef Full Text | Google Scholar

Oppenheimer, J. R., and Volkoff, G. M. (1939). On massive neutron cores. Phys. Rev. 55, 374–381. doi:10.1103/physrev.55.374

CrossRef Full Text | Google Scholar

Pal, S., Podder, S., and Chaudhuri, G. (2025). Is the central compact object in HESS j1731-347 a hybrid star with a quark core? An analysis with the constant speed of sound parameterization. ApJ 983, 24. doi:10.3847/1538-4357/adbc6b

CrossRef Full Text | Google Scholar

Pereira, J. P., Flores, C. V., and Lugones, G. (2018). Phase transition effects on the dynamical stability of hybrid neutron stars. Astrophys. J. 860, 12. doi:10.3847/1538-4357/aabfbf

CrossRef Full Text | Google Scholar

Qi, L., Zheng, S., Zhang, J., Ge, M., Li, A., Zhang, S. N., et al. (2025). PSR j1231-1411 revisited: pulse profile analysis of X-ray observation. Astrophys. J. 981, 99. doi:10.3847/1538-4357/adb42f

CrossRef Full Text | Google Scholar

Rather, I. A., Marquez, K. D., Backes, B. C., Panotopoulos, G., and Lopes, I. (2024). Radial oscillations of hybrid stars and neutron stars including delta baryons: the effect of a slow quark phase transition. JCAP 05 2024, 130. doi:10.1088/1475-7516/2024/05/130

CrossRef Full Text | Google Scholar

Rather, I. A., Panotopoulos, G., and Lopes, I. (2023). Quark models and radial oscillations: decoding the HESS J1731-347 compact object’s equation of state. Eur. Phys. J. C 83, 1065. doi:10.1140/epjc/s10052-023-12223-1

CrossRef Full Text | Google Scholar

Reed, B. T., Fattoyev, F. J., Horowitz, C. J., and Piekarewicz, J. (2021). Implications of PREX-2 on the equation of state of neutron-rich matter. Phys. Rev. Lett. 126, 172503. doi:10.1103/physrevlett.126.172503

PubMed Abstract | CrossRef Full Text | Google Scholar

Reed, B. T., Fattoyev, F. J., Horowitz, C. J., and Piekarewicz, J. (2024). Density dependence of the symmetry energy in the post–PREX-CREX era. Phys. Rev. C 109, 035803. doi:10.1103/physrevc.109.035803

CrossRef Full Text | Google Scholar

Reinhard, P. G., Roca-Maza, X., and Nazarewicz, W. (2023). Combined theoretical Analysis of the parity-violating Asymmetry for48Ca and208Pb. Phys. Rev. Lett. 129, 232501. doi:10.1103/physrevlett.129.232501

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, J., and Zhang, C. (2020). Quantum nucleation of up-down quark matter and astrophysical implications. Phys. Rev. D. 102, 083003. doi:10.1103/physrevd.102.083003

CrossRef Full Text | Google Scholar

Riley, T. E., Watts, A. L., Bogdanov, S., Ray, P. S., Ludlam, R. M., Guillot, S., et al. (2019). A NICER view of PSR J0030+0451: millisecond pulsar parameter estimation. Astrophys. J. Lett. 887, L21. doi:10.3847/2041-8213/ab481c

CrossRef Full Text | Google Scholar

Romani, R. G., Kandel, D., Filippenko, A. V., Brink, T. G., and Zheng, W. (2022). PSR j0952-0607: the fastest and heaviest known galactic neutron star. Astrophys. J. Lett. 934, L17. doi:10.3847/2041-8213/ac8007

CrossRef Full Text | Google Scholar

Routaray, P., Das, H. C., Amrit Pattnaik, J., and Kumar, B. (2024). Constraining neutron star properties and dark matter admixture with the NITR-I equation of state: insights from observations and universal relations. Int. J. Modern Phy. E, 33, 11. doi:10.1142/s0218301324500526

CrossRef Full Text | Google Scholar

Sagun, V., Giangrandi, E., Dietrich, T., Ivanytskyi, O., Negreiros, R., and Providencia, C. (2023). What is the nature of the HESS j1731-347 compact object? Astrophys. J. 958 (49), 49. doi:10.3847/1538-4357/acfc9e

CrossRef Full Text | Google Scholar

Salinas, M., and Piekarewicz, J. (2024). Impact of tensor couplings with scalar mixing on covariant energy density functionals. Phys. Rev. C 109, 045807. doi:10.1103/physrevc.109.045807

CrossRef Full Text | Google Scholar

Salmi, T., Choudhury, D., Kini, Y., Riley, T. E., Vinciguerra, S., Watts, A. L., et al. (2024b). The radius of the high-mass pulsar PSR J0740+6620 with 3.6 yr of NICER data. Astrophys. J. Lett. 974, 294. doi:10.3847/1538-4357/ad5f1f

CrossRef Full Text | Google Scholar

Salmi, T., Deneva, J. S., Ray, P. S., Watts, A. L., Choudhury, D., Kini, Y., et al. (2024a). A NICER view of PSR j1231-1411: a complex case. Astrophys. J. 976, 58. doi:10.3847/1538-4357/ad81d2

CrossRef Full Text | Google Scholar

Samsing, J., and Hotokezaka, K. (2021). Populating the black hole mass gaps in stellar clusters: general relations and upper limits. ApJ 923, 126. doi:10.3847/1538-4357/ac2b27

CrossRef Full Text | Google Scholar

Shao, W. (2022). On the neutron star/black hole mass gap and black hole searches. Astron. Astrophys. 22, 122002. doi:10.1088/1674-4527/ac995e

CrossRef Full Text | Google Scholar

Tagami, S., Wakasa, T., and Yahiro, M. (2022). Slope parameters determined from CREX and PREX2. Results Phys. 43, 106037. doi:10.1016/j.rinp.2022.106037

CrossRef Full Text | Google Scholar

Tewari, S., Chatterjee, S., Kumar, D., and Mallick, R. (2024). Analyzing the dense matter equation of states in the light of the compact object HESS J1731-347. Phys. Rev. D 111, 103009 (2025). doi:10.1103/PhysRevD.111.103009

CrossRef Full Text | Google Scholar

Thompson, T. A., Kochanek, C. S., Stanek, K. Z., Badenes, C., Post, R. S., Jayasinghe, T., et al. (2019). A noninteracting low-mass black hole–giant star binary system. Science 366, 637–640. doi:10.1126/science.aau4005

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsaloukidis, L., Koliogiannis, P. S., Kanakis-Pegios, A., and Moustakidis, Ch.C. (2023). Twin stars as probes of the nuclear equation of state: effects of rotation through the PSR J0952-0607 pulsar and constraints via the tidal deformability from the GW170817 event. Phys. Rev. D. 107, 023012. doi:10.1103/physrevd.107.023012

CrossRef Full Text | Google Scholar

Typel, S., Oertel, M., Klähn, T., Chatterjee, D., Dexheimer, V., Ishizuka, C., et al. (2022). CompOSE core team collaboration, CompOSE reference manual. Eur. Phys. J. 58, 221. doi:10.1140/epja/s10050-022-00847-y

CrossRef Full Text | Google Scholar

Vinciguerra, S., Salmi, T., Watts, A. L., Choudhury, D., Riley, T. E., Ray, P. S., et al. (2024). An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451. Astrophys. J. 961, 62. doi:10.3847/1538-4357/acfb83

CrossRef Full Text | Google Scholar

Weber, F. (2005). Strange quark matter and compact stars. Prog. Part. Nucl. Phys. 54 (1), 193–288. doi:10.1016/j.ppnp.2004.07.001

CrossRef Full Text | Google Scholar

Witten, E. (1984). Cosmic separation of phases. Phys. Rev. D. 30, 272–285. doi:10.1103/physrevd.30.272

CrossRef Full Text | Google Scholar

Yang, S.-H., and Pi, C.-M. (2024). Color-flavor locked strange stars admixed with mirror dark matter and the observations of compact stars, JCAP09 2024, 052.

Google Scholar

Yang, S.-H., Pi, C.-M., and Weber, F. (2025). Strange stars admixed with mirror dark matter: confronting observations of XTE J1814-338. Phys. Rev. D. 111, 043037. doi:10.1103/physrevd.111.043037

CrossRef Full Text | Google Scholar

Yang, S.-H., Pi, C.-M., Zheng, X.-P., and Weber, F. (2021). Constraints from compact star observations on non-Newtonian gravity in strange stars based on a density dependent quark mass model. Phys. Rev. D. 103, 043012. doi:10.1103/physrevd.103.043012

CrossRef Full Text | Google Scholar

Yang, S.-H., Pi, C.-M., Zheng, X.-P., and Weber, F. (2023). Confronting strange stars with compact-star observations and new physics. Universe 9, 202. doi:10.3390/universe9050202

CrossRef Full Text | Google Scholar

Zhang, C., Luo, Y., Li, H. B., Shao, L., and Xu, R. (2024). Radial and nonradial oscillations of inverted hybrid stars. Phys. Rev. D. 109, 063020. doi:10.1103/physrevd.109.063020

CrossRef Full Text | Google Scholar

Zhang, C., and Ren, J. (2023). Hybrid stars may have an inverted structure. Phys. Rev. D. 108, 063012. doi:10.1103/physrevd.108.063012

CrossRef Full Text | Google Scholar

Zhang, N.-B., and Li, B.-A. (2025). Impact of the nuclear equation of state on the formation of twin stars. Eur. Phys. J. A 61, 31. doi:10.1140/epja/s10050-025-01497-6

CrossRef Full Text | Google Scholar

Zhou, T., and Huang, C. (2025). Hidden twin star solutions from an agnostic speed-of-sound model: confronting XTE j1814-338’s extreme compactness. arXiv:2504, 08662. doi:10.48550/arXiv.2504.08662

CrossRef Full Text | Google Scholar

Keywords: neutron stars, quark stars, color-flavor locked matter, hadron-quark phase transition, hybrid stars, equation of state

Citation: Kourmpetis K, Laskos-Patkos P and Moustakidis CC (2025) Confronting recent light compact star observations with color-flavor locked quark matter. Front. Astron. Space Sci. 12:1600563. doi: 10.3389/fspas.2025.1600563

Received: 26 March 2025; Accepted: 11 June 2025;
Published: 22 July 2025.

Edited by:

Armen Sedrakian, University of Wrocław, Poland

Reviewed by:

Malte Albrecht, Jefferson Lab (DOE), United States
Giuseppe Pagliara, University of Ferrara, Italy
Milva Orsaria, National Scientific and Technical Research Council (CONICET), Argentina

Copyright © 2025 Kourmpetis, Laskos-Patkos and Moustakidis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ch. C. Moustakidis, bW91c3Rha2lAYXV0aC5ncg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.