Next Article in Journal
Analog Particle Production Model for General Classes of Taub-NUT Black Holes
Next Article in Special Issue
Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers
Previous Article in Journal
Superconformal Line Defects in 3D
Previous Article in Special Issue
Inclination Estimates from Off-Axis GRB Afterglow Modelling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties

1
Institute of Astronomy, National Tsing Hua University, Hsinchu 300044, Taiwan
2
SUPA, School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, UK
3
School of Physics and Astronomy, University of Leicester, Leicester LE1 7RH, UK
*
Author to whom correspondence should be addressed.
Universe 2021, 7(9), 349; https://doi.org/10.3390/universe7090349
Submission received: 2 September 2021 / Revised: 15 September 2021 / Accepted: 15 September 2021 / Published: 17 September 2021
(This article belongs to the Special Issue Waiting for GODOT—Present and Future of Multi-Messenger Astronomy)

Abstract

:
In the field of multi-messenger astronomy, Bayesian inference is commonly adopted to compare the compatibility of models given the observed data. However, to describe a physical system like neutron star mergers and their associated gamma-ray burst (GRB) events, usually more than ten physical parameters are incorporated in the model. With such a complex model, likelihood evaluation for each Monte Carlo sampling point becomes a massive task and requires a significant amount of computational power. In this work, we perform quick parameter estimation on simulated GRB X-ray light curves using an interpolated physical GRB model. This is achieved by generating a grid of GRB afterglow light curves across the parameter space and replacing the likelihood with a simple interpolation function in the high-dimensional grid that stores all light curves. This framework, compared to the original method, leads to a ∼90× speedup per likelihood estimation. It will allow us to explore different jet models and enable fast model comparison in the future.

1. Introduction

The gravitational wave emitted from a merging binary neutron star, GW170817 [1], followed swiftly by a Gamma-ray burst (GRB) of short duration known as GRB170817A [2,3] commenced the era of multi-messenger astrophysics (MMA). This is the first direct evidence showing neutron star mergers as the progenitor of short GRBs [4]. The union between electromagnetic (EM) and gravitational wave (GW) observations allow for astrophysical objects to be investigated with alternative probes, enabling the source properties to be studied under a new light.
The coalescence of compact objects with at least one neutron star is of particular interest. By the end of the third LIGO observing run (O3), one more binary neutron star (BNS) event, GW190425 [5], and two NS-BH mergers are reported [6]. However, GW170817 remains the only event with a confirmed GW-EM detection and has allowed a major leap in terms of its scientific implications.
A broadband search of the EM counterpart of GW170817 shows both thermal kilonova emission (AT2017gfo) generated from the radioactive decay of newly synthesized elements in the ejecta and a non-thermal synchrotron component from the relativistic jet [7]. The latter is the GRB afterglow that typically shines in X-ray, optical, and radio wavelengths and can last for months to years following the GRB. For the case of GRB170817A, the first confirmed X-ray counterpart was observed by the Chandra X-ray Observatory and announced at around 9 days post-merger [8]. Late-time monitoring of this event unveils the X-ray afterglow emission reaching its peak luminosity at around 155 days [9,10] and continues to shine at more than 1000 days after the GRB [11,12,13].
Modeling the afterglow light curves would shed lights on the physical environment of the system, as the afterglow peak time is dependent on the observers’ viewing angle and the angular profile of the jet. Studies of the temporal behavior of GRB170817A afterglow suggest a structured jet seen off-axis, see, e.g., in [13,14,15,16,17,18]. Understanding of the jet structure is essential to explain the observations of off-axis BNS events such as GRB170817A and uncovering the underlying physical mechanisms behind their production. With better modelling of the jet structure, we are able to more confidently infer the rate of GRB production [19], and the inference of the viewing angle of individual MMA events improves, allowing for tighter constraints on inference of the Hubble constant [20,21].
In this paper, we describe a semi-analytical GRB afterglow model that takes into account the effects of lateral spreading in Section 2. Based on this model, we develop a Bayesian framework that would allow quick parameter estimation and explain it in Section 3. We present the results with this methodology and discuss its future application in Section 4.

2. GRB Afterglow

A broadband afterglow lasting hours-to-days typically follows a GRB, where this GRB afterglow is produced in the shock system that forms as the ultra-relativistic jet that emitted the GRB decelerates. Using energy conservation for a spherical blastwave, the dynamical behaviour of the decelerating system can be modeled, and the change in Lorentz factor, d Γ , with the change in swept-up mass found, see, e.g., in [22]. The instantaneous Γ and swept-up mass of the blastwave can be used to find the synchrotron emission responsible for the broadband GRB afterglow. By following the method in [23], a single jetted outflow with a half opening angle, θ j , is split into multiple emission components and the dynamics and synchrotron flux, relative to the observers line-of-sight from each segment, can be found. Summing across the multiple components at a given observer time, t, is equivalent to integrating across the equal arrival time surface for the jet, and by allowing each component to have a unique energy and Γ , then the afterglow from complex jet structures can be found, see, e.g., in [24].
This method naturally accounts for the edge-effect in a GRB afterglow, where the jet edge becomes visible for an on-axis observer when Γ ( t ) < 1 / θ j , and the resulting jet-break in the lightcurve. We further add the effects of lateral spreading in the jet by considering the maximum sound speed of each component and the change in radius due to this lateral spreading (see in [23,25]). With this method, the afterglow from spreading and non-spreading jets, with a defined structure profile can be found for observers at any viewing angle i.e., within the jet opening angle for cosmological GRBs, see, e.g., in [26], and at higher viewing angles for gravitational wave counterparts to neutron star mergers, see e.g., in [15,16].
The semi-analytic approach used in generating the lightcurves divides the jet emission area over its radial and angular components into a grid of a given resolution. The dynamics and flux are solved for each resolution element within the jet individually before summing to give the resultant lightcurve. The duration of each afterglow iteration depends on the resolution, where for off-axis events the resolution can be set to low values > 50, while for events with low viewing angle, wide opening angles or high Γ , the required resolution becomes much higher, and typically >10,000 to maintain the same level of accuracy. Fitting algorithms can require 100s of thousands of individual model lightcurves for a reliable fit. In our test model with four free parameters, it requires at least ≥1 × 10 4 iterations for a single analysis. A model that takes several seconds, or even minutes, per iteration is therefore not ideal.

Jet Structuring

For a power-law jet, the energy and Lorentz factor distribution of the ejecta with respect to jet central axis is scaled as
y ( θ ) = 1 if 0 θ θ c , θ θ c k if θ c θ θ j , 0 if θ θ j .
In this scenario, the jet is parameterized by two characteristic angles, θ c and θ j . Starting from the center of the jet, the energy is uniformly distributed within the core, θ c . Outside of θ c , the energy distribution follows a power-law decay with a sharp decline at the edge, θ j . This structure was proposed to account for the observed power-law decay in GRB afterglow light curves. Here, we assume the power-law index k = 2 as used in Lamb and Kobayashi [24].

3. Parameter Estimation

The model described in Section 2 is used to simulate X-ray afterglow light curves from 0.15 to 1000 days after the GRB event. Based on the X-ray fluxes, we generate random noise assuming a minimum signal-to-nose ratio (SNR) of 3. The injection data are the sum of both the signal and the noise components. We perform the parameter estimation using Bayesian inference, from which the posterior distribution of jet parameters is given by, according to Bayes theorem,
p ( θ | d ) = L ( d | θ ) π ( θ ) Z ,
where L ( d | θ ) is the likelihood function of the data given a set of model parameters θ , and π ( θ ) is the prior distribution for θ and Z is the evidence, or the observed data. A nested sampling algorithm is used to calculate the evidence and the posterior densities. We carry out the analysis with Bilby, a Python-based Bayesian inference library [27], and the Dynesty sampler [28]. However, one of the issues with stochastic sampling over a complicated physical system is that the likelihood evaluation in the process are computationally expensive.
In order to increase the efficiency in investigating different jet models, we employ an interpolated GRB model and substitute the original likelihood function with it. This is done by first simulating each model parameter across a certain range in the parameter space and storing the resulting light curves in a multidimensional grid. The detailed configuration of the grid is specified below:
  • l o g 10 ϵ 0 ( 50 , 53 )
  • θ j ( 0 , 30 ° )
  • θ c ( 0 , 30 ° )
  • θ o b s ( 0 , 45 ° ),
where θ j is the jet half-opening angle, θ c is the core width of a structured jet, θ o b s is the angle between the observer and the jet central axis, and ϵ 0 is the isotropic equivalent energy of the jet. For each dimension, the parameter range is linearly cut into subsets. We increase the density of splits for parameters that have higher impacts to the resulting light curves. The total number of entries in our 4D Cartesian grid is then ϵ 0 × θ j × θ c × θ o b s = 40 × 50 × 40 × 50 = 4,000,000. Figure 1 shows how afterglow light curves distribute in the three-dimensional space for an observer at different viewing angles. When the inclination angle increases, the afterglow flux peaks at later times. Jet parameters in this dataset are θ j = 0.2 rad (∼ 11.5 ° ), θ c = 0.15 rad (∼ 8.6 ° ) and ϵ 0 = 2 × 10 52 ergs 1 .
To put the emphasis on the jet structure and save computational power, we keep the rest of the afterglow parameters constant throughout the simulation. The fixed parameters include the initial bulk Lorentz factor Γ = 100 , the ambient particle number density n = 0.001 cm 3 , the shock accelerated election index p = 2.15 , the microphysical parameters ϵ B = 0.01 and ϵ e = 0.1 [24], and the luminosity distance d = 40 Mpc [1] for a GW170817-like source. The resolution in the jet model is set to 50 for all light curves.
We adopt uniform priors with upper and lower bounds the same as the simulation range for every parameter of interest. At each Monte Carlo sampling point, instead of calculating the likelihood from the original model, it performs a linear interpolation between the adjacent parameter values in the grid. A comparison between the interpolated and theoretical light curves with same parameter values is shown in Figure 2. The largest difference appears at observing angles near the jet edge, i.e., slightly off-axis case with ~10% deviation. As the geometry between the jet and observers affects the light curves the most, we set the numbers of splits for θ j and θ o b s larger than other parameters. Increasing the density of splits in the grid would raise the accuracy but the simulation time also dramatically grows. For post-peak observing epochs and the rest of off-axis scenarios, the interpolated light curves well resemble the theoretical values with a <10% deviation. As our input data has a SNR value of 3, the error caused by the lack of entries in the interpolation process is not the dominant source of uncertainty.

4. Results and Discussion

The input data, produced with the method described in Section 3, are the X-ray light curve generated from a GRB jet with a power-law profile. The width of the jet is set to θ j = 0.3 rad with the core angle being θ c = 0.15 rad. Considering an off-axis observer, we set the viewing angle θ o b s = 0.36 rad ( 21 ° ) and the isotropic energy is ϵ 0 = 2 × 10 52 ergs 1 . The remaining parameters have the same values as those used to simulate the 4D grid. This light curve is composed of the 0.4 keV fluxes at observing times from 0.16 days to 1000 days post-event and visualized as blue dots in Figure 3. After incorporating the interpolated model in the sampling process, the time required for calculating the likelihood per iteration reduces from 7.78 s to ∼ 0.09 s with only 1 CPU. The overall run time with this method decreases by ∼90 times. Figure 4 shows the joint posterior of the four parameters mentioned above. The injection values are well recovered within the 68% confidence interval. This serves as a proof-of-principle that this method is not only efficient, but also capable of capturing the jet features if the afterglow emission is intrinsically produced by such a system.
In this work, we perform an analysis on simulated light curves with a uniform signal-to-noise ratio throughout all observing epochs. In reality, for the same patch of sky and the same observational instrument, the background fluence would remain at a similar scale. The temporal evolution of a GRB afterglow, which arises until reaching the peak luminosity at t p e a k and then declines with time, yields a varying observational significance. The X-ray emission from the position of GW170817, except for the first observation at 2.3 days, reports a significance of ≥ 3 σ at all epochs. We set all SNR equaling to the detection limit for a conservative estimate. Recent studies show that the X-ray counterpart of GW170817 is still bright after 1200 days [11,12]. The 3 σ excess compared to the prediction from an off-axis structured jet poses a challenge to previous agreements from multi-wavelength studies. It also makes modeling the afterglow of GRB170817 worthwhile even after four years as new evidence is still emerging.
On the other hand, we use relatively wide uniform priors for jet structure parameters in the analysis. One of the strengths of multi-messenger astronomy is that we can combine the information obtained from various approaches and place tighter constraints on the source properties. As gravitational wave observation provides an independent measure of viewing angle and luminosity distance, the posteriors from GW parameter estimation can be passed to the afterglow analysis as prior function and thus increase the accuracy in search of the GRB afterglow parameters.
In this paper, we present the validity of a Bayesian framework combined with an interpolated GRB model. It effectively lessens the computational cost while being capable of constraining jet structure parameters in low SNR circumstances. In the upcoming era of a global gravitational wave detector network (LIGO-VIRGO-KAGRA) and new EM observatories like the Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor (GECAM), it is foreseeable that we will find more GW-EM counterparts in the near future. A fast, functional parameter estimation will be beneficial for modeling the GRB afterglow light curves. With this methodology, more jet models, e.g., Gaussian or Top-hat jets, will be investigated and its application on the real data of GW170817, the only event so far with GW-EM detection, will be included in our future works.

Author Contributions

Conceptualization, E.-T.L., F.H. and I.S.H.; Formal analysis/Investigation/ visualization, E.-T.L.; Methodology, E.-T.L., F.H. and G.P.L.; Resources, M.J.W.; Writing—original draft, E.-T.L., F.H. and G.P.L.; Writing—review and editing, A.K.H.K. and S.S. Supervision, J.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Ministry of Science and Technology of Taiwan through grant 109-2628-M-007-005-RSP.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abbott, B.P.; Birney, R.; Jawahar, S.; Lockerbie, N.A.; Reid, S.; Tokmakov, K.V. GW170817: Observation of gravitational waves from a binary neutron star inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [Green Version]
  2. Goldstein, A.; Veres, P.; Burns, E.; Briggs, M.; Hamburg, R.; Kocevski, D.; Wilson-Hodge, C.; Preece, R.; Poolakkil, S.; Roberts, O.; et al. An ordinary short gamma-ray burst with extraordinary implications: Fermi-GBM detection of GRB 170817A. ApJL 2017, 848, L14. [Google Scholar] [CrossRef] [Green Version]
  3. Savchenko, V.; Ferrigno, C.; Kuulkers, E.; Bazzano, A.; Bozzo, E.; Brandt, S.; Chenevez, J.; Courvoisier, T.L.; Diehl, R.; Domingo, A.; et al. INTEGRAL detection of the first prompt gamma-ray signal coincident with the gravitational-wave event GW170817. ApJL 2017, 848, L15. [Google Scholar] [CrossRef] [Green Version]
  4. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. Gravitational Waves and Gamma-rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A. ApJL 2017, 848, L13. [Google Scholar] [CrossRef]
  5. Abbott, B.; Abbott, R.; Abbott, T.; Abraham, S.; Acernese, F.; Ackley, K.; Adams, C.; Adhikari, R.; Adya, V.; Affeldt, C.; et al. GW190425: Observation of a compact binary coalescence with total mass 3.4 M⊙. ApJL 2020, 892, L3. [Google Scholar] [CrossRef]
  6. Abbott, R.; Abbott, T.D.; Abraham, S.; Acernese, F.; Ackley, K.; Adams, A.; Adams, C.; Adhikari, R.X.; Adya, V.B.; Affeldt, C.; et al. Observation of gravitational waves from two neutron star–black hole coalescences. ApJL 2021, 915, L5. [Google Scholar] [CrossRef]
  7. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. Multi-messenger observations of a binary neutron star merger. ApJL 2017, 848, L12. [Google Scholar] [CrossRef]
  8. Troja, E.; Piro, L.; van Eerten, H.; Wollaeger, R.T.; Im, M.; Fox, O.D.; Butler, N.R.; Cenko, S.B.; Sakamoto, T.; Fryer, C.L.; et al. The X-ray counterpart to the gravitational-wave event GW170817. Nature 2017, 551, 71–74. [Google Scholar] [CrossRef] [Green Version]
  9. Dobie, D.; Kaplan, D.L.; Murphy, T.; Lenc, E.; Mooley, K.P.; Lynch, C.; Corsi, A.; Frail, D.; Kasliwal, M.; Hallinan, G. A Turnover in the Radio Light Curve of GW170817. ApJL 2018, 858, L15. [Google Scholar] [CrossRef]
  10. Lin, E.T.; Yu, H.F.; Kong, A.K. Bayesian analysis on the X-ray spectra of the binary neutron star merger GW170817. JHEAP 2019, 21, 1–5. [Google Scholar] [CrossRef]
  11. Hajela, A.; Margutti, R.; Bright, J.S.; Alexander, K.D.; Metzger, B.D.; Nedora, V.; Kathirgamaraju, A.; Margalit, B.; Radice, D.; Berger, E.; et al. The emergence of a new source of X-rays from the binary neutron star merger GW170817. arXiv 2021, arXiv:2104.02070. [Google Scholar]
  12. Troja, E.; O’Connor, B.; Ryan, G.; Piro, L.; Ricci, R.; Zhang, B.; Piran, T.; Bruni, G.; Cenko, S.B.; van Eerten, H. Accurate flux calibration of GW170817: Is the X-ray counterpart on the rise? arXiv 2021, arXiv:2104.13378. [Google Scholar]
  13. Troja, E.; van Eerten, H.; Zhang, B.; Ryan, G.; Piro, L.; Ricci, R.; O’Connor, B.; Wieringa, M.H.; Cenko, S.B.; Sakamoto, T. A thousand days after the merger: Continued X-ray emission from GW170817. MNRAS 2020, 498, 5643–5651. [Google Scholar] [CrossRef]
  14. Ghirlanda, G.; Salafia, O.S.; Paragi, Z.; Giroletti, M.; Yang, J.; Marcote, B.; Blanchard, J.; Agudo, I.; An, T.; Bernardini, M.G.; et al. Compact radio emission indicates a structured jet was produced by a binary neutron star merger. Science 2019, 363, 968–971. [Google Scholar] [CrossRef] [Green Version]
  15. Lyman, J.D.; Lamb, G.P.; Levan, A.J.; Mandel, I.; Tanvir, N.R.; Kobayashi, S.; Gompertz, B.; Hjorth, J.; Fruchter, A.S.; Kangas, T.; et al. The optical afterglow of the short gamma-ray burst associated with GW170817. Nat. Astron. 2018, 2, 751–754. [Google Scholar] [CrossRef]
  16. Lamb, G.P.; Lyman, J.D.; Levan, A.J.; Tanvir, N.R.; Kangas, T.; Fruchter, A.S.; Gompertz, B.; Hjorth, J.; Mandel, I.; Oates, S.R.; et al. The Optical Afterglow of GW170817 at One Year Post-merger. ApJL 2019, 870, L15. [Google Scholar] [CrossRef] [Green Version]
  17. Resmi, L.; Schulze, S.; Ishwara-Chandra, C.H.; Misra, K.; Buchner, J.; De Pasquale, M.; Sánchez-Ramírez, R.; Klose, S.; Kim, S.; Tanvir, N.R.; et al. Low-frequency View of GW170817/GRB 170817A with the Giant Metrewave Radio Telescope. ApJL 2018, 867, 57. [Google Scholar] [CrossRef] [Green Version]
  18. Troja, E.; Piro, L.; Ryan, G.; van Eerten, H.; Ricci, R.; Wieringa, M.H.; Lotti, S.; Sakamoto, T.; Cenko, S.B. The outflow structure of GW170817 from late-time broad-band observations. MNRAS 2018, 478, L18–L23. [Google Scholar] [CrossRef]
  19. Tan, W.W.; Yu, Y.W. The jet structure and the intrinsic luminosity function of short gamma-ray bursts. ApJ 2020, 902, 83. [Google Scholar] [CrossRef]
  20. Wang, H.; Giannios, D. Multimessenger Parameter Estimation of GW170817: From Jet Structure to the Hubble Constant. ApJ 2021, 908, 200. [Google Scholar] [CrossRef]
  21. Nakar, E.; Piran, T. Afterglow constraints on the viewing angle of binary neutron star mergers and determination of the Hubble constant. ApJ 2021, 909, 114. [Google Scholar] [CrossRef]
  22. Peer, A. Dynamical Model of an Expanding Shell. ApJL 2012, 752, L8. [Google Scholar] [CrossRef]
  23. Lamb, G.P.; Mandel, I.; Resmi, L. Late-time evolution of afterglows from off-axis neutron star mergers. MNRAS 2018, 481, 2581–2589. [Google Scholar] [CrossRef] [Green Version]
  24. Lamb, G.P.; Kobayashi, S. Electromagnetic counterparts to structured jets from gravitational wave detected mergers. MNRAS 2017, 472, 4953–4964. [Google Scholar] [CrossRef]
  25. Lamb, G.P.; Kann, D.A.; Fernández, J.J.; Mandel, I.; Levan, A.J.; Tanvir, N.R. GRB jet structure and the jet break. arXiv 2021, arXiv:2104.11099. [Google Scholar]
  26. Lamb, G.P.; Tanvir, N.R.; Levan, A.J.; de Ugarte Postigo, A.; Kawaguchi, K.; Corsi, A.; Evans, P.A.; Gompertz, B.; Malesani, D.B.; Page, K.L.; et al. Short GRB 160821B: A Reverse Shock, a Refreshed Shock, and a Well-sampled Kilonova. ApJ 2019, 883, 48. [Google Scholar] [CrossRef]
  27. Ashton, G.; Hübner, M.; Lasky, P.D.; Talbot, C.; Ackley, K.; Biscoveanu, S.; Chu, Q.; Divakarla, A.; Easter, P.J.; Goncharov, B.; et al. BILBY: A User-friendly Bayesian Inference Library for Gravitational-wave Astronomy. ApJs 2019, 241, 27. [Google Scholar] [CrossRef] [Green Version]
  28. Speagle, J.S. DYNESTY: A dynamic nested sampling package for estimating Bayesian posteriors and evidences. MNRAS 2020, 493, 3132–3158. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Three-dimensional illustration of afterglow light curves of a power-law jet seen from θ o b s = 0 ° to θ o b s = 45 ° . Gray dots are data stored in our 4D grid. The open-angle of the jet is θ j = 0.2 (rad) and the core width θ c = 0.15 (rad). Other fixed parameters include ϵ 0 = 2 × 10 52 ergs 1 , n = 0.001 cm 3 , p = 2.15 , ϵ B = 0.01 , ϵ e = 0.1 and d = 40 Mpc.
Figure 1. Three-dimensional illustration of afterglow light curves of a power-law jet seen from θ o b s = 0 ° to θ o b s = 45 ° . Gray dots are data stored in our 4D grid. The open-angle of the jet is θ j = 0.2 (rad) and the core width θ c = 0.15 (rad). Other fixed parameters include ϵ 0 = 2 × 10 52 ergs 1 , n = 0.001 cm 3 , p = 2.15 , ϵ B = 0.01 , ϵ e = 0.1 and d = 40 Mpc.
Universe 07 00349 g001
Figure 2. A comparison between the interpolated (dark blue) and simulated light curves (red) at observation angles of θ o b s = 20 ° ,   21 ° and 22 ° . These observing angles lie in the midpoints of the 4D-grid of simulated light curves (light blue) used to perform the interpolation. The inset shows the fractional difference between the interpolated light curves and the simulated light curves. Note that all parameter values are same as Figure 1 except for a slightly wider jet with θ j = 0.3 rad (∼ 17 ° ).
Figure 2. A comparison between the interpolated (dark blue) and simulated light curves (red) at observation angles of θ o b s = 20 ° ,   21 ° and 22 ° . These observing angles lie in the midpoints of the 4D-grid of simulated light curves (light blue) used to perform the interpolation. The inset shows the fractional difference between the interpolated light curves and the simulated light curves. Note that all parameter values are same as Figure 1 except for a slightly wider jet with θ j = 0.3 rad (∼ 17 ° ).
Universe 07 00349 g002
Figure 3. Injection data points (blue) overlaid with the best-fit power-law jet model (black solid line). Red lines are 1000 random draws from the posterior samples.
Figure 3. Injection data points (blue) overlaid with the best-fit power-law jet model (black solid line). Red lines are 1000 random draws from the posterior samples.
Universe 07 00349 g003
Figure 4. Probability density distribution of parameters in a power-law structured jet model. Orange vertical lines indicate the injected values from which the input data was simulated.
Figure 4. Probability density distribution of parameters in a power-law structured jet model. Orange vertical lines indicate the injected values from which the input data was simulated.
Universe 07 00349 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, E.-T.; Hayes, F.; Lamb, G.P.; Heng, I.S.; Kong, A.K.H.; Williams, M.J.; Saha, S.; Veitch, J. A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties. Universe 2021, 7, 349. https://doi.org/10.3390/universe7090349

AMA Style

Lin E-T, Hayes F, Lamb GP, Heng IS, Kong AKH, Williams MJ, Saha S, Veitch J. A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties. Universe. 2021; 7(9):349. https://doi.org/10.3390/universe7090349

Chicago/Turabian Style

Lin, En-Tzu, Fergus Hayes, Gavin P. Lamb, Ik Siong Heng, Albert K. H. Kong, Michael J. Williams, Surojit Saha, and John Veitch. 2021. "A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties" Universe 7, no. 9: 349. https://doi.org/10.3390/universe7090349

APA Style

Lin, E. -T., Hayes, F., Lamb, G. P., Heng, I. S., Kong, A. K. H., Williams, M. J., Saha, S., & Veitch, J. (2021). A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties. Universe, 7(9), 349. https://doi.org/10.3390/universe7090349

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop