Feast and Famine: regulation of black hole growth in low-redshift galaxies

We analyse the observed distribution of Eddington ratios  ( L / L _Edd)  as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation [   M _*/star formation rate (SFR) ∼  a Hubble time] in their central kiloparsec regions, and is characterized by a broad lognormal distribution of accretion rates peaked at a few per cent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows little dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations (   M _*/SFR ≫  a Hubble time), and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate on to black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power law and the decrease in the accretion rate on to black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.     

The growth and assembly of a massive galaxy at z∼ 2

We study the stellar mass assembly of the Spiderweb galaxy  (MRC 1138−262)  , a massive   z = 2.2  radio galaxy in a protocluster and the probable progenitor of a brightest cluster galaxy. Nearby protocluster galaxies are identified and their properties are determined by fitting stellar population models to their rest-frame ultraviolet to optical spectral energy distributions. We find that within 150 kpc of the radio galaxy the stellar mass is centrally concentrated in the radio galaxy, yet most of the dust-uncorrected, instantaneous star formation occurs in the surrounding low-mass satellite galaxies. We predict that most of the galaxies within 150 kpc of the radio galaxy will merge with the central radio galaxy by   z = 0  , increasing its stellar mass by up to a factor of ≃2. However, it will take several hundred Myr for the first mergers to occur, by which time the large star formation rates are likely to have exhausted the gas reservoirs in the satellite galaxies. The tidal radii of the satellite galaxies are small, suggesting that stars and gas are being stripped and deposited at distances of tens of kpc from the central radio galaxy. These stripped stars may become intracluster stars or form an extended stellar halo around the radio galaxy, such as those observed around cD galaxies in cluster cores.     

Evolution of globular cluster systems in three galaxies of the Fornax cluster

We studied and compared the radial profiles of globular clusters and of the stellar bulge component in three galaxies of the Fornax cluster observed with the WFPC2 of the Hubble Space Telescope ( HST ). The stars are more concentrated toward the galactic centres than globular clusters, in agreement with what has already been observed in many other galaxies: if the observed difference is the result of evolution of the globular cluster systems starting from initial profiles similar to those of the halo–bulge stellar components, a relevant fraction of their initial mass (74, 47 and 52 per cent for NGC 1379, 1399 and 1404, respectively) should have disappeared in the inner regions. This mass has probably contributed to the nuclear field population, local dynamics and high-energy phenomena in the primeval life of the galaxy. An indication in favour of the evolutionary interpretation of the difference between the globular cluster system and stellar bulge radial profiles is given by the positive correlation we found between the value of the mass lost from the globular cluster system and the central galactic black hole mass in the set of seven galaxies for which these data are available.     

Rates of tidal disruption of stars by massive central black holes

There is strong evidence for some kind of massive dark object in the centres of many galaxy bulges. The detection of flares from tidally disrupted stars could confirm that these objects are black holes (BHs). Here we present calculations of the stellar disruption rates in detailed dynamical models of real galaxies, taking into account the refilling of the loss cone of stars on disruptable orbits by two-body relaxation and tidal forces in non-spherical galaxies. The highest disruption rates (one star per 10⁴ yr) occur in faint ( L ≲10¹⁰ L_⊙) galaxies, which have steep central density cusps. More luminous galaxies are less dense and have much longer relaxation times and more massive BHs. Dwarf stars in such galaxies are swallowed whole by the BH and hence do not emit flares; giant stars could produce flares as often as every 10⁵ yr, although the rate depends sensitively on the shape of the stellar distribution function. We discuss the possibility of detecting disruption flares in current supernova searches. The total mass of stars consumed over the lifetime of the galaxy is of the order of 10⁶ M_⊙, independent of galaxy luminosity; thus, disrupted stars may contribute significantly to the present BH mass in galaxies fainter than ∼10⁹ L_⊙.     

Evidence for recent star formation in BCGs: a correspondence between blue cores and UV excess

We present a joint analysis of near-ultraviolet ( NUV ) data from the GALEX ( Galaxy Evolution Explorer ) mission and (optical) colour profiles for a sample of seven brightest cluster galaxies (BCGs) in the Canadian Cluster Comparison Project. We find that every BCG, which has a blue rest-frame UV colour, also shows a blue core in its optical colour profile. Conversely, BCGs that lack blue cores and show monotonic colour gradients typical of old elliptical galaxies are red in the UV. We interpret this as evidence that the NUV enhancement in the blue BCGs is driven by recent star formation and not from old evolved stellar populations such as horizontal branch stars. Furthermore, the UV enhancement cannot be from an active galactic nuclei (AGN) because the spatial extent of the blue cores is significantly larger than the possible contamination region due to a massive black hole. The recent star formation in the blue BCGs typically has an age less than 200 Myr and contributes mass fractions of less than a per cent. Although the sample studied here is small, we demonstrate, for the first time , a one-to-one correspondence between blue cores in elliptical galaxies (in particular BCGs) and a NUV enhancement observed using GALEX . The combination of this one-to-one correspondence and the consistently young age of recent star formation, coupled with additional correlations with the host cluster's X-ray properties, strongly suggests that the star formation is fuelled by gas cooling out of the intracluster medium. In turn, this implies that any AGN heating of the intracluster medium in massive clusters only acts to reduce the magnitude of the cooling flow and once this flow starts, it is nearly always active. Collectively, these results suggest that AGN feedback in present-day BCGs, while important, cannot be as efficient as suggested by the recent theoretical model by proposed by De Lucia et al.     

Black hole demographics from the relation

We analyse a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M _bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M _•, using the tight empirical correlation between M _• and bulge stellar velocity dispersion. The frequency function     is reasonably well described as a Gaussian with     and standard deviation ∼0.45; the implied mean ratio of black hole mass to bulge mass is a factor of ∼5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be ∼     consistent with that inferred from high-redshift     active galactic nuclei.     

Radiative feedback from quasars and the growth of massive black holes in stellar spheroids

We discuss the importance of feedback via photoionization and Compton heating on the co-evolution of massive black holes (MBHs) at the centre of spheroidal galaxies, and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar medium (ISM). We then demonstrate that the observed   M _BH–σ  relation could be established following the conversion of most of the gas of an elliptical progenitor into stars, specifically when the gas-to-stars mass ratio in the central regions has dropped to a low level  ∼0.01  or less, so that gas cooling is no longer able to keep up with the radiative heating by the growing central massive black hole (MBH). A considerable amount of the remaining gas will be expelled and both MBH accretion and star formation will proceed at significantly reduced rates thereafter, in agreement with observations of present-day ellipticals. We find further support for this scenario by evolving over an equivalent Hubble time a simple, physically based toy model that additionally takes into account the mass and energy return for the spheroid evolving stellar population, a physical ingredient often neglected in similar approaches.     

Pre-enriched, not primordial ellipticals

We follow the chemical evolution of a galaxy through star formation and its feedback into the interstellar medium (ISM), starting from primordial gas and allowing for gas to inflow into the region being modelled. We attempt to reproduce observed spectral line strengths for early-type galaxies in order to constrain their star formation histories (SFH). The efficiencies and times of star formation are varied, as are the amount and duration of inflow. We evaluate the chemical enrichment and the mass of stars made with time. Single stellar population (SSP) data are then used to predict line strengths for composite stellar populations. The results are compared with observed line strengths in 10 ellipticals, including some features which help to break the problem of age–metallicity degeneracy in old stellar populations. We find that the elliptical galaxies modelled require high metallicity SSPs (> 3 Z⊙) at later times. In addition, the strong lines observed cannot be produced by an initial starburst in primordial gas, even if a large amount of inflow is allowed for during the first few × 10⁸ yr. This is because some pre-enrichment is required for lines in the bulk of the stars to approach the observed line strengths in ellipticals. These strong lines are better modelled by a system with a delayed burst of star formation, following an early SFH which can be a burst or more steady star formation. Such a model is representative of star formation in normal ellipticals or spirals, respectively, followed by a starburst and gas inflow during a merger or strong interaction with a gas-rich galaxy. Alternatively, a single initial burst of normal stars with a Salpeter initial mass function could produce the observed strong lines if it followed some pre-enrichment process which did not form long-lived stars (e.g. population III stars).     

High-redshift galaxies, their active nuclei and central black holes

We demonstrate that the luminosity function of the recently detected population of actively star-forming galaxies at redshift z  = 3 and the B -band luminosity function of quasi-stellar objects (QSOs) at the same redshift can both be matched with the mass function of dark matter haloes predicted by standard variants of hierarchical cosmogonies for lifetimes of optically bright QSOs anywhere in the range 10⁶ to 10⁸ yr. There is a strong correlation between the lifetime and the required degree of non-linearity in the relation between black hole and halo mass. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M _bh ∝  v ⁵_halo ∝  M ^5/3_halo and a typical duration of the optically bright QSO phase of a few times 10⁷ yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection-dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.     

On the evolution of the black hole: spheroid mass ratio

We present the results of a study which uses the 3C RR sample of radio-loud active galactic nuclei to investigate the evolution of the black hole:spheroid mass ratio in the most massive early-type galaxies from  0 < z < 2  . Radio-loud unification is exploited to obtain virial (linewidth) black hole mass estimates from the 3C RR quasars, and stellar mass estimates from the 3C RR radio galaxies, thereby providing black hole and stellar mass estimates for a single population of early-type galaxies. At low redshift  ( z ≲ 1)  , the 3C RR sample is consistent with a black hole:spheroid mass ratio of   M _bh/ M _sph≃ 0.002  , in good agreement with that observed locally for quiescent galaxies of similar stellar mass  ( M _sph≃ 5 × 10¹¹ M_⊙)  . However, over the redshift interval  0 < z < 2  the 3C RR black hole:spheroid mass ratio is found to evolve as   M _bh/ M _sph∝ (1 + z )^2.07±0.76  , reaching   M _bh/ M _sph≃ 0.008  by redshift   z ≃ 2  . This evolution is found to be inconsistent with the local black hole:spheroid mass ratio remaining constant at a moderately significant level (98 per cent). If confirmed, the detection of evolution in the 3C RR black hole:spheroid mass ratio further strengthens the evidence that, at least for massive early-type galaxies, the growth of the central supermassive black hole may be completed before that of the host spheroid.     

A unified model for the evolution of galaxies and quasars

We incorporate a simple scheme for the growth of supermassive black holes into semi-analytic models that follow the formation and evolution of galaxies in a cold dark matter-dominated Universe. We assume that supermassive black holes are formed and fuelled during major mergers. If two galaxies of comparable mass merge, their central black holes coalesce and a few per cent of the gas in the merger remnant is accreted by the new black hole over a time-scale of a few times 10⁷ yr. With these simple assumptions, our model not only fits many aspects of the observed evolution of galaxies, but also reproduces quantitatively the observed relation between bulge luminosity and black hole mass in nearby galaxies, the strong evolution of the quasar population with redshift, and the relation between the luminosities of nearby quasars and those of their host galaxies. The strong decline in the number density of quasars from z ∼2 to z =0 is a result of the combination of three effects: (i) a decrease in the merging rate; (ii) a decrease in the amount of cold gas available to fuel black holes, and (iii) an increase in the time-scale for gas accretion. The predicted decline in the total content of cold gas in galaxies is consistent with that inferred from observations of damped Ly α systems. Our results strongly suggest that the evolution of supermassive black holes, quasars and starburst galaxies is inextricably linked to the hierarchical build-up of galaxies.     

The evolution of barium and europium in local dwarf spheroidal galaxies

By means of a detailed chemical evolution model, we follow the evolution of barium (Ba) and europium (Eu) in four Local Group Dwarf Spheroidal (dSph) galaxies, in order to set constraints on the nucleosynthesis of these elements and on the evolution of this type of galaxies compared with the Milky Way. The model, which is able to reproduce several observed abundance ratios and the present-day total mass and gas mass content of these galaxies, adopts up-to-date nucleosynthesis and takes into account the role played by supernovae (SNe) of different types (II, Ia) allowing us to follow in detail the evolution of several chemical elements (H, D, He, C, N, O, Mg, Si, S, Ca, Fe, Ba and Eu). By assuming that Ba is a neutron-capture element produced in low-mass asymptotic giant branch stars by s-process but also in massive stars (in the mass range 10–30 M_⊙) by r-process, during the explosive event of SNe of Type II, and that Eu is a pure r-process element synthesized in massive stars also in the range of masses 10–30 M_⊙, we are able to reproduce the observed [Ba/Fe] and [Eu/Fe] as functions of [Fe/H] in all four galaxies studied. We confirm also the important role played by the very low star formation (SF) efficiencies (ν= 0.005–0.5 Gyr⁻¹) and by the intense galactic winds (6–13 times the star formation rate) in the evolution of these galaxies. These low SF efficiencies (compared to the one for the Milky Way disc) adopted for the dSph galaxies are the main reason for the differences between the trends of [Ba/Fe] and [Eu/Fe] predicted and observed in these galaxies and in the metal-poor stars of our Galaxy. Finally, we provide predictions for Sagittarius galaxy for which data of only two stars are available.     

X-ray active galactic nuclei in the core of the Perseus cluster

We present a study of the X-ray emission from the nuclei of galaxies observed in the core of the Perseus cluster in a deep exposure with Chandra . Point sources are found coincident with the nuclei of 13 early-type galaxies, as well as the central galaxy NGC 1275. This corresponds to all galaxies brighter than M _B > −18 in the Chandra field. All of these sources have a steep power-law spectral component and four have an additional thermal component. The unabsorbed power-law luminosities in the 0.5–7.0 keV band range from 8 × 10³⁸ to 5 × 10⁴⁰ erg s⁻¹. We find no simple correlations between the K -band luminosity, or the FUV and NUV AB magnitudes of these galaxies and their X-ray properties. We have estimated the black hole masses of the nuclei using the K -band   M _BH– L _{K bol}  relation and again find no correlation between black hole mass and the X-ray luminosity. Bondi accretion on to the black holes in the galaxies with minihaloes should make them much more luminous than observed.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

The K-band Hubble diagram for the brightest cluster galaxies: a test of hierarchical galaxy formation models

We analyse the K -band Hubble diagram for a sample of brightest cluster galaxies (BCGs) in the redshift range 0< z <1. In good agreement with earlier studies, we confirm that the scatter in the absolute magnitudes of the galaxies is small (0.3 mag). The BCGs exhibit very little luminosity evolution in this redshift range: if q ₀=0.0, we detect no luminosity evolution; for q ₀=0.5, we measure a small negative evolution (i.e., BCGs were about 0.5 mag fainter at z =1 than today). If the mass in stars of these galaxies had remained constant over this period of time, substantial positive luminosity evolution would be expected: BCGs should have been brighter in the past, since their stars were younger. A likely explanation for the observed zero or negative evolution is that the stellar mass of the BCGs has been assembled over time through merging and accretion, as expected in hierarchical models of galaxy formation. The colour evolution of the BCGs is consistent with that of an old stellar population ( z _for>2) that is evolving passively. We can thus use evolutionary population synthesis models to estimate the rate of growth in stellar mass for these systems. We find that the stellar mass in a typical BCG has grown by a factor ≃2 since z ≃1 if q ₀=0.0, or by factor ≃4 if q ₀=0.5. These results are in good agreement with the predictions of semi-analytic models of galaxy formation and evolution set in the context of a hierarchical scenario for structure formation. The models predict a scatter in the luminosities of the BCGs that is somewhat larger than the observed one, but that depends on the criterion used to select the model clusters.     

Can gas dynamics in centres of galaxies reveal orbiting massive black holes?

If supermassive black holes in centres of galaxies form by merging of black hole remnants of massive Population III stars, then there should be a few black holes of mass one or two orders of magnitude smaller than that of the central ones, orbiting around the centre of a typical galaxy. These black holes constitute a weak perturbation in the gravitational potential, which can generate wave phenomena in gas within a disc close to the centre of the galaxy. Here, we show that a single orbiting black hole generates a three-arm spiral pattern in the central gaseous disc. The density excess in the spiral arms in the disc reaches values of 3–12 per cent when the orbiting black hole is about 10 times less massive than the central black hole. Therefore, the observed density pattern in gas can be used as a signature in detecting the most massive orbiting black holes.     

An elementary model for the dust cycle in galaxies

We give an elementary model for the evolution of dust in galaxies, based on abundance arguments. The model takes account of grain core production in both supernovae and giant stars, and includes mantle growth in the interstellar medium. Destruction of grain cores does not appear to be a dominant effect. We show that a self-consistent picture can be made in which the interstellar dust mass is an approximately constant fraction of the heavy element mass in the interstellar medium. This result is demonstrated to be essentially independent of outflow or inflow of interstellar material.     

The relationship between X-ray variability and the central black hole mass

We assembled a sample of Seyfert 1 galaxies, quasi-stellar objects (QSOs) and low-luminosity active galactic nuclei (LLAGNs) observed by ASCA , the central black hole masses of which have been measured. We found that the X-ray variability (which is quantified by the 'excess variance' σ _rms²) is significantly anti-correlated with the central black hole mass, and it is likely that a linear relationship of σ _rms²∝ M _bh⁻¹ exists. It can be interpreted that the short time-scale X-ray variability is caused by some global coherent variations in the X-ray emission region, which is scaled by the size of the central black hole. Hence the central black hole mass is the driving parameter of the previously established relation between X-ray variability and luminosity. Our findings favour the hypothesis that the narrow-line Seyfert 1 galaxies and QSOs harbour smaller black holes than the broad-line objects, and can also easily explain the observational fact that high-redshift QSOs have greater variability than local AGNs at a given luminosity. Further investigations are needed to confirm our findings, and a large sample X-ray variability investigation can give constraints on the physical mechanisms and evolution of AGNs.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

A catalogue of potential adaptive optics survey fields from the UKIRT archive

We present a multicolour catalogue of faint galaxies situated close to bright stars,   V ≲ 15  , with the aim of identifying high-redshift galaxies suitable for study with adaptive optics-equipped near-infrared imagers and spectrographs. The catalogue is constructed from archival calibration observations of the United Kingdom Infrared Telescope (UKIRT) Faint Standard stars with the UKIRT Fast Track Imager (UFTI) camera on UKIRT. We have analysed the deepest 16 fields from the archive to provide a catalogue of galaxies brighter than   K ∼ 20.3  lying between 3 and 25 arcsec of the guide stars. We identify 111 objects in a total survey area of  8.7 arcmin²  . Of these, 87 are classified as galaxies based on their light profiles in our ∼0.5 arcsec median seeing K -band images. 12 of the galaxies have  ( J − K ) ≥ 2.0  consistent with them lying at high redshifts,   z ≳ 2  . These 12 very red galaxies have K -band magnitudes of   K = 18.1–20.1  and separations from the guide stars of 4–20 arcsec and hence are very well suited to adaptive optics studies to investigate their morphologies and spectral properties on sub-kpc scales. We provide coordinates and JHK photometry for all catalogued objects.