skip to main content

Title: Shock Breakout in Three-dimensional Red Supergiant Envelopes

UsingAthena++, we perform 3D radiation-hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) that has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density “halo” of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to ≈1–2 hr. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to ≈3–6 hr and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminosities more » by a factor of 3–10 (Lbol∼ 1044erg s−1), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects.

« less
; ;
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 164
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this

    We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure.more »Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

    « less
  2. Abstract

    A growing number of core-collapse supernovae (SNe) that show evidence for interaction with dense circumstellar medium (CSM) are accompanied by “precursor” optical emission rising weeks to months prior to the explosion. The precursor luminosities greatly exceed the Eddington limit of the progenitor star, implying that they are accompanied by substantial mass loss. Here, we present a semi-analytic model for SN precursor light curves, which we apply to constrain the properties and mechanisms of the pre-explosion mass loss. We explore two limiting mass-loss scenarios: (1) an “eruption” arising from shock breakout following impulsive energy deposition below the stellar surface; and (2) a steady “wind,” due to sustained heating of the progenitor envelope. The eruption model, which resembles a scaled-down version of Type IIP SNe, can explain the luminosities and timescales of well-sampled precursors, for ejecta masses ∼ 0.1–1Mand velocities ∼ 100–1000 km s−1. By contrast, the steady wind scenario cannot explain the highest precursor luminosities ≳ 1041erg s−1, under the constraint that the total ejecta mass does not exceed the entire progenitor mass (though the less luminous SN 2020tlf precursor can be explained by a mass-loss rate ∼ 1Myr−1). However, shock interaction between the wind and pre-existing (earlier ejected) CSMmore »may boost its radiative efficiency and mitigate this constraint. In both the eruption and wind scenarios, the precursor ejecta forms compact (≲1015cm) optically thick CSM at the time of core collapse; though only directly observable via rapid post-explosion spectroscopy (≲ a few days before being overtaken by the SN ejecta), this material can boost the SN luminosity via shock interaction.

    « less
  3. Abstract

    We present Atacama Large Millimeter/submillimeter Array observations with a 800 au resolution and radiative-transfer modeling of the inner part (r≈ 6000 au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultracompact Hiiregion G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The Markov Chain Monte Carlo fits to the data give total stellar massesMfrom 120 to 200M, with much smaller ionized-gas massesMion-gas= 0.2–0.25M. The stellar mass is distributed within the gravitational radiusRg≈ 1000 to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation isi= 49°–56°. Radial motions at radiir>Rgconverge tovr,0≈ 8.7 km s−1, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars withMstar= 32–60M, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. Themore »inferred cluster density is very large, comparable to that reported at similar scales in the Galactic center. Stellar interactions are likely to occur within the next million years.

    « less
  4. Abstract

    For the first ∼3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium withn< 0.01 cm−3. We report on observational evidence for an excess of X-ray emission atδt> 900 days after the merger. WithLx≈ 5 × 1038erg s−1at 1234 days, the recently detected X-ray emission represents a ≥3.2σ(Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context ofJetFitafterglow models, current data represent a departure with statistical significance ≥3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the mergermore »ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.

    « less
  5. Abstract

    Gravitational-wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate,RBBH(z). We make predictions forRBBH(z) as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations,COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterized by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about 30Mand short delay times (tdelay≲ 1 Gyr), while the stable RLOF channel primarily forms systems with BH masses above 30Mand long delay times (tdelay≳ 1 Gyr). We provide a new fit for the metallicity-dependent specific star formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction ofRBBH(z). This leads to a distinct redshift evolution ofRBBH(z) for high and low primary BH masses. We furthermore find that, at high redshift,RBBH(z) is dominated by the CE channel, while at low redshift, it contains a large contribution (∼40%) from the stable RLOF channel.more »Our results predict that, for increasing redshifts, BBHs with component masses above 30Mwill become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution ofRBBH(z) for different BH masses can be tested with future detectors.

    « less