The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end
Galaxies fall into two clearly distinct types: ‘blue-sequence’ galaxies that are rapidly forming young stars, and ‘red-sequence’ galaxies in which star formation has almost completely ceased. Most galaxies more massive than follow the red-sequence while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy’s centre. We develop a simple analytic model for this interaction. In galaxies less massive than young stars and supernovae drive a high entropy outflow which is more buoyant than any tenuous corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by an increasingly hot corona. Above a halo mass of , the outflow ceases to be buoyant and star formation is unable to prevent the build up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy’s corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations. So long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice.
keywords:black hole physics, galaxies: formation, galaxies: active, methods: hydrodynamical simulations, quasars: general.
Galaxies fall into two clearly distinct types: active, ‘blue-sequence’ galaxies that are rapidly forming young stars, and passive ‘red-sequence’ galaxies in which star formation has almost completely ceased. The two sequences are clearly seen when galaxy colours or star formation rates are plotted as a function of galaxy mass (eg., Kauffmann et al. (2003); Baldry et al. (2006)). Low-mass galaxies generally follow the ‘blue-sequence’ with a tight, almost linear, relationship between star formation rate and stellar mass (eg., Brinchmann et al., 2004), while massive galaxies follow the ‘red-sequence’ with almost undetectable levels of star formation (eg., Bower et al., 1992). At the present day, the transition between the two types occurs at a stellar mass scale of . Galaxies less massive than the transition-scale grow through star formation, doubling their stellar mass on a timescale comparable to the age of the Universe; above the transition mass, galaxy growth slows and is driven primarily by galaxy mergers (eg., De Lucia et al., 2006; Parry et al., 2009; Qu et al., 2016; Rodriguez-Gomez et al., 2016). The existence of the transition mass is closely related to the form of the galaxy stellar mass function, creating the exponential break at high masses (eg., Benson et al., 2003; Peng et al., 2010). The transition mass is sometimes referred to as the ‘quenching’ mass scale. In this paper we will focus on the properties of central galaxies (we will not consider the environmental effects that may suppress star formation in lower mass satellite galaxies, see Trayford et al. (2016)) and show that the transition mass scale arises from a competition between star formation driven outflows and black hole accretion.
Most of the stars in the Universe today were, however, formed when the Universe was less than half its present age. Recently, deep redshift surveys have been able to convincingly demonstrate the existence of a transition mass at higher redshifts (eg., Peng et al., 2010; Ilbert et al., 2015; Darvish et al., 2016). It is useful to illustrate the balance of galaxies on the two sequences by revisiting the analysis of Kauffmann et al. (2003). This is illustrated in Fig. 1, where contours show the star formation rate growth timescale of galaxies at using observational data from the COSMOS high-redshift galaxy survey (Ilbert et al., 2015). So that passive galaxies appear in the figure, we assign galaxies without detectable star formation a growth timescale of with a scatter of 0.2 dex. The relative contributions of blue-and red-sequence galaxies are then computed using the luminosity functions presented by Ilbert et al. (2013). The figure clearly shows that the division into sequences seen in present-day galaxies was already established 8 billion years ago and that the transition mass scale of has changed little over the intervening time.
A simple way to begin to understand galaxy formation is to view galaxies as equilibrium systems in which the star formation rate must balance the gas inflow rate, either by converting the inflowing gas into stars or, more importantly, by driving an outflow (eg., White & Frenk, 1991; Finlator & Davé, 2008; Schaye et al., 2010; Bouché et al., 2010; Davé et al., 2012). Such a model broadly explains many aspects of galaxy evolution, such as the almost linear correlation between stellar mass and star formation rate, and the rate of evolution of this sequence. In order to explain the flat faint-end slope of the galaxy mass function, such models require that the mass loading of the outflow depends strongly on galaxy mass. Low-mass galaxies lose most of their mass in the outflow (and form few stars) while galaxies around the transition mass scale consume much of the inflowing mass in star formation.
Such models do not generally consider the role of the nuclear supermassive black hole, however. Yet the energy liberated when one solar mass of gas is accreted by a black hole is 10,000 times the supernova energy released by forming the same mass of stars. Observational measurements of energetic outflows and radio jets from galaxy nuclei (eg., Harrison et al., 2012; Maiolino et al., 2012; Fabian et al., 2003) indeed suggest that black holes play an important role in galaxy formation, most likely by heating the surrounding gas corona, offsetting cooling losses and disrupting the gas inflow (eg., Binney & Tabor, 1995; Silk & Rees, 1998; Dubois et al., 2013). Such observations have motivated the inclusion of black hole feedback in cosmological galaxy formation models (eg., Di Matteo et al., 2005; Croton et al., 2006; Bower et al., 2006; Sijacki et al., 2007; Booth & Schaye, 2009; Dubois et al., 2013; Sijacki et al., 2015; Dubois et al., 2016).
The success of these models results from including two different modes of black hole feedback or accretion depending on halo mass or Eddington accretion ratio. In the most extreme implementation, the effect of black hole feedback is implemented by switching cooling off in massive haloes (Gabor & Davé, 2015). There are three possible arguments that may be advanced to support a halo mass or accretion rate dependence: (i) (in the absence of disk instabilities) black holes are only able to accrete efficiently if the surrounding gas is hot and pressure supported (Croton et al., 2006); (ii) black hole feedback is only effective if the black hole is surrounded by a hot halo that is able to capture the energy of the black hole jet (Bower et al., 2006); (iii) only black holes accreting well below the Eddington limit (perhaps in the ADAF regime) are able to generate mechanical outflows (Meier, 1999; Nemmen et al., 2007). In the latter case, a halo mass dependence of the feedback efficiency emerges because of the strong dependence of black hole mass on halo mass. In each case, however, the models provide only a qualitative explanation of the transition mass scale. The first case assumes angular momentum prevents accretion of cold gas; the second argument does not explain why black holes do not undergo run-away growth in low-mass haloes; the third links galaxy properties to the highly uncertain AU-scale physics of black hole accretion disks.
In contrast to models that include an explicit mass or accretion rate dependence, the MassiveBlackII (Khandai et al., 2015) and EAGLE (Schaye et al., 2015) simulations adopt a simpler description in which AGN feedback power is a fixed fraction of the rest mass accretion rate. Despite this simple proportionality, the EAGLE simulations reproduce the observed galaxy transition mass scale well. For example, the star formation growth time scales of EAGLE galaxies, shown by the coloured points in Fig 1, match the observation data well (for further discussion, see Furlong et al. (2015b)). Although a transition mass is not so clearly evident in the MassiveBlackII simulations, the results are broadly consistent with a variation of the EAGLE simulation in which star formation driven outflow are weak (see §3.4).
The success of the simple black hole feedback scheme in EAGLE suggests that the transition mass scale emerges from the interaction between star formation feedback and black hole fueling. This has motivated us to explore a scenario in which star formation-driven outflows themselves regulate the density of gas reaching the black hole. This has previously been considered in the context of high resolution simulations of individual high redshift galaxies (eg., Dubois et al., 2015; Habouzit et al., 2016). In low-mass galaxies, such outflows are efficient, making it difficult for high gas densities to build up in the central regions of the galaxy. The critical and novel component of our model is the recognition that outflows from more massive galaxies are not buoyant. As the mass of the galaxy’s dark matter halo increases, a hot corona begins to form (eg., White & Frenk, 1991; Birnboim & Dekel, 2003) and the outflow becomes trapped. Star formation can then no longer prevent the gas density increasing in the central regions of the galaxy, triggering a strongly non-linear response from the black hole, heating the corona and disrupting the inflow of fresh fuel for star formation. In this picture, the falling effectiveness of feedback from star formation leads to an increase in accretion on to the black hole. Black hole feedback will take over just as star formation driven feedback fails.
The structure of the paper is as follows. In §2, we develop a simple analytic model that captures the key physics of the problem. We begin by exploring the strongly non-linear growth rates of black holes accreting from a constant density medium. We show that in the Bondi accretion regime black holes initially grow slowly and then abruptly switch to a rapid accretion phase. In the absence of feedback, the black hole grows to infinite mass in a finite time. We go on to consider how the density in the central parts of the galaxy will evolve as the halo mass grows, and develop a model in which the gas density is determined by the buoyancy of the star formation driven outflow. We show that a critical mass scale emerges: below this scale galaxies balance gas inflow with star formation driven outflows, but above this mass the galaxy and its gas corona are regulated by black hole accretion and star formation is strongly suppressed. We use this simple model to explore the growth of black holes in a cosmological context. In §3, we validate the analytic model by comparing it to galaxies forming in the EAGLE hydrodynamic simulation suite. We show that the approximations of the analytic model are well supported, and explore the effects of varying the sub-grid parameters of the simulations. We show that the galaxy transition mass scale is robust to the choice of subgrid parameters, but that it disappears entirely if star formation driven outflow are absent. We present a summary of the results and a discussion of the model’s wider implications in §4.
2 An Analytic Model For the Growth of Black Holes
2.1 Black hole growth in a constant density medium
We begin by considering the growth of a black hole in a constant density medium. Black holes accreting at their theoretical maximum rate (the Eddington rate) grow exponentially with a timescale, (since ) that is much shorter than the present-age of the Universe (Salpeter, 1964) . In practice, however, the accretion rates of black holes usually lie below this rate because the gas density is low. If the gas surrounding the black hole has density and effective sound speed (including turbulent pressure of the interstellar medium), , the black hole accretion disk is fed at a rate (Bondi, 1952):
(where is Newton’s gravitational constant). The relevant scale on which the density is measured is given by the Bondi radius, corresponding to around 500 pc for a black hole of mass accreting from the diffuse interstellar medium (assuming an effective sound speed of ). We have inserted a factor, , to account for the suppression of accretion by the gas motion relative to the black hole and angular momentum (Bondi & Hoyle, 1944); on average, in the EAGLE simulations irrespective of the system mass (Rosas-Guevara et al., 2015). We will later validate our analytic model by comparing to the EAGLE simulations that assume star forming gas follows an effective equation of state when averaged over the kpc-scale areas of the galaxy. In this case, . Assuming, instead, an isothermal equation of state (ie., constant) would result in a stronger density dependence, , which would strengthen our conclusions.
In the Bondi regime, black hole accretion is highly non-linear: as the black hole mass doubles, the rate of accretion increases four-fold, so that the growth timescale decreases rapidly as the black hole grows. If the density of the external medium remains constant, the black hole will grow to an infinite mass in a finite time given by:
where depends on the effective equation of state of the interstellar medium (ISM), the relative motion of the surrounding gas, and the initial black hole seed mass. This behaviour is illustrated in the inset to Fig. 2. In practice, of course, the black hole does not grow indefinitely, since it will eventually exhaust or expel the surrounding gas. This sets an upper limit to the rapid growth of the black hole, and we assume that further black hole growth is limited by the binding energy within the cooling radius of the dark matter halo. We estimate this limiting mass in §2.3. Because the coefficient depends on the initial seed mass, we might expect the growth of black holes in the EAGLE simulations to be very sensitive to this parameter. We find, however, that a reduction in the initial black hole mass is compensated by an increase in the density around the black hole (§3.4). As a result, the halo mass at which black holes enter the rapid growth phase is unchanged. As we explain below, this mass-scale is determined by the ability of feedback from star formation to regulate the gas density within the galaxy, particularly in its centre, and is weakly dependent on the details of the black hole accretion physics.
2.2 The buoyancy of star formation driven outflows
The critical factor is therefore the ability of the star formation-driven outflow to effectively carry mass away from the galaxy, regulating its density and limiting black hole growth. Recent numerical simulations (Scannapieco
& Brüggen, 2015) have shown that supernova cannot efficiently launch a ballistic wind and that most mass escapes the galaxy in a hot diffuse form.
In the absence of a corona, gas hotter than the virial temperature will flow out as a rarefraction wave by converting its thermal energy into
a bulk flow. This type of escape is seen in simulations of isolated disk patches (eg., Creasey
et al., 2013). The outflow leaves behind a diffuse
corona which slowly cools. As the halo mass increases, and the characteristic cooling time of the halo becomes
longer, the density of the corona increases (White &
Frenk, 1991; Birnboim &
Dekel, 2003). Once the mass of the corona becomes comparable to that of the outflow, supernova heated gas must be buoyant compared to the surrounding gas corona in order to escape. (Although a sufficiently energetic outflow may sweep up and expel all of the halo material, this places more stringent requirements on the outflow.)
If the corona is hot and dense (ie., formed of high entropy material), the outflow will not be buoyant and will stall, cool and fall back onto the galaxy. The buoyancy of the outflow is determined by the ratio of its ‘adiabat’
To illustrate this, we begin by comparing the adiabat of the outflow to that of the halo. To determine the adiabat of the outflow, we assume that roughly half the energy available from supernovae is used to drive an outflow with a mass loading .
Specifically, we assigning the outflow a specific energy of (60% of the supernova energy available from
a Chabrier initial stellar mass function, Crain
et al. (eg., 2015), and equivalent to a temperature of ).
To estimate the density of the entrained/heated gas, we assume that galaxy disk sizes are set by the angular momentum of
the dark matter halo (Mo
et al., 1998) so that and that disk mass is proportional to the halo mass. (We use the notation : at high redshift, where the effect of dark energy
In order for the outflow to be a significant loss of mass from the galaxy disk we require that , and define a critical adiabat by setting so that .
In contrast to the outflow, the adiabat of the galaxy’s gas corona increases strongly with the mass of the halo:
assuming the inner parts of the corona have a temperature of
in hydrostatic equilibrium, and a density of 30 times the times the baryonic virial density
Equating the critical outflow abiabat () with that of the halo (Eq. 4) picks out a critical halo mass scale of
In systems less massive than , the outflow is buoyant, or escapes as a rarefraction wave, carrying away a significant fraction of the star forming gas, but in systems more massive than , the outflow is either trapped by the surrounding corona (if ) or does not cause a significant loss of mass (if the outflow adjusts so that this requires ). The gas that is not consumed by star formation then builds up in the central regions until rapid black hole growth is triggered. Our analysis does not consider the angular momentum carried away by the outflow: we assume that tidal torques mediated by spiral density waves and clumping instabilities will carry gas into the centre of the system to feed the black hole.
The halo mass dependence of the mass loading factor
The argument we have presented above explains the existence of a critical halo mass scale, above which the star formation is unable to drive an outflow with significant mass loading (see also Keller et al., 2016). In haloes less massive than , the outflow will be buoyant if , but we can take the argument a step further by assuming that adjusts so that the outflow is only just buoyant (ie., ). This is plausible. When a parcel of gas starts to be heated it remains trapped until its adiabat is comparable to that of the halo. Once the adiabat exceeds , however, the gas element will start to rise, so that it leaves the star forming region and is not heated further. As it leaves, it is replaced by another element, so that the energy output is spread over a greater gas mass. This is an overly simple interpretation of the complex physics of the interstellar medium, but it suggests that lower mass galaxies will drive outflows with a greater mass loading.
Using this argument, we can estimate the dependence of the mass loading on . Equating this to the halo adiabat, we can estimate the halo mass dependence of the outflow:
Note that the redshift dependence of the mass loading in this equation is encapsulated in the redshift dependence of .
A simple equilibrium model
The mass scaling we have derived from the buoyancy argument above can be combined with a simple equilibrium model to investigate the evolution of the host galaxy. The model we present is deliberately as simplified as possible and is presented solely to motivate the choice of fitting functions used in the following section. A much more complete overview of equilibrium models can be found in Dekel & Mandelker (2014).
The basis of equilibrium models is that the inflow of gas from the halo (roughly) balances star formation and outflow, ie., Since we are concerned with low mass haloes, we assume that the mass inflow rate of gas is proportional to the halo growth rate, and that most of this mass flows directly onto the galaxy in cold streams. We also assume that the outflow is the dominant sink for this material
where is the baryon fraction, and is the recycled stellar fraction. We write the outflow mass loading as a function of halo mass, , and note that power-law solutions of Eq. 7 must be of the form:
Due to the logarithmic nature of the problem, these solutions are attractors, and if star formation and gas accretion are sufficiently rapid, memory of the initial stellar mass is quickly forgotten. Our analysis of the buoyancy of the outflow suggests , and that the – relation should have a low mass slope close to , as is indeed seen in observational data (eg., Moster et al., 2013; Behroozi et al., 2013). This lends credence to our anzatz that the outflow will be only just buoyant.
Relation to the central gas density
In order to complete the analysis, we need to relate the outflow efficiency to the build-up of gas in the central regions of the galaxy. We cannot expect to relate the global mass-loading factor directly to the density of gas on sub-kiloparsec scales around the black hole. Indeed, in §3.4, we show that a decrease in the normalisation factor in the black hole accretion timescale is compensated by an increase in . A full description of this self-regulation would be complex (see §3.4 for a discussion of the numerical simulations).
We can, however, motivate a choice of scaling by assuming that the fraction of gas that is able to reach the central regions of the galaxy is inversely proportional to . Using the gas density from Eq. 3, and rounding to , this suggests
where is the gas density around the blackhole, expressed in Hydrogen atoms per cm, and is a normalisation constant.
We need to turn to the reference EAGLE simulation to see if the scaling suggested above is reasonable in practice, and to determine the normalisation constant. In §3.3, Fig. 5 shows that the density surrounding the black hole is indeed reasonably well described by Eq. 9. In the simulation, the slope of the mass dependence appears to increase at with redshift. This can be qualitatively understood if gas collects rapidly at in the central regions once the outflow becomes in efficient. The effect of adopting a redshift dependent slope is to trigger rapid growth in slightly lower mass haloes at high redshift, making BH trajectories (in the – plane) more similar.
Fig. 5 shows that a suitable choice of normalisation is . With this choice of normalisation, black holes in present-day haloes have growth time, .
2.3 Black hole growth in high mass haloes
The focus of this paper is to explain the rapid growth of black holes at the transition mass scale, , and their low accretion rates below it. Above this mass scale, we assume that the mass of the black hole is limited to a fraction of the halo binding energy. Using the total halo binding energy, however, does not take into account the long cooling times of massive haloes in the present-day Universe. We therefore compute the binding energy within a halo mass dependent radius. For the relevant halo masses, the cooling rates of metal enriched plasmas are not particularly sensitive to temperature, so that the cooling time will roughly scale as
where and are the halo virial radius and temperature respectively. Equating the cooling time to the age of the Universe at redshift suggests that the radius within which we should calculate the binding energy scales as roughly relative to the virial radius. We then calculate the binding energy within this radius (allowing for the variation in concentration with halo mass (Duffy et al., 2008)) to compute the upper limit to the black hole mass. Clearly this estimate is approximate, and our approach is deliberately simple. A more advanced analytic treatment would account for low entropy gas expelled from the system (McCarthy et al., 2011) and the impact of the subsequent rearrangement on the halo cooling time. The treatment would also need to account for the growth of black holes due to black hole – black hole mergers. All of these effects are, however, accounted for in the EAGLE numerical simulations.
2.4 Linking black hole growth and halo growth
Today’s supermassive black holes are thought to originate from black hole ‘seeds’ of mass to a few times that are created from the direct collapse of metal free gas clouds (eg., Begelman et al., 2006; Regan & Haehnelt, 2009; Mayer et al., 2010). We translate the growth of black holes into a relation between black hole mass and halo mass by assuming that black holes are created as seeds in haloes capable of hosting a small galaxy, and then comparing the black hole growth rate with that of its host dark matter halo. In a -CDM cosmology, the halo growth rate is well approximated by Correa et al. (2015)
The equilibrium model allows us to link the halo growth rate, the effectiveness of star formation driven outflows and the gas density in the galaxy. The density of gas surrounding the black hole depends critically on the effectiveness of this feedback from star formation. If star formation-driven outflows are not effective at limiting the density of the ISM, the black hole will be driven into its rapid growth phase. A 100-fold increase in gas density translates into an order of magnitude decrease in the timescale for black hole growth, and the timescale would then fall from being longer than the age of the Universe to being a small fraction of it. Black hole growth and the effectiveness of star formation-driven outflows are thus intertwined.
Combining Eqs. 11 and 9 enables us to plot the trajectory of a black hole in the black hole mass – halo mass plane. To match the EAGLE simulations, we assume that black holes have an initial seed mass of . Seeding black holes at different times creates a variety of trajectories, which we illustrate by coloured lines in Fig. 2. The transition halo mass at which this rapid growth phase occurs is rather weakly dependent on the seeding time, and the spread in trajectories is further compressed if the trajectory is plotted as a function of stellar mass. By connecting the final black hole masses, we predict the black hole mass – halo mass relation at the present day (shown as a black dashed line). The model is idealised, and we would expect scatter in both the growth rates of haloes and the density of gas surrounding the black hole. This will lead to the populations of red and blue galaxies overlapping at the transition mass scale: red systems having significantly more massive black holes than blue galaxies of equal halo mass (as seen in Fig. 1).
The model we have presented makes some other remarkable observable predictions. Firstly, it suggests that the transition mass scale should be weakly dependent on redshift. This is indeed confirmed by observations that measure the fraction of red sequence galaxies as a function of mass (as can be seen by comparing Fig. 1 with the equivalent plot for present-day galaxies in Kauffmann et al. (2003)). Secondly, the model predicts that the relation between black hole mass and galaxy mass will be strongly non-linear, rising steeply below a galaxy mass of , before flattening to an almost linear relation.
We compare the prediction of the analytic model to observational data in Fig. 4, using the observationally derived relation between halo mass and galaxy mass (Moster et al., 2013) to convert the black hole mass – halo mass relation predicted by the analytic model into a relation between black hole mass and galaxy mass. The prediction of the analytic model is shown as a thick dashed line. The observational data are shown as red and blue points, depending on the host galaxy morphology (red and blue corresponding to early and late-type morphologies, respectively). A key point in this plot is that recent observational developments have made it possible to measure the masses of black holes over a wide range of galaxy mass and morphology (eg., McConnell & Ma, 2013; Savorgnan et al., 2016). To plot total galaxy mass, rather than than bulge mass alone, we convert the mid-IR photometry given by Savorgnan et al. (2016) to stellar mass using the calibration of Meidt et al. (2014). The error bars on the point show the uncertainty in this conversion. The lower limit shows the bulge mass, while the upper limit assumes there is no non-stellar contribution to the emission. The data points correspond to the assumption that 50% of the disk emission is from diffuse gas in late-type systems, and 10% in early-type systems. When plotted in this way, the non-linear relation between black hole mass and galaxy mass becomes clear, as does the connection between black hole mass and galaxy morphology. The observations confirm the model’s prediction that early-type (red sequence) galaxies harbour more massive black holes than late-type (blue sequence) galaxies of equivalent mass.
Clearly the analytic model makes the prediction that isolated galaxies with masses below the current limit of the observation data will have black hole masses less than . Unfortunately, it is extremely hard to directly measure the masses of such small black holes. Further more, it is critical to correctly account for non-detections and to avoid contamination of the sample by tidally stripped satellite galaxies (Barber et al., 2016). By estimating black hole masses from the emission line widths, Reines & Volonteri (2015) give us a glimpse of the relation at lower masses. The results are fully consistent with the predictions of the analytic model.
3 Validating the Analytic Model with the EAGLE Hydrodynamic Simulations
3.1 The Eagle simulation suite
In order to validate our simple analytic model, we have analysed the growth of black holes in the state-of-the-art EAGLE cosmological simulations (Schaye et al., 2015; Crain et al., 2015). The simulation suite consists of a large number of hydrodynamical simulations of galaxy formation, including different resolutions, simulated volumes and physical models. The simulations use advanced smoothed particle hydrodynamics (ANARCHY, Schaller et al. (2015)), based on the Gadget code (Springel, 2005), and phenomenological subgrid models to capture unresolved physics including cooling, metal enrichment and energy input from stellar feedback and black hole growth (Dalla Vecchia & Schaye, 2008; Wiersma et al., 2009b; Wiersma et al., 2009a). A complete description of the code can be found in Schaye et al. (2015). The simulations successfully reproduce key observational properties of galaxies, including the evolution of the galaxy stellar mass function, galaxy sizes and the transition mass scale for galaxy star formation rates (Schaye et al., 2015; Furlong et al., 2015b, a; Trayford et al., 2016). The simulations also provide a reasonable description of the observed black hole X-ray luminosity function (Rosas-Guevara et al., 2016). The simulations uses -CDM cosmological parameters taken from the Planck Collaboration 2014 measurements (Planck Collaboration et al., 2014), specifically , , , and .
The largest EAGLE simulation, which we label REFERENCE and use in most plots, uses a simulation carried out in a comoving volume of using baryon particles of initial mass and the same number of dark matter particles of mass . The simulation uses a gravitational softening that is no larger than (we use the notation to clearly distinguish proper length scales from comoving length scales). A convergence study is shown in Appendix A: we find that the results we present are not changed by a factor of 8 decrease in the particle mass. We use additional simulation models to explore the impact of different subgrid physics assumptions. These are listed in Table 1. The simulation outputs were analysed using the SUBFIND programme to identify bound sub-structures (Springel et al., 2001; Dolag et al., 2009) within each dark matter halo. We identify these substructures as galaxies and measure stellar masses within a radius of . The black hole mass we quote corresponds to the central black hole of each galaxy.
3.2 Black holes and galaxies in the EAGLE simulations
In the EAGLE simulations, the accretion on to a black hole is determined by a sub-grid model that accounts for the mean density, effective sound speed and relative motion and angular momentum of the surrounding gas as detailed in Rosas-Guevara et al. (2015) (with the exception that, in the EAGLE simulations, we do not increase the accretion rate to account for an unresolved clumping factor). The model assumes that once gas reaches high densities on sub-kpc scales around the black hole, it will be accreted on to the black hole at the Bondi rate unless its angular momentum is sufficiently high to prevent this. Our simulations allow us to understand how the density on sub-kpc scales around the black hole is determined by the interaction of star formation, feedback and gas accretion (on scales of 1 kpc and greater) without making the simplifications adopted in the analytic model. Importantly, the simulation does not impose a galaxy transition mass scale by varying the black hole feedback efficiency as a function of halo mass or accretion rate. We always assume a that the energy generated by feedback is 1.5% of the rest mass energy of the accreted material.
The reference EAGLE simulation accurately predicts the bimodality of galaxy types and the existence of a tight and rapidly evolving sequence of star formation growth timescales in blue-sequence galaxies (Furlong et al., 2015b; Trayford et al., 2016). This is illustrated in Fig. 1, where coloured points show a selection of galaxies from the simulation. The simulation also reproduces many aspects of the observed AGN population well, including the creation of a population of X-ray bright AGN in the first billion years of the Universe (Rosas-Guevara et al., 2016).
In Fig. 4, we show galaxies sampled from the simulation as contour lines and a background image coloured by the average star formation growth timescale, . The simulation follows both the analytic model and the observations well, reproducing the red-/blue-sequence dichotomy. Furthermore, the colouring of the image shows that the transition is driven by black hole growth: at a given mass, galaxies lying on the blue (star forming) sequence have black holes an order of magnitude smaller than equivalent red-sequence (passive) galaxies, showing the importance of black hole growth in this transition region. The agreement on the transition mass scale between the analytic model and the simulation is encouraging since the choice of the parameter was based on the gas density around black holes (in haloes) in the simulation, and not fit to reproduce the transition mass seen in the simulation. The greatest value of the hydrodynamic simulations is, however, that we can vary the parameters that control feedback from star formation and the accretion rates of black holes, and thus experiment with their impact on the resulting galaxy and black hole correlations. In the following section, we will explore how the value of changes in response to variations in the simulation parameters, gaining deeper insight into the connection between the galaxies corona and the density of gas around the black hole.
3.3 Density of gas around the black hole
The ability of the galaxy to launch an outflow is part of a complex story that cannot be completely captured in a simple analytic model, and we turn to the full hydrodynamic simulation to support the functional form adopted for the density of gas around the black hole on scales (the resolution of the reference simulation). Fig. 5 shows that the density of material surrounding the black hole is a strong function of both redshift and halo mass. The high densities reached at high redshift bring the average black hole accretion rates close to the Eddington limit, and allows the model to explain the presence of high redshift quasars (Mortlock et al., 2011; Rosas-Guevara et al., 2016).
The functional form suggested by the analytic model (dashed lines) matches that in the simulation reasonably well, capturing the main features, particularly at lower redshift. At higher redshift, the simple analytic model underestimates the build-up of the gas density in haloes more massive than seen in the simulation. This has the effect of reducing the spread in the halo mass at which black holes enter their rapid growth phase, leading to a situation in which the most rapidly growing black holes are found in haloes of mass , regardless of redshift. Black holes in lower-mass haloes have yet to enter the non-linear growth phase, while black holes in higher-mass haloes have injected so much energy into the host halo, re-arranging its density and entropy structure (eg., McCarthy et al., 2011), that inflowing gas is unable to remain cool and provide further fuel for star formation (eg., Dubois et al., 2013).
3.4 Parameter variations
The EAGLE simulation suite includes many variations of model parameters. In the model we have presented, the accretion rate onto the black hole undergoes a transition that is driven by the hot halo of the galaxy. This predicts that the transition mass scale should be insensitive to the star formation law and the details of the AGN accretion and feedback model. In Fig. 6, we show three models in which the AGN accretion parameter and black hole seed mass are varied. These models are discussed in more detail in Crain et al. (2015) and Salcido et al. (2016), and were run in smaller, , volumes compared to the REFERENCE simulation (but with the same resolution). In models EagleVariation_ALPHA1p0e8 and EagleVariation_ALPHA1p0e4, the parameter controlling the angular momentum dependence of the black hole accretion rate (and thus the coefficient in Eq. 2) is increased or decreased by a factor 100, respectively. In the left hand panel, we show the relation between halo mass and the central black hole mass. Despite varying the normalisation of the black hole accretion rate parameter by 4 orders of magnitude, the models show remarkable similarity to the reference EAGLE simulation. The initial black hole masses in the EagleVariation_SEED1p0e4 simulation are a factor 10 smaller than in the reference run. Despite this, massive black holes are able to grow, following a similar relation to that seen in the REFERENCE simulation. The robustness of the model is remarkable, confirming our analytic model in which the growth of the black hole, the steep transition in the black hole mass-halo mass relation, and the transition in galaxy star formation rates are driven by the properties of the halo, and are not a coincidence of the details of the black hole accretion parameters.
The right-hand panel of Fig. 6 provides an explanation of this behaviour, showing the density of the gas surrounding the black hole as a function of halo mass in the same models. Models in which the black hole accretion rate normalisation is lower, or the seed black hole mass is smaller, reach higher densities, compensating for the lower normalisation of the initial accretion rate. This is the behaviour expected in our model: when the star formation-driven outflow ceases to be effective in higher-mass haloes, the gas density in the centre of the galaxy builds up until efficient black hole accretion and feedback occur.
In Fig. 7, we show the central black hole mass and black hole gas density relations for three further model variations. The simulations EagleVariation_EOS1.666 and EagleVariation_AGNdT9_ALPHA1p0e4 show the effects of varying the exponent of the equation of state assumed for the ISM and the heating temperature used in the AGN feedback model. These simulations were run in and boxes, respectively, using the same resolution as the REFERENCE simulation. The increase in the equation of state exponent, from to , reduces the compressibility of the ISM. Despite its impact on the internal properties of galaxies, the assumed equation of state has relatively little impact on the transition mass scale since the gas reaches similar densities around the black hole. The higher AGN heating temperature in simulation EagleVariation_AGNdT9_ALPHA1p0e4 was introduced in order to increase the abiabat of gas heated by AGN feedback (but not its overall energy), but this also has little impact on the masses of black holes or the density of the gas that builds up around them. This model does, however, result in more realistic X-ray properties in galaxy groups (Schaye et al., 2015).
Fig. 7 shows a third model, in which stellar feedback is switched off. This model is also shown as a dot-dashed line in Fig. 4. The resulting black hole masses and densities are markedly different in this case. Since stellar feedback is no longer removing gas from small galaxies, the density around the black hole builds up regardless of the presence or absence of a hot, high-entropy corona. As a result, no transition in the masses or activity of the black hole is seen. The model does not match the observed stellar (or black hole) properties of galaxies, clearly demonstrating the need for a transition between stellar and AGN feedback in setting the properties of galaxies and establishing the existence of the characteristic mass scale in galaxy formation.
4 Discussion and Conclusions
Galaxy properties show an abrupt change in behaviour around a halo mass of . In lower-mass haloes, central galaxies are rapidly star forming, doubling their stellar mass through star formation on a timescale that is shorter than the age of the Universe. In higher-mass haloes, galaxies have much longer star formation growth timescales (Fig. 1) and grow primarily through galaxy mergers and the accretion of external stars. This dichotomy is often referred to as the blue and red galaxy sequences, respectively. The change in star formation rate at the transition mass scale leads to a change in slope of the galaxy mass – halo mass relation, creating a break in the galaxy mass function at a galaxy mass scale of . The decline in the importance of on-going star formation and the rise in the importance of growth by mergers may also account for the shift from late to early-type galaxy morphology at this mass scale. Explaining the origin of this mass scale is of fundamental importance to explaining the structure of the observable Universe.
In this paper, we have developed a simple analytic model in which black holes and star formation compete to regulate the gas content of a galaxy as it grows by accretion from the cosmic web. The main elements of the model are the non-linear growth of black holes in the Bondi accretion regime and the buoyancy of star formation - driven outflow relative to galaxies hot coronae.
In the Bondi regime, black holes grow slowly before abruptly switching to a rapid growth phase (Fig. 2). The onset of this phase is highly dependent on the surrounding gas density, and thus on the ability of star formation driven outflows to prevent gas accumulating in the centre of the galaxy. In our model, the effectiveness of stellar feedback depends on the buoyancy of the outflow relative to the corona (in the absence of a corona, the outflow escapes as a rarefraction wave). We show that this leads to a critical halo mass above which star formation driven outflows are unable to prevent the build up of gas. Making simple assumptions about the evolution of the gas density of star forming galaxies, we compare the buoyancy of the star formation driven outflow to that of the halo (Fig. 3) in order to determine the dependence of the central gas density on halo mass and redshift. The model links the build up of the gas density in the central regions of the galaxy, the on-set of rapid black hole growth, the galaxy’s eventual transition to the red sequence and the build-up of the hot corona.
In order to test the simplified analytic model, we compare the model to the growth of black holes in the EAGLE hydrodynamic simulations. The EAGLE simulations assume a constant efficiency of black hole feedback, but at the same time obtain a remarkable match to the observed properties of galaxies, including the abrupt change in galaxy star formation rates in the transition regime (Fig. 1). Although no halo mass dependence of black hole accretion rates is imposed in the simulations, the transition in galaxy properties emerges as we would expect from the analytic model (Fig. 4).
The greatest value of the hydrodynamic simulations is, however, that we can vary the parameters that control feedback from star formation and the accretion rates of black holes, and thus experiment with their impact on the resulting galaxy and black hole correlations. This enables us to confirm the causal connections implied by the simple analytic model. For a fixed density, varying the implicit accretion disk viscosity, or the seed black hole mass, alters , and hence the onset of the rapid growth phase, . In practice, however, reductions in are compensated by an increasingly steep relation between the gas density around the black hole and halo mass, and the galaxy transition mass is remarkably insensitive to the parameter choice (Fig. 6). Variations in the effective equation of state assumed for star-forming gas or the heating temperature of AGN feedback (Fig. 7) also have little impact on the transition mass scale. If we reduce or eliminate feedback from supernovae, however, we find that black holes are able to grow effectively in haloes of all masses, as illustrated by the dot-dashed line in Fig. 4 and the cyan line in Fig. 7. This change eliminates the transition mass scale and consequently does not reproduce the properties of observed galaxies. This confirms that, below the transition mass, black hole growth is suppressed by stellar feedback.
Galaxies fall into two distinct galaxy types, characterised by the presence or absence of significant on-going star formation. This dichotomy is established in the early universe and driven by a transition in galaxy properties at a mass of . We have presented a new way of understanding this transition and, thus, the origin of the distinct red and blue galaxy sequences. Rapid growth of black holes is triggered at a halo mass of by the development of a sufficiently hot diffuse gas corona which confines the star formation-driven outflow by preventing it from rising buoyantly. This simple analytic model is supported by observational data and confirmed in numerical experiments. It makes predictions for the relations between the growth of galaxies and their black holes which will be tested by forthcoming observations. In particular, the model predicts that the black holes present in isolated low-mass galaxies will be small, with most objects having black holes smaller than .
We gratefully acknowledge Oliver Ilbert and COSMOS team for assistance with the star formation rate data for young galaxies, and Alistair Graham and Sandra Savorgnan for assistance with the stellar and black holes masses of local galaxies. This work was supported by STFC grant ST/L00075X/1 and used the DiRAC Data Centric system at Durham University, operated by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant ST/K00042X/1, STFC capital grant ST/H008519/1, and STFC DiRAC Operations grant ST/K003267/1 and Durham University. DiRAC is part of the National E-Infrastructure. We also gratefully acknowledge PRACE for awarding us access to the resource Curie based in France at Trs Grand Centre de Calcul. This work was sponsored by the Dutch National Computing Facilities Foundation (NCF) for the use of supercomputer facilities, with financial support from the Netherlands Organization for Scientific Research (NWO). The research was supported in part by the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013) / ERC Grant agreements 278594-GasAroundGalaxies, GA 267291 Cosmiway, and 321334 dustygal, the Interuniversity Attraction Poles Programme initiated by the Belgian Science Policy OWNce ([AP P7/08 CHARM]), the National Science Foundation under Grant No. NSF PHY11-25915, the UK Science and Technology Facilities Council (grant numbers ST/F001166/1 and ST/I000976/1), Rolling and Consolidated Grants to the ICC, Marie Curie Reintegration Grant PERG06-GA-2009-256573, RAC is a Royal Society University Research Fellow. RGB, JS and RAC thank the Simons Foundation for their hospitality at the ‘Galaxy Superwinds’ symposium where this work was completed. The data used in this project is available from the EAGLE public data base at http://icc.dur.ac.uk/Eagle/database.php (McAlpine et al., 2016).
Appendix A Numerical Convergence
For the reference EAGLE simulation, we are able to investigate the effect of repeating simulations at higher numerical resolution. We compare REFERENCE to two simulations in which the particle mass is 8 times smaller (and the gravitational softening a factor of 2 smaller) in Fig. 8, reproducing the black hole mass and black hole density panels discussed for the model variations in §3.4. The higher resolution simulations were carried out in a volume. Using the nomenculture introduced in Schaye et al. (2015), the simulations test the ‘strong’ convergence (in which we keep all parameters fixed at their fiducial values) in HIRES-REFERENCE and the ‘weak’ convergence of the simulation in HIRES-RECALIBRATED. As discussed by Schaye et al., a calculation with higher resolution resolves additional physical processes and that may require recalibration of the parameters encapsulating the subgrid processes. In the left hand panel, which shows the black hole mass – halo mass relation, both approaches to convergence yield similar results, and agree extremely well with the REFERENCE calculation on which our parameter investigations are based. In the right hand panel, which shows the gas density around the black hole as a function of halo mass, the two higher resolution simulations differ. This is expected since the parameter, controlling the angular momentum dependence of the accretion rate, is a factor of 100 smaller in HIRES-RECALIBRATED than in the reference run. This confirms that the parameter dependencies investigated at the resolution of the REFERENCE simulation also hold in the higher-resolution simulations.
- pubyear: 2016
- pagerange: The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end–LABEL:lastpage
- We use for the Boltzmann constant. The mean molecular weight of the plasma is where is the mass of hydrogen atom and for a fully ionized plasma of primordial composition. Thermodynamic entropy is proportional the logarithm of the adiabat.
- and are the cosmological density parameters for matter and dark energy respectively.
- We assume a density of 30 times where is the baryon density parameter and is the critical density. The temperature we assume is twice the virial temperature of the halo since the radius of interest is close to the peak of the halo’s rotation curve. We define the virial temperature as , where is the radius enclosing an overdensity of 200 times and is the mass within it.
- Baldry I. K., Balogh M. L., Bower R. G., Glazebrook K., Nichol R. C., Bamford S. P., Budavari T., 2006, \mnras, 373, 469
- Barber C., Schaye J., Bower R. G., Crain R. A., Schaller M., Theuns T., 2016, \mnras, 460, 1147
- Begelman M. C., Volonteri M., Rees M. J., 2006, \mnras, 370, 289
- Behroozi P. S., Marchesini D., Wechsler R. H., Muzzin A., Papovich C., Stefanon M., 2013, \apjl, 777, L10
- Benson A. J., Bower R. G., Frenk C. S., Lacey C. G., Baugh C. M., Cole S., 2003, \apj, 599, 38
- Binney J., Tabor G., 1995, \mnras, 276, 663
- Birnboim Y., Dekel A., 2003, \mnras, 345, 349
- Bondi H., 1952, \mnras, 112, 195
- Bondi H., Hoyle F., 1944, \mnras, 104, 273
- Booth C. M., Schaye J., 2009, \mnras, 398, 53
- Bouché N., et al., 2010, \apj, 718, 1001
- Bower R. G., Lucey J. R., Ellis R. S., 1992, \mnras, 254, 601
- Bower R. G., Benson A. J., Malbon R., Helly J. C., Frenk C. S., Baugh C. M., Cole S., Lacey C. G., 2006, \mnras, 370, 645
- Brinchmann J., Charlot S., White S. D. M., Tremonti C., Kauffmann G., Heckman T., Brinkmann J., 2004, \mnras, 351, 1151
- Chang Y.-Y., van der Wel A., da Cunha E., Rix H.-W., 2015, \apjs, 219, 8
- Correa C. A., Wyithe J. S. B., Schaye J., Duffy A. R., 2015, \mnras, 450, 1521
- Crain R. A., et al., 2015, ArXiv e-prints, arXiv:1501.01311,
- Creasey P., Theuns T., Bower R. G., 2013, \mnras, 429, 1922
- Croton D. J., et al., 2006, \mnras, 365, 11
- Dalla Vecchia C., Schaye J., 2008, \mnras, 387, 1431
- Darvish B., Mobasher B., Sobral D., Rettura A., Scoville N., Faisst A., Capak P., 2016, preprint, (arXiv:1605.03182)
- Davé R., Finlator K., Oppenheimer B. D., 2012, \mnras, 421, 98
- De Lucia G., Springel V., White S. D. M., Croton D., Kauffmann G., 2006, \mnras, 366, 499
- Dekel A., Mandelker N., 2014, \mnras, 444, 2071
- Di Matteo T., Springel V., Hernquist L., 2005, \nat, 433, 604
- Dolag K., Borgani S., Murante G., Springel V., 2009, \mnras, 399, 497
- Dubois Y., Pichon C., Devriendt J., Silk J., Haehnelt M., Kimm T., Slyz A., 2013, \mnras, 428, 2885
- Dubois Y., Volonteri M., Silk J., Devriendt J., Slyz A., Teyssier R., 2015, \mnras, 452, 1502
- Dubois Y., Peirani S., Pichon C., Devriendt J., Gavazzi R., Welker C., Volonteri M., 2016, preprint, (arXiv:1606.03086)
- Duffy A. R., Schaye J., Kay S. T., Dalla Vecchia C., 2008, \mnras, 390, L64
- Fabian A. C., Sanders J. S., Allen S. W., Crawford C. S., Iwasawa K., Johnstone R. M., Schmidt R. W., Taylor G. B., 2003, \mnras, 344, L43
- Finlator K., Davé R., 2008, \mnras, 385, 2181
- Furlong M., et al., 2015a, preprint, (arXiv:1510.05645)
- Furlong M., et al., 2015b, \mnras, 450, 4486
- Gabor J. M., Davé R., 2015, \mnras, 447, 374
- Habouzit M., Volonteri M., Dubois Y., 2016, preprint, (arXiv:1605.09394)
- Harrison C. M., et al., 2012, \mnras, 426, 1073
- Ilbert O., et al., 2013, \aap, 556, A55
- Ilbert O., et al., 2015, \aap, 579, A2
- Kauffmann G., et al., 2003, \mnras, 341, 54
- Keller B. W., Wadsley J., Couchman H. M. P., 2016, preprint, (arXiv:1604.08244)
- Khandai N., Di Matteo T., Croft R., Wilkins S., Feng Y., Tucker E., DeGraf C., Liu M.-S., 2015, \mnras, 450, 1349
- Maiolino R., et al., 2012, \mnras, 425, L66
- Mayer L., Kazantzidis S., Escala A., Callegari S., 2010, \nat, 466, 1082
- McAlpine S., et al., 2016, Astronomy and Computing, 15, 72
- McCarthy I. G., Schaye J., Bower R. G., Ponman T. J., Booth C. M., Dalla Vecchia C., Springel V., 2011, \mnras, 412, 1965
- McConnell N. J., Ma C.-P., 2013, \apj, 764, 184
- Meidt S. E., et al., 2014, \apj, 788, 144
- Meier D. L., 1999, \apj, 522, 753
- Mo H. J., Mao S., White S. D. M., 1998, \mnras, 295, 319
- Mortlock D. J., et al., 2011, \nat, 474, 616
- Moster B. P., Naab T., White S. D. M., 2013, \mnras, 428, 3121
- Nemmen R. S., Bower R. G., Babul A., Storchi-Bergmann T., 2007, \mnras, 377, 1652
- Parry O. H., Eke V. R., Frenk C. S., 2009, \mnras, 396, 1972
- Peng Y.-j., et al., 2010, \apj, 721, 193
- Planck Collaboration et al., 2014, \aap, 571, A16
- Qu Y., Di Matteo P., Lehnert M. D., van Driel W., Jog C. J., 2016, \mnras, 000, 000
- Regan J. A., Haehnelt M. G., 2009, \mnras, 393, 858
- Reines A. E., Volonteri M., 2015, \apj, 813, 82
- Rodriguez-Gomez V., et al., 2016, \mnras, 458, 2371
- Rosas-Guevara Y. M., et al., 2015, \mnras, 454, 1038
- Rosas-Guevara Y., Bower R. G., Schaye J., McAlpine S., Dalla-Vecchia C., Frenk C. S., Schaller M., Theuns T., 2016, preprint, (arXiv:1604.00020)
- Salcido J., Bower R. G., Theuns T., McAlpine S., Schaller M., Crain R. A., Schaye J., Regan J., 2016, preprint, (arXiv:1601.06156)
- Salpeter E. E., 1964, \apj, 140, 796
- Savorgnan G. A. D., Graham A. W., Marconi A., Sani E., 2016, \apj, 817, 21
- Scannapieco E., Brüggen M., 2015, \apj, 805, 158
- Schaller M., Dalla Vecchia C., Schaye J., Bower R. G., Theuns T., Crain R. A., Furlong M., McCarthy I. G., 2015, \mnras, 454, 2277
- Schaye J., et al., 2010, \mnras, 402, 1536
- Schaye J., et al., 2015, \mnras, 446, 521
- Sijacki D., Springel V., Di Matteo T., Hernquist L., 2007, \mnras, 380, 877
- Sijacki D., Vogelsberger M., Genel S., Springel V., Torrey P., Snyder G. F., Nelson D., Hernquist L., 2015, \mnras, 452, 575
- Silk J., Rees M. J., 1998, \aap, 331, L1
- Springel V., 2005, \mnras, 364, 1105
- Springel V., White S. D. M., Tormen G., Kauffmann G., 2001, \mnras, 328, 726
- Trayford J. W., Theuns T., Bower R. G., Crain R. A., Lagos C. d. P., Schaller M., Schaye J., 2016, \mnras,
- White S. D. M., Frenk C. S., 1991, \apj, 379, 52
- Wiersma R. P. C., Schaye J., Smith B. D., 2009a, \mnras, 393, 99
- Wiersma R. P. C., Schaye J., Theuns T., Dalla Vecchia C., Tornatore L., 2009b, \mnras, 399, 574