An official website of the United States government
Abstract We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25M⊙star in the early phase of the main sequence. We use the explicit gas-dynamics codePPMstar, which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. We perform very long simulations on 8963grids accelerated by luminosity boost factors of 1000, 3162 and 10,000. In these simulations, the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a 1D method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. The 1D model is based on a spherically averaged reduced entropy equation that takes the turbulent dissipation as input from the 3D hydrodynamic simulation and takes buoyancy and all other energy sources and sinks into account. This 1D method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale.
more »
« less
Shock Breakout in Three-dimensional Red Supergiant Envelopes
Abstract 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 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.
more »
« less
- Award ID(s):
- 1663688
- PAR ID:
- 10368997
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 933
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 164
- Size(s):
- Article No. 164
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT 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. 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.more » « less
-
Abstract We explore the three-dimensional properties of convective, luminous (L≈ 104.5–105L⊙), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry usingAthena++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3M⊙of the simulatedM≈ 15M⊙stars. This work reveals a radius,Rcorr, around which the nature of the convection changes. Forr>Rcorr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well studied in less luminous stars, but in RSGs, the near- (or above-)Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outward rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and early lightcurves of Type IIP supernovae. Inside ofRcorr, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length theory (MLT). Radiation pressure provides ≈1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length ofα= 3–4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes (≈300 days), with amplitudes consistent with those observed photometrically.more » « less
-
Abstract Superluminous supernovae (SLSNe) radiate ≳10–100 times more energy than ordinary stellar explosions, implicating a novel power source behind these enigmatic events. One frequently discussed source, particularly for hydrogen-poor (Type I) SLSNe, is a central engine such as a millisecond magnetar or accreting black hole. Both black hole and magnetar engines are expected to channel a fraction of their luminosity into a collimated relativistic jet. Using 3D relativistic hydrodynamical simulations, we explore the interaction of a relativistic jet, endowed with a luminosityLj≈ 1045.5erg s−1and durationteng≈ 10 days compatible with those needed to power SLSNe, launched into the envelope of the exploding star. The jet successfully breaks through the expanding ejecta, and its shocked cocoon powers ultraviolet/optical emission lasting several days after the explosion and reaching a peak luminosity ≳1044erg s−1, corresponding to a sizable fraction ofLj. This high radiative efficiency is the result of the modest adiabatic losses the cocoon experiences owing to the low optical depths of the enlarged ejecta at these late times, e.g., compared to the more compact stars in gamma-ray bursts. The luminosity and temperature of the cocoon emission match those of the “bumps” in SLSN light curves observed weeks prior to the optical maximum in many SLSNe. Confirmation of jet breakout signatures by future observations (e.g., days-long to weeks-long internal X-ray emission from the jet for on-axis observers, spectroscopy confirming large photosphere velocitiesv/c≳ 0.1, or detection of a radio afterglow) would offer strong evidence for central engines powering SLSNe.more » « less
-
Abstract We present UV and/or optical observations and models of SN 2023ixf, a type II supernova (SN) located in Messier 101 at 6.9 Mpc. Early time (flash) spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of Hi, Hei/ii, Civ, and Niii/iv/vwith a narrow core and broad, symmetric wings arising from the photoionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for ∼8 days with respect to first light, at which time Doppler broadened the features from the fastest SN ejecta form, suggesting a reduction in CSM density atr≳ 1015cm. The early time light curve of SN 2023ixf shows peak absolute magnitudes (e.g.,Mu= −18.6 mag,Mg= −18.4 mag) that are ≳2 mag brighter than typical type II SNe, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light-curve and multiepoch spectral models from the non-LTE radiative transfer codeCMFGENand the radiation-hydrodynamics codeHERACLESsuggests dense, solar-metallicity CSM confined tor= (0.5–1) × 1015cm, and a progenitor mass-loss rate of yr−1. For the assumed progenitor wind velocity ofvw= 50 km s−1, this corresponds to enhanced mass loss (i.e.,superwindphase) during the last ∼3–6 yr before explosion.more » « less
