A hierarchical Bayesian approach for reconstructing the Initial Mass Function of Single Stellar Populations

Abstract

Recent studies based on the integrated light of distant galaxies suggest that the initial mass function (IMF) might not be universal. Variations of the IMF with galaxy type and/or formation time may have important consequences for our understanding of galaxy evolution. We have developed a new stellar population synthesis (SPS) code specifically designed to reconstruct the IMF. We implement a novel approach combining regularization with hierarchical Bayesian inference. Within this approach we use a parametrized IMF prior to regulate a direct inference of the IMF. This direct inference gives more freedom to the IMF and allows the model to deviate from parametrized models when demanded by the data. We use Markov Chain Monte Carlo sampling techniques to reconstruct the best parameters for the IMF prior, the age, and the metallicity of a single stellar population. We present our code and apply our model to a number of mock single stellar populations with different ages, metallicities, and IMFs. When systematic uncertainties are not significant, we are able to reconstruct the input parameters that were used to create the mock populations. Our results show that if systematic uncertainties do play a role, this may introduce a bias on the results. Therefore, it is important to objectively compare different ingredients of SPS models. Through its Bayesian framework, our model is well-suited for this.

Accepted for publication in MNRAS, 30 pages, 28 figures

References in corpus (27)

• emcee: The MCMC Hammer
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• MILES: A Medium resolution INT Library of Empirical Spectra
• Evolution of asymptotic giant branch stars II. Optical to far-infrared isochrones with improved TP-AGB models
• Importance Nested Sampling and the MultiNest Algorithm
• Evolutionary Stellar Population Synthesis with MILES. Part I: The Base Models and a New Line Index System
• A systematic variation of the stellar initial mass function in early-type galaxies
• ULySS: A Full Spectrum Fitting Package
• The initial mass function of early-type galaxies
• Evolutionary stellar population synthesis with MILES - II. Scaled-solar and α-enhanced models
• Population Synthesis in the Blue IV: Accurate Model Predictions for Lick Indices and UBV Colors in Single Stellar Populations
• The galaxy stellar mass-star formation rate relation: Evidence for an evolving stellar initial mass function?
• CosmoSIS: modular cosmological parameter estimation
• MIUSCAT: extended MILES spectral coverage. I. Stellar populations synthesis models
• SPIDER VIII - Constraints on the Stellar Initial Mass Function of Early-type Galaxies from a Variety of Spectral Features
• Medium-resolution Isaac Newton Telescope library of empirical spectra - II. The stellar atmospheric parameters
• Counting Low-Mass Stars in Integrated Light
• Evidence of Cosmic Evolution of the Stellar Initial Mass Function
• Systematic variation of the stellar Initial Mass Function with velocity dispersion in early-type galaxies
• Evidence For A Mild Steepening And Bottom-Heavy IMF In Massive Galaxies From Sodium And Titanium-Oxide Indicators
• The stellar IMF determined in early-type galaxies from a non-degenerate set of optical line indices
• The atmospheric parameters and spectral interpolator for the stars of MILES
• From star clusters to dwarf galaxies: The properties of dynamically hot stellar systems
• A recalibration of optical photometry: Tycho-2, Strömgren, and Johnson systems
• Dissecting the Red Sequence - III. Mass-to-Light Variations in 3D Fundamental Plane Space
• The X-shooter Spectral Library (XSL): I. DR1. Near-ultraviolet through optical spectra from the first year of the survey

Cited by in corpus (7)

• The Extended IRTF Spectral Library: Expanded coverage in metallicity, temperature, and surface gravity
• The Stellar Initial Mass Function in Early-Type Galaxies From Absorption Line Spectroscopy. IV. A Super-Salpeter IMF in the center of NGC 1407 from Non-Parametric Models
• A Comprehensive Bayesian Discrimination of the Simple Stellar Population Model, Star Formation History and Dust Attenuation Law in the Spectral Energy Distribution Modeling of Galaxies
• Exploring the IMF of star clusters: a joint SLUG and LEGUS effort
• Detecting episodes of star formation using Bayesian model selection
• Hierarchical Bayesian inference of the Initial Mass Function in Composite Stellar Populations
• Counting stars from the integrated spectra of galaxies