Arxiv:1905.05573V3 [Astro-Ph.HE] 5 Aug 2019 O H Xlso Oe Htte Assume, They That Model Explosion the for Iharltvl Atoto of (CSM) Outflow Matter Fast Curve
Total Page:16
File Type:pdf, Size:1020Kb
Draft version August 6, 2019 Typeset using LATEX twocolumn style in AASTeX62 The strongly interacting binary scenarios of the enigmatic supernova iPTF14hls Roni Anna Gofman1 and Noam Soker1,2 1Department of Physics, Technion, Haifa, 3200003, Israel; [email protected]; [email protected] 2Guangdong Technion Israel Institute of Technology, Shantou 515069, Guangdong Province, China ABSTRACT We argue that some scenarios for the enigmatic supernova (SN) iPTF14hls and its progenitor require a strong binary interaction. We examine scenarios that attribute the extra power of iPTF14hls to a magnetar, to a late fallback on to the neutron star (NS) that launches jets, to an interaction of the ejecta with a circumstellar matter (CSM), or to a common envelope jets SN (CEJSN). For each of these four scenarios, we study the crucial process that supplies the extra energy and conclude that a binary companion to the progenitor must be present. For the magnetar scenario and late jets we claim that a companion should spin-up the pre-collapse core, in the ejecta-CSM scenario we find that the formation of the equatorial CSM requires a companion, and in the CEJSN where a NS spirals-in inside the giant envelope of the progenitor and launches jets the strong binary interaction is built-in. We argue that these types of strong binary interactions make the scenarios rare and explain the enigmatic nature of iPTF14hls. We further study processes that might accompany the binary interaction, in particular, the launching of jets before, during and after the explosion and their observational consequences. We do not consider the difficulties of the different scenarios and neither do we determine the best scenario for iPTF14hls. We rather focus on the binary nature of these scenarios that might as well explain other rare types of SNe. Keywords: Supernovae — stars: jets — stars: variables: general — binaries: general 1. INTRODUCTION (2017) find the kinetic energy of the absorbing gas to be E 1052 erg. However, some scenarios for Arcavi et al. (2017) report the discovery and evo- CSM ≈ lution and Sollerman et al. (2019) describe the late iPTF14hls do not need such a large kinetic energy, e.g., evolution of the extraordinary Type II supernova Andrews & Smith (2017). (SN) iPTF14hls (AT 2016bse; Gaia16aog). It is even Sollerman et al. (2019) estimate that iPTF14hls emit- ted E = 3.6 1050 erg from discovery to their last not clear whether iPTF14hls is a core collapse su- rad × pernova (CCSN), a pair instability supernova, or a observation. As we argue later, this energetic emis- common envelope jets supernova (CEJSN). It might sion hints that iPTF14hls could not have been driven turn out that iPTF14hls is not an extremely rare by neutrinos, as the explosion energy most likely was E > 2 1051 erg. Sollerman et al. (2019) argue that type of a SN, as Arcavi et al. (2018) suggest that exp × SN 2018aad (ASASSN-18eo; Hosseinzadeh et al. 2018; the fast luminosity decline in the third year excludes the arXiv:1905.05573v3 [astro-ph.HE] 5 Aug 2019 Nicholls et al. 2018) is similar in many aspects to late mass accretion scenario and the magnetar scenario iPTF14hls. Milisavljevic & Margutti (2018) in their for iPTF14hls. They find the ejecta-CSM interaction review of peculiar SNe mention that some SNe that are to better fit their observations, but like Woosley (2018) initially classified as peculiar are later incorporated into they find that none of the scenarios can fit all properties the spectrum of standard events. of iPTF14hls. The peculiar properties of iPTF14hls include the fol- The peculiar nature of iPTF14hls motivated sev- lowing. (1) The light curve evolution is about an or- eral theoretical scenarios. Examples include scenarios der of magnitude slower than that of typical type II-P that build on a rapidly rotating magnetic neutron star SNe. (2) There are at least five peaks in the light curve. (NS), i.e., a magnetar (e.g., Arcavi et al. 2017; Dessart (3) There is an absorbing circumstellar matter (CSM) 2018; Woosley 2018), some that attribute the pro- with a relatively fast outflow of v 6000 km s−1. longed powering to fallback accretion (e.g., Arcavi et al. CSM ≈ For the explosion model that they assume, Arcavi et al. 2017; Wang et al. 2018; Liu et al. 2019) or interac- tion of the ejecta with a circumstellar matter (CSM; 2 Andrews & Smith 2017), while others explicitly dis- tar to the CEJSN scenario. In section 6 we summarise cuss a binary companion (e.g., Soker & Gilkis 2018). A our conclusion that some of the scenarios for peculiar different class of models consider the possibility that SNe that might account also for iPTF14hls require a iPTF14hls is a pair instability supernova (e.g., Woosley strong binary interaction. 2018; Vigna-G´omez et al. 2019). Jets play significant roles in some scenarios, either 2. THE MAGNETAR SCENARIO only in the explosion itself (e.g., Chugai 2018), or both in 2.1. The properties of the magnetar the explosion and a long-lasting powering of iPTF14hls, As Arcavi et al. (2017) already mentioned, one pos- as in the CEJSN scenario (Soker & Gilkis 2018) and in sible long-lasting power source for iPTF14hls is the the scenario of late accretion of hydrogen-rich gas with rotational energy of the central NS as it spins down stochastic angular momentum (Quataert et al. 2019), or by dipole radiation, i.e., a magnetar. Let RNS be only in a late powering. Liu et al. (2019) attribute the the radius of the magnetar, INS its moment of in- prolonged activity of iPTF14hls to the launching of jets ertia, B its magnetic field, and P its rotation pe- by a black hole (BH) that continues to accrete mass for riod. The initial rotational energy of such a magnetar a long time, hundreds of days and more. 50 45 2 −2 is Ep 2 10 INS/10 g cm (P/10 ms) erg. As well, it might be that in all cases where magne- ≈ × Its magnetic dipole radiation power is Pp 4.54 tars supply a non-negligible amount of energy, jets play 50 14 2 6 −4 ≈ −1 × 10 (B/10 G) (RNS/12 km) (P/10 ms) erg yr , a more significant role than the magnetar in the total and it can power a supernova to have a maximum lumi- energy budget (Soker 2016a; Soker & Gilkis 2017a). nosity of (e.g., Kasen, & Bildsten 2010) Wang et al. (2018) propose that intermittent fallback accretion of 0.2M⊙ powered iPTF14hls. In their sce- E t t t ≈ L f p p ln 1+ d d , (1) nario the total ejecta mass in the explosion is about peak 2 ≈ td tp − td + tp 51 21M⊙ and the explosion energy is 2.2 10 erg, that × they argue can be driven by the delayed neutrino mech- where f is a correction parameter, td is the photon dif- anism. They do not mention any binary model, and it fusion time, and tp is the spin-down time scale of the is not clear why this explosion should be different from magnetar. By comparing eq. (1) to numerical simula- other CCSNe. They could not fit the third peak in the tions Kasen, & Bildsten (2010) find f =3/2. The diffu- sion time depends on the the ejecta properties according light curve with a fallback model, and attribute it to a 1/2 magnetic outburst on the central NS. to td = (3κMej/4πvf c) where κ is the ejecta opacity, In the present study we examine some aspects of sce- Mej is the ejecta mass, vf = (Ep + ESN) /2Mej is the final characteristic ejecta velocity, and E is the kinetic narios for iPTF14hls that, we argue, require strong p SN stellar binary interaction. We consider only cases energy of the SN explosion itself. The spin-down time where the explosion energy comes from gravitational scale is given by energy, and we do not consider thermonuclear mod- − E P 2 B 2 els. Vigna-G´omez et al. (2019) claim for pair instabil- t = p 0.44 p P ≈ 10 ms 1014 G ity supernova resulting from the merger of two giant p −4 (2) stars of about equal mass. This is another scenario for INS RNS yr. iPTF14hls of a strongly interacting binary system, but × 1045 g cm2 12 km we do not consider it in the present study. At late times after peak luminsoity the SN light curve We do not rank the scenarios we study by how well follows the spin-down power L (t)= E t−1 (1 + t/t )−2. they fit the observations of iPTF14hls. We limit our- p p p p By fitting the late time light curve of iPTF14hls to selves to examine the implications of the scenarios on L (t) Arcavi et al. (2017) find that a magnetar with binary interaction, as even if a scenario does not fit p intial spin period of P = 5 ms and an initial magnetic iPTF14hls it might fit newly discovered peculiar SNe field of B = 0.7 1014 G can produce the observed in the future. × average luminosity and timescale. In section 2 we discuss the magnetar model that Dessart (2018) shows that a magnetar with P 7 ms, Woosley (2018) proposes, in section 3 we discuss the ≈ I = 1045 g cm2, R = 10 km and B =0.7 1014 G shocked CSM scenario that Andrews & Smith (2017) NS NS × that powers a typical Type II SN ejecta with Mej = propose and which Milisavljevic & Margutti (2018) sup- 51 13.35M⊙ and E = 1.32 10 erg can produce most port based on new observations (unpublished), and in SN × of the observed properties of SN iPTF14hls.