On the unique evolutionary mechanisms of massive quiescent galaxies in the epoch of reionisation

Abstract: We investigate the evolutionary histories of a population of high mass, high redshift, quiescent galaxies in the cosmohydrodynamical simulation Thesan; studying the characteristic properties of their haloes and environments over the epoch of reionisation. Thesan employs a modified version of the Arepo moving-mesh code utilised in IllustrisTNG, whose explicit handling of ionising radiation couples haloes and galaxies with cosmic reionisation. This results in a regime of rapid growth and quenching, producing nine massive quiescent galaxies at $z = 5.5$, in a $95.5\text{cMpc}^3$ volume, with no counterpart in IllustrisTNG. We find that the rapid assembly of stellar mass is attributed to smooth accretion of mass onto their haloes in dense environments, particularly from massive neighbouring structures; while these galaxies exhibit fast-growing potential wells hosting massive central black holes. With insignificant merger activity, we find AGN feedback to be the primary source of quenching. We find megaparsec-scale densities and halo masses to continue growing after quenching, suggesting that massive quiescent galaxies will be found in some of the largest haloes and densest regions of space by observations. With an unprecedented plethora of massive quiescent galaxies found in JWST data, these reionisation simulations may provide the paradigm for the evolutionary histories of these galaxies, while constraining the relatively unknown physics of reionisation. Despite the rarity of simulated massive quiescent galaxies, their identification in Thesan enables identification of halo and environmental properties most conducive to their quenching, hopefully guiding future efforts into the study of analogous systems in deep observational surveys and large-volume N-body simulations.