Abstract
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 . 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; ). 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 and radius of the central compact object (CCO) in the HESS J1731-347 remnant. More precisely, the authors reported that, at the level, and 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; ; Huang et al., 2024; Li and Sedrakian, 2023; Kubis et al., 2023; Routaray et al., 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) (; 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 .
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 and km. In addition, by employing an uninformative prior, and constraining the radius between 10 and 14 km, they reported a similar value for and a higher value for . 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 and 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 (; ; ) 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) (; 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 (; ; 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 level. In addition, the fact that PSR J1231-1411 has a slightly higher mass and radius indicates that the curve may exhibit positive slope . 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 (; Lugones and Horvath, 2002; ), assuming the absolute stability of strange matter, in light of the aforementioned measurements (Doroshenko et al., 2022; Salmi et al., 2024a; ; Salmi et al., 2024b; Vinciguerra et al., 2024; Riley et al., 2019; Miller et al., 2019; ; ; Romani et al., 2022). In any case, this characteristic property 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 (; ; ; ; ; ; 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 ) 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 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 and , where represents the mass of the strange quark and the quark chemical potential, the pressure and energy density can be expressed as follows (Lugones and Horvath, 2002; )
where is a parameter that mimics the impact of perturbative QCD (pQCD) corrections (), represents the bag constant and . Notably, an analytical expression for can be obtained by combining Equations 1, 2:
Equation 3 reduces to the one presented in Lugones and Horvath (2002) for , 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, 56, at zero pressure and temperature (Farhi and Jaffe, 1984). This condition establishes an upper limit for the bag constant as a function of and , for a fixed value of :
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 , this condition is expressed as (Farhi and Jaffe, 1984):
This value decreases when considering lower values for . In this study, we use the constraint derived with (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 ():
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), , , 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 and values to identify the curves that are compatible with the light neutron star measurements for HESS J1731-347 and PSR J1231-1411 (at the 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 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 parameter we employed two distinct values, and . The latter value was selected as it results to EOSs similar to the quark matter EOS from pQCD [e.g., see , 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 (). 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 pairs that are compatible with the above measurements for two different values of , 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 considering Case I, while Figure 1b shows the same graph for . As is evident, although a smaller 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 and respectively.
FIGURE 1
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 pairs (above ). The latter is demonstrated in Figure 2b, where the maximum masses for all combinations of the purple region of Figure 1b are calculated (for the case).
FIGURE 2
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 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 estimation for PSR J0437-4715 (see Figure 2a, where we have plotted the 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 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 case). The resulting maximum masses span the range of approximately , 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 , points toward the possible existence of quark stars with masses exceeding . 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 (), and gravitational wave detections from binary mergers (; ; ; ; ).
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 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 (; )]. 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 (; Reed et al., 2021). Then, the subsequent CREX experiment extracted a puzzling value for the neutron skin thickness of calcium (), 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; ). 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; ). 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 (fixing ). Notably, the pairs of and 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 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 . The motivation to do that relies on the knowledge that reducing reduces the value of pressure for a given baryon chemical potential and, as a consequence, for the same parameter values of 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 (), leading to EOS stiffening. While the constructed EOSs achieve higher masses (still not the conservative 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 , where is the coupling constant and 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
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 , compared to the period of radial oscillations that the stellar configurations may undergo. In particular, if , 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 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 ), 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 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 . 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 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 . While a value of around 400 MeV is potentially compatible to their results, larger values may be disfavored [see Figure 2 in ].
FIGURE 4
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 (where is the nuclear saturation density), which correspond to values of baryon chemical potential that most likely exceed 2 GeV (; Kurkela et al., 2010; Fraga et al., 2014; Kurkela et al., 2014; ; ; 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 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 ) 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 relations with maximum masses potentially exceeding , 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 . While altering the pQCD related parameter 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 , to produce a hybrid branch that initiates at 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.
Statements
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
1
AbacA. G.AbbottR.AbouelfettouhI.AcerneseF.AckleyK.AdhicaryS.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. 10.3847/2041-8213/ad5beb
2
AbbotB. P.AbbottR.AbbottT.AcerneseF.AckleyK.AdamsC.et al (2017). (LIGO scientific, virgo), GW170817: Observation of gravitational Waves from a binary neutron star inspiral. Phys. Rev. Lett.119, 161101. 10.1103/physrevlett.119.161101
3
AbbottB. P.AbbottR.AbbottT. D.AbrahamS.AcerneseF.AckleyK.et al (2020a). (LIGO scientific, Virgo), GW190425: observation of a compact binary coalescence with total mass ∼ 3.4M⊙. Astrophysical J. Lett.892. 10.3847/2041-8213/ab75f5
4
AbbotR.AbbottT.AcerneseF.AckleyK.AdamsC.AdhikariN.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. X13, 041039. 10.1103/physrevx.13.041039
5
AbbottR.AbbottT. D.AbrahamS.AcerneseF.AckleyK.AdamsC.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. 10.3847/2041-8213/ab960f
6
AbbottR.DetmoldW.IllaM.ParreñoA.PerryR. J.Romero-LópezF.et al (2025). NPLQCD collaboration, QCD Constraints on isospin-dense Matter and the nuclear Equation of state. Phys. Rev. Lett.134, 011903. 10.1103/physrevlett.134.011903
7
AdhikariD.AlbatainehH.AndroicD.AniolK.ArmstrongD. S.AverettT.et al (2021). Accurate Determination of the neutron skin Thickness of 208Pb through parity-Violation in electron scattering. Phys. Rev. Lett.126, 172502. 10.1103/PhysRevLett.126.172502
8
AdhikariD.AlbatainehH.AndroicD.AniolK. A.ArmstrongD. S.AverettT.et al (2022). Precision Determination of the neutral weak form Factor of 48Ca. Phys. Rev. Lett.129, 042501. 10.1103/PhysRevLett.129.042501
9
AlfordJ. A. J.HalpernJ. P. (2023). Do central compact objects have carbon atmospheres?Astrophys. J.944, 36. 10.3847/1538-4357/acaf55
10
AlfordM.BrabyM.ParisM.ReddyS. (2005a). Hybrid stars that masquerade as neutron stars. APJ629, 969–978. 10.1086/430902
11
AlfordM.KouvarisC.RajagopalK. (2005b). Evaluating the gapless color-flavor locked phase. Phys. Rev. D.71, 054009. 10.1103/physrevd.71.054009
12
AlfordM.RajagopalK.ReddyS.WilczekF. (2001). Minimal color-flavor-locked–nuclear interface. Phys. Rev. D.64, 074017. 10.1103/physrevd.64.074017
13
AlfordM.RajagopalK.WilczekF. (1999). Color-flavor locking and chiral symmetry breaking in high density QCD. Nucl. Phys. B537, 443–458. 10.1016/s0550-3213(98)00668-3
14
AlfordM.SedrakianA. (2017). Compact stars with sequential QCD phase transitions. Phys. Rev. Lett.119, 161104. 10.1103/physrevlett.119.161104
15
AlfordM. G. (2001). Color-superconducting quark matter, ann. Rev. Of nucl. And part. Phys.51, 131. 10.1146/annurev.nucl.51.101701.132449
16
AlfordM. G.SchmittA.RajagopalK.SchäferT. (2008). Color superconductivity in dense quark matter. Rev. Mod. Phys.80, 1455–1515. 10.1103/revmodphys.80.1455
17
Alvarez-CastilloD. E.BlaschkeD. B.GrunfeldA. G.PaguraV. P. (2019). Third family of compact stars within a nonlocal chiral quark model equation of state. Phys. Rev. D.99, 063010. 10.1103/physrevd.99.063010
18
AnnalaE.GordaT.KurkelaA.NättiläJ.VuorinenA. (2020). Evidence for quark-matter cores in massive neutron stars. Nat. Phys.16 (9), 907–910. 10.1038/s41567-020-0914-9
19
AnnalaE.GordaT.KurkelaA.VuorinenA. (2018). Gravitational-wave constraints on the neutron-star-matter equation of state. Phys. Rev. Lett.120, 172703. 10.1103/physrevlett.120.172703
20
AntoniadisJ.FreireP.WexN.TaurisT. M.LynchR. S.van KerkwijkM. H.et al (2013). A massive pulsar in a compact relativistic binary. Science340, 448, 1233232. 10.1126/science.1233232
21
BarrE. D.DuttaA.FreireP. C. C.CadelanoM.GautamT.KramerM.et al (2024). A pulsar in a binary with a compact object in the mass gap between neutron stars and black holes. Science383, 275–279. 10.1126/science.adg3005
22
BaymG.HatsudaT.KojoT.PowellP. D.SongY.TakatsukaT. (2018). From hadrons to quarks in neutron stars: a review. Rep. Prog. Phys.81, 056902. 10.1088/1361-6633/aaae14
23
BaymG.PethickP.SutherlandP. (1971). The ground state of matter at high densities: equation of state and stellar models. Astrophys. J.170, 299. 10.1086/151216
24
BenicS.BlaschkeD.Alvarez-CastilloD. E.FischerT.TypelS. (2015). A new quark-hadron hybrid equation of state for astrophysics — I. High-mass twin compact stars. Astron. Astrophys.577, A40. 10.1051/0004-6361/201425318
25
BlaschkeD.AyriyanA.Alvarez-CastilloD. E.GrigorianH. (2020). Was GW170817 a canonical neutron star merger? Bayesian analysis with a third family of compact stars. Universe6, 81. 10.3390/universe6060081
26
BombaciI.PandaP. K.ProvidênciaC.VidañaI. (2008). Metastability of hadronic compact stars. Phys. Rev. D.77, 083002. 10.1103/physrevd.77.083002
27
BombaciI.LogotetaD.PandaP. K.ProvidênciaC.VidañaI. (2009). Quark matter nucleation in hot hadronic matter. Phys. Lett.680B, 448. 10.1016/j.physletb.2009.09.039
28
BrodieL.HaberA. (2023). Nuclear and hybrid equations of state in light of the low-mass compact star in HESS J1731-347. Phys. Rev. C108, 025806. 10.1103/physrevc.108.025806
29
BurgioG. F.DasH. C.VidañaI. (2024). The nuclear symmetry energy and the neutron skin thickness in nuclei. Front. Astron. Space Sci.11, 1505560. 10.3389/fspas.2024.1505560
30
CharP.BiswasB. (2024). The compact object of HESS J1731-347 and its implication on neutron star matter. arXiv:2408.15220. 10.48550/arXiv.2408.15220
31
ChoudhuryD.SalmiT.VinciguerraS.RileyT. E.KiniY.WattsA. L.et al (2024). A NICER view of the nearest and brightest millisecond pulsar: PSR j0437–4715. Astrophys. J. Lett.971, L20. 10.3847/2041-8213/ad5a6f
32
ChristianJ.-E.Schaffner-BielichJ. (2021). Supermassive neutron stars rule out twin stars. Phys. Rev. D.103, 063042. 10.1103/physrevd.103.063042
33
ChristianJ.-E.Schaffner-BielichJ. (2022). Confirming the existence of twin stars in a NICER way. Astrophys. J.935, 122. 10.3847/1538-4357/ac75cf
34
CromartieH.FonsecaE.RansomS.DemorestP. B.ArzoumanianZ.BlumerH.et al (2020). Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar. Nat. Astron.4, 72–76. 10.1038/s41550-019-0880-2
35
DanielewiczP.LeeJ. (2009). Symmetry energy I: semi-infinite matter. Nucl. Phys. A818, 36–96. 10.1016/j.nuclphysa.2008.11.007
36
DasH. C.LopesL. L. (2023). Anisotropic strange stars in the spotlight: unveiling constraints through observational data. Mon. Not. Roy. Astron. Soc.525 (3), 3571–3575. 10.1093/mnras/stad2554
37
Di ClementeF.DragoA.PagliaraG. (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. 10.3847/1538-4357/ad445b
38
DoroshenkoV.SuleimanovV.PühlhoferG.SantangeloA. (2022). A strangely light neutron star within a supernova remnant. Nat. Astron.6, 1444–1451. 10.1038/s41550-022-01800-1
39
FarhiE.JaffeR. L. (1984). Strange matter. Phys. Rev. D.30, 2379–2390. 10.1103/physrevd.30.2379
40
FragaE. S.KurkelaA.VuorinenA. (2014). Interacting quark matter equation of state for compact stars. Astrophys. J. Lett.781, L25. 10.1088/2041-8205/781/2/l25
41
FragaE. S.PisarskiR. D.Schaffner-BielichJ. (2001). Small, dense quark stars from perturbative QCD. Phys. Rev. D.63, 121702(R). 10.1103/physrevd.63.121702
42
GaoB.YanY.HaradaM. (2024). Reconciling constraints from the supernova remnant HESS J1731-347 with the parity doublet model. Phys. Rev. C109, 065807. 10.1103/physrevc.109.065807
43
GholamiH.RatherI. A.HofmannM.BuballaM.Schäffner-BielichJ. (2024). Astrophysical constraints on color-superconducting phases in compact stars within the RG-consistent NJL model. Phys. Rev. D, 111, 103034. 10.1103/PhysRevD.111.103034
44
GlendenningN. K.KettnerC. (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.
45
GulminelliF.RadutaAd. R. (2015). Unified treatment of subsaturation stellar matter at zero and finite temperature. Phys. Rev. C92, 055803. 10.1103/physrevc.92.055803
46
HeiselbergH.Hjorth-JensenM. (2000). Phases of dense matter in neutron stars. Phys. Rep.328, 237–327. 10.1016/s0370-1573(99)00110-6
47
HorvathJ. E.BenvenutoO. G.VucetichH. (1992). Nucleation of strange matter in dense stellar cores. Phys. Rev. D.45, 3865–3868. 10.1103/physrevd.45.3865
48
HorvathJ. E.RochL. S.de SáL. M.MoraesP. H. R. S.BarãoL. G.de AvellarM. G. B.et al (2023). A light strange star in the remnant HESS J1731-347: minimal consistency checks. Astron. Astrophys.672, L11. 10.1051/0004-6361/202345885
49
HuangK.ShenH.HuJ.ZhangY. (2024). Hadronic equation of state of low-mass neutron stars from a relativistic mean-field model with tensor couplings. Phys. Rev. D.109, 043036. 10.1103/physrevd.109.043036
50
IidaK.SatoK. (1998). Effects of hyperons on the dynamical deconfinement transition in cold neutron star matter. Phys. Rev. C58, 2538–2559. 10.1103/physrevc.58.2538
51
JayasingheT.StanekK. Z.ThompsonT. A.KochanekC. S.RowanD. M.VallelyP. J.et al (2021). A Unicorn in Monoceros: the 3M⊙dark 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. 10.1093/mnras/stab907
52
JiménezJ. C.LazzariL.GonçalvesV. P. (2024). How the QCD trace anomaly behaves at the core of twin stars?Phys. Rev. D.110, 114014. 10.1103/physrevd.110.114014
53
Kanakis-PegiosA.PetousisV.VeselskýM.LejaJ.MoustakidisCh.C. (2024). Constraints for the X17 boson from compact objects observations. Phys. Rev. D.109, 043028. 10.1103/physrevd.109.043028
54
KiniY.SalmiT.VinciguerraS.WattsA. L.BilousA.GallowayD. 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. 10.1093/mnras/stae2398
55
KöhlerH. S. (1976). Skyrme force and the mass formula. Nucl. Phys. A258 (2), 301–316. 10.1016/0375-9474(76)90008-7
56
KubisS.WojcikW.CastilloD. A.ZabariN. (2023). Relativistic mean-field model for the ultracompact low-mass neutron star HESS J1731-347. Phys. Rev. C108, 045803. 10.1103/physrevc.108.045803
57
KumarM.KumarS.ThakurV.KumarR.AgrawalB. K.DhimanS. K. (2023). CREX- and PREX-II-motivated relativistic interactions and their implications for the bulk properties of nuclear matter and neutron stars. Phys. Rev. C107, 055801. 10.1103/physrevc.107.055801
58
KumarM.KumarS.ThakurV.KumarR.AgrawalB. K.DhimanS. K. (2024). Trends of nuclear matter and neutron-star observables post CREX-PREX-II results within the relativistic mean field model. Phys. Rev. C110, 055802. 10.1103/physrevc.110.055802
59
KurkelaA.FragaE. S.Schaffner-BielichJ.VuorinenA. (2014). Constraining neutron star matter with quantum chromodynamics. Astrophys. J.789, 127. 10.1088/0004-637x/789/2/127
60
KurkelaA.RajagopalK.SteinhorstR. (2024). Astrophysical equation-of-state constraints on the color-superconducting gap. Phys. Rev. Lett.132, 262701. 10.1103/physrevlett.132.262701
61
KurkelaA.RomatschkeP.VuorinenA. (2010). Cold quark matter. Phys. Rev. D.81, 105021. 10.1103/physrevd.81.105021
62
Laskos-PatkosP.KoliogiannisP. S.MoustakidisCh. C. (2024). Hybrid stars in light of the HESS J1731-347 remnant and the PREX-II experiment. Phys. Rev. D.109, 063017. 10.1103/physrevd.109.063017
63
Laskos-PatkosP.LalazissisG. A.WangS.MengJ.RingP.MoustakidisCh. C. (2025). Speed of sound bounds and first-order phase transitions in compact stars. Phys. Rev. C111, 025801. 10.1103/physrevc.111.025801
64
Laskos-PatkosP.MoustakidisCh.C. (2025). XTE J1814-338: a potential hybrid star candidate. Phys. Rev. D.111, 063058. 10.1103/physrevd.111.063058
65
LattimerJ. M.PrakashM. (2001). Neutron star structure and the equation of state. Astrophys. J.550, 426–442. 10.1086/319702
66
LiJ. J.SedrakianA. (2023). Baryonic models of ultra-low-mass compact stars for the central compact object in HESS J1731-347. Phys. Lett.844B, 138062. 10.1016/j.physletb.2023.138062
67
LiJ. J.SedrakianA.AlfordM. (2020). Relativistic hybrid stars with sequential first-order phase transitions and heavy-baryon envelopes. Phys. Rev. D.101, 063022. 10.1103/physrevd.101.063022
68
LiJ. J.SedrakianA.AlfordM. (2023a). Relativistic hybrid stars with sequential first-order phase transitions in light of multimessenger constraints. Astrophys. J.944, 206. 10.3847/1538-4357/acb688
69
LiJ. J.SedrakianA.AlfordM. (2023b). Ultracompact hybrid stars consistent with multimessenger astrophysics. Phys. Rev. D.107, 023018. 10.1103/physrevd.107.023018
70
LiJ. J.SedrakianA.AlfordM. (2024). Hybrid star models in the light of new multimessenger data. Astrophys. J.967, 116. 10.3847/1538-4357/ad4295
71
LiJ. J.SedrakianA.AlfordM. (2025). Confronting new NICER mass-radius measurements with phase transition in dense matter and twin compact stars. JCAP02, 002. 10.1088/1475-7516/2025/02/002
72
LopesL. L.IssifuA. (2025). XTE J1814-338 as a dark matter admixed neutron star. Phys. Dark. Univ.48, 101922. 10.1016/j.dark.2025.101922
73
LugonesG.HorvathJ. E. (2002). Color-flavor locked strange matter. Phys. Rev. D.66, 074017. 10.1103/physrevd.66.074017
74
M. VeselskýPetousisV.KoliogiannisP.MoustakidisC.LejaJ. (2025). Simultaneous explanation of XTE J1814-338 and HESS J1731-347 objects using K− and condensates. Phys. Rev. D.111, L061308. 10.1103/physrevd.111.l061308
75
MarianiM.Ranea-SandovalI. F.LugonesG.OrsariaM. G. (2024). Could a slow stable hybrid star explain the central compact object in HESS J1731-347?Phys. Rev. D.110, 043026. 10.1103/physrevd.110.043026
76
MillerM. C.LambF. K.DittmannA. J.BogdanovS.ArzoumanianZ.GendreauK. 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. 10.3847/2041-8213/ab50c5
77
MiyatsuT.CheounM.-K.KimK.SaitoK. (2023). Can the PREX-2 and CREX results be understood by relativistic mean-field models with the astrophysical constraints?Phys. Lett. B843, 138013. 10.1016/j.physletb.2023.138013
78
NaseriM.BozzolaG.PaschalidisV. (2024). Exploring pathways to forming twin stars. Phys. Rev. D.110, 044037. 10.1103/physrevd.110.044037
79
NegreirosR.ZhangC.XuR. (2025). Rotational properties of inverted hybrid stars. Phys. Rev. D.111, 063026. 10.1103/physrevd.111.063026
80
OikonomouP. T.MoustakidisCh.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. 10.1103/physrevd.108.063010
81
OlesenM. L.MadsenJ. (1994). Nucleation of quark matter bubbles in neutron stars. Phys. Rev. D.49, 2698–2702. 10.1103/physrevd.49.2698
82
OppenheimerJ. R.VolkoffG. M. (1939). On massive neutron cores. Phys. Rev.55, 374–381. 10.1103/physrev.55.374
83
PalS.PodderS.ChaudhuriG. (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. ApJ983, 24. 10.3847/1538-4357/adbc6b
84
PereiraJ. P.FloresC. V.LugonesG. (2018). Phase transition effects on the dynamical stability of hybrid neutron stars. Astrophys. J.860, 12. 10.3847/1538-4357/aabfbf
85
QiL.ZhengS.ZhangJ.GeM.LiA.ZhangS. N.et al (2025). PSR j1231-1411 revisited: pulse profile analysis of X-ray observation. Astrophys. J.981, 99. 10.3847/1538-4357/adb42f
86
RatherI. A.MarquezK. D.BackesB. C.PanotopoulosG.LopesI. (2024). Radial oscillations of hybrid stars and neutron stars including delta baryons: the effect of a slow quark phase transition. JCAP 052024, 130. 10.1088/1475-7516/2024/05/130
87
RatherI. A.PanotopoulosG.LopesI. (2023). Quark models and radial oscillations: decoding the HESS J1731-347 compact object’s equation of state. Eur. Phys. J. C83, 1065. 10.1140/epjc/s10052-023-12223-1
88
ReedB. T.FattoyevF. J.HorowitzC. J.PiekarewiczJ. (2021). Implications of PREX-2 on the equation of state of neutron-rich matter. Phys. Rev. Lett.126, 172503. 10.1103/physrevlett.126.172503
89
ReedB. T.FattoyevF. J.HorowitzC. J.PiekarewiczJ. (2024). Density dependence of the symmetry energy in the post–PREX-CREX era. Phys. Rev. C109, 035803. 10.1103/physrevc.109.035803
90
ReinhardP. G.Roca-MazaX.NazarewiczW. (2023). Combined theoretical Analysis of the parity-violating Asymmetry for48Ca and208Pb. Phys. Rev. Lett.129, 232501. 10.1103/physrevlett.129.232501
91
RenJ.ZhangC. (2020). Quantum nucleation of up-down quark matter and astrophysical implications. Phys. Rev. D.102, 083003. 10.1103/physrevd.102.083003
92
RileyT. E.WattsA. L.BogdanovS.RayP. S.LudlamR. M.GuillotS.et al (2019). A NICER view of PSR J0030+0451: millisecond pulsar parameter estimation. Astrophys. J. Lett.887, L21. 10.3847/2041-8213/ab481c
93
RomaniR. G.KandelD.FilippenkoA. V.BrinkT. G.ZhengW. (2022). PSR j0952-0607: the fastest and heaviest known galactic neutron star. Astrophys. J. Lett.934, L17. 10.3847/2041-8213/ac8007
94
RoutarayP.DasH. C.Amrit PattnaikJ.KumarB. (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. 10.1142/s0218301324500526
95
SagunV.GiangrandiE.DietrichT.IvanytskyiO.NegreirosR.ProvidenciaC. (2023). What is the nature of the HESS j1731-347 compact object?Astrophys. J.958 (49), 49. 10.3847/1538-4357/acfc9e
96
SalinasM.PiekarewiczJ. (2024). Impact of tensor couplings with scalar mixing on covariant energy density functionals. Phys. Rev. C109, 045807. 10.1103/physrevc.109.045807
97
SalmiT.ChoudhuryD.KiniY.RileyT. E.VinciguerraS.WattsA. 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. 10.3847/1538-4357/ad5f1f
98
SalmiT.DenevaJ. S.RayP. S.WattsA. L.ChoudhuryD.KiniY.et al (2024a). A NICER view of PSR j1231-1411: a complex case. Astrophys. J.976, 58. 10.3847/1538-4357/ad81d2
99
SamsingJ.HotokezakaK. (2021). Populating the black hole mass gaps in stellar clusters: general relations and upper limits. ApJ923, 126. 10.3847/1538-4357/ac2b27
100
ShaoW. (2022). On the neutron star/black hole mass gap and black hole searches. Astron. Astrophys.22, 122002. 10.1088/1674-4527/ac995e
101
TagamiS.WakasaT.YahiroM. (2022). Slope parameters determined from CREX and PREX2. Results Phys.43, 106037. 10.1016/j.rinp.2022.106037
102
TewariS.ChatterjeeS.KumarD.MallickR. (2024). Analyzing the dense matter equation of states in the light of the compact object HESS J1731-347. Phys. Rev. D111, 103009 (2025). 10.1103/PhysRevD.111.103009
103
ThompsonT. A.KochanekC. S.StanekK. Z.BadenesC.PostR. S.JayasingheT.et al (2019). A noninteracting low-mass black hole–giant star binary system. Science366, 637–640. 10.1126/science.aau4005
104
TsaloukidisL.KoliogiannisP. S.Kanakis-PegiosA.MoustakidisCh.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. 10.1103/physrevd.107.023012
105
TypelS.OertelM.KlähnT.ChatterjeeD.DexheimerV.IshizukaC.et al (2022). CompOSE core team collaboration, CompOSE reference manual. Eur. Phys. J.58, 221. 10.1140/epja/s10050-022-00847-y
106
VinciguerraS.SalmiT.WattsA. L.ChoudhuryD.RileyT. E.RayP. S.et al (2024). An updated mass-radius analysis of the 2017-2018 NICER data set of PSR J0030+0451. Astrophys. J.961, 62. 10.3847/1538-4357/acfb83
107
WeberF. (2005). Strange quark matter and compact stars. Prog. Part. Nucl. Phys.54 (1), 193–288. 10.1016/j.ppnp.2004.07.001
108
WittenE. (1984). Cosmic separation of phases. Phys. Rev. D.30, 272–285. 10.1103/physrevd.30.272
109
YangS.-H.PiC.-M. (2024). Color-flavor locked strange stars admixed with mirror dark matter and the observations of compact stars, JCAP09 2024, 052.
110
YangS.-H.PiC.-M.WeberF. (2025). Strange stars admixed with mirror dark matter: confronting observations of XTE J1814-338. Phys. Rev. D.111, 043037. 10.1103/physrevd.111.043037
111
YangS.-H.PiC.-M.ZhengX.-P.WeberF. (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. 10.1103/physrevd.103.043012
112
YangS.-H.PiC.-M.ZhengX.-P.WeberF. (2023). Confronting strange stars with compact-star observations and new physics. Universe9, 202. 10.3390/universe9050202
113
ZhangC.LuoY.LiH. B.ShaoL.XuR. (2024). Radial and nonradial oscillations of inverted hybrid stars. Phys. Rev. D.109, 063020. 10.1103/physrevd.109.063020
114
ZhangC.RenJ. (2023). Hybrid stars may have an inverted structure. Phys. Rev. D.108, 063012. 10.1103/physrevd.108.063012
115
ZhangN.-B.LiB.-A. (2025). Impact of the nuclear equation of state on the formation of twin stars. Eur. Phys. J. A61, 31. 10.1140/epja/s10050-025-01497-6
116
ZhouT.HuangC. (2025). Hidden twin star solutions from an agnostic speed-of-sound model: confronting XTE j1814-338’s extreme compactness. arXiv:2504, 08662. 10.48550/arXiv.2504.08662
Summary
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
Volume
12 - 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
Updates
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, moustaki@auth.gr
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.