arXiv:2101.12220v2 [astro-ph.HE] 13 Apr 2021 10 al se[6 7 o ueia iuain)b adopting by simulations) analyti- numerical estimated for by often 17] star [16, is the (see hole of cally revisited rate black accretion been endoparasitic subsequent recently the The have the 13]. star, near [12, PBH neutron by the the of down via of settling center gravitational energy the to and kinetic lead surface that PBH waves accretion, of cap- friction, dissipation the dynamical focusing, the for gravitational via by Estimates stars followed time-scale neutron by neutron-star reliable PBHs processes. of oldest on ture these the hinges for of therefore estimates age it the time- on populations; than accretion- proceeds subsequent shorter star neutron and scales the capture of the destruction if driven only PBHs on for discussed, 15]. been [14, have in process observa-example, this Other to of Universe. contribution signatures the their tional of on content hence dark-matter and the the on PBHs, that limit so of a [13]), density provides also the This see populations but neutron-star neu- 12], of 11]). by [6, existence [1, capture e.g., (see, for e.g., stars efficient tron (see, particularly ac- be stars then may these would process and swallow stars, and by fact captured the crete be from results can 10]) PBHs [4, that e.g., see, observations; other swl s[]frarcn eiw.Cntanso the Sufficiently estab- ( on observations. been PBHs Constraints different have of small number matter review). a dark recent by the lished a to for contribution [3], [4] PBH e.g., as (see, matter well dark accounting as elusive not and if mysterious con- the to, been for, long contributing have for 2]) [1, candidates e.g., sidered (see, Universe early the in otedr atri h asrne10 range mass the in matter dark the to detection possible a for also [9] gravitational- NANOGrav). see as from by well therein; arising as constraints references observations, of wave and review 7] a mass for [6, different [8] e.g., for limits (see, types in ranges Hubble different resulted a masses, have observations than black-hole of larger less for within while, [5] time, radiation Hawking to due yee yREVT by Typeset vdnl,teaoeagmn a rvd constraints provide can argument above the Evidently, rmrilbakhls(Bs htmyhv formed have may that (PBHs) holes black Primordial n opligcntan ntePHcontribution PBH the on constraint compelling One − 10 M eto tr abrn rmrilBakHl:MxmmSu Maximum Hole: Black Primordial a Harboring Stars Neutron ⊙ 3 eateto srnm n CA nvriyo liosat Illinois of University NCSA, and Astronomy of Department wihcretyi o elcntandby well-constrained not is currently (which h lc oe eas hwta hsmxmmsria ied time survival Γ maximum with this EOSs that st whic show polytropic for after also for We time earlier survival hole. found black maximum we the the that analytically flow derive Correspon adiabatic we center. for their rate at accretion hole black primordial, possibly akmte otn fteUniverse. the of content dark-matter h otiuino rmrilbakhlsi h asrange mass the in holes black s Establishing primordial alone. of mass contribution hole black the initial the on depends M eepoei eea eaiiytesria ieo neutro of time survival the relativity general in explore We 2 eateto hsc,Uiest fIlni tUrbana-Ch at Illinois of University Physics, of Department . E 1 X eateto hsc n srnm,BwonClee Brun College, Bowdoin Astronomy, and Physics of Department 5 × 10 14 )wudhv evaporated have would g) hmsW Baumgarte W. Thomas ≥ 5 − / 16 ,s htti uvvltm sue eryuieslvalu universal nearly a assumes time survival this that so 3, M ⊙ . M . 1 n tatL Shapiro L. Stuart and h peia,sed-tt,Bniaceinformula, accretion Bondi steady-state, spherical, the et nldn t eaiitcgnrlzto) Here treat- generalization). textbook relativistic a for its [19] including also ment, see ([18]; flow adiabatic for sue ob elgbei oprsnwt h black the with parameter is The comparison radius accretion in still mass. the negligible hole holes, within be gas black to the of small assumed mass central for The homogeneous nearly core. but, star’s neutron distances hole, the inside large well black at the measured from respectively, speed, sound ic h rdc ftetoterms two the of product the Since and infinity, at and where inegnau” hc,fr1 for which, eigenvalue”, tion o sgvnb h lpe function “lapse” the fac- to by redshift given due This rate is infinity tor core. a the at at from redshift measured advances gravitational the that time than This slower so- [20]). somewhat flow spacetime footnote Bondi the see relativistic of the derivation lution; of full applicability a and (see for core metric [17] star in neutron the A inside mea- Appendix the well as from again far time but observer hole dot to static black the asymptotic” respect star, “local with a neutron derivative by a sured a of denotes center (1) the at in hole black a onto unity. order esrdb 1 sfnaetlya(baryon) a fundamentally is (1) by measured oi bevr yial 0 typically observer, totic dM mass ino diinlitra energy, internal additional of tion se[7) conigfrbt fteefcos emay we as factors, (1) these rate of accretion both the for rewrite Accounting [17]). (see eut nadmninesnme eycoet unity, to close very number dimensionless αβ a in results nagnrlrltvsi omlto,adfraccretion for and formulation, relativistic general a In rest ≈ a M igt h iiu taysaeBondi steady-state minimum the to ding ,w ilinr hspouti h following. the in product this ignore will we 1, t hc nacstebakhl’ gravitational hole’s black the enhances which , c au sipratfrconstraining for important is value a uch r h aypoi”rs-asdniyadthe and density rest-mass “asymptotic” the are o esrstetm sosre ya observer an by observed as time the measures now raaCapin raa lios61801 Illinois Urbana, Urbana-Champaign, h niesa ilb osmdby consumed be will star entire the h ffncereutoso tt (EOSs), state of equations nuclear iff yasmwa agraon u oteaccre- the to due amount larger somewhat a by 10 tr hths nendoparasitic, an host that stars n mag,Ubn,Ilni 61801 Illinois Urbana, ampaign, pnsol ekyo h stiffness the on weakly only epends − 16 M M ⊙ dM stebakhl’ rvttoa mass. gravitational hole’s black the is dt wc,Mie04011 Maine swick, . M ,3 2, ˙ 4 = M 4 = . . παβλ 6 πλ . λ 10 ≤ M sadmnines“accre- dimensionless a is < α − a 10 M Γ 2 3 ρ M a ≤ vvlTime rvival 2 0 α .Teacee mass accreted The 1. 3 β ρ ⊙ 5 0 o h oa asymp- local the for α ≡ / , othe to that e ,i osatof constant a is 3, . dM/dM and 1 rest rest β mass & & (1) (2) ρ 1 1 0 2

For accretion onto black holes, a relativistic treatment for the maximum time that a can survive for Schwarzschild black holes shows that the requirement accretion by an endoparasitic . that the sound speed be less than the speed of light de- We first observe that, for Γ = 5/3, the accretion rate mands that the flow pass through a critical point, yield- (6) becomes independent of the sound speed a. A rel- ing a unique value for λ [19]. Adopting a polytropic equa- ativistic treatment shows that this value represents a tion of state (EOS) minimum accretion rate, since relativistic corrections for higher sound speeds a . 1 and correspondingly larger Γ P = Kρ0 , (3) densities lead to higher accretion rates. As shown in [22], Bondi accretion in general relativity exhibits such where K is a constant, Γ = 1+1/n is the adiabatic a minimum accretion rate for stiffer EOSs as well. This index and n is the polytropic index, that unique value minimum rate is given by for 1 ≤ Γ ≤ 5/3, assuming a ≪ 1, is dM 4πλ¯ M 2 (5−3Γ)/2(Γ−1) = , (7) 1 2  dt  Γ1/(Γ−1) K1/(Γ−1) λ = , (4) min 4 5 − 3Γ for all 5/3 ≤ Γ ≤ 3. Here λ¯ is defined as which is the value found in the Newtonian formulation. Note that here and throughout we adopt geometrized − 1/(Γ−1) 3/2 ¯ (5−3Γ)/(2Γ−2) Γ 1 (1+3¯xIII) units with G = 1 = c. λ =x ¯III , Γ − 1 − x¯III  4 While the accretion timescale may be small compared (8) to the total capture and settling time for many black with hole masses, there is some uncertainty in determining 1/2 the actual accretion timescale when applying the Bondi x¯III =Γ − 7/6 − (12Γ − 11) /6, (9) accretion rate given by Eq. (2). We first observe that (2) depends on both ρ and a, which may vary from star to where we have adopted the notation of [22]. The exis- star (even though their order of magnitude is probably tence of the minimum accretion rate (7) results from the similar for all neutron stars). More importantly, the ac- fact that, for Γ > 5/3 and to leading order in a, the accre- cretion eigenvalues λ are easily derived only for soft EOSs tion eigenvalues λ can be written as λ = λa¯ (3Γ−5)/(Γ−1), with 1 ≤ Γ ≤ 5/3, as in Eq. (4). As a result, some au- so that the dependence on a in (6) cancels out, thereby thors (e.g. [12, 21]) have resorted to approximating the leaving M˙ finite but non-zero. Note that (7) reduces to accretion rate by adopting values for Γ such as Γ = 4/3, (6) for Γ=5/3, in which case λ = λ¯ =1/4. even though one expects that the neutron star interior is Since the black hole accretes matter at a rate greater governed by a significantly stiffer EOS with Γ & 2, for than the minimum rate (7), we can find the maximum which (4) clearly does not apply. Finally, to make mat- time by integrating [23] ters worse, consider that the polytropic EOS (3) implies that the sound speed is given by Mtot dM Mtot dM tmax = ≤ ZM0 dM/dt ZM0 (dM/dt)min a2 ≃ ΓKρΓ−1 (5) 0 Mtot 1/(Γ−1) dM = κK 2 whenever P ≪ ρ, where ρ is the total mass-energy den- ZM0 M sity. This is the standard Newtonian relation. Inserting 1/(Γ−1) 1 1 (5) into (2) results in = κK − . (10) M0 Mtot  dM ≃ 4πλM 2a(5−3Γ)/(Γ−1). (6) Here M0 is the initial black-hole mass, Mtot is the total dt gravitational mass of the system, and we have defined

→ − We now observe that, for a 0, (6) suggests that dM/dt Γ1/(Γ 1) becomes infinite whenever Γ > 5/3, seemingly making κ = ¯ . (11) this expression totally unreliable for stiff EOSs. 4πλ In this short paper we clarify this issue. Building on Assuming that M0 ≪ Mtot we may neglect the last term our relativistic treatment of Bondi accretion for stiff poly- in (10) and obtain tropic EOSs (see [22]) we use our finding there that there exists a finite, nonzero, minimum accretion rate whenever K1/(Γ−1) Γ ≥ 5/3 to show that, under quite general conditions, tmax = κ . (12) M0 there exists a corresponding maximum accretion time for a PBH to swallow a neutron star. This maximum accre- This maximum time is dominated by the initial phase tion time depends only weakly on the stiffness for EOSs of accretion, when M and the accretion rate given by with Γ ≥ 5/3, thereby providing a nearly universal es- Eq. (7) assume their lowest values. During this phase timate, depending only on the initial black-hole mass, the accretion proceeds in a quasistationary fashion and 3

¯ ¯ ¯ 2 n Γ λ κ Mmax κ/Mmax We next observe that, according to Fig. 4 of [22], ac- 1.5 1.67 0.25 0.685 0.265 9.75 cretion rates for Γ > 5/3 depend only very weakly on the 1.4 1.71 0.432 0.392 0.237 6.98 asymptotic sound speed, meaning that the actual accre- 1.3 1.77 0.640 0.261 0.213 5.75 tion rate is comparable to its minimum value, and hence 1.2 1.83 0.889 0.185 0.194 4.92 the actual accretion time is comparable to its maximum value (16). 1.1 1.91 1.17 0.138 0.177 4.40 For concreteness, we shall take the maximum mass of 1.0 2.0 1.49 0.107 0.164 3.96 a nonrotating neutron star to be Mmax = 2.3M⊙. This 0.9 2.11 1.83 0.085 0.152 3.69 mass is consistent with the likely range of values esti- 0.8 2.25 2.15 0.071 0.142 3.46 mated by several groups (see, e.g., [26–29]) using data 0.7 2.43 2.45 0.060 0.135 3.32 from the binary event GW170817 de- 0.6 2.67 2.69 0.053 0.129 3.21 tected by LIGO/Virgo in gravitational waves [30], as well 0.5 3.0 2.83 0.049 0.125 3.12 as from the counterpart gamma-ray burst GRB17081A and AT 2017gfo. It is also consistent with the TABLE I. Polytropic parameters. Here n is the polytropic measurements of [31], which revealed the highest index, Γ = 1+1/n the adiabatic index, the coefficients λ¯ and mass to date. We then may rewrite (16) as −2/n κ are defined in Eqs. (8) and (11), and M¯ max = K Mmax −10 22 is the maximum dimensionless gravitational mass allowed 10 M⊙ 10 g ≃ × 6 ≃ × 7 for non-rotating neutron stars. Note that the combination tmax 1 10 s 2 10 s , (17) 2  M0   M0  κ/M¯ max, which appears in the maximum accretion time (16), depends only weakly on Γ for stiff EOSs (with Γ & 2, say). resulting in values that are remarkably close to those adopted by, for example, [12] (see their Eq. 39). This close agreement, though reassuring, is rather coinciden- the neutron star maintains its initial state (e.g., Γ and tal, as our value results from a detailed treatment of K) to good approximation. Bondi accretion for stiff EOSs in full general relativity. We can now determine the value of the constant K for Regardless, results quite similar to (17) have been in- any polytropic EOS by matching its dimensionless max- voked by previous investigators (see, e.g., [6, 12] and ref- −1/(2Γ−2) imum mass M¯ max = K Mmax to the best cur- erences therein) to constrain PBHs in the mass range −16 −10 rent value value for the maximum mass of a nonrotating, 10 M⊙ . M . 10 M⊙ as dark-matter candidates. isolated neutron star, Mmax. The quantity M¯ max is ob- Corresponding to Eqs. (14) and (17), the minimum tained by setting K equal to unity and integrating the accretion rate given by Eq. (7) can be evaluated to yield Tolman-Oppenheimer-Volkoff equations (see [24, 25]) to initially find the mass at the turning point along an equilibrium ¯ 2 2 2 sequence of stars parametrized by their central density. 5 M⊙ Mmax M⊙ M0 M˙ min =2 × 10 We can therefore express K as s κ Mmax  M⊙  2 2(Γ−1) −9 M⊙ M0 Mmax . 3 × 10 − , (18) K = . (13) yr 10 10M⊙  M¯ max  where, in the last step, we have again assumed (15) and Inserting (13) into (12) we now obtain Mmax =2.3M⊙. To summarize, we provide estimates for the accretion 2 κ M M⊙ times of black holes residing in neutron stars. While our t =5 × 10−6s max , (14) max ¯ 2 results are very similar to previously adopted values, they Mmax  M⊙   M0  are based on a rigorous, relativistic treatment of spherical where we have used that, in geometrized units, 1M⊙ ≃ Bondi accretion for stiff EOSs. We also demonstrate that 5 × 10−6 s. there exists a maximum accretion time, depending chiefly We list values for the polytropic parameters appear- on the initial black hole mass alone, by which time the ing in (14) in Table I. Note that, for moderately stiff black hole will have consumed the entire neutron star. ¯ 2 EOSs with n . 1.0 and Γ & 2, say, the factor κ/Mmax We further argue that actual accretion times will not varies only very moderately, with values between 3 and differ much from the maximum accretion time, so that 4. Adopting the maximum of these values, the latter provides an approximate universal value for the lifetime of a neutron star with a small black hole κ ≃ 4, (15) residing at its center. M¯ 2 While our discussion here centers on PBHs and their max max contribution to the dark matter, the same estimates ap- we can rewrite (12) as ply to some alternative scenarios as well. Specifically, the neutron star could also capture other dark matter 2 −10 5 Mmax 10 M⊙ particles which, under sufficiently favorable conditions, tmax ≃ 2 × 10 s . (16)  M⊙   M0  could form a high density object that collapses to form 4 a small black hole in the neutron star interior and ulti- could, in principle, be phase transitions and other effects mately consume the entire star. Invoking observations of that might occur in the accretion flow, but we note that neutron star populations, several authors have used these the Bondi critical radius occurs rather close to the black arguments to derive constraints and bounds on dark mat- hole for typical cases (see Fig 2 in [22]). This means ter particles (see, e.g., [32–36]). that the density and temperature in the accreting gas do Our arguments build on a number of assumptions, of not increase substantially above their core values before course. We assume that the accretion process is dom- the gas is captured. In spite of these considerations, we inated by spherical, steady-state, adiabatic Bondi ac- believe that our result provides an interesting limit on cretion (which may not hold for extremely stiff EOSs, the timescale for the demise of a neutron star by an en- e.g. polytropes with Γ > 3; see, for example, [22, 37]), doparasitic black hole, not least because it provides more and we ignore effects of rotation (which seems justified rigorous justification for a key assumption used in con- according to the findings of [16]; see also [11, 21]). We straining the contribution of PBHs to the dark matter also ignore radiation and radiation pressure, although the content of the Universe. trapping radius for photons [38] is likely very large and radiation may be totally ineffective in holding back the hyper-Eddington accretion that can arise here (see, e.g., ACKNOWLEDGMENTS [39]). Radiation processes, including the possible role of neutrinos, as well as the role of conduction, need to It is a pleasure to thank Gordon Baym and Chloe be probed further and the effects of fermion degeneracy Richards for interesting discussions. This work was sup- must be taken into account. The possible effects of mag- ported in part by National Science Foundation (NSF) netic fields also have been ignored and should be investi- grants No. PHY-1707526 and No. PHY-2010394 to Bow- gated. We implicitly assume that the weak dependence doin College, and NSF grants No. PHY-1662211 and of the accretion rates and times on the polytropic in- No. PHY-2006066 and National Aeronautics and Space dex of the EOS indicate a similarly weak dependence on Administration (NASA) grant No. 80NSSC17K0070 to the detailed properties of realistic nuclear EOSs. There the University of Illinois at Urbana-Champaign.

[1] S. Hawking, Gravitationally collapsed objects of very low [10] B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, mass, Mon. Not. R. Astron. Soc. 152, 75 (1971). Constraints on Primordial Black Holes, [2] B. J. Carr and S. W. Hawking, Black holes in the early arXiv:2002.12778 [astro-ph.CO] (2020). Universe, Mon. Not. R. Astron. Soc. 168, 399 (1974). [11] D. Markovic, Evolution of a primordial [3] G. F. Chapline, Cosmological effects of primordial black black hole inside a rotating solar-type star, holes, Nature (London) 253, 251 (1975). Mon. Not. R. Astron. Soc. 277, 25 (1995). [4] B. Carr and F. K¨uhnel, Primordial Black Holes [12] Y. G´enolini, P. D. Serpico, and P. Tinyakov, as Dark Matter: Recent Developments, An- Revisiting capture into nual Review of Nuclear and Particle Science neutron stars, Phys. Rev. D 102, 083004 (2020), 70, 10.1146/annurev-nucl-050520-125911 (2020), arXiv:2006.16975 [astro-ph.HE]. arXiv:2006.02838 [astro-ph.CO]. [13] P. Montero-Camacho, X. Fang, G. Vasquez, [5] S. W. Hawking, Black hole explosions?, M. Silva, and C. M. Hirata, Revisiting con- Nature (London) 248, 30 (1974). straints on asteroid-mass primordial black holes as [6] F. Capela, M. Pshirkov, and P. Tinyakov, dark matter candidates, JCAP 2019 , 031 (2019), Constraints on primordial black holes as arXiv:1906.05950 [astro-ph.CO]. dark matter candidates from capture by neu- [14] G. M. Fuller, A. Kusenko, and V. Takhistov, tron stars, Phys. Rev. D 87, 123524 (2013), Primordial Black Holes and r -Process Nucle- arXiv:1301.4984 [astro-ph.CO]. osynthesis, Phys. Rev. Lett. 119, 061101 (2017), [7] F. K¨uhnel and K. Freese, Constraints on arXiv:1704.01129 [astro-ph.HE]. primordial black holes with extended mass [15] V. Takhistov, G. M. Fuller, and A. Kusenko, A functions, Phys. Rev. D 95, 083508 (2017), Test for the Origin of Solar Mass Black Holes, arXiv:1701.07223 [astro-ph.CO]. arXiv:2008.12780 [astro-ph.HE] (2020). [8] M. Sasaki, T. Suyama, T. Tanaka, [16] W. E. East and L. Lehner, Fate of a neutron and S. Yokoyama, Primordial black star with an endoparasitic black hole and implica- holes—perspectives in gravitational wave astronomy, tions for dark matter, Phys. Rev. D 100, 124026 (2019), Classical and Quantum Gravity 35, 063001 (2018), arXiv:1909.07968 [gr-qc]. arXiv:1801.05235 [astro-ph.CO]. [17] C. B. Richards, T. W. Baumgarte, and S. L. Shapiro, Ac- [9] V. Vaskonen and H. Veerm¨ae, Did NANOGrav cretion onto a small black hole at the center of a neutron see a signal from primordial black hole for- star, arXiv:2102.09574 [astro-ph.HE] (2021). mation?, Phys. Rev. Lett. 126, 051303 (2021), [18] H. Bondi, On spherically symmetrical accretion, arXiv:2009.07832 [astro-ph.CO]. Mon. Not. R. Astron. Soc. 112, 195 (1952). 5

[19] S. L. Shapiro and S. A. Teukolsky, Black holes, white imum mass, Phys. Rev. D 97, 021501 (2018), dwarfs, and neutron stars : the physics of compact objects arXiv:1711.00473 [astro-ph.HE]. (Wiley Interscience, 1983). [30] B. P. Abbott, LIGO Scientific Collaboration, and [20] The underlying reasons for this result are that (1) the Virgo Collaboration, GW170817: Observation self-gravity of the neutron star gas that is bound to the of Gravitational Waves from a Binary Neutron 2 black hole inside ra ∼ M/a ≪ R can be neglected in Star Inspiral, Phys. Rev. Lett. 119, 161101 (2017), comparison to the gravitational field of the black hole, arXiv:1710.05832 [gr-qc]. and (2) the (spherical) stellar matter exterior to ∼ ra has [31] H. T. Cromartie, E. Fonseca, S. M. Ransom, P. B. De- no influence on the accretion flow by virtue of Birkhoff’s morest, Z. Arzoumanian, H. Blumer, P. R. Brook, M. E. theorem. DeCesar, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. [21] C. Kouvaris and P. Tinyakov, Growth of Ferrara, N. Garver-Daniels, P. A. Gentile, M. L. Jones, black holes in the interior of rotating neu- M. T. Lam, D. R. Lorimer, R. S. Lynch, M. A. McLaugh- tron stars, Phys. Rev. D 90, 043512 (2014), lin, C. Ng, D. J. Nice, T. T. Pennucci, R. Spiewak, I. H. arXiv:1312.3764 [astro-ph.SR]. Stairs, K. Stovall, J. K. Swiggum, and W. W. Zhu, Rela- [22] C. B. Richards, T. W. Baumgarte, and S. L. Shapiro, tivistic Shapiro delay measurements of an extremely mas- Relativistic Bondi accretion for stiff equations of sive , Nature Astronomy 4, 72 (2020), state, Mon. Not. R. Astron. Soc. 502, 3003 (2021), arXiv:1904.06759 [astro-ph.HE]. arXiv:2101.08797 [astro-ph.HE]. [32] I. Goldman and S. Nussinov, Weakly inter- [23] For sufficiently small initial black holes, M/Mtot ≪ 1, acting massive particles and neutron stars, the longest time is spent undergoing quasistatic Bondi Phys. Rev. D 40, 3221 (1989). accretion, while the final phase, when M/Mtot ∼ 1, only [33] A. de Lavallaz and M. Fairbairn, Neutron stars as lasts a much shorter dynamical timescale; see [17]. dark matter probes, Phys. Rev. D 81, 123521 (2010), [24] J. R. Oppenheimer and G. M. Volkoff, On Massive Neu- arXiv:1004.0629 [astro-ph.GA]. tron Cores, Physical Review 55, 374 (1939). [34] J. Bramante and T. Linden, Detecting Dark [25] R. C. Tolman, Static Solutions of Ein- Matter with Imploding in the Galac- stein’s Field Equations for Spheres of Fluid, tic Center, Phys. Rev. Lett. 113, 191301 (2014), Physical Review 55, 364 (1939). arXiv:1405.1031 [astro-ph.HE]. [26] B. Margalit and B. D. Metzger, Con- [35] J. Bramante and F. Elahi, Higgs portals to straining the Maximum Mass of Neutron pulsar collapse, Phys. Rev. D 91, 115001 (2015), Stars from Multi-messenger Observations of arXiv:1504.04019 [hep-ph]. GW170817, Astrophys. J. Lett. 850, L19 (2017), [36] J. Bramante, T. Linden, and Y.-D. Tsai, Search- arXiv:1710.05938 [astro-ph.HE]. ing for dark matter with neutron star mergers [27] M. Shibata, S. Fujibayashi, K. Hotokezaka, K. Ki- and quiet kilonovae, Phys. Rev. D 97, 055016 (2018), uchi, K. Kyutoku, Y. Sekiguchi, and M. Tanaka, arXiv:1706.00001 [hep-ph]. Modeling GW170817 based on numerical relativity [37] M. C. Begelman, Accretion of γ > 5/3 gas by a and its implications, Phys. Rev. D 96, 123012 (2017), Schwarzschild black hole, Astronomy and arXiv:1710.07579 [astro-ph.HE]. 70, 583 (1978). [28] L. Rezzolla, E. R. Most, and L. R. Weih, Using [38] M. J. Rees, Accretion and the phenomenon, Gravitational-wave Observations and Quasi-universal Physica Scripta 17, 193 (1978). Relations to Constrain the Maximum Mass of Neu- [39] K. Inayoshi, Z. Haiman, and J. P. Ostriker, Hyper- tron Stars, Astrophys. J. Lett. 852, L25 (2018), Eddington accretion flows on to massive black arXiv:1711.00314 [astro-ph.HE]. holes, Mon. Not. R. Astron. Soc. 459, 3738 (2016), [29] M. Ruiz, S. L. Shapiro, and A. Tsokaros, arXiv:1511.02116 [astro-ph.HE]. GW170817, general relativistic magnetohydrody- namic simulations, and the neutron star max-