The nature and evolution of early massive quenched galaxies in the simba-C simulation

We examine the nature, origin, and fate of early ( z >= 2) massive ( M-* > 10(10) M-circle dot) quenched galaxies (EQGs) in a new (100 h(-1) Mpc(3)) run of the SIMBA-C galaxy formation model. We define 'quenched' to be > 4 sigma below an iterative polynomial fit to the star-forming sequence (SFS), and find that SIMBA-C produces EQGs as early as z similar to 5 and number densities agreeing with observations at z less than or similar to 3 (though slightly low at z > 4). Using a photometric-based EQG selection or a fixed specific star formation rate cut of 10(-10) yr(-1) yields similar results. EQGs predominantly arise in central galaxies with stellar mass M * similar to 10(10 . 5 -11.3) M-circle dot, not necessarily the most massive systems. A Uniform Manifold Approximation and Projection shows that quenched galaxies have notably large black hole-to-stellar mass ratios, lower rotational support, and less dust, but are not atypical versus similar-mass non-EQGs in their environments, halo mass, or halo gas temperatures at the time of quenching. However, via galaxy tracking we show that the progenitor environments of EQGs are significantly more overdense than that of non-EQGs, which drives higher black hole mass fractions and stellar-to-halo mass ratios. This results in the Eddington ratio dropping sufficiently low for SIMBA-C's jet mode feedback to turn on, which quickly quenches the host galaxies. EQGs thus seem to be galaxies that grow their black holes quickly within highly dense environments, but end up in moderately dense environments where black hole feedback can quench effectively. We find that greater than or similar to 30 per cent of EQGs rejuvenate, but the rejuvenating fraction drops quickly below z less than or similar to 2. By z = 0, it is difficult to distinguish the descendants of EQGs versus non-EQGs.