Abstract
White dwarfs represent the final state of evolution for most stars1,2,3. Certain classes of white dwarfs pulsate4,5, leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution6. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. Here we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments7,8, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure–density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars9.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
Source data are provided with this paper. Additional data are available upon request.
Code availability
Owing to its complexity, the data analysis algorithm that supports the findings of this study is available from D.C.S. upon reasonable request and will be outlined in more detail in a supporting publication.
References
Giammichele, N. et al. A large oxygen-dominated core from the seismic cartography of a pulsating white dwarf. Nature 554, 73–76 (2018).
Winget, D. E. et al. An independent method for determining the age of the Universe. Astrophys. J. 315, L77–L81 (1987).
Fontaine, G., Brassard, P. & Bergeron, P. The potential of white dwarf cosmochronology. Publ. Astron. Soc. Pacif. 113, 409–435 (2001).
Córsico, A. H., Althaus, L. G., Miller Bertolami, M. M. & Kepler, S. O. Pulsating white dwarfs: new insights. Astron. Astrophys. Rev. 27, 7–98 (2019).
Fontaine, G., Brassard, P. & Dufour, P. Might carbon-atmosphere white dwarfs harbour a new type of pulsating star? Astron. Astrophys. 483, L1–L4 (2008).
Winget, D. E. & Kepler, S. O. Pulsating white dwarf stars and precision asteroseismology. Annu. Rev. Astron. Astrophys. 46, 157–199 (2008).
Hurricane, O. A. et al. Fuel gain exceeding unity in an inertially confined fusion implosion. Nature 506, 343–348 (2014).
Gaffney, J. A. et al. Review of equation-of-state models for inertial confinement fusion materials. High Energ. Dens. Phys. 28, 7–24 (2018).
Fontaine, G. et al. The physics of pulsating white dwarf stars. ASP Conf. Series 462, 477–490 (2012).
Kraus, D. et al. X-ray scattering measurements on imploding CH spheres at the National Ignition Facility. Phys. Rev. E 94, 011202(R) (2016).
Zhang, S. et al. Path integral Monte Carlo simulations of dense carbon–hydrogen plasmas. J. Chem. Phys. 148, 102318 (2018).
Swift, D. C. et al. Sensitivity of stellar physics to the equation of state. Preprint at https://arxiv.org/abs/2006.16208 (2020).
Dufour, P., Liebert, J., Fontaine, G. & Behara, N. Rare white dwarf stars with carbon atmospheres. Nature 450, 522–524 (2007).
Dufour, P. et al. White dwarfs: something different. Astrophys. J. 683, 978–989 (2008).
Fontaine, G. & Van Horn, H. M. Convective white-dwarf envelope model grids for H-, He-, and C-rich compositions. Astrophys. J. Suppl. Ser. 31, 467–487 (1976).
Clark, D. S. et al. Plastic ablator ignition capsule design for the National Ignition Facility. Phys. Plasmas 17, 052703 (2010).
Hammel, B. A. et al. High-mode Rayleigh–Taylor growth in NIF ignition capsules. High Energ. Dens. Phys. 6, 171–178 (2010).
Vladimirov, A. S., Voloshin, N. P., Nogin, V. N., Petrovtsev, A. V. & Simonenko, V. A. Shock compressibility of aluminum at p ≥ 1 Gbar. JETP Lett. 39, 82–85 (1984).
Döppner, T. et al. Absolute equation-of-state measurement for polystyrene from 25–60 Mbar using a spherically converging shock wave. Phys. Rev. Lett.121, 025001 (2018).
Swift, D. C. et al. Gigabar Material Properties Experiments on NIF and Omega. Report No. LLNL-PROC-492261 (LLNL, 2011).
Swift, D. C. et al. Absolute Hugoniot measurements from a spherically convergent shock using X-ray radiography. Rev. Sci. Instrum. 89, 053505 (2018).
Kritcher, A. L. et al. Probing matter at Gbar pressures at the NIF. High Energ. Dens. Phys. 10, 27–34 (2014).
Smith, R. F. et al. Ramp compression of diamond to five terapascals. Nature 511, 330–333 (2014).
Smith, R. F. et al. Equation of state of iron under core conditions of large rocky exoplanets. Nat. Astron. 2, 452–458 (2018).
Swift, D. C. et al. Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker. Preprint at https://arxiv.org/abs/2006.13342 (2020).
Moses, E., Boyd, R., Remington, B., Keane, C. & Al-Ayat, R. The National Ignition Facility: ushering in a new age for high energy density science. Phys. Plasmas 16, 041006 (2009).
Bachmann, B. et al. Using penumbral imaging to measure micrometer size plasma hot spots in Gbar equation of state experiments on the National Ignition Facility. Rev. Sci. Instrum. 85, 11D614 (2014).
Scott, H. A. & Hansen, S. B. Advances in NLTE modeling for integrated simulations. High Energ. Dens. Phys. 6, 39–47 (2010).
Cauble, R. et al. Absolute equation-of-state data in the 10–40 Mbar (1–4 TPa) regime. Phys. Rev. Lett. 80, 1248–1251 (1998).
Ozaki, N. et al. Shock Hugoniot and temperature data for polystyrene obtained with quartz standard. Phys. Plasmas 16, 062702 (2009).
Barrios, M.A. et al. High-precision measurements of the equation of state of hydrocarbons at 1–10 Mbar using laser-driven shock waves. Phys. Plasmas 17, 056307 (2010).
Hu, S. X., Boehly, T. R. & Collins, L. A. Properties of warm dense polystyrene plasmas along the principal Hugoniot. Phys. Rev. E 89, 063104 (2014).
Sterne, P. A. et al. Equations of state for ablator materials in inertial confinement fusion simulations. J. Phys. Conf. Ser. 717, 012082 (2016).
Wilson, B., Sonnad, V., Sterne, P. & Isaacs, W. PURGATORIO—a new implementation of the INFERNO algorithm. J. Quant. Spectrosc. Radiat. Transf. 99, 658–679 (2006).
Benedict, L. X. et al. Multiphase equation of state for carbon addressing high pressures and temperatures. Phys. Rev. B 89, 224109 (2014).
Kerley, G. I. Sandia National Laboratories Report SAND2003, 3613 (2003).
Lyon, S. & Johnson, J. SESAME: The Los Alamos National Laboratory Equation of State Database. Technical Report No. LA-UR-92-3407 (Table 7592) (LANL, 1995).
Marinak, M. M. et al. Three-dimensional HYDRA simulations of National Ignition Facility targets. Phys. Plasmas 8, 2275–2280 (2001).
Berzak Hopkins, L. F. et al. First high-convergence cryogenic implosion in a near-vacuum hohlraum. Phys. Rev. Lett. 114, 175001 (2015).
Hohenberger, M. et al. Time-resolved measurements of the hot-electron population in ignition-scale experiments on the National Ignition Facility. Rev. Sci. Instrum. 85, 11D501 (2014).
Hinkel, D.E. et al. Development of improved radiation drive environment for high foot implosions at the National Ignition Facility. Phys. Rev. Lett. 117, 225002-1 (2016).
Dewald, E. L. et al. Dante soft X-ray power diagnostic for National Ignition Facility. Rev. Sci. Instrum. 75, 3759–3761 (2004).
Hamel, S. et al. Equation of state of CH1.36: first-principles molecular dynamics simulations and shock-and-release wave speed measurements. Phys. Rev. B 86, 094113 (2012).
Córsico, A. H., Romero, A. D., Althaus, L. G. & Garca-Berro, E. Hot C-rich white dwarfs: testing the DB–DQ transition through pulsations. Astron. Astrophys. 506, 835–843 (2009).
Montgomery, M. H. et al. SDSS J142625.71+575218.3: a prototype for a new class of variable white dwarf. Astrophys. J. 678, L51–L54 (2008).
Parsons, S. G. et al. A pulsating white dwarf in an eclipsing binary. Nat. Astron. 4, 690–696 (2020).
Acknowledgements
This work was performed under the auspices of the US Department of Energy under contract number DE-AC52-07NA27344 and contract number 89233218CNA000001. This work has been supported by Laboratory Directed Research and Development (LDRD) award number 13-ERD-073, and also by the University of California, Office of the President, Lab Fee Grant LFR-17-449059, by the Department of Energy, National Nuclear Security Administration, award DE-NA0003842, and by the Department of Energy, Office of Science, Office of Fusion Energy Sciences, awards DE-SC0018298, DE-SC0019269 and FWP 100182. We also acknowledge the NIF Discovery Science programme for providing access to the facility. This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favouring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes.
Author information
Authors and Affiliations
Contributions
A.L.K. wrote the manuscript and performed the design calculations. D.C.S. developed the analysis method and analysed the experimental data. B.B., B.M., M.J.M., M.E.M., J.N. and N.K. provided input and feedback to improve the data analysis method. T.D. fielded the experiments and follow-on supporting experiments together with B.B., D.K. and A.L. L.X.B., J.L.D., S.H., P.A.S., A.A.C. and H.D.W. provided theoretical calculations of equation of state models included in this paper and/or provided theoretical models for understanding the data. G.F. and D.S. performed calculations of white dwarf stars and worked closely with S.H.G. and A.L.K. to understand the impact of this work for white dwarfs. W.R.J., J.A.G. and J.N. provided calculations of theoretical opacities. W.S., F.E. and A.N. fabricated and characterized the targets for these experiments. D.C.S., S.H.G, R.W.F., P.N., B.A.R. and G.W.C. proposed the experiments and/or aided in obtaining experimental beamtime at NIF. All co-authors provided input on interpretation of the data and results and/or on their impact for white dwarf modelling.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature thanks Stephanie Hansen and Don Winget for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Schematic of the radiographic analysis method.
Normalized transmission intensity of an X-ray radiograph versus changes in compression and opacity. A given transmission corresponds to a range of possible compressions and opacities (black dashed curves). Density at the shock front is further constrained through knowledge of the initial material density ahead of the shock front and known mass contained within a region on the radiograph corresponding to a Ge marker layer in the sample. This further constrains the opacity and is included in the analysis as the shock traverses the sample. An example shock trajectory is denoted by a red dashed curve.
Extended Data Fig. 2 Comparison of two C9(H)10 shock Hugoniot measurements.
Measured pressure versus mass density (ρ) normalized to the initial density (ρ0) along the shock Hugoniot. Green, data from shot N130103-009-999 fielded at 24 K; red (from Fig. 3), shot N130701-002-999 fielded at ambient temperature. Also plotted are the theoretical Hugoniots for AA-DFT and AA-TFD at 24 K (ρ0 = 1.136 g cm−3), indicating that the initial density conditions are predicted to access approximately the same Hugoniot states at high pressure.
Extended Data Fig. 3 Sensitivity of the theoretical Hugoniot to fluorine.
Calculations show insensitivity of the theoretical Hugoniot (AA-DFT33,34,35) to fluorobenzene solvent (C6H5F) for concentrations up to 20%, corresponding to 1% atomic fraction of fluorine (green curve). Concentrations of 0.5% F (red curve) and 0% F (blue curve) are also shown, but not visibly distinguishable.
Extended Data Fig. 5 Extended C9(H)10 shock Hugoniot measurements.
Measured pressure versus mass density (ρ) normalized to the initial density (ρ0) along the shock Hugoniot from this work (red and purple curves and shaded region). The purple curve corresponds to the extended dataset that may be impacted by radiative shock-front preheat. Also plotted are previous experimental data and theoretical modelling of the Hugoniot (see Fig. 3).
Source data
Rights and permissions
About this article
Cite this article
Kritcher, A.L., Swift, D.C., Döppner, T. et al. A measurement of the equation of state of carbon envelopes of white dwarfs. Nature 584, 51–54 (2020). https://doi.org/10.1038/s41586-020-2535-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-020-2535-y
This article is cited by
-
Production of diamond using intense heavy ion beams at the FAIR facility and application to planetary physics
Scientific Reports (2023)
-
Observing the onset of pressure-driven K-shell delocalization
Nature (2023)
-
Probing atomic physics at ultrahigh pressure using laser-driven implosions
Nature Communications (2022)
-
Nonequilibrium band occupation and optical response of gold after ultrafast XUV excitation
Scientific Reports (2022)
-
Accurate temperature diagnostics for matter under extreme conditions
Nature Communications (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
