Draft version September 6, 2021 A Typeset using L TEX twocolumn style in AASTeX631 Fallback Accretion Halted by R-process Heating in Neutron Star Mergers and Gamma-Ray Bursts Wataru Ishizaki,1 Kenta Kiuchi,1,2 Kunihito Ioka,1 and Shinya Wanajo2,3 1Center for Gravitational Physics, Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 2Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Am M¨uhlenberg, Potsdam-Golm, 14476, Germany 3Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), RIKEN, Wako, Saitama 351-0198, Japan ABSTRACT The gravitational wave event GW170817 with a macronova/kilonova shows that a merger of two neutron stars ejects matter with radioactivity including r-process nucleosynthesis. A part of the ejecta inevitably falls back to the central object, possibly powering long-lasting activities of a short gamma- ray burst (sGRB), such as extended and plateau emissions. We investigate the fallback accretion with the r-process heating by performing one-dimensional hydrodynamic simulations and developing a semi-analytical model. We show that the usual fallback rate dM/dt t−5/3 is halted by the heating ∝ because pressure gradients accelerate ejecta beyond an escape velocity. The suppression is steeper than Chevalier’s power-law model through Bondi accretion within a turn-around radius. The characteristic halting timescale is 104–108 sec for the GW170817-like r-process heating, which is longer than ∼ the typical timescale of the long-lasting emission of sGRBs. The halting timescale is sensitive to the uncertainty of the r-process. Future observation of fallback halting could constrain the r-process heating on the year scale. Keywords: hydrodynamics — accretion, accretion discs — nucleosynthesis — gamma-ray burst: gen- eral 1. INTRODUCTION itational waves and electromagnetic counterparts are The discovery of the short-duration gamma-ray broadly consistent with the predictions from a series of burst (sGRB), GRB 170817A, coincided with the theoretical studies including numerical relativity calcu- detection of GW170817 by the Laser Interferome- lations (Shibata et al. 2017; Cowperthwaite et al. 2017; ter Gravitational-Wave Observatory (LIGO) and the Kasen et al. 2017; Kasliwal et al. 2017; Villar et al. Virgo Consortium (LVC), is a direct evidence that 2017; Ioka & Nakamura 2018, 2019) and these si- binary neutron star (BNS) mergers are one of the multaneous detection has established the paradigm sources of sGRBs (Abbott et al. 2017a,b). The si- that sGRBs originate from compact binary mergers multaneous detection of macronova/kilonova emis- (Paczynski 1986; Goodman 1986; Eichler et al. 1989). sion is a strong indication for the ejection of Although sGRBs are classified as GRBs with the arXiv:2104.04708v3 [astro-ph.HE] 3 Sep 2021 neutron-rich matter and r-process element synthesis in duration of prompt emission being less than 2 sec this BNS merger (Arcavi et al. 2017; Chornock et al. (Kouveliotou et al. 1993), their central engines are 2017; Coulter et al. 2017; Cowperthwaite et al. 2017; thought to remain active for longer time (Burrows et al. Drout et al. 2017; Kasen et al. 2017; Kasliwal et al. 2005; Ioka et al. 2005). This was suggested from the 2017; Kilpatrick et al. 2017; McCully et al. 2017; observations of extended emission which lasts for about 48 49 −1 Nicholl et al. 2017; Shappee et al. 2017; Smartt et al. 100 sec with a luminosity of 10 –10 erg s and pos- 2017; Soares-Santos et al. 2017; Tanaka et al. 2017; sibly from the observations of plateau emission which 4 45 Tanvir et al. 2017). The observed properties of the grav- lasts for about 10 sec with a luminosity of 10 – 1046 erg s−1 (Barthelmy et al. 2005; Rowlinson et al. 2013; Kisaka et al. 2017; Kagawa et al. 2019). These Corresponding author: Wataru Ishizaki long-lasting emissions are fainter than the prompt [email protected] emission in luminosity. However, because of the long duration of these emissions, they are compara- 2 Ishizaki et al. ble with or greater than the prompt emission in flu- more sophisticated calculation based on the model in ence (Kisaka et al. 2017). Some of the extended emis- Metzger et al. (2010), using the ejecta profiles obtained sions have been observed to darken abruptly, which are from numerical relativity simulations and the radioac- not expected in the sGRB afterglow (Ioka et al. 2005). tive heating rates obtained from nucleosynthesis calcu- Therefore, a central engine activity is very likely to ex- lations. They showed that the mass accretion stops and plain these long-lasting emission1. There are two mod- resumes after a finite time, the so-called “gap”, because els for the energy source of the long-lasting emissions, the marginally bound fluid elements become unbound which have been extensively discussed in previous works. by the heating due to the radioactive decay of r-process One is the release of rotational or magnetic field en- elements. Furthermore, they argued that the timescale ergy from a strongly magnetized massive neutron star of this resumption is (100) sec, which is in agreement O that forms after a BNS merger (Zhang & M´esz´aros 2001; with the timescale of the extended emission. Metzger et al. 2008; Bucciantini et al. 2012; Fan et al. In the test-particle model of Metzger et al. (2010) and 2013; Murase et al. 2018). Currently, this scenario is Desai et al. (2019), they assumed that all the radioac- not conclusive because the expected observational fea- tive energy is converted to the kinetic energy. In reality, tures, e.g., late radio emission, have not been detected the energy from the radioactive heating is converted to (Metzger & Bower 2014; Horesh et al. 2016). the kinetic energy through the pressure gradient force The other scenario for the energy source of resulting from the increased internal energy. It is un- long-lasting activity is the fallback accretion clear whether the pressure gradient is large enough to of ejecta (Rosswog 2007; Lee & Ramirez-Ruiz convert all the internal energy of the fluid element into 2007; Rossi & Begelman 2009; Kisaka & Ioka 2015; the kinetic energy, so that the validity of their assump- Kisaka et al. 2017). Numerical calculations of com- tion is also unclear. Because their assumption cannot be pact binary mergers show that a part of the ejecta verified within the framework of the test-particle model, is still gravitationally bound (e.g., Rosswog et al. 1999; it is necessary to solve the hydrodynamic equations in- Bauswein et al. 2013; Kyutoku et al. 2015; Radice et al. corporating the effects of radioactive heating. 2016; Kiuchi et al. 2017). Kyutoku et al. (2015) calcu- In this paper, we reconsider the effect of the radioac- lated the coalescence of a black hole and a neutron star tive heating due to decaying r-process elements on the (BH–NS) based on numerical relativity and showed that fallback accretion in BNS mergers by numerically solv- a part of the ejecta falls back to the merger remnant. ing one-dimensional fluid equations. We also construct While the gravitational energy released by fallback ac- a semi-analytical model that reproduces the hydrody- cretion is large enough to explain the extended emission namical simulation results and explores it over a large and the plateau emission, the simple theory of fallback parameter space including the realistic r-process heat- accretion, which assumes a zero-temperature fluid, pre- ing. This paper is organized as follows. In Section 2, dicts that the mass accretion rate is proportional to we describe the method of the numerical calculation of the power of time as t−5/3 with no typical timescale the fluid equations and show the results. In Section 3 (Rees 1988; Michel 1988). However, the observed light and Section 4, the semi-analytical models of the accre- curve of the long-lasting emission of sGRB clearly has tion flow are developed for a constant heating rate and a certain timescale, which is not compatible with the a heating rate of a broken power-law form, respectively. simple theory (Kagawa et al. 2019). In Section 5, the semi-analytical model is applied to ex- The coincidence of the macronova/kilonova emission plore the parameters relevant to the radioactive heating. with the gravitational wave source GW170817 indi- In Section 6, we summarize this work, and discuss the cates that the ejecta of the BNS merger is heated by scope of application of the spherically symmetric mod- the radioactive decay of r-process elements. There- eling used in this study and future prospects. fore, the assumption of zero-temperature fluid is in- appropriate (Smartt et al. 2017; Waxman et al. 2018; 2. HYDRODYNAMICAL SIMULATION OF Kawaguchi et al. 2018). Metzger et al. (2010) discussed FALLBACK ACCRETION the effect of ejecta heating by the radioactive decay 2.1. Numerical method of r-process elements on the mass accretion using a In order to investigate the effect of the radioactive test-particle model. Desai et al. (2019) performed a heating due to r-process nuclei on the fallback accre- tion in the BNS merger, we have performed long-term 1 Note that, although such a rapid fade-out has not been reported one-dimensional hydrodynamic simulations of the mat- for the plateau emission, it may also reflect the late-time activity ter ejected during the merger. The ejecta profiles of the of the central engine (see also Matsumoto et al. 2020). velocity and the mass density are derived from the nu- Fallback Halted by R-process Heating 3 merical relativity simulations performed by Kiuchi et al. 1 1012 (2017). The hydrodynamical equations for spherically symmetric and purely radial flow, including the heating 1010 and the point source gravity, are written as follows: 108 ∂ρ 1 ∂ vr (Model) 2 ρ (Model) + r ρv =0, (1) θ 2 [c] 0.1 =0 deg r ∂t r ∂r v 6 θ=30 deg 10 θ=60 deg θ=90 deg vesc ∂v ∂v 1 ∂P GM 4 + v = , (2) 10 rest mass density [g/cc] ∂t ∂r −ρ ∂r − r2 102 ∂ 1 1 ∂ 1 0.01 ρv2 + ǫ + r2v ρv2 + ǫ + P 100 1000 ∂t 2 int r2 ∂r 2 int Radius [km] GM = ρv + Q , (3) − r2 heat Figure 1.