Near-Gaussian distributions for modelling discrete stellar velocity data with heteroskedastic uncertainties
arXiv:2009.07858 · doi:10.1093/mnras/staa2860
Abstract
The velocity distributions of stellar tracers in general exhibit weak non-Gaussianity encoding information on the orbital composition of a galaxy and the underlying potential. The standard solution for measuring non-Gaussianity involves constructing a series expansion (e.g. the Gauss-Hermite series) which can produce regions of negative probability density. This is a significant issue for the modelling of discrete data with heteroskedastic uncertainties. Here, we introduce a method to construct positive-definite probability distributions by the convolution of a given kernel with a Gaussian distribution. Further convolutions by observational uncertainties are trivial. The statistics (moments and cumulants) of the resulting distributions are governed by the kernel distribution. Two kernels (uniform and Laplace) offer simple drop-in replacements for a Gauss-Hermite series for negative and positive excess kurtosis distributions with the option of skewness. We demonstrate the power of our method by an application to real and mock line-of-sight velocity datasets on dwarf spheroidal galaxies, where kurtosis is indicative of orbital anisotropy and hence a route to breaking the mass-anisotropy degeneracy for the identification of cusped versus cored dark matter profiles. Data on the Fornax dwarf spheroidal galaxy indicate positive excess kurtosis and hence favour a cored dark matter profile. Although designed for discrete data, the analytic Fourier transforms of the new models also make them appropriate for spectral fitting, which could improve the fits of high quality data by avoiding unphysical negative wings in the line-of-sight velocity distribution.
21 pages, 12 figures, accepted for publication in MNRAS, slight adjustment of parameterization to match normalized definition of Gauss-Hermite coefficients (results unchanged)
References in corpus (18)
- The NumPy array: a structure for efficient numerical computation
- The Gaia mission
- Improving the full spectrum fitting method: accurate convolution with Gauss-Hermite functions
- Evolutionary Stellar Population Synthesis with MILES. Part I: The Base Models and a New Line Index System
- The DART imaging and CaT survey of the Fornax Dwarf Spheroidal Galaxy
- Stellar Velocities in the Carina, Fornax, Sculptor and Sextans dSph Galaxies: Data from the Magellan/MMFS Survey
- The SAMI Galaxy Survey: Revisiting Galaxy Classification Through High-Order Stellar Kinematics
- The mass and velocity anisotropy of the Carina, Fornax, Sculptor and Sextans dwarf spheroidal galaxies
- Clean Kinematic Samples in Dwarf Spheroidals: An Algorithm for Evaluating Membership and Estimating Distribution Parameters When Contamination is Present
- Systemic Proper Motions of Milky Way Satellites from Stellar Redshifts: the Carina, Fornax, Sculptor and Sextans Dwarf Spheroidals
- Cores and Cusps in the Dwarf Spheroidals
- Action-based distribution functions for spheroidal galaxy components
- The distance to the Fornax Dwarf Spheroidal Galaxy
- Three Mechanisms for Bar Thickening
- 3D motions in the Sculptor dwarf galaxy as a glimpse of a new era
- Made-to-Measure Dark Matter Haloes, Elliptical Galaxies and Dwarf Galaxies in Action Coordinates
- Models of Bars I: Flattish Profiles for Early-Type Galaxies
- Models of Bars with Exponential Density Profiles