Follow-up of Neutron Star Mergers with CTA and Prospects for Joint Detection with Gravitational-Wave Detectors

Abstract: The joint gravitational wave (GW) and electromagnetic observations of the binary neutron star (BNS) merger GW170817 marked a giant leap in multi-messenger astrophysics. The extensive observation campaign of the associated Gamma-Ray Burst (GRB) and its afterglow has strengthened the hypothesis associating GRBs with BNS mergers and provided insights on mass ejection, particularly the relativistic outflow launched in BNS mergers. In this paper, we investigate the joint detection probabilities of BNS mergers by GW detectors and the upcoming ground-based very-high-energy (VHE) $\gamma$-ray instrument, the Cherenkov Telescope Array (CTA). Using an empirical relation that constrains the distance-inclination angle plane, we simulated BNS mergers detectable in the O5 run of the LIGO/Virgo/Kagra (LVK) network with $300$~Mpc BNS horizon. Assuming Gaussian structured jets and ignoring large sky localization challenges of GW detectors, we estimated VHE afterglow detection probability by CTA. We have explored the afterglow parameter space to identify conditions favourable for detecting synchrotron self-Compton emission by CTA. Our study reveals that events viewed at angles $\lesssim3$ times the jet core angle are detectable by CTA when the initial bulk Lorentz factor at the jet axis ranges between 100 and 800. We find high kinetic energy ($E_k>10^{50}$ erg), ambient density ($n_0>10^{-1}$ $cm^{-3}$), and energy content in non-thermal electrons significantly enhance the likelihood of CTA detection within 300 Mpc. The joint detection rate varies significantly with afterglow parameter distributions, ranging from $0.003$ to $0.5$ per year.