We use a large suite of carefully controlled full hydrodynamic simulations to study the ram pressure stripping of the hot gaseous haloes of galaxies as they fall into massive groups and clusters. The ...sensitivity of the results to the orbit, total galaxy mass, and galaxy structural properties is explored. For typical structural and orbital parameters, we find that ∼30 per cent of the initial hot galactic halo gas can remain in place after 10 Gyr. We propose a physically simple analytic model that describes the stripping seen in the simulations remarkably well. The model is analogous to the original formulation of Gunn & Gott, except that it is appropriate for the case of a spherical (hot) gas distribution (as opposed to a face-on cold disc) and takes into account that stripping is not instantaneous but occurs on a characteristic time-scale. The model reproduces the results of the simulations to within ≈10 per cent at almost all times for all the orbits, mass ratios, and galaxy structural properties we have explored. The one exception involves unlikely systems where the orbit of the galaxy is highly non-radial and its mass exceeds about 10 per cent of the group or cluster into which it is falling (in which case the model underpredicts the stripping following pericentric passage). The proposed model has several interesting applications, including modelling the ram pressure stripping of both observed and cosmologically simulated galaxies and as a way to improve present semi-analytic models of galaxy formation. One immediate consequence is that the colours and morphologies of satellite galaxies in groups and clusters will differ significantly from those predicted with the standard assumption of complete stripping of the hot coronae.
Abstract
We present a MUSE (Multi-Unit Spectroscopic Explorer) and KMOS (K-band Multi-Object Spectrograph) dynamical study 405 star-forming galaxies
at redshift z = 0.28–1.65 (median redshift
...$\bar{z}$
= 0.84). Our sample is representative of the star-forming ‘main sequence’, with star formation rates of SFR = 0.1–30 M⊙ yr−1 and stellar masses M
⋆ = 108–1011 M⊙. For 49 ± 4 per cent of our sample, the dynamics suggest rotational support, 24 ± 3 per cent are unresolved systems and 5 ± 2 per cent appear to be early-stage major mergers with components on 8–30 kpc scales. The remaining 22 ± 5 per cent appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination-corrected rotational velocities and show that these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j
⋆ = J/
$M_\star \propto M_\star ^{2/3}$
but with a redshift evolution that scales as j
⋆ ∝ M
$_\star ^{2/3}(1+z)^{-1}$
. We also identify a correlation between specific angular momentum and disc stability such that galaxies with the highest specific angular momentum (log(j
⋆/M
$_\star ^{2/3}$
) > 2.5) are the most stable, with Toomre Q = 1.10 ± 0.18, compared to Q = 0.53 ± 0.22 for galaxies with log(j
⋆/M
$_\star ^{2/3}$
) < 2.5. At a fixed mass, the Hubble Space Telescope morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas discs: at z ∼ 1, massive galaxies that have discs with low specific angular momentum are globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular momentum have evolved into stable discs with spiral structure where star formation is a local (rather than global) process.
Breaking the hierarchy of galaxy formation Bower, R. G.; Benson, A. J.; Malbon, R. ...
Monthly notices of the Royal Astronomical Society,
August 2006, Letnik:
370, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Recent observations of the distant Universe suggest that much of the stellar mass of bright galaxies was already in place at z > 1. This presents a challenge for models of galaxy formation because ...massive haloes are assembled late in the hierarchical clustering process intrinsic to the cold dark matter (CDM) cosmology. In this paper, we discuss a new implementation of the Durham semi-analytic model of galaxy formation in which feedback due to active galactic nuclei (AGN) is assumed to quench cooling flows in massive haloes. This mechanism naturally creates a break in the local galaxy luminosity function at bright magnitudes. The model is implemented within the Millennium N-body simulation. The accurate dark matter merger trees and large number of realizations of the galaxy formation process enabled by this simulation result in highly accurate statistics. After adjusting the values of the physical parameters in the model by reference to the properties of the local galaxy population, we investigate the evolution of the K-band luminosity and galaxy stellar mass functions. We calculate the volume-averaged star formation rate density of the Universe as a function of redshift and the way in which this is apportioned amongst galaxies of different mass. The model robustly predicts a substantial population of massive galaxies out to redshift z∼ 5 and a star formation rate density which rises at least out to z∼ 2 in objects of all masses. Although observational data on these properties have been cited as evidence for 'antihierarchical' galaxy formation, we find that when AGN feedback is taken into account, the fundamentally hierarchical CDM model provides a very good match to these observations.
Abstract
We present sensitive 870 μm continuum measurements from our ALMA programmes of 114 X-ray selected active galactic nuclei (AGN) in the Chandra Deep Field-South and Cosmic Evolution Survey ...fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at z = 1.5–3.2 with stellar mass >2 × 1010 M⊙. We constructed broad-band spectral energy distributions in the infrared band (8–1000 μm) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above ≈1010 M⊙.
The case for AGN feedback in galaxy groups McCarthy, I. G.; Schaye, J.; Ponman, T. J. ...
Monthly notices of the Royal Astronomical Society,
August 2010, Letnik:
406, Številka:
2
Journal Article
Recenzirano
Odprti dostop
The relatively recent insight that energy input from supermassive black holes (BHs) can have a substantial effect on the star formation rates (SFRs) of galaxies motivates us to examine the effects of ...BH feedback on the scale of galaxy groups. At present, groups contain most of the galaxies and a significant fraction of the overall baryon content of the Universe and, along with massive clusters, they represent the only systems for which it is possible to measure both the stellar and gaseous baryonic components directly. To explore the effects of BH feedback on groups, we analyse two high-resolution cosmological hydrodynamic simulations from the OverWhelmingly Large Simulations (OWLS) project. While both include galactic winds driven by supernovae, only one of the models includes feedback from accreting BHs. We compare the properties of the simulated galaxy groups to a wide range of observational data, including the entropy and temperature profiles of the intragroup medium, hot gas mass fractions, the luminosity–temperature and mass–temperature scaling relations, the K-band luminosity of the group and its central brightest galaxy (CBG), SFRs and ages of the CBG, and gas- and stellar-phase metallicities. Both runs yield entropy distributions similar to the data, while the run without active galactic nucleus (AGN) feedback yields highly peaked temperature profiles, in discord with the observations. Energy input from supermassive BHs significantly reduces the gas mass fractions of galaxy groups with masses less than a few × 1014 M⊙, yielding a gas mass fraction and X-ray luminosity scaling with system temperature that is in excellent agreement with the data, although the detailed scatter in the L–T relation is not quite correct. The run without AGN feedback suffers from the well-known overcooling problem – the resulting stellar mass fractions are several times larger than observed and present-day cooling flows operate uninhibitedly. By contrast, the run that includes BH feedback yields stellar mass fractions, SFRs and stellar age distributions in excellent agreement with current estimates, thus resolving the long-standing ‘cooling crisis’ of simulations on the scale of groups. Both runs yield very similar gas-phase metal abundance profiles that match X-ray measurements, but they predict very different stellar metallicities. Based on the above, galaxy groups provide a compelling case that feedback from supermassive BHs is a crucial ingredient in the formation of massive galaxies.
The diffuse plasma that fills galaxy groups and clusters (the intracluster medium) is a by-product of galaxy formation. The present thermal state of this gas results from a competition between gas ...cooling and heating. The heating comes from two distinct sources: gravitational heating associated with the collapse of the dark matter halo and additional thermal input from the formation of galaxies and their black holes. A long-term goal of this research is to decode the observed temperature, density and entropy profiles of clusters and to understand the relative roles of these processes. However, a long-standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh-based codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near power-law entropy profiles, while those formed in Eulerian simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas distribution, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealized merger of two clusters and that the differences are not the result of the lower effective resolution of Eulerian cosmological simulations. The difference is not sensitive to the minimum mesh size (in Eulerian codes) or the number of particles used (in SPH codes). We investigate whether the difference is the result of the different gravity solvers, the Galilean non-invariance of the mesh code or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of vortices in the two codes. Particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. The origin of this difference lies in the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the Eulerian simulations that very efficiently mix the fluid and erase the low-entropy gas. We discuss the potentially profound implications of these results.
Galaxy groups are not scaled down versions of massive galaxy clusters - the hot gas in groups known as the intragroup medium (IGrM) is, on average, less dense than the intracluster medium, implying ...that one or more non-gravitational processes (e.g. radiative cooling, star formation and/or feedback) has had a relatively larger effect on groups. In the present study, we compare a number of cosmological hydrodynamic simulations that form part of the OverWhelmingly Large Simulations project to isolate and quantify the effects of cooling and feedback from supernovae (SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing Lagrangian thermal histories of the gas in the different runs, which were all started from identical initial conditions. While radiative cooling, star formation and SN feedback are all necessary ingredients, only runs that also include AGN feedback are able to successfully reproduce the optical and X-ray properties of groups and low-mass clusters. We isolate how, when and exactly what gas is heated by AGN. Interestingly, we find that the gas that constitutes the present-day IGrM is that which was not strongly heated by AGN. Instead, the low median density/high median entropy of the gas in present-day groups is achieved by the ejection of lower entropy gas from low-mass progenitor galaxies at high redshift (primarily 2 ≲z≲ 4). This corresponds to the epoch when supermassive black holes accreted most of their mass, typically at a rate that is close to the Eddington limit (i.e. when the black holes are in a 'quasar mode').
Only ∼10 per cent of baryons in the Universe are in the form of stars, yet most models of luminous structure formation have concentrated on the properties of the luminous stellar matter. Such models ...are now largely successful at reproducing the observed properties of galaxies, including the galaxy luminosity function and the star formation history of the universe. In this paper we focus on the ‘flip side’ of galaxy formation and investigate the properties of the material that is not presently locked up in galaxies. This ‘by-product’ of galaxy formation can be observed as an X-ray emitting plasma the intracluster medium (ICM) in groups and clusters. Since much of this material has been processed through galaxies, observations of the ICM represent an orthogonal set of constraints on galaxy formation models. In this paper, we attempt to self-consistently model the formation of galaxies and the heating of the ICM. We set out the challenges for such a combined model and demonstrate a possible means of bringing the model into line with both sets of constraints. In this paper, we present a version of the Durham semi-analytic galaxy formation model galform that allows us to investigate the properties of the ICM. As we would expect on the basis of gravitational scaling arguments, the previous model fails to reproduce even the most basic observed properties of the ICM. We present a simple modification to the model to allow for heat input into the ICM from the active galactic nucleus (AGN) ‘radio-mode’ feedback. This heating acts to expel gas from the X-ray luminous central regions of the host halo. With this modification, the model reproduces the observed gas mass fractions and luminosity–temperature (L–T) relation of groups and clusters. In contrast to simple ‘pre-heating’ models of the ICM, the model predicts mildly positive evolution of the L–T relation, particularly at low temperatures. The model is energetically plausible, but seems to exceed the observed heating rates of intermediate-temperature clusters. Introducing the heating process into the model requires changes to a number of model parameters in order to retain a good match to the observed galaxy properties. With the revised parameters, the best-fitting luminosity function is comparable to that presented in Bower et al. The new model makes a fundamental step forward, providing a unified model of galaxy and cluster ICM formation. However, the detailed comparison with the data is not completely satisfactory, and we highlight key areas for improvement.
We analyse a z < 0.1 galaxy sample from the Sloan Digital Sky Survey focusing on the variation in the galaxy colour bimodality with stellar mass and projected neighbour density Σ, and on measurements ...of the galaxy stellar mass functions. The characteristic mass increases with environmental density from about 1010.6 to (Kroupa initial mass function, H0= 70) for Σ in the range 0.1–10 Mpc−2. The galaxy population naturally divides into a red and blue sequence with the locus of the sequences in colour–mass and colour–concentration indices not varying strongly with environment. The fraction of galaxies on the red sequence is determined in bins of 0.2 in log Σ and bins). The red fraction fr generally increases continuously in both Σ and such that there is a unified relation: . Two simple functions are proposed which provide good fits to the data. These data are compared with analogous quantities in semi-analytical models based on the Millennium N-body simulation: the Bower et al. and Croton et al. models that incorporate active galactic nucleus feedback. Both models predict a strong dependence of the red fraction on stellar mass and environment that is qualitatively similar to the observations. However, a quantitative comparison shows that the Bower et al. model is a significantly better match; this appears to be due to the different treatment of feedback in central galaxies.