Is there an upper limit to black hole masses?

We make a case for the existence for ultra-massive black holes (UMBHs) in the Universe, but argue that there exists a likely upper limit to black hole (BH) masses of the order of   M ∼ 10¹⁰ M_⊙  . We show that there are three strong lines of argument that predicate the existence of UMBHs: (i) expected as a natural extension of the observed BH mass bulge luminosity relation, when extrapolated to the bulge luminosities of bright central galaxies in clusters; (ii) new predictions for the mass function of seed BHs at high redshifts predict that growth via accretion or merger-induced accretion inevitably leads to the existence of rare UMBHs at late times; (iii) the local mass function of BHs computed from the observed X-ray luminosity functions of active galactic nuclei predict the existence of a high-mass tail in the BH mass function at   z = 0  . Consistency between the optical and X-ray census of the local BH mass function requires an upper limit to BH masses. This consistent picture also predicts that the slope of the   M _bh–σ  relation will evolve with redshift at the high-mass end. Models of self-regulation that explain the co-evolution of the stellar component and nuclear BHs naturally provide such an upper limit. The combination of multiwavelength constraints predicts the existence of UMBHs and simultaneously provides an upper limit to their masses. The typical hosts for these local UMBHs are likely the bright, central cluster galaxies in the nearby Universe.     

Iron line profiles including emission from within the innermost stable orbit of a black hole accretion disc

We propose a model for the source of the X-ray background (XRB) in which low-luminosity active nuclei ( L  ∼ 10⁴³ erg s⁻¹) are obscured ( N  ∼ 10²³ cm⁻²) by nuclear starbursts within the inner ∼ 100 pc. The obscuring material covers most of the sky as seen from the central source, rather than being distributed in a toroidal structure, and hardens the averaged X-ray spectrum by photoelectric absorption. The gas is turbulent with velocity dispersion ∼ few × 100 km s⁻¹ and cloud–cloud collisions lead to copious star formation. Although supernovae tend to produce outflows, most of the gas is trapped in the gravity field of the star-forming cluster itself and the central black hole. A hot ( T  ∼ 10⁶ − 10⁷ K) virialized phase of this gas, comprising a few per cent of the total obscuring material, feeds the central engine of ∼ 10⁷ M⊙ through Bondi accretion, at a sub-Eddington rate appropriate for the luminosity of these objects. If starburst-obscured objects give rise to the residual XRB, then only 10 per cent of the accretion in active galaxies occurs close to the Eddington limit in unabsorbed objects.     

Mass and spin co-evolution during the alignment of a black hole in a warped accretion disc

In this paper, we explore the gravitomagnetic interaction of a black hole (BH) with a misaligned accretion disc to study BH spin precession and alignment jointly with BH mass M _BH and spin parameter a evolution, under the assumption that the disc is continually fed, in its outer region, by matter with angular momentum fixed on a given direction     . We develop an iterative scheme based on the adiabatic approximation to study the BH–disc co-evolution: in this approach, the accretion disc transits through a sequence of quasi-steady warped states (Bardeen–Petterson effect) and interacts with the BH until the spin   J _BH  aligns with     . For a BH aligning with a corotating disc, the fractional increase in mass is typically less than a few per cent, while the spin modulus can increase up to a few tens of per cent. The alignment time-scale     is of  ∼10⁵–10⁶ yr  for a maximally rotating BH accreting at the Eddington rate. BH–disc alignment from an initially counter-rotating disc tends to be more efficient compared to the specular corotating case due to the asymmetry seeded in the Kerr metric: counter-rotating matter carries a larger and opposite angular momentum when crossing the innermost stable orbit, so that the spin modulus decreases faster and so the relative inclination angle.     

Ultraluminous X-ray source 1E 0953.8+6918 (M81 X-9): an intermediate-mass black hole candidate and its environs

We present a ROSAT and ASCA study of the Einstein source X-9 and its relation to a shock-heated shell-like optical nebula in a tidal arm of the M81 group of interacting galaxies. Our ASCA observation of the source shows a flat and featureless X-ray spectrum well described by a multicolour disc blackbody model. The source most likely represents an optically thick accretion disc around an intermediate-mass black hole  ( M ∼10² M_⊙)  in its high/soft state, similar to other variable ultraluminous X-ray sources observed in nearby disc galaxies. Using constraints derived from both the innermost stable orbit around a black hole and the Eddington luminosity, we find that the black hole is fast-rotating and that its mass is between ∼80 M_⊙–1.5×10² M_⊙. The inferred bolometric luminosity of the accretion disc is ∼(1.1×10⁴⁰ erg s⁻¹)/(cos  i ). Furthermore, we find that the optical nebula is very energetic and may contain large amounts of hot gas, accounting for a soft X-ray component as indicated by archival ROSAT PSPC data. The nebula is apparently associated with X-9; the latter may be powering the former and/or they could be formed in the same event (e.g. a hypernova). Such a connection, if confirmed, could have strong implications for understanding both the birth of intermediate-mass black holes and the formation of energetic interstellar structures.     

Numerical simulations of hot halo gas in galaxy mergers

Galaxy merger simulations have explored the behaviour of gas within the galactic disc, yet the dynamics of hot gas within the galaxy halo have been neglected. We report on the results of high-resolution hydrodynamic simulations of colliding galaxies with metal-free hot halo gas. To isolate the effect of the halo gas, we simulate only the dark matter halo and the hot halo gas over a range of mass ratios, gas fractions and orbital configurations to constrain the shocks and gas dynamics within the progenitor haloes. We find that (i) a strong shock is produced in the galaxy haloes before the first passage, increasing the temperature of the gas by almost an order of magnitude to   T ∼ 10^6.3 K  . (ii) The X-ray luminosity of the shock is strongly dependent on the gas fraction; it is  ≳10³⁹ erg s⁻¹  for halo gas fractions larger than 10 per cent. (iii) The hot diffuse gas in the simulation produces X-ray luminosities as large as  10⁴² erg s⁻¹  . This contributes to the total X-ray background in the Universe. (iv) We find an analytic fit to the maximum X-ray luminosity of the shock as a function of merger parameters. This fit can be used in semi-analytic recipes of galaxy formation to estimate the total X-ray emission from shocks in merging galaxies. (v) ∼10–20 per cent of the initial gas mass is unbound from the galaxies for equal-mass mergers, while 3–5 per cent of the gas mass is released for the 3:1 and 10:1 mergers. This unbound gas ends up far from the galaxy and can be a feasible mechanism to enrich the intergalactic medium with metals.     

ROSAT ASCA observations of X-ray luminous starburst galaxies: NGC 3310 and 3690

We present ROSAT [High Resolution Imager (HRI) and Position Sensitive Proportional Counter (PSPC)] and ASCA observations of the two luminous ( L _x ∼ 10⁴¹⁻⁴² erg s⁻¹) star-forming galaxies NGC 3310 and 3690. The HRI shows clearly that the sources are extended with the X-ray emission in NGC 3690 coming from at least three regions. The combined 0.1–10 keV spectrum of NGC 3310 can be described by two components, a Raymond–Smith plasma with temperature kT  = 0.81^+0.09_−0.12 keV and a hard power law, Γ = 1.44^−0.20_−0.11 (or alternatively a harder Raymond–Smith plasma with kT  ∼ 15 keV), while there is no substantial excess absorption above the Galactic column value. The soft component emission is probably a super wind while the nature of the hard emission is more uncertain with the likely origins being X-ray binaries, inverse Compton scattering of infrared photons, an active galactic nucleus or a very hot gas component (∼10⁸ K). The spectrum of NGC 3690 is similar, with kT  = 0.83^+0.02_−0.04 keV and Γ = 1.56^+0.11_−0.11. We also employ more complicated models such as a multi-temperature thermal plasma, a non-equilibrium ionization code or the addition of a third softer component, which improve the fit but not at a statistically significant level (2σ). These results are similar to recent results on the archetypal star-forming galaxies M82 and NGC 253.     

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.     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Radiation from early black holes – I. Effects on the neutral intergalactic medium

In the pre-reionization Universe, the regions of the intergalactic medium (IGM) which are far from luminous sources are the last to undergo reionization. Until then, they should be scarcely affected by stellar radiation; instead, the X-ray emission from an early black hole (BH) population can have much larger influence. We investigate the effects of such emission, looking at a number of BH model populations (differing for the cosmological density evolution of BHs, the BH properties, and the spectral energy distribution of the BH emission). We find that BH radiation can easily heat the IGM to  10³–10⁴ K  , while achieving partial ionization. The most interesting consequence of this heating is that BHs are expected to induce a 21-cm signal (  δ T _b∼ 20–30 mK  at   z ≲ 12  ) which should be observable with forthcoming experiments (e.g. LOFAR). We also find that at   z ≲ 10  BH emission strongly increases the critical mass separating star-forming and non-star-forming haloes.     

A variable Quasi-Periodic Oscillation in M82 X-1. Timing and spectral analysis of XMM–Newton and RossiXTE observations

We report results from a spectral and timing analysis of M82 X-1, one of the brightest known ultraluminous X-ray sources. Data from a new 105-ks XMM–Newton observation of M82 X-1, performed in 2004 April, and of archival RossiXTE observations are presented. A very soft thermal component is present in the XMM spectrum. Although it is not possible to rule out a residual contamination from the host galaxy, modelling it with a standard accretion disc would imply a black hole (BH) mass of  ≈10³ M_⊙  . An emission line was also detected at an energy typical for fluorescent Fe emission. The power density spectrum of the XMM observation shows a variable Quasi-Periodic Oscillation (QPO) at frequency of 113 mHz with properties similar to those discovered by Strohmayer and Mushotzky. The QPO was also found in seven archival RXTE observations, that include those analysed by Strohmayer and Mushotzky, and Fiorito and Titarchuk. A comparison of the properties of this QPO with those of the various types of QPOs observed in Galactic black hole candidates strongly suggests an association with the type-C, low-frequency QPOs. Scaling the frequency inversely to the BH mass, the observed QPO frequency range (from 50 to 166 mHz) would yield a BH mass anywhere in the interval few tens to  1000  M_⊙  .     

Chandra High-Resolution Camera observations of the luminous X-ray source in the starburst galaxy M82

We analyse Chandra High Resolution Camera observations of the starburst galaxy M82, concentrating on the most luminous X-ray source. We find a position for the source of         (J2000) with a 1 σ radial error of 0.7 arcsec. The accurate X-ray position shows that the luminous source is neither at the dynamical centre of M82 nor coincident with any suggested radio AGN candidate. The source is highly variable between observations, which suggests that it is a compact object and not a supernova or remnant. There is no significant short-term variability within the observations. Dynamical friction and the off-centre position place an upper bound of 10⁵–10⁶ M_⊙ on the mass of the object, depending on its age. The X-ray luminosity suggests a compact object mass of at least 500 M_⊙. Thus the luminous source in M82 may represent a new class of compact object with a mass intermediate between those of stellar-mass black hole candidates and supermassive black holes.     

Supermassive Black Holes in Early-Type Galaxies: Relationship with Radio Emission and Constraints on the Black Hole Mass Function

Using recently published estimates — based on high spatial resolution spectroscopy — of the mass M _BH of nuclear black holes for a sample of nearby galaxies, we explore the dependence of galaxy nucleus emissivity at various wavelengths on M _BH. We confirm an almost linear scaling of the black hole mass with the baryonic mass of the host spheroidal galaxy. A remarkably tight relationship is also found with both nuclear and total radio centimetric flux, with a very steep dependence of the radio flux on M _BH ( P  ∝  M ^2.5_BH). The high-frequency radio power is thus a very good tracer of a supermassive black hole, and a good estimator of its mass. This, together with the lack of significant correlations with the low-energy X-ray and far-IR flux, supports the view that advection-dominated accretion is ruling the energy output in the low accretion rate regime. Using the tight dependence of total radio power on M _BH and the rich statistics of radio emission of galaxies, we derive an estimate of the mass function of remnants in the nearby Universe. This is compared with current models of quasar and active galactic nucleus (AGN) activity and of the origin of the hard X-ray background (HXRB). As for the former, continuous long-lived AGN activity is excluded by the present data with high significance, whereas the assumption of a short-lived, possibly recurrent, activity pattern gives remarkable agreement. The presently estimated black hole mass function also implies that the HXRB has been produced by a numerous population (∼ 10⁻² Mpc⁻³) of moderately massive ( M _BH ∼ 10⁷ M⊙) black holes.     

The growth of supermassive black holes in pseudo-bulges, classical bulges and elliptical galaxies

Using results from structural analysis of a sample of nearly 1000 local galaxies from the Sloan Digital Sky Survey, we estimate how the mass in central black holes is distributed amongst elliptical galaxies, classical bulges and pseudo-bulges, and investigate the relation between their stellar masses and central stellar velocity dispersion σ. Assuming a single relation between elliptical galaxy/bulge mass, M _Bulge, and central black hole mass, M _BH, we find that  55⁺⁸₋₄  per cent of the mass in black holes in the local universe is in the centres of elliptical galaxies,  41⁺⁴₋₂  per cent in classical bulges and  4^+0.9_−0.4  per cent in pseudo-bulges. We find that ellipticals, classical bulges and pseudo-bulges follow different relations between their stellar masses and σ, and the most significant offset occurs for pseudo-bulges in barred galaxies. This structural dissimilarity leads to discrepant black hole masses if single   M _BH– M _Bulge  and   M _BH–σ  relations are used. Adopting relations from the literature, we find that the   M _BH–σ  relation yields an estimate of the total mass density in black holes that is roughly 55 per cent larger than if the   M _BH– M _Bulge  relation is used.     

Quasi-stars: accreting black holes inside massive envelopes

We study the structure and evolution of 'quasi-stars', accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate on to the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass,   M _*+ M _BH≈ M _*  . This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasi-stars. We show that the photospheric temperature of the envelope scales as   T _ph∝ M ^−2/5_BH M ^7/20_*  , and decreases with time while the black hole mass increases. Once   T _ph < 10⁴ K  , the photospheric opacity drops precipitously and T _ph hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities, this limiting temperature is approximately 4000 K. After a quasi-star reaches this limiting temperature, it is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 10³ and  10⁴ M_⊙  could form via this mechanism in less than a few Myr.     

Fields and filaments in the core of the Centaurus cluster

We present high-resolution images of the Faraday rotation measure (RM) structure of the radio galaxy PKS 1246−410 at the centre of the Centaurus cluster. Comparison with Hα-line and soft X-ray emission reveals a correspondence between the line-emitting gas, the soft X-ray emitting gas, regions with an excess in the RM images and signs of depolarization. Magnetic field strengths of 25 μG, organized on scales of ∼1 kpc and intermixed with gas at a temperature of 5 × 10⁶ K with a density of ∼0.1 cm⁻³, can reproduce the observed RM excess, the depolarization and the observed X-ray surface brightness. This hot gas may be in pressure equilibrium with the optical line-emitting gas, but the magnetic field strength of 25 μG associated with the hot gas provides only 10 per cent of the thermal pressure and is therefore insufficient to account for the stability of the line-emitting filaments.     

The dynamical state of RX J1347.5−1145 from a combined strong lensing and X-ray analysis

We perform a combined X-ray and strong lensing analysis of RX J1347.5−1145, one of the most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong lensing alone, based on published Very Large Telescope (VLT) and new Hubble Space Telescope ( HST ) data, strongly argues in favour of a complex structure. The analysis takes into account arc positions, shapes and orientations, and is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal mass distribution, with a total projected mass of   M _tot= (9.9 ± 0.3) × 10¹⁴ M_⊙  h ⁻¹  within a radius of 360 kpc  h ⁻¹ (1.5 arcmin). Such a complex structure could be a signature of a recent major merger as further supported by X-ray data. A temperature map of the cluster, based on deep Chandra observations, reveals a hot front located between the first main component and an X-ray emitting south-eastern subclump. The map also unveils a filament of cold gas in the innermost regions of the cluster, most probably a cooling wake caused by the motion of the cD inside the cool core region. A merger scenario in the plane of the sky between two dark matter subclumps is consistent with both our lensing and X-ray analyses, and can explain previous discrepancies with mass estimates based on the virial theorem.     

The BH mass of nearby QSOs: a comparison of the bulge luminosity and virial methods

We report on the analysis of the photometric and spectroscopic properties of a sample of 29 low-redshift  ( z < 0.6)  QSOs for which both Hubble Space Telescope ( HST ) WFPC2 images and ultraviolet HST   FOS spectra are available. For each object we measure the R -band absolute magnitude of the host galaxy, the C  iv (1550 Å) linewidth and the 1350 Å continuum luminosity. From these quantities we can estimate the black hole (BH) mass through the   M _BH– L _bulge  relation for inactive galaxies, and from the virial method based on the kinematics of the regions emitting the broad-lines. The comparison of the masses derived from the two methods yields information on the geometry of the gas emitting regions bound to the massive BH. The cumulative distribution of the linewidths is consistent with that produced by matter laying in planes with inclinations uniformly distributed between ∼10° and ∼50°, which corresponds to a geometrical factor   f ∼ 1.3  . Our results are compared with those of the literature and discussed within the unified model of active galactic nuclei.     

The host galaxy of a narrow-line Seyfert 1 galaxy, RE J1034+396, with X-ray quasi-periodic oscillations

Using simple stellar population synthesis, we model the bulge stellar contribution in the optical spectrum of a narrow-line Seyfert 1 galaxy, RE J1034+396. We find that its bulge stellar velocity dispersion is  67.7 ± 8 km s⁻¹  . The supermassive black hole (SMBH) mass is about  (1–4) × 10⁶ M_⊙  if it follows the well-known   M _BH–σ_*  relation found in quiescent galaxies. We also derive the SMBH mass from the Hβ second moment, which is consistent with that from its bulge stellar velocity dispersion. The SMBH mass of (1–4)  × 10⁶ M_⊙  implies that the X-ray quasi-periodic oscillation (QPO) of RE J1034+396 can be scaled to a high-frequency QPO at 27–108 Hz found in Galactic black hole binaries with a  10-M_⊙  black hole. With the mass distribution in different age stellar populations, we find that the mean specific star formation rate (SSFR) over the past 0.1 Gyr is  0.0163 ± 0.0011  Gyr⁻¹, the stellar mass in the logarithm is  10.155 ± 0.06  in units of solar mass and the current star formation rate is  0.23 ± 0.016 M_⊙ yr⁻¹  . For RE J1034+396, there is no relation between the Eddington ratio and the SSFR as suggested by Chen et al., despite a larger scatter in their relation. We also suggest that about 7.0 per cent of the total Hα luminosity and 50 per cent of the total [O  ii ] luminosity come from the star formation process.     

Unresolved emission and ionized gas in the bulge of M31

We study the origin of unresolved X-ray emission from the bulge of M31 based on archival Chandra and XMM–Newton observations. We demonstrate that three different components are present. (i) Broad-band emission from a large number of faint sources – mainly accreting white dwarfs and active binaries, associated with the old stellar population, similar to the Galactic ridge X-ray emission of the Milky Way. The X-ray to K -band luminosity ratios are compatible with those for the Milky Way and for M32; in the 2–10 keV band, the ratio is  (3.6 ± 0.2) × 10²⁷ erg s⁻¹ L⁻¹_⊙  . (ii) Soft emission from ionized gas with a temperature of about ∼300 eV and a mass of  ∼2 × 10⁶ M_⊙  . The gas distribution is significantly extended along the minor axis of the galaxy, suggesting that it may be outflowing in the direction perpendicular to the galactic disc. The mass and energy supply from evolved stars and Type Ia supernovae is sufficient to sustain the outflow. We also detect a shadow cast on the gas emission by spiral arms and the 10-kpc star-forming ring, confirming significant extent of the gas in the 'vertical' direction. (iii) Hard extended emission from spiral arms, most likely associated with young stellar objects and young stars located in the star-forming regions. The   L _X/SFR  (star formation rate) ratio equals  ∼9 × 10³⁸ (erg s⁻¹)(M_⊙ yr⁻¹)⁻¹  , which is about ∼1/3 of the high-mass X-ray binary contribution, determined earlier from Chandra observations of other nearby galaxies.