FBOTs and AT2018cow following electron-capture collapse of merged white dwarfs
We suggest that fast-rising blue optical transients (FBOTs), and the brightest event of the class AT2018cow, result from electron-capture collapse following a merger of a massive ONeMg white dwarf (WD) with another WD. Two distinct evolutionary channels lead to the disruption of the less massive WD during the merger and formation of a shell burning non-degenerate star incorporating the ONeMg core. During the shell burning stage a large fraction of the envelope is lost to the wind, while mass and angular momentum are added to the core. As a result, the electron-capture collapse occurs with a small envelope mass, after years. During the formation of a neutron star as little as of the material is ejected at the bounce-off with mildly relativistic velocities and total energy few ergs. This ejecta becomes optically thin on time scales of days - this is the FBOT. During the collapse the neutron star is spun up and magnetic field is amplified. The ensuing fast magnetically-dominated relativistic wind from the newly formed neutron star shocks against the ejecta, and later against the wind. The radiation-dominated forward shock produces the long-lasting optical afterglow, while the termination shock of the relativistic wind produces the high energy emission in a Pulsar Wind Nebulae-like manner. The late X-ray variability is akin to magnetically-driven Crab flares (and possibly late GRB flares). If the secondary WD was of the DA type - the most frequent - the wind will have hydrogen, of the order of : this explains appearance of hydrogen late in the afterglow spectrum. The model explains many of the puzzling properties of FBOTs/AT2018cow: host galaxies, fast and light anisotropic ejecta producing bright optical peak, afterglow high energy emission of similar luminosity to optical, hard X-ray and infra-red features, presence of dense wind environment, late powerful radio emission.
AT2018cow (Smartt et al., 2018; Perley et al., 2018; Ho et al., 2018; Margutti et al., 2018) is a mysterious astrophysical event. It is likely to be the brightest member of the class of fast-rising blue optical transient (FBOT, Drout et al., 2014). AT2018cow seems to be at a cross-road of supernova explosion (and associated complicated nuclear reactions, neutrino transports physics), pulsars/magnetars, (early) pulsar wind nebulae (PWNe), possibly, GRBs and, possibly, white dwarf physics.
AT2018cow had a few surprising features, most important in our view being:
optical rise-time of days - this is order of magnitude shorter than the supernovae
peak optical luminosity erg s - this exceeds the typical peak power of supernovae
The X-ray emission of initial power erg s had an extra component at days, peaking at keV (Fig. 6 in Margutti et al., 2018)
There is a clear change of properties of the emission at 20 days.
There is an indication of the rising IR component at days (Fig. 5 in Perley et al., 2018)
There is bright radio emission day peaking at Hz (Fig. 11 in Margutti et al., 2018)
These properties exclude normal Ni-powered supernovae and require a separate formation channel. We discuss one such possible channel in the present paper, trying to build a coherent model of AT2018cow.
2 The model: a specific channel of WD mergers
In this paper we discuss a scenario where fast-rising blue optical transients (FBOTs) are powered by the electron-capture collapse following a merger of a massive ONeMg white dwarf (WD) with another WD. (Previously, the electron-capture collapse was mostly used in the accretion models, e.g., Canal & Schatzman, 1976, hence the name Accretion Induced Collapse, AIC.); see also Timmes & Woosley (1992); Miyaji et al. (1980); Nomoto & Kondo (1991); Kitaura et al. (2006); Jones et al. (2016) Some details of the binary evolution and of the collapse in such systems were previously discussed by Lyutikov & Toonen (2017), see also Brooks et al. (2017); Ruiter et al. (2018); Yungelson & Kuranov (2017).
Let us outline the main stages (see Lyutikov & Toonen, 2017, for a more detailed discussion); also Ruiter et al. (2018). An initial system with the primary mass and the secondary mass forms a double degenerate CO-ONeMg WD system, via two distinct evolutionary channels (Lyutikov & Toonen, 2017),. For sufficiently large mass ratio (Marsh et al., 2004) the ensuing gravitational wave-driven mass transfer is unstable, whereby the less massive CO WD is disrupted on a few orbital time scales and forms a disk around the primary. Disk accretion at high rates creates a spreading layer - a belt-like structure on the surface of the primary Inogamov & Sunyaev (1999); Balsara et al. (2009); Inogamov & Sunyaev (2010); Belyaev et al. (2013); Philippov et al. (2016). After the spreading is complete (on viscose time scale of seconds, e.g., Shen et al., 2012) the resulting star of consists of a slowly rotating degenerate ONeMg core, and a fast rotating, with period in the hundreds of seconds, non-degenerate envelope. After the removal of the degeneracy the envelope expands to few cm. The star will emit near the Eddington limit and drive powerful winds. Angular momentum contained in the shell will be both lost to the wind, and transported to the core through the (turbulent) boundary layer.
The merger product ignites shell CO burning, adding mass to the degenerate core; at the same time mass and angular momentum is lost due to powerful winds. For an ONeMg WD sufficiently close to the Chandrasekhar mass an electron-capture/accretion induced collapse follows after years later. During the collapse the magnetic field is amplified, Mösta et al. (2015) and the neutron star is spun to millisecond periods. As a result it produces long lasting relativistic wind that shocks against the ejecta material and, later on, against material lost during the shell-burning stage. The highly magnetized relativistic wind produced by a central neutron star will interact with the fairly dense newly ejected material and dense pre-AIC wind, producing -ray afterglow in the highly magnetized reverse shock, in a way similar to the case of afterglows from long GRBs, as suggested by Lyutikov & Camilo Jaramillo (2017).
In Fig. 1 we picture the immediate surrounding of a FBOT at times days-weeks after the collapse. This picture is our working model.
The key point of the model is that a merger of a heavy ONeMg WD with another WD, and ensuring mass loss during the shell burning stage, results in a collapse of the core surrounded by a fairly light envelope, tens of percent of at most. Depending on the particular collapse time one expects few cases when the AIC happens right before the shell burning is about to end. This would produce fast ejecta with small mass and leads to AT2018cow-like events. In other cases, envelopes of few tenths of solar masses are ejected, producing longer and less bright transients, which are still fast and luminous if compared with conventional supernovae.
In addition, to spinning up the core large amount of angular momentum of the shell (disrupted secondary WD) is lost to the wind. As a result the AIC can be direct, without formation of the accretion disk. In this case the newly born neutron star looses most if its rotational energy to the fast, relativistic wind. For shorter shell-burning stage, there is too much angular momentum in the shell, that leads to the formation of the accretion disk, that spins down the central neutron star in a propeller stage. As a result, a large fraction of the energy is deposited into relatively slow, matter-dominated wind with low radiation efficiency.
Regarding the origin of Type Ia SNe, our model is both somewhat independent of the active discussions on single versus double degenerate/mixed origin of SNIa (in a sense that it can work in both cases), but at the same time the model is closely related to the controversy. In any case, the merger of ONeMg-CO WDs will likely proceed in a regime different to the more common CO WDs merger.
3 The pre-collapse wind and the ejecta
3.1 Pre-collapse wind
Let’s assume that the post-merger/shell-burning star launches a wind with velocity . After time the edge of wind reaches radius . For a total mass in the wind we find
For mass loss rate of yr and total mass in the envelope of this stage lasts yrs; we adopt notation .
Optical depth trough the wind is at
As the shell is accreted onto the proto-neutron star, a narrow outer layer will be ejected. Previously, studies of the AIC ejected mass are predict ejected (Woosley & Baron, 1992; Fryer et al., 1999; Abdikamalov et al., 2010).
These studies were mostly concerned with single degenerate scenario of AIC, basically with no envelope. It is not clear how a presence of a tenuous envelope, of few , would affect the ejecta mass. Let us assume, for a particular case of AT2018cow, a bounced ejected mass of , with maximal velocity . For homologous expansion with , the energy in the ejected part is
Before the ejecta slows down due to the interaction with the pre-existing wind, its density evolves according to
where is the maximum velocity of the ejecta.
For numerical estimates we chose . (Many of the final relations depend weakly on ejecta mass, , e.g., Eq. (16) and 18)). The velocity of the ejecta will be related to the escape velocity from the proto-neutron star,
for a proto-neutron star radius of km. The ejecta energy is then erg. Thus, our model naturally produces fast and light ejecta.
The Thomson optical depth through ejecta is
where is time in days since the explosion. The surface evolves with time according to
It reaches maximum at
This explains the short rise time of the transient. Also note that radius (8) is larger than the radius (2) when the wind becomes optically thin. Thus, the maximum of the ejecta emission is not affected by scattering in the pre-explosion wind.
In fact, for the expanding optically thick plasma the peak in luminosity will appears before the shell becomes optically thin due to the fact that photons diffuse out somewhat faster (Arnett, 1982). The peak time is
The free-free emission of the optically thin part of the ejecta is fairly small,
where is free-free emissivity (Lang, 1999). Thus, we associated the early fast optical transient with the emission of the ejected shell.
3.3 Ejecta-wind interaction
The ejecta will launch a forward shock into the pre-explosion wind. As the shock propagates through the wind the accumulated mass is then
It becomes equal the ejecta mass at
Thus, the external wind has little effect on the ejecta until a approximately a month after the explosion. Before that the ejecta is in free expansion stage.
This will be approximately the time that the NS-driven shock exits the ejecta and enters the wind. We associate this with the transition of the afterglow properties at days (Margutti et al., 2018).
4 Non-thermal emission in FBOTs: pulsar-like termination shock in fast NS wind
4.1 Wind power
The newly created neutron star (NS) is spun up to short periods and magnetic field is amplified. The central NS will produce a highly magnetized wind that shocks against the ejecta (and, later on, against the pre-explosion wind). The NS wind-ejecta interaction will produce two shocks: forwards shock in the ejecta and termination shock in the wind. It is the wind termination shock that produces the X-ray emission, while the radiation-dominated forward shock produces the optical transient, see §4.4. In the termination shock the accelerated particles will produce synchrotron emission in the fast cooling regime, so that a large fraction of the wind power will be emitted as radiation. (See Lyutikov & Camilo Jaramillo, 2017, for discussion of emission produced at the highly magnetized termination shock in GRBs).
In the fast cooling regime most of the power given to the accelerated particles is emitted. Let us then identify the observed X -ray luminosity with NS wind power. It is nearly constant at erg s until 20 days and then falls (Fig. 5 in Margutti et al., 2018). Using the initial pulsar spin-down power and spin-down time , we find the magnetic field and the initial spin
At we have . Thus we suggest that a magnetar-type object is formed; its initial spin is not very high - milliseconds. The wind power at time is then
(In fact, since observed X-ray luminosity is a lower limit on wind power, estimates (13) are upper limit on magnetic field and lower limit on the initial spin.
4.2 Propagation of NS wind-driven shock through ejecta
The newly formed neutron star generates powerful wind that first propagates within the ejecta, and later on through the pre-explosion wind. Let us consider the dynamic of the wind-driven shock propagating through ejecta with density (4). In the Kompaneets approximation (Kompaneets, 1960) 111Kompaneets approximation assumes supersonic driving, while Sedov scaling assumes subsonic driving of the expansion. At early times, when the termination shock is close to the contact discontinuity, the Kompaneets approximation is more justified that the Sedov’s., the relativistic wind with power will produce a cavity expanding according to
(function has limits for and for ). For example, at one day the neutron star produces a cavity if a size .
The relative velocity of the shock with respect to the ejecta is
Thus, it changes only slowly with time.
The corresponding equipartition post-termination shock magnetic field (in the highly magnetized wind) is
( is time measured in days.) Thus, at early times magnetic field , while later .
4.3 Propagation of NS wind-driven shock through pre-explosion ejecta
Neglecting the slow-varying logarithm (the Kompaneets approximation is expected to less applicable at later times)
Thus, after break-out the expansion accelerates mildly, cf. Eq. (17).
The equipartition magnetic field is then
It is similar (both in value and in time-dependence) to the magnetic field before the NS wind breaks out of the ejecta.
4.4 Optical emission: forward shock from the NS-driven wind
The NS wind is shocked at the termination shock and will also produce a forward shock. Similar amount of energy will be dissipated in both shocks. Particles accelerated at the termination shock produces the non-thermal X-ray emission, while the forward shock will produce the long lasting optical emission (in addition to the emission produced by the ejecta).
Let us first consider the forward shock emission. The forward shock initially propagates through the ejecta, and later-on through the pre-explosion wind. The shock driven by the NS wind will be modified by radiation pressure, (Weaver, 1976; Budnik et al., 2010; Ofek et al., 2010; Lyutikov, 2018). The observed properties of mildly relativistic shocks are fairly complicated and often have steep dependence on the underlying parameters due to steep dependence of photon production rates, phone escape and pair production on the plasma properties. Qualitatively, radiation-dominated shocks can temporarily reach temperature exceeding the shock jump conditions (even taking into account radiation pressure). For fast photon production this may results in isothermal jumps (Landau & Lifshitz, 1959; Lyutikov, 2018), in which case the peak temperature exceeds the final temperature by a factor of a few. For slow photon production the temperature peak may exceed the final temperature by a large amount (Ito et al., 2018; Ioka et al., 2018).
Let us give here an order-of-magnitude estimate of the final post-shock temperature in strongly radiationally-dominated shock. The post shock pressure is
where and . Comparing with (17) we conclude that the forward shock is radiation-dominated with post-shock temperature (far downstream!)
This matches the observed temperatures at early times, both in value and (presumably) the temporal decrease (Margutti et al., 2018).
On the other hand, radiation-dominated shocks require sufficiently high optical depth, at least of the order of , while final stationary configurations may be reached at optical depths of thousands (Ioka et al., 2018). This is not really satisfied in the particular case: from , with given by (17) and density by (4), the condition is satisfied for
Thus, only at very early times the shock is highly radiation-dominated.
At the moment of shock breakout, at day, we expect an X-ray flash with duration seconds, Eq. (25) and luminosities erg s (e.g., Svirski & Nakar, 2014). Later on, after the shock is no longer radiation dominated, the post-shock temperature evaluates to
The corresponding free-free luminosity is far too small erg s.
Thus, optical emission is puzzling - it is hard to see how requirement of radiation-modified shocks (and hence large optical depth) can be reconciled with short transient time scales (and hence small optical depth). A possible answer is anisotropy. All the ingredients - pre-collapse wind, ejecta and the NS wind are expected to be anisotropic, see Fig. 1. First, the wind is launched by star that rotates with nearly critical velocity. Even Solar wind is highly anisotropic: slow and dense in the equatorial sector and fast in the polar regions (with approximately constant ). Second, the newly formed NS can, under certain parameters, be nearly critically spinning, so that the ejecta is anisotropic as well (see Lyutikov & Toonen, 2017). Third, the remaining NS’s wind is equatorially collimated, with power (Michel, 1973).
There are strong observational arguments in favor of anisotropy. First, hydrogen and helium lines show spectral asymmetry, with a tail towards longer wavelengths (Margutti et al., 2018). This can be explained if the line of sight samples unevenly the ejecta.
Anisotropic ejecta can also reconcile the requirements of low ejecta mass, (and hence early transparency) and the requirement of larger ejecta mass, to keep the NS wind-generated forward shock to be radaiatively dominated, which would nicely explain the optical temperature, Eq. (24).
An alternative possibility is that the later optical emission originates as a synchrotron emission from particles accelerated at the forward shock. One can construct a model of particle acceleration at the forward shock following the standard GRB parametrization (e.g., Piran, 2004), with , , in the fast cooling regime (fast cooling for the forward shock is important - X-ray luminosity for the termination shock and optical from the forward shock are comparable, (Margutti et al., 2018). But the expected spectrum will be non-thermal.
5 The X-ray transient: synchrotron emission from the NS-driven termination shock
5.1 Typical frequencies
We associate the X-ray emission with the particles accelerated at the termination shock. Let us consider its properties next. Suppose the neutron star-launched wind is propagating with Lorentz factor . The peak synchrotron frequency of particles accelerated by the wind termination shock is then, using the estimate of the magnetic field (18),
where is time since explosion in days. It is sharply decreases with time. Thus, after two weeks the injection peak passes through keV energy band. This is the hard X-ray component seen before 2 weeks, (Fig. 6 in Margutti et al., 2018)
At days . If at this moment the wind Lorentz factor falls to , then the injection falls in the near-IR. This explains the rising near-IR light curve after days (Fig. 5 in Perley et al., 2018).
The cooling energy
is well below - the emission is in the fast cooling regime (Sari et al., 1996). Thus, a large fraction of the energy given to the accelerated particles is radiated. The high radiation efficiency of the termination shock in the post-explosion wind was previously stressed by Lyutikov & Camilo Jaramillo (2017).
5.2 Evolution of particle distribution in expanding post termination shock flow
The magnetic field in the post termination shock flow, Eq. (18), decreases with time. As a result, particles will first cool quickly, but later on the radiative losses will decrease. This will lead to effective pile of cooled previously injected particles at some low Lorentz factor. We suggest this as an origin of late radio peak. Let us first consider the case of constant wind luminosity, . Neglecting the adiabatic losses (radiative losses are much more important), the evolution of the particles’ Lorentz factor is given by
At times , when the magnetic field evolves differently with time, Eq (18)
For particles injected at time with distribution function at , at later times the distribution function is given by
see Fig. (2)
At any moment the emission is contributed to by a distribution of particles with different injection energies, different injection times and different cooling history in the changing magnetic field. The outcome is bound to be complicated - the total distribution at any time is given by integration of (32) over and injection rate. Qualitatively, particles that are injected early on cool quickly and do not contribute much to the later emission. On the other hand, particles that are injected late, find themselves in low magnetic field and also do not emit efficiently.
As a simple order-of-magnitude estimate, we can assume that all particles injected at some cool quickly to the minimal energy , Eq (32). For the terminal Lorentz factors are
Then for constant injection rate we expected the distribution near the lowest energy range to be
This will produce fluxes and correspondingly. This exercise is not meant to produce a fit to data, but to demonstrate that cooling in decreasing magnetic field can produce rising spectral flux at the lower end of the spectrum.
6 Low frequency emission - free-free absorption in the ejecta and the wind
At lower frequencies, radio and IR waves can experience free-free absorption (Lang, 1999, Eq. 1.223) both within the ejecta and in the pre-explosion wind.
Ejecta contribute to free-free absorption a lot
It is optically thin for very high frequencies for a long time
Thus, in the radio and far IR the ejecta remains mostly opaque until the shock breakout from the ejecta, after approximately a month, Eq. (12).
The free-free optical depth through the wind, with density given by (1), becomes unity at
The effects of free-free absorption explain the evolution of the radio and IR luminosities, see (Ho et al., 2018, Fig. 1). High frequencies, 341 and 230 GHz, are transparent all along, while lower frequency, 34 traces the expanding radius of the corresponding surface.
7 Population synthesis
7.1 Pre-merger evolutionary channels
Most calculations of WD-WD mergers are aimed at explaining the Type Ia SNe, thus looking for detonation (see Maoz et al., 2014, for a recent review). Less attention has been given to models that fail to detonate. As we argue, failed SN Ia, that collapse via electron capture, may be related to the FBOTs. Dan et al. (2014) discussed the results of the WD-WD mergers and argued that there is large phase space available for WD-WD mergers to produce an accretion induced collapse (AIC). Nomoto & Iben (1985) stressed the role of carbon ignition during WD mergers in order to produce a Type Ia SN. Thus, in order to avoid explosion, there should be little carbon in the system. We suggest then that the primary is a heavy ONeMg WD. In this section we calculate possible evolutionary scenarios and rates for the corresponding mergers.
We use the binary population synthesis (BPS) method to predict the properties of the binary mergers, that is the merger rates, host galaxies and formation channels. Using the BPS code SeBa (Portegies Zwart & Verbunt, 1996; Toonen et al., 2012; Toonen & Nelemans, 2013), we simulate the evolution of a large number of binaries following in detail those that lead to the merger of an ONeMg and CO WD. Processes such as wind mass loss, stable & unstable mass transfer, accretion, angular momentum loss, and gravitational wave emission are taken into account. It was shown by Toonen et al. (2014) that the main source of uncertainty in the BPS outcomes come from the uncertainty in the input assumptions, in particular the CE-phase. For this reason, we follow Toonen et al. (2012), in performing two sets of population synthesis calculations using their model and . For full details on the models, see Toonen et al. (2012). In short, these models differ from one another with respect to the modeling of the CE-phase. Despite the importance of this phase for the formation of compact binaries and the enormous effort of the community, the CE-phase is still poorly constrained (see Ivanova et al., 2013, for a review). Commonly the CE-phase is modeled in BPS codes by energy conservation (Webbink, 1984), with a parameter that describes the efficiency with which orbital energy is consumed to unbind the CE. This recipe is used in model for every CE-phase. An alternative model has been proposed by (Nelemans et al., 2000) in order to reproduce the observed population of double white dwarfs. This model is based on a balance of angular momentum with an efficiency parameter . In our model , the -recipe is used unless the binary contains a compact object or the CE is triggered by a tidal instability (rather than dynamically unstable Roche lobe overflow, as proposed by Nelemans et al. (2000) for detached double white dwarfs.
Figures 3-4 show the initial parameters of binaries leading to mergers between ONeMg and CO WDs in our simulations. Every point represents a single system in the BPS simulations. The figures show that there are different evolutionary paths that can lead to an ONe-CO WD merger, however the dominant channels involve initially compact systems (in blue circles) i.e. the ’inverted channel’ and initially wide systems (in green squares) i.e the ’direct’ channel. For single stars, the initial mass of the progenitor of an ONeMg WD ranges between approximately 6.5-8M according to SeBa. This is similar to the range of initial masses in the direct channel where the primary forms the ONeMg WD (majority of green points in Fig. 3). The progenitors of ONeMg WDs in the ’inverted’ channel, i.e., the secondaries denoted in blue, have lower masses as these stars accrete a significant amount of mass from their companion stars.
In Fig. 5 we show the final masses of the ONeMg and CO WD that merge according to model and respectively. The masses of the ONeMg WDs are in the range , while the majority of CO WDs have masses in the range 0.5-0.8M. As described above it is possible that the ONeMg WD forms before the other WD in the system (channel ’direct’), or it forms afterwards (channel ’inverted’). In model , 48% of merging ONe-CO DWDs go through the ’direct’ channel, whereas for model the fraction goes up to 69%. The masses of the CO WDs in the ’inverted’ channel are systematically higher than those of the ’direct’ channel.
In Fig. 6 we show the distribution of the mass ratio as a function of the primary mass. For donor masses in the range Marsh et al. (2004) (see their Fig. 1) find that mass transfer is always unstable if the companion mass is above . It is always stable for . The blue systems are well above the limit for unstable mass transfer. The green systems occupy a larger part of parameter space. The far majority of the systems have a mass ratio that make stable mass transfer unlikely. Also note that given the ’optimistic’ stability limits of Marsh et al. (2004) the AM CVn rate is overestimated by orders of magnitude, indicating that mass transfer is likely less stable than the ’optimistic scenario’. In addition, the results from Marsh et al. (2004) do not take into account the effect of novae outbursts on the evolution of the systems. As shown by Shen et al. (2012) these outburst have a destabilizing effect on the mass transfer.
Using the BPS simulations, we also estimate the rate of CO-ONeMg WD mergers. Assuming a constant star formation history of yr for a Hubble time (roughly alike the Milky Way), the current merger rate ranges from yr for model and yr for model . This is in agreement with the BPS calculations of Yungelson & Kuranov (2017) and Ruiter et al. (2018). The CO-ONeMg are less common than mergers between CO-CO WDs for which we find a merger rate of yr. Mergers of CO-CO WDs with a combined mass above Chandrasekhar - the classical supernova type Ia progenitors (Webbink, 1984; Iben & Tutukov, 1984) have a merger rate of yr in our simulations, about an order of magnitude above that of CO-ONeMg WDs.
Integrated over time, the total number of CO-ONeMg WD mergers ranging between per Solar mass of created stars222This is independent of the assumed star formation history.. The rate of FBOTs has been estimated by Drout et al. (2014) to be 4%-7% of core-collapse supernova. Assuming the core-collapse rate is about 0.0025-0.010M (Maoz & Graur, 2017; Horiuchi et al., 2011; Graur et al., 2017)333The upper limit is calculated by Maoz & Graur (2017). The lower limit is based on the observed supernova type Ia rate (M Maoz & Graur, 2017) and the ratio of core-collapse to type Ia supernova (0.25-4 Graur et al., 2017), the estimated CO-ONeMg WD merger rate is consistent with the lower limit of the FBOT rate.
7.3 Host galaxies
In Fig. 7 we show the distributions of delay times of the CO-ONeMg mergers after a single burst of starformation. Their merger rates peaks at short delay times of about 50-100Myr, with a long tail to long delay times. The peak occurs significantly earlier than expected from the classical type Ia supernovae progenitors consisting of superChandrasekhar mergers of CO-CO WDs (consistent with Yungelson & Kuranov, 2017; Ruiter et al., 2018). The typical delay time of the CO-ONeMg mergers is closer to that of core-collapse supernovae, which peaks sharply at 50Myr (see e.g. Fig.3 of Zapartas et al., 2017). As a result, we expect the host galaxies of CO-ONeMg mergers to be more similar to those of core-collapse supernovae instead of Type Ia supernovae. This is consistent with the observed host galaxies of FBOTs (Drout et al., 2014).
An alternative to the WD merger scenario is the possibility of ultra-stripped envelopes in close binaries (Tauris et al., 2013, 2015). Binary interactions may strip the primary star of the envelope and also affect the mass of the collapsing core (Langer, 2012; Podsiadlowski et al., 2004). Müller et al. (2018a, b) estimate that ultra-stripped SNe produce normal slowly rotating pulsars, which are unlikely to produce fast spinning central neutron star. Müller et al. (2018a) studied the case when the angular momentum is implanted onto the neutron star only due to accretion of a tenuous envelope, without accounting for the progenitor’s initial spins. In contrast, in the case of double WD merger the envelope has a lot of angular momentum and has time to implant it to the core during the shell burning stage.
In a possibly related line of research, Piran & Shaviv (2005), see also Dall’Osso et al. (2014), argued that the pulsar J0737-3039B was born through non-standard SN mechanism (presumably via AIC), similar to the ultra-stripped case considered by Tauris et al. (2015). Piran & Shaviv (2005) also argued for slow initial spin (and slow kick velocity). Thus, ultra-stripped binary cores produce slowly rotating remnants, while AT2018cow needs a powerful central source.
In this paper we discuss a channel for transient emission after electron-capture collapse following a merger of two WDs. Qualitatively, this channels allows a collapse into a neutron star to occur with a small envelope mass. As a result the ejecta is light, have high velocity and becomes optically transparent much earlier. This early transparency allows higher radiation efficiency, as energy is not spent on adiabatic expansion of the envelope. In addition, light envelopes allow more efficient conversion of the rotational energy of the central engine - presumably a neutron star - into radiation. (If the envelope remains optically thick for along time any wind energy is spent on matter acceleration, not production of radiation. Also, large envelope can slow down the central neutron star in a propeller regime.)
Thus, small envelope mass, and small ejecta mass, lead to brighter and shorter SNe. In AT2018cow the ejecta was the lightest, with only ejected. In this case the AIC occurred at the time when most of the envelope was already lost to the wind. Other FBOTs may have larger remaining envelopes at the moment of AIC and larger ejected masses, but all smaller than . In our model the envelop mass depends on the mass of the primary ONeMg WD (how close it is to the Chandrasekhar limit), and the mass of the companion (how quickly mass is added to the core, mass loss rate - how long the shell burning continues).
We also advocate the Pulsar Wind paradigm - the reverse shock, as opposed to the forward shock - as the origin of the X-ray emission, similar to the GRB model of Lyutikov & Camilo Jaramillo (2017); see also Usov (1992); Lyutikov (2006b); Metzger et al. (2011).
Let us highlight how the key observational results discussed by Margutti et al. (2018) are explained in our model
A very short rise time to peak, few days - optical transient is generated by an envelope ejected during the bounce from the proto-neutron star. Ejected mass is small, while velocity is nearly relativistic (of the order of the escape velocity from the surface of a proto-neutron star).
Large bolometric peak luminosity - as the ejecta becomes optically thin early on, a large fraction of the internal energy is emitted (as opposed to been spent on adiabatic expansion in conventional SN explosion).
Persistent blue colors, with lack of evidence for cooling at days - later-on the emission starts to become dominated by the non-thermal particles accelerators at the termination shock.
Large blackbody radius cm - wind-driven cavity expands to these scales on time scales of few days, Eq. (16).
Persistent optically thick UV/optical emission with no evidence for transition into a nebular phase - emission is dominated by the radiation-modified forward shock.
Abrupt change of the velocity of the material which dominates the emission at times days - the jet breaks through the ejecta and enters the pre-explosion wind, after time given by (12).
NIR excess of emission - injection energy at the termination shock moved to IR both due to decreasing magnetic field in the post-shock flow and falling Lorentz factor of the wind, Eq, (27) and discussion afterwards.
The spectra evolve from a hot, blue, and featureless continuum around the optical peak, to very broad features - this is a transition from radiation-modified forward shock at early time to regular matter-dominated forward shock, combined with emission from the termination shock.
“Late-time optical spectra at days show line widths of 4000 km (, indicating substantially lower outflow velocities than at earlier times (when v = 0.1c), and an abrupt transition from very high velocity to lower velocity emitting material” (Margutti et al., 2018). We associate this transition with the moment that the NS-driven shock plows through the ejecta and enters the pre-collapse wind, Eq. (12) - this is due to slowing down of the ejecta.
There are indication of hydrogen later-on (Margutti et al., 2018) - if the disrupted WD was of the DA type - one expects some hydrogen in the pre-collapse wind, of the order of total. This explains the late appearance of hydrogen lines, presumably when the NS-driven forward shock exits the ejecta and enters the pre-existing wind. (We have no way of telling what the atmosphere of a post interaction WD will be. However, it is probable that the majority will be a DA, and therefore the rates would not be wildly different from the ones presented here.)
Similar X-ray and optical luminosities are naturally explained as emission from forward and NS wind termination shocks (the latter in the fast cooling regime).
Early X-ray spectral bump is also due to the passage of the peak frequency (late similar effect will produce IR increase)
Erratic inter-day variability of the X-ray emission (Ho et al., 2018) is hard to reproduce within the forward shock scenario (Lyutikov, 2006b), since the forward shock emission properties depend on the integrated quantities - central engine total energy and total matter swept. But they can be reproduced within the internal shock paradigm - in the fast cooling regime - either by the changes of the wind properties of the central source (Lyutikov & Camilo Jaramillo, 2017), minijets in the outflow (Lyutikov, 2006a), or by Crab flare-like reconnection processes in the shocked pulsar wind (Clausen-Brown & Lyutikov, 2012; Lyutikov et al., 2018).
Also note, that AIC with a formation of a neutron star is probably responsible for formation of young pulsars in globular clusters Lyne et al. (1996). This is consistent with the present model.
The present model, connecting FBOTs to the merger of WDs, is related to the possibility that some short GRBs come from a similar channel of WD mergers (Lyutikov & Toonen, 2017). The detection of gravitational waves associated with a GRB (Abbott et al., 2017) identifies merger of neutron stars as the central engine. It is not clear at the moment whether this identification is generic to the whole class of short GRBs. As discussed by Lyutikov (2009); Lyutikov & Toonen (2017), there is a number of observational contradictions to the binary neutron star merger paradigm (like extended emission and late flares - both not seen in GW/GRB170817). One possibility, that is still viable, is that some short GRBs originate from WD mergers. Several parameters may separate outcomes of WD mergers (e.g., FBOTs and short GRBs): masses of the merging WDs, the amount of the material lost to the wind, spin right before the AIC (so that AIC can either occur directly to a neutron star or with a formation of an accretion disk). For preferential intrinsic parameters and viewing angles (e.g., observer along the spin axis of the collapsing WD) we may see a short GRB and later on an FBOT.
This work had been supported by DoE grant DE-SC0016369 and NASA grant 80NSSC17K0757.
We would like to thank Maxim Barkov, Ori Fox, Raffaella Margutti, Danny Milisavljevic, Amir Levinson, Eran Ofek, Amir Sharon, Nir Shaviv and Thomas Tauris for discussions.
- Abbott et al. (2017) Abbott, B. P., Abbott, R., Abbott, T. D., Acernese, F., Ackley, K., Adams, C., Adams, T., Addesso, P., Adhikari, R. X., Adya, V. B., & et al. 2017, Physical Review Letters, 119, 161101
- Abdikamalov et al. (2010) Abdikamalov, E. B., Ott, C. D., Rezzolla, L., Dessart, L., Dimmelmeier, H., Marek, A., & Janka, H.-T. 2010, Phys. Rev. D, 81, 044012
- Arnett (1982) Arnett, W. D. 1982, ApJ, 253, 785
- Balsara et al. (2009) Balsara, D. S., Fisker, J. L., Godon, P., & Sion, E. M. 2009, ApJ, 702, 1536
- Belyaev et al. (2013) Belyaev, M. A., Rafikov, R. R., & Stone, J. M. 2013, ApJ, 770, 67
- Brooks et al. (2017) Brooks, J., Schwab, J., Bildsten, L., Quataert, E., Paxton, B., Blinnikov, S., & Sorokina, E. 2017, ApJ, 850, 127
- Budnik et al. (2010) Budnik, R., Katz, B., Sagiv, A., & Waxman, E. 2010, ApJ, 725, 63
- Canal & Schatzman (1976) Canal, R., & Schatzman, E. 1976, A&A, 46, 229
- Chevalier (1998) Chevalier, R. A. 1998, ApJ, 499, 810
- Clausen-Brown & Lyutikov (2012) Clausen-Brown, E., & Lyutikov, M. 2012, MNRAS, 426, 1374
- Dall’Osso et al. (2014) Dall’Osso, S., Piran, T., & Shaviv, N. 2014, MNRAS, 438, 1005
- Dan et al. (2014) Dan, M., Rosswog, S., Brüggen, M., & Podsiadlowski, P. 2014, MNRAS, 438, 14
- Drout et al. (2014) Drout, M. R., Chornock, R., Soderberg, A. M., Sanders, N. E., McKinnon, R., Rest, A., Foley, R. J., Milisavljevic, D., Margutti, R., Berger, E., Calkins, M., Fong, W., Gezari, S., Huber, M. E., Kankare, E., Kirshner, R. P., Leibler, C., Lunnan, R., Mattila, S., Marion, G. H., Narayan, G., Riess, A. G., Roth, K. C., Scolnic, D., Smartt, S. J., Tonry, J. L., Burgett, W. S., Chambers, K. C., Hodapp, K. W., Jedicke, R., Kaiser, N., Magnier, E. A., Metcalfe, N., Morgan, J. S., Price, P. A., & Waters, C. 2014, ApJ, 794, 23
- Fryer et al. (1999) Fryer, C., Benz, W., Herant, M., & Colgate, S. A. 1999, ApJ, 516, 892
- Graur et al. (2017) Graur, O., Bianco, F. B., Huang, S., Modjaz, M., Shivvers, I., Filippenko, A. V., Li, W., & Eldridge, J. J. 2017, ApJ, 837, 120
- Ho et al. (2018) Ho, A. Y. Q., Phinney, E. S., Ravi, V., Kulkarni, S. R., Petitpas, G., Emonts, B., Bhalerao, V., Blundell, R., Cenko, S. B., Dobie, D., Howie, R., Kamraj, N., Kasliwal, M. M., Murphy, T., Perley, D. A., Sridharan, T. K., & Yoon, I. 2018, ArXiv e-prints
- Horiuchi et al. (2011) Horiuchi, S., Beacom, J. F., Kochanek, C. S., Prieto, J. L., Stanek, K. Z., & Thompson, T. A. 2011, ApJ, 738, 154
- Iben & Tutukov (1984) Iben, Jr., I., & Tutukov, A. V. 1984, ApJS, 54, 335
- Inogamov & Sunyaev (1999) Inogamov, N. A., & Sunyaev, R. A. 1999, Astronomy Letters, 25, 269
- Inogamov & Sunyaev (2010) —. 2010, Astronomy Letters, 36, 848
- Ioka et al. (2018) Ioka, K., Levinson, A., & Nakar, E. 2018, arXiv e-prints
- Ito et al. (2018) Ito, H., Levinson, A., Stern, B. E., & Nagataki, S. 2018, MNRAS, 474, 2828
- Ivanova et al. (2013) Ivanova, N., Justham, S., Chen, X., De Marco, O., Fryer, C. L., Gaburov, E., Ge, H., Glebbeek, E., Han, Z., Li, X.-D., Lu, G., Marsh, T., Podsiadlowski, P., Potter, A., Soker, N., Taam, R., Tauris, T. M., van den Heuvel, E. P. J., & Webbink, R. F. 2013, A&A Rev., 21, 59
- Jones et al. (2016) Jones, S., Röpke, F. K., Pakmor, R., Seitenzahl, I. R., Ohlmann, S. T., & Edelmann, P. V. F. 2016, A&A, 593, A72
- Kitaura et al. (2006) Kitaura, F. S., Janka, H.-T., & Hillebrandt, W. 2006, A&A, 450, 345
- Kompaneets (1960) Kompaneets, A. S. 1960, Soviet Physics Doklady, 5, 46
- Landau & Lifshitz (1959) Landau, L. D., & Lifshitz, E. M. 1959, Fluid mechanics
- Lang (1999) Lang, K. R. 1999, Astrophysical formulae
- Langer (2012) Langer, N. 2012, ARA&A, 50, 107
- Lyne et al. (1996) Lyne, A. G., Manchester, R. N., & D’Amico, N. 1996, ApJ, 460, L41
- Lyutikov (2006a) Lyutikov, M. 2006a, MNRAS, 369, L5
- Lyutikov (2006b) —. 2006b, New Journal of Physics, 8, 119
- Lyutikov (2009) —. 2009, ArXiv e-prints
- Lyutikov (2018) —. 2018, MNRAS, 477, 816
- Lyutikov & Camilo Jaramillo (2017) Lyutikov, M., & Camilo Jaramillo, J. 2017, ApJ, 835, 206
- Lyutikov et al. (2018) Lyutikov, M., Komissarov, S., & Sironi, L. 2018, Journal of Plasma Physics, 84, 635840201
- Lyutikov & Toonen (2017) Lyutikov, M., & Toonen, S. 2017, ArXiv e-prints
- Maoz & Graur (2017) Maoz, D., & Graur, O. 2017, ApJ, 848, 25
- Maoz et al. (2014) Maoz, D., Mannucci, F., & Nelemans, G. 2014, ARA&A, 52, 107
- Margutti et al. (2018) Margutti, R., Metzger, B. D., Chornock, R., Vurm, I., Roth, N., Grefenstette, B. W., Savchenko, V., Cartier, R., Steiner, J. F., Terreran, G., Migliori, G., Milisavljevic, D., Alexander, K. D., Bietenholz, M., Blanchard, P. K., Bozzo, E., Brethauer, D., Chilingarian, I. V., Coppejans, D. L., Ducci, L., Ferrigno, C., Fong, W., GÖtz, D., Guidorzi, C., Hajela, A., Hurley, K., Kuulkers, E., Laurent, P., Mereghetti, S., Nicholl, M., Patnaude, D., Ubertini, P., Banovetz, J., Bartel, N., Berger, E., Coughlin, E. R., Eftekhari, T., Frederiks, D. D., Kozlova, A. V., Laskar, T., Svinkin, D. S., Drout, M. R., Macfadyen, A., & Paterson, K. 2018, ArXiv e-prints
- Marsh et al. (2004) Marsh, T. R., Nelemans, G., & Steeghs, D. 2004, MNRAS, 350, 113
- Metzger et al. (2011) Metzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, N., & Quataert, E. 2011, MNRAS, 413, 2031
- Michel (1973) Michel, F. C. 1973, ApJ, 180, L133
- Miyaji et al. (1980) Miyaji, S., Nomoto, K., Yokoi, K., & Sugimoto, D. 1980, PASJ, 32, 303
- Mösta et al. (2015) Mösta, P., Ott, C. D., Radice, D., Roberts, L. F., Schnetter, E., & Haas, R. 2015, Nature, 528, 376
- Müller et al. (2018a) Müller, B., Gay, D. W., Heger, A., Tauris, T. M., & Sim, S. A. 2018a, MNRAS, 479, 3675
- Müller et al. (2018b) Müller, B., Tauris, T. M., Heger, A., Banerjee, P., Qian, Y.-Z., Powell, J., Chan, C., Gay, D. W., & Langer, N. 2018b, ArXiv e-prints
- Nelemans et al. (2000) Nelemans, G., Verbunt, F., Yungelson, L. R., & Portegies Zwart, S. F. 2000, A&A, 360, 1011
- Nomoto & Iben (1985) Nomoto, K., & Iben, Jr., I. 1985, ApJ, 297, 531
- Nomoto & Kondo (1991) Nomoto, K., & Kondo, Y. 1991, ApJ, 367, L19
- Ofek et al. (2010) Ofek, E. O., Rabinak, I., Neill, J. D., Arcavi, I., Cenko, S. B., Waxman, E., Kulkarni, S. R., Gal-Yam, A., Nugent, P. E., Bildsten, L., Bloom, J. S., Filippenko, A. V., Forster, K., Howell, D. A., Jacobsen, J., Kasliwal, M. M., Law, N., Martin, C., Poznanski, D., Quimby, R. M., Shen, K. J., Sullivan, M., Dekany, R., Rahmer, G., Hale, D., Smith, R., Zolkower, J., Velur, V., Walters, R., Henning, J., Bui, K., & McKenna, D. 2010, ApJ, 724, 1396
- Perley et al. (2018) Perley, D. A., Mazzali, P. A., Yan, L., Cenko, S. B., Gezari, S., Taggart, K., Blagorodnova, N., Fremling, C., Mockler, B., Singh, A., Tominaga, N., Tanaka, M., Watson, A. M., Ahumada, T., Anupama, G. C., Ashall, C., Becerra, R. L., Bersier, D., Bhalerao, V., Bloom, J. S., Butler, N. R., Copperwheat, C., Coughlin, M. W., De, K., Drake, A. J., Duev, D. A., Frederick, S., Jesús González, J., Goobar, A., Heida, M., Ho, A. Y. Q., Horst, J., Hung, T., Itoh, R., Jencson, J. E., Kasliwal, M. M., Kawai, N., Kulkarni, S. R., Kumar, B., Kumar, H., Kutyrev, A. S., Khanam, T., Lee, W. H., Maeda, K., Mahabal, A., Murata, K. L., Neill, J. D., Ngeow, C.-C., Penprase, B., Pian, E., Quimby, R., Ramirez-Ruiz, E., Richer, M., Román-Zúñiga, C. G., Srivastava, S., Socia, Q., Sollerman, J., Tachibana, Y., Taddia, F., Tinyanont, K., Troja, E., Ward, C., & Wee, J. 2018, ArXiv e-prints
- Philippov et al. (2016) Philippov, A. A., Rafikov, R. R., & Stone, J. M. 2016, ApJ, 817, 62
- Piran (2004) Piran, T. 2004, Reviews of Modern Physics, 76, 1143
- Piran & Shaviv (2005) Piran, T., & Shaviv, N. J. 2005, Physical Review Letters, 94, 051102
- Podsiadlowski et al. (2004) Podsiadlowski, P., Langer, N., Poelarends, A. J. T., Rappaport, S., Heger, A., & Pfahl, E. 2004, ApJ, 612, 1044
- Portegies Zwart & Verbunt (1996) Portegies Zwart, S. F., & Verbunt, F. 1996, A&A, 309, 179
- Ruiter et al. (2018) Ruiter, A. J., Ferrario, L., Belczynski, K., Seitenzahl, I. R., Crocker, R. M., & Karakas, A. I. 2018, arXiv e-prints
- Sari et al. (1996) Sari, R., Narayan, R., & Piran, T. 1996, ApJ, 473, 204
- Shen et al. (2012) Shen, K. J., Bildsten, L., Kasen, D., & Quataert, E. 2012, ApJ, 748, 35
- Smartt et al. (2018) Smartt, S. J., Clark, P., Smith, K. W., McBrien, O., Maguire, K., O’Neil, D., Fulton, M., Magee, M., Prentice, S., Colin, C., Tonry, J., Denneau, L., Stalder, B., Heinze, A., Weiland, H., Flewelling, H., & Rest, A. 2018, The Astronomer’s Telegram, 11727
- Svirski & Nakar (2014) Svirski, G., & Nakar, E. 2014, ApJ, 788, 113
- Tauris et al. (2013) Tauris, T. M., Langer, N., Moriya, T. J., Podsiadlowski, P., Yoon, S.-C., & Blinnikov, S. I. 2013, ApJ, 778, L23
- Tauris et al. (2015) Tauris, T. M., Langer, N., & Podsiadlowski, P. 2015, MNRAS, 451, 2123
- Timmes & Woosley (1992) Timmes, F. X., & Woosley, S. E. 1992, ApJ, 396, 649
- Toonen et al. (2014) Toonen, S., Claeys, J. S. W., Mennekens, N., & Ruiter, A. J. 2014, A&A, 562, A14
- Toonen & Nelemans (2013) Toonen, S., & Nelemans, G. 2013, A&A, 557, A87
- Toonen et al. (2012) Toonen, S., Nelemans, G., & Portegies Zwart, S. 2012, A&A, 546, A70
- Usov (1992) Usov, V. V. 1992, Nature, 357, 472
- Weaver (1976) Weaver, T. A. 1976, ApJS, 32, 233
- Webbink (1984) Webbink, R. F. 1984, ApJ, 277, 355
- Weiler et al. (1986) Weiler, K. W., Sramek, R. A., Panagia, N., van der Hulst, J. M., & Salvati, M. 1986, ApJ, 301, 790
- Woosley & Baron (1992) Woosley, S. E., & Baron, E. 1992, ApJ, 391, 228
- Yungelson & Kuranov (2017) Yungelson, L. R., & Kuranov, A. G. 2017, MNRAS, 464, 1607
- Zapartas et al. (2017) Zapartas, E., de Mink, S. E., Izzard, R. G., Yoon, S.-C., Badenes, C., Götberg, Y., de Koter, A., Neijssel, C. J., Renzo, M., Schootemeijer, A., & Shrotriya, T. S. 2017, A&A, 601, A29