Arxiv:1805.05262V1

Home , Bright giant

arXiv:1805.05262v1 [astro-ph.HE] 14 May 2018 lysgetdt eascae ihavreyo astrophysi long of as variety such a with phenomena, associated be gene to are suggested magnetars) millisecond ally as (dubbed (NSs) stars tron eso h el onmgeas uha h ioemagnetic dipole the parame- as important such magnetars, fields some born and newly LCs the masses of the ters ejected fitting the Through SLSNe, [3–5]. ejecta brilliant of (SN) quite supernova SN the con- the heat and making to extracted energy be heating can into magnetar verted born newly energ central rotational the huge 1 of model, 5, magnetar-powered [4, the Refs. In e.g., (see, 14]). repro scenario well magnetar millisecond be the can within by (LCs) duced powered curves light be their SLSNe to because type-I magnetars suggested as generally SLSNe (classified are of [10] [11]) subclass a SLSNe as hydrogen-poor Especially, the [9]. bursts [ radio transients fast luminous pre- and and x-ray evolving an rapidly by some accompanied [7], cursor possibly are that mergernovae [6] bright emissions [3–5], (SLSNe) supernovae perluminous ae,temgeasaerqie opossess to required 5 are magnetars the cases, therein). × rbn h hsc fnwybr antr hog observa through magnetars born newly of physics the Probing el onmlieodrttn ihymgeie neu- magnetized highly rotating millisecond born Newly 10 14 B d ,and G, nta pnperiods spin initial , antr novd eso ht eedn nteinternal the on depending that, show curves We light involved. SLSNe magnetars magnetar-powered the calculate mag We through mainly sion. down spin which magnetars, rotating fast T quantity nlsgo to grow angles o el onmgeaswt ioemgei fields magnetic dipole with magnetars born newly For r obmsi h C if LCs the in bumps no are upi h Cdsper o h antrwt smaller with magnetar the for disappears LC the in bump id fbmsaeosre,b tigteeLswt u mode our with LCs these fitting by observed, are bumps of kinds ihsc id fbmshtet a led u oe( some put already may hitherto bumps of kinds such with ASnumbers: PACS magnetars. T u lotecuilpyia quantities physical crucial the also but P B ¯ c c i t P n h radii the and , hrfr,osraino LN C a rvd e appro new a provide may LCs SLSNe of observation Therefore, . h eta nie fsm uelmnu uenve(SLSN supernovae superluminous some of engines central The illa obt rgtrpa n rgtrbm naL.Wh LC. a in bump brighter a and peak brighter a both to lead will for i ∼ 3 .INTRODUCTION I. – s(e es 1,1]adreferences and 15] [14, Refs. (see ms 1–8 P R 2 DU eto uefliiy up ol peri h LN C af LCs SLSNe the in appear could bumps superfluidity, neutron unCheng Quan / a a can hr am-a uss[,2,su- 2], [1, bursts gamma-ray short π/ R P .Tevleof value The 2. DU ff i 4 c h rsn ieadtelmnst fabm,a ela th as well as bump, a of luminosity the and time arising the ect nvriyo hns cdm fSine,Biig100049, Beijing Sciences, of Academy Chinese of University a edtrie.I most In determined. be can fdrc ra(U oe fnwybr antr,a ela t as well as magnetars, born newly of cores (DU) Urca direct of 5 colo hsc n Mechanical and Physics of School 1 eta hn omlUiest,Whn407,China 430079, Wuhan University, Normal China Central , 2 ue nvriyo dcto,Whn400,China 430205, Wuhan Education, of University Hubei hns cdm fSine,Biig104,China 100049, Beijing Sciences, of Academy Chinese ∗ T hagNnZhang Shuang-Nan , c = 2 ua nvriy ua 302 China 430072, Wuhan University, Wuhan 2 ∗ e aoaoyo atceAstrophysics, Particle of Laboratory Key × lcrncades [email protected] address: Electronic T B 1 c d 10 colo hsc n Technology, and Physics of School nttt fHg nryPhysics, Energy High of Institute eemnsteaiigtm n h eaieapiueo b a of amplitude relative the and time arising the determines ∼ B 9 ¯ t ,or K, , 5 3 nttt fAstrophysics, of Institute R × DU Dtd o 2017) Nov (Dated: 10 and , R DU 13 cal 8], 2– r- – y - 2 , , = , T 4 B u-e Yu Yun-Wei , c 1 d . ol edtrie.Nnbevto fSSeLCs SLSNe of Nonobservation determined. be could 5 = adaMIdie ubln yao fe oebuc fthe if bounce core of after field dynamo) turbulent toroidal MRI-driven with a (and PNS a that ∼ showed CCSNe of lations ierwnig ttoayaceinsokisaiiy M instability, shock an compression, accretion and flux stationary magnetic as winding, amplified such be linear ways, can of PNSs series of a fields through magnetic the [24 SLSNe here, with NS interest associated rotating of (CCSNe) fast the supernovae a and collapse in core [23], In act stars massive instabilities Tayler of and collapse r-mode core the of and merger [22] the NSs during arises mag (MRI) or instability [21], torotational instability Kelvin-Helmholtz [20], (PNS) star di a an [19], in progenitor dynamo magnetized of highly fields a of magnetic collapse during core strong conservation flux of magnetic instance, formation for magnetars; the explain to posed 0[6 l niaeta h oodlfil ol eafwto few a be could field ing toroidal the that ∼ indicate all [16] 20 h lwpaemdlto ntexryeiso fMagnetar of emission x-ray 0142 the 4U in [17], Observationall modulation afterglows phase burst slow gamma-ray 17]). the short [16, some toroid of LCs Refs. interior x-ray the e.g., stronger (see, even fields possess magnetic to considered erally P × huhvr rough very though i 5 esgetta,oc h LN C ihsuch with LCs SLSNe the once that, suggest We . 0 ie ihrta h ioefil famgea,reach- magnetar, a of field dipole the than higher times 100 eie h togsraedpl ed,mgeasaegen- are magnetars fields, dipole surface strong the Besides 10 10 & ei ioerdainadgaiainlwv emis- wave gravitational and radiation dipole netic oodlmgei fields magnetic toroidal × Ls ihtetl nl vlto fnwyborn newly of evolution angle tilt the with (LCs) ∼ 5 15 lcrclEngineering, Electrical 10 ,ntonly not l, m oevr ti neetn htastronger a that interesting is it Moreover, cm. 10 14 )aegnrlysgetdt enwyborn newly be to suggested generally are e) o even (or G α + 16 l epn te uniisucagd the unchanged, quantities other keeping ile G, − 1[8,adtebih in aefo G 1806- SGR from flare giant bright the and [18], 61 3 c opoetepyiso el born newly of physics the probe to ach rhge.Vrosmcaim aebe pro- been have mechanisms Various higher. or G n ioPn Zheng Xiao-Ping and , B ω ¯ t e h aiu ihswe h tilt the when lights maximum the ter = yao[5.Seicly ueia simu- numerical Specifically, [25]. dynamo ff ino uelmnu supernovae superluminous of tion B 4 rnilyrttn ilscn protoneutron millisecond rotating erentially d osrit on constraints ) . 6 and ∼ × 10 10 P eklmnst faLC. a of luminosity peak e China i 16 16 B ¯ fnwybr magnetars born newly of t ,and G, h nta pnperiods spin initial the , 2] a epoue i MRI via produced be can [23]) G eciia temperature critical he B P ¯ t i , 3 = P , 5 i s there ms, 1 , m.The ump. R DU and , α two − ne- RI, the ω al y, ]. 2

2 precollapse iron core is highly magnetized and rapidly rotat- perature Tc adopted and the cooling mechanism of the newly ing [25–28]. Meanwhile, the PNSs may probably have ini- born magnetar. Consequently, the coupled evolutions of the tial spin periods of the order of milliseconds [28–32], which tilt angle, spin, and stellar temperature of a magnetar should could further trigger the turbulent dynamo that is driven by be taken into account while calculating the magnetar-powered differential rotation and convection and amplify the interior SLSNe LCs. Furthermore, the effects of physical quantities B¯ t 16 toroidal fields to ∼ 10 G [20]. Hence, newly born rapidly and Tc should also be involved. 16 rotating magnetars with toroidal fields of a few×10 G can 3 In fact, the critical temperature for P2 neutron superflu- possibly be formed in CCSNe with highly magnetized and idity is a function of density, Tcn(ρ), which has a parabolic fast rotating precollapse cores. It also seems that the strong shape and peaks at a certain density ρ between the core-crust magnetic fields of newly born magnetars are tightly related to boundary and the stellar center. The maximum value of Tcn(ρ) their fast rotations. However, an estimation of what portionof 3 is usually dubbed as the critical temperature Tc for P2 neu- newborn NSs are magnetars with such strong toroidal fields tron superfluid transition. The specific value of Tc is still un- is impossible at present because we still know little about the certain since the roles of medium effects and complicated in- properties of the progenitor cores and the evolution process teractions on the result of Tc are not clearly known [42, 43]. from PNSs to newborn NSs. Strong magnetic fields can in- Previous result suggests that Tc may be within the range from duce nonaxisymmetric quadrupole deformation, which man- a rather small value to ∼ 1010 K (see Ref. [44] and refer- ifests the newly born magnetars as strong gravitational wave ences therein). However, observations of the cooling behavior (GW) sources [16, 33–37]. However, initially, the tilt angle of the NS in Cassiopeia A [45, 46] and surface temperatures between the spin and magnetic axes of a magnetar may be 8−9 of other isolated NSs [47] suggest Tc ∼ 10 K. Following very tiny,1 the gravitational wave emission would be strongly Refs. [34, 45], in this paper, Tc is left to be a parameter within suppressed consequently [34]. 8 9 the range 5 × 10 –2 × 10 K. We will show the value of Tc The tilt angle evolution of a newly born magnetar with can obviously affect the shape of a SLSN LC, which suggests strong toroidal field was first investigated in Ref. [34] and that observation of SLSNe may provide another approach to then involved in the calculation of gravitational wave back- determine Tc. ground from newly born magnetars [38]. The tilt angle trends The effect of magnetically induced GW emission on the to increase to π/2 in order to minimize the NS’s spin energy. magnetar-powered SLSNe has been studied in Refs. [15, 41, Generally, the angle evolution can be divided into two stages 48]. The results show that strong toroidal fields can over- [34]. The first stage is when the stellar temperature is so high 9 all reduce the emitted luminosities of SLSNe because a large (& 10 K) that the whole star is in the liquid state without amount of rotational energy of the central magnetars is re- superfluidity. The tilt angle evolution is determined by the leased in the GW channel [41, 48]. To require that the ma- competition between the bulk viscosity (BV) of stellar matter jority of rotational energy can be used to power the SLSNe, and gravitational radiation reaction (GRR) of the magnetar. the GW emission must be weakened, and an upper limit for As the newly born magnetar cools down due to intense neu- 16 toroidal fields is derived as B¯ t . several × 10 G [15]. How- trino emission, a solid crust can form for stellar temperatures 9 3 ever, all these results are obtained under the assumption of lower than ∼ 10 K [39]. Moreover, the P2 neutron superflu- constant tilt angles as mentioned before. A more detailed idity will occur in the core when the stellar temperature drops investigation about the role of B¯ t on the SLSNe LCs is still below a critical value Tc. After the formations of a solid crust needed as the angles should indeed evolve with time. and neutron superfluidity, the tilt angle evolution goes into the second stage, in which the angle evolution is driven by viscous On the other hand, the cooling mechanism of newly born dissipation of the free-body precession due to core-crust cou- NSs is still an open issue. It is hard to address this issue both pling [40]. Evolution of the tilt angle can lead to a change in theoretically and observationally because of the poor knowl- the magnetic dipole luminosity that is ejected into a SN; thus, edge of dense matter properties and obscuration of thermal the SLSN LC may be changed accordingly. This represents emissions of NSs by the surrounding dense ejected materi- the fundamental starting point of this paper, and the essential als. As generally considered, in the early period, the classi- 3 difference compared to previous work, in which a constant an- cal NSs composed purely of neutrons, protons, and electrons gle is generally assumed (e.g., Refs. [3–5, 14, 15, 41]). cool down mainly through the modified Urca (MU) process [49] if the proton fraction in stellar interiors is below a thresh- The most crucial quantities that determine the angle evo- old of about 11% [50, 51]. However, some dense matter equa- ¯ lution in the first stage are toroidal field Bt (volume-averaged tions of state (EOSs) (e.g., Akmal-Pandharipande-Ravenhall strength), spin period P, and stellar temperature T. In the later (APR) [52] and Prakash-Ainsworth-Lattimer (PAL) [53]) pre- stage, since the angle could increase to π/2 in a very short dict that in the central region of some NSs with sufficient time if the fast rotating magnetar has a strong toroidal field masses [54, 55] the proton fraction can surpass the thresh- [16, 40], the start time of this stage becomes crucial. The start time of the second stage strongly depends on the critical tem-

2 The crust formation temperature [39] may be a little higher than (or ap- proximate to) Tc; we therefore expect that the second evolution stage will 1 This seems to be a direct consequence of the field amplification due to MRI begin when Tc is reached. and the α − ω dynamo. 3 In this paper, we focus on this type of NSs only as in Ref. [34]. 3 old, leading to the occurrence of the direct Urca (DU) process born magnetars. For the toroidal-dominated interior mag- [50, 51]. The DU neutrino emission in the central region of netic field configuration,4 the quadrupole ellipticity of mag- ¯ 2 4 2 NSs can greatly expedite the cooling of NSs. The size of the netic deformation is ǫB = −5Bt R /(6GM ) [33, 34]. Follow- DU core dependson the mass of a NS and EOS (see Ref. [56] ing Dall’Osso et al. [34], at the first stage, the tilt angle χ of a for a review), both of which are uncertain for NSs embedded liquid NS evolves under the combined effect of GRR and BV, in SNe. Since the stellar temperature evolution could affect which can be written as the tilt angle evolution, and further the shape of SLSNe LCs, cos χ 2G observation of SLSNe may give some clues on two crucial is- χ˙ = − Iǫ2 Ω4 sin χ cos χ(15 sin2 χ + 1). (2) 5 B sues: (i) could the DU process occur in a NS, and (ii) if it τd sin χ 5c occurred, how large is the DU core? These could help one to The first term on the rhs of Eq. (2) is related to the effect of understand the NS interior structures and constrain the EOS BV of stellar matter on damping of the free-body precession, of dense matter. which can essentially increase χ of a prolate star (ǫB < 0). The The paper is organized as follows. We show the evolution correspondingdamping timescale of free-body precession cal- of newly born magnetars in Sec. II. The model for magnetar- culated for a classical NS with only the MU process involved powered SLSNe is briefly introduced in Sec. III. Our results is [34] are presented in Sec. IV. Finally, a conclusion and discussions 2 2 2 −6 are given in Sec. V. cot χ B¯ t P T τd ≃ 3.9 s 1 + 3cos2χ 1016 G! 1 ms 1010 K −1 II. EVOLUTION OF NEWLY BORN MAGNETARS M R 3 × . (3) 1.4M⊙ ! 12 km The collapse of a massive progenitor core may give rise to a The second term on the rhs of Eq. (Ref. 2) (taken from [63]) differential rotating PNS with strong convective motions in its represents the damping of the free-body precession due to interior. Initially, the ultrahot (with a temperature of a few tens GRR, which can actually lead to alignment of the spin and MeV) PNS is opaque to neutrinos and may have a radius of magnetic axes even for a NS with ǫ < 0. several tens of kilometers. Subsequently, as the PNS becomes B It should be stressed that in the presence of the DU process transparent to neutrinos, it will contract and become a newly the BV of stellar matter is several orders of magnitude higher born NS with a radius of ∼ 10 km at ∼ 10 s after the core than that with MU process involved [64]. Hence, if the DU bounce [29, 31]. Contraction of the PNS can lead to spin-up process could occur in a NS core, the damping timescale τ of the newly born NS because of angular momentum conser- d should be modified accordingly in principle. However, in the vation; thus, the newly born NS possibly has a spin period of simple phenomenological NS structure models we consider & 1 ms at birth [29, 57]. Simulations of the collapse of mas- hereinafter, the DU process, if it occurred in the core region, sive progenitor cores have demonstrated that the initial spins the radius of the DU core, R , is assumed to be much smaller of newly born NSs could indeed reach of the order of 1 ms DU than R. This is reasonable because, depending on EOSs, the [30, 58]. It should be noted that fallback accretion probably DU process can marginally occur or be quenched completely exists at early periods, which may further spin up the central in a 1.4M NS (see, e.g., Refs. [54–56, 65]). The DU core remnant and result in a newly born magnetar with an initial ⊙ contains a rather small portion of the total precession energy period very close to 1 ms. Theoretically, fast spin (. 5 ms) is due to the small R and (thus) the small DU core mass M . suggested to be an indispensable condition so as to avoid the DU DU As a very rough estimation, taking the NS model derived in early unstable phase that newly born magnetars will undergo Ref. [56] for a 1.4M NS as an example, the precession en- [59]. In this paper, we take 1 ms as the possibly minimum ⊙ ergy of the DU core is [34] initial spin periods of newly born magnetars. We also note that newly born magnetars may have different initial spin pe- 1 E = − I Ω2ǫ cos2 χ riods due to different progenitor properties, and more detailed pre,DU 2 DU B simulations are still needed in order to know how rapidly such M R 2 Ω 2 magnetars can rotate. ≃ × 46 DU DU 5.3 10 erg 4 A newly born magnetar is generally considered to be spun 0.023M⊙ ! 2.4 km 10 rad/s ! down via magnetic dipole radiation (MDR) and magnetically |ǫ | × B cos2 χ, (4) induced GW emission, which can be expressed as follows 10−3 ! (e.g., Ref. [38]), 2 where IDU = 2MDUR /5 is the moment of inertia of the DU 2 6 3 2 5 DU BdR Ω 2GǫBIΩ core. While the precession energy of the remaining part of Ω=˙ − sin2χ − sin2 χ(15 sin2 χ + 1), (1) 6Ic3 5c5 where Ω is the angular frequency, Bd is the surface dipole magnetic field at the magnetic pole, and I = 0.35MR2 is the 4 Though some magnetohydrodynamics simulations show that a twisted- moment of inertia of the NS with M and R representing the torus magnetic configuration composed of both poliodal and toroidal fields stellar mass and radius, respectively [60]. From now on, we may naturally form in NS interiors [61], the dominant one is still the take typical values M = 1.4M⊙ and R = 12 km for newly toroidal component [62]. 4

49 2 2 the star is Epre,MU ≃ 8.1 × 10 erg Ω4 ǫB,−3 cos χ ≫ Epre,DU, DU and MU processes of the stellar matter, and on the other where we have adopted the notation Q = Q/10x here and hand provide new channels for neutrino emission that are as- x hereinafter. Thus, it can be seen that, though the damping rate sociated with the pair breaking and formation processes [51]. of the precession energy E˙pre is much larger in the DU core For simplicity, below Tc, the thermal evolution of magnetars due to its larger BV ([34, 64]), dissipation of most of the pre- is ignored since it is trivial to our results. cession energy still occurs in the MU region. Consequently, the occurrence of the DU process in a small core region will not modify the form of τd given in Eq. (3) significantly. III. MODEL FOR MAGNETAR-POWERED SLSNE If the orthogonal configuration between the two axes is not reached in the first stage, when the stellar temperature cools Newly born magnetars lose rotational energy mainly via down to Tc, the second evolution stage will begin, in which MDR and GW emission, of which only theenergyin the MDR the viscosity due to core-crust coupling plays a dominant role channel can be used to thermalize the SNe ejecta. The MDR in damping of the free-body precession [40]. The tilt angle of luminosity emitted by a newly born magnetar is a prolate star could increase to π/2 on a timescale τ ≃ nP/ |ǫB| with n ∼ 102–104 representing the number of precession cy- B2R6Ω4 L = d sin2χ. (6) cles [16, 33, 40, 66]. For a magnetar with P ∼ 30 ms and m 6c3 16 B¯ t = 4.6 × 10 G as shown in Fig. 2, at the beginning of the second stage, the maximal orthogonal timescale is esti- During the expansion of SN ejecta, it loses internal energy 16.66 −2 due to adiabatic expansion and thermal radiation of the ejecta; mated to be τmax ≃ 0.49d(P/30 ms)(B¯ t/10 G) by tak- 4 meanwhile, the ejecta can also be heated by the energies that ing n = 10 , which is still much shorter than the evolution 56 timescale (of the order of 100 d) of the magnetar. Hence, it stem from MDR, Ni cascade decay, and ejecta-circumstellar is reasonable to assume that χ = π/2 can be realized immedi- medium interaction. Since in the magnetar-powered scenario the dominant energy source is MDR from magnetars, the evo- ately when Tc is reached (see, e.g., Ref. [38]). Since the tilt angle evolution of a newly born magnetar is lution formula of the internal energy Eint can thus be written tightly related to the stellar temperature T, the thermal evolu- as [4, 8, 41, 68] tion is thus an important issue to be addressed. For the poor dE dV int = −P − L + L , (7) knowledge of dense matter EOS and the NS interior struc- dt ej dt th m tures, here we adopt a phenomenological NS model, in which = the NSs are comprisedof a small DU core with radius RDU and where Pej Eint/(3V) is the pressure dominated by radiation a large MU shell of radius R − RDU. Moreover, for simplicity, with V denoting the volume of the SN ejecta and Lth is the the isothermal assumption is adopted with temperature dis- thermal radiation luminosity. At early times, the ejecta are op- 2 tributing uniformly in a NS. The “DU core+MU shell” model tically thick, namely, the optical depth τo = 3κMej/(4πRej) ≫ was used previously in Refs. [54, 65] while discussing the 1 with κ, Mej, and Rej representing the opacity, ejecta mass, cooling of the 2.21M⊙ NS with all calculations based on the and radius of the SN ejecta, respectively. In this case, Lth has realistic EOS APR. For a 1.4M⊙ newly born magnetar, if the the following form [8]: DU process could take place in the core region, the evolution Eintc of T would follow the formula below, Lth = . (8) τoRej dT CV = −Lν,DU − Lν,MU While at late times the ejecta will become optically thin (τ ∼ dt o 1), one then has [8] 4π 3 4π 3 3 = − QDURDU − QMU R − RDU , (5) 3 3 Eintc Lth = . (9) 39 −1 Rej where CV ≈ 10 T9 erg K is the total heat capacity of 5 27 6 −3 −1 the NS and QDU ≈ 10 T9 erg cm s and QMU ≈ Moreover, the dynamical evolution of the SN ejecta can gen- 21 8 −3 −1 10 T9 erg cm s are the DU and MU neutrino emissivities, erally be determined by the following equations [4, 8], respectively [51, 56]. From Eq. (5), one can see, by setting R = 0, the thermal evolutionreturns to the pure MU cooling dRej DU = vej, (10) case, and the evolution equation is consistent with the analyt- dt − − ical expression T(t) = 109 K(t/τ + 10 6) 1/6 with τ ≃ 1 dvej c c = 4πR2 P /M , (11) yr [34, 51, 67]. It should be stressed that below the critical dt ej ej ej temperature T the presence of neutron pairing can on one c where v is the expansion velocity of the SN ejecta. hand suppress the heat capacity, neutrino emissivities of the ej

IV. RESULTS

5 For diﬀerent values of RDU, the expression for CV remains unchanged because the heat capacities of the stellar matter are the same for the DU and Combining Eqs. (1), (2), (5), (7), (10), and (11) and taking MU processes [44, 51, 65]. into account the tilt angle evolution in the second stage, we 5

can determine the LCs (evolution of Lth) of magnetar-powered 44 10

SLSNe. In all calculations below, typical values Mej = 5M⊙ T =2.0E9 K and κ = 0.2 cm2 g−1 are taken for the ejecta mass and opacity, c T =1.5E9 K

c 43 respectively [4, 41]. Moreover, the initial values for the SNe 10

T =1.2E9 K parameters, i.e., radius, velocity, and internal energy, are taken c

8 9 −1 51 T =1.0E9 K as Rej,i = 3 × 10 cm, vej,i = 10 cm s , and Eint,i = 10 c ◦ 42 = 10 = 10 erg, respectively. We take Ti 10 K and χi 1 for the (erg / s )

initial temperatures and tilt angles of newly born magnetars, th ¯ L respectively. The strength of Bt and the values of RDU and Tc 41 are set as free parameters in order to investigate their eﬀects 10 on the LCs. However, as discussed in Sec. I, B¯ t are taken to be of the order of 1016 G, and no more than 100B as inferred 40 d 10 from observations of magnetars.

Assuming diﬀerent critical temperatures Tc, the evolutions 80 of thermal radiation luminosities Lth (upper panel) of SLSNe χ and tilt angles (lower panel) of the central magnetars are 60 14 shown in Fig. 1. We take dipole ﬁelds Bd = 5 × 10 G,

16 toroidal ﬁelds B¯ t = 4.6×10 G, and initial spin periods Pi = 1 40 ms for the magnetars, which are assumed to cool down via (degree) only the MU process (RDU = 0) here. Since the growth of χ

20 of a newly born magnetar can be suppressed by its strong B¯ t during the ﬁrst evolution stage [see Eqs. (2) and (3)], χ will

0 increase to π/2 only when the magnetar cools down to Tc so that the viscosity due to core-crust coupling becomes eﬀec- 0 100 200 300 400

t (d) tive. The higher Tc is, the earlier χ = π/2 can be achieved, 9 as shown in the lower panel of Fig. 1. For Tc ≤ 1.5 × 10 K, the rapid growth of χ in the second evolution stage leads FIG. 1: Evolutions of the radiated luminosity Lth and tilt angle χ. The curves are obtained by assuming different critical temperatures to an enhancement in the injected power Lm, resulting in a T for 3P neutron superfluidity as shown in the legends. bump in the LCs after the maxima. With the decrease of Tc, c 2 the bump arises at later times; however, its relative amplitude gradually grows. This can be understood as follows. The GW emission is more sensitive to χ; thus, more rotational energy of a magnetar will be released in the GW channel if χ = π/2 is achieved earlier. The energy in the MDR channel that can (black curves) and a constant angle χ = π/2 [thick light-gray be used to energize the SN is therefore reduced, leading to a curves in panels (a), (b), and (c)]. The magnetar parameters relatively small bump. In contrast, for T = 2 × 109 K, no c B , B¯ , and P are taken the same as in Fig. 1, while the critical apparent bumps appear in the resultant LC. The reason is the d t i temperature is fixed to be T = 109 K. orthogonal configuration (χ = π/2) can be achieved at the c very beginning of the evolution if T is high enough. In this c The occurrence of the DU process in the core greatly ac- case, the rapid increase in L cannot lead to a remarkable in- m celerates the cooling of a magnetar [panel (d) of Fig. 2]; the crease in L , the evolution at the very early period of which th orthogonal configuration can thus be realized earlier in com- is essentially determined by the initial internal energy E of int,i parison with the result of only MU cooling [see panel (b)]. a SN because very little energy of the central magnetar has For R ≤ 1.5 × 104 cm, with the increase of R , the arising been injected into the SN. The very short timescale needed DU DU time of the bumps in the LCs is gradually brought forward. for χ = π/2 to be realized makes the resultant LC have little However, the luminosities of the bumps as well as the peak differences with that calculated by assuming a constant angle luminosities are gradually enhanced. The reason is a larger χ = π/2 (see Fig. 2). RDU will result in a lower stellar temperature T, which can The results above show that with the tilt angle evolution hinder the growth of χ of a magnetar in the first evolution involved the arising of a bump after the maximum light in a stage. A small χ can, on one hand, reduce the rotational en- 3 SLSN LC can be related to the occurrence of P2 neutron su- ergy loss due to GW emission of the central magnetar and thus perfluidity in the magnetar core. Given a specific Tc, while increase the energy in the MDR channel. On the other hand, other quantities remain unchanged, the arising time of the in the case of MDR-dominated spin-down,6 the peak luminos- bumps in LCs is determined by the thermal evolution of cen- ity of a magnetar-powered SLSN can be approximately given tral magnetars. This can be seen from Fig. 2, in which we show the results calculated by assuming that the magnetars cool down via both the DU and MU processes; however, the 6 radii RDU of the DU cores are different. For a comparison, we For the values of Bd, B¯ t, and Pi taken in this paper, the braking effect of also show the results obtained by assuming only MU cooling MDR overwhelms that of GW emission. 6

100 10

(erg / s)

40 (degree) th L

(a) (b)

10 1000

10 R =0.0 cm

R =1.0E4 cm

R =1.5E4 cm

R =1.5E5 cm

9 (Hz) 100

T (K) 10

(d) (c)

10 10

0 100 200 300 400 0 100 200 300 400

t (d) t (d)

FIG. 2: Evolutions of the radiated luminosity Lth, tilt angle χ, spin frequency ν, and stellar temperature T, calculated by assuming different DU core radii RDU for the central magnetars as shown in the legends. The thick light-gray lines show the results derived for a constant angle χ = π/2. The dotted line in panel (d) represents the critical temperature Tc. as [4, 48] different from the one obtained with a constant angle χ = π/2 (see, e.g., Refs. [41, 48]), and (ii) the thermal evolution ofa E τ L ∼ rot,i sd , (12) newly born magnetar, which is related to its interior structure, th,p 2 ff tdiff can obviously a ect the shape of the SLSN LC that is powered by it. where Erot,i, τsd, and tdiff are the initial rotational energy, spin- In Fig. 3, we show the evolution curves derived by setting 8 down timescale, and the photon diffusion timescale, respec- Tc = 5 × 10 K, as inferred from the cooling behavior of the tively. The spin-down timescale of a magnetar with tilt angle χ NS in Cassiopeia A [45]. One can see that for Tc of this value 3 2 2 2 6 2 ¯ is τsd = 3Ic Pi /(2π BdR sin χ). The peak luminosity is thus bumps still exist in the LCs for the values of Bt and RDU taken 2 3 2 6 2 2 B P Lth,p ∼ 6I c /(BdR tdiffsin χ), which shows that a smaller χ (see the legends), while d and i are kept the same as in Fig. will result in a higher Lth,p, as seen in Fig. 2. Furthermore, for 1. Again, the effect of the NS structure (the DU core radius 4 RDU ≤ 1.5 × 10 cm, at the points when the tilt angles grow to RDU) on the shape of a SLSN LC is clearly shown. Of the π/2, the magnetars have been considerably spun down, with most important, in contrast to previous results [15, 41, 48],a spin frequencies ν(=Ω/2π)in the range ∼ 24–50 Hz [see panel higher B¯ t in a magnetar does not necessarily lead to a lower (c)]. Lth,p of the resultant LC. Instead, with the tilt angle evolution ¯ If RDU is large enough that the DU cooling becomes dom- involved, a higher Bt can more obviously suppress the growth 5 inant, as is the case for RDU = 1.5 × 10 cm, the magnetar of χ [see panel (b) of Fig. 3]. Following the analysis above, cools down to Tc and has its tilt angle enlarged to π/2 almost the smaller χ will naturally lead to a higher peak luminosity immediately after its birth. The increase of χ doesnotleadto a and a brighter bump, as seen in panel (a). The close relation 9 bump in the LC, which is similar to the case of Tc = 2×10 K, between the strength of the toroidal field in a magnetar and 5 as shown in Fig. 1. The resultant LC for RDU = 1.5 × 10 cm the shape of the magnetar-powered SLSN LC may enable us differs little from that obtained with a constant angle χ = π/2, to probe the internal toroidal fields of newly born magnetars except that the former has a higher Lth,p, which results from through observation of SLSNe LCs. the suppression of GW emission and the promotion of Lth,p With χ = π/2 assumed, by fitting the SLSNe LCs in the before the orthogonal configuration is realized, as analyzed magnetar engine scenario, some central magnetars were found 14 above. In short, from Fig. 2, two important conclusionscan be to have Bd . 10 G and Pi & 1 ms [14, 15]. Moreover, as drawn: (i) with the tilt angle evolution involved, under some suggested in Refs. [29, 30], newly born magnetars may have conditions, the presence of strong B¯ t in a newly born magne- Pi & 1 ms. We therefore show in Fig. 4 the resultant LCs tar will result in a magnetar-powered SLSN LC that is quite powered by magnetars with lower dipole magnetic fields of 7

100

(erg / s)

40 (degree) rad L B =3.5E16 G, R =4.9E4 cm

t DU

B =3.4E16 G, R =4.9E4 cm 20

t DU

B =3.5E16 G, R =4.8E4 cm (a) (b)

t DU

9 40

1.0x10 10

1000

8.0x10

6.0x10 (Hz)

T (K)

100

4.0x10

2.0x10

0 50 100 150 200 250 300 0 50 100 150 200 250 300

t (d) t (d)

8 FIG. 3: The evolution results obtained by setting Tc = 5 × 10 K. Bumps in the LCs can be produced for the values of B¯ t and RDU taken as shown in the legends. The dotted line in panel (d) represents the critical temperature Tc.

14 Bd = 2 × 10 G and larger initial spin periods of Pi = 2 8 44 ms. The critical temperature is taken to be Tc = 5 × 10 10 K. The magnetars are assumed to have a DU core of radius 4 RDU = 5 × 10 cm, while their toroidal magnetic ﬁelds B¯ t are

10 diﬀerent. For comparison, we also show the results calculated 16 by taking Pi = 1 ms and B¯ t = 2 × 10 G, but other quantities

are kept unchanged (red solid lines). Obviously, for magne- 10 (erg / s )

rad tars with lower Bd and larger Pi, bumps still appear in the L resultant LCs after the maxima if the magnetars have toroidal 41

10 16 ﬁelds B¯ t > 1.9 × 10 G. The amplitudes of the bumps are

very sensitive to the strengths of B¯ t, and a positive correlation

10 exists between the former and the latter. Thus, bumps appear in the LCs only for magnetars with strong enough toroidal ﬁelds. The most interesting result is that compared to the no 80 toroidal ﬁeld case (B¯ t = 0), the strong B¯ t in a magnetar could enhance the peak luminosity and the emission after the peak

B =2.00E16 G, P =2 ms of a SLSN, rather than reduce them due to GW emission as t i

B =1.95E16 G, P =2 ms formerly considered [41, 48], if the tilt angle evolution is in-

t i 40 (degree) volved. Furthermore, with other quantities kept the same, a B =1.90E16 G, P =2 ms

t i smaller initial spin period of Pi = 1 ms does not lead to a

B =0.00 G, P =2 ms

20 t i bump in the LC, as seen in Fig. 4. This is because fast rota-

B =2.00E16 G, P =1 ms t i tion of the magnetar can reduce the damping timescale τd [Eq.

0 (3)], so χ = π/2 can be achieved very soon.

0 50 100 150 200 250 300

t (d) V. CONCLUSION AND DISCUSSIONS

FIG. 4: Evolutions of the radiated luminosity Lth and tilt angle χ 14 ¯ calculated for central magnetars with Bd = 2 × 10 G, while their Bt Newly born magnetars are generally considered to possess and P are different as shown in the legends. The light-gray curves i small initial spin periods (Pi & 1 ms), strong internal toroidal show the results derived by assuming a constant angle χ = π/2. magnetic fields, and initially small tilt angles. The last can 8 grow to π/2 due to damping of the free-body procession of a We note that, in addition to MDR and magnetically in- magnetar by internal viscosities. In this paper, by involving duced GW emission, newly born rapidly rotating magnetars the tilt angle evolution, we calculated the magnetar-powered may also spin down due to some instabilities that are driven SLSNe LCs. We find that, depending on Tc, RDU, B¯ t, and Pi, by the emission of GWs [69, 70], for instance, secular in- at the point when the tilt angle of a magnetar grows to π/2, stability [71], f-mode instability [72], and r-mode instability a bump could appear in the resultant LC after the maximum [73]. The first two are inclined to arise in more massive stars light. The arising of the bump can be associated with the oc- [57, 74], while the last one may arise in NSs with various 3 currence of a P2 neutron superfluid in the magnetar interior, masses. If a newly born magnetar can be considerably spun thus furthering the critical temperature Tc and the cooling of down through GW emissions associated with these instabili- the magnetar. We find that for newly born magnetars with ties in an extremely short time after birth, a bump may appear 14 16 Bd = 5 × 10 G, B¯ t = 4.6 × 10 G, and Pi = 1 ms there in the resultant LC due to a larger spin period of the magne- 9 will be no bumps in the LCs if Tc is as high as 2 × 10 K, or tar, as inferred from Fig. 4 by increasing Pi from 1 to 2 ms. 5 they have large DU cores with radii RDU = 1.5 × 10 cm. A Hence, whether bumps could appear depends on the braking lower Tc can result in a bump with larger relative amplitude at effects of the GW emissions that related to these instabilities. later times. Similarly, a smaller RDU leads to a later but more Specifically, for the f-mode and r-mode, their rather uncertain dim bump, and a lower peak luminosity. The most interesting saturation amplitudes are one of the key ingredients that deter- result is that the presence of strong toroidal magnetic field B¯ t mine the braking effects [74]. Detailed calculations with the in a newly born magnetar does not lower the peak luminos- effects of these instabilities involved are the aim of our future ity of the LC when the tilt angle evolution is involved. In- work. ¯ stead, a stronger Bt actually leads to both a brighter peak and In present paper, we adopted canonical values M = 1.4M⊙ a brighter bump. Moreover, for newly born magnetars with and R = 12 km for newly born magnetars. For more com- 14 Bd = 2 × 10 G and Pi = 2 ms, bumps still arise in the re- pact (massive) stars, larger DU cores may exist in their interi- ¯ sultant LCs when certain values for Bt, RDU, and Tc are taken. ors, especially when a realistic EOS is considered [54]. This While keeping other quantities unchanged, a smaller Pi does could expedite the cooling of the newly born magnetars, and not actually result in a bump in the LC because of a reduction possibly no bumps will appear in the LCs. Moreover,for more in the damping timescale τ . We therefore suggest that if the d massive NSs with larger DU cores, the damping timescale τd SLSNe LCs with such kinds of bumps could be observed in should be recalculated, though we expect that τd may be con- the future, by fitting these LCs with the model presented in siderably reduced due to stronger BV of stellar matter in the this paper, one can determine not only Bd and Pi of the newly presence of DU process. In future work, by adopting some re- ¯ born magnetars but also Bt, RDU, and Tc that are not easy to alistic EOS, we will calculate the structures and thermal evo- probe in other ways. The latter three quantities are crucial for lutions of NSs with various masses. Based on these results, our understanding of internal magnetic fields of magnetars, the magnetar-powered SLSNe LCs will be revisited, and by EOSs, and neutron superfluidity at supranuclear densities. comparing them with the observed ones, we may get some in- Until now, none of the observed SLSNe shows such kinds formation about the physical parameters of newly born mag- of bumps in its LCs. This may be due to the following rea- netars as well as EOS of dense matter. 9 sons: (i) A high critical temperature of Tc > 1.5 × 10 K is 16 favored. (ii) The central newly born magnetars have B¯ t . 10 G, which may be attributed to a less effective field amplification process related to MRI or an α − ω dynamo. (iii) The Acknowledgments DU process plays a dominant role in the cooling of the central 5 magnetars (RDU & 10 cm), and therefore EOS PAL is more favorable as compared to EOS APR. (iv) The magnetars have We gratefully thank the anonymous referee for constructive 16 small Pi (of 1 ms) but with relatively weak B¯ t (of 2 × 10 G), and helpful comments and suggestions for improving this pa- and fast growths of the tilt angles in extremely early periods per. We also thank Shao-Ze Li for useful discussions. Quan are thus unavoidable, as presented in Fig. 4. Consequently, in Cheng acknowledges funding support by China Postdoctoral the current status, only very rough constraints may be set on Science Foundation under grant No. 2018M632907. This these key physical quantities of magnetars. In order to obtain work is also supported by the National Natural Science Foun- rigorous constraints, more observations of SLSNe with higher dation of China (Grants No. 11773011, No. 11373036, No. precision are needed so as to find such bumps in the LCs. On 11133002, No. 11473008, and No. 11622326), the National the other hand, further improvements of our model are also Program on Key Research and Development Project (Grants necessary (see the discussion in the last paragraph of this sec- No. 2016YFA0400802, and No. 2016YFA0400803), and the tion). In a word, observation of SLSNe LCs may provide a Key Research Program of Frontier Sciences, CAS (Grant No. new approach to probe the physics of newly born magnetars. QYZDY-SSW-SLH008).

[1] Usov, V. V. 1992, Nature, 357, 472; Dai, Z. G., Lu, T. 1998, 552, L35; Metzger, B. D., Giannios, D., Thompson, T. A., Buc- Phys. Rev. Lett., 81, 4301; Zhang, B., M´esz´aros, P. 2001, ApJ, 9

ciantini, N., Quataert, E. 2011, MNRAS, 413, 2031 [33] Cutler, C. 2002, Phys. Rev. D, 66, 084025 [2] Campana, S., et al. 2006, A&A, 454, 113; Metzger, B. D., [34] Dall’Osso, S., Shore, S. N., Stella, L. 2009, MNRAS, 398, 1869 Quataert, E., Thompson, T. A. 2008, MNRAS, 385, 1455; [35] Dall’Osso, S., Giacomazzo, B., Perna, R., Stella, L. 2015, ApJ, Rowlinson, A., et al. 2013, MNRAS, 430, 1061; Gompertz, B. 798, 25 P., O’Brien, P. T., Wynn, G. A., Rowlinson, A. 2013, MNRAS, [36] Marassi, S., Ciolfi, R., Schneider, R., Stella, L. Ferrari, V. 2011, 431, 1745 MNRAS, 411, 2549 [3] Woosley, S. E. 2010, ApJL, 719, L204 [37] Cheng, Q., Zhang, S. N., Zheng, X. P. 2017, Phys. Rev. D, 95, [4] Kasen, D., Bildsten, L. 2010, ApJ, 717, 245 083003 [5] Inserra, C., Smartt, S. J., Jerkstrand, A., et al. 2013, ApJ, 770, [38] Cheng, Q., Yu, Y. W., Zheng, X. P. 2015, MNRAS, 454, 2299 128 [39] Chamel, N., Haensel, P. 2008, Living Rev. Relativ., 11, 10 [6] Yu, Y. W., Zhang, B., Gao, H. 2013, ApJL, 776, L40 [40] Alpar, A., Sauls, J. A. 1988, ApJ, 327, 723 [7] Li, S. Z., Yu, Y. W. 2016, ApJ, 819, 120 [41] Ho, W. C. G. 2016, MNRAS, 463, 489 [8] Yu, Y. W., Li, S. Z., Dai, Z. G. 2015, ApJL, 806, L6 [42] Baldo, M., et al. 1998, Phys. Rev. C, 58, 1921 [9] Zhang, B. 2014, ApJL, 780, L21; Ravi, V., Lasky, P. D. 2014, [43] Schwenk, A., Friman, B. 2004, Phys. Rev. Lett. 92, 082501; MNRAS, 441, 2433; Piro, A. L. 2016, ApJ, 824, L32; Metzger, Khodel, V. A., et al. 2004, Phys. Rev. Lett., 93, 151101 B. D., Berger, E., Margalit, B. 2017, ApJ, 841, 14; Kashiyama, [44] Page, D., Lattimer, J. M., Prakash, M., Steiner, A. W. 2004, K., Murase, K. 2017, ApJ, 839, L3 ApJS, 155, 623 [10] Quimby, R. M., et al. 2011, Nature, 474, 487 [45] Page, D., Prakash, M., Lattimer, J. M., Steiner, A. W. 2011, [11] Gal-Yam, A. 2012, Science, 337, 927 Phys. Rev. Lett., 106, 081101; [12] Chatzopoulos, E., Wheeler, J. C., Vinko, J., Horvath, Z. L., [46] Shternin, P. S., Yakovlev, D. G., Heinke, C. O., Ho, W. C. G., Nagy, A. 2013, ApJ, 773, 76; Patnaude, D. J. 2011, MNRAS, 412, L108 [13] Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2013, Nature, [47] Beloin, S., Han, S., Steiner, A. W., Page, D. 2018, Phys. Rev. 502, 346; Howell, D. A., et al. 2013, ApJ, 779, 98; Nicholl, C, 97, 015804 M., Jerkstrand, A., Inserra, C., et al. 2014, MNRAS, 444, 2096; [48] Kashiyama, K., Murase, K., Bartos, I., Kiuchi, K., Margutti, R. Vreeswijk, P. M., Savaglio, S, Gal-Yam, A., et al. 2014, ApJ, 2016, ApJ, 818, 94 797, 24; Wang, S. Q., Wang, L. J., Dai, Z. G., Wu, X. F. [49] Friman, B. L., Maxwell, O. V. 1979, ApJ, 232, 541 2015, ApJ, 799, 107; Metzger, B. D., Margalit, B., Kasen, D., [50] Lattimer, J. M., Pethick, C. J., Prakash, M., Haensel, P. 1991, Quataert, E. 2015, MNRAS, 454, 3311; Kasen, D., Metzger, B. Phys. Rev. Lett., 66, 2701 D., Bildsten, L. 2016, ApJ, 821, 36; Lunnan, R., et al. 2016, [51] Page, D., Geppert, U., Weber, F. 2006, Nucl. Phys. A, 777, 497 ApJ, 831, 144; Chen, K. J., Woosley, S. E., Sukhbold, T. 2016, [52] Akmal, A., Pandharipande, V. R., Ravenhall, D. G. 1998, Phys. ApJ, 832, 73; Nicholl, M., et al. 2017, ApJL 835, L8 Rev. C, 58, 1804 [14] Yu, Y. W., Zhu, J. P., Li, S. Z., Lü, H. J., Zou, Y. C. 2017, ApJ, [53] Prakash, M., Ainsworth, T. L., Lattimer, J. M. 1988, Phys. Rev. 840, 12 Lett. 61, 2518 [15] Moriya, T. J., Tauris, T. M. 2016, MNRAS, 460, L55 [54] Alford, M. G., Mahmoodifar, S., Schwenzer, K. 2012, Phys. [16] Stella, L., DallOsso, S., Israel, G. L., Vecchio, A. 2005, ApJ, Rev. D, 85, 024007 634, L165 [55] Page, D., Applegate, J. H. 1992, ApJL, 394, L17 [17] Fan, Y. Z., Wu, X. F., Wei, D. M. 2013, Phys. Rev. D, 88, [56] Yakovlev, D. G., Pethick, C. J. 2004, Ann. Rev. Astron. Astro- 067304 phys., 42, 169 [18] Makishima, K., et al. 2014, Phys. Rev. Lett., 112, 171102 [57] Strobel, K., Schaab, C., Weigel, M. K. 1999, A&A, 350, 497 [19] Ferrario, L., Wickramasinghe, D. 2006, MNRAS, 367, 1323 [58] Heger, A., Langer, A., Woosley, S. E. 2000, ApJ, 528, 368; Ster- [20] Duncan, R. C., Thompson, C. 1992, ApJ, 392, L9; Thompson, gioulas, N. 2003, Living Rev. Relativ., 6, 3; Kissin, Y., Thomp- C., Duncan, R. C. 1993, ApJ, 408, 194 son, C. 2017, arXiv: 1705.07906 [21] Price, D. J., Rosswog, S. 2006, Science 312, 719 [59] Geppert, U., Rheinhardt, M. 2006, A&A, 456, 639 [22] Duez, M. D., Liu, Y. T., Shapiro, S. L., Shibata, M., Stephens, [60] Lattimer, J. M., Prakash, M. 2001, ApJ, 550, 426 B. C. 2006, Phys. Rev. D, 73, 104015; Guilet, J., Bauswein, A., [61] Braithwaite, J., Spruit, H. C. 2004, Nature, 431, 819; Braith- Just, O., Janka, H.-T. 2017, MNRAS, 471, 1879 waite, J., Spruit, H. C. 2006, A&A, 450, 1097; [23] Akiyama, S., Wheeler, J. C., Meier, D. L., Lichtenstadt, I. 2003, [62] Braithwaite, J. 2009, MNRAS, 397, 763; Akgün, T., Reiseneg- ApJ, 584, 954 ger, A., Mastrano, A., Marchant, P. 2013, MNRAS, 433, 2445 [24] Cheng, Q., Yu, Y. W. 2014, ApJL, 786, L13; Cheng, Q., Zhang, [63] Cutler, C., Jones, D. I. 2000, Phys. Rev. D, 63, 024002 S. N., Zheng, X. P. 2017, Research in Astron. and Astophys., [64] Zdunik, J. L. 1996, A&A, 308, 828 17, 54 [65] Mahmoodifar, S., Strohmayer, T. 2013, ApJ, 773, 140 [25] Rembiasz, T., Guilet, J., Obergaulinger, M., Cerdá-Durán, P., [66] Jones, P. B. 1976, Astrophys. Space Sci., 45, 369 Aloy, M. A., Müller, E. 2016, MNRAS, 460, 3316 [67] Shapiro, S. L., Teukolsky, S. A. 1983, Black holes, white [26] Mösta, P., Ott, C. D., Radice, D., Roberts, L. F., Schnetter, E., dwarfs, and neutron stars: The physics of compact objects Haas, R. 2015, Nature, 528, 376 (New York: Wiley-Interscience) [27] Guilet, J., Müller, E. 2015, MNRAS, 450, 2153 [68] Arnett, W. D. 1979, ApJ, 230, L37 [28] Halevi, G., Mösta, P. 2018, arXiv:1801.08943 [69] Chandrasekhar, S. 1970, Phys. Rev. Lett. 24, 611 [29] Metzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, [70] Friedman, J. L., Schutz, B. F. 1978, ApJ, 222, 281 N., Quataert, E. 2011, MNRAS, 413, 2031 [71] Lai, D., Shapiro, S. L. 1995, ApJ, 442, 259 [30] Ott, C. D., Burrows, A., Thompson, T. A., Livne, E., Walder, [72] Andersson, N., Kokkotas, K. D. 1998a, MNRAS, 299, 1059 R. 2006, ApJS, 164, 130 [73] Andersson, N. 1998b, ApJ, 502, 708 [31] Obergaulinger, M., Aloy, M. A. 2017, MNRAS, 469, L43 [74] Glampedakis, K., Gualtieri, L. 2017, arXiv:1709.07049v2 [32] Gilkis, A. 2018, MNRAS, 474, 2419