Arxiv:2005.06103V2
Total Page:16
File Type:pdf, Size:1020Kb
View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by MPG.PuRe Draft version July 10, 2020 Typeset using LATEX twocolumn style in AASTeX63 Intermediate Luminosity Red Transients by Black Holes Born from Erupting Massive Stars Daichi Tsuna,1, 2 Ayako Ishii,1, 3 Naoto Kuriyama,1, 4 Kazumi Kashiyama,1, 2 and Toshikazu Shigeyama1, 4 1Research Center for the Early Universe (RESCEU), Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 2Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 3Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Am M¨uhlenberg 1, Potsdam-Golm 14476, Germany 4Department of Astronomy, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan (Received |; Revised |; Accepted |) Submitted to ApJL ABSTRACT We consider black hole formation in failed supernovae when a dense circumstellar medium (CSM) is present around the massive star progenitor. By utilizing radiation hydrodynamical simulations, we calculate the mass ejection of blue supergiants and Wolf-Rayet stars in the collapsing phase and the radiative shock occurring between the ejecta and the ambient CSM. We find that the resultant emission is redder and dimmer than normal supernovae (bolometric luminosity of 1040{1041 erg s−1, effective temperature of ∼ 5 × 103 K, and timescale of 10{100 days) and shows a characteristic power- law decay, which may comprise a fraction of intermediate luminosity red transients (ILRTs) including AT 2017be. In addition to searching for the progenitor star in the archival data, we encourage X-ray follow-up observations of such ILRTs ∼ 1-10 yr after the collapse, targeting the fallback accretion disk. Keywords: black holes|high energy astrophysics; transient sources|high energy astrophysics 1. INTRODUCTION formation. During the protoneutron star phase before The main channel by which black holes (BHs) form BH formation, neutrinos carry away a significant frac- is considered to be the gravitational collapse of mas- tion of the energy of the core, around 10% of its rest sive stars. Numerical studies agree that a progenitor mass energy (e.g., O'Connor & Ott 2013). This results with a compact inner core can fail to revive the bounce in a decrease in the core's gravitational mass, and gen- shock (e.g. O'Connor & Ott 2011). These failed super- erates a sound pulse, which can eventually steepen into novae will be observed as massive stars suddenly vanish- a shock and unbind the outer envelope of the star upon ing upon BH formation, but their observational signa- shock breakout (Nadezhin 1980; Lovegrove & Woosley tures are not completely known. Previous studies imply 2013; Fern´andezet al. 2018; Coughlin et al. 2018a,b). that the outcome depends on the star's angular momen- This work investigates the emission from mass ejec- tum. Failed supernovae of stars with moderate or rapid tion of blue supergiant (BSG) and Wolf-Rayet (WR) rotation are expected to form accretion disks around progenitors that fail to explode. As shown in Figure1, arXiv:2005.06103v2 [astro-ph.HE] 9 Jul 2020 the BHs, from which energetic transients are generated we particularly consider emission from the interaction (Bodenheimer & Woosley 1983; Woosley 1993; Mac- between the ejecta and a dense circumstellar medium Fadyen & Woosley 1999; Kashiyama & Quataert 2015; (CSM), which is commonly introduced to explain Type Kashiyama et al. 2018; see also Quataert et al. 2019). IIn supernovae (e.g. Grasberg & Nadezhin 1986). In For the dominant slowly rotating case, the collapsing our case, where the ejecta mass is expected to be much star can still leave behind a weak transient after BH lighter than normal supernovae, CSM interaction is still, Corresponding author: Daichi Tsuna [email protected] 2 Tsuna et al. Weak explosion by Pre-explosion outburst Intermediate luminosity red transient? Fallback accretion disk? core neutrino emission ~ years before black hole formation after ~10-100 days after ~1-10 years Figure 1. Schematic figure of the emission we consider in this work. The ejecta created upon BH formation collides with the CSM created from a mass eruption ∼ years before core collapse. The kinetic energy of the ejecta is efficiently converted to radiation, being observable as intermediate luminosity red transients. The fallback of the outer layers of the envelope may form an accretion disk, that can be observable by X-rays 1{10 yr after core collapse. and even more, important to efficiently convert the ki- the huge mass loss implies energy injection greatly ex- netic energy into radiation 1. ceeding the Eddington rate, or instantaneous injection We find that the resulting emission is similar to what with timescales shorter than the outer envelope's dy- has been classified as intermediate luminosity red tran- namical timescale. sients (ILRTs) observed in the previous decades (Kulka- Kuriyama & Shigeyama(2020) studied the CSM re- rni et al. 2007; Botticella et al. 2009; Smith et al. 2009; sulting from energy injection, for various progenitors Berger et al. 2009; Bond et al. 2009; Cai et al. 2018; Jenc- while being agnostic of the injection model. A notable son et al. 2019; Williams et al. 2020; Stritzinger et al. finding was that the inner part has a profile of roughly 2020). The origin of ILRTs is unknown, with several ρ(r) / r−1:5, shallower than the commonly adopted interpretations such as electron-capture supernovae or wind profile (ρ(r) / r−2). An example is in Figure2, luminous blue variable-like mass eruptions. We propose where we plot the density profile of a BSG erupting a an intriguing possibility that BH formation of massive mass of 0:02 M at 8:6 yr before core collapse. Moti- stars can explain at least a fraction, if not all, of ILRTs. vated by this, we assume the CSM profile to be a power This Letter is constructed as follows. In Section 2 we law ρ(r) = qr−s, and take s = 1:5 as a representative present the details of our emission model and demon- value. The profile may be even shallower inwards due strate that an ILRT AT 2017be can be naturally ex- to the stronger gravitational pull from the central star. plained with our model. In Section 3 we estimate the However, our assumption should not severely affect the detectability of these signals by present and future op- dynamics of CSM interaction (given that the ejecta is tical surveys, and suggest ways to distinguish this from heavier than the inner CSM), since the CSM mass is other transients. concentrated at where s = 1:5. As neither the velocity nor mass-loss rate is constant, 2. OUR EMISSION MODEL the standard parameterization of using wind velocity 2.1. The Dense CSM and mass-loss rate is not appropriate. Instead, we pa- rameterize the dense CSM with its mass M , and Past observations of Type IIn SNe found that the mas- CSM time t between its eruption and core collapse. Since sive star's final years can be dramatic, with mass-loss CSM the outer edge of the CSM is v t , where v is the rates of 10−4 { 1 M yr−1 (Kiewe et al. 2012; Taddia out CSM out velocity of the outer edge, the CSM mass is related to et al. 2013). Such mass loss may be common also for q by q ≈ M (3 − s)(v t )s−3=4π. Because the Type IIP supernovae ( 70%; Morozova et al. 2018), CSM out CSM & CSM is created from the marginally bound part, v which comprises about half of all core-collapse super- out is of order the escape speed at the stellar surface. For novae. Although the detailed mechanism is unknown, s = 1:5, 1 For red supergiants, hydrogen recombination in the ejecta is ex- pected to dominantly power the emission (Lovegrove & Woosley 2013), and a strong candidate is already found (Adams et al. −1:5 2017). 7 MCSM vouttCSM q ∼ 10 cgs −2 15 : (1) 10 M 5 × 10 cm Red Transients from Newborn Black Holes 3 Pre-explosion Outburst Weak explosion by core neutrino emission 10 6 10 16 17 10 10 8 ) ) 18 3 3 10 19 m m 10 10 c 10 c g g 20 ( ( 10 y y t t 12 21 i i 10 10 s s n n 22 e e 10 d BSG-0.3-rad d 10 14 10 23 BSG model BSG-0.3 7 10 24 r 1.5(q = 5 × 106 cgs) v 10 16 7 8 9 1015 1016 10 10 10 radius (cm) velocity (cm s 1) Figure 2. (Left panel): Example density profile at core collapse of the unbound CSM that was erupted by energy injection 8.6 yr before core collapse. For this calculation, energy of 3:0 × 1048 erg is injected in the base of the envelope for a timescale of 4 −2 1:9 × 10 seconds, and the resulting total mass ejected is 1:9 × 10 M . (Right panel): Ejecta density profile for the \BSG{ 0.3" (light blue line, from the hydrodynamic simulation) and \BSG{0.3{rad" (dark blue line, from the radiation hydrodynamic simulation). We use the data at 105 s after core collapse, when the ejecta have become nearly homologous. Kuriyama & Shigeyama(2020) finds for their BSG and WR (their WR-1) models q ∼ 106{109 cgs and (n−3)=2 1=n 15 1 [2(5 − δ)(n − 5)Eej] vouttCSM ∼ 5 × 10 cm. g = ; (4) (n−5)=2 The optical depth is given as a function of radius 4π(n − δ) [(3 − δ)(n − 3)Mej] 1−s 1=2 by τ(r) ≈ κqr =(s − 1), where κ is the opacity.