%Finished 7.10.00 \documentclass[12pt]{article} \usepackage{latexsym} \usepackage{amsmath} % \usepackage{amssym} \usepackage{graphicx} \setlength{\textwidth}{6.5truein} \setlength{\oddsidemargin}{0truein} \setlength{\topmargin}{-0.7truein} \setlength{\textheight}{9.5truein} \setlength{\parindent}{0.0truein} \setlength{\parskip}{0.15truein} \def\lesssim{\mathrel{\hbox{\rlap{\hbox{\lower5pt\hbox{$\sim$}}}\hbox{$<$}}}} \def\gtrsim{\mathrel{\hbox{\rlap{\hbox{\lower5pt\hbox{$\sim$}}}\hbox{$>$}}}} \def\ni{\noindent} \def\bfn{\mbox{\boldmath$\nabla$}} \begin{document} \addtocounter{section}{17} \addtocounter{subsection}{2} \addtocounter{figure}{2} \addtocounter{page}{118} \subsection{Stable and Unstable orbits around a Schwarzchild black hole} In Newtonian dynamics the equation of motion of a particle in a central potential is $${1 \over 2} \left ({dr \over dt}\right )^2 + V(r) = E, $$ where $V(r)$ is an ``effective potential''. For an orbit around a point mass, the effective potential is $$V(r) = {h^2 \over 2r^2} - {GM \over r} , $$ where $h$ is the specific angular momentum of the particle. The concept of an effective potential is, I am sure, familiar to you. If you sketch the effective potential, you can see easily that bound orbits have two turning points and that a circular orbit corresponds to the special case where the particle sits at the minimum of the effective potential. \begin{figure}[h] \begin{center} {\scalebox{0.80}[0.80]{\includegraphics{L11fig2.eps}}} \end{center} \caption{The Newtonian effective potential showing how an angular momentum barrier prevents particles reaching $r=0$.} \end{figure} Furthermore, you can see that in Newtonian dynamics, a finite angular momentum provides an {\it angular momentum barrier} preventing a particle reaching $r=0$. This is not true in General Relativity. In GR, the equation of motion of a particle around a central mass point is $$ \beta(r) \left ({dr \over dp}\right )^2 + {J^2 \over r^2} - {1 \over \alpha(r)} = - E, \eqno(1) $$ where $$ \alpha (r) = {1 \over \beta (r)} = \left (1-{2\mu \over r}\right ). $$ Since $d\tau^2 = E dp^2$ we can rewrite equation (1) as $$ \left ({dr \over d\tau}\right )^2 + V(r) = {1 \over E} \eqno (2) $$ where $h = J/\sqrt{E}$ and $V(r)$ is an effective potential $$ V(r) = \left (1- { 2\mu \over r} \right ) \left ( 1+ {h^2 \over r^2} \right ). \eqno(3) $$ If we transform to the dimensionless variables \begin{eqnarray*} x & = & {2 r \over r_s} = {rc^2 \over GM}, \\ h^{\prime 2} & =& { 4 h^2 \over r_s^2}, \end{eqnarray*} we can write the effective potential as $$ V(x) = \left ( 1- { 2 \over x} \right ) \left ( 1+ { h^{\prime 2} \over x^2} \right ) = 1 - { 2 \over x} - { 2 h^{\prime 2} \over x^3} + {h^{\prime 2} \over x^2}. \eqno (4) $$ Differentiate this expression, $${dV \over dx} = { 2 \over x^2} + {6h^{\prime 2} \over x^4} - {2h^{\prime 2} \over x^3}, $$ and so the extrema of the effective potential are located at the solutions of the quadratic equation $$x^2 - h^{\prime 2} x + 3h^{\prime 2} = 0,$$ {\it i.e.} at $$x = {h^\prime \over 2} \left \{ h^\prime \pm \sqrt{h^{\prime 2} - 12} \right \}. $$ IF $h^\prime = \sqrt{12} = 2\sqrt{3}$ {\it there is only one extremum}, and there are no turning points in the orbit for lower values of $h^\prime$. Figure 4 shows the effective potential for several values of $h^\prime$. The dots show the locations of stable circular orbits. The maxima in the potential are the locations of {\it unstable} circular orbits. \begin{figure}[h] \begin{center} \rotatebox{-90}{\scalebox{0.75}[0.75]{\includegraphics{pg_pot.ps}}} \end{center} \caption{The effective potential (4) plotted for several values of the angular momentum parameter $h^\prime$.} \end{figure} What is the physical significance of this result? The smallest stable circular orbit has $$x_{\rm min} = 6,\ \ {\it i.e.} \ \ r_{\rm min} = 6 {GM \over c^2}.$$ Gas in an accretion disc settles into circular orbits around the compact object as shown schematically in Figure 2. However, the gas slowly loses angular momentum because of turbulent viscosity (the turbulence is thought to be generated by magnetohydrodynamic instabilities). As the gas loses angular momentum it moves slowly towards the black hole, gaining gravitational potential energy and heating up. Eventually it loses enough angular momentum that it can no longer follow a stable circular orbit and so it falls into the black hole. We can, therefore, make a rough energy of the efficiency of energy radiation in an accretion disc. The maximum efficiency is of order the gravitational binding energy at the smallest stable circular orbit divided by the rest mass enery of the gas $$\epsilon_{\rm acc} \approx {1 \over 2} { G M m \over r_{\rm min}} {1 \over mc^2} \simeq {1 \over 12} \sim 8\%.$$ An accretion disc can convert perhaps a few percent of the rest mass energy of the gas into radiation. Compare this with the efficiency of nuclear burning of hydrogen to helium ($26$ MeV per He nucleus), $$\epsilon _{\rm nuclear} \sim 0.7 \%. $$ Accretion discs are capable of converting rest mass energy into radiation with an efficiency that is about 10 times greater than the efficiency of nuclear burning of hydrogen. The `accretion power' of black holes cause some of the most energetic phenomena known in the Universe. \subsection{Supermassive black holes} The first quasar\footnote {{\it Quasi-stellar radio sources}. We now know that the majority of quasars are radio quiet, and so they are often called QSOs for {\it quasi-stellar object}. } (3C273) was discovered in $1963$ by Maarten Schmidt. He measured a cosmological redshift of $z = 0.15$ for this object -- unprecedently high at the time. (Quasars have since been discovered with redshifts as high as $z = 5.8$) Quasars are very luminous, typically $100$-- $1000$ times brighter than a large galaxy. However, they are {\it compact}, so compact in fact that quasars look like stars in photographs. In fact, from variability and other studies one can infer that the size of the continuum emitting region of a quasar is of order a few parsecs. How can we explain such a phenomenon? Imagine an object radiating many times the luminosity of an entire galaxy from a region smaller than the Solar System. Donald Lynden-Bell was one of the first to suggest that the quasar phenomenon is caused by accretion of gas on to a {\it supermassive} black hole residing at the centre of a galaxy. The black hole masses required to explain the high luminosities of quasars are truly spectacular -- we require black holes with a few million to a few billion times the mass of the Sun. Do such supermassive black holes exist? The evidence in recent years has become extremely strong. Using the Hubble Space Telescope it is possible to probe the velocity dispersions of stars in the very central regions of galaxies. According to Newtonian dynamics, we would expect the characteristic velocities to vary as $$ v^2 \sim {GM \over r}. \eqno (5) $$ If the central mass is dominated by a supermassive black hole, then we expect the typical velocities of stars to {\it increas} as we go closer to the centre. This is indeed what is found in a number of galaxies. From the rate of increase of the velocities with radius, we can estimate the mass of the central object which seems to be correlated with the mass of the bulge component of the galaxy: $$ M_{bh} \approx 0.006 M_{bulge}. $$ It seems as though at the time of galaxy formation, about half a percent of the mass of the bulge collapses right to the very centre of a galaxy to form a supermassive black hole. During this phase the infall gas radiates efficiently producing a quasar. When the gas supply is used up, the quasar quickly fades away leaving a dormant massive black hole that is starved of fuel. Nobody has yet developed a convincing theory of how this happens, or of what determined the masses of the central black holes. \begin{center} \begin{tabular}{rlll} \multicolumn{4}{c}{\Large{Table 2: Supermassive Black Holes}}\\ \hline & & \\ & &$M_{bh}/M_\odot$ & Evidence\\ $\ast\ast\ast$ & M87 & $2 \times 10^9$ & stars \& optical disc\\ $\ast\ast$ & NGC 3115 & $1 \times 10^9$ & stars \\ $\ast\ast$ & NGC 4594 (Sombrero) & $5 \times 10^8$ & stars \\ $\ast\ast$ & NGC 3377 & $1 \times 10^9$ & stars \\ $\ast\ast\ast\ast\ast$ & NGC 4258 & $4 \times 10^7$ & masing $H_2 0$ disc \\ $\ast\ast$ & M31 (Andromeda) & $3 \times 10^7$ & stars \\ $\ast\ast$ & M32 & $3 \times 10^6$ & stars \\ $\ast\ast\ast\ast$ & Galactic Centre & $2.5 \times 10^6$ & stars \& 3D motions\\ \hline \end{tabular} \end{center} A skeptic might argue that these observations merely prove that a dense compact object exists at the centres of galaxies, not necessarily a black hole. But there are two beautiful observational results that probe compact objects on parsec scales -- making it almost certain that the central objects are black holes. In our own Milky Way Galaxy it is possible to measure the {\it proper motions} of stars in the galactic centre (using infrared wavelengths to penetrate through the dense dust that obscures optical light). This has allowed astronomers to measure actually see the stars moving and so infer their three dimensional motions. These observations imply that there exists a black hole of mass $2.5 \times 10^6 M_\odot$ at the centre of our Galaxy. In a remarkable set of observations, a disc of $H_2 O$ masars has been detected in the galaxy NGC 4258 using VLBI. The VLBI observations measure the velocities of the masing clouds on scales of $\sim 0.3$ -- $2$ parsec and are well fitted by a thin (actually slightly warped) disc in circular motion. The mass of the central black hole is estimated to be $4 \times 10^7 M_\odot$. \begin{figure}[h] \begin{center} \rotatebox{0}{\scalebox{0.55}[0.55]{\includegraphics{ngc4258.ps}}} \end{center} \caption{The masing $H_20$ disc in the centre of NGC4258.} \end{figure} Table 2 lists the masses of some supermassive black holes. I have given the observations a $5$ star rating. The masing disc of NGC 4528 gets a full five stars -- this is the strongest observational evidence for a supermassive black hole. The stellar motions in the Galactic Centre get four stars, though some astronomers might argue that this evidence is so strong that it should rate five stars. Most of the other observations are based on measurements of stellar velocity dispersions. This is fairly strong evidence, but not completely convincing\footnote {The interpretation of velocity dispersion measurements requires some assumptions about the degree of velocity anisotropy.} and so rate only two stars. Finally, to end this section, I will show you some recent results from X-ray spectroscopy. The following picture shows spectral line of iron at X-ray wavelengths measured in the Seyfert galaxy MCG-6-30-15. The iron lines come from the inner parts of the accretion disc around the black hole and the line profile allows one to probe the strong gravitational regime. The lines are assymetric as a result of special relativistic beaming and general relativity. The detailed shape of the line profile depends on the metric and hence on whether the black hole is rotating (see the notes on the Kerr metric). The hope is that the line profiles can be measured and modelled in detail and so be used to infer the {\it angular momenta} of black holes as well as their masses. \begin{figure}[h] \begin{center} {\scalebox{0.80}[0.80]{\includegraphics{mcg6_prof.cps}}} \end{center} \caption{ The line profile of iron K-alpha from MCG-6-30-15 observed by the ASCA satellite (Tanaka et al, 1995, Nature, 375, 659). The emission line is extremely broad, with a width indicating velocities of order one-third of the speed of light. There is a marked asymmetry towards energies lower than the rest-energy of the emission line (6.4 keV). This asymmetry is most likely caused by gravitational and relativistic-Doppler shifts near the black hole at the center of the galaxy. The solid line shows the model profile expected from a disk of matter orbiting the hole, extending between 3 and 10 Schwarzschild radii.} \end{figure} \end{document}