The impact of quenching on galaxy profiles in the SIMBA simulation

We study specific star formation rate (sSFR) and gas profiles of star-forming (SF) and green valley (GV) galaxies in the SIMBA cosmological hydrodynamic simulation. SF galaxy half-light radii (R_half) at z = 0 and their evolution (∝(1 + z)^-0.78) agree with observations. Passive galaxy R_half agree with observations at high redshift, but by z = 0 are too large, owing to numerical heating. We compare SIMBAz = 0 sSFR radial profiles for SF and GV galaxies to observations. SIMBA shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in GV galaxies and massive SF systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, SIMBA yields GV galaxies with strongly dropping sSFR profiles beyond ≳R_half, regardless of active galactic nucleus feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. SIMBA's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At z = 2, SIMBA does not show central depressions in massive SF galaxies, suggesting SIMBA's X-ray feedback should be more active at high-z. High-resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in z = 0 GV galaxies represents a unique success of SIMBA. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.