skip to main content

Title: Massive central galaxies of galaxy groups in the Romulus simulations: an overview of galaxy properties at z  = 0

Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central galaxies in group-scale haloes at z = 0 to observed BGGs. The comparison encompasses the stellar mass–halo mass relation, various kinematic properties and scaling relations, morphologies, and the star formation rates. Generally, we find that Romulus reproduces the full spectrum of diversity in the properties of the BGGs very well, albeit with a tendency toward lower than the observed fraction of quenched BGGs. We find both early-type S0 and elliptical galaxies as well as late-type disc galaxies; we find Romulus galaxies that are fast-rotators as well as slow-rotators; and we observe galaxies transforming from late-type to early-type following strong dynamical interactions with satellites. We also carry out case studies of selected Romulus galaxies to explore the link between their properties, and the recent evolution of the stellar system as well as the surrounding intragroup/circumgalactic medium. In general, mergers/strong interactions quench star-forming activity and disrupt the more » stellar disc structure. Sometimes, however, such interactions can also trigger star formation and galaxy rejuvenation. Black hole feedback can also lead to a decline of the star formation rate but by itself, it does not typically lead to complete quenching of the star formation activity in the BGGs.

« less
; ; ; ; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 22-47
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The kinematic morphology–density relation of galaxies is normally attributed to a changing distribution of galaxy stellar masses with the local environment. However, earlier studies were largely focused on slow rotators; the dynamical properties of the overall population in relation to environment have received less attention. We use the SAMI Galaxy Survey to investigate the dynamical properties of ∼1800 early and late-type galaxies with log (M⋆/M⊙) > 9.5 as a function of mean environmental overdensity (Σ5) and their rank within a group or cluster. By classifying galaxies into fast and slow rotators, at fixed stellar mass above log (M⋆/M⊙) > 10.5, we detect a higher fraction (∼3.4σ) of slow rotators for group and cluster centrals and satellites as compared to isolated-central galaxies. We find similar results when using Σ5 as a tracer for environment. Focusing on the fast-rotator population, we also detect a significant correlation between galaxy kinematics and their stellar mass as well as the environment they are in. Specifically, by using inclination-corrected or intrinsic $\lambda _{R_{\rm {e}}}$ values, we find that, at fixed mass, satellite galaxies on average have the lowest $\lambda _{\, R_{\rm {e}},\rm {intr}}$, isolated-central galaxies have the highest $\lambda _{\, R_{\rm {e}},\rm {intr}}$, and group and clustermore »centrals lie in between. Similarly, galaxies in high-density environments have lower mean $\lambda _{\, R_{\rm {e}},\rm {intr}}$ values as compared to galaxies at low environmental density. However, at fixed Σ5, the mean $\lambda _{\, R_{\rm {e}},\rm {intr}}$ differences for low and high-mass galaxies are of similar magnitude as when varying Σ5 ($\Delta \lambda _{\, R_{\rm {e}},\rm {intr}} \sim 0.05$, with σrandom = 0.025, and σsyst < 0.03). Our results demonstrate that after stellar mass, environment plays a significant role in the creation of slow rotators, while for fast rotators we also detect an independent, albeit smaller, impact of mass and environment on their kinematic properties.« less

    Misalignments between the rotation axis of stars and gas are an indication of external processes shaping galaxies throughout their evolution. Using observations of 3068 galaxies from the SAMI Galaxy Survey, we compute global kinematic position angles for 1445 objects with reliable kinematics and identify 169 (12 per cent) galaxies which show stellar-gas misalignments. Kinematically decoupled features are more prevalent in early-type/passive galaxies compared to late-type/star-forming systems. Star formation is the main source of gas ionization in only 22 per cent of misaligned galaxies; 17 per cent are Seyfert objects, while 61 per cent show Low-Ionization Nuclear Emission-line Region features. We identify the most probable physical cause of the kinematic decoupling and find that, while accretion-driven cases are dominant, for up to 8 per cent of our sample, the misalignment may be tracing outflowing gas. When considering only misalignments driven by accretion, the acquired gas is feeding active star formation in only ∼1/4 of cases. As a population, misaligned galaxies have higher Sérsic indices and lower stellar spin and specific star formation rates than appropriately matched samples of aligned systems. These results suggest that both morphology and star formation/gas content are significantly correlated with the prevalence and timescales of misalignments. Specifically, torques on misaligned gas discs are smaller for more centrallymore »concentrated galaxies, while the newly accreted gas feels lower viscous drag forces in more gas-poor objects. Marginal evidence of star formation not being correlated with misalignment likelihood for late-type galaxies suggests that such morphologies in the nearby Universe might be the result of preferentially aligned accretion at higher redshifts.

    « less

    We study the alignments of galaxy spin axes with respect to cosmic web filaments as a function of various properties of the galaxies and their constituent bulges and discs. We exploit the SAMI Galaxy Survey to identify 3D spin axes from spatially resolved stellar kinematics and to decompose the galaxy into the kinematic bulge and disc components. The GAMA survey is used to reconstruct the cosmic filaments. The mass of the bulge, defined as the product of stellar mass and bulge-to-total flux ratio Mbulge = M⋆ × (B/T), is the primary parameter of correlation with spin–filament alignments: galaxies with lower bulge masses tend to have their spins parallel to the closest filament, while galaxies with higher bulge masses are more perpendicularly aligned. M⋆ and B/T separately show correlations, but they do not fully unravel spin–filament alignments. Other galaxy properties, such as visual morphology, stellar age, star formation activity, kinematic parameters, and local environment, are secondary tracers. Focussing on S0 galaxies, we find preferentially perpendicular alignments, with the signal dominated by high-mass S0 galaxies. Studying bulge and disc spin–filament alignments separately reveals additional information about the formation pathways of the corresponding galaxies: bulges tend to have more perpendicular alignments, while discs showmore »different tendencies according to their kinematic features and the mass of the associated bulge. The observed correlation between the flipping of spin–filament alignments and the growth of the bulge can be explained by mergers, which drive both alignment flips and bulge formation.

    « less
  4. It is now well established that galaxies have different morphologies, gas contents, and star formation rates (SFR) in dense environments like galaxy clusters. The impact of environmental density extends to several virial radii, and galaxies appear to be pre-processed in filaments and groups before falling into the cluster. Our goal is to quantify this pre-processing in terms of gas content and SFR, as a function of density in cosmic filaments. We have observed the two first CO transitions in 163 galaxies with the IRAM-30 m telescope, and added 82 more measurements from the literature, thus forming a sample of 245 galaxies in the filaments around the Virgo cluster. We gathered HI-21cm measurements from the literature and observed 69 galaxies with the Nançay telescope to complete our sample. We compare our filament galaxies with comparable samples from the Virgo cluster and with the isolated galaxies of the AMIGA sample. We find a clear progression from field galaxies to filament and cluster galaxies for decreasing SFR, increasing fraction of galaxies in the quenching phase, an increasing proportion of early-type galaxies, and decreasing gas content. Galaxies in the quenching phase, defined as having a SFR below one-third of that of the main sequencemore »(MS), are only between 0% and 20% in the isolated sample, according to local galaxy density, while they are 20%–60% in the filaments and 30%–80% in the Virgo cluster. Processes that lead to star formation quenching are already at play in filaments; they depend mostly on the local galaxy density, while the distance to the filament spine is a secondary parameter. While the HI-to-stellar-mass ratio decreases with local density by an order of magnitude in the filaments, and two orders of magnitude in the Virgo cluster with respect to the field, the decrease is much less for the H 2 -to-stellar-mass ratio. As the environmental density increases, the gas depletion time decreases, because the gas content decreases faster than the SFR. This suggests that gas depletion precedes star formation quenching.« less
  5. ABSTRACT We examine the 1.4 GHz radio luminosities of galaxies arising from star formation and active galactic nuclei (AGNs) within the state-of-the-art cosmological hydrodynamic simulation Simba. Simba grows black holes via gravitational torque limited accretion from cold gas and Bondi accretion from hot gas, and employs AGN feedback including jets at low Eddington ratios. We define a population of radio loud AGNs (RLAGNs) based on the presence of ongoing jet feedback. Within RLAGN, we define high and low excitation radio galaxies (HERGs and LERGs) based on their dominant mode of black hole accretion: torque limited accretion representing feeding from a cold disc, or Bondi representing advection-dominated accretion from a hot medium. Simba predicts good agreement with the observed radio luminosity function (RLF) and its evolution, overall as well as separately for HERGs and LERGs. Quiescent galaxies with AGN-dominated radio flux dominate the RLF at $\gtrsim 10^{22-23}$ W Hz−1, while star formation dominates at lower radio powers. Overall, RLAGNs have higher black hole accretion rates and lower star formation rates than non-RLAGN at a given stellar mass or velocity dispersion, but have similar black hole masses. Simba predicts an LERG number density of 8.53 Mpc−3, ∼10× higher than for HERGs, broadly as observed. While LERGs dominate amongmore »most massive galaxies with the largest black holes and HERGs dominate at high specific star formation rates, they otherwise largely populate similar-sized dark matter haloes and have similar host galaxy properties. Simba thus predicts that deeper radio surveys will reveal an increasing overlap between the host galaxy demographics of HERGs and LERGs.« less