Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

A spectroscopic quadruple as a possible progenitor of sub-Chandrasekhar type Ia supernovae

Abstract

Binaries have received much attention as possible progenitors of type Ia supernova explosions1, but long-term gravitational effects2,3 in tight triple or quadruple systems could also play a key role in producing type Ia supernovae. Here we report on the properties of a spectroscopic quadruple found within a star cluster: the 2 + 2 hierarchical system HD 744384. Its membership in the open cluster IC 2391 makes it the youngest (43 Myr) spectroscopic quadruple known to us, and among the quadruple systems with the shortest outer orbital period. The eccentricity of the 6 yr outer period is 0.46, and the two inner orbits, with periods of 20.6 d and 4.4 d and eccentricities of 0.36 and 0.15, are not coplanar. Using an innovative combination of ground-based high-resolution spectroscopy5,6,7 and Gaia/Hipparcos astrometry8,9,10,11, we show that this system is undergoing secular interaction, which probably pumped the eccentricity of one of the inner orbits higher than expected for the spectral types of its components. We compute the future evolution of HD 74438 and show that this system is an excellent candidate progenitor of sub-Chandrasekhar type Ia supernovae through white dwarf mergers. Taking into account the contribution of this specific type of type Ia supernova accounts for the chemical evolution of iron-peak elements in the Galaxy better than considering only near Chandrasekhar-mass type Ia supernovae12.

This is a preview of subscription content, access via your institution

Access options

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Spectral fitting of the HRS/SALT spectra (taken on 2018 December 31).
Fig. 2: Location of HD 74438 in the HR diagram.
Fig. 3: Comparison of the orbital parameters of HD 74448 with those of other quadruple and binary systems.
Fig. 4: MSE simulations of the HD 74438 future.
Fig. 5: Distribution of merging WDs in our simulations.

Similar content being viewed by others

Data availability

Supplementary data are provided with this paper: the measured RVs from which all the orbital solutions are computed are available in a separate csv file. In addition, raw spectroscopic data are available in the associated observatory archives:

http://archive.eso.org/scienceportal for ESO/GES and GIRAFFE spectra.

https://ssda.saao.ac.za/ for SALT/HRS spectra. The spectra from the monitoring programme 2018-1-MLT-009 are publicly available except for the 2020-2-MLT-003 programme (three spectra, available from mid-2024 on).

• The HERCULES/UCMJO raw spectra will be made publicly available on the VizieR/CDS repository.

Code availability

The MSE code used to perform the simulation of the HD 74438 system is available upon request to A.S.H.

References

  1. Maoz, D., Manucci, P. & Nelemans, G. Observational clues to the progenitors of type Ia supernovae. Annu. Rev. Astron. Astrophys. 52, 107–170 (2014).

    Article  ADS  Google Scholar 

  2. Lidov, M. L. The evolution of orbits of artificial satellites of planets under the action of gravitational perturbations of external bodies. Planet. Space Sci. 9, 719–759 (1962).

    Article  ADS  Google Scholar 

  3. Kozai, Y. Secular perturbations of asteroids with high inclination and eccentricity. Astron. J. 67, 591–598 (1962).

    Article  ADS  MathSciNet  Google Scholar 

  4. Merle, T. et al. The Gaia-ESO Survey: double-, triple-, and quadruple-line spectroscopic binary candidates. Astron. Astrophys. 608, A95 (2017).

    Article  Google Scholar 

  5. Gilmore, G. et al. The Gaia-ESO Public Spectroscopic Survey. Messenger 147, 25–31 (2012).

    ADS  Google Scholar 

  6. Crause, L. A. et al. Performance of the Southern African Large Telescope (SALT) High Resolution Spectrograph (HRS). Proc. SPIE 9147, 6–20 (2014).

    Google Scholar 

  7. Hearnshaw, J. B. et al. The HERCULES Echelle Spectrograph at Mt John. Exp. Astron. 13, 59–76 (2002).

    Article  ADS  Google Scholar 

  8. Gaia Collaboration et al. Gaia Data Release 2. Summary of the contents and survey properties. Astron. Astrophys. 616, A1 (2018).

  9. Gaia Collaboration. Gaia Early Data Release 3. Summary of the contents and survey properties. Astron. Astrophys. 649, A1 (2021).

    Article  Google Scholar 

  10. ESA. The HIPPARCOS and TYCHO Catalogues. Astrometric and Photometric Star Catalogues Derived from the ESA HIPPARCOS Space Astrometry Mission. SP-1200 (ESA Publications Division, 1997).

  11. Høg, E. et al. The Tycho-2 catalogue of the 2.5 million brightest stars. Astron. Astrophys. 355, L27–L30 (2000).

    ADS  Google Scholar 

  12. Eitner, P. et al. Observational constraints on the origin of the elements III. Evidence for the dominant role of sub-Chandrasekhar SN Ia in the chemical evolution of Mn and Fe in the Galaxy. Astron. Astrophys. 635, A38 (2020).

    Article  Google Scholar 

  13. Gaia Collaboration et al. Gaia Data Release 2. Observational Hertzsprung–Russell diagrams. Astron. Astrophys. 616, A10 (2018).

    Article  ADS  Google Scholar 

  14. Randich, S. et al. The Gaia-ESO Survey: open clusters in Gaia-DR1. A way forward to stellar age calibration. Astron. Astrophys. 612, A99 (2018).

    Article  Google Scholar 

  15. Spina, L. et al. The Gaia-ESO Survey: the present-day radial metallicity distribution of the Galactic disc probed by pre-main-sequence clusters. Astron. Astrophys. 601, A70 (2017).

    Article  Google Scholar 

  16. Platais, I. et al. WIYN open cluster study. XXVI. Improved kinematic membership and spectroscopy of IC 2391. Astron. Astrophys. 461, 509–522 (2007).

    Article  ADS  Google Scholar 

  17. Bressan, A. et al. PARSEC: stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code. Mon. Not. R. Astron. Soc. 427, 127–145 (2012).

    Article  ADS  Google Scholar 

  18. Andrae, R. et al. Gaia Data Release 2. First stellar parameters from Apsis. Astron. Astrophys. 616, A8 (2018).

    Article  Google Scholar 

  19. Bravi, L. et al. The Gaia-ESO Survey: a kinematical and dynamical study of four young open clusters. Astron. Astrophys. 615, A37 (2018).

    Article  Google Scholar 

  20. Pineda, J. E. et al. The formation of a quadruple star system with wide separation. Nature 518, 213–215 (2015).

    Article  ADS  Google Scholar 

  21. Lee, A. T., Offner, S. S. R., Kratter, K. M., Smullen, R. A. & Li, P. S. The formation and evolution of wide-orbit stellar multiples in magnetized clouds. Astrophys. J. 887, 232 (2019).

    Article  ADS  Google Scholar 

  22. Tokovinin, A. The updated Multiple Star Catalog. Astrophys. J. Suppl. Ser. 235, 6 (2018).

    Article  ADS  Google Scholar 

  23. Muterspaugh, M. W. et al. Masses, luminosities, and orbital coplanarities of the μ Orionis quadruple-star system from PHASES differential astrometry. Astron. J. 135, 766–776 (2008).

    Article  ADS  Google Scholar 

  24. Tokovinin, A. From binaries to multiples. II. Hierarchical multiplicity of F and G dwarfs. Astron. J. 147, 87 (2014).

    Article  ADS  Google Scholar 

  25. Naoz, S. The eccentric Kozai–Lidov effect and its applications. Annu. Rev. Astron. Astrophys. 54, 441–489 (2016).

    Article  ADS  Google Scholar 

  26. Cochran, W. D., Hatzes, A. P., Butler, R. P. & Marcy, G. W. Detection of a planetary companion to 16 Cygni B. Bull. Am. Astron. Soc. 28, 1111 (1996).

    ADS  Google Scholar 

  27. Baron, F. et al. Imaging the Algol triple system in the H band with the CHARA interferometer. Astrophys. J. 752, 20 (2012).

    Article  ADS  Google Scholar 

  28. Pejcha, O., Antognini, J. M., Shappee, B. J. & Thompson, T. A. Greatly enhanced eccentricity oscillations in quadruple systems composed of two binaries: implications for stars, planets and transients. Mon. Not. R. Astron. Soc. 431, 943–951 (2013).

    Article  ADS  Google Scholar 

  29. Hamers, A. S. & Lai, D. Secular chaotic dynamics in hierarchical quadruple systems, with applications to hot Jupiters in stellar binaries and triples. Mon. Not. R. Astron. Soc. 470, 1657–1672 (2017).

    Article  ADS  Google Scholar 

  30. Hamers, A. S. On the rates of Type Ia supernovae originating from white dwarf collisions in quadruple star systems. Mon. Not. R. Astron. Soc. 478, 620–637 (2018).

    Article  ADS  Google Scholar 

  31. Fang, X., Thompson, T. A. & Hirata, C. M. Dynamics of quadruple systems composed of two binaries: stars, white dwarfs, and implications for Ia supernovae. Mon. Not. R. Astron. Soc. 476, 4234–4262 (2018).

    Article  ADS  Google Scholar 

  32. Pourbaix, D. et al. SB9: the Ninth Catalogue of Spectroscopic Binary Orbits. Astron. Astrophys. 424, 727–732 (2004).

    Article  ADS  Google Scholar 

  33. Zasche, P. et al. Doubly eclipsing systems. Astron. Astrophys. 630, A128 (2019).

    Article  Google Scholar 

  34. Geller, A. M., Hurley, J. R. & Mathieu, R. D. Direct N-body modeling of the old open cluster NGC 188: a detailed comparison of theoretical and observed binary star and blue straggler populations. Astron. J. 145, 8 (2013).

    Article  ADS  Google Scholar 

  35. Hamers, A. S., Rantala, A., Neunteufel, P., Preece, H. & Vynatheya, P. Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems. Mon. Not. R. Astron. Soc. 502, 4479–4512 (2021).

    Article  ADS  Google Scholar 

  36. Wang, B. & Han, Z. Progenitors of type Ia supernovae. N. Astron. Rev. 56, 122–141 (2012).

    Article  ADS  Google Scholar 

  37. Hamers, A. S., Pols, O. R., Claeys, J. S. W. & Nelemans, G. Population synthesis of triple systems in the context of mergers of carbon–oxygen white dwarfs. Mon. Not. R. Astron. Soc. 430, 2262–2280 (2013).

    Article  ADS  Google Scholar 

  38. Toonen, S., Perets, H. B. & Hamers, A. S. Rate of WD-WD head-on collisions in isolated triples is too low to explain standard type Ia supernovae. Astron. Astrophys. 610, A22 (2018).

    Article  ADS  Google Scholar 

  39. Moe, M. & Di Stefano, R. Mind your Ps and Qs: the interrelation between period (P) and mass-ratio (q) distributions of binary stars. Astrophys. J. Suppl. Ser. 230, 15 (2017).

    Article  ADS  Google Scholar 

  40. Flörs, A. et al. Sub-Chandrasekhar progenitors favoured for Type Ia supernovae: evidence from late-time spectroscopy. Mon. Not. R. Astron. Soc. 491, 2902–2918 (2020).

    ADS  Google Scholar 

  41. Ochsenbein, F., Bauer, P. & Marcout, J. The VizieR database of astronomical catalogues. Astron. Astrophys. Suppl. Ser. 143, 23–32 (2000).

    Article  ADS  Google Scholar 

  42. Sacco, G. G. et al. The Gaia–ESO Survey: processing FLAMES–UVES spectra. Astron. Astrophys. 565, A113 (2014).

    Article  Google Scholar 

  43. Crawford, S. M. et al. PySALT: the SALT science pipeline. Proc. SPIE 7737, 773725 (2010).

    Article  Google Scholar 

  44. Kurucz, R. L. in Stellar Atmospheres: Beyond Classical Models (eds L. Crivellari et al.) 441 (NATO ASI Series C Vol. 341, Springer, 1991)..

  45. Plez, B. Turbospectrum: code for spectral synthesis. Astrophysics Source Code Library ascl:1205.004 (2012).

  46. Ryabchikova, T. et al. A major upgrade of the VALD database. Phys. Scr. 90, 054005 (2015).

    Article  ADS  Google Scholar 

  47. Heiter, U. et al. Atomic and molecular data for optical stellar spectroscopy. Phys. Scr. 90, 054010 (2015).

    Article  ADS  Google Scholar 

  48. Moya, A., Zuccarino, F., Chaplin, W. J. & Davies, G. R. Empirical relations for the accurate estimation of stellar masses and radii. Astrophys. J. 237, 21 (2018).

    Article  Google Scholar 

  49. Jancart, S., Jorissen, A., Babusiaux, C. & Pourbaix, D. Astrometric orbits of SB9 stars. Astron. Astrophys. 442, 365–380 (2005).

    Article  ADS  Google Scholar 

  50. Tokovinin, A. & Moe, M. Formation of close binaries by disc fragmentation and migration, and its statistical modelling. Mon. Not. R. Astron. Soc. 491, 5158–5171 (2020).

    Article  ADS  Google Scholar 

  51. van den Berk, J., Portegies Zwart, S. F. & McMillan, S. L. W. The formation of higher order hierarchical systems in star clusters. Mon. Not. R. Astron. Soc. 379, 111–122 (2007).

    Article  ADS  Google Scholar 

  52. Wall, J. E., McMillan, S. L. W., Mac Low, M.-M., Klessen, R. S. & Portegies Zwart, S. Collisional N-body dynamics coupled to self-gravitating magnetohydrodynamics reveals dynamical binary formation. Astrophys. J. 887, 62 (2019).

    Article  ADS  Google Scholar 

  53. Hamers, A. S., Glanz, H. & Neunteufel, P. A statistical view on stable and unstable Roche lobe overflow of a tertiary star onto the inner binary in triple systems. Astrophys. J. Suppl. Ser. 259, 25 (2022).

    Article  ADS  Google Scholar 

Download references

Acknowledgements

T.M. and S.V.E. are supported by a grant from the Fondation ULB. A.S.H. thanks the Max Planck Society for support through a Max Planck Research Group. T.M. acknowledges M.G. Perderson and T. Van Reeth for investigating the TESS data. T.M. acknowledges A. Tokovinin for explanations about the Multiple Star Catalog. T.Z. and G.T. acknowledge financial support from the Slovenian Research Agency (research core funding P1-0188) and from the European Space Agency (PRODEX Experiment Arrangement C4000127986). G.T. acknowledges support by the Swedish strategic research programme eSSENCE, project grant ‘The New Milky Way’ from the Knut and Alice Wallenberg foundation and grant 2016-03412 from the Swedish Research Council. Based on data products from observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF/Osservatorio Astrofisico di Arcetri. These data have been obtained from the GES Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. This work was partly supported by the European Union FP7 programme through ERC grant 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge support from INAF and Ministero dell’Istruzione, dell’Università e della Ricerca (MIUR) in the form of the grant ‘Premiale VLT 2012’. The results presented here have benefited from discussions held during the Gaia–ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. Some of the observations reported in this paper were obtained with SALT under programmes 2018-1-MLT-009 (principal investigator R.S.) and 2020-2-MLT-03 (principal investigator R. Manick). Polish participation in SALT is funded by grant MNiSW DIR/WK/2016/07. This work has made use of data from the European Space Agency mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This work has made use of the VALD database, operated at Uppsala University, the Institute of Astronomy RAS in Moscow and the University of Vienna, the SIMBAD database, operated at CDS, Strasbourg, France, and the VizieR catalogue access tool, CDS, Strasbourg, France (https://doi.org/10.26093/cds/vizier). The original description of the VizieR service was published in ref. 41. This research has made use of Python3 and IPython, and modules NumPy, SciPy and Pandas. All the graphics were generated with Matplotlib except Extended Data Fig. 5.

Author information

Authors and Affiliations

Authors

Contributions

T.M. initiated the project, performed the follow-up with HRS/SALT and the analysis and interpretation of the data and designed most figures. A.S.H. performed the simulations of the evolution of the quadruple and designed some of the figures. S.V.E. and A.J. contributed to the analysis and interpretation of the results. T.M., S.V.E., A.J. and A.S.H. wrote the manuscript, with input from all authors. M.V.d.S. and D.P. contributed to the analysis of the data. K.P. performed the acquisition and reduction of HERCULES/UCMJO spectra. R.S. contributed to the follow-up with HRS/SALT and to the interpretation of the data. T.Z. and G.T. contributed to the analysis and interpretation of the data as well as to the writing. G.G. and S.R. are the principal investigators of the GES in which the quadruple was discovered. A.G., A.H., G.S. and C.C.W. contributed to the acquisition and reduction of the GES data. All authors provided critical feedback.

Corresponding author

Correspondence to Thibault Merle.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Astronomy thanks Andrei Tokovinin and Tamas Borkovits for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Cross-correlation functions computed from high-resolution spectra of HD 74438 at different epochs as labeled (black lines) and multi-Gaussian fits (red lines).

The vertical dotted line is the mean velocity of IC 2391. Left: HERCULES/UCMJO CCFs (with multiple spectra observed in some of the nights). The time series covers one week and the vertical position of the spectra scales with time. Right: same as left for HRS/SALT CCFs. The time series covers 7 months. The vertical shift of the spectra is arbitrary. Spectral fitting was performed using the two spectra marked with an asterisk and plotted in bold.

Extended Data Fig. 2 RV solutions using HERCULES/UCMJO and HRS/SALT data points for the brightest AB pair (left) and for the faintest CD pair (right).

The GES data points are not included in the computation of the orbital solutions because taken 5 y before, showing the gravitational effect between the two inner pairs. The uncertainties on the Gaussian mixture model fit to the CCFs provide the error bars on the RVs.

Extended Data Fig. 3 Orbital parameters of the two inner orbits (A-B and C-D) and the outer one (AB–CD) of the quadruple 2 + 2 system HD 74438.

P is the orbital period, e the eccentricity, ω the periastron argument, T0 the julian date at periastron, v0 the center-of-mass velocity, K1 and K2 the radial velocity amplitudes of the primary and secondary in each orbit. qspec and qdyn are the spectroscopic and dynamical mass ratios. i is the inclination of the orbit on the sky, a the semi-major axis of the orbit around the center of mass, \(\mu _{{{{\mathrm{phot}}}}}^{\prime\prime}\) the proper motion of the photocentre and Ω the argument of the ascending node.

Extended Data Fig. 4 Astrophysical parameters of the quadruple system HD 74438 and its components.

Teff is the effective temperature, M, L and R are the mass, the luminosity and the radius in solar units.

Extended Data Fig. 5 RV solutions of the wide pair AB–CD using center of mass RVs of the AB and CD pairs.

The error bars on the center-of-mass velocities of each pair come from the error propagation of the individual component RVs.

Extended Data Fig. 6 Orbital solutions from Extended Data Fig. 3 and measured RVs for the four components of the SB4.

Top: GIRAFFE data, middle top: UVES/VLT from GES data, middle bottom: HERCULES/UCMJO data, bottom: HRS/SALT data. Not all the RVs are presented here.

Extended Data Fig. 7 The two possible orbits of the photocentre of the AB pair around the centre of mass of the AB - CD system (in units of the semi-major axis a″AB-CBphot).

The black solid line corresponds to the orbit of inclination 73.2°, whereas the red solid line corresponds to the orbit with inclination 106.8°. On these orbits, following the direction of orbital motion, the solid circle marks the periastron, the open circle marks the epoch 2015.3 (April 24) when the predicted orbital velocity projected on the plane of the sky (solid-line arrow) matches the Gaia DR2 proper motion (differential with respect to the cluster, and expressed in units of the proper-motion modulus μAB-CDphot) averaged over the Gaia DR2 time span (represented by the thick part of the orbit), the triangle marks the average Gaia DR2 epoch (2015.5), the open square marks the Gaia eDR3 epoch (2016.0). For this epoch, the solid-line arrow marks the predicted orbital motion projected on the plane of the sky, whereas the dashed arrow corresponds to the Gaia eDR3 proper motion (differential with respect to the cluster). The counter-clockwise evolution of the observed proper motion (Gaia DR2 and Gaia eDR3) favours the orbit with 73.2° inclination (black line).

Extended Data Fig. 8 Probability distributions of merger times in the simulations.

The distribution of times of all merger events is displayed with the blue dotted lines, and the times of the first merger in the system (if applicable) with the red solid line. The distributions of the maximum age reached by the system for the systems in which the maximum wall time was exceeded is shown with the green dashed line.

Extended Data Fig. 9 Probability distributions of the current mutual inclinations (left panel: Φ AB / AB–CD; right panel: Φ CD / AB–CD) leading to different events in our Monte Carlo simulations (as described in the legends).

The curve labelled ICs corresponds to the distribution of the initial mutual inclinations (encompassing all possible outcomes). This distribution is obtained from the observed values of the individual orbital inclinations on the sky, complemented by flat distributions for the longitudes of the unknown ascending nodes (equation 12).

Extended Data Fig. 10 One possible future evolution of HD 74438.

Each panel shows an event of interest as labelled at the top, with the time indicated in each subpanel and the system represented schematically in a so-called mobile diagram35, with orbital parameters (semimajor axes a and eccentricities e) and masses (in units of M) indicated. The meaning of the dot colours (representing the stars) is indicated in the legend at the top. The quadruple experiences unstable RLOF, triple CE and merger events, leaving a WD remnant with a sub-Chandrasekhar mass of 1.3 M.

Supplementary information

Supplementary Information

Supplementary Notes and Tables 1–4.

Supplementary Data 1

RV measurements in kilometres per second of (i) each component (A, B, C and D) and (ii) centre of mass of pairs AB and CD, together with their error bars. The associated epochs are given in Gregorian years.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Merle, T., Hamers, A.S., Van Eck, S. et al. A spectroscopic quadruple as a possible progenitor of sub-Chandrasekhar type Ia supernovae. Nat Astron 6, 681–688 (2022). https://doi.org/10.1038/s41550-022-01664-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41550-022-01664-5

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing