Sub-millimeter galaxies in hierarchical models: revisiting the need for a top-heavy stellar initial mass function with Bayesian optimisation

TL;DR

This study uses Bayesian optimisation to efficiently explore galaxy formation models, demonstrating that a top-heavy stellar initial mass function in starbursts is essential to match both high-redshift sub-millimeter galaxy observations and local galaxy properties.

Contribution

It introduces an extensive Bayesian parameter search in GALFORM, confirming the necessity of a top-heavy IMF in starbursts for consistent galaxy evolution modeling.

Findings

Universal IMF models fail to fit all data simultaneously.

Allowing a top-heavy IMF improves fit to both high-redshift and local observations.

Optimal IMF slope found to be approximately 0.7, more top-heavy than recent models.

Abstract

The properties of high-redshift sub-millimetre galaxies (SMGs) remain controversial within hierarchical structure formation models. We revisit whether a top-heavy stellar initial mass function (IMF) in starbursts is required to reproduce both SMG observations and local galaxy properties. Using Bayesian optimisation, we perform an extensive search of the 15-dimensional parameter space of the GALFORM semi-analytical model. This efficient approach converges to optimal parameter values in fewer than 200 model evaluations, representing orders of magnitude fewer runs than traditional methods. We test whether GALFORM can simultaneously match three key observational constraints: the $z=0$ $K$-band luminosity function, the SMG number counts at 850~$\mu$m, and the SMG redshift distribution. We consider two model variants: one with a universal solar neighbourhood IMF for all star formation, and another allowing the IMF slope in starbursts to vary as a free parameter. When assuming a universal Chabrier IMF, we find no parameter combination that simultaneously reproduces all three datasets. The model either matches the SMG constraints while grossly overpredicting the local $K$-band luminosity function, or matches the local luminosity function while severely underpredicting SMG counts by factors of 3--100. In contrast, allowing a top-heavy IMF in starbursts enables excellent simultaneous fits to all constraints. The best-fitting model prefers an IMF slope parameter $x \approx 0.7$ (where d$n$/dlog$m \propto m^{-x}$), somewhat more top-heavy than recent models but less extreme than early proposals. Our comprehensive parameter space exploration definitively confirms that, within the GALFORM framework, a top-heavy IMF in starbursts is necessary to reconcile high-redshift dusty star-forming galaxies with local galaxy populations.