No Arabic abstract
We analyse the tidal disruption probability of potential neutron star--black hole (NSBH) merger gravitational wave (GW) events, including GW190426_152155, GW190814, GW200105_162426 and GW200115_042309, detected during the third observing run of the LIGO/Virgo Collaboration, and the detectability of kilonova emission in connection with these events. The posterior distributions of GW190814 and GW200105_162426 show that they must be plunging events and hence no kilonova signal is expected from these events. With the stiffest NS equation of state allowed by the constraint of GW170817 taken into account, the probability that GW190426_152155 and GW200115_042309 can make tidal disruption is $sim24%$ and $sim3%$, respectively. However, the predicted kilonova brightness is too faint to be detected for present follow-up search campaigns, which explains the lack of electromagnetic (EM) counterpart detection after triggers of these GW events. Based on the best constrained population synthesis simulation results, we find that disrupted events account for only $lesssim20%$ of cosmological NSBH mergers since most of the primary BHs could have low spins. The associated kilonovae for those disrupted events are still difficult to be discovered by LSST after GW triggers in the future, because of their low brightness and larger distances. For future GW-triggered multi-messenger observations, potential short-duration gamma-ray bursts and afterglows are more probable EM counterparts of NSBH GW events.
We present detailed simulations of black hole-neutron star (BH-NS) mergers kilonova and gamma-ray burst (GRB) afterglow and kilonova luminosity function, and discuss the detectability of electromagnetic (EM) counterpart in connection with gravitational wave (GW) detections, GW-triggered target-of-opportunity observations, and time-domain blind searches. The predicted absolute magnitude of the BH-NS kilonovae at $0.5,{rm days}$ after the merger falls in $[-10,-15.5]$. The simulated luminosity function contains the potential viewing-angle distribution information of the anisotropic kilonova emission. We simulate the GW detection rates, detectable distances and signal duration, for the future networks of 2nd/2.5th/3rd-generation GW detectors. BH-NSs tend to produce brighter kilonovae and afterglows if the BH has a higher aligned-spin, and a less massive NS with a stiffer EoS. The detectability of kilonova is especially sensitive to the BH spin. If BHs typically have low spins, the BH-NS EM counterparts are hard to discover. For the 2nd generation GW detector networks, a limiting magnitude of $m_{rm limit}sim23-24,{rm mag}$ is required to detect the kilonovae even if BH high spin is assumed. Thus, a plausible explanation for the lack of BH-NS associated kilonova detection during LIGO/Virgo O3 is that either there is no EM counterpart (plunging events), or the current follow-ups are too shallow. These observations still have the chance to detect the on-axis jet afterglow associated with an sGRB or an orphan afterglow. Follow-up observations can detect possible associated sGRB afterglows, from which kilonova signatures may be studied. For time-domain observations, a high-cadence search in redder filters is recommended to detect more BH-NS associated kilonovae and afterglows.
We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted before the outflow is launched, which depends on the disk position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We estimate the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for GW190814.
We investigate the orbital dynamics of hierarchical three-body systems containing a double neutron star system orbiting around a massive black hole. These systems show complex dynamical behaviour because of relativistic coupling between orbits of the neutron stars in the double neutron star system and the orbit of the double neutron star system around the black hole. The orbital motion of the neutron stars around each other drives a loop mass current, which gives rise to gravito-magnetism. Generally, gravito-magnetism involves a rotating black hole. The hierarchical three-body system that we consider is an unusual situation in which black hole rotation is not required. Using a gravito-electromagnetic formulation, we calculate the orbital precession and nutation of the double neutron star system. These precession and nutation effects are observable, thus providing probes to the spacetime around black holes as well as tests of gravito-electromagnetism in the framework of general relativity.
Recent detailed 1D core-collapse simulations have brought new insights on the final fate of massive stars, which are in contrast to commonly used parametric prescriptions. In this work, we explore the implications of these results to the formation of coalescing black-hole (BH) - neutron-star (NS) binaries, such as the candidate event GW190426_152155 reported in GWTC-2. Furthermore, we investigate the effects of natal kicks and the NSs radius on the synthesis of such systems and potential electromagnetic counterparts linked to them. Synthetic models based on detailed core-collapse simulations result in an increased merger detection rate of BH-NS systems ($sim 2.3$ yr$^{-1}$), 5 to 10 times larger than the predictions of standard parametric prescriptions. This is primarily due to the formation of low-mass BH via direct collapse, and hence no natal kicks, favored by the detailed simulations. The fraction of observed systems that will produce an electromagnetic counterpart, with the detailed supernova engine, ranges from $2$-$25$%, depending on uncertainties in the NS equation of state. Notably, in most merging systems with electromagnetic counterparts, the NS is the first-born compact object, as long as the NSs radius is $lesssim 12,mathrm{km}$. Furthermore, core-collapse models that predict the formation of low-mass BHs with negligible natal kicks increase the detection rate of GW190426_152155-like events to $sim 0.6 , $yr$^{-1}$; with an associated probability of electromagnetic counterpart $leq 10$% for all supernova engines. However, increasing the production of direct-collapse low-mass BHs also increases the synthesis of binary BHs, over-predicting their measured local merger density rate. In all cases, models based on detailed core-collapse simulation predict a ratio of BH-NSs to binary BHs merger rate density that is at least twice as high as other prescriptions.
We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public MOSFiT package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly ($B>10^{10}$) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is $q=0.92pm0.07$ (90% credible interval). We place tight constraints on the maximum stable NS mass, $M_{rm TOV}=2.17^{+0.08}_{-0.11}$ M$_odot$. For a uniform prior in tidal deformability, the radius of a 1.4 M$_odot$ NS is $R_{1.4}sim 10.7$ km. Re-weighting with a prior based on equations of state that support our credible range in $M_{rm TOV}$, we derive a final measurement $R_{1.4}=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally-detected neutron star merger, GW190425, we estimate that an associated kilonova would have been fainter (by $sim0.7$ mag at one day post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.