Bayesian Disk-Fitting with Irregular Flux Distributions


Mathew R. Varidel, Scott M. Croom, Geraint F. Lewis, Brendon J. Brewer, Enrico M. Di Teodoro, Joss Bland-Hawthorn, Julia J. Bryant, Christoph Federrath, Caroline Foster, Karl Glazebrook, Michael Goodwin, Brent Groves, Andrew M. Hopkins, Jon S. Lawrence, Angel R. Lopez-Sanchez, Anne M. Medling, Matt S. Owers, Samuel N. Richards, Richard Scalzo, Nicholas Scott, Sarah M. Sweet, Dan S. Taranu and Jesse van de Sande


We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H-alpha gas velocity dispersion for our sample to be in the range [7, 30] km/s. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is -0.29 +- 0.18 for the H-alpha flux-weighted mean velocity dispersion.

Update as of March 10, 2019: Accepted for MNRAS (view).

Publication Date: 
March 2019
Paper PDF: