1. Introduction
White dwarfs (WDs) are astrophysical objects that originate from the remnants of stars whose initial mass was below approximately
[
1]. These stars, after depleting their reserves of nuclear fuel, enter a phase during which their cores contract because nuclear reactions can no longer counteract gravitational forces, while their outer layers expand. This collapse is halted only when the electrons within the star core become degenerate, providing the necessary pressure to counteract further gravitational collapse. Depending on the mass of the progenitor, the stellar outcome can be different. Indeed, the nuclear fusion reactions that occur during the star’s evolution can result in the formation of different nuclei, ultimately influencing the nature of the resulting WDs, such as helium (He), carbon–oxygen (C-O), and oxygen–neon–magnesium (O-Ne-Mg) WDs. It is essential to note that the maximum mass that a WD can attain, referred to as the Chandrasekhar mass and calculated to be approximately
[
2], varies depending on the composition of the WD. In practice, the majority of observed WDs are of the C-O type.
In work by Glendenning et al. [
3,
4], it was proposed that WDs could possess an inner core composed of absolutely stable strange quark matter. This is a consequence of the Bodmer–Witten hypothesis [
5,
6]. What makes this idea even more interesting is that the presence of this stable strange quark matter core has the potential to make some of these compact objects stable, while the corresponding configuration without the strange quark core would be unstable.
These objects, named “strange dwarfs” (SDs), exhibit characteristics distinct from those of conventional WDs. Specifically, they can have different radii, masses, and astrophysical evolution. It was conjectured that SDs could form either by accumulating normal nuclear matter on the surface of a strange quark star (SQS) or by collecting clusters of strange quark matter, commonly referred to as “strangelets”, onto WDs.
Glendenning et al. [
3] studied the radial oscillations of SDs to assess their stability, finding that they can remain stable even if the density of the nuclear matter envelope surpasses the maximum density observed in typical WDs.
The question concerning the stability of SDs underwent a thorough reexamination in the work of Alford et al. [
7]. By studying the fundamental radial mode, they found those objects to be unstable under radial oscillations.
It was suggested that the previous works by Glendenning et al. [
3,
4] may have inadvertently misinterpreted their findings by confusing the second–lowest eigenmode with the lowest one. However, upon closer examination and analysis, it became evident that these two research studies were built upon different underlying hypotheses [
8]. Contrary to the initial belief that the two studies were grounded on the same assumptions, it became clear that they operated within slightly distinct theoretical frameworks, each with its own validity. This realization effectively solved the apparent contradiction between the results obtained by the two works. Specifically, Di Clemente et al. [
8] and Di Clemente [
9] focused on the boundary conditions between the nuclear matter and the quark core. The analysis is based on the formalism established in work by Pereira et al. [
10] and Di Clemente et al. [
11]. These studies provide insights into the boundary conditions that ought to be applied in the context of rapid (and slow) conversions between nuclear matter and quark matter and in the context of phase transitions in general. Crucially, the specific boundary conditions employed can exert a substantial influence on the eigenvalues governing radial oscillations, and, by extension, they can have a profound impact on the overall stability of the star. In addition, Di Clemente et al. [
8] also addressed the applicability of the stability criterion based on the analysis of the extrema of the MR curve [
12,
13]. However, in this specific case, a crucial refinement is added to this criterion, by emphasizing the need for explicit specification regarding whether the quark content of the star remains constant or undergoes changes during the radial oscillations.
In this review, based on our previous work [
8], and in part of the first author’s PhD thesis [
9], we build upon prior research by incorporating a comprehensive analysis of the equations of state (EoSs) relevant to SDs. This involves a detailed examination of the mass–radius relationship for these objects. Furthermore, we give an explicit mathematical framework that yields a formula for fixing the quark content within the core of an SD.
Additionally, we have broadened our discussion to include potential astrophysical signatures of SDs. This extension is crucial for observational astrophysics, as it may provide insights into identifying and verifying the existence of SDs through observable phenomena.
2. Equation of State
One important consideration in assessing the stability of SDs lies in the nature of their EoSs. Historically, when examining this aspect in previous works [
3,
4,
7], the EoS was formulated as follows:
In this expression,
represents the Baym–Pethick–Sutherland (BPS) EoS [
14], while
denotes an EoS characterizing strange quark matter, which can be based on the MIT bag model [
15]. The most important parameter in this formulation is the transition pressure
(corresponding to the energy density
), which is the pressure at the boundary between quark matter and nuclear matter.
As the EoS for nuclear matter, we use the BPS EoS, which represents the ideal “limit” for WDs in which all elements up to Fe have been produced. This EoS displays a Chandrasekhar mass of approximately
, a value lower than the typical Chandrasekhar mass for C-O WDs or O-Mg-Ne WDs (
and
). It has been pointed out in work by Benvenuto and Althaus [
16] that the use of the BPS EoS is not realistic for WDs. Nevertheless, as will be explained in
Section 5.1, it represents the limit in compactness for WDs since it provides the smallest radius for a given mass.
In our work, we use a fit of the BPS EoS in order to avoid artifacts due to the numerical differentiation of a piecewise interpolation. The form of the fit of the BPS EoS is as follows:
where the function
f is as follows:
Note that the inclusion of decimal digits is essential for achieving adequate precision across a broad spectrum of pressures and densities. The fit ranges in energy density from ∼7 g/cm
3 to ∼
g/cm
3 and it is visible in
Figure 1. The percentage error is normally smaller than about
, except for a small region at densities of about 20 g/cm
3 where the error reaches about
.
When dealing with SDs, it is important to understand that one can use any value for as long as it is less than g/cm3.
Unlike regular WDs, where one typically only needs to specify the central pressure
to define a star’s configuration when solving the TOV equation [
17], SDs require two parameters. As clear from Equation (
1), the first parameter is the transition pressure
, which represents the pressure at the interface where the quark core meets the outer nuclear matter envelope, the second one is indeed the central pressure as in the normal case
.
What allows the formation of SDs is the existence of the Coulomb barrier that separates the outer nuclear matter from the inner core of quark matter. This separation occurs when the maximum density of nuclear matter is lower than the neutron drip density . Beyond this density, free neutrons start to appear. Importantly, since they are not subject to the constraints of the Coulomb barrier, they can readily penetrate the core of quark matter. Upon entering the core, they are absorbed, which leads to the deconfinement of their constituent quarks.
Given that for SDs, the solutions of the TOV depend on two parameters, a question arises about the suitability of choosing the pair of parameters
for characterizing these configurations. Choosing a value of
does not account for the fact that below the neutron drip density, the baryonic content of the core remains constant, despite changing the central pressure
. The studies by Vartanyan et al. [
18] and Vartanyan et al. [
19] discuss the case in which nuclear matter cannot transition into quark matter, allowing for the definition of sequences of SDs with the same quark content in the core, which we define as
. Consequently, one can solve the TOV equation with an alternative parameter pair, namely
.
The quark baryon number, represented as
, can be expressed as follows:
Here, is the baryon density within the quark core. It is important to note that these two parameter choices are not interchangeable. If one opts to keep constant while varying , this leads to changes in , implying that hadrons can deconfine into quarks because the change in the core is encoded in the fact that the central pressure is changing by fixing the external pressure of the core (transition pressure). Conversely, when is maintained at a constant value, one necessitates an increase in with higher values of . Notice that depends on and in the first case, and on and in the second case.
In order to correctly consider an SD EoS at its equilibrium, we want to choose the parameter pair . is a function of and ; therefore, we need to find the inverse relation that, given a choice of , gives back the value of .
In our analysis, since the core of the system is relatively small, we can reasonably approximate its impact on the gravitational force by using Newtonian physics in a first–order approximation and we will later add a general relativistic correction.
For the EoS governing a small quantity of quark matter within the core, we utilize a parametric expression, as follows:
Here, the variable
a serves as a multiplicative parameter that encompasses various factors, including the bag constant. The term
corresponds to the Witten density, defined as
. When considering the Newtonian hydrostatic equilibrium and Equation (
5), we obtain the following:
Upon integrating both sides and combining all constants into a single parameter, denoted as
K, we obtain the following equation:
In this equation,
represents the radius of the core, and
denotes the energy density at the center of the system. The solution to both sides of this equation leads us to the following expression:
This equation shows how the central energy density is connected to the Witten density , considering the parameter K and the core radius .
We can modify Equation (
8) to replace the core radius
with its baryon content, denoted as
, since
. Additionally, considering that our core experiences slight compression due to the surrounding nuclear matter, we can replace
with the effective energy density specific to the quark core, denoted as
, which is slightly greater than
. Indeed, the transition density at the boundary of the quark core satisfies
. The new form of the equation reads as follows:
Now, it is reasonable to incorporate some higher–order corrections into Equation (
9) to account for general relativistic effects:
In this modified equation, we introduce two numerical parameters, and , to account for these higher–order relativistic effects. These parameters are determined numerically and they depend on the quark matter EoS.
For the quark matter, we utilize the following thermodynamic potential [
20,
21,
22]:
Here,
is a parametrization of perturbative QCD corrections, the gap parameter
= 80 MeV, the strange quark mass
= 120 MeV, and the bag constant
MeV
4, as in work by Bombaci et al. [
23]. The specific parameter choices allow us to reach maximum masses for an SQS of about 2.6
.
It is important to note that the parameters
and
are primarily influenced by the bag constant, and any reasonable adjustments made to the other parameters in Equation (
11) have negligible impacts on
and
.
Mass–Radius Relation
The relationship between mass and radius exhibits notable differences depending on whether we consider the parameter pairs or . When constructing a mass–radius (MR) diagram, it is essential to vary one parameter, typically the central pressure (or central energy density), while keeping the other parameters constant.
If we opt to fix the transition pressure
, we are essentially exploring configurations where quark matter consistently appears at the same energy density threshold. We begin with a configuration, where
, denoted as point
b (or
and
depending on the transition energy density), as illustrated in
Figure 2. Then, we progress along the bottom branch in a clockwise direction, increasing the central pressure. Point
a indicates the Chandrasekhar mass of the WD, while
c is the maximum mass for SDs with
. The mechanically stable configurations are the ones to the right of
c.
In
Figure 2, it becomes clear how the final point of the MR of the BPS EoS (WD), defined as the point at which
, joins with the curve obtained by fixing
. Indeed, the SD curves, constructed by choosing
as the transition pressure, join exactly with the WDs’ MR at the point where, for the BPS,
(or
). When the transition density is relatively low, the point where the curves join falls before reaching the WDs’ maximum mass on the MR diagram (point
a in
Figure 2). This is evident in
Figure 3, where a low energy transition density (
g/cm
3) is represented by the dashed black curve, which has its maximum allowed mass at approximately 0.5
; this is where it joins the WD curve, which is not shown in the figure.
On the other hand, the MR relation for SDs tends to converge with that of an SQS when considering small radii. In particular, if the value of
is significantly smaller than the neutron drip density, it implies that there is insufficient matter above the quark core (which now constitutes the majority of the star) to exert significant compression on it. For the high values of
, specifically the neutron drip density, the radii of SDs are slightly smaller compared to those of an SQS. This occurs because there is a broader range of pressure that nuclear matter must cover in these cases. This behavior is shown in
Figure 4, with detailed representations of the maxima in
Figure 5 and the low–mass region in
Figure 6.
When we fix , we implicitly let vary. Conversely, when we fix the baryon content of the core (), it is the transition density that changes along the diagram.
When we establish a fixed value for
, the corresponding configuration contains a specific quantity of quarks in its core. If we start from a point at which
, there is no nuclear matter situated above the core to exert compression. In other words, it corresponds to the extreme point on the left side of
Figure 3 (red dashed curve) and satisfies Equation (
8).
As we increase the central energy density (which means adding nuclear matter on top of the quark core), we move in a counter-clockwise direction on the MR diagram. During this progression, we intersect curves in
Figure 3 that correspond to increasing values of
. This means that a curve representing a constant
is comprised of configurations with varying
. The initial point on this curve has
(or equivalently
), while the final point corresponds to
.
In
Figure 7, we can observe specific behavior where, if the value of
is too large, the condition
is achieved at relatively small radii. The extreme points belong to the curve representing the highest possible density of nuclear matter within the star, which corresponds to the neutron drip density (the blue dashed one).
4. SD Collapse
The transition from ordinary matter to strange quark matter can release huge amounts of energy and it can be associated with violent phenomena such as gamma–ray bursts [
28,
29,
30] and extremely energetic supernovae [
31]. In this review, we will only discuss the possible transition of ordinary matter in an SD to strange quark matter and we will discuss its possible signatures. In a binary system where a WD orbits a main sequence star, mass transfer occurs as the WD accretes material from its companion. This typical scenario culminates in a type Ia supernova (SN) event. However, a different outcome known as accretion–induced collapse (AIC) is theoretically possible. It is important to note that while the concept of AIC has been explored, actual observations of such events are notably absent [
32]. This absence can be attributed to the substantial difference in timescales between the collapse process and the nuclear reactions responsible for igniting the WD deflagration.
The presence of the strange quark matter core in SDs becomes important when the object faces significant perturbations, like in the early stages of a type Ia supernova (SN) event. Specifically, if the quark matter core is sufficiently large, it can potentially cause the object to collapse rather than undergo the usual deflagration process (see the illustration in
Figure 12). The difficulties of achieving an AIC arise from the fact that nuclear reactions occur when the star is near the Chandrasekhar limit. This phase, characterized by marginal mechanical stability (
), leads to the star’s disruption before AIC can take place [
33].
The mechanical instability of SDs is strictly related to the rapid conversion of hadrons into quarks. This process is a crucial mechanism that allows the star to undergo a collapse by becoming mechanically unstable.
As long as
, the object remains mechanically stable. However, the system becomes unstable if a fluctuation causes the density to exceed
in a small region near the core or if free neutrons are produced and fall into the quark core. In this analysis, temperature can play an important role. From Haensel et al. [
34] and Hempel and Schaffner-Bielich [
35], one can notice that at temperatures exceeding about 0.5 MeV, a significant fraction of free neutrons appears, already at densities of the order of
g/cm
3. To assess this instability, we calculate the fundamental eigenvalue of a star at the Chandrasekhar limit, assuming
remains constant. While for a slow transition,
, in the case of a rapid transition, for the same point in the MR diagram,
becomes significantly negative. It is important to remark that each point at constant
corresponds to a point at constant
; therefore, one can go from a situation in which the transition is physically slow to a situation in which the star’s internal boundary is in a rapid transition regime and the baryon content of the core is not constant anymore. In any case, it is illogical to apply a slow transition scenario when
remains constant or to employ a fast transition scenario when
is held constant because
cannot increase in that case (and, therefore, the transition must not occur). However, exceptions may arise when the star is in close proximity to the Chandrasekhar mass. In such situations, perturbations could potentially drive a small region of the star, located near the strange core, to exceed the neutron drip density or generate some free neutrons. This, in turn, could trigger the phase transition, going from having a constant
to a situation in which
increases because of the neutron flux.
Figure 13 displays the e–folding time, defined as
. When
, the e–folding, which is the typical collapse timescale, falls below 1 second, implying that the collapse can be more rapid than the full development of a deflagration, which is of the order of several seconds [
36]. It can be argued that the amount of quark matter that triggers the collapse depends on the EoS. Therefore, we calculated the e–folding time for the set of parameters for the quark EoS presented in work by Bombaci et al. [
23]. The results remained consistent with our previous findings. In the same figure, the maximum density reached by the nuclear matter component at the boundary (
) is also shown. From the behavior of
, it is possible to determine when the static structure of a 1
SD remains similar to the one of a WD. When
, the structure of the star changes, and the boundary density
deviates from ∼
g/cm
3, which is the typical central density of a WD at the Chandrasekhar mass. This suggests that the presence of the quark core does not influence the static properties of the star unless the value of
is large enough to exert a noticeable gravitational pull.
An essential query regarding SDs pertains concerns how they accumulate the strange quark matter in their cores. The most straightforward explanation lies in the idea that WDs gradually accumulate strangelets over their lifespan. This idea is linked to the possibility that dark matter is made, at least in part, of strangelets [
6]. In a few papers, it has been shown that this scenario is compatible with the most recent data from cosmology and astrophysics [
37,
38,
39,
40]. In work by Di Clemente et al. [
41], we have shown that an astrophysical path leading to the formation of subsolar compact stars in electron-capture supernovae can be based on the hypothesis of dark matter made of strangelets. The existence of subsolar-mass compact objects has been suggested in work by Doroshenko et al. [
42] and it cannot be explained by solely considering standard equations of state [
43].
Another mechanism to produce subsolar-mass compact objects is based on AIC. If AIC takes place in an SD instead of a WD, an SQS is produced instead of an NS. Since this phenomenon is very energetic and the collapsed object is more bound than NSs, the final object can be a subsolar-mass km-sized compact star [
8,
44].