The formation of galaxies has been the subject of intensive theoretical and observational work during the past few decades. Observational data reveal the existence of supermassive black holes lurking at the centers of galaxies formed at a period very close to the recombination era. However, current observations remain incapable of providing a decisive answer as to the origin of these supermassive black holes. Recent development in the detection techniques of gravitational waves may provide a clue to the processes leading to the formation of supermassive black holes at the center of galaxies. Gravitational waves emanating from these processes in the early universe are extremely weak. Nevertheless, future generations of gravitational wave detectors may shed light on the formation of supermassive black holes in the early universe.
The demography of black holes at the center of galaxies is a promising channel for a better understanding of galaxy formation . The Hubble Space Telescope and the Chandra X-ray Observatory detected supermassive black holes with masses in excess of one billion solar masses in quasars at redshifts corresponding to only a few hundred million years after the Big Bang . The existence of supermassive black holes imposes important constraints on the formation mechanism of galaxies . Furthermore, the gas physics involved in the formation of supermassive black holes is not fully understood yet . So far there is no decisive answer concerning the formation of the earliest black holes, primarily because their growth process masks the origin and properties of the initial progenitor. The logical question that might arise is to what extent modern astronomical observations can set the stage for probing the evolutionary phase of galaxy formation. In this paper, we review some theoretical and observational studies of black holes in the early universe, their possible origin, formation, and fate. In the first part we outline the observational evidence connecting the existence of supermassive black holes in the early universe, followed by a brief introduction to the technique employed in the determination of their masses. In the second part, we present the various models for the formation of supermassive black holes in the early universe, including the collapse of population III stars, dynamical instabilities, the collapse of gas due to dynamical instabilities, the collapse of supermassive stars, the dynamical processes in enriched halos, and the formation and fate of primordial black holes.
2. Masses of Supermassive Black Holes
The determination of the masses of black holes is an essential requirement in verifying the degree of accuracy of various galactic formation and black hole growth models and the extent of their adherence to observations. Two important approaches are used to estimate the masses of black holes lurking at the centers of galaxies. The first technique is known as reverberation mapping. It involves the measurement of the structure of the broad emission line region (BLR) around a supermassive black hole. The equation describing this process is given by 
where G is the gravitational constant, is the rms velocity of gas moving near the black hole as determined from the Doppler broadening of the gaseous emission lines, is the radius of the broad line region, and f is a form factor that depends on the shape of the BLR. The measurement of is considered a serious challenge  , and the measurement of the f factor is also difficult .
The other technique to estimate the masses of supermassive black holes is the use of the M-sigma relation. It represents the correlation between the mass of the supermassive black hole and the velocity dispersion of its host galaxy’s bulge . It is used to estimate black hole masses in faraway galaxies, and is given by:
where A is a constant of order 3, σ is the stellar velocity dispersion of the galaxy’s bulge, and α is a constant of order 5 representing the slope of the M-σ relation. Ferrarese and Merritt  found A = 3.1 and α = 4.8 ± 0.5. Other studies gave values very close to the above results . The tight nature of the M-σ relation suggests that a feedback mechanism is operating between the growth of supermassive black holes and the growth of the galaxy’s bulge.
3. Processes in Black Hole Formation
In this section we introduce some basic processes leading to the formation of supermassive black holes. Accretion around a black hole is one of the key ingredients of its growth. The second element is the existence of a dark matter halo that forms and grows from primordial density fluctuations characterized by a virial radius, a mass over-density, and a virial temperature. This halo serves as a host for the pre-galactic disk, which usually grows via gas dynamic processes.
3.1. Accretion around Supermassive Black Holes
Accretion around supermassive black holes is considered the only mechanism capable of producing the observed luminosities produced by supermassive black holes in quasars . If accretion is an acceptable mechanism for black hole growth, the progenitor remains controversial, with primordial black holes  , dark stars  , and collapsing clouds of gas  being the main candidates. The existence of an efficient mechanism for transporting angular momentum outward will enable the accretion material to approach a marginally stable orbit. The existence of magnetic fields in the matter flowing into the disk, as well as turbulent motions, is one such mechanism, since it leads to the transfer of angular momentum outward . The Eddington accretion rate is a characteristic scale for accretion, and is given by 
The total energy released in the disk is equal to the Eddington luminosity
This is a critical luminosity for any given mass M, beyond which the radiation force overcomes gravity. Luminosities ranging from 1042 to 1048 erg/s have been observed for active galactic nuclei (AGNs), corresponding to black hole masses ranging from 105 to 109 solar masses . Even though a good agreement has been found between the observation and the theory of the spectral distribution of radiation, these theories are mainly concerned with mass accretion rates and the luminosity of the accretion disk irrespective of the origin of the accreting supermassive black hole.
3.2. Primordial Dark Matter Halos
Galaxies are thought to be formed from baryonic matter in dark matter halos born out of small primordial density fluctuations . There are three important parameters that can be inferred from these halos. The first is the virial mass Mvir that can be calculated directly from the virial theorem. The second is the circular velocity Vc which can be calculated from the relation
where rvir is the virial radius. The third parameter is the virial temperature, which is given by
where mp is the proton mass, μ is the mean molecular weight, and kB is the Boltzmann constant. The gravitational collapse of the baryon component can proceed when the mass of the over dense region reaches the Jeans mass MJ. At masses in excess of MJ, baryons are captured and are then shock-heated by the subsequent collapse and virilization of dark matter. Gas dynamics processes predict that low-mass objects are less efficient in dissipating energy and cool rather slowly, whereas more massive objects can cool at a faster rate . The collapsing halos in the early universe exhibit a virial temperature smaller than 104 K and are referred to as mini-halos. A necessary condition for the gas to cool down and form the first stars is that the halos should rely on the less efficient H2 cooling .
4. Formation of Supermassive Black Holes
During the past few decades, several models have been proposed to explain the presence of massive black holes (MBHs) at redshifts corresponding to the era when the universe was less than one billion years old. Important questions to answer are when did the seeding black holes at the center of galaxies form and what mechanism was involved in their growth? Several possible formation channels have been investigated to understand the MBH seed, as outlined below.
4.1. Collapse of Population III Stars
Population III stars are massive metal-free objects comprising the first generation of stars after the Big Bang. These stars are postulated to have formed in mini-halos with masses of the order 106 solar masses and to have collapsed from the highest primordial density field. For Tvir > 103 K, the cooling process is mediated by molecular hydrogen . Atomic hydrogen cooling takes place in the larger halos with a total mass of 108 solar masses and . Simulation of the collapse of molecular clouds suggests that massive stars with M > 100 solar masses can form . The fate of population III stars depends primarily on their masses. The collapse of 40 - 140 solar masses low metallicity stars is predicted to directly form a black hole. When the mass of the population III star is in the range 140 to 260 solar masses, the fate of the star is determined by the electron-positron pair production instability that leads to supernova explosions. Supernovae predicted by this model for certain ranges of massive stars will release a colossal amount of energy that can be detected by current observatories. No such events have been recorded so far. On the other hand, the above model has large uncertainties concerning the final mass of the population III stars.
4.2. Gas Dynamic Instabilities
Metal-free or metal-poor proto-galaxies are efficient nurseries where black holes can for and grow. In these systems, supermassive black holes can also be formed directly out of a dense gas cloud . On the other hand, enriched halos exhibit an efficient cooling process which favors fragmentation and star formation rather than direct black hole formation. In metal-free gas clouds that characterize the very first proto-galaxies, the collapse is expected to occur only in massive halos with virial temperatures , where the formation of molecular hydrogen is inhibited . At these temperatures H2 formation is inhibited and atomic hydrogen becomes an efficient agent for cooling down the tenuous gas until it reaches 4000 K . At , the line-trapping of Lyman-α photons in isothermally collapsing gas causes the equation of state to stiffen with the consequence that fragmentation becomes harder to achieve provided that the metallicity does not exceed about 10−4 of the solar metallicity . The dissociation of H2 in these systems is brought about by Lyman-α trapping. In such halos, gas cooling and contraction proceed gradually with no fragmentation until rotational support halts the collapse, which usually occurs before reaching densities that allow the formation of a massive black hole (MBH).
Local, rather than global, instabilities in a self-gravitating galactic disk can be calculated using the Toomre stability parameter formalism. The Toomre parameter Q is defined as
where Σ is the surface mass density, cs is the speed of sound, and is the epicyclical frequency, and V is the circular velocity of the disk. Gravitational instabilities occur when Q approaches a critical value QC. Instabilities might lead to mass in-fall rather than fragmentation and star formation, provided that the destabilization of the system is kept below a threshold value. This happens when the inflow rate is below a critical threshold
where αc is the viscosity parameter. This process continues until the mass accumulated at the center (Ma) is enough to make the disk marginally stable. The upper limit of the mass that can contribute to MBH formation is determined by the mass and spin parameter of the halo.
4.3. Collapse of Supermassive Stars
Gas dynamical processes can also lead to the formation of supermassive stars (SMS) that may collapse, under certain conditions, to form an MBH. Gas accumulated at the few parsecs around the center of the halo, by processes described in the previous section, can reach 104 to 106 solar masses. For efficient gas accumulation, an SMS ( ) may form, which eventually collapses to form a black hole . In systems where mass accumulation is fast enough, the outer layers of the SMS are not thermally relaxed during much of the lifetime of a main sequence star . These stars exhibit complex structures with a convective core surrounded by a convectively stable envelope containing most of the star’s mass. Hydrogen burning in the core of these stars is relatively low and continues throughout most of its massive stages. When hydrogen is exhausted, the SMS will contract and suffer catastrophic neutrino losses that lead to its collapse to an initial black hole with a mass of a few solar masses that grows subsequently via accretion from the resulting bloated envelope. This object is referred to as a quasistar  , and it consists of a low-mass central black hole surrounded by a massive radiation-pressure-supported envelope. The black hole grows gradually at the expense of the massive envelope until the resulting MBH is unveiled.
4.4. Dynamical Processes in Enriched Halos
Star formation can proceed in mini-halos characterized by a virial temperature Tvir < 104 K  . The halos will be enriched with metals by the first generation of population III stars, and thus fragmentation and formation of low mass stars will be a natural outcome of this enrichment . This process sets the stage for new horizons of MBH formation. Stellar dynamical processes may lead to the formation of compact star clusters   , resulting from collisions. These collisions arise from dynamical interactions and may play a major role in the formation of very massive stars (VMS) leading to the formation of MBH remnants in the range 102 - 104 solar masses . In an attempt to reach equilibrium, the compact core cluster initially contracts and then starts to decouple thermally from its outer region. Energy transfer from the central dense core will cause a rapid core collapse . Dynamical friction causes a segregation of more massive stars in the center. If these massive stars remain in the main sequence stage, then a subsystem will be developed and will decouple from the cluster. In this subsystem, star collisions can proceed in a runaway manner eventually leading to the growth of VMSs . The fate of VMSs depends essentially on their metal enrichment. Metal enriched VMSs will lose much of their mass and end their life as less massive objects (≈150 MO) . The final fate is either a low-mass black hole or a pair-instability supernova. For low metallicity, VMSs may have a different fate. Stars with masses ≥ 40MO and sub-solar metallicity may collapse directly into a black hole without a supernova explosion.
4.5. Dark Stars in the Early Universe
Dark stars are a new line of research that proposes that the first stars in the universe were fueled by dark matter heating rather than by nuclear fusion . Weakly interacting massive particles (WIMPs) are considered among the best dark matter candidates . It is assumed that in the early universe the density of dark matter was sufficiently high to trigger dark matter annihilation . The annihilation products of WIMPs inside a star can be trapped to produce enough energy to heat its core and prevent its collapse. The first stars are postulated to form inside dark matter halos of masses of the order of 106 solar masses  , with one single star per halo. It is also argued that these stars set the stage for many important processes like reionization, the seeding of supermassive black holes, and the production of heavy elements in subsequent generations of stars. The lightest neutralino is motivated by supersymmetry (SUSY) arguments and is considered the best WIMP candidate in the Minimal Supersymmetric Standard Model .
In the dark star model, authors assumed a mass of 100 GeV for the annihilating WIMPS. So far, WIMPs in general and neutralinos in particular have not been detected despite intensive searches over the past few decades. Furthermore, no trace of supersymmetric particles has been found in the Large Hadron Collider (LHC), even though it attains energies of seven tera-electron volts, which is far in excess of the 100 GeV postulated for annihilating DM particles in dark stars.
5. Primordial Black Holes
Theoretical studies of the possibility of formation of primordial black holes (PBH) in the early universe date back to the original work of Hawking . He argued that extreme densities and inhomogeneities in the early universe can lead to the local collapse of matter resulting in the formation of black holes. More recently, Choptuik  and Kim  demonstrated the formation of PBHs in the inflationary era, during which the energy density of the universe experienced a dramatic decrease leading to a cosmological phase transition. Hawking  argues that PBHs formed in a wide spectrum of masses in the early universe ranging from 10−5 g, corresponding to the Planck mass, to 1017 solar masses. His upper limit for mass exceeds, by many orders of magnitude, even the greatest masses of supermassive black holes observed today in galactic centers. On the other hand, the formation of very small black holes may arise either from the softening of the equation of state, phase transitions, or from the collapse of hypothetical cosmic strings.
Over dense regions in the early universe may collapse to a black hole if the gravitational attraction overcomes the pressure forces and the velocity of expansion . This condition is fulfilled when the potential energy for self-gravitation
exceeds the kinetic energy of expansion
where R is the radius of a region in the early universe, and is the energy density. The units are such that . In a Friedman universe the sum of these energies is zero, therefore
Furthermore, Hawking assumed that the equation of state relating the pressure P and the energy density has the form , and that is proportional to . Thus,
A necessary condition for the collapse to occur is that the gravitational energy, , should exceed the internal energy U. Taking and , the condition for collapse to occur becomes
for , and , the condition for collapse reduces to
Once a black hole is formed, it will grow by accreting nearby matter. The rate of accretion was calculated by Zeldovitch and Navikov 
where here is the density of the background universe. But since (see above), then
where is the initial mass of the black hole and is the time of its formation. Thus, if is small compared to , that is, if the black hole is small compared to the particle horizon, then remains small and there would be almost no accretion. However, if is of the same order as , then the Zeldovitch-Navikov argument leads to . In this case, the accretion would cause the black hole to grow at the same rate as the particle horizon, producing black holes of the order of the Hubble radius if the growth continued to the present time, or it would reach a mass of 1015 to 1017 solar masses if the growth was at the same rate as that of the particle horizon.
6. Gravitational Waves from the Relic of Galaxy Formation
Gravitational waves may be produced extensively during the initial phase of galaxy formation in the early universe. Several models have been proposed to explain the formation of massive black holes at the center of galaxies at redshifts corresponding to the era when the universe was less than one billion years old. Masses of black holes at the center of faraway galaxies can be estimated using the M-sigma relation . An important question to answer is at what epoch the seeding black holes at the center of galaxies were formed, and what are the possible mechanisms involved in their growth. Many hypotheses have been advanced to elucidate the origin of the formation of massive or supermassive black holes at the center of galaxies. As mentioned earlier, several possible formation channels have been proposed for the formation of massive black holes seed  , collapse of supermassive stars  , merger of small and intermediate mass black holes. The collision of massive or supermassive black holes at the center of these galaxies will produce energetic GWs that may lie within a sphere of radius again corresponding to the sensitivity of the actual or future GWs detectors . As the sensitivity of gravitational wave detectors improves, they will be capable of detecting GW events at redshifts beyond the era of star formation, and it is highly probable that these detections will provide evidence to support that primordial black holes are involved  .
The formation of supermassive black holes at the centers of galaxies in the early universe is one of the most challenging problems in astrophysics. In this paper we introduced the various models of supermassive black hole formation and we compared these models with observational findings. This approach reveals a gap in our understanding of the most efficient processes leading to the formation of supermassive black holes. Gravitational waves may be associated with the formation of supermassive black holes that can be detected by future generations of gravitational waves detectors. Alternatively, primordial black holes could be considered as a good candidate for seeding new born galaxies in the very early universe and in an era when no supermassive black holes could have had the time to form. An analysis of the cosmic background radiation may give a clue about the existence of primordial black holes. Future missions such as the James Webb Space Telescope will reveal the structure of new born galaxies in the very early universe and will provide a wealth of information regarding the nature of the initial seeding black holes at the centers of galaxies.