Constraining the Outflow Structure of the Binary Neutron Star Merger Event GW170817 with Markov-Chain Monte Carlo Analysis

TL;DR

This study uses a detailed MCMC analysis with hydrodynamic simulations to constrain the outflow structure of GW170817, favoring a relativistic structured jet over other models.

Contribution

It introduces a physically motivated two-parameter 'boosted fireball' model and applies an efficient MCMC fitting method to observational data, advancing understanding of BNS merger outflows.

Findings

Relativistic structured jet model is favored with a ~5° jet opening angle.

Off-axis viewing angle estimated at ~27°.

Parameter degeneracies suggest a wide range of explosion energies and densities.

Abstract

The multi-wavelength non-thermal emission from the binary neutron star (BNS) merger GW170817 has raised a heated debate concerning the post-merger outflow structure. Both a relativistic structured jet viewed off-axis and a mildly relativistic quasi-spherical outflow can explain the observational data of GW170817 up to ~260 days. We utilize a physically motivated analytic two-parameter model called the "boosted fireball" for the outflow structure after it has expanded far from the merger site and has entered the self-similar coasting phase. This model consists of a family of outflows with a structure varying smoothly between a highly collimated ultra-relativistic jet and an isotropic fireball. We simulate the dynamical evolution, starting with "boosted fireball" initial conditions, of 240 outflows using the moving-mesh relativistic hydrodynamics code JET to follow their evolution through the afterglow phase. We compute nearly 2,000,000 synchrotron spectra from the hydrodynamic simulations using the standard synchrotron radiation model. By making use of scaling relations in the hydrodynamic and radiation equations, we develop a synthetic light curve generator with an efficient sampling speed. This allows us to fit the observational data by performing Markov-Chain Monte Carlo (MCMC) analysis in a 8-dimensional parameter space, consisting of hydrodynamic parameters, radiation parameters and observational parameters. Our results favor the relativistic structured jet, with a jet opening angle ~5 deg and Lorentz factor ~175, viewed from an off-axis angle of 27(+9-3) deg. Due to parameter degeneracies, we find broad distributions for the explosion energy E_0, the circumburst density n_0, epsilon_e and epsilon_B. The combination of a high n_0 and a low epsilon_B can also produce a good fit, indicating that an extremely low n_0 may not be required for GW170817.