arXiv:2002.11975v2 [astro-ph.HE] 9 Apr 2020 rglci ngas a ( far so enigma tragalactic 1 2 2013 2016 2020 ( bursts repeating 2016 into divided nomenally ( rarely peaks 2016 multiple or mainly luminosity peak peak the with of pulses limits radio upper coherent average millisecond-duration are by discovered were they since 2018 ua ta.2019 al. et Kumar al. et Keane atojcshv enpooe ( com- proposed in involved been have models literature progenitor objects the of pact dozens in Catalogue FRB ena, reported the been in collected have and FRBs 100 × × 1 atrdobrt FB)hv eanda ex- an remained have (FRBs) bursts radio Fast yee sn L using 2020 Typeset 10, April version Draft http://www.frbcat.org 10 2 1 10 .Mawie oepanti ai phenom- radio ; this explain over to Meanwhile, date, To ). bursts. non-repeating and ) ; ; ; e aoaoyo oenAtooyadAtohsc (Nanji Astrophysics and Astronomy Modern of laboratory Key colo srnm n pc cec,NnigUniversity, Nanjing Science, Space and Astronomy of School 42 aa ta.2018 al. et Farah 42 erffe l 2019 al. et Petroff ihnalrernei oa niomnsadhs galaxies. host Keywords: F and environments and local 180924-like in FRB range of large properties a NSWD within observational This the channels. hos both BNS/BWD/AIC 180924- and of explain FRB environments, those of local resemble subset mergers ejecta, a NSWD the least since at me magnetars, NSWD produce flaring If a could field. rotation objects fossil differential such a to of then scenario t due conservation in mechanism flux result magnetic dynamo could the the transfer via mass field unstable with magnetic mergers NSWD that find NW)mres oiae yFB102-iefs ai uss( bursts radio by fast suggested channels 180924-like BNS/BWD/AIC in FRB prod born be by also motivated could mergers, magnetars whether (NSWD) investigate accretion-ind we or paper, mergers, this (BWD) dwarf In white binary or (BNS) star antr rmNurnSa–ht wr egr:Applicat Mergers: Dwarf Star–White Neutron from Magnetars ace&Rzol 2014 Rezzolla & Falcke − HM/R olbrto ta.2019a al. et Collaboration CHIME/FRB r ( erg ti ieyblee htmgeascudb oni oeclas sup core-collapse in born be could magnetars that believed widely is It 8 × ( 2012 hn 2018 Zhang 10 A T 1. E 44 X ,and ), (1339) ; opc iaysas(8) rvttoa ae 68;Magnet (678); waves Gravitational (283); stars binary Compact INTRODUCTION r s erg twocolumn HM/R olbrto tal. et Collaboration CHIME/FRB ; − ; hrtne al. et Thornton rcak ta.2019 al. et Prochaska 1 ,caatrzdb single a by characterized ), ods&Catre 2019 Chatterjee & Cordes n energy and az2018 Katz tl nAASTeX63 in style oie tal. et Lorimer ; ybrk 2014 Lyubarsky ahym tal. et Kashiyama h-igZhong Shu-Qing hmine al. et Champion 1 E ; ( ( ∼ 2013 erffe al. et Petroff pte tal. et Spitler oo tal. et Popov 7 × .They ). ,phe- ), ( 10 L 2007 p 39 ABSTRACT , b ∼ − ), ) aglte al. et Margalit ; ; ,2 1, ajn 103 hn;szoghtalcm [email protected] [email protected], China; 210093, Nanjing gUiest) iityo dcto,Nnig209,Chi 210093, Nanjing Education, of Ministry University), ng and enpooe ogv iet o-eetn rre- or non-repeating ( to respectively have rise flares give bursts, clustered to peating or proposed single been with magnetar flaring dr eid”btenbrt n lsee us ar- burst clustered and bursts between periods” “dark measure rotation polarization, as properties such observational (RM; all FRBs nearly of explain cessfully a onepr ( counterpart cal ( 2019 2017 al. et Kumar al. et Platts 2019 Beloborodov 2017 al. et Metzger 2016 2015 Huang & Geng 2017 2016 2018 2016 eetn Rs ulmdlls a ee to refer can list model full A FRBs. repeating esl ta.2019 al. et Hessels iGoDai Zi-Gao 2 n neetn ru fmdl eeatt young a to relevant models of group interesting One http://frbtheorycat.org ,cru-us iprinmaue(M,adthe and (DM), measure dispersion circum-burst ), suc- to developed been has them of One ). ; ; ; ; ihlie l 2018 al. et Michilli aglt&Mtgr2018 Metzger & Margalit eoooo 2017 Beloborodov age l 2016 al. et Wang ( hn 2016 Zhang 2019 aaiso h nlrmat from remnants final the of galaxies t cdclas AC fwiedwarfs. white of (AIC) collapse uced cdfo eto trwiedwarf star–white neutron from uced ieFB ihntefaeokof framework the within FRBs like ( ,2 1, .B rlmnr aclto,we calculation, preliminary a By ). gr a nedcet magnetars, create indeed can rgers 2019 B101.05+5lk FRBs 180916.J0158+65-like RB Rs osbyfo magnetars from possibly FRBs) eN curn nultra-strong an acquiring NS he htejee l 2017 al. et Chatterjee roa Se,bnr neutron binary (SNe), ernovae ,acutn o o-eetn and/or non-repeating for accounting ), dcneto rpsil via possibly or convection nd ) ; 2 o oFs ai Bursts Radio Fast to ion ,pritn ai oreadopti- and source radio persistent ), . ; u&Kmr2018 Kumar & Lu hne sas bet well to able also is channel , ; eoooo 2017 Beloborodov r 92;Rdobursts Radio (992); ars 2017 ; , ,feunydwwr drift downward frequency ), a ta.2016 al. et Dai hn 2014 Zhang 2019 ; ybrk 2014 Lyubarsky ; ; ytkve l 2016 al. et Lyutikov eze ta.2017 al. et Metzger eze ta.2019 al. et Metzger ; edla tal. et Tendulkar ; ; eze tal. et Metzger ; uaee al. et Murase ; ege al. et Deng ue al. et Gu .cn na ; Katz ; ; ; 2 Zhong & Dai rival times appearing in FRB 121102 and its hosted low tron star–white dwarf (NSWD) mergers. This channel metallicity dwarf star forming galaxy and its surround- has also been briefly mentioned and/or discussed by ing dense, highly magnetized, and dynamic plasma envi- Liu (2018, 2020), Khokhriakova & Popov (2019), and ronment (Chatterjee et al. 2017; Tendulkar et al. 2017) Beloborodov (2019) previously. To answer this question, within the framework of synchrotron maser emission we need to study (1) whether this channel can form mag- from decelerating relativistic blast waves produced by netars or not, and (2) if it can, whether the formed mag- flare ejecta from young magnetars (Metzger et al. 2019). netars can account for the observations of FRB 180924- In this framework, repeating FRBs similar to FRB like bursts in the flaring magnetar framework. If this 121102 are arise from young and very active millisec- channel can indeed form magnetars, it could increase ond magnetars quite possibly connected with superlumi- the magnetar event rate to some extent and contribute nous supernovae (SLSNe) and long-duration gamma-ray to at least a subset of FRBs similar to FRBs 180924, bursts (LGRBs; Metzger et al. 2017). Therefore, young 190523, 181112, and even 180916.J0158+65. In the fol- magnetars giving rise to FRB 121102-like FRBs might lowing, we organize the structure of the paper: §2 intro- be formed during the core-collapse of massive stars as- duces the possibility and speculation that NSWD merg- sociated with SLSNe or LGRBs (SLSNe/LGRBs chan- ers could form magnetars; whether or not the NSWD nels). channel can explain the observations of FRB 180924-like On the other hand, young millisecond magne- cases is discussed in §3; and a summary and discussion tars could also be born in binary neutron star are presented in §4. (BNS) mergers (Rosswog et al. 2003; Price & Rosswog 2006; Giacomazzo & Perna 2013), binary white dwarf 2. MAGNETARS FROM NSWD MERGERS (BWD) mergers (King et al. 2001; Yoon et al. 2007; The explosive outcomes of NSWD mergers have Schwab et al. 2016), or accretion-induced collapse been explored in the literature (Metzger 2012; (AIC) of white dwarfs (WD; Nomoto & Kondo 1991; Margalit & Metzger 2016, 2017; Zenati et al. 2019a,b; Tauris et al. 2013; Schwab et al. 2015). These mag- Fern´andez et al. 2019), but the final remnants netars could produce FRBs analogous to FRB 180924 of these events have been little investigated (see (Bannister et al. 2019), as suggested by Margalit et al. Paschalidis et al. 2011a,b; Margalit & Metzger 2016). (2019). Compared with FRBs created from magnetars Generally, there are two evolutionary pathways for born in SLSNe/LGRBs channels, the FRBs produced NSWD binaries, which depend on the critical mass from magnetars born in BNS/BWD/AIC channels could ratio q = M /M , where M is the crit- have a distinct observational properties. Just like FRBs crit WD,crit NS WD,crit ical WD mass and M is the NS mass. The first 180924 and 190523 (Ravi et al. 2019), likely as well as NS pathway is that the WD fills its Roche lobe and its FRB 181112 (Prochaska et al. 2019), they host an old matter undergoes stable mass transfer to the NS if massive galaxy with a relatively low rate of star for- q < q , evolving into an ultra-compact X-ray bi- mation and relatively high metallicity, lie in a large crit nary. The second pathway is that the WD is tidally spatially offset location relative to the central contain- disrupted by the NS via unstable mass transfer on a ment region of the galaxy, and have low DM and RM rather short dynamical timescale for q > q , lead- contributions from the host galaxy and no bright per- crit ing to an NSWD merger (Hjellming & Webbink 1987; sistent radio source (Bannister et al. 2019; Ravi et al. Hurley et al. 2002). The critical mass ratio is related 2019). If this is the case, FRBs could be divided into two to the critical WD mass M , which is found to populations: FRB 121102-like bursts stem from young WD,crit be M = 0.37 M (van Haaften et al. 2012) or magnetars born in SLSNe/LGRBs channels, while FRB WD,crit ⊙ M = 0.2 M (Bobrick et al. 2017). Thus an 180924-like cases come from those young magnetars WD,crit ⊙ NSWD merger with q > q is in the case of un- born in BNS/BWD/AIC channels. Additionally, most crit stable mass transfer. Toonen et al. (2018) pointed FRB 180924-like bursts should also be repeating in out that over 99.9% of semi-detached NSWD systems the flaring magnetar framework due to the event rate would merge when M =0.2 M , which indicates comparison between magnetars and total FRB events WD,crit ⊙ that the NSWD merger is a prevalent fate of semi- (Nicholl et al. 2017; Margalit et al. 2019; Ravi 2019), detached NSWD binaries. After an NSWD merger, as which is supported by FRB 171019 followed by faint shown by Paschalidis et al. (2011a,b), the final rem- bursts (Kumar et al. 2019). nant both in the inspiraling case and in the head-on In this paper, we investigate whether or not FRB case is a spinning Thorne-Zytkow-like object (TZIO; 180924-like bursts are also likely to be generated by Thorne & Zytkow 1977) surrounded by a massive ex- magnetars born in an alternative possible channel: neu- tended hot disk composed of WD debris without ex- FRBs from Magnetars Formed from NSWD Mergers 3 plosive outcomes. Considering the disk winds and nu- ical timescale ∼ 3 × 104 s. In this case, the accretion clear burning, Margalit & Metzger (2016) suggested an of the final NS surrounded by a massive disk may let NSWD merger likely forms an isolated millisecond pul- it differentially and rapidly rotate, as possibly shown sar surrounded by an accretion disk at the final stage. by the simulation results of Paschalidis et al. (2011b), Whether or not these final remnants evolve into mag- even if Paschalidis et al. (2011b) did not take into ac- netars has received less attention. How the magnetic count explosive outcomes. Moreover, during this short fields of the final remnants evolve remains unknown. dynamical timescale, transient ohmic dissipation of the Fortunately, it has been suggested that the ultra-strong final NS could be ignored, see an estimate in Sun et al. magnetic fields in magnetars may result from two main (2019) and Equation (9) of Urpin & Konenkov (1997). scenarios (for a review see Turolla et al. 2015). Al- Therefore, we just need to focus on the magnetic field though these two scenarios are mainly used for nascent amplification of the final NS that accretes the WD de- NSs born in SLSNe/LGRBs/BNS/BWD/AIC channels, bris material from the disk during the merger. Owing we guess that they might also be used for “renascent” to the lack of previous studies of the magnetic field (magnetic field undergoes amplification) NSs formed in evolution of NSWD mergers with unstable and rapid the NSWD channel. mass transfer, we perform only a preliminary analysis on the α − ω dynamo induced by possible differential 2.1. α–ω Dynamo rotation and/or convection that can amplify the NS’ magnetic field during the NSWD mergers. In the fi- The first scenario that we consider is magnetic field nal paragraph of this subsection, we would also discuss amplification by a vigorous dynamo action at the the dynamo induced by the magneto-rotational instabil- early, highly convective stage of the NSs after merg- ity (MRI; Balbus & Hawley 1998) in the disk that also ers: the α dynamo arising from the coupling of con- might contribute to the NS’ magnetic field amplification. vective motions and rotation, and the ω dynamo driven We assume that the final NS with differential rotation by differential rotation (Duncan & Thompson 1992; induced by accretion has two components: the core and Thompson & Duncan 1993). For an NSWD merger, ∼ the shell divided by a boundary at the radius Rc = 0.5R the magnetic field of the NS remnant surrounded by a (where R is the NS radius), as done by Dai et al. (2006), massive disk may either increase via a dynamo winding- its accretion can be generally determined by the Alfv´en up process (as suggested by Paschalidis et al. 2011b) or radius decrease through an enhanced ohmic dissipation of ac- 3 4/7 −1/7 −2/7 creted matter in the NS’ crust (Urpin & Konenkov 1997; rm = (BsR ) (GM) M˙ Konar & Bhattacharya 1997; Cumming et al. 2004), 4/7 12/7 −1/7 Bs R M somewhat similar to the finding of Sun et al. (2019). =5.3 1012 G 12 km 1.4 M We discuss whether this finding is true in the follow- ⊙ −2/7 ing. Bobrick et al. (2017) showed that NSWD merg- M˙ × km, (1) ers exhibit an exponentially increasing rate of mass −6 −1 10 M⊙ s ! transfer during different phase transitions (see Figure
12 in Bobrick et al. 2017), and the NSWD mergers where Bs, R, M, and M˙ are the surface magnetic dipole in which the WDs have a higher mass would have a field strength, radius, mass, and accretion rate of the shorter dynamical timescale between the onset of sig- NS, respectively, and the corotation radius nificant mass transfer and the final merger (e.g., only 1/3 −3 −5 4 2 1/3 2 tdyn ∼ 10 (10 ) yr ∼ 3 × 10 (3 × 10 ) s for WD GM GMPs rc = = mass ∼ 0.75(1.2) M , see Figure 13 in Bobrick et al. Ω2 4π2 ⊙ s 2017). However, it is not true that the maximal rate 1/3 2/3 3 M Ps of disk accretion onto the final NS in Figures 11 and =7.8 × 10 km, (2) 1.4 M⊙ 10 s 12 of Bobrick et al. (2017) is limited by the Edding- ton rate since the disk accretion of NS can be highly where Ωs = 2π/Ps and Ps are the angular velocity and super-Eddington. If the majority of mass is lost via a spin period of the NS’ shell, respectively. One additional disk wind or possible a jet in the mass transfer pro- key radius is the light cylinder radius, cess and only 0.05 M can be accreted onto the NS, ⊙ P as discussed in Margalit & Metzger (2016) for a WD r = c/Ω =4.8 × 105 s km. (3) lc s 10 s with mass 0.6 M⊙ (close to 0.75 M⊙), the NS would accrete the disk material onto its surface with an aver- One expects that in the case of rm < rc < rlc for a −6 −1 12 age rate M˙ ∼ 10 M⊙ s during the short dynam- normal NS in an NSWD merger with Bs = 10 G, 4 Zhong & Dai
M = 1.4 M⊙, R = 12 km, and Ps = 10 s, as well where Ic is the moment of inertia of the core. Through −6 −1 as an accretion rate M˙ = 10 M⊙ s , disk mate- Equations (4), (8), and (9), we can solve Ωs, Ωc, ∆Ω, rial is column-accreted onto the NS’ surface and leads and Bφ as illustrated in Figure 1 via numerical cal- to the NS’ shell to spin up (Frank et al. 1992). Follow- culation, combining Equations (5) and (7). To obtain ing Piro & Ott (2011) and Dai et al. (2006), the time- these results, we have also employed: (1) typical values dependent angular velocity for the NS’ shell can be for a normal NS in an NSWD merger: M = 1.4 M⊙, 12 solved by R = 12 km, initial period Ps,0 = 10 s, Bs = 10 G, I =∼ I ∼ 1045 g cm2, ǫ = 0.3 (Dai et al. 2006), as dΩs s c I = N − N − N , (4) ˙ −6 −1 s dt acc dip mag well as an accretion rate M = 10 M⊙ s , such that the term N can be ignored in comparison with N where I is the moment of inertia of the shell, and (1) dip acc s even if P possibly reaches down to its break-up limit N is the accretion torque described by, when r < R s acc m P = 0.96 ms (Lattimer & Prakash 2004); (2) ini- for a normal NS from Equation (1), s,min tial conditions: Ωs,0 = 2π/Ps,0, Ωc,0 = (1+ A0)Ωs,0 Ωs (the initial angular velocity of the core Ω should be N = 1 − (GMR)1/2M,˙ (5) c,0 acc Ω larger than that of the shell Ω for a normal NS), K s,0 A = 10−3 related to a small residual differential ro- 3 1/2 0 where ΩK = GM/R ; (2) Ndip is the magnetic 8 tation, and Bφ,0 ∼ 10 G; (3) boundary conditions: dipole radiation torque for accreting NSs with rm < rlc, Ωs ≤ Ωs,max = 6541 due to Ps,min = 0.96 ms, Ωc < enhanced over the standard dipole torque by a factor Ω = c/R , and B < B = 1017 G because of 2 c,max c φ φ,max of (rlc/rm) > 1 due to the enhanced open magnetic the buoyancy effect. From Figure 1, we can acquire the field lines via the compression of the magnetosphere following: (Parfrey et al. 2016; Metzger et al. 2018), read as
2 6 3 2 • The top panel shows that the angular velocity 2B R Ω rlc N = s s of the shell gradually increases up to its limit at dip 3c3 r m about 200 s and lasts this lasts until the end of 2Ωs 3 6/7 2/7 4/7 = (BsR ) (GM) M˙ ; (6) the dynamical process. While the angular velocity 3c of the core reverses (rotating in an opposite direc- (3) Nmag is the magnetic torque acting between the core tion) at about 50 s and its absolute value rapidly 4 and shell, written as rises to a very large value c/Rc at around 1000 2 s. This is because the toroidal magnetic field Bφ N = R3B B , (7) mag 3 c r φ reverses, as in the bottom panel. The evolution of the differential angular velocity ∆Ω follows the where Br = Bs/ǫ (here, ǫ is defined by the ratio of the effective surface dipole field strength to the radial field angular velocity of the core. strength) and Bφ are the radial magnetic field compo- • nent and toroidal field component, respectively. They In the bottom panel, the toroidal magnetic field can be related to each other by Bφ rapidly declines to zero and then reverses be- fore 0.1 s, its absolute value continues going up to dB 17 φ = (∆Ω)B ≡ (Ω − Ω )B , (8) the buoyancy limit 10 G at about 150 s. dt r c s r where ∆Ω is the differential angular velocity and Ωc is In short, these results manifest the toroidal magnetic the NS’ core angular velocity. On the other hand, the field can be enhanced during the dynamical timescale time-dependent angular velocity for the NS’ core can be as long as the initial remnant NS in an NSWD merger solved by3 has a small residual differential rotation. Additionally, during the field amplification, the spin-down torque Ndip dΩc Ic = Nmag, (9) responds to the magnetic dipole radiation luminosity dt
Ldip = NdipΩs 3 Please note that the right term of Equation (9) has no negative 2 sign, differing from Equation (3) in Dai et al. (2006) because the 2Ω = s (B R3)6/7(GM)2/7M˙ 4/7, (10) magnetic field is amplified via the differential angular velocity 3c s resulting from the angular momentum transport of the accreting material onto the NS’ shell rather than the angular momentum 4 loss of the core for an accreting NS system we consider here. The The core angular velocity Ωc that exceeds the break-up limit of same reason that the right term of Equation (4) in this paper the shell is reasonable since this break-up limit should not be that differs from Equation (2) in Dai et al. (2006). of the core angular velocity. Instead, its limit should be c/Rc. FRBs from Magnetars Formed from NSWD Mergers 5
1 1 1 1 1 1 Δ
1 1 1 1 1 1 1 1 1 1