Sensitivity analysis of a galaxy formation model

Abstract

We present the first application of a variance-based sensitivity analysis (SA) to a model that aims to predict the evolution and properties of the whole galaxy population. SA is a well-established technique in other quantitative sciences, but is a relatively novel tool for the evaluation of astrophysical models. We perform a multi-parameter exploration of the GALFORM semi-analytic galaxy formation model, to compute how sensitive the present-day K-band luminosity function is to varying different model parameters. The parameter space is scanned using a low-discrepancy sampling technique proposed by Saltelli. We first demonstrate the usefulness of the SA approach by varying just two model parameters, one which controls supernova feedback and the other the heating of gas by AGN. The SA analysis matches our physical intuition regarding how these parameters affect the predictions for different parts of the galaxy luminosity function. We then use SA to compute Sobol' sensitivity indices varying seven model parameters, connecting the variance in the model output to the variance in the input parameters. The sensitivity is computed in luminosity bins, allowing us to probe the origin of the model predictions in detail. We discover that the SA correctly identifies the least- and most important parameters. Moreover, the SA also captures the combined responses of varying multiple parameters at the same time. Our study marks a much needed step away from the traditional "one-at-a-time" parameter variation often used in this area and improves the transparency of multi-parameter models of galaxy formation.

16 pages, 10 figures

References in corpus (10)

• The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations
• The Role of Pressure in GMC Formation II: The H_2 - Pressure Relation
• A primer on hierarchical galaxy formation: the semi-analytical approach
• Monte Carlo Markov Chain Parameter Estimation in Semi-Analytic Models of Galaxy Formation
• Comparing galaxy formation in semi-analytic models and hydrodynamical simulations
• Why So Many Published Sensitivity Analyses Are False. A Systematic Review of Sensitivity Analysis Practices
• Building a Predictive Model of Galaxy Formation - I: Phenomenological Model Constrained to the $z=0$ Stellar Mass Function
• Constraints on galaxy formation models from the galaxy stellar mass function and its evolution