Hierarchical Bayesian Inference of Globular Cluster Properties

Abstract

We present a hierarchical Bayesian inference approach to estimating the structural properties and the phase space center of a globular cluster (GC) given the spatial and kinematic information of its stars based on lowered isothermal cluster models. As a first step towards more realistic modelling of GCs, we built a differentiable, accurate emulator of the lowered isothermal distribution function using interpolation. The reliable gradient information provided by the emulator allows the use of Hamiltonian Monte Carlo methods to sample large Bayesian models with hundreds of parameters, thereby enabling inference on hierarchical models. We explore the use of hierarchical Bayesian modelling to address several issues encountered in observations of GC including an unknown GC center, incomplete data, and measurement errors. Our approach not only avoids the common technique of radial binning but also incorporates the aforementioned uncertainties in a robust and statistically consistent way. Through demonstrating the reliability of our hierarchical Bayesian model on simulations, our work lays out the foundation for more realistic and complex modelling of real GC data.

16 pages, 12 figures, and 2 tables

References in corpus (11)

• Hubble Space Telescope proper motion (HSTPROMO) catalogs of Galactic globular clusters. II. Kinematic profiles and maps
• The eye of Gaia on globular clusters kinematics: internal rotation
• Overview of the LAMOST survey in the first decade
• Globular cluster number density profiles using Gaia DR2
• Mass models of NGC 6624 without an intermediate-mass black hole
• The orbital anisotropy profiles of nearby globular clusters from Gaia Data Release 2
• Spherical models of star clusters with potential escapers
• Multimass modelling of Milky Way globular clusters -- I. Implications on their stellar initial mass function above 1 M${}_{\odot }$
• Structure and Internal Kinematics of Nine Inner Milky Way Globular Clusters
• Truncated $\gamma$-exponential models for tidal stellar systems
• Investigating Dynamical Properties of Globular Clusters through a Family of Lowered Isothermal Models