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.
arXiv comment21 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 (53)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Astropy: A Community Python Package for Astronomy
- emcee: The MCMC Hammer
- The NumPy array: a structure for efficient numerical computation
- A Universal Density Profile from Hierarchical Clustering
- Gaia Data Release 2. Summary of the contents and survey properties
- The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
- The Gaia mission
- Dwarf Galaxies of the Local Group
- Parametric Recovery of Line-of-Sight Velocity Distributions from Absorption-Line Spectra of Galaxies via Penalized Likelihood
- The observed properties of dwarf galaxies in and around the Local Group
- Improving the full spectrum fitting method: accurate convolution with Gauss-Hermite functions
- The Atlas3D project -- I. A volume-limited sample of 260 nearby early-type galaxies: science goals and selection criteria
- Evolutionary Stellar Population Synthesis with MILES. Part I: The Base Models and a New Line Index System
- A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies
- Gaia Data Release 2: Kinematics of globular clusters and dwarf galaxies around the Milky Way
- AGAMA: Action-based galaxy modelling architecture
- The SAURON project - III. Integral-field absorption-line kinematics of 48 elliptical and lenticular galaxies
- Dynamical Modeling of Velocity Profiles: The Dark Halo Around the Elliptical Galaxy NGC2434
- Galaxia: a code to generate a synthetic survey of the Milky Way
- The DART imaging and CaT survey of the Fornax Dwarf Spheroidal Galaxy
- Expansions for nearly Gaussian distributions
- Stellar Velocities in the Carina, Fornax, Sculptor and Sextans dSph Galaxies: Data from the Magellan/MMFS Survey
- Improved evidence for a black hole in M32 from HST/FOS spectra - II. Axisymmetric dynamical models
- Dark matter heats up in dwarf galaxies
- Dark Matter Cores and Cusps: The Case of Multiple Stellar Populations in Dwarf Spheroidals
- The SAMI Galaxy Survey: Revisiting Galaxy Classification Through High-Order Stellar Kinematics
- Dark matter distribution in dwarf spheroidal galaxies
- 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
- The Kinematic Signature of Face-On Peanut-Shaped Bulges
- Dark matter distribution in the Draco dwarf from velocity moments
- The core size of the Fornax dwarf Spheroidal
- 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
- Model comparison of the dark matter profiles of Fornax, Sculptor, Carina and Sextans
- Action-based distribution functions for spheroidal galaxy components
- Galaxy structure from multiple tracers: II. M87 from parsec to megaparsec scales
- Phase Space Models of the Dwarf Spheroidal Galaxies
- The distance to the Fornax Dwarf Spheroidal Galaxy
- Three Mechanisms for Bar Thickening
- Mass profiles and galaxy orbits in nearby galaxy clusters from the analysis of the projected phase-space
- 3D motions in the Sculptor dwarf galaxy as a glimpse of a new era
- Action-based dynamical models of dwarf spheroidal galaxies: application to Fornax
- Analytical Solutions to the Mass-Anisotropy Degeneracy with Higher Order Jeans Analysis: A General Method
- Dynamical Models of Elliptical Galaxies -- II. M87 and its Globular Clusters
- The morphology and structure of stellar populations in the Fornax dwarf spheroidal galaxy from Dark Energy Survey Data
- Made-to-Measure Dark Matter Haloes, Elliptical Galaxies and Dwarf Galaxies in Action Coordinates
- Inner bars also buckle. The MUSE TIMER view of the double-barred galaxy NGC 1291
- Stellar Density Profiles of Dwarf Spheroidal Galaxies
- Line Profiles from Discrete Kinematic Data
- Models of Bars I: Flattish Profiles for Early-Type Galaxies
- Models of Bars with Exponential Density Profiles