Hao Wang,Ranadeep G. Dastidar,Dimitrios Giannios,Paul C. Duffell
Abstract:Gamma-ray burst (GRB) afterglows are emissions from ultrarelativistic blast waves produced by a narrow jet interacting with surrounding matter. Since the first multimessenger observation of a neutron star merger, hydrodynamic modeling of GRB afterglows for structured jets with smoothly varying angular energy distributions has gained increased interest. While the evolution of a jet is well described by self-similar solutions in both ultrarelativistic and Newtonian limits, modeling the transitional phase remains challenging. This is due to the nonlinear spreading of a narrow jet to a spherical configuration and the breakdown of self-similar solutions. Analytical models are limited in capturing these nonlinear effects, while relativistic hydrodynamic simulations are computationally expensive, which restricts the exploration of various initial conditions. In this work, we introduce a reduced hydrodynamic model that approximates the blast wave as an infinitely thin two-dimensional surface. Further assuming axial symmetry, this model simplifies the simulation to one dimension and drastically reduces the computational costs. We have compared our modeling to relativistic hydrodynamic simulations and semianalytic methods, and applied it to fit the light curve and flux centroid motion of GRB 170817A. These comparisons demonstrate good agreement and validate our approach. We have developed this method into a numerical tool, \texttt{jetsimpy}, which models the synchrotron GRB afterglow emission from a blast wave with arbitrary angular energy and Lorentz factor distribution. Although the code is built with GRB afterglow in mind, it applies to any relativistic jet. This tool is particularly useful in Markov Chain Monte Carlo studies and is provided to the community.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve several key problems in gamma - ray burst (GRB) afterglow simulations:
1. **Modeling of non - linear effects**:
- **Problem description**: In the simulation of GRB afterglows, the non - linear expansion of a narrow jet to a spherical configuration is a complex process. Existing analytical models have difficulty capturing these non - linear effects, and while relativistic hydrodynamical simulations are accurate, they are computationally expensive, limiting the exploration of different initial conditions.
- **Solution**: The authors propose a simplified hydrodynamical model. By approximating the shock region as an infinitely thin two - dimensional surface and assuming axial symmetry, the simulation is simplified to one - dimensional, thereby greatly reducing the computational cost.
2. **Simulation of arbitrary angular energy distributions**:
- **Problem description**: Existing GRB afterglow models usually assume that the energy distribution of the jet is a "top - hat" model, that is, the energy is evenly distributed within a narrow core and drops sharply at the edges. This model works well when the observation direction is within the jet core, but when observing off - axis, a smoothly varying angular energy distribution needs to be considered.
- **Solution**: The authors develop a numerical tool `jetsimpy` that can simulate jets with arbitrary angular energy distributions and Lorentz factor distributions and generate synthetic afterglow observables. This makes the tool applicable not only to GRB afterglows but also to the study of other relativistic jets.
3. **Computational efficiency and accuracy**:
- **Problem description**: Although the fully numerical hydrodynamical method is accurate, due to its high computational cost, its application in parameter estimation is limited, especially when using the Markov Chain Monte Carlo (MCMC) method.
- **Solution**: Through the above - mentioned simplified model, `jetsimpy` reduces the computational cost by thousands of times while maintaining an accuracy comparable to that of full - numerical simulations. It is especially suitable for situations requiring a large amount of computation, such as MCMC studies.
### Main contributions of the paper
1. **Proposing an efficient hydrodynamical model**: By approximating the shock region as an infinitely thin two - dimensional surface and assuming axial symmetry, the simulation is simplified to one - dimensional, significantly reducing the computational cost.
2. **Developing the numerical tool `jetsimpy`**: This tool can simulate jets with arbitrary angular energy distributions and Lorentz factor distributions and generate synthetic afterglow observables, applicable to the study of various relativistic jets.
3. **Verifying the accuracy of the model**: By comparing with relativistic hydrodynamical simulations and semi - analytical methods, the validity of the model is verified, especially in cases where non - linear hydrodynamical effects are significant.
### Application prospects
This model and tool have important application prospects in the following aspects:
1. **Multi - messenger astronomy**: GRB afterglows can be used as probes to detect the expansion rate of the universe, helping to solve the Hubble constant tension problem.
2. **Standard siren research**: By constraining the observation angle, it helps to break the degeneracy of distance in gravitational wave data and provides more accurate distance measurements.
3. **MCMC studies**: Greatly reducing the computational cost makes large - scale parameter estimation possible, especially when dealing with complex jet structures.
In conclusion, this paper solves multiple key problems in GRB afterglow simulations by proposing an efficient and accurate hydrodynamical model, providing powerful tools and support for research in related fields.