Sub-millimetre galaxies in hierarchical models: revisiting the need for a top-heavy stellar initial mass function with Bayesian optimization
Year: 2026
Type: article
Abstract: 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 optimization, 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 data sets. 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 dn/dln$m \propto m^{-x}$). Our comprehensive parameter space exploration confirms that, given the implementation of galaxy formation physics within galform, a top-heavy IMF in starbursts is necessary to reconcile high-redshift dusty star-forming galaxies with local galaxy populations.
Institution: Durham University
Language: English
Subfield: Astronomy and Astrophysics
Field: Physics and Astronomy
Domain: Physical Sciences
Open Access status: gold
Estimated APC paid: $3,046
Awards: Durham DiRAC-2 petascale facility, G1468651527, G1805294933, DiRAC-2.5 - the pathway to DiRAC Phase 3, G2309215924, G2368581224, G2843074508, G4251892262, G4831232377, G629320885, DiRAC-2.5 DC - Operations 2017-2020, Astronomy at Durham 2023-2026, DiRAC: Memory Intensive 2.5y, Astronomy at Durham 2020-2023, STFC Centre for Doctoral Training in Data Intensive Science, DiRAC: Memory Intensive 2.5x
OpenAlex ID: W7140230605
