No Arabic abstract
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 present radiative transfer simulations for blue kilonovae hours after neutron star (NS) mergers by performing detailed opacity calculations for the first time. We calculate atomic structures and opacities of highly ionized elements (up to the tenth ionization) with atomic number Z = 20 - 56. We find that the bound-bound transitions of heavy elements are the dominant source of the opacities in the early phase (t < 1 day after the merger), and that the ions with a half-closed electron shell provide the highest contributions. The Planck mean opacity for lanthanide-free ejecta (with electron fraction of Ye = 0.30 - 0.40) can only reach around kappa ~ 0.5 - 1 cm^2 g^-1 at t = 0.1 day, whereas that increases up to kappa ~ 5 - 10 cm^2 g^-1 at t = 1 day. The spherical ejecta model with an ejecta mass of Mej = 0.05Msun gives the bolometric luminosity of ~ 2 x 10^42 erg s^-1 at t ~ 0.1 day. We confirm that the existing bolometric and multi-color data of GW170817 can be naturally explained by the purely radioactive model. The expected early UV signals reach 20.5 mag at t ~ 4.3 hours for sources even at 200 Mpc, which is detectable by the facilities such as Swift and the Ultraviolet Transient Astronomy Satellite (ULTRASAT). The early-phase luminosity is sensitive to the structure of the outer ejecta, as also pointed out by Kasen et al. (2017). Therefore, the early UV observations give strong constraints on the structure of the outer ejecta as well as the presence of a heating source besides r-process nuclei.
We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgos third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)
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 report here the non-detection of gravitational waves from the merger of binary neutron star systems and neutron-star--black-hole systems during the first observing run of Advanced LIGO. In particular we searched for gravitational wave signals from binary neutron star systems with component masses $in [1,3] M_{odot}$ and component dimensionless spins $< 0.05$. We also searched for neutron-star--black-hole systems with the same neutron star parameters, black hole mass $in [2,99] M_{odot}$ and no restriction on the black hole spin magnitude. We assess the sensitivity of the two LIGO detectors to these systems, and find that they could have detected the merger of binary neutron star systems with component mass distributions of $1.35pm0.13 M_{odot}$ at a volume-weighted average distance of $sim$ 70Mpc, and for neutron-star--black-hole systems with neutron star masses of $1.4M_odot$ and black hole masses of at least $5M_odot$, a volume-weighted average distance of at least $sim$ 110Mpc. From this we constrain with 90% confidence the merger rate to be less than 12,600 Gpc$^{-3}$yr$^{-1}$ for binary-neutron star systems and less than 3,600 Gpc$^{-3}$yr$^{-1}$ for neutron-star--black-hole systems. We find that if no detection of neutron-star binary mergers is made in the next two Advanced LIGO and Advanced Virgo observing runs we would place significant constraints on the merger rates. Finally, assuming a rate of $10^{+20}_{-7}$Gpc$^{-3}$yr$^{-1}$ short gamma ray bursts beamed towards the Earth and assuming that all short gamma-ray bursts have binary-neutron-star (neutron-star--black-hole) progenitors we can use our 90% confidence rate upper limits to constrain the beaming angle of the gamma-ray burst to be greater than ${2.3^{+1.7}_{-1.1}}^{circ}$ (${4.3^{+3.1}_{-1.9}}^{circ}$).