Massive Central Galaxies Research Articles

Nearly all Massive Quiescent Disk Galaxies Have a Surprisingly Large Atomic Gas Reservoir

Abstract The massive galaxy population above the characteristic Schechter mass M * ≈ 1010.6 contributes to about half of the total stellar mass in the local universe. These massive galaxies usually reside in hot dark matter halos above the critical shock-heating mass ∼1012 , where the external cold gas supply to these galaxies is expected to be suppressed. When galaxies run out of their cold gas reservoir, they become dead and quiescent. Therefore, massive quiescent galaxies living in hot halos are commonly believed to be gas-poor. Based on the data from SDSS, ALFALFA, GASS, and COLD GASS surveys, here we show that the vast majority of the massive, quiescent, central disk galaxies in the nearby universe have a remarkably large amount of cold atomic hydrogen gas, surprisingly similar to star-forming galaxies. Both star-forming and quiescent disk galaxies show identical symmetric double-horn H i spectra, indicating similar regularly rotating H i disks. Relative to their star-forming counterparts, massive quiescent central disk galaxies are quenched because of their significantly reduced molecular gas content, lower dust content, and lower star formation efficiency. Our findings reveal a new picture, which clearly demonstrates the detailed star formation quenching process in massive galaxies and provides a stringent constraint on the physical mechanism of quenching.

Black hole mass of central galaxies and cluster mass correlation in cosmological hydro-dynamical simulations

Context. The correlations between the properties of the brightest cluster galaxy (BCG) and the mass of its central super-massive black hole (SMBH) have been extensively studied from a theoretical and observational angle. More recently, relations connecting the SMBH mass and global properties of the hosting cluster, such as temperature and mass, were observed. Aims. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical MBH − MBCG relation. Moreover, we study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. Methods. We employed 135 groups and clusters with a mass range 1.4 × 1013 M⊙ − 2.5 × 1015 M⊙ extracted from a set of 29 zoom-in cosmological hydro-dynamical simulations where the baryonic physics is treated with various sub-grid models, including feedback by active galactic nuclei. Results. In our simulations we find that MBH correlates well with M500 and T500, with the scatter around these relations compatible within 2σ with the scatter around MBH − MBCG at z = 0. The MBH − M500 relation evolves with time, becoming shallower at lower redshift as a direct consequence of hierarchical structure formation. On average, in our simulations the contribution of gas accretion to the total SMBH mass dominates for the majority of the cosmic time (z > 0.4), while in the last 2 Gyr the BH-BH mergers become a larger contributor. During this last process, substructures hosting SMBHs are disrupted in the merger process with the BCG and the unbound stars enrich the diffuse stellar component rather than increase BCG mass. Conclusions. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the MBH − T500 relation is comparable to the scatter around the MBH − MBCG relation and that, given the observational difficulties related to the estimation of the BCG mass, clusters temperature and mass can be a useful proxy for the SMBHs mass, especially at high redshift.

3C 294 revisited: Deep Large Binocular Telescope AO NIR images and optical spectroscopy

Context. High redshift radio galaxies are among the most massive galaxies at their redshift, are often found at the center of protoclusters of galaxies, and are expected to evolve into the present day massive central cluster galaxies. Thus they are a useful tool to explore structure formation in the young Universe. Aims. 3C 294 is a powerful FR II type radio galaxy at z = 1.786. Past studies have identified a clumpy structure, possibly indicative of a merging system, as well as tentative evidence that 3C 294 hosts a dual active galactic nucleus (AGN). Due to its proximity to a bright star, it has been subject to various adaptive optics imaging studies. Methods. In order to distinguish between the various scenarios for 3C 294, we performed deep, high-resolution adaptive optics near-infrared imaging and optical spectroscopy of 3C 294 with the Large Binocular Telescope. Results. We resolve the 3C 294 system into three distinct components separated by a few tenths of an arcsecond on our images. One is compact, the other two are extended, and all appear to be non-stellar. The nature of each component is unclear. The two extended components could be a galaxy with an internal absorption feature, a galaxy merger, or two galaxies at different redshifts. We can now uniquely associate the radio source of 3C 294 with one of the extended components. Based on our spectroscopy, we determined a redshift of z = 1.784 ± 0.001, which is similar to the one previously cited. In addition we found a previously unreported emission line at λ6749.4 Å in our spectra. It is not clear that it originates from 3C 294. It could be the Ne [IV] doublet λ2424/2426 Å at z = 1.783, or belong to the compact component at a redshift of z ∼ 4.56. We thus cannot unambiguously determine whether 3C 294 hosts a dual AGN or a projected pair of AGNs.

Kelvin–Helmholtz instability in self-gravitating streams

ABSTRACT Self-gravitating gaseous filaments exist on many astrophysical scales, from sub-pc filaments in the interstellar medium to Mpc scale streams feeding galaxies from the cosmic web. These filaments are often subject to Kelvin–Helmholtz Instability (KHI) due to shearing against a confining background medium. We study the non-linear evolution of KHI in pressure-confined self-gravitating gas streams initially in hydrostatic equilibrium, using analytic models and hydrodynamic simulations, not including radiative cooling. We derive a critical line mass, or mass per unit length, as a function of the stream Mach number and density contrast with respect to the background, μcr(Mb, δc) ≤ 1, where μ = 1 is normalized to the maximal line mass for which initial hydrostatic equilibrium is possible. For μ < μcr, KHI dominates the stream evolution. A turbulent shear layer expands into the background and leads to stream deceleration at a similar rate to the non-gravitating case. However, with gravity, penetration of the shear layer into the stream is halted at roughly half the initial stream radius by stabilizing buoyancy forces, significantly delaying total stream disruption. Streams with μcr < μ ≤ 1 fragment and form round, long-lived clumps by gravitational instability (GI), with typical separations roughly eight times the stream radius, similar to the case without KHI. When KHI is still somewhat effective, these clumps are below the spherical Jeans mass and are partially confined by external pressure, but they approach the Jeans mass as μ → 1 and GI dominates. We discuss potential applications of our results to streams feeding galaxies at high redshift, filaments in the ISM, and streams resulting from tidal disruption of stars near the centres of massive galaxies.

From ‘bathtub’ galaxy evolution models to metallicity gradients

We model gas phase metallicity radial profiles of galaxies in the local Universe by building on the `bathtub' chemical evolution formalism - where a galaxy's gas content is determined by the interplay between inflow, star formation and outflows. In particular, we take into account inside-out disc growth and add physically-motivated prescriptions for radial gradients in star formation efficiency (SFE). We fit analytical models against the metallicity radial profiles of low-redshift star-forming galaxies in the mass range $\log(M_\star/M_\odot)$ = [9.0-11.0] derived by Belfiore et al. 2017, using data from the MaNGA survey. The models provide excellent fits to the data and are capable of reproducing the change in shape of the radial metallicity profiles, including the flattening observed in the centres of massive galaxies. We derive the posterior probability distribution functions for the model parameters and find significant degeneracies between them. The parameters describing the disc assembly timescale are not strongly constrained from the metallicity profiles, while useful constrains are obtained for the SFE (and its radial dependence) and the outflow loading factor. The inferred value for the SFE is in good agreement with observational determinations. The inferred outflow loading factor is found to decrease with stellar mass, going from nearly unity at $\log(M_\star/M_\odot) = 9.0$ to close to zero at $\log(M_\star/M_\odot) =11.0$, in general agreement with previous empirical determinations. These values are the lowest we can obtain for a physically-motivated choice of initial mass function and metallicity calibration. We explore alternative choices which produce larger loading factors at all masses, up to order unity at the high-mass end.

PyUltraLight: a pseudo-spectral solver for ultralight dark matter dynamics

PYULTRALIGHT simulates the dynamics of ultralight dark matter in a non-expanding background. PYULTRALIGHT can describe the evolution of several interacting ultralight dark matter halos or one or more halos orbiting a central, fixed Newtonian potential, the latter scenario corresponding to dwarf galaxies orbiting a massive central galaxy. We verify PYULTRALIGHT by showing that it reproduces qualitative dynamical features of previously published simulations and demonstrate that it has excellent energy-conservation properties. PYULTRALIGHT is implemented in a Python-based Jupyter notebook, solving the Schrödinger-Poisson equation governing ultralight scalar field dark matter dynamics in the non-relativistic regime using a symmetrised split-step pseudospectral algorithm. The notebook interface makes it simple to specify simulation parameters and visualise the resulting output but performance-critical routines are managed via calls to computationally efficient compiled libraries. PYULTRALIGHT runs on standard desktop hardware with support for shared memory mutlithreading and is available on GitHub.

Barlenses in the CALIFA survey: Combining photometric and stellar population analyses

Aims: It is theoretically predicted that, at low galaxy inclinations, boxy/peanut bar components have a barlens appearance of a round central component embedded in the narrow bar. We investigate barlenses in the Calar Alto Legacy Integral Field Area (CALIFA) survey galaxies, studying their morphologies, stellar populations, and metallicities. We show that, when present, barlenses account for a significant portion of light of photometric bulges, i.e., the excess light on top of the disks, which highlights the importance of bars in accumulating central galaxy mass concentrations in the cosmic timescale. Methods: We made multi-component decompositions for a sample of 46 barlens galaxies drawn from the CALIFA survey, where M⋆/M⊙ = 109.7 − 1011.4 and z = 0.005 − 0.03. Unsharp masks of the Sloan Digital Sky Survey (SDSS) r′-band mosaics were used to identify the boxy/peanut or X-shaped features. Barlenses are identified in the images using our simulation snapshots as an additional guide. Our decompositions with GALFIT include bulges, disks, and bars as well as barlenses as a separate component. For 26 of the decomposed galaxies the CALIFA DR2 V500 grating data cubes were used to explore stellar ages and metallicities at the regions of various structure components. Results: We find that 25 ± 2% of the 1064 galaxies in the whole CALIFA sample show either X-shaped or barlens features. In the decomposed galaxies with barlenses, on average 13% ± 2% of the total galaxy light belongs to this component, leaving less than 10% for possible separate bulge components. Most importantly, bars and barlenses are found to have similar cumulative stellar age and metallicity distributions. The metallicities in barlenses are on average near solar, but exhibit a large range. In some of the galaxies barlenses and X-shaped features appear simultaneously, in which case the bar origin of the barlens is unambiguous. Conclusion: This is the first time that a combined morphological and stellar population analysis is used to study barlenses. We show that their stars are accumulated in a prolonged time period concurrently with the evolution of the narrow bar.

The impact of assembly bias on the halo occupation in hydrodynamical simulations

We investigate the variations in galaxy occupancy of the dark matter haloes with the large-scale environment and halo formation time, using two state-of-the-art hydrodynamical cosmological simulations, EAGLE and Illustris. For both simulations, we use three galaxy samples with a fixed number density ranked by stellar mass. For these samples we find that low-mass haloes in the most dense environments are more likely to host a central galaxy than those in the least dense environments. When splitting the halo population by formation time, these relations are stronger. Hence, at a fixed low halo mass, early-formed haloes are more likely to host a central galaxy than late-formed haloes since they have had more time to assemble. The satellite occupation shows a reverse trend where early-formed haloes host fewer satellites due to having more time to merge with the central galaxy. We also analyse the stellar mass -- halo mass relation for central galaxies in terms of the large-scale environment and formation time of the haloes. We find that low mass haloes in the most dense environment host relatively more massive central galaxies. This trend is also found when splitting the halo population by age, with early-formed haloes hosting more massive galaxies. Our results are in agreement with previous findings from semi-analytical models, providing robust predictions for the occupancy variation signature in the halo occupation distribution of galaxy formation models.

Pushing back the limits: detailed properties of dwarf galaxies in a ΛCDM universe

We present the results of a set of high-resolution chemo-dynamical simulations of dwarf galaxies in a ΛCDM cosmology. Out of an original (3.4 Mpc/h)3 cosmological box, a sample of 27 systems are re-simulated from z = 70 to z = 0 using a zoom-in technique. Gas and stellar properties are confronted to the observations in the greatest details: in addition to the galaxy global properties, we investigated the model galaxy velocity dispersion profiles, half-light radii, star formation histories, stellar metallicity distributions, and [Mg/Fe] abundance ratios. The formation and sustainability of the metallicity gradients and kinematically distinct stellar populations are also tackled. We show how the properties of six Local Group dwarf galaxies, NGC 6622, Andromeda II, Sculptor, Sextans, Ursa Minor and Draco are reproduced, and how they pertain to three main galaxy build-up modes. Our results indicate that the interaction with a massive central galaxy could be needed for a handful of Local Group dwarf spheroidal galaxies only, the vast majority of the systems and their variety of star formation histories arising naturally from a ΛCDM framework. We find that models fitting well the local Group dwarf galaxies are embedded in dark haloes of mass between 5 × 108 to a few 109 M⊙, without any missing satellite problem. We confirm the failure of the abundance matching approach at the mass scale of dwarf galaxies. Some of the observed faint however gas-rich galaxies with residual star formation, such as Leo T and Leo P, remain challenging. They point out the need of a better understanding of the UV-background heating.

A detection of the environmental dependence of the sizes and stellar haloes of massive central galaxies

We use ~100 square deg of deep (>28.5 mag arcsec$^{-2}$ in i-band), high-quality (median 0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) survey to reveal the halo mass dependence of the surface mass density profiles and outer stellar envelopes of massive galaxies. The i-band images from the HSC survey reach ~4 magnitudes deeper than Sloan Digital Sky Survey and enable us to directly trace stellar mass distributions to 100 kpc without requiring stacking. We conclusively show that, at fixed stellar mass, the stellar profiles of massive galaxies depend on the masses of their dark matter haloes. On average, massive central galaxies with $\log M_{\star, 100\ \mathrm{kpc}}>11.6$ in more massive haloes at 0.3 < z < 0.5 have shallower inner stellar mass density profiles (within ~10-20 kpc) and more prominent outer envelopes. These differences translate into a halo mass dependence of the mass-size relation. Central galaxies in haloes with $\log M_{\rm{Halo}}>14.0$ are ~20% larger in $R_{\mathrm{50}}$ at fixed stellar mass. Such dependence is also reflected in the relationship between the stellar mass within 10 and 100 kpc. Comparing to the mass--size relation, the $\log M_{\star, 100\ \rm{kpc}}$-$\log M_{\star, 10\ \rm{kpc}}$ relation avoids the ambiguity in the definition of size, and can be straightforwardly compared with simulations. Our results demonstrate that, with deep images from HSC, we can quantify the connection between halo mass and the outer stellar halo, which may provide new constraints on the formation and assembly of massive central galaxies.

Supermassive black holes and their feedback effects in the IllustrisTNG simulation

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.

Study of central light concentration in nearby galaxies

We propose a novel technique to estimate the masses of super massive black holes (SMBHs) residing at the centres of massive galaxies in the nearby Universe using simple photometry. Aperture photometry using SEXTRACTOR is employed to determine the central intensity ratio (CIR) at the optical centre of the galaxy image for a sample of 49 nearby galaxies with SMBH mass estimations. We find that the CIR of ellipticals and classical bulges is strongly correlated with SMBH masses whereas pseudo bulges and ongoing mergers show significant scatter. Also, the CIR of low luminosity AGNs in the sample shows significant connection with the 5 GHz nuclear radio emission suggesting a stronger link between the former and the SMBH evolution in these galaxies. In addition, it is seen that various structural and dynamical properties of the SMBH host galaxies are correlated with the CIR making the latter an important parameter in galaxy evolution studies. Finally, we propose the CIR to be an efficient and simple tool not only to distinguish classical bulges from pseudo bulges but also to estimate the mass of the central SMBH.

First results from the IllustrisTNG simulations: the stellar mass content of groups and clusters of galaxies

The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the Arepo code and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light. The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( $-3.5 < \alpha_{\rm 3D} < -3$). Total halo mass is a very good predictor of stellar mass, and vice versa: at $z=0$, the 3D stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of TNG less massive galaxies ($< 10^{11}M_{\rm sun}$ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, with $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.

Millimetre-wave emission from an intermediate-mass black hole candidate in the Milky Way

It is widely accepted that black holes (BHs) with masses greater than a million solar masses (Msun) lurk at the centres of massive galaxies. The origins of such `supermassive' black holes (SMBHs) remain unknown (Djorgovski et al. 1999), while those of stellar-mass BHs are well-understood. One possible scenario is that intermediate-mass black holes (IMBHs), which are formed by the runaway coalescence of stars in young compact star clusters (Portagies Zwart et al. 1999), merge at the centre of a galaxy to form an SMBH (Ebisuzaki et al. 2001). Although many candidates for IMBHs have been proposed to date, none of them are accepted as definitive. Recently we discovered a peculiar molecular cloud, CO-0.40-0.22, with an extremely broad velocity width near the centre of our Milky Way galaxy. Based on the careful analysis of gas kinematics, we concluded that a compact object with a mass of ~1E5 Msun is lurking in this cloud (Oka et al. 2016). Here we report the detection of a point-like continuum source as well as a compact gas clump near the center of CO-0.40-0.22. This point-like continuum source (CO-0.40-0.22*) has a wide-band spectrum consistent with 1/500 of the Galactic SMBH (Sgr A*) in luminosity. Numerical simulations around a point-like massive object reproduce the kinematics of dense molecular gas well, which suggests that CO-0.40-0.22* is the most promising candidate for an intermediate-mass black hole.

Are Fossil Groups Early-forming Galaxy Systems?

Abstract Using the Illustris cosmological simulation, we investigate the origin of fossil groups in the mass regime. We examine the formation of the two primary features of fossil groups: the large magnitude gap between their two brightest galaxies and their exceptionally luminous brightest group galaxy (BGG). For fossils and nonfossils identified at z = 0, we find no difference in their halo mass assembly histories at early times, departing from previous studies. However, we do find a significant difference in the recent accretion history of fossil and nonfossil halos; in particular, fossil groups show a lack of recent accretion and have in majority assembled 80% of their mass before . For fossils, massive satellite galaxies accreted during this period have enough time to merge with the BGG by the present day, producing a more massive central galaxy. In addition, the lack of recent group accretion prevents replenishment of the bright satellite population, allowing for a large magnitude gap to develop within the past few Gyr. We thus find that the origin of the magnitude gap and overmassive BGG of fossils in Illustris depends on the recent accretion history of the groups and merger history of the BGGs after their collapse at . This indicates that selecting galaxy groups by their magnitude gap does not guarantee obtaining either early-forming galaxy systems or undisturbed central galaxies.

Physics of Galactic Metals: Evolutionary Effects due to Production, Distribution, Feedback, and Interaction with Black Holes

Abstract We ask how the inclusion of various physical heating processes due to the metal content of gas affects the evolution of central massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion, and cooling rate based on individual element abundances. The feedback models include (1) active galactic nucleus (AGN) feedback via high-velocity broad absorption line winds and Compton/photoionization heating; (2) stellar feedback from multiple processes, including powerful winds from supernovae, stellar winds from young massive stars, and AGB stars, as well as radiative heating within Str mgren spheres; and (3) additional heating effects due to the presence of metals, including grain photoelectric heating and metallicity-dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background. With a suite of zoom-in simulations of 30 halos with , we show that energy and momentum budgeted from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal-induced heating reduces the fraction of accreted stellar material but overall has a relatively minor effect on the massive central galaxies. The inclusion of AGN feedback significantly improves the ability of our simulations to yield realistic gas and stellar properties of massive galaxies with a reasonable accreted star fraction from other galaxies.

CANDELS Sheds Light on the Environmental Quenching of Low-mass Galaxies

Abstract We investigate the environmental quenching of galaxies, especially those with stellar masses (M *) &lt; 109.5 M ⊙, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test whether or not such a dwarf QG-massive central galaxy connection exists beyond the local universe. For this purpose, we only need a statistically representative, rather than complete, sample of low-mass galaxies, which enables our study to z ≳ 1.5. For each low-mass galaxy, we measure the projected distance (d proj) to its nearest massive neighbor (M * &gt; 1010.5 M ⊙) within a redshift range. At a given z and M *, the environmental quenching effect is considered to be observed if the d proj distribution of QGs ( ) is significantly skewed toward lower values than that of star-forming galaxies ( ). For galaxies with 108 M ⊙ &lt; M * &lt; 1010 M ⊙, such a difference between and is detected up to z ∼ 1. Also, about 10% of the quenched galaxies in our sample are located between two and four virial radii (R Vir) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, , decreases with satellite M * at M * ≲ 109.5 M ⊙, but increases with satellite M * at M * ≳ 109.5 M ⊙. This trend suggests a smooth, if any, transition of the quenching timescale around M * ∼ 109.5 M ⊙ at 0.5 &lt; z &lt; 1.0.

Star Formation Quenching Timescale of Central Galaxies in a Hierarchical Universe

Abstract Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at . Thus, the mechanism(s) responsible for quenching central galaxies play a crucial role in galaxy evolution as whole. We combine a high-resolution cosmological N-body simulation with observed evolutionary trends of the “star formation main sequence,” quiescent fraction, and stellar mass function at to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range we infer quenching e-folding times that span 1.5–0.5 Gyr with more massive central galaxies quenching faster. For , this implies a total migration time of from the star formation main sequence to quiescence. Compared to satellites, central galaxies take longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this remain unclear.

Colours, star formation rates and environments of star-forming and quiescent galaxies at the cosmic noon

We analyse the star formation rates (SFRs), colours, and dust extinctions of galaxies in massive (10^12.5-10^13.5 Msun) haloes at z~2 in high-resolution, cosmological zoom-in simulations as part of the Feedback in Realistic Environments (FIRE) project. The simulations do not model feedback from active galactic nuclei (AGN) but reproduce well the observed relations between stellar and halo mass and between stellar mass and SFR. About half (a third) of the simulated massive galaxies (massive central galaxies) at z~2 have broad-band colours classifying them as 'quiescent', and the fraction of quiescent centrals is steeply decreasing towards higher redshift, in agreement with observations. The progenitors of z~2 quiescent central galaxies are, on average, more massive, have lower specific SFRs, and reside in more massive haloes than the progenitors of similarly massive star forming centrals. The simulations further predict a morphological mix of galaxies that includes disk-dominated, irregular, and early-type galaxies. However, our simulations do not reproduce the reddest of the quiescent galaxies observed at z~2. We also do not find evidence for a colour bimodality, but are limited by our modest sample size. In our simulations, the star formation activity of central galaxies of moderate mass (Mstar~10^10-10^11 Msun) is affected by a combination of two distinct physical processes. Outflows powered by stellar feedback result in a short-lived (<100 Myr), but almost complete, suppression of star formation activity after which many galaxies quickly recover and continue to form stars at normal rates. In addition, galaxies residing in slowly growing haloes tend to experience a moderate reduction of their SFRs ('cosmological starvation'). The relative importance of these processes and AGN feedback is uncertain and will be explored in future work.

The H α luminosity-dependent clustering of star-forming galaxies from z ∼ 0.8 to ∼2.2 with HiZELS

We present clustering analyses of identically-selected star-forming galaxies in 3 narrow redshift slices (at z=0.8, z=1.47 and z=2.23), from HiZELS, a deep, near-infrared narrow-band survey. The HiZELS samples span the peak in the cosmic star-formation rate density, identifying typical star-forming galaxies at each epoch. Narrow-band samples have well-defined redshift distributions and are therefore ideal for clustering analyses. We quantify the clustering of the three samples, and of H-alpha luminosity-selected subsamples, initially using simple power law fits to the two-point correlation function. We extend this work to link the evolution of star-forming galaxies and their host dark matter halos over cosmic time using sophisticated dark matter halo models. We find that the clustering strength, r0, and the bias of galaxy populations relative to the clustering of dark matter increase linearly with H-alpha luminosity (and, by implication, star-formation rate) at all three redshifts, as do the host dark matter halo masses of the HiZELS galaxies. The typical galaxies in our samples are star-forming centrals, residing in halos of mass M_halo ~ a few times 10^12M_solar. We find a remarkably tight redshift-independent relation between the H-alpha luminosity scaled by the characteristic luminosity, L(H-alpha)/L(H-alpha)*(z), and the minimum host dark matter halo mass of central galaxies. This reveals that the dark matter halo environment is a strong driver of galaxy star-formation rate and therefore of the evolution of the star-formation rate density in the Universe.