The SAMI Galaxy Survey: Bayesian inference for gas disc kinematics using a hierarchical Gaussian mixture model. Issue 3 (8th March 2019)

• Record Type: Journal Article
• Authors: Varidel, Mathew R; Croom, Scott M; Lewis, Geraint F; Brewer, Brendon J; Di Teodoro, Enrico M; Bland-Hawthorn, Joss; Bryant, Julia J; Federrath, Christoph; Foster, Caroline; Glazebrook, Karl; Goodwin, Michael; Groves, Brent; Hopkins, Andrew M; Lawrence, Jon S; López-Sánchez, Ángel R; Medling, Anne M; Owers, Matt S; Richards, Samuel N; Scalzo, Richard; Scott, Nicholas; Sweet, Sarah M; Taranu, Dan S; van de Sande, Jesse
• Is Part Of: Monthly notices of the Royal Astronomical Society. Volume 485:Issue 3(2019)
• Page Start: 4024
• Page End: 4044
• Publication Date: 2019-03-08
• Subjects: methods: data analysis -- methods: statistical -- techniques: imaging spectroscopy -- galaxies: kinematics and dynamics; Astronomy -- Periodicals; Periodicals; 520.5
• DOI: 10.1093/mnras/stz670
• Languages: English
• ISSNs: 0035-8711

Abstract

Abstract: 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. The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disc 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 α gas velocity dispersion for our sample to be in the range $\bar{\sigma }_v \sim$ [7, 30] km s −1 . The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $\Delta \bar{\sigma }_v / \bar{\sigma }_v = - 0.29 \pm 0.18$ for the H α flux-weighted mean velocity dispersion.