Signatures of Exploding Supermassive PopIII Stars at High Redshift in JWST, EUCLID and Roman Space Telescope

Cédric Jockel, Kyohei Kawaguchi, Sho Fujibayashi, Masaru Shibata

First listed 2025-07-21 | Last updated 2025-07-21

Abstract

Recently discovered supermassive black holes with masses of $\sim10^8\,M_\odot$ at redshifts $z\sim9$-$11$ in active galactic nuclei (AGN) pose severe challenges to our understanding of supermassive black hole formation. One proposed channel are rapidly accreting supermassive PopIII stars (SMSs) that form in large primordial gas halos and grow up to $<10^6\,M_\odot$. They eventually collapse due to the general relativistic instability and could lead to supernova-like explosions. This releases massive and energetic ejecta that then interact with the halo medium via an optically thick shock. We develop a semi-analytic model to compute the shock properties, bolometric luminosity, emission spectrum and photometry over time. The initial data is informed by stellar evolution and general relativistic SMS collapse simulations. We find that SMS explosion light curves reach a brightness $\sim10^{45\mathrm{-}47}\,\mathrm{erg/s}$ and last $10$-$200$ years in the source frame - up to $250$-$3000$ years with cosmic time dilation. This makes them quasi-persistent sources which vary indistinguishably to little red dots and AGN within $0.5$-$9\,(1+z)$ yrs. Bright SMS explosions are observable in long-wavelength JWST filters up to $z\leq20$ ($24$-$26$ mag) and pulsating SMSs up to $z\leq15$. EUCLID and the Roman space telescope (RST) can detect SMS explosions at $z<11$-$12$. Their deep fields could constrain the SMS rate down to $10^{-11}$Mpc$^{-3}$yr$^{-1}$, which is much deeper than JWST bounds. Based on cosmological simulations and observed star formation rates, we expect to image up to several hundred SMS explosions with EUCLID and dozens with RST deep fields.

Short digest

Semi-analytic models, anchored to stellar-evolution and GR collapse simulations, follow ejecta–CSM shocks from exploding supermassive Pop III stars to predict light curves, spectra, and photometry. The events peak at ~10^45–10^47 erg/s and last 10–200 yr in the source frame (observed 250–3000 yr), yielding quasi-persistent variability indistinguishable from little red dots/AGN on 0.5–9(1+z) yr baselines. Bright explosions should be detectable in long-wavelength JWST filters to z≤20 at 24–26 mag, with pulsating SMSs to z≤15; Euclid and Roman reach z<11–12. Deep fields could constrain the event rate to ~10^-11 Mpc^-3 yr^-1 and yield several hundred (Euclid) to dozens (Roman) detections, directly testing heavy-seed SMBH pathways.

Key figures to inspect

• Figure 1: Use the phase cartoon to track how a GR-unstable SMS forms a BH, ejects a core-driven outflow, sweeps up the bloated atmosphere, and transitions to an ejecta–CSM shock; note how core/atmosphere fractions set the eventual ejecta mass budget.
• Figure 2: Inspect the ejecta, shocked shell, and CSM density profiles with forward/reverse shock jump conditions to read off diffusion timescales and where radiation should remain thermalized versus become non-thermal.
• Figure 3: Bolometric light curves (thermal vs non-thermal phases) show quasi-persistence and compare directly to two high‑z AGN luminosities; use these panels to infer detectability windows and observer-frame durations from the optically thick/thin transition.
• Figure 4: Parameter-space map of peak optically thick-phase luminosity versus ejecta mass and kinetic energy; identify which model combinations exceed survey limits and note the gray region where thin-shock luminosity dominates (model validity boundary).