ORIGINAL RESEARCH article

Front. Astron. Space Sci., 18 December 2023

Sec. Cosmology

Volume 10 - 2023 | https://doi.org/10.3389/fspas.2023.1334642

Ultra-low-mass and small-radius white dwarfs made of heavy elements

  • 1. Center for Gravitation and Cosmology, College of Physical Science and Technology, Yangzhou University, Yangzhou, China

  • 2. School of Astronomy and Space Science, Nanjing University, Nanjing, China

  • 3. Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing, China

  • 4. School of Physics, Peking University, Beijing, China

  • 5. Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing, China

  • 6. National Astronomical Observatories, Chinese Academy of Sciences, Beijing, China

Abstract

Seven possible ultra-low-mass and small-radius white dwarfs have been recently identified, with masses ranging from ∼0.02 M to ∼0.08 M and radii ranging from ∼ 4,270 km to 10670 km. The mass–radius measurements of these white dwarfs pose challenges to traditional white dwarf models, assuming they are mostly made of nuclei lighter than 56Fe. In this work, we consider the possibility that those white dwarfs are made of heavier elements. Due to the small charge-to-mass ratios in heavy elements, the electron number density in white dwarf matter is effectively reduced, which reduces the pressure with additional contributions of lattice energy and electron polarization corrections. This consequently leads to white dwarfs with much smaller masses and radii, which coincide with the seven ultra-low-mass and small-radius white dwarfs. The mass of the most massive white dwarfs is effectively reduced and could possibly account for the sub-Chandrasekhar progenitors in underluminous Type Ia supernovae. The corresponding equation of state and matter contents of dense stellar matter with and without reaching the cold-catalyzed ground state are presented, which are obtained using the latest Atomic Mass Evaluation (AME 2020). Further observations are necessary to unveil the actual matter contents in those white dwarfs via, e.g., spectroscopy, asteroseismology, and the discoveries of other ultra-low-mass and small-radius white dwarfs.

1 Introduction

White dwarfs represent the final destiny of the vast majority of stars, which may reach a temperature of ∼100 eV and a density of g/cm3 in their centers (). The matter contents of typical white dwarfs are 12C and 16O, covered by a thin envelope of 4He (and 1H, if not burned entirely). If the progenitor star approaches the 10M limit, white dwarfs are thought to be made of 16O and 20Ne (); however, He white dwarfs are also possible if the progenitors are of very low mass and in a binary system (; ). As a white dwarf slowly cools down, a crystallized core will be formed, releasing latent heat that delays the cooling process (). Most white dwarfs are expected to go through at least one pulsation phase during their evolution, displaying periodic variations in their brightness that arise from global oscillations (; ). Additionally, various oscillation modes can be excited during the late inspiral or merger of white dwarf binaries, which may emit gravitational waves that are detectable for future space-borne gravitational wave detectors ().

Combined with the measurements of the distance and surface temperature Teff of a white dwarf (), its radius R can be inferred according to the observational flux. The surface gravity of a white dwarf can also be measured according to the gravitational redshift of the spectrum lines from the atmosphere (; ; ), which can be used to fix the mass with additional information on its radius. Throughout the available data on the masses and radii of white dwarfs in the Montreal White Dwarf Database (MWDD; ), as indicated in Table 1, seven ultra-low-mass and small-radius white dwarfs have been identified with masses ranging from ∼0.02 M to ∼0.08 M and radii ranging from ∼ 4,270 km to 10,670 km (), which deviate from the MR relation of typical white dwarfs. It should be noted that no spectral lines have been identified in those stars, and the surface gravity needs to be inferred theoretically (), which results in large uncertainties. Further confirmation of those white dwarfs may pose challenges to traditional white dwarf models, which are considered to be made of light elements such as 12C, 16O, 4He, and 20Ne. Consequently, traditional white dwarf models predict much larger radii than those indicated in Table 1. In addition to these anomalous white dwarfs, there are underluminous Type Ia supernovae whose progenitor masses appear to be well below 1.4 M (; ), which again challenges the large Chandrasekhar mass limit (∼1.4 M) predicted by traditional white dwarf models.

TABLE 1

MWDD IDMRTeff
MkmK
LSPM J0815 + 16330.082 ± 0.03113,563.23 ± 1,024.764,655 ± 35
LP 240-300.081 ± 0.01613,542.5 ± 626.64,680 ± 25
BD+20 5125B0.08 ± 0.03813,046.72 ± 1,124.184,395 ± 90
LP 462-120.054 ± 0.02411,999.23 ± 1,552.784,800 ± 20
WD J1257 + 54280.032 ± 0.0312,403.13 ± 3,561.257,485 ± 85
2MASS J13453297 + 42004370.031 ± 0.049,186.23 ± 3,405.514,270 ± 75
SDSS J085557.46 + 053524.50.02 ± 0.24514,688.29 ± 767.0710,670 ± 1,677

Masses, radii, and surface temperatures of seven ultra-low-mass and small-radius white dwarfs (; ; ).

To understand the physical origin of such low-mass and small-radius white dwarfs, possible candidates made of various exotic materials were considered. For example, proposed that they are, in fact, strange dwarfs comprised of a strange quark matter (SQM) core and a thick normal matter crust (). It was suggested that the intermittent fractional collapses of the crust induced by the refilling of materials accreted from its low-mass companion lead to repeating fast radio bursts (). Additionally, there are possibilities that these white dwarfs may be strangelet dwarfs () or ud quark matter (udQM) dwarfs (; ) comprised of strangelets or udQM nuggets immersed in a sea of electrons.

In this work, we consider the possibility that those low-mass and small-radius white dwarfs are made of heavy elements such as 56Fe, 62Ni, 108Pd, and 208Pb, which significantly reduces the mass and radius of a white dwarf in comparison with those made of light elements. In fact, there exist extensive observations of (DZ) white dwarfs contaminated by heavy elements, which were identified by the characteristic spectral lines (; ; ). It is, thus, reasonable to consider the possibility that white dwarfs can also be made entirely of heavy elements, where the shell of light elements is either stripped by its companion object in a binary system or by a nuclear explosion, e.g., in a thermonuclear electron-capture supernova event ().

The remainder of this paper is organized as follows. Section 2 presents the theoretical framework for obtaining the properties of white dwarf matter under various constraints on the mass numbers of nuclei, including either light or heavy elements. The obtained equation of state (EOS) and matter contents that minimize the energy density of cold white dwarf matter are presented in Section 3, while the corresponding white dwarf structures are investigated and compared with the low-mass and small-radius white dwarfs. Section 4 presents our conclusion.

2 Theoretical framework

For cold white dwarf matter comprised of crystallized nuclei immersed in a sea of electrons, the energy density can be divided into three parts (), i.e.,where ne is the average electron number density and α = 1/137.03599911 is the fine structure constant. The first term in Eq. 1 represents the energy density of nuclei, with MN(Z, A) being the mass of a nucleus with Z protons and A nucleons, which is determined bywhere MA(Z, A) is the measured atomic mass (AME 2020; ; ), me = 510998.95 eV is the electron mass, and eV is the electron binding energy (). The baryon number density is then given by nb = ne/fZ with the charge-to-mass ratio fZ = Z/A. The second term in Eq. 1 corresponds to the energy density of electrons, including exchange contributions. Ee is the energy density of free electron gas, which is given bywith . The third term in Eq. 1 represents the lattice energy density, including electron polarization corrections. For a body-centered cubic lattice, the Madelung constant is KM = −0.895929255682 (), and the function σ(Z) (; ) is given by

For nuclei to stably exist inside white dwarfs, they should be stable against electron capture and β-decay reactions, i.e.,

This indicates the following stability condition, i.e.,

with the energy density

E

fixed by Eq.

1

. It should be noted that the electron number density changes to

ne

(

Z

± 1)/

Z

as we vary the charge number of a nucleus from

Z

to

Z

± 1. A vast number of nuclei species fulfilling the stability condition (7) are obtained, which could all stably exist inside white dwarfs if there are no other decay channels. At a fixed baryon number density

nb

, we search for the nucleus that minimizes the energy density of white dwarf matter under three different considerations, i.e.,

  • 1. Massive white dwarfs comprised of light elements (A ≤ 16, 24, 28)

  • 2. Catalyzed ones with all possible nuclear species

  • 3. Ultra-low-mass and small-radius white dwarfs made of heavy elements (A ≥ 62, 108, 208)

Once the nucleus is fixed, the energy density can then be obtained using Eq. 1. A similar procedure is carried out in a vast density range with nb ≈ 10–13–10−4 fm−3. At larger densities, the nucleus becomes too neutron-rich, causing neutrons to start dripping out and forming a neutron gas, a phenomenon beyond the scope of the current study. According to basic thermodynamic relations, the pressure of white dwarf matter is then determined bywith .

At smaller densities with pressure P ≲ 0.4Pe ≈ 10–17 MeV/fm3, nuclei are not fully ionized as few electrons start to bind to them. In such cases, the EOS predicted by Eqs 1, 8 is no longer valid. We follow the treatment outlined by and employ the results obtained by , where the pressure of white dwarf matter now becomes

It should be noted that a dampening factor f(Z, ne) is introduced, where Eq. 9 represents the pressure of non-interacting electrons in the non-relativistic limit when f = 1. The exact value of f(Z, ne) is determined by interpolating the results presented in Figure 1 of the work of .

FIGURE 1

; ). For comparison, the nuclear species from the BPS EOS is presented (). The open circles indicate the central densities of the most massive white dwarfs, corresponding to the values in the sixth column of Table 2.

3 Results and discussion

In Figure 1, we present the proton number Z, mass number A, and charge-to-mass ratio fZ = Z/A of nuclei in white dwarf matter as functions of energy density, where various constraints on the mass numbers of nuclei are adopted with A ≤ 16, 24, 28 or A ≥ 62, 108, 208. Similar to the BPS model (), the catalyzed one indicated by the black solid curve is obtained by searching for all possible nuclear species without any restriction on A, which generally agrees with the BPS EOS at E ≲ 0.1 MeV/fm3. The latest data from the Atomic Mass Evaluation (AME 2020) have been adopted in our calculations (; ). It is found that the nuclear species remains unchanged at E ≲ 3 × 10−4 MeV/fm3 except for the catalyzed one, which covers most of the density range in white dwarfs. As indicated by the symbols in the bottom panel of Figure 1 and in the second column of Table 2, at small densities, the nuclei 16O, 24Mg, 28Si, 56Fe, 62Ni, 108Pd, and 208Pb minimize the energy density of white dwarf matter under various constraints on A. The corresponding charge-to-mass ratio fZ starts to decrease with A at A ≥ 62, which reduces the electron number density with ne = fZnb. As will be shown later, this will consequently reduce the pressure and make white dwarfs more compact. At larger densities, the nuclear species starts to change, which typically takes place at the center of the most massive white dwarfs (indicated by the open circles in Figure 1), except for the catalyzed one. If we further increase the density, the nuclear species continues to change, which decreases the charge-to-mass ratio fZ.

TABLE 2

CriterionAX0.03R0.03MmaxRmaxEcPcAXc
kmMkmeV/fm3eV/fm3
A ≤ 1616O21,3001.381,4224,790.24.28316O
A ≤ 2424Mg19,7671.332,201976.40.48924Mg
A ≤ 2828Si19,1231.302,522477.90.23828Si
Catalyzed56Fe15,1371.002,080830.30.38466Ni
A ≥ 6262Ni14,3001.002,064811.40.38366Ni
A ≥ 108108Pd11,7600.892,087601.50.238108Ru
A ≥ 208208Pb9,1300.751,3912,040.71.098208Hg

Summary of white dwarf properties obtained under different constraints on nuclear mass numbers. The radii (R0.03) of 0.03M white dwarfs and their matter contents (AX0.03) are indicated in the third and second columns. The masses (Mmax) and radii (Rmax) of the most massive white dwarfs are presented, along with the energy density Ec, pressure Pc, and nuclei species (AXc) at the center of these white dwarfs.

Based on the nuclear species indicated in Figure 1, the pressure of white dwarf matter can then be fixed using Eqs 8, 9. The corresponding EOSs under various constraints are then presented in Figure 2. At E ≳ 10–4 MeV/fm3, the EOSs of white dwarf matter generally coincide with each other, while there are slight variations due to the sudden changes in nuclear species causing mild first-order phase transitions. At E ≲ 10–4 MeV/fm3, the effects of chemical composition become evident, and the pressure is effectively reduced if white dwarf matter is made of heavy elements. This is attributed to the reduction in the charge-to-mass ratio fZ at A ≥ 62, which decreases the electron number density and, consequently, the pressure. The pressure reduction becomes more evident at smaller densities, where the variations in lattice energy density and electron polarization corrections become sizable even for cases with the same fZ value.

FIGURE 2

To better illustrate the effects of heavy elements on the EOSs of white dwarf matter, in Figure 3, we present the corresponding adiabatic index as a function of energy density, which is determined by

FIGURE 3

For electron gas in the non-relativistic limit, one expects Γ = 5/3, which turns into 4/3 in the extreme relativistic limit with a softer EOS. This is indeed the case for white dwarf matter at E ≳ 10–4 MeV/fm3, where the EOSs in Figure 2 generally coincide with each other with Γ = 4/3. Nevertheless, there are few exceptions during the first-order phase transitions with sudden changes in nuclear species, which reduces the adiabatic index to Γ = 0. At smaller densities, however, Γ easily exceeds the non-relativistic limit of 5/3. In particular, as density decreases, Γ quickly increases and becomes larger if white dwarf matter is made of heavy elements with a larger A. This can be attributed to the additional contributions of lattice energy and electron polarization corrections, which are sizable at small densities.

Based on the EOSs presented in Figure 2, the structures of white dwarfs are fixed by solving the TOV equation. where G = 6.707 × 10−45 MeV−2 is the gravitational constant. In Figure 4, we present the mass–radius relations of white dwarfs under various constraints on their matter contents. Typical white dwarfs with M > 0.5M () and seven ultra-low-mass and small-radius white dwarfs are indicated by the dots with the corresponding error bars (). It is evident that the typical white dwarfs are made of light elements with A ≤ 28, where the mass and radius could become slightly larger if lighter elements (12C, 4He, and 1H), magnetic field, and temperature effects are accounted for. However, compared with the seven ultra-low-mass and small-radius white dwarfs, the radii of these normal white dwarfs are too large. To resolve this, it was suggested that the low-mass and small-radius white dwarfs are made of exotic matter such as SQM or udQM. As illustrated by , a strange dwarf with a mass similar to that of a normal white dwarf could harbor an extremely dense SQM core, resulting in a smaller radius that coincides with that of the seven white dwarfs. For absolute stable SQM or udQM with sufficiently small surface tension, it was shown that strangelets or udQM nuggets of a certain size can be more stable than others (); consequently, the surface of a quark star fragments into a crystalline crust comprised of strangelets or udQM nuggets immersed in a sea of electrons (). Further investigations have revealed that the crustal material could also form strangelet () or udQM dwarfs (; ), where the radii are typically smaller than those of normal white dwarfs, and can, thus, accommodate the seven small-radius white dwarfs.

FIGURE 4

), while the seven dots below correspond to the ultra-low-mass and small-radius white dwarfs, namely, LSPM J0815 + 1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257 + 5428, 2MASS J13453297 + 4200437, and SDSS J085557.46 + 053524.5 ().

In this work, we resort to the more familiar explanation that the white dwarfs are made of heavier elements. As indicated in Figure 4, the masses and radii of white dwarfs decrease significantly if they are made of heavy elements, which coincide with the observed seven ultra-low-mass and small-radius white dwarfs (). In particular, to accommodate the extreme small radius of 2MASS J13453297 + 4200437, the white dwarf should be comprised of elements with mass numbers A ≥ 108. Similar cases are observed for other small-radius white dwarfs LSPM J0815 + 1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257 + 5428, and SDSS J085557.46 + 053524.5, whose radii are too small if they are made of light elements with A ≤ 28. This is attributed to the reduction in pressure at small densities if white dwarfs are made of heavy elements, causing them to become more compact than typical white dwarfs. The maximum mass of white dwarfs is also reduced, reaching 0.75 M if they are made of heavy elements, which could be the progenitor for underluminous Type Ia supernovae.

To show this explicitly, in Figure 5, we present the internal energy density profiles of 0.03M white dwarfs, where the horizontal axis corresponds to the total mass m(r) enclosed in a sphere of radius r in Eq. 12. The corresponding nuclear species that makes up the white dwarf and its radius are indicated in the second and third columns of Table 2. It is evident that the energy density of white dwarf matter increases with the nuclear mass number A, leading to more compact white dwarfs with smaller radii.

FIGURE 5

Finally, it is worth mentioning that we considered only ideal scenarios in which white dwarfs are made entirely of elements with fixed upper or lower limits on the mass number A, while in reality, we expect they are comprised of elements with various mass numbers and fractions. In that case, the EOS can be obtained as a combination of the EOSs indicated in Figure 2, connected at certain critical pressures, where the corresponding mass–radius relation lies between those indicated in Figure 4. Meanwhile, the possible existence of a strong magnetic field could alter significantly the mass–radius relations of white dwarfs (), while the temperature effects are expected to have a minor impact.

4 Conclusion

The recent observations of the seven possible ultra-low-mass and small-radius white dwarfs, namely, LSPM J0815 + 1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257 + 5428, 2MASS J13453297 + 4200437, and SDSS J085557.46 + 053524.5, have posed challenges to traditional white dwarf models assuming that they are mostly made of nuclei lighter than 56Fe. To resolve this, previous investigations suggest that they may be made of exotic matter such as SQM or udQM (; ; ; ). In this work, we consider the possibility that these white dwarfs are made of heavier elements, which effectively reduces the pressure and leads to white dwarfs with much smaller masses and radii. In such cases, we can accommodate the low-mass and small-radius white dwarfs without introducing any exotic matter. The maximum mass of white dwarfs is also reduced if they are made of heavier elements, which could possibly account for the sub-Chandrasekhar progenitors in underluminous Type Ia supernovae.

It should be noted that we assumed white dwarfs are made entirely of elements with fixed upper or lower limits on the mass numbers, where the predicted radii should be viewed as lower limits for white dwarfs containing lighter elements. Additional theoretical uncertainties are also expected under the influence of a strong magnetic field and high temperatures. These issues should be considered in our future study, where their impacts on various properties of white dwarfs aside from mass–radius relations need to be addressed. Meanwhile, the uncertainty in the mass and radius measurements of the seven anomalous white dwarfs is still large, which may originate from possible unresolved binary systems or extrapolation errors at small masses. Further observations are, thus, necessary to unveil the actual structures and matter contents in those white dwarfs via, e.g., their cooling processes (), characteristic spectral lines emitted by the heavy elements (; ; ), pulsation that arises from their global oscillations (; ), gravitational waves excited during the late inspiral or merger of white dwarf binaries (), and searching for other white dwarfs with similar low-mass and small-radius characteristics ().

Statements

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

C-JX: conceptualization, funding acquisition, investigation, and writing–original draft. Y-FH: conceptualization, funding acquisition, and writing–review and editing. H-BL: conceptualization, validation, and writing–review and editing. LS: funding acquisition, investigation, and writing–review and editing. R-XX: conceptualization, funding acquisition, and writing–review and editing.

Funding

The authors declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Grant Nos. 12275234, 12342027, 12233002, and 12041306), the National SKA Program of China (Grant Nos. 2020SKA0120300 and 2020SKA0120100), and the National Key R&D Program of China (Grant No. 2021YFA0718500).

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.

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

Summary

Keywords

white dwarfs, stars, equation of state, heavy elements, Type Ia supernovae

Citation

Xia C-J, Huang Y-F, Li H-B, Shao L and Xu R-X (2023) Ultra-low-mass and small-radius white dwarfs made of heavy elements. Front. Astron. Space Sci. 10:1334642. doi: 10.3389/fspas.2023.1334642

Received

07 November 2023

Accepted

27 November 2023

Published

18 December 2023

Volume

10 - 2023

Edited by

She-Sheng Xue, International Center for Relativistic Astrophysics, Italy

Reviewed by

Fridolin Weber, San Diego State University, United States

Xinjian Wen, Shanxi University, China

Updates

Copyright

*Correspondence: Cheng-Jun Xia, ; Yong-Feng Huang, ; Hong-Bo Li, ; Lijing Shao, ; Ren-Xin Xu,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics