Many core-collapse supernova (SN) progenitors show indications of enhanced pre-SN mass loss and outbursts, some of which could be powered by wave energy transport within the progenitor star. Depending on the star's structure, convectively excited waves driven by late-stage nuclear burning can carry substantial energy from the core to the envelope, where the wave energy is dissipated as heat. We examine the process of wave energy transport in single-star SNe progenitors with masses between 11 and 50 M e . Using MESA stellar evolution simulations, we evolve stars until core collapse and calculate the wave power produced and transmitted to the stars' envelopes. These models improve upon prior efforts by incorporating a more realistic wave spectrum and nonlinear damping effects, reducing our wave-heating estimates by ∼1 order of magnitude compared to prior work. We find that waves excited during oxygen/neon burning typically transmit ∼10 46 -10 47 erg of energy at 0.1-10 yr before core collapse in typical (M<30 M e ) SN progenitors. High-mass progenitors can often transmit ∼10 47 -10 48 erg of energy during oxygen/neon burning, but this tends to occur later, at about 0.01-0.1 yr before core collapse. Pre-SN outbursts may be most pronounced in low-mass SN progenitors (M12 M e ) undergoing semidegenerate neon ignition and in high-mass progenitors (M30 M e ) exhibiting convective shell mergers.

After the tidal disruption of a star by a massive black hole, disrupted stellar debris can fall back to the hole at a rate significantly exceeding its Eddington limit. To understand how black hole mass affects the duration of super-Eddington accretion in tidal disruption events, we first run a suite of simulations of the disruption of a Solarlike star by a supermassive black hole of varying mass to directly measure the fallback rate onto the hole, and we compare these fallback rates to the analytic predictions of the "frozen-in" model. Then, adopting a Zero-Bernoulli Accretion flow as an analytic prescription for the accretion flow around the hole, we investigate how the accretion rate onto the black hole evolves with the more accurate fallback rates calculated from the simulations. We find that numerically-simulated fallback rates yield accretion rates onto the hole that can, depending on the black hole mass, be nearly an order of magnitude larger than those predicted by the frozen-in approximation. Our results place new limits on the maximum black hole mass for which super-Eddington accretion occurs in tidal disruption events.

General-relativistic magnetohydrodynamic (GRMHD) simulations are an important tool to study a variety of astrophysical systems such as neutron star mergers, core-collapse supernovae, and accretion onto compact objects. A conservative GRMHD scheme numerically evolves a set of conservation equations for 'conserved' quantities and requires the computation of certain primitive variables at every time step. This recovery procedure constitutes a core part of any conservative GRMHD scheme and it is closely tied to the equation of state (EOS) of the fluid. In the quest to include nuclear physics, weak interactions, and neutrino physics, state-of-the-art GRMHD simulations employ finitetemperature, composition-dependent EOSs. While different schemes have individually been proposed, the recovery problem still remains a major source of error, failure, and inefficiency in GRMHD simulations with advanced microphysics. The strengths and weaknesses of the different schemes when compared to each other remain unclear. Here we present the first systematic comparison of various recovery schemes used in different dynamical spacetime GRMHD codes for both analytic and tabulated microphysical EOSs. We assess the schemes in terms of (i) speed, (ii) accuracy, and (iii) robustness. We find large variations among the different schemes and that there is not a single ideal scheme. While the computationally most efficient schemes are less robust, the most robust schemes are computationally less efficient. More robust schemes may require an order of magnitude more calls to the EOS, which are computationally expensive. We propose an optimal strategy of an efficient three-dimensional Newton-Raphson scheme and a slower but more robust one-dimensional scheme as a fall-back. * NASA Einstein Fellow arXiv:1712.07538v2 [astro-ph.HE]

The discovery of rapidly rising and fading supernovae powered by circumstellar interaction has suggested the pre-supernova mass eruption phase as a critical phenomenon in massive star evolution. It is important to understand the mass and radial extent of the circumstellar medium (CSM) from theoretically predicted mass ejection mechanisms. In this work, we study the wave heating process in massive hydrogen-poor stars, running a suite of stellar models in order to predict the wave energy and pre-explosion timescale of surface energy deposition. We survey stellar models with main-sequence progenitor masses from 20–70 M ⊙ and metallicity from 0.002–0.02. Most of these models predict that less than ∼1047 erg is deposited in the envelope, with the majority of the energy deposited in the last week of stellar evolution. This translates to CSM masses less than ∼10−2 M ⊙ that extend to less than ∼1014 cm, too small to greatly impact the light curves or spectra of the subsequent supernovae, except perhaps during the shock breakout phase. However, a few models predict somewhat higher wave energy fluxes, for which we perform hydrodynamical simulations of the mass ejection process. Radiative transfer simulations of the subsequent supernovae predict a bright but brief shock-cooling phase that could be detected in some Type Ib/c supernovae if they are discovered within a couple days of explosion.

Most massive stars exchange mass with a companion, leading to evolution which is altered drastically from that expected of stars in isolation. Such systems result from unusual binary evolution pathways and can place stringent constraints on the physics of these interactions. We use the R4 binary system's B[e] supergiant, which has been postulated to be the product of a stellar merger, to guide our understanding of such outcomes by comparing observations of R4 to the results of simulating a merger with the 3D hydrodynamics code FLASH. Our approach tailors the simulation initial conditions to observed properties of R4 and implements realistic stellar profiles from the 1D stellar evolution code MESA onto the 3D grid, resolving the merger inspiral to within 0.02 R e. We map the merger remnant into MESA to track its evolution on the H-R diagram over a period of 10 4 yr. This generates a model for a B[e] supergiant with stellar properties, age, and nebula structure in qualitative agreement with those of the R4 system. Our calculations provide evidence to support the idea that R4ʼs B[e] supergiant was originally a member of a triple system in which the inner binary merged after its most massive member evolved off the main sequence, producing a new object of similar mass but significantly more luminosity than the A supergiant companion. The code framework presented in this paper, which was constructed to model tidal encounters, can be used to generate accurate models of a wide variety of merger stellar remnants.

Many core-collapse supernovae (SNe) with hydrogen-poor and low-mass ejecta, such as ultra-stripped SNe and type Ibn SNe, are observed to interact with dense circumstellar material (CSM). These events likely arise from the core collapse of helium stars that have been heavily stripped by a binary companion and have ejected significant mass during the last weeks to years of their lives. In helium star models run to days before core collapse we identify a range of helium core masses ≈2.5–3 M ⊙ whose envelopes expand substantially due to the helium shell burning while the core undergoes neon and oxygen burning. When modeled in binary systems, the rapid expansion of these helium stars induces extremely high rates of late-stage mass transfer ( M ̇ ≳ 10 − 2 M ⊙ yr − 1 ) beginning weeks to decades before core collapse. We consider two scenarios for producing CSM in these systems: either mass transfer remains stable and mass loss is driven from the system in the vicinity of the accreting companion, or mass transfer becomes unstable and causes a common envelope event (CEE) through which the helium envelope is unbound. The ensuing CSM properties are consistent with the CSM masses (∼10−2–1 M ⊙) and radii (∼1013–1016 cm) inferred for ultra-stripped SNe and several type Ibn SNe. Furthermore, systems that undergo a CEE could produce short-period neutron star binaries that merge in less than 100 Myr.

The ∼100 tidal disruption events (TDEs) observed so far exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs radiate predominantly at X-ray energies, while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. In particular, for optical TDEs, while their UV and optical emissions decline somewhat following the fallback pattern, some events can greatly rebrighten in X-rays at late time. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates in the regime of super-Eddington accretion. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al., which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation. Furthermore, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect that more X-rays will leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.