Testing Bayesian inference of GRMHD model parameters from VLBI data

A. I. Yfantis,S. Zhao,R. Gold,M. Mościbrodzka,A. E. Broderick

2024-09-24

Abstract:Recent observations by the Event Horizon Telescope (EHT) of supermassive black holes M87* and Sgr A* offer valuable insights into their spacetime properties and astrophysical conditions. Utilizing a library of model images (~2 million for Sgr A*) generated from general-relativistic magnetohydrodynamic (GRMHD) simulations, limited and coarse insights on key parameters such as black hole spin, magnetic flux, inclination angle, and electron temperature were gained. The image orientation and black hole mass estimates were obtained via a scoring and an approximate rescaling procedure. Lifting such approximations, probing the space of parameters continuously, and extending the parameter space of theoretical models is both desirable and computationally prohibitive with existing methods. To address this, we introduce a new Bayesian scheme that adaptively explores the parameter space of ray-traced, GRMHD models. The general relativistic radiative transfer code \ipole is integrated with the EHT parameter estimation tool THEMIS. The pipeline produces a ray-traced model image from GRMHD data, computes predictions for VLBI observables from the image for a specific VLBI array configuration and compares to data thereby sampling the likelihood surface via an MCMC scheme. At this stage we focus on four parameters: accretion rate, electron thermodynamics, inclination, and source position angle. Our scheme faithfully recovers parameters from simulated VLBI data and accommodates time-variabibility via an inflated error budget. We highlight the impact of intrinsic variability on model fitting approaches. This work facilitates more informed inferences from GRMHD simulations and enables expansion of the model parameter space in a statistically robust and computationally efficient manner.

High Energy Astrophysical Phenomena,Instrumentation and Methods for Astrophysics

What problem does this paper attempt to address?

The paper aims to address the following issues:

1. **Limitations of existing methods**: Current methods for inferring General Relativistic Magnetohydrodynamics (GRMHD) model parameters using Event Horizon Telescope (EHT) data have some limitations, such as sparse parameter sampling, fixed partial parameters, and high computational costs.

2. **Exploration of continuous parameter space**: Existing methods find it difficult to explore the continuous parameter space in detail, especially when dealing with the expansion of the model parameter space, which appears computationally inefficient.

3. **Variability challenges**: The variability (spatiotemporal changes and structural changes) in GRMHD simulations makes it difficult for traditional methods to accurately match observational data.

To address these issues, the authors propose a new Bayesian inference scheme that integrates the General Relativistic Radiative Transfer code IPOLE with the EHT parameter estimation tool THEMIS, enabling adaptive exploration of GRMHD model parameters. This method not only improves computational efficiency but also better handles model variability and statistically robustly expands the model parameter space.

Specifically, the focus of the paper is on developing a new method that can efficiently and accurately infer GRMHD model parameters from VLBI data, and the effectiveness and accuracy of this new method are validated through synthetic data.