2018/06/26 by J. J. Argyle, Joshua Argyle, J. Méndez-Abreu +5 · 1 citation
Mathematics · Physics and Astronomy · #Advanced Statistical Methods and Models #Algorithm #Artificial intelligence #Astrophysics #Bayesian probability #Bulge #Computer science #Galaxies: Formation, Evolution, Phenomena #Galaxy #Markov chain Monte Carlo #Physics #Statistical Methods and Inference #Surface brightness #astro-ph.GA
paper · pdf · doi:10.1093/mnras/sty1691
21 pages, 16 figures. Accepted for publication in MNRAS
openalex publication_date 2018/06/26 · arxiv created 2018/07/05 · arxiv updated 2018/07/11 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
We introduce phi, a fully Bayesian Markov chain Monte Carlo algorithm designed for the structural decomposition of galaxy images. phi uses a triple layer approach to effectively and efficiently explore the complex parameter space. Combining this with the use of priors to prevent non-physical models, phi offers a number of significant advantages for estimating surface brightness profile parameters over traditional optimization algorithms. We apply phi to a sample of synthetic galaxies with Sloan Digital Sky Survey (SDSS)-like image properties to investigate the effect of galaxy properties on our ability to recover unbiased and well-constrained structural parameters. In two-component bulge+disc galaxies, we find that the bulge structural parameters are recovered less well than those of the disc, particularly when the bulge contributes a lower fraction to the luminosity, or is barely resolved with respect to the pixel scale or point spread function (PSF). There are few systematic biases, apart from for bulge+disc galaxies with large bulge Sérsic parameter, n. On application to SDSS images, we find good agreement with other codes, when run on the same images with the same masks, weights, and PSF. Again, we find that bulge parameters are the most difficult to constrain robustly. Finally, we explore the use of a Bayesian information criterion method for deciding whether a galaxy has one or two components.