Abstract

There is a lot of current astrophysical evidence and interest in intermediate mass black holes (IMBH), ranging from a few hundred to several thousand solar masses. The active galaxy M82 and the globular cluster G1 in M31, for example, are known to host such objects. Here, we discuss several aspects of IMBH such as their expected luminosity, spectral nature of radiation, and associated jets. We also discuss possible scenarios for their formation including the effects of dynamical friction, and gravitational radiation. We also consider their formation in the early universe and also discuss the possibility of supermassive black holes forming from mergers of several IMBH and compare the relevant time scales involved with other scenarios.

1. Introduction

Intermediate mass black holes (IMBH) are those black holes having mass between that of stellar black holes and supermassive black holes, that is, in the range of 500 to 104. Recent observations indicate an intermediate mass black hole in the elliptical galaxy NGC 4472, with an X-ray luminosity of 4 × 1032 J/s [1].

The Schwarzschild radius corresponding to these masses is of the order of

An accretion disc is formed by material falling into a gravitational source. Conservation of angular momentum requires that, as a large cloud of material collapses inward, any small rotation it may have will increase. Centrifugal force causes the rotating cloud to collapse into a disc, and tidal effects will tend to align this disc's rotation with the rotation of the gravitational source in the middle.

The Bondi accretion rate for a black hole (of mass M moving with a velocity in a medium of density), is given by [2]: where is the cross section and it is given by .

The velocity is given by where is the velocity of sound in the medium of density and it is given by , and is the ratio of specific heats and is the universal gas constant.

The equation describing the velocity of isothermal winds has many solutions depending on the initial conditions at the base of the wind. There is only one critical solution for which the velocity increases from subsonic at the base to supersonic far out. This velocity passes through the critical point and implies one particular value of the initial velocity at the lower boundary of the isothermal region.

If the density at this point is fixed, the mass loss rate is fixed by (2) as . The total energy increases from negative at the base of the wind to positive in the supersonic region, so the flow requires the input of energy into the wind. This energy input is needed to keep the flow isothermal and it is this energy that is transferred into kinetic energy of the wind through the gas pressure.

As we will see in the next section, for an intermediate mass black hole, the temperature is typically of the order of  K. In this case, the velocity of sound is of the order of  m/s.

Also we have ; hence, the accretion rate is given by

For a typical number density of , and the IMBH mass of , the accretion rate is: .

The above expression implies that the higher the density of the medium through which the black hole is travelling, the more the accretion rate. Also, the accretion rate is inversely related to the velocity with which the black hole is travelling in the medium. Although there is not much conclusive evidence for the IMBH, there is indirect evidence.

For a given mass of luminous object, there is a maximum value for its luminosity (Eddington luminosity), as the radiation pressure would tend to push the matter apart exceeding the gravitational force supporting it. This is the luminosity a body would have to have so that the force generated by radiation pressure exceeds the gravitational force. Thus, observed luminosity can set a lower limit on the mass of an accreting black hole.

Chandra and XMM-Newton observations in the nearby spiral galaxy have detected X-ray sources of luminosities of the order of , with the source away from the centre [3]. In Section 4, we will discuss the dynamics of the IMBH.

The X-ray luminosity corresponding to the accretion rate of  kg/s is given by . This matches with recent observed results [4].

The Eddington luminosity (for an IMBH of mass ) is given by:

where, the Thomson cross section is given by:

A possible candidate for the IMBH is the globular cluster G1 in M31, which is the most massive stellar cluster in the local group, with a mass of the order of .

Due to a large number of field stars contained within the accretion radius the Bondi accretion by an IMBH is complicated [4]. The accretion radius is given by

In the case of globular cluster G1 in M31, there are more than stars within 0.4 pc of the centre [5]. The dynamical effects of these stars should also be considered.

Many more such dense stellar clusters are known to exist. For instance, a recently recognised super star cluster in our own galaxy is the Westerlund I. Most of its estimated half a million stars are crowded in a region hardly 3 parsecs wide. Several dozen of these stars are among the most massive and luminous superhot Wolf-Rayet stars, LBV, red supergiants, yellow hypergiants, and so forth. Also, near the Milky Way centre, we have the well-known Arches and Quintuplet clusters.

In the centre of M15 there are approximately five million stars per cubic parsec, having hundred million times more stellar density than the solar neighbourhood. Also galaxies like M31, M33, and the Milky Way itself have comparable central stellar densities. However M32 (satellite galaxy of Andromeda) has thirty million stars in a cubic parsec at its core. Even HST cannot resolve individual stars in this region.

More and more massive binary stars are being found, like, for instance, WR20a (>80  each, with 3.7-day period). As discussed above, IMBH is likely to form in dense clusters containing young massive stars. IMBH is believed to power ultraluminous X-ray sources (ULXs). IMBH can capture companion stars, which provide accreting material to sustain ULX. These stars may be blue giants, white dwarfs, and so forth.

Tidal forces can rip giants and the material can fall into the IMBH. To rip apart these stars the mass of the black hole is around the mass of the IMBH [6]. The energy released during this process will be the binding energy of the stars which is given by

Several ULXs are identified in the Antenna, a pair of colliding galaxies, producing several stars in dense cluster. Stellar collisions lead to formation of so-called megastars (~103), which collapses on short time scales to form IMBH. Clusters in M82 have such stellar densities. Binary and multiple IMBH can also form.

2. Black Body Considerations of IMBH

We have already obtained the Eddington luminosity as

By considering the black hole emission to obey black body radiation, we can use Stefan’s law to relate the luminosity to the temperature, , as where is Stefan’s constant and area of the black hole and it is given by where is the Schwarzschild radius and indicates the size of the ambient gas around the black hole.

Making use of this in (10) we get

This implies that the temperature of the IMBH (accretion disk) of mass of the order of 500 to  MΘ and ranging from 10 to 100 is

From Wien's law the corresponding wavelength is given by

The corresponding wavelength is

This lies in the soft X-ray region of the spectrum [7].

During the earlier epochs of the universe, the density was much larger; hence, the ambient density is also larger by the same factor of .

The present density of the universe is of the order of one proton per cubic metre. As we will see in Section 5, the maximum redshift up to which we can detect supermassive black hole is of the order of .

The number density at this epoch is given by

The temperature corresponding to this redshift and this number density is of the order of

And the corresponding wavelength is of the order of .

This falls in the X-ray region of the spectrum.

This wavelength is further red-shifted by a factor of . Hence, the observed wavelength will be of the order of .

This lies in the UV region of the spectrum.

3. Jets from the Black Hole

One of the manifestations of this accretion energy release is the production of so-called jets: the collimated beams of matter that are expelled from the innermost regions of accretion disks. These jets shine particularly brightly at radio frequencies. In rotating black holes, the matter forms a disk due to the mechanical forces present. In a Schwarzschild black hole, the matter would be drawn in equally from all directions and thus would form an omnidirectional accretion cloud rather than disk. Jets form in Kerr black holes (rotating black holes) that have an accretion disk [8].

Black holes convert a specific fraction of accretion energy into radiation, which is traced by the X-ray luminosity and jet kinetic energy, which is traced by the radio-emission luminosity [9]. The matter is funnelled into a disk-shaped torus by the black hole's spin and magnetic fields, but, in the very narrow regions over the black hole's poles, matter can be energized to extremely high temperatures and speeds, escaping the black hole in the form of high-speed jets. Inferred jet velocities close to the speed of light suggest that jets are formed within a few gravitational radii of the event horizon of the black hole.

The horizon for the Kerr black hole is given by

Here, is the geometric mass and is the geometric angular momentum.

From the condition that should be real, the limiting case is given by .

From this, the maximum angular momentum is given by

From the classical expression for the angular momentum associated with a jet of length , assuming the particles to be travelling at near speed of light, the expression becomes

Considering a conical jet with base radius r and density , the mass of the jet is given by

Then the angular momentum becomes

From the geometry of the jet, we can relate the length of the jet to the radius r as . Here, we have assumed the small opening angle of the jet to be .

Length of the jet is

For typical densities of the ambient gas and for the IMBH of mass , the length is of the order of 20 pc.

4. Evolution of a Star Cluster and Possible Scenario for IMBH Formation

One of the possible models for the formation of an IMBH is the collapse of a cluster of stars [10]. The collapsed core can accrete matter ejected during the formation. If a collection of a thousand 10 solar mass stars in a volume of a parsec cube collapses, ejecting 30% of its mass, then the total mass of the ambient gas is given by

The accretion rate is given by (4) as

If 30% of the mass of each star is ejected, then this implies an accretion rate of . And the corresponding (Eddington) luminosity is

Dynamical friction is related to loss of momentum and kinetic energy of moving bodies through a gravitational interaction with surrounding matter in space. The effect must exist if the principle of conservation of energy and momentum is valid since any gravitational interaction between two or more bodies corresponds to elastic collisions between those bodies. For example, when a heavy body moves through a cloud of lighter bodies, the gravitational interaction between this body and the lighter bodies causes the lighter bodies to accelerate and gain momentum and kinetic energy.

Since energy and momentum are conserved, this body has to lose a part of its momentum and energy equal to the sums of all momenta and energies gained by the light bodies. Because of the loss of momentum and kinetic energy of the body under consideration, the effect is called dynamical friction. Of course the mechanism works the same way for all masses of interacting bodies and for any relative velocities between them.

However, while in the above case the most probable outcome is the loss of momentum and energy by the body under consideration, in the general case it might be either loss or gain. In a case when the body under consideration is gaining momentum and energy, the same physical mechanism is called sling effect.

The full Chandrasekhar dynamical friction formula for the change in velocity of the object involves integrating over the phase space density of the field of matter [11]. By assuming a constant density, though, a simplified equation for the force from dynamical friction, , is given as where is the gravitational constant, is the mass of the moving object, is the density, and is the velocity of the object in the frame in which the surrounding matter was initially at rest.

In this equation, is not a constant but depends on how compares to the velocity dispersion of the surrounding matter. The greater is the density of the surrounding media, the stronger will be the force from dynamical friction. Similarly, the force is proportional to the square of the mass of the object. The force is also proportional to the inverse square of the velocity. This means the fractional rate of energy loss drops rapidly at high velocities.

Dynamical friction is, therefore, unimportant for objects that move relativistically, such as photons. Dynamical friction is particularly important in the formation of planetary systems and interactions between galaxies.

During the formation of planetary systems, dynamical friction between the protoplanet and the protoplanetary disk causes energy to be transferred from the protoplanet to the disk. This results in the inward migration of the protoplanet.

When galaxies interact through collisions, dynamical friction between stars causes matter to sink toward the centre of the galaxy and for the orbits of stars to be randomised. The dynamical friction comes into effect in the evolution of cluster between the ambient gas and dust and the central IMBH.

If an IMBH of mass is moving with a velocity of in a uniform background of ‘‘fixed’’ lighter stars of equal masses . Then as the IMBH moves, a star approaching with impact parameter , will have a velocity change given by: where is the encounter duration and is the acceleration. They are given by

The change in velocity becomes

The kinetic energy gained by the star corresponding to this change in velocity is given by

If is the number density of the stars, then the number of encounters with impact parameter between and is given by

The total change in velocity is given by where , the lower bound on , is given for the case where the gravitational energy is of the order of the kinetic energy of the black hole. That is,

And the upper limit is the size of the system. Let .

Then the total change in velocity of the black hole is given by

In the above discussion we have assumed that the stars are stationary. This need not be true. If the stars have a velocity dispersion of and the black hole is moving at a slower velocity than this, that is, , then the dynamical friction force is given by where

In the case of the velocity of the black hole being faster than , the dynamical friction will be reduced and the black hole will slide through the cluster.

The black hole is also subjected to a gravitational force due to the star cluster, which is given by where is the density distribution of stars and is the typical mass of the star (assuming all stars are of the same mass).

For a constant density, , the potential is given by where is the radius of the cluster.

The gravitational force on the black hole is given by

(The particle inside a homogenous gravitational system performs simple harmonic motion!)

The equation of motion of the black hole in the star cluster is given by

This corresponds to a damped oscillator. And the solution is given by [12]

The black hole undergoes a damped oscillation in the star cluster. The damping time corresponding to the system is given by where

In the case of the system M82, which harbours an IMBH of mass in the range of 500 to  MΘ, the black hole is not found at the centre but displaced by about one kilo-parsec from the centre.

The period corresponding to the oscillation is given by

From (37), we can work out the time taken for the BH to shift by this distance. Using this equation along with (32) and (35), we get the time of the order of years.

For the system under consideration, the number density of the stars in the cluster is , typical mass of the star in the cluster is about one solar mass, and the velocity dispersion is of the order of  km/s.

We get the damping time for the system as

This works out to be of the order of years.

Taking the effects of dynamical friction into consideration, the relaxation time for the system to form the IMBH is given by

For a denser core, that is, about , the time taken to form the IMBH is given by

5. Formation of Supermassive Black Hole (SBH): Merger of IMBH

Zwart et al. [13] propound that dense star densities near galactic centres can lead to runaway stellar mergers thus efficiently producing IMBH. Indeed they suggest that the inner most ten parsecs of our galaxy can contain about 50 IMBH each of about thousand solar masses. These can sink towards the core as they interact with the stars. Every time an IMBH has a stellar encounter, the stars’ velocity is boosted. This reduces the potential energy making the IMBH sink to the galactic core. These IMBH can thus ultimately merge with the galaxy’s supermassive black hole.

However if there are not enough stars very near to the galaxy’s SBH, then the IMBH may stop falling inward, halting at about 0.01 light years from the centre. However if there are several IMBH near the galactic core, they can interact with one another and once every million years an IMBH can merge with the SBH.

It has been suggested that the presence of several IMBH near the galactic centre might explain why there are so many clusters of young stars (like S2, the B type star) in the galactic core, where tidal forces should rip apart the gas clouds from which the stars form.

These stars could have formed much further away and the IMBH could have shepherded them inwards, much like the thin rings of Uranus or Saturn that are kept in place by being “shepherded” by tiny satellites.

This IMBH can merge with the stars in the surrounding volume to give a supermassive black hole in the time scale given by

This is the time scale for relaxing in a cluster of objects of average mass moving with a mean velocity of . If is the average cross section for interaction, gives the “radius” of the sphere of influence around a given individual object and the mean free path is , and is the number density of objects. Then the “collision” or “interaction” time scale is given by

The term comes from many body effects (i.e., , the so-called Coulomb logarithm in plasma physics). This leads to (48). Note the sharp dependence as well as .

The maximum redshift observed till now is of the order of . The age of the universe corresponding to this redshift is given by

The above equation is the standard formula for the age of the universe, corresponding to any redshift . Higher implies younger objects. So corresponds to ~109 years. Now, recent supermassive black holes are claimed to have been observed at , hardly ~5 × 108 years after the big bang [14].

We estimate the time to form the black hole from the accretion time scale, corresponding to Eddington luminosity. This gives a logarithmic time scale of years. These are typical time scales for growth of massive black holes in the early universe.

Essentially we have (with as the initial mass) accretion power , with , is a numerical parameter, and integrating we get . This time can be translated into corresponding redshift.

According to the model suggested above, the time taken for the formation of SBH is of the order of years. The corresponding redshift is of the order of . For redshifts above this limit, as per this model, we should not be able to detect any supermassive black holes.

Other possible ways in which SBH can form are the following. The first is by the merger of two or more IMBH. Moving masses like black holes produce gravitational waves in the fabric of space-time. A more massive moving object will produce more powerful waves, and objects that move very quickly will produce more waves over a certain time period. The second is by accretion of matter by the IMBH in systems such as AGNs and quasars.

Gravitational waves are usually produced in an interaction between two or more compact masses. Such interactions include the binary orbit of two black holes orbiting each other. As the black holes orbit each other, they send out waves of gravitational radiation that reaches the Earth; however, once the waves do get to the Earth, they are extremely weak.

This is because gravitational waves decrease in strength as they move away from the source. Even though they are weak, the waves can travel unobstructed within the fabric of space-time.

From Kepler’s third law, the period is related to the separation by

Knowing the period, we can determine the orbital velocity from .

Power lost by gravitational waves (quadrupole formula for gravitational waves) is given by (for objects of mass with separation of )

This emission results in the objects coming closer. The merger time is obtained by integrating the rate of binding energy change as follows:

This implies a merger time for two objects separated by an initial distance of

With the constant, , for equal mass objects.

So a larger initial separation would imply a longer merger time. So for the merger time to be years or less, the initial separation can be precisely calculated.

For two IMBH of mass  MΘ each, to merge in years, the distance of separation should be of the order of where

And the energy loss due to the gravitational waves emission is given by

Therefore, the distance of separation is given by

The corresponding orbital frequency is given by

The orbital velocity is of the order of

This implies that the frequency is  Hz, and hence the period will be given by .

For two IMBH with separation of about parsecs, the time taken to merge is of the order of Hubble time. For such a model to produce an SBH, we need about IMBH merging together.

Hence, the model discussed earlier, with the IMBH merging with the surrounding stars, gives a much more efficient way of generating a supermassive black hole.

In the case of accretion of matter by IMBH to form SBH, the increase in mass is exponential with time where years is the characteristic time required for the mass to increase -fold, with the accreting disk emitting at maximum luminosity.

For the IMBH to accrete enough matter to become SBH of mass say , we have

The corresponding time scale is of the order of

This implies that the mass will have to increase -fold by accretion for the IMBH to become an SBH. Hence, even this model does not provide an efficient way of formation of SBH from an IMBH.

IMBH would have formed in the early universe. Owing to low metal content, the earliest stars would have been very massive, a few hundred solar masses [15]. Such stars would end up in a pair-instability supernova (around oxygen-neon burning temperature of 2 billion degrees) and would collapse into a black hole (if their mass exceeds 250 solar masses).

Oxygen-neon burning occurs at ~2 × 109 K. At this stage, temperatures are high enough for the electron-positron pair production processes . At ~6 × 109 K, it is maximal. Now if the nature of the black body radiation changes from photons to electron-positron pairs (which are fermions) the energy density is no longer but . So once the pairs are forming, the radiation pressure drops by . So there is less support against gravitational collapse or gravitational contraction.

The balance between the two (or the “stability condition”) can be expressed as

Now,  K,  m (appropriate density for oxygen-neon burning); this implies a limiting mass of 250 solar masses, beyond which the gravitational energy density is higher. This result in the context of population III supermassive stars is discussed by Heger and Woosley [16].

Signatures of such explosions of supermassive stars at (when the universe was only half a billion years old and eleven times smaller) could be sought with future space telescopes [17, 18].

Observations of galaxy core show correlation between black hole masses and the spheroidal (bulge) component of the host galaxy. Thus, the primordial low mass galaxies (blue galaxies) would host low mass central black holes, which just correspond to the IMBH mass (about to the galaxy mass). There is a well-known relation between black hole mass and spheroidal bulge component mass of the host galaxy [19].

Typically, the black hole mass is the bulge mass. Primordial galaxies (so-called blue galaxies as seen, e.g., in Hubble deep field) have a lower mass ~108 solar mass. This would give a central black hole mass of ~104 solar mass, corresponding to a typical IMBH.

Merger of these primeval galaxies would lead to larger galaxies and cluster and the IMBH would also merge and sink to the core forming an SBH (like, e.g., M87 has a black hole).

6. Concluding Remarks

We have discussed several aspects of the expected characteristics of IMBH like their luminosity, accretion rate, formation, and so forth. They are likely to have formed in the early universe, around . This redshift corresponds to a minimal time scale for formation (and growth) of a massive black hole. Beyond , the mass would not be large enough, as not enough time would have elapsed for the growth of black holes by accretion or merger. Supermassive black holes are not likely to form from merger of IMBH. Other possible scenarios are also discussed.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.