Towards an understanding of long gamma-ray burst environments through circumstellar medium population synthesis predictions

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock, separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a ‘CSM population synthesis’ with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind, and ISM densities are compared with the largest sample of afterglow derived parameters yet compiled, which we make available for the community. We find that high ISM densities of n ∼ 1000 cm−3 best reproduce observations. If long GRBs instead occur in typical ISM densities of n ∼ 1 cm−3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius. © 2022 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society.
BPG and MN are supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 948381). DAK acknowledges support from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). AJvM is supported by the ANR-19-CE31-0014GAMALO project. AJL has received funding from the European Research Council (ERC) under the European Union’s Seventh Framework Programme (FP7-2007-2013) (Grant agreement No. 725246). ERS has been supported by STFC consolidated grant ST/P000495/1. PJG is supported by NRF SARChI Grant 111692. We gratefully acknowledge the use of GOTOHEAD, the computing cluster of the Gravitational-wave Optical Transient Observer (GOTO), as well as support from Joe Lyman and Krzysztof Ulaczyk.
With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.