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.

  • Article
  • Published:

Shock cooling of a red-supergiant supernova at redshift 3 in lensed images

Abstract

The core-collapse supernova of a massive star rapidly brightens when a shock, produced following the collapse of its core, reaches the stellar surface. As the shock-heated star subsequently expands and cools, its early-time light curve should have a simple dependence on the size of the progenitor1 and therefore final evolutionary state. Measurements of the radius of the progenitor from early light curves exist for only a small sample of nearby supernovae2,3,4,5,6,7,8,9,10,11,12,13,14, and almost all lack constraining ultraviolet observations within a day of explosion. The several-day time delays and magnifying ability of galaxy-scale gravitational lenses, however, should provide a powerful tool for measuring the early light curves of distant supernovae, and thereby studying massive stellar populations at high redshift. Here we analyse individual rest-frame exposures in the ultraviolet to the optical taken with the Hubble Space Telescope, which simultaneously capture, in three separate gravitationally lensed images, the early phases of a supernova at redshift z ≈ 3 beginning within 5.8 ± 3.1 hours of explosion. The supernova, seen at a lookback time of approximately 11.5 billion years, is strongly lensed by an early-type galaxy in the Abell 370 cluster. We constrain the pre-explosion radius to be \(53{3}_{-119}^{+154}\) solar radii, consistent with a red supergiant. Highly confined and massive circumstellar material at the same radius can also reproduce the light curve, but because no similar low-redshift examples are known, this is unlikely.

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: A multiply imaged SN in the archival HST imaging.
Fig. 2: Posterior distributions of the progenitor radius R, envelope mass M and shock velocity v from the RSG model for RV = 3.1.
Fig. 3: Two examples of the reconstructed light curve of the SN with fixed time delays and magnifications.
Fig. 4: Volumetric CCSN rate (RCC) as a function of redshift.

Similar content being viewed by others

Data availability

The HST data used for this study can be retrieved from the NASA Mikulski Archive for Space Telescopes (http://archive.stsci.edu). The supernova is found in the HST imaging of the Abell 370 field acquired from programme GO-11591 (PI J.-P. Kneib). The LBT spectroscopy data are available from the LBT archive (http://archive.lbto.org). The Keck MOSFIRE data can be retrieved from the Keck Observatory Archive (https://www2.keck.hawaii.edu/koa/public/koa.php). The HFF data and models can be downloaded from https://archive.stsci.edu/prepds/frontier/lensmodels/#modelsandinput. Additional data including the MMT spectroscopic data, the HST coaddition and image differencing data, the GALFIT scripts and resulting models, the HST photometry data of the SN host galaxy, the SN light curve fitting script and resulting MCMC data, and our best-fit GLAFIC lens model are available from https://doi.org/10.5281/zenodo.6725770.

References

  1. Waxman, E. & Katz, B. in Handbook of Supernovae (eds Alsabti, A. & Murdin, P.) 967–1015 (Springer, 2017).

  2. Garnavich, P. M. et al. Shock breakout and early light curves of type II-P supernovae observed with Kepler. Astrophys. J. 820, 23 (2016).

    Article  ADS  Google Scholar 

  3. Arcavi, I. et al. SN 2011dh: discovery of a type IIb supernova from a compact progenitor in the nearby galaxy M51. Astrophys. J. Lett. 742, L18 (2011).

    Article  ADS  Google Scholar 

  4. Morales-Garoffolo, A. et al. SN 2011fu: a type IIb supernova with a luminous double-peaked light curve. Mon. Not. R. Astron. Soc 454, 95–114 (2015).

    Article  ADS  CAS  Google Scholar 

  5. Ben-Ami, S. et al. Discovery and early multi-wavelength measurements of the energetic type Ic supernova PTF12gzk: a massive-star explosion in a dwarf host galaxy. Astrophys. J. Lett. 760, L33 (2012).

    Article  ADS  Google Scholar 

  6. Valenti, S. et al. The first month of evolution of the slow-rising type IIP SN 2013ej in M74. Mon. Not. R. Astron. Soc. Lett. 438, L101–L105 (2014).

    Article  ADS  CAS  Google Scholar 

  7. Bersten, M. C. et al. A surge of light at the birth of a supernova. Nature 554, 497–499 (2018).

    Article  ADS  CAS  PubMed  Google Scholar 

  8. Tartaglia, L. et al. The progenitor and early evolution of the type IIb SN 2016gkg. Astrophys. J. Lett. 836, L12 (2017).

    Article  ADS  Google Scholar 

  9. Arcavi, I. et al. Constraints on the progenitor of SN 2016gkg from its shock-cooling light curve. Astrophys. J. Lett. 837, L2 (2017).

    Article  ADS  Google Scholar 

  10. Piro, A. L., Haynie, A. & Yao, Y. Shock cooling emission from extended material revisited. Astrophys. J. 909, 209 (2021).

    Article  ADS  CAS  Google Scholar 

  11. Szalai, T. et al. The type II-P supernova 2017eaw: from explosion to the nebular phase. Astrophys. J. 876, 19 (2019).

    Article  ADS  CAS  Google Scholar 

  12. Rui, L. et al. Probing the final-stage progenitor evolution for type IIP supernova 2017eaw in NGC 6946. Mon. Not. R. Astron. Soc 485, 1990–2000 (2019).

    Article  ADS  CAS  Google Scholar 

  13. Xiang, D. et al. Observations of SN 2017ein reveal shock breakout emission and a massive progenitor star for a type Ic supernova. Astrophys. J. 871, 176 (2019).

    Article  ADS  CAS  Google Scholar 

  14. Soumagnac, M. T. et al. SN 2018fif: the explosion of a large red supergiant discovered in its infancy by the Zwicky Transient Facility. Astrophys. J. Letters 902, 6 (2020).

    Article  CAS  Google Scholar 

  15. Lotz, J. M. et al. The Frontier Fields: survey design and initial results. Astrophys. J. 837, 97 (2017).

    Article  ADS  Google Scholar 

  16. Struble, M. F. & Rood, H. J. A compilation of redshifts and velocity dispersions for Abell clusters (Epoch 1991.2). Astrophys. J. Suppl. Ser. 77, 363–377 (1991).

    Article  Google Scholar 

  17. Benítez, N. Bayesian photometric redshift estimation. Astrophys. J. 536, 571–583 (2000).

    Article  Google Scholar 

  18. Brammer, G. B., van Dokkum, P. G. & Coppi, P. EAZY: a fast, public photometric redshift code. Astrophys. J. 686, 1503–1513 (2008).

    Article  Google Scholar 

  19. Kawamata, R. et al. Size–luminosity relations and UV luminosity functions at z = 6–9 simultaneously derived from the complete Hubble Frontier Fields data. Astrophys. J. 855, 4 (2018).

    Article  ADS  Google Scholar 

  20. Kawamata, R., Oguri, M., Ishigaki, M., Shimasaku, K. & Ouchi, M. Precise strong lensing mass modeling of four Hubble Frontier Field clusters and a sample of magnified high-redshift galaxies. Astrophys. J. 819, 114 (2016).

    Article  ADS  Google Scholar 

  21. Oguri, M. The mass distribution of SDSS J1004.4112 revisited. Publ Astron. Soc. Jpn 62, 1017–1024 (2010).

    Article  Google Scholar 

  22. Oguri, M. Fast calculation of gravitational lensing properties of elliptical Navarro–Frenk–White and Hernquist density profiles. Publ. Astron. Soc. Pacif. 133, 074504 (2021).

    Article  ADS  Google Scholar 

  23. Nakar, E. & Sari, R. Early supernovae light curves following the shock breakout. Astrophys. J. 725, 904–921 (2010).

    Article  Google Scholar 

  24. Rabinak, I. & Waxman, E. The early UV/optical emission from core-collapse supernovae. Astrophys. J. 728, 63 (2011); erratum 770, 81 (2013).

    Article  ADS  Google Scholar 

  25. Nakar, E. & Piro, A. L. Supernovae with two peaks in the optical light curve and the signature of progenitors with low-mass extended envelopes. Astrophys. J. 788, 193 (2014).

    Article  ADS  Google Scholar 

  26. Piro, A. L. Using double-peaked supernova light curves to study extended material. Astrophys. J. Lett. 808, L51 (2015).

    Article  ADS  Google Scholar 

  27. Sapir, N. & Waxman, E. UV/optical emission from the expanding envelopes of type II supernovae. Astrophys. J. 838, 130 (2017).

    Article  ADS  Google Scholar 

  28. Morozova, V., Piro, A. L. & Valenti, S. Unifying type II supernova light curves with dense circumstellar material. Astrophys. J. 838, 28 (2017).

    Article  ADS  Google Scholar 

  29. Morozova, V., Piro, A. L. & Valenti, S. Measuring the progenitor masses and dense circumstellar material of type II supernovae. Astrophys. J. 858, 15 (2018).

    Article  ADS  Google Scholar 

  30. Margalit, B. Analytic light curves of dense CSM shock breakout and cooling. Astrophys. J. 933, 238 (2022).

    Article  ADS  Google Scholar 

  31. Kass, R. E. & Raftery, A. E. Bayes factors. J. Am. Stat. Assoc. 90, 773–795 (1995).

    Article  MathSciNet  MATH  Google Scholar 

  32. Li, W. et al. Nearby supernova rates from the Lick Observatory Supernova Search – II. The observed luminosity functions and fractions of supernovae in a complete sample. Mon. Not. R. Astron. Soc 412, 1441–1472 (2011).

    Article  ADS  Google Scholar 

  33. Kelly, P. L. et al. Multiple images of a highly magnified supernova formed by an early-type cluster galaxy lens. Science 347, 1123–1126 (2015).

    Article  Google Scholar 

  34. Goobar, A. et al. iPTF16geu: a multiply imaged, gravitationally lensed type Ia supernova. Science 356, 291–295 (2017).

    Article  ADS  CAS  PubMed  Google Scholar 

  35. Rodney, S. A. et al. A gravitationally lensed supernova with an observable two-decade time delay. Nat. Astron. 5, 1118–1125 (2021).

  36. Oguri, M. Strong gravitational lensing of explosive transients. Rep. Prog. Phys. 82, 126901 (2019).

    Article  ADS  MathSciNet  CAS  PubMed  Google Scholar 

  37. Foxley-Marrable, M. et al. Observing the earliest moments of supernovae using strong gravitational lenses. Mon. Not. R. Astron. Soc 495, 4622–4637 (2020).

    Article  ADS  CAS  Google Scholar 

  38. Grogin, N. A. et al. CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. Astrophys. J. Suppl. Ser. 197, 35 (2011).

    Article  ADS  Google Scholar 

  39. Postman, M. et al. The cluster lensing and supernova survey with Hubble: an overview. Astrophys. J. Suppl. Ser. 199, 25 (2012).

    Article  ADS  Google Scholar 

  40. Strolger, L. G. et al. The rate of core collapse supernovae to redshift 2.5 from the CANDELS and CLASH supernova surveys. Astrophys. J. 813, 93 (2015).

    Article  ADS  Google Scholar 

  41. Schmidt, B. P. et al. The distance of five type II supernovae using the expanding photosphere method and the value of H0. Astrophys. J. 432, 42–48 (1994).

    Article  Google Scholar 

  42. Hamuy, M. & Pinto, P. A. Type II supernovae as standardized candles. Astrophys. J. 566, L63–L66 (2002).

    Article  Google Scholar 

  43. Kelly, P. L. & Kirshner, R. P. Core-collapse supernovae and host galaxy stellar populations. Astrophys. J. 759, 107 (2012).

    Article  ADS  Google Scholar 

  44. Drout, M. R. et al. The first systematic study of type Ibc supernova multi-band light curves. Astron. J. 741, 97 (2011); erratum 753, 180 (2012).

    Article  Google Scholar 

  45. Prentice, S. J. et al. The bolometric light curves and physical parameters of stripped-envelope supernovae. Mon. Not. R. Astron. Soc 458, 2973–3002 (2016).

    Article  ADS  CAS  Google Scholar 

  46. Madau, P. & Dickinson, M. Cosmic star-formation history. Annu. Rev. Astron. Astrophys. 52, 415–486 (2014).

    Article  ADS  Google Scholar 

  47. Salpeter, E. D. The luminosity function and stellar evolution. Astrophys. J. 121, 161–167 (1955).

    Article  Google Scholar 

  48. Fruchter, A. S., Hack, W., Dencheva, M., Droettboom, M., & Greenfield, P. BetaDrizzle: a redesign of the MultiDrizzle package. In 2010 Space Telescope Science Institute Calibration Workshop (eds Deustua, S. & Oliveira, C.) 382–387 (Space Telescope Science Institute, 2010).

  49. Jones, D. O., Scolnic, D. M., & Rodney, S. A. PythonPhot: simple DAOPHOT-type photometry in Python. Astrophysics Source Code Library https://www.ascl.net/1501.010 (2015).

  50. Stetson, P. B. DAOPHOT: a computer program for crowded-field stellar photometry. Publ. Astron. Soc. Pacif. 99, 191–222 (1987).

    Article  Google Scholar 

  51. Peng, C. Y., Ho, L. C., Impey, C. D. & Rix, H.-W. Detailed structural decomposition of galaxy images. Astron. J. 124, 266–293 (2002).

    Article  Google Scholar 

  52. Dong, S. et al. ASASSN-15lh: a highly super-luminous supernova. Science 351, 257–260 (2016).

    Article  ADS  CAS  PubMed  Google Scholar 

  53. Bradač, M. et al. Hubble Frontier Field photometric catalogues of Abell 370 and RXC J2248.7-4431: multiwavelength photometry, photometric redshifts, and stellar properties. Mon. Not. R. Astron. Soc. 489, 99–107 (2019).

    Google Scholar 

  54. Maraston, C. & Strömbäck, G. Stellar population models at high spectral resolution. Mon. Not. R. Astron. Soc 418, 2785–2811 (2011).

    Article  Google Scholar 

  55. Schulze, S. et al. The Palomar Transient Factory core-collapse supernova host-galaxy sample. I. Host-galaxy distribution functions and environment dependence of core-collapse supernovae. Astron. J. Suppl. Ser. 255, 29 (2021).

    Article  ADS  CAS  Google Scholar 

  56. Schmidt, G. D., Weymann, R. J. & Foltz, C. B. A moderate-resolution, high-throughput CCD channel for the MMT spectrograph. Publ. Astron. Soc. Pacif. 101, 713–724 (1989).

    Article  Google Scholar 

  57. Rothberg, B. et al. Current status of the facility instrumentation suite at the Large Binocular Telescope Observatory. In Proc. SPIE 10702: Ground-based and Airborne Instrumentation for Astronomy VII (eds Evans, C. J. et al.) 1070205 (SPIE, 2018).

  58. McLean, I. S. et al. MOSFIRE, the multi-object spectrometer for infrared exploration at the Keck Observatory. In Proc. SPIE 8446, Ground-based and Airborne Instrumentation for Astronomy IV (eds McLean, I. S. et al.) 84460J (SPIE, 2012).

  59. Prochaska, J. X. et al. PypeIt: the Python spectroscopic data reduction pipeline. Preprint at https://arxiv.org/abs/2005.06505 (2020).

  60. Prochaska, J. X. et al. pypeit/PypeIt: release 1.0.0. Zenodo https://doi.org/10.5281/zenodo.3743493 (2020).

  61. Konidaris, N. & Steidel, C. MOSFIRE DRP https://keck-datareductionpipelines.github.io/MosfireDRP/#mosfire-drp (2018).

  62. Wilkinson, D. M., Maraston, C., Goddard, D., Thomas, D. & Parikh, T. FIREFLY (Fitting IteRativEly For Likelihood analYsis): a full spectral fitting code. Mon. Not. R. Astron. Soc 472, 4297–4326 (2017).

    Article  Google Scholar 

  63. Reddy, N. A. et al. The MOSDEF Survey: significant evolution in the rest-frame optical emission line equivalent widths of star-forming galaxies at z = 1.4–3.8. Astrophys. J. 869, 92 (2018).

    Article  ADS  CAS  Google Scholar 

  64. Keeton, C. R. On modeling galaxy-scale strong lens systems. Gen. Rel. Grav. 42, 2151–2176 (2010).

    Article  MathSciNet  MATH  Google Scholar 

  65. McCully, C., Keeton, C. R., Wong, K. C. & Zabludoff, A. I. A new hybrid framework to efficiently model lines of sight to gravitational lenses. Mon. Not. R. Astron. Soc 443, 3631–3642 (2014).

    Article  Google Scholar 

  66. Ammons, S. M., Wong, K. C., Zabludoff, A. I. & Keeton, C. R. Mapping compound cosmic telescopes containing multiple projected cluster-scale halos. Astron. J. 781, 2 (2014).

    Article  Google Scholar 

  67. Johnson, T. L. et al. Lens models and magnification maps of the six Hubble Frontier Fields clusters. Astron. J. 797, 48 (2014).

    Article  Google Scholar 

  68. Jullo, E. et al. A Bayesian approach to strong lensing modelling of galaxy clusters. New J. Phys. 9, 447 (2007).

    Article  ADS  Google Scholar 

  69. Faber, S. M. & Jackson, R. E. Velocity dispersions and mass-to-light ratios for elliptical galaxies. Astrophys. J. 204, 668–683 (1976).

    Article  ADS  CAS  Google Scholar 

  70. Witt, H. J. & Mao, S. On the magnification relations in quadruple lenses: a moment approach. Mon. Not. R. Astron. Soc 311, 689–697 (2000).

    Article  ADS  Google Scholar 

  71. Drout, M. R. et al. Rapidly evolving and luminous transients from Pan-STARRS1. Astrophys. J. 794, 23 (2014).

    Article  ADS  Google Scholar 

  72. Ho, A. Y. Q. et al. AT2018cow: a luminous millimeter transient. Astrophys. J. 871, 73 (2019); erratum 935, 62 (2022).

    Article  ADS  CAS  Google Scholar 

  73. Margutti, R. et al. An embedded X-ray source shines through the aspherical AT2018cow: revealing the inner workings of the most luminous fast-evolving optical transients. Astrophys. J. 872, 18 (2019).

    Article  ADS  CAS  Google Scholar 

  74. Foreman-Mackey, D., Hogg, D. W., Lang, D. & Goodman, J. emcee: the MCMC hammer. Publ. Astron. Soc. Pacif. 125, 306–312 (2013).

    Article  Google Scholar 

  75. Chen, W. Additional data for ‘Shock cooling of a red-supergiant supernova at redshift 3 in lensed images’. Zenodo https://doi.org/10.5281/zenodo.6725770 (2022).

  76. Kriek, M. et al. An ultra-deep near-infrared spectrum of a compact quiescent galaxy at z = 2.2. Astrophys. J. 700, 221–231 (2009).

    Article  Google Scholar 

  77. Schreiber, C. & Dickinson, H. FAST++ v1.3.1. GitHub https://github.com/cschreib/fastpp (2021).

  78. Morishita, T. et al. Characterizing intracluster light in the Hubble Frontier Fields. Astrophys. J. 846, 139 (2017).

    Article  ADS  Google Scholar 

  79. Diego, J. M. et al. Dark matter under the microscope: constraining compact dark matter with caustic crossing events. Astrophys. J. 857, 25 (2018).

    Article  ADS  Google Scholar 

  80. Menzies, J. W. et al. Spectroscopic and photometric observations of SN 1987a: the first 50 days. Mon. Not. R. Astron. Soc 227, 39P–49P (1987).

    Article  Google Scholar 

  81. Hosseinzaden, G. et al. Weak mass loss from the red supergiant progenitor of the type II SN 2021yja. Preprint at https://arxiv.org/abs/2203.08155 (2022).

  82. Bullivant, C. et al. SN 2013fs and SN 2013fr: exploring the circumstellar-material diversity in type II supernovae. Mon. Not. R. Astron. Soc 476, 1497–1518 (2018).

    Article  ADS  CAS  Google Scholar 

  83. Valenti, S. et al. The diversity of type II supernova versus the similarity in their progenitors. Mon. Not. R. Astron. Soc 459, 3939–3962 (2016).

    Article  ADS  CAS  Google Scholar 

  84. Ho, A. Y. Q. et al. SN 2020bvc: a broad-line type Ic supernova with a double-peaked optical light curve and a luminous X-ray and radio counterpart. Astrophys. J. 902, 86 (2020).

    Article  ADS  CAS  Google Scholar 

  85. Yang, Y. et al. The young and nearby normal type Ia supernova 2018gv: UV–optical observations and the earliest spectropolarimetry. Astrophys. J. 902, 46 (2020).

    Article  ADS  CAS  Google Scholar 

  86. Ho, A. Y. Q. et al. The Koala: a fast blue optical transient with luminous radio emission from a starburst dwarf galaxy at z = 0.27. Astrophys. J. 895, 49 (2020).

    Article  ADS  CAS  Google Scholar 

  87. Richardson, D., Jenkins, R. L. III, Wright, J. & Maddox, L. Absolute-magnitude distributions of supernovae. Astrophys. J. 147, 118 (2014).

    Google Scholar 

  88. Taylor, M. et al. The core collapse supernova rate from the SDSS-II supernova survey. Astrophys. J. 792, 135 (2014).

    Article  ADS  Google Scholar 

  89. Hatano, K., Branch, D. & Deaton, J. Extinction and radial distribution of supernova properties in their parent galaxies. Astrophys. J. 502, 177–181 (1998).

    Article  Google Scholar 

  90. Calzetti, D. et al. The dust content and opacity of actively star-forming galaxies. Astrophys. J. 533, 682–695 (2000).

    Article  ADS  Google Scholar 

  91. Jeffreys, H. An invariant form for the prior probability in estimation problems. Proc. R. Soc. Lond. A 186, 453–461 (1946).

    Article  ADS  MathSciNet  CAS  MATH  Google Scholar 

  92. Kelly, P. L. & Refsdal, S. N. et al. Classification as a luminous and blue SN 1987A-like type II supernova. Astrophys. J. 831, 205 (2016).

    Article  ADS  Google Scholar 

  93. Pastorello, A. et al. SN 2009E: a faint clone of SN 1987A. Astron. Astrophys. 537, A141 (2012).

    Article  Google Scholar 

  94. Taddia, F. et al. Long-rising type II supernovae from Palomar Transient Factory and Caltech Core-Collapse Project. Astron. Astrophys. 588, A5 (2016).

    Article  Google Scholar 

  95. Coe, D., Benítez, N., Broadhurst, T. & Moustakas, L. A. A high-resolution mass map of galaxy cluster substructure: LensPerfect analysis of A1689. Astrophys. J. 723, 1678 (2010).

    Article  ADS  Google Scholar 

Download references

Acknowledgements

This work was supported by the HST Cycle 27 Archival Research programme (grant AR-15791), as well as by GO-15936 and GO-16278. We utilize gravitational lensing models produced by the GLAFIC group. The lens modelling was partially funded by the HST Frontier Fields programme conducted by STScI. The lens models were obtained from the Mikulski Archive for Space Telescopes. Some of the observations reported here were obtained at the MMT Observatory, a joint facility of the Smithsonian Institution and the University of Arizona. The data presented here were obtained in part at the LBT Observatory. The LBT is an international collaboration among institutions in the USA, Italy and Germany. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration (NASA). The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. We wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the Indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. P.L.K. is supported by US National Science Foundation (NSF) grant AST-1908823. J.M.D. acknowledges support from projects PGC2018-101814-B-100 and MDM-2017-0765. M.O. acknowledges support from World Premier International Research Center Initiative, MEXT, Japan, and JSPS KAKENHI grants JP20H00181, JP20H05856, JP22H01260 and JP18K03693. A.Z. acknowledges support by grant 2020750 from the USA–Israel Binational Science Foundation (BSF) and grant 2109066 from the US NSF, and by the Ministry of Science & Technology, Israel. A.V.F. is grateful for assistance from the Christopher R. Redlich Fund, the UC Berkeley Miller Institute for Basic Research in Science (where he was a Miller Senior Fellow), and many individual donors. We acknowledge the help of W. Zheng with the Keck MOSFIRE observations.

Author information

Authors and Affiliations

Authors

Contributions

W.C. analysed the HST, Keck, LBT and MMT data, wrote the manuscript and developed simulations. P.L.K. aided the interpretation of the events, planned the Keck, LBT and MMT observations, edited the manuscript and cross-checked the results. M.O. constructed the gravitational-lensing model for this work. T.J.B., J.M.D. and A.Z. helped with the HST proposal and reviewed the manuscript. P.L.K., N.E. and A.Z. acquired the Keck observation. P.L.K. and N.E. acquired the LBT observation. A.Z. constructed independent lens models for the cluster to verify the time-delay results. A.V.F. obtained the Keck time, helped acquire the Keck MOSFIRE observations, and edited the HST proposal and this manuscript. T.L.T. contributed to the manuscript and interpretation.

Corresponding author

Correspondence to Wenlei Chen.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

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 Photometry of the multiply imaged SN: image differencing and the blackbody fitting.

We used GALFIT51 to fit bright sources in the lensing system and subtracted the best-fit model from the field prior to measuring the flux. ac, The template image (a), the best-fit GALFIT model (b), and the residual through the HST ACS-WFC F814W filter (c). d, The WFC3-IR F160W difference image by subtracting the template from the event image directly. e, The difference image using our GALFIT-based method, where we subtracted the PSF-convolved best-fit GALFIT source models from the images, and then calculated the difference from GALFIT residuals between the coadded images and the event images. We can see that this method can reduce the significant residual from the mismatched PSFs around bright sources as shown in d. fh, Distributions of effective temperature and luminosity from the MCMC samples from fitting the blackbody emission into the photometry of each SN image.

Extended Data Fig. 2 Photometry and photometric redshift of the SN host galaxy.

ac, Broadband spectral energy distributions and best-fitting BPZ17 templates for G1, G2 and G3. The coloured circles designate the observed AB magnitudes with ±1σ uncertainties measured from ACS-WFC (green) and WFC3-IR (red) images, and arrows correspond to 95% upper limits. Dark-blue curves plot the best-fit spectral templates. Grey rectangles mark the magnitudes calculated from the best-fitting BPZ templates (with approximate uncertainties) for those filters95. df, Broadband spectral energy distributions and best-fitting EAZY17 templates for G1, G2 and G3. The data points are the observed flux with ±1σ uncertainties. Blue curves show the best-fit spectral template. The three panels on the right display the posterior probability distribution of the photometric redshift. g,h, Photometric redshift probability distributions for G1 (green dashed line), G2 (blue dashed line), G3 (red dashed line), and the joint analysis (solid line) derived from BPZ and EAZY fitting. i, The distribution of CCSN host-galaxy stellar mass (M) and star-formation rate (SFR)55, where the red star marks M and SFR of the host galaxy of the SN reported in this work.

Extended Data Fig. 3 Modelling the strongly lensed images.

a, Best-fit multiple-image geometry from the published Keeton V4 (left), GLAFIC V4 (middle) and Sharon V4 (right) models, by minimizing the r.m.s. angular separation of the predicted images from the real images in the image plane, for z = 3. b, Best-fit r.m.s. angular separation of the predicted images as a function of z from the published Keeton V4, GLAFIC V4 and Sharon V4 models. c, Best-fit multiple-image geometry from our revised GLAFIC model. d, Distribution of μS1, μS2 and μS3 from the MCMC samples given by our revised GLAFIC model, where μ1, μ2 and μ3 are predicted magnifications of images S1, S2 and S3, respectively. Cyan, yellow and magenta markers show the distribution projected on the μS1μS3, μS1μS2 and μS2μS3 planes, where red-dashed lines are the best-fit lines from linear regressions on the projected distribution. e, Correlations between Δt13 and Δt23 from the MCMC samples, where Δt13 and Δt23 are the time delay of S1 and S2 relative to S3, respectively.

Extended Data Fig. 4 Distribution of parameters from fitting the RSG model to the observations.

Here R is the extended envelope radius, M is the envelope mass, v is the shock velocity, t0 is the initial time, and Δt13 and μS1 are, respectively, the time delay and magnification of the S1 image. We set t = 0 at the observation time of the S3 image in the ACS-WFC F814W band. The time delay of each image is defined as the difference of its light-travel time from that of S3. The result is for RV = 3.1.

Extended Data Fig. 5 Distribution of parameters from fitting the CSM-homologous model to the observations.

The result is for a prior on the CSM mass of M < 1M. The parameters are the same as in Extended Data Fig. 5. Data of the MCMC samples for all models with all our choices of RV and priors are available from an online repository75.

Extended Data Fig. 6 Reconstructed light curves from the best-fit model parameters from all the light curve models for RV = 3.1.

af, Data shown are from the RSG model (a), the BSG model (b), the CSM-homologous model for a prior on the CSM mass of M < 0.1M (c), the CSM-homologous model for a prior on the CSM mass of M < 1M (d), the CSM-planar model for a prior on the CSM mass of M < 0.1M (e), and the CSM-planar model for a prior on the CSM mass of M < 1M (f). Data points and error bars correspond to the measured values and ±1σ uncertainties of the flux densities, respectively. Reconstructed light curves with all our choices of RV are available from an online repository75.

Extended Data Fig. 7 Comparisons of the reconstructed SN light curve with early UV light curves from other SNe.

a, The early-time evolution of the effective blackbody temperature of the SN reported here (green data points, where error bars are 68% confidence intervals), where the effective temperatures of the SN images are obtained by independently fitting the blackbody emission into the photometry of each SN image. The green solid line is from the best-fit RSG model. Red and blue dashed lines are two examples from the RSG model with progenitor radius of 200R and 1,000R, respectively, for the same envelope mass and shock velocity from our best-fit RSG model for the SN reported here. Red, yellow, blue and cyan data points show the early-time evolution of SN 2018fif14 (type IIP), SN 2013ej6 (type IIP/L), SN 2017eaw11,12 (type IIP), SN 1987A80 and SN 2016gkg7,8,9,10 (type IIb). bf, Stars and solid curves show the absolute magnitude of the SN reported here and reconstructed light curves from the best-fit RSG model, respectively, with arbitrary magnitude offsets for better visualization. For this SN, orange, teal and purple points and lines are for ACS-WFC F814W, WFC3-IR F110W and WFC3-IR F160W, respectively. In the rest frame of the SN, central wavelengths of these three filters fall into the B, UVW1, and UVW2 bands of Swift-UVOT. In b, we compare the SN’s early light curve to the Swift-UVOT observations of SN 2018fif14 and SN 2021yja81—two type-II SNe whose early light curves indicate only small amounts of CSM around their progenitors14,81. c,d, Comparison to the two type-II SNe, SN 2013fs82 and ASASSN-14gm83, which are believed to have dense CSM shells around their progenitors. In c,d, we plot the light curve from our best-fit CSM-homologous model for the SN of this work (dash-dotted lines) and light curves from the CSM-homologous model for the CSM radius and mass given by the analyses of CSM-rich CCSNe28,29 (dotted lines) for the two type-II SNe. e,f, Comparisons to the early UV light curves of SN 2016gkg8 (type IIb), SN 2020bvc84 (type Ib/c), SN 2018gv85 (type Ia), and ZTF18abvkwla86 (fast blue optical transient). ZTF-g band light curve is shown for ZTF18abvkwla at z = 0.27, corresponding to the rest-frame wavelength of 3,820 Å. All magnitudes are AB magnitudes. BB, blackbody.

Extended Data Fig. 8 Estimating the CCSN rate based on observations of the multiply imaged SNe.

a, Cumulative distribution of CCSN host galaxies as a function of the host-galaxy extinction AV43. The black solid line shows the cumulative distribution function of the exponential distribution P(AV) = λVexp(−λVAV) with best-fit λV = 0.98. Black, red and blue dashed lines show the cumulative distribution function of the exponential distribution with λV = 1, λV = 0.187 and λV = 2, respectively. b, Differential comoving volume (dV/dz) for the strongly lensed sources as a function of redshift z. Blue curve shows the lensed volume for all the sources in the six Hubble Frontier Fields within a 0.03° × 0.03° search window for each cluster field. Orange curve shows the effective lensing volume for sources that can be significantly (>5σ) detected by HST in the last decade. c, Differential number of detectable multiply imaged core-collapse supernovae (dN/dz) above the 5σ signal-to-noise-ratio level by HST in the last decade as a function of z. d, Volumetric CCSN rate (RCC) as a function of redshift. Contents are the same as in Fig. 4, but for λV = 0.187. e. Volumetric CCSN rate (RCC) as a function of redshift for different choices of λV and RV. Red, green and blue shaded regions are for λV = 0.187, λV = 1 and λV = 2, respectively, for RV = 4.05. Red, green and blue lines are from the MLE kCC for the three λV options. The green dotted line is the same as the green solid line, but for RV = 3.1.

Extended Data Table 1 Results of photometry, gravitational-lensing modelling and the Bayesian test of light curve models
Extended Data Table 2 Best-fit values of free model parameters from light curve models

Supplementary information

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, W., Kelly, P.L., Oguri, M. et al. Shock cooling of a red-supergiant supernova at redshift 3 in lensed images. Nature 611, 256–259 (2022). https://doi.org/10.1038/s41586-022-05252-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-022-05252-5

This article is cited by

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.

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