Filip Huško

Recent observations with James Webb Space Telescope have revealed a population of UV-bright galaxies at $z\gtrsim 10$. This discovery naturally raises the question: What do such early galaxies evolve into by the present day? In this work, we address this descendant question using the new-generation COLIBRE (COLd ISM and Better REsolution) cosmological hydrodynamical simulations to trace bright galaxies selected at $z=10$ and follow their descendants to the present day. Most of the high-redshift galaxies do not survive as distinct, self-bound objects to $z=0$; instead, the majority are incorporated into more massive systems through merging or disruption. The surviving descendants span a broad range of present-day stellar masses, although they are most commonly intermediate- to high-mass, $M_\star \sim 10^{10}$–$10^{11}\, {\rm M_\odot }$. They typically reside in galaxy groups and clusters, with host halo masses, $M_{200{\rm c}}\sim 10^{13}$–$10^{14}\, {\rm M_\odot }$. The large scatter in descendant stellar mass shows that present-day outcomes retain only a weak memory of the stellar mass of the high-redshift progenitor. We show that the evolution of descendant host halo masses is consistent with the forward conditional distribution predicted by extended Press–Schechter (EPS) theory, both in the median growth and in the large scatter in descendant mass. In particular, EPS confirms that massive present-day galaxies typically do not originate from the most massive objects at high redshift. A galaxy observed at $z\gtrsim 10$ therefore cannot be interpreted as the direct progenitor of a single class of $z=0$ galaxies.

James Webb Space Telescope (JWST) has enabled the detection of galaxies in the earliest stages of cosmic history. We compare the ultraviolet luminosity functions (UVLFs) at redshifts $z=7\mathrm{-}20$ predicted by a set of new cosmological hydrodynamical simulations, colibre (COLd ISM and Better REsolution), with observations, including those from JWST. The UV luminosities of colibre galaxies are derived using the radiative transfer code skirt, which tracks stellar emission and its processing through the multiphase interstellar medium and dust distribution predicted by colibre. We find that although colibre is consistent with the observed evolution of the stellar mass function up to $z=12$, its dust-attenuated UVLFs fall systematically below the observations at the bright end: at the number density of $10^{-6}\, \mathrm{Mpc^{-3}\, mag^{-1}}$, the brightest galaxies are underluminous by $\approx 1\, \rm mag$ at $z=7$, increasing to $\approx 2.5\, \rm mag$ at $z=15$. Accounting for observational uncertainties brings the colibre UVLFs closer to the observational data, but does not fully resolve the discrepancy. Ignoring dust attenuation allows colibre to produce sufficiently bright galaxies at $7\lesssim z \lesssim 12$, but at $z=15$, colibre still underpredicts the luminosities of the brightest galaxies, indicating the need for additional physical mechanisms to boost the UV luminosities at the earliest cosmic epochs, such as a ‘top-heavy’ stellar initial mass function. We fit the colibre UVLFs with Schechter functions and calculate the evolution of the best-fitting parameters. We find that the galaxy number density decreases, the characteristic luminosity becomes fainter and the faint-end slope becomes steeper towards higher redshifts. The UV luminosity density decreases by a factor of several hundred from $z = 7$ to $15$.

Early ($z \ge 2$) Massive ($M_{\star } \gt 10^{10}\, \mathrm{{\rm M}_{\odot }}$), Quenched Galaxies (MQGs) challenge current galaxy formation models. In this series, we study these systems using the new colibre cosmological hydrodynamical simulations. Following the broad agreement between their predictions and observations found in the first paper, this second paper explores the processes driving galaxies to become massive and quenched in colibre, identifying Active Galactic Nucleus (AGN) feedback as the primary quenching mechanism in both the thermal (L200m6 simulation) and hybrid (thermal + jet, L200m7h simulation) AGN feedback models implemented. However, the two models behave differently: while the thermal model efficiently quenches massive galaxies at $z\gt 3$, the hybrid model is less effective because black holes (BHs) grow more slowly in the early Universe, and the jet component, which dominates the feedback energy, acts on longer time-scales to impact galaxies. Both models predict quasar-like MQGs (AGN with $L_{\rm bol}\gtrsim 10^{45}\, \mathrm{erg\, s^{-1}}$), with the most luminous systems associated with more recently quenched galaxies. Compared to star-forming galaxies of similar mass, MQGs host more massive BHs and exhibit higher star formation efficiencies. These differences arise primarily from their environments before quenching, particularly at local ($\rm 0.3\, cMpc$) to intermediate scales ($\rm 1.0\, cMpc$), where overdense regions are associated with enhanced gas inflows, higher BH accretion and, hence, feedback power. We find that about 54 per cent (20 per cent) of the $z=3$ MQGs survive as the main progenitors of $z=0$ galaxies, although up to 56 per cent (60 per cent) experience rejuvenation episodes at a given redshift $z\lt 3$ in L200m6 (L200m7h). Our results highlight the central role of BH growth, AGN feedback and environment in driving rapid quenching in the early Universe.

We present predictions from the colibre cosmological hydrodynamical simulations of galaxy formation for the present-day galaxy luminosity functions (LFs) at wavelengths ranging from the far-ultraviolet (FUV) to the submillimetre. The simulations are post-processed with the radiative transfer code skirt, accounting for dust attenuation and emission using the distribution and properties of dust grains predicted directly by colibre. Results from simulations varying in mass resolution by a factor of $\sim 10^2$ ($\sim 10^5\text {-} 10^7\, \mathrm{{\rm M}_{\odot }}$) show very good convergence over most luminosity ranges. The colibre-skirt LFs match the data remarkably well from the FUV to $3.4\, \rm{\mu m}$ and also in the far-infrared and submillimetre wavelength range ($70\text {-}850\, \rm{\mu m}$). In the mid-infrared (MIR; $8\text {-}24\, \rm{\mu m}$), colibre-skirt matches the data well at low luminosities but significantly underpredicts the luminosities of MIR-bright galaxies, with the discrepancy increasing towards longer wavelengths. The total infrared LF, obtained by integrating the spectral energy distributions over $8\text {-}1000\, \rm{\mu m}$, also matches observations well at the faint end but underpredicts the number of very bright galaxies. The unprecedented agreement at all other wavelengths indicates that colibre, coupled with this calibration-free skirt post-processing framework, successfully predicts the properties of stellar populations at the present day and the amount and distribution of interstellar dust.

The diversity of galaxy morphologies and their relations with galaxy and halo properties is fundamental to understanding galaxy formation. Cosmological simulations of representative volumes can help disentangle the origin of observed correlations, but most suffer from two main limitations that affect morphologies: an overpressurized interstellar medium and spurious interactions between stellar and dark matter particles. We present an overview of galaxy morphologies in the COLIBRE simulations, which address these limitations and reproduce many observed galaxy scaling relations. To quantify galaxy morphology, we use four (strongly correlated) theory-space metrics, three kinematic and one spatial. We explore how different choices and limitations affect these indicators, including luminosity- versus mass-weighting, aperture size and shot noise. Overall, we find good convergence in present-day morphologies across two orders of magnitude in mass resolution. COLIBRE predicts that kinematic morphology correlates strongly with stellar mass and colour, and that galaxies with stellar masses of $\approx (1-2)\times 10^{10}\, \, \mathrm{M}_{\odot }$ tend to be the most rotationally dominated. At fixed stellar mass, the morphology of central galaxies correlates weakly with the properties of their host halo. Morphology correlates more strongly with internal galaxy properties, with more discy galaxies being more gas-rich, having higher star formation rates and exhibiting younger and more extended stellar populations. Other properties, like the mass of the most massive black hole, the fraction of stars that are accreted and stellar metallicity, also correlate with morphology, but with correlation strengths sensitive to the stellar mass of the galaxy and whether it is a central or satellite.

The James Webb Space Telescope has uncovered a substantial population of Massive ($M_{\star } \gt 10^{10}\, \mathrm{M_{\odot }}$), Quenched Galaxies (MQGs) in the early Universe ($z \ge 2$), whose properties challenge current galaxy formation models. In this series, we examine this population of MQGs within the new colibre cosmological hydrodynamical simulations, which introduce key innovations in their sub-grid physics. In this first paper, we find a dependence of MQG number densities on both mass resolution and the Active Galactic Nucleus feedback implementation, as well as a significant impact from potential observational uncertainties. Using the fiducial $(200\, \rm cMpc)^3$ volume L200m6 simulation, which provides adequate volume, mass, and spatial resolution to study these systems, we report number densities and stellar mass functions in broad agreement with the latest observations. The predicted quenching and formation time-scales are qualitatively consistent with observational inferences, indicating extended formation (medians $t_{50}\approx 0.5{-}1.5\, \mathrm{Gyr}$) followed by rapid quenching (medians $t_{\mathrm{q}}\lesssim 0.6\, \mathrm{Gyr}$) with strong starburst episodes. Leveraging the state-of-the-art physics in colibre, the model predicts that MQGs have dust and $\rm H_{2}$ fractions more than 1 dex lower than their massive star-forming counterparts; generally consistent with the (scarce) observational estimates. MQGs and massive star-forming systems show broadly similar stellar sizes and kinematics, suggesting that size or morphological transformations occur after quenching in colibre. Our results provide robust predictions for MQGs and show that tensions with observations are reduced when an effective observational uncertainty is forward-modelled.

We investigate the co-evolution of the stellar and dark matter mass of satellite galaxies using the colibre cosmological hydrodynamical simulations with subhaloes resolved by the history-based hbt-herons subhalo finder. We identify a universal tidal track connecting stellar mass-loss to subhalo mass-loss characterized by two distinct phases, which can be well described by the two-parameter model. The initial phase consists primarily of dark matter stripping, whereas stellar stripping becomes significant only after the subhalo total bound mass fraction drops below a critical value ($\sim 0.057$). We find a bimodal mass-loss rate distribution of subhaloes. In satellites with modest mass-loss rates, the stellar mass is largely frozen. By contrast, the galaxy quickly becomes unresolved, along with the dark matter component for the extreme mass-loss population, naturally explaining the lack of ‘orphan’ galaxies in previous hydrodynamical simulations. Our model also predicts the formation condition for dark-matter-deficient galaxies (DMDGs), whose abundance peaks at $m_{*}\sim 10^{9.5}\, \rm {M}_{\odot }$. The abundance of DMDGs can be very sensitive to numerical effects, with colibre resolving a much larger DMDG population than previous hydrodynamical simulations. We also estimate the influence of artificial disruption on the satellite stellar mass function, which can amount to 20 (50) per cent at $m_* \sim 10^{9} (10^{8}) \, \rm M_\odot$, given a baryonic mass resolution of $\sim 10^{6}\, \rm {M}_{\odot }$.

We present the COLIBRE galaxy formation model and the COLIBRE suite of cosmological hydrodynamical simulations. COLIBRE includes new models for radiative cooling, dust grains, star formation, stellar mass loss, turbulent diffusion, pre-supernova stellar feedback, supernova feedback, supermassive black holes, and active galactic nucleus (AGN) feedback. The multiphase interstellar medium is explicitly modelled without a pressure floor. Hydrogen and helium are tracked in non-equilibrium, with their contributions to the free electron density included in metal-line cooling calculations. The chemical network is coupled to a dust model that tracks three grain species and two grain sizes. In addition to the fiducial thermally driven AGN feedback, a subset of simulations uses black hole spin-dependent hybrid jet/thermal AGN feedback. To suppress spurious transfer of energy from dark matter to stars, dark matter is supersampled by a factor 4, yielding similar dark matter and baryonic particle masses. The subgrid feedback model is calibrated to match the observed $z\approx 0$ galaxy stellar mass function, galaxy sizes, and black hole masses in massive galaxies. The COLIBRE suite includes three resolutions, with particle masses of $\sim 10^5$, $10^6$, and $10^7\, \text{M}_\odot$ in cubic volumes of up to 100, 200, and 400 cMpc on a side, respectively. The largest runs use 136 billion ($5\times 3008^3$) particles. We describe the model, assess its strengths and limitations, and present both visual impressions and quantitative results. Comparisons with various low-redshift galaxy observations generally show very good numerical convergence and excellent agreement with the data.

We present the calibration of stellar and active galactic nucleus (AGN) feedback in the subgrid model for the new colibre hydrodynamical simulations of galaxy formation. colibre directly simulates the multiphase interstellar medium and the evolution of dust grains, which is coupled to the chemistry. colibre is calibrated at three resolutions: particle masses of $m_{\rm gas} \approx m_{\rm dm} \sim 10^7$ (m7), $10^6$ (m6), and $10^5\,\mathrm{M_\odot }$ (m5). To calibrate the colibre feedback at m7 resolution, we run Latin hypercubes of $\approx 200$ simulations that vary up to four subgrid parameters in cosmological volumes of ($50\,\mathrm{cMpc})^{3}$. We train Gaussian process emulators on these simulations to predict the $z=0$ galaxy stellar mass function (GSMF) and size–stellar mass relation (SSMR) as functions of the model parameters, which we then fit to observations. The trained emulators not only provide the best-fitting parameter values but also enable us to investigate how different aspects of the prescriptions for supernova and AGN feedback affect the predictions. In particular, we demonstrate that while the observed $z=0$ GSMF and SSMR can be matched individually with a relatively simple supernova feedback model, simultaneously reproducing both necessitates a more sophisticated prescription. We show that the calibrated m7 colibre model not only reproduces the calibration target observables, but also matches various other galaxy properties to which the model was not calibrated. Finally, we apply the calibrated m7 model to the m6 and m5 resolutions and, after slight manual adjustments of the subgrid parameters, achieve a similar level of agreement with the observed $z=0$ GSMF and SSMR.