Jay V. Kalinani, Riccardo Ciolfi, Manuela Campanelli, Bruno Giacomazzo, Andrea Pavan, Allen Wen, Yosef Zlochower

We present general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers, exploring the conditions to launch relativistic outflows compatible with gamma-ray burst (GRB) jets. Employing an unprecedentedly low numerical density floor decreasing as the sixth power of radial distance, we obtain for the first time an evolution where such a floor has no effect on any outflow produced during and after merger. In models forming an accreting black hole (BH) 25 ms after merger, incipient jets are launched with terminal Lorentz factors and Poynting-flux luminosities compatible with GRBs. The rather high jet velocities reached in a few tens of ms, when compared to previous simulations adopting a much higher and uniform density floor, show the importance of removing any influence of the latter as the jet evolves toward increasingly large time and spatial scales. In absence of collapse, we find the emergence of much heavier and slower polar outflows, incapable of producing a GRB. We thus favor a BH origin for GRB jets.