Free Access
Issue
A&A
Volume 640, August 2020
Article Number L11
Number of page(s) 4
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202038879
Published online 06 August 2020

© ESO 2020

1. Introduction

As a white dwarf cools down, the fully ionized plasma that makes up its core eventually becomes so correlated that a first-order phase transition occurs, leading to the formation of a solid core. Predicted more than 50 years ago (van Horn 1968), the unequivocal observational signature of this phenomenon has recently been brought to light using data from the Gaia satellite (Tremblay et al. 2019). This signature comes in the form of a pile-up of objects in the cooling sequence of evolving white dwarfs (see also Bergeron et al. 2019). The core crystallization of white dwarfs is accompanied by the release of latent heat and gravitational energy from the change in the C/O abundance profile (Mochkovitch 1983; García-Berro et al. 1988; Isern et al. 1997; Fontaine et al. 2001; Althaus et al. 2010). These two phenomena temporarily slow down the evolution of solidifying white dwarfs, leading to a pile-up of objects at the luminosities where the phase transition takes place.

Detecting this pile-up for normal-mass white dwarfs (∼0.6 M) is a delicate exercise as crystallization occurs at the same time as convective coupling, another important event in the evolution of a white dwarf. Convective coupling refers to the contact of the superficial convection zone with the degenerate and highly conductive interior. This temporarily slows down the evolution of the white dwarf (Fontaine et al. 2001, Fig. 5), thereby masking the delay caused by crystallization. For more massive objects, however, crystallization starts well before convective coupling, which allowed Tremblay et al. (2019) to unambiguously attribute to core crystallization the pile-up structure detected in the luminosity function of DA white dwarfs with masses between 0.9 and 1.1 M.

Using white dwarf population simulations, Tremblay et al. (2019) showed that the observed pile-up of massive white dwarfs is roughly consistent with the predictions of theoretical evolution sequences and that both latent heat release and C/O phase separation are needed to explain the observations1. However, discrepancies between the models and the observations were noticed: (1) the crystallization bump is predicted to start too early and (2) the amplitude of the pile-up is significantly underestimated (see also Kilic et al. 2020). These problems are a source of concern for the use of white dwarfs as cosmochronometers (Winget et al. 1987; Oswalt et al. 1996; Fontaine et al. 2001; Kalirai 2012; Hansen et al. 2013; Tremblay et al. 2014; Kilic et al. 2017; Isern 2019; Fantin et al. 2019). Core crystallization is a significant event in the evolution of a white dwarf and an accurate description of this phenomenon is needed to generate reliable theoretical cooling sequences.

In this Letter we present a new C/O phase diagram aimed at improving the modeling of the pile-up structure discovered by Tremblay et al. (2019). We describe our phase diagram in Sect. 2 and our white dwarf population simulations in Sect. 3, where we show the impact of our improved description of the phase transition on the luminosity function of massive hydrogen-atmosphere white dwarfs.

2. C/O phase diagram

The exact shape of the C/O phase diagram is crucial in determining the impact of core crystallization on white dwarf cooling (e.g., Althaus et al. 2012). The position of the liquidus dictates the temperature at the liquid–solid transition, and the separation between the liquidus and the solidus, ΔxO, determines the importance of O sedimentation during crystallization. Several versions of the C/O phase diagram already exist in the literature, each with its particular limitations. Early calculations used density-functional methods (Segretain & Chabrier 1993), which are intrinsically more approximate than modern simulation techniques. Other studies relied on analytic fits to Monte Carlo (MC) simulations (Ogata et al. 1993; Medin & Cumming 2010). This approach can be delicate due to the use of approximate linear mixing rules and to the sensitivity of the phase diagram on the somewhat arbitrary choices that are made to construct the analytic functions used to interpolate the MC data (DeWitt et al. 1996). Finally, molecular dynamics (MD) simulations were used to directly simulate the phase transition (Horowitz et al. 2010). Challenges associated with those MD methods include the existence of finite-size effects and their extreme computational cost, which prohibits the detailed well-sampled mapping of the phase diagram that is needed for white dwarf modeling.

We have computed a new accurate phase diagram by adapting the Gibbs–Duhem integration technique coupled to MC simulations (Kofke 1993a,b) to plasmas. This advanced approach, previously only used for mixtures of neutral particles, was specifically designed to address the limitations of other phase diagram mapping techniques outlined above. It consists in calculating the phase diagram by integrating at constant pressure the Clapeyron equation dT/dξ along the liquid–solid coexistence curve in the temperature (T) – fugacity fraction (ξ) space. The Clapeyron equation depends on the enthalpies and the compositions of the liquid and solid phases, which we calculate with Monte Carlo simulations in the semi-grand canonical ensemble: at constant pressure, temperature, number of ions, and fugacity fraction. The details of our adaptation of the method to plasmas and analytical fits to our C/O phase diagram for an accurate implementation in white dwarf evolution codes will be published elsewhere.

Our new C/O phase diagram is shown in Fig. 1, where we compare it to previous calculations. Qualitatively, it is close to those of Medin & Cumming (2010) and Horowitz et al. (2010), with a similar azeotropic form. However, the quantitative differences are important in the context of white dwarf modeling. For instance, the narrow separation between the liquidus and the solidus in the phase diagram of Horowitz et al. (2010) would lead to a slight underestimation of O sedimentation. We are confident that our new phase diagram is the most accurate as it is free of many of the limitations and approximations of previous studies. In particular, (1) we directly integrate the MC simulations along the coexistence curve, removing the need for arbitrary fit functions to a sparse set of simulations; (2) the relatively low cost of MC simulations allows a very fine sampling of the diagram; (3) finite-size effects are easily mitigated since there is no need to simulate a liquid–solid interface as with MD methods; (4) all calculations are performed at constant pressure and not at constant volume as in all other approaches (phase transitions occur at constant pressure); (5) the relativistic electron jellium is explicitly included in the MC simulations; (6) screening of the ion–ion interactions by relativistic electrons is accounted for; and (7) the numerical precision of the MC simulations and Gibbs–Duhem integration is revealed by the smoothness of the final coexistence curve and the recovery of the exact melting temperature at the end point of the integration at xO = 1 (see Fig. 1). Finally, our approach provides a full description of the thermodynamics of the phase transition.

thumbnail Fig. 1.

C/O phase diagram. The liquidus (above which the plasma is entirely liquid) and the solidus (below which it is entirely solid) are both shown. The horizontal axis gives the number fraction of O, xO = NO/(NC + NO), and the vertical axis corresponds to the ratio of the temperature to the melting temperature of a pure carbon plasma. The results from this work (in red) are compared to those of Segretain & Chabrier (1993), Medin & Cumming (2010), and Horowitz et al. (2010). The region where an azeotrope is predicted is enlarged in the upper left corner.

3. Pile-up in the cooling sequence

We have computed new evolution sequences for massive hydrogen-atmosphere white dwarfs in order to test our updated phase diagram against the observed core crystallization pile-up. To do so we used STELUM, the Montreal white dwarf evolution code (Brassard & Fontaine 2018). The constitutive physics is identical to that described in Fontaine et al. (2001), except that (1) we now use the Cassisi et al. (2007) conductive opacities, which we correct in the moderately coupled and moderately degenerate regime (Blouin et al. 2020) following the new theory of Shaffer & Starrett (2020); (2) the plasma coupling parameter at the phase transition is given by our new C/O phase diagram; (3) the release of gravitational energy due to O sedimentation is implemented following Isern et al. (1997, 2000); and (4) diffusion of 22Ne in the liquid phase is included. We assume an initially chemically homogeneous core with X(C) = X(O) = 0.49 and X(22Ne) = 0.02 (consistent with the results of Salaris et al. 2010 for M ≈ 1 M) and we use an envelope stratification given by MH/M = 10−4 and MHe/M = 10−2 (the canonical values for DA stars). We find that the cooling delay imposed by C/O phase separation is of 1.0 Gyr at log L/L = −4.5 for our 0.9 M sequence. This is close to the value obtained by Althaus et al. (2012, Fig. 5) using the Horowitz et al. (2010) phase diagram, which is unsurprising given the similarity between the two phase diagrams (Fig. 1).

In order to compare our cooling sequences to the luminosity function given in Tremblay et al. (2019), we developed our own MC population synthesis code. We use the initial mass function of Salpeter (1955), main-sequence lifetimes from Hurley et al. (2000), the initial–final mass relation of Cummings et al. (2018), and synthetic photometry from state-of-the-art atmosphere models2 (Bergeron et al. 1995; Kowalski & Saumon 2006; Tremblay et al. 2011; Blouin et al. 2018). Our slightly different choice of main-sequence lifetimes and initial–final mass relation from those of Tremblay et al. (2019) has no effect on the results presented below.

Figure 2 compares our theoretical luminosity function for hydrogen-atmosphere white dwarfs between 0.9 and 1.1 M (black solid line) to the data of Tremblay et al. (2019, in red). We assume a constant stellar formation rate and a 10 Gyr age for the Galactic disk. The new C/O phase diagram leads to a luminosity function that is very close to that obtained by Tremblay et al. (2019). The magnitude of the crystallization pile-up is very similar, and the fit to the low-luminosity cutoff is nearly identical. We note that the corrections to the Cassisi et al. (2007) opacities in the regime of moderate Coulomb coupling and partial electron degeneracy (Blouin et al. 2020), which significantly increase the conductivity of the H envelope, were crucial to obtaining a cutoff close to the observations. Tremblay et al. (2019) were able to reproduce the cutoff without these corrections because they were relying on older calculations (Hubbard & Lampe 1969; Itoh et al. 1983; Mitake et al. 1984) that predict higher conductivities than Cassisi et al. (2007) in the core and in the He envelope (Salaris et al. 2013).

thumbnail Fig. 2.

Lower panel: luminosity functions for hydrogen-atmosphere white dwarfs with masses between 0.9 and 1.1 M. Observational data from Tremblay et al. (2019) are shown in red and their population synthesis is represented by the dotted curve. Results from our own simulations are shown as dashed and solid lines (with and without an attempt to simulate the effect of 22Ne phase separation, respectively). The normalization is arbitrary, and the Gaia G magnitude limit from Gentile Fusillo et al. (2019) is used. Upper panel: evolution of the fraction of the core that is crystallized for the 1.1, 1.0, and 0.9 M sequences.

One notable difference between our luminosity function and that of Tremblay et al. (2019) is that the bump associated with crystallization starts later in our simulations (log L/L ≈ −2.6 instead of −2.3), bringing the theoretical luminosity function closer to the data. This shift occurs because the Coulomb coupling parameter at the liquid–solid transition for an equimass C/O mixture is Γ ≈ 215 according to our phase diagram, while Tremblay et al. (2019) used the one-component plasma result of Γ = 175 (Potekhin & Chabrier 2000). We note that this problem does not affect evolution codes that already include a detailed treatment of crystallization (e.g., Althaus et al. 2012), but a comparison with only Tremblay et al. (2019) is possible at this point since no other studies have yet tried to model the Gaia crystallization pile-up.

We performed the most accurate calculation of the classical C/O phase diagram to date. Because it is in good agreement with previous calculations, we can conclude that the physics of the solidification of the C/O plasma in white dwarfs is well understood. This implies that the excess in the luminosity function at log L/L ≈ −3.3 above the prediction from C/O crystallization is likely due to a separate mechanism. The gravitational settling of 22Ne in the C/O liquid phase is another source of cooling delay in white dwarfs (Bildsten & Hall 2001; García-Berro & Oswalt 2016; Camisassa et al. 2016). It is already included in our sequences, but we note that its importance may be underestimated given the remaining uncertainties on the initial X(22Ne) profile and on the diffusion coefficients (e.g., Cheng et al. 2019). However, we checked that arbitrarily increasing the importance of 22Ne settling does not lead to a more prominent crystallization bump, and worsens the overall agreement with the observations by reducing the number of objects that have had the time to evolve to lower luminosities. At log L/L ≈ −3.3 a significant portion of the core has already solidified (see upper panel of Fig. 2), meaning that most of the 22Ne diffusion has already stopped. Therefore, any change in the treatment of 22Ne settling, whether from the initial X(22Ne) profile or the diffusion coefficients, is unlikely to solve the discrepancy at log L/L ≈ −3.3.

Another possibly important cooling delay may arise from the phase separation of 22Ne during crystallization (Isern et al. 1991; Althaus et al. 2010). Our current best understanding is that at the small 22Ne concentrations typical of C/O white dwarfs (∼1% by number), the presence of 22Ne should not affect the phase diagram, except near the azeotropic point of the C/O/Ne phase diagram. Thus, the crystallization of the C/O core initially proceeds as in the case without 22Ne with no redistribution of neon ions between the solid and liquid phases. After a significant fraction of the core has crystallized, the temperature approaches the azeotropic point and the existing calculations indicate that the liquid phase is enriched in 22Ne relative to the solid (Segretain 1996; García-Berro et al. 2008). The 22Ne-poor solid is lighter than the surrounding liquid and floats upward where it eventually melts. This gradually displaces the 22Ne-rich liquid downward toward the solid–liquid interface until the azeotropic composition is reached, thereby releasing a considerable amount of gravitational energy. Given our very limited knowledge of the ternary C/O/Ne phase diagram (Segretain 1996; Hughto et al. 2012), this effect cannot be quantitatively implemented in our evolution models. However, we note that our current understanding of 22Ne phase separation is remarkably consistent with the missing cooling delay. In Fig. 2 we show the luminosity function obtained by adding an artificial 0.6 Gyr delay when 60% of the core is crystallized. These parameters are entirely consistent with those found in preliminary studies (Segretain 1996; García-Berro et al. 2008) and yield an excellent fit to the crystallization pile-up3. Based on the current (albeit limited) knowledge of the C/O/Ne phase diagram, we propose that the phase separation of 22Ne in the advanced stage of crystallization significantly contributes to the pile-up in the luminosity function of 0.9−1.1 M white dwarfs (Fig. 2).

Finally, we speculate that the cooling anomaly for very massive white dwarfs (1.08−1.23 M) identified by Cheng et al. (2019), where roughly 6% of objects are affected by an unexplained ≈8 Gyr cooling delay, may also be at least partially explained by 22Ne phase separation rather than only by 22Ne diffusion in the liquid phase, as originally suggested. The energy source responsible for the unexplained cooling delay has an effect that is highly peaked on the crystallization sequence. We found that such a peaked effect is unlikely to occur from simple diffusion alone (which is inhibited by crystallization), while it can realistically be expected from 22Ne phase separation. A 0.6 Gyr delay is too small to explain the findings of Cheng et al. (2019), but this delay could be much more important if the initial 22Ne concentration were higher. A significant fraction of the massive objects of the Cheng et al. (2019) sample must come from double white dwarf mergers and, interestingly, additional 22Ne is expected to be formed during merger events (Staff et al. 2012). This would mean that those objects have a higher 22Ne abundance than the usual X(22Ne) = 0.01−0.02 (and a different distribution throughout the core), possibly leading to a longer cooling delay due to the phase separation of neon during crystallization.

By removing any remaining uncertainties on the classical C/O phase diagram, we have shown that the pile-up detected in the Gaia cooling sequence cannot be explained by latent heat release and O sedimentation alone. 22Ne phase separation appears to play a crucial role in the formation of the excess of massive white dwarfs observed at log L/L ≈ −3.3. Our results highlight the need for a complete and accurate ternary C/O/Ne phase diagram to establish quantitatively the importance of 22Ne phase separation in white dwarf evolution. We plan to generalize our newly developed Gibbs–Duhem integration method to three-component mixtures to address this problem.


1

Earlier results already supported the occurrence of O sedimentation in white dwarfs (García-Berro et al. 2010).

3

The additional cooling delay from 22Ne phase separation worsens the fit to the low-luminosity cutoff, but this may simply be due to the unrealistic assumption that 22Ne phase separation has the same importance for all stars in our simulation. In particular, the old stars that form the cutoff of the luminosity function likely contain less 22Ne than the younger ones that form the crystallization bump since the 22Ne abundance of a white dwarf increases for higher metallicity progenitors. Alternatively, this mismatch could be due to our assumption on the age of the Galactic disk.

Acknowledgments

Research presented in this article was supported by the Laboratory Directed Research and Development program of Los Alamos National Laboratory under project number 20190624PRD2. This work was performed under the auspices of the U.S. Department of Energy under Contract No. 89233218CNA000001.

References

  1. Althaus, L. G., Córsico, A. H., Isern, J., & García-Berro, E. 2010, A&ARv, 18, 471 [NASA ADS] [CrossRef] [Google Scholar]
  2. Althaus, L. G., García-Berro, E., Isern, J., Córsico, A. H., & Miller Bertolami, M. M. 2012, A&A, 537, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bergeron, P., Wesemael, F., & Beauchamp, A. 1995, PASP, 107, 1047 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bergeron, P., Dufour, P., Fontaine, G., et al. 2019, ApJ, 876, 67 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bildsten, L., & Hall, D. M. 2001, ApJ, 549, L219 [NASA ADS] [CrossRef] [Google Scholar]
  6. Blouin, S., Dufour, P., & Allard, N. F. 2018, ApJ, 863, 184 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blouin, S., Shaffer, N. R., Saumon, D., & Starrett, C. E. 2020, ApJ, submitted [arXiv:2006.16390] [Google Scholar]
  8. Brassard, P., & Fontaine, G. 2018, https://doi.org/10.5281/zenodo.1755253 [Google Scholar]
  9. Camisassa, M. E., Althaus, L. G., Córsico, A. H., et al. 2016, ApJ, 823, 158 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cassisi, S., Potekhin, A. Y., Pietrinferni, A., Catelan, M., & Salaris, M. 2007, ApJ, 661, 1094 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cheng, S., Cummings, J. D., & Ménard, B. 2019, ApJ, 886, 100 [NASA ADS] [CrossRef] [Google Scholar]
  12. Cummings, J. D., Kalirai, J. S., Tremblay, P. E., Ramirez-Ruiz, E., & Choi, J. 2018, ApJ, 866, 21 [NASA ADS] [CrossRef] [Google Scholar]
  13. DeWitt, H., Slattery, W., & Chabrier, G. 1996, Phys. B, 228, 21 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fantin, N. J., Côté, P., McConnachie, A. W., et al. 2019, ApJ, 887, 148 [CrossRef] [Google Scholar]
  15. Fontaine, G., Brassard, P., & Bergeron, P. 2001, PASP, 113, 409 [NASA ADS] [CrossRef] [Google Scholar]
  16. García-Berro, E., & Oswalt, T. D. 2016, New Astron. Rev., 72, 1 [NASA ADS] [CrossRef] [Google Scholar]
  17. García-Berro, E., Hernanz, M., Mochkovitch, R., & Isern, J. 1988, A&A, 193, 141 [Google Scholar]
  18. García-Berro, E., Althaus, L. G., Córsico, A. H., & Isern, J. 2008, ApJ, 677, 473 [NASA ADS] [CrossRef] [Google Scholar]
  19. García-Berro, E., Torres, S., Althaus, L. G., et al. 2010, Nature, 465, 194 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  20. Gentile Fusillo, N. P., Tremblay, P.-E., Gänsicke, B. T., et al. 2019, MNRAS, 482, 4570 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hansen, B. M. S., Kalirai, J. S., Anderson, J., et al. 2013, Nature, 500, 51 [NASA ADS] [CrossRef] [Google Scholar]
  22. Horowitz, C. J., Schneider, A. S., & Berry, D. K. 2010, Phys. Rev. Lett., 104, 231101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  23. Hubbard, W. B., & Lampe, M. 1969, ApJS, 18, 297 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hughto, J., Horowitz, C. J., Schneider, A. S., et al. 2012, Phys. Rev. E, 86, 066413 [CrossRef] [Google Scholar]
  25. Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 [NASA ADS] [CrossRef] [Google Scholar]
  26. Isern, J. 2019, ApJ, 878, L11 [NASA ADS] [CrossRef] [Google Scholar]
  27. Isern, J., Hernanz, M., Mochkovitch, R., & Garcia-Berro, E. 1991, A&A, 241, L29 [NASA ADS] [Google Scholar]
  28. Isern, J., Mochkovitch, R., García-Berro, E., & Hernanz, M. 1997, ApJ, 485, 308 [NASA ADS] [CrossRef] [Google Scholar]
  29. Isern, J., García-Berro, E., Hernanz, M., & Chabrier, G. 2000, ApJ, 528, 397 [NASA ADS] [CrossRef] [Google Scholar]
  30. Itoh, N., Mitake, S., Iyetomi, H., & Ichimaru, S. 1983, ApJ, 273, 774 [NASA ADS] [CrossRef] [Google Scholar]
  31. Kalirai, J. S. 2012, Nature, 486, 90 [NASA ADS] [Google Scholar]
  32. Kilic, M., Munn, J. A., Harris, H. C., et al. 2017, ApJ, 837, 162 [NASA ADS] [CrossRef] [Google Scholar]
  33. Kilic, M., Bergeron, P., Kosakowski, A., et al. 2020, ApJ, 898, 84 [CrossRef] [Google Scholar]
  34. Kofke, D. A. 1993a, Mol. Phys., 78, 1331 [CrossRef] [Google Scholar]
  35. Kofke, D. A. 1993b, J. Chem. Phys., 98, 4149 [CrossRef] [Google Scholar]
  36. Kowalski, P. M., & Saumon, D. 2006, ApJ, 651, L137 [Google Scholar]
  37. Medin, Z., & Cumming, A. 2010, Phys. Rev. E, 81, 036107 [NASA ADS] [CrossRef] [Google Scholar]
  38. Mitake, S., Ichimaru, S., & Itoh, N. 1984, ApJ, 277, 375 [NASA ADS] [CrossRef] [Google Scholar]
  39. Mochkovitch, R. 1983, A&A, 122, 212 [NASA ADS] [Google Scholar]
  40. Ogata, S., Iyetomi, H., Ichimaru, S., & van Horn, H. M. 1993, Phys. Rev. E, 48, 1344 [NASA ADS] [CrossRef] [Google Scholar]
  41. Oswalt, T. D., Smith, J. A., Wood, M. A., & Hintzen, P. 1996, Nature, 382, 692 [NASA ADS] [CrossRef] [Google Scholar]
  42. Potekhin, A. Y., & Chabrier, G. 2000, Phys. Rev. E, 62, 8554 [NASA ADS] [CrossRef] [Google Scholar]
  43. Salaris, M., Cassisi, S., Pietrinferni, A., Kowalski, P. M., & Isern, J. 2010, ApJ, 716, 1241 [NASA ADS] [CrossRef] [Google Scholar]
  44. Salaris, M., Althaus, L. G., & García-Berro, E. 2013, A&A, 555, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
  46. Segretain, L. 1996, A&A, 310, 485 [NASA ADS] [Google Scholar]
  47. Segretain, L., & Chabrier, G. 1993, A&A, 271, L13 [NASA ADS] [Google Scholar]
  48. Shaffer, N. R., & Starrett, C. E. 2020, Phys. Rev., E, 101 [Google Scholar]
  49. Staff, J. E., Menon, A., Herwig, F., et al. 2012, ApJ, 757, 76 [NASA ADS] [CrossRef] [Google Scholar]
  50. Tremblay, P. E., Bergeron, P., & Gianninas, A. 2011, ApJ, 730, 128 [NASA ADS] [CrossRef] [Google Scholar]
  51. Tremblay, P. E., Kalirai, J. S., Soderblom, D. R., Cignoni, M., & Cummings, J. 2014, ApJ, 791, 92 [NASA ADS] [CrossRef] [Google Scholar]
  52. Tremblay, P.-E., Fontaine, G., Fusillo, N. P. G., et al. 2019, Nature, 565, 202 [NASA ADS] [CrossRef] [Google Scholar]
  53. van Horn, H. M. 1968, ApJ, 151, 227 [NASA ADS] [CrossRef] [Google Scholar]
  54. Winget, D. E., Hansen, C. J., Liebert, J., et al. 1987, ApJ, 315, L77 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

C/O phase diagram. The liquidus (above which the plasma is entirely liquid) and the solidus (below which it is entirely solid) are both shown. The horizontal axis gives the number fraction of O, xO = NO/(NC + NO), and the vertical axis corresponds to the ratio of the temperature to the melting temperature of a pure carbon plasma. The results from this work (in red) are compared to those of Segretain & Chabrier (1993), Medin & Cumming (2010), and Horowitz et al. (2010). The region where an azeotrope is predicted is enlarged in the upper left corner.

In the text
thumbnail Fig. 2.

Lower panel: luminosity functions for hydrogen-atmosphere white dwarfs with masses between 0.9 and 1.1 M. Observational data from Tremblay et al. (2019) are shown in red and their population synthesis is represented by the dotted curve. Results from our own simulations are shown as dashed and solid lines (with and without an attempt to simulate the effect of 22Ne phase separation, respectively). The normalization is arbitrary, and the Gaia G magnitude limit from Gentile Fusillo et al. (2019) is used. Upper panel: evolution of the fraction of the core that is crystallized for the 1.1, 1.0, and 0.9 M sequences.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.