%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{page}{113} \section*{Lecture Seventeen: Theory of Relativity} \subsection{Collapse of a Massive Star to form a Black Hole} Do black holes exist in Nature? A star like the Sun is held up by radiation pressure. It is burning hydrogen to make helium releasing about $26\; {\rm Mev}$ for each atom of He that is formed. What happens when all of the nuclear fuel is used up? The star must collapse to a high density. In fact, we expect that the Sun will collapse to a form a white dwarf -- a star with a radius of about $5000\; {\rm km}$ and a spectacularly high mean density of about $10^{6}\; {\rm g}\; {\rm cm}^{-3}$. Astronomers knew about white dwarfs a long time ago. For example, in 1915 Adams discovered that the companion of the bright star Sirius had a `white' spectrum from which he deduced an effective temperature of $8000$ K. The mass of the companion (known as Sirius B) was estimated to be about $0.8 M_\odot$ from observations of the binary orbit. Given its observed flux and temperature, the blackbody radiation law could be used to infer a radius of $\sim 20,000 \; {\rm km}$ and so, given its mass, there could be no escape from the conclusion that Sirius B was a very compact star of extraordinarily high density. In 1925, Adams showed further that the spectral lines of Sirius B were gravitationally redshifted -- as expected from General Relativity. Here is what Eddington said about this discovery in his classic book {\it The Internal Constitution of the Stars} written in 1926: \begin{quotation} ``Prof. Adams has killed two birds with one stone; he has carried out a new test of Einstein's general theory of relativity and he has confirmed our suspicion that matter 2000 times denser than platinum is not only possible, but is present in the Universe.'' \end{quotation} This, incidentally, is a good example of why astronomy is so interesting. Astronomy turns up strange objects in which conditions are much more extreme than anything we encounter here on Earth. This allows us to study physics in unusual regimes, where we can hope to learn something fundamentally new . Examples include neutron stars, quasars, relativistic jets, accretion discs, gamma ray bursts and cosmological fluctuations in the early Universe. Nobody knew how to explain white dwarfs. What holds them up? What provides the pressure? The answer had to await the development of quantum mechanics and the formulation of Fermi-Dirac statistics. Fowler realised in 1926 that white dwarfs were held up by electron degeneracy pressure. The electrons in a white dwarf behave like the free electrons in a metal. Because of the Pauli exclusion principle, the electrons completely fill phase space up to a characteristic Fermi-energy. It is the Pauli exclusion principle that holds up a white dwarf. Chandrasekhar, in 1930, was on his way to Cambridge to start work as a research student under Eddington's supervision. In transit, Chandrasekhar realized that the more massive a white dwarf, the denser it must be and so the stronger the gravitational field. For white dwarfs about a critical mass of $1.4M_\odot$ (now called the Chandrasekhar limit), gravity would overwhelm degeneracy pressure and no stable solution was possible. The white dwarf would collapse to a point. Eddington was not impressed (though Chandrasekhar later won the Nobel prize for this work). Here is another quote from Eddington, \begin{quotation} ``The star apparantly has to go on radiating and radiating and contracting and contracting until, I suppose, it gets down to a few kilometres radius when gravity becomes strong enough to hold the radiation and the star at last can find peace.... I think that there should be a law of Nature to prevent the star from behaving in this absurd way.'' \end{quotation} After the discovery of the neutron, people realized that at extremely high densities the electrons would react with the protons to form neutrons. A new stable configuration for a star was possible -- neutron stars. A neutron star of one solar mass would have a radius of only $10\;{\rm km}$. But as with white dwarfs, neutron stars have a maximum mass above which no stable configuration is possible. This maximum mass is believed to be about $3 M_\odot$ or less -- the exact value is uncertain because of uncertainties in the equation of state of matter at such high densities. Nevertheless, we now believe that Eddington was wrong -- there is no law of Nature to prevent the gravitational collapse of a massive star. We believe that massive stars collapse to form {\it black holes}. \begin{figure}[h] \begin{center} \scalebox{0.8}[0.8]{\includegraphics{L11fig1.eps}} \end{center} \caption{Collapse of a star to form a black hole.} \end{figure} Figure 1 shows a space time diagram of a collapsing star. Eventually, the radius of the star collapses to less than the Schwarzchild radius, at which point the formation of a black hole becomes inevitable. Some theorists were very skeptical about the formation of black holes. The Schwarzchild solution is very special -- it is exactly spherically symmetric by construction. In reality, a star will not be perfectly symmetric and so perhaps, as it collapses, the assymmetries amplify and avoid the formation of an event horizon. Perhaps black holes never form or are exceedingly rare in Nature. Penrose in the early $1960$'s applied global geometrical techniques to prove a famous series of `singularity theorems'. These show that in realistic situations an event horizon (a closed trapped surface) will be formed and that there must exist a singularity within this surface, {\it i.e.} a point at which the curvature diverges and General Relativity ceases to be valid. I will touch on Penrose's approach in a later lecture. The singularity theorems were important in convincing people that black holes must form in Nature. In the rest of this section, I will look at some of the observational evidence for the existence of black holes. As we will see, there is compelling evidence that black holes do indeed exist. Furthermore, it should become possible within the next few years not only to measure the masses of black holes, but also to measure their angular momenta using powerful X-ray telescopes! Direct experimental probes of the strong gravity regime are now possible. \subsection{Compact binary systems} Although radiation cannot escape from black holes, one of the best way of finding candidate black holes is to search for luminous compact X-ray sources. The reason for this is that if a black hole has a stellar companion, the intense tidal field can pull gas from the companion producing an accretion disc around the black hole. A schematic picture is shown below. I will show you that accretion discs can radiate very efficiently. \begin{figure}[h] \begin{center} \scalebox{0.7}[0.7]{\includegraphics{artcv.ps}} \end{center} \caption{Schematic picture of a compact binary system.} \end{figure} The following table summarizes the common classes of compact binaries. The compact object can be a white dwarf, neutron star, or a black hole. \begin{center} \begin{tabular}{|l||l|l|l|}%[here] \hline \multicolumn{4}{|c|}{}\\ \multicolumn{4}{|c|}{\bf{{\Large Compact accreting binary systems}}}\\ \multicolumn{4}{|c|}{}\\ \hline & \multicolumn{3}{c|}{}\\ & \multicolumn{3}{c|}{{\large Compact object}} \\ \cline{2-4} % & \multicolumn{3}{c|}{}\\ {{\large Companion Star}} & White Dwarf& Neutron Star & Black Hole \\ \hline Early type, massive & None known & Massive X-ray binaries & Cyg X-I, LMC X-3\\ & & & \\ Late-type, low mass & Cataclysmic variables & Low mass X-ray binaries & A0620-00\\ & (e.g., Dwarf novae) & & \\ \hline \end{tabular} \end{center} If you find a compact binary system, then you can set limits on the mass of the compact object from the dynamics of the binary orbit. If you find evidence for a compact object that is more massive than the Chandrasekhar limit, then you have good evidence that the object might be a black hole. In fact it is not so straightforward. What observers actually measure is the {\it mass function} $$ f(M)={P K^3 \over 2\pi G}, $$ where $P$ is the orbital period, and $K$ is the radial velocity amplitude. For example, for the low mass X-ray binary A0620-00, the period is $P=7.7$ hours and $K = 457\; {\rm km}\;{s}^{-1}$. From Kepler's laws we can show that the mass function is related to the masses of the compact object ($M_1$), the companion star ($M_2$) and the inclination angle $i$ of the orbit to the plane of the sky by $$ f(M)={M_1^3 \sin^3 i \over ( M_1 + M_2)^2}. $$ You can see from this equation that the mass function is a strict lower limit on the mass of the compact object. It is equal to the mass, $f = M_1$, {\it only} if $M_2 = 0$ and the orbit is viewed edge on ($\sin i = 1$). For example, for A0620-00, the lower limit on the mass of the compact object is $2.9 M_\odot$, and this makes it a very good black hole candidate because this mass limit is very close to the theoretical {\it upper} limit for the mass of a neutron star. In fact, it is possible to make reasonable estimates \footnote{An estimate of the mass $M_2$ can be made by identified the spectral type and luminosity of the companion star and the inclination angle can be estimated from the shape of the stars light curve -- searching for evidence of eclipsing by the compact object.} for $M_2$ and $\sin i$ in this system leading to a probable mass of $\approx 10 M_\odot$ for the compact object -- well into the black hole regime. The following table from a review by Phil Charles summarizes the dynamical mass limits on some good black hole candidates (so called short X-ray transients). As you can see in several systems, V404 Cyg, G2000+25 and N Oph 77, the minimum mass inferred from the mass function is {\it well above} the theoretical maximum mass limit for a neutron star. As we understand things at present there can be no other explanation other than that the compact objects are black holes. \begin{table}[t] %\hrulefill %\medskip \centerline{\large Dynamical Mass Estimates of Binaries} \medskip \begin{center} \medskip \begin{tabular}{|lccccccl|} \cline{1-8} & & & & & & &\\ \multicolumn{8}{|c|}{Table 1. Derived Parameters and Dynamical Mass Measurements of SXTs}\\ & & & & & & &\\ \cline{1-8} & & & & & & &\\ {\it Source} & $f(M)$ & $\rho$ & $q$ & $i$ & $M_1$ & $M_2$ & \\ & $(M_\odot)$ & $({\rm g}\; {\rm cm}^{-3}$) & $(=M_1/M_2)$ & & $(M_\odot)$ & $(M_\odot)$ & Ref. \\ & & & & & & &\\ \cline{1-8} & & & & & & &\\ V404 Cyg & $6.08\pm0.06$ & $0.005$ & $17\pm1$ & $55\pm4$ & $12\pm2$ & $0.6$ & [1-2] \\ G2000+25 & $5.01\pm0.12$ & $1.6$ & $24\pm10$ & $56\pm15$ & $10\pm4$ & $0.5$ & [3-5]\\ N Oph 77 & $4.86\pm0.13$ & $0.7$ & $>19$ & $60\pm10$ & $6\pm2$ & $0.3$ & [6-9]\\ N Mus 91 & $3.01\pm0.15$ & $1.0$ & $8\pm2$ & $54^{+20}_{-15}$ &$6^{+5}_{-2}$ & $0.8$ & [13-15]\\ A0620-00 & $2.91\pm0.08$ & $1.8$ & $15\pm1$ & $37\pm5$ & $10\pm5$ & $0.6$ & [16-18]\\ J0422+32 & $1.21\pm0.06$ & $4.2$ & $>12$ & $20-40$ & $10\pm5$ & $0.3$ & [19-20]\\ & & & & & & &\\ J1655-40 & $3.24\pm0.14$ & $0.03$ & $3.6\pm0.9$ & $67\pm3$ & $6.9 \pm1$ & $2.1$ & [10-12]\\ 4U1543-47 & $0.22\pm0.02$ & $0.2$ & -- & $20-40$ & $5.0\pm2.5$ & $2.5$ & [21]\\ & & & & & & &\\ Cen X-4 & $0.21\pm0.08$ & $0.5$ & $5\pm1$ & $43\pm11$ & $1.3 \pm 0.6$ & $0.4$ & [22-23]\\ & & & & & & &\\ \cline{1-8} & & & & & & &\\ \multicolumn{8}{|c|}{\emph{References}}\\ & & & & & & &\\ \multicolumn{8}{|c|}{\parbox{15.5cm}{ [1] Casares \& Charles 1994; [2] Shahbaz et. al. 1994b; [3] Filippenko et. al. 1995a; \newline [4] Beekman et. al. 1996; [5] Harlaftis et. al. 1996; [6] Filippenko et. al. 1997; \newline [7] Remillard et. al. 1996; [8] Martin et. al. 1995; [9] Harlaftis et. al. 1997; \newline [10] Orosz \& Bailyn 1997; [11,\ 12] van der Hooft 1997, 1998; [13] Orosz et. al. 1996; [14] Casares et. al. 1997; [15] Shahbaz et. al. 1997; [16] Orosz et. al. 1994; \newline [17] Marsh et al 1994; [18] Shahbaz et. al. 1994a; [19] Filippenko et. al. 1995b; \newline [20] Beekman et. al. 1997; [21] Orosz et. al. 1998; [22] McClintock \& Remillard 1990; [23] Shahbaz et. al. 1993.}}\\ & & & & & & & \\ \cline{1-8} \end{tabular} \end{center} \end{table} \end{document} \newpage \subsection{Stable and Unstable orbits around a Schwarzchild black hole} In Newtonian dynamics $${1 \over 2} \left ({dr \over dt}\right )^2 + V(r) = E $$ $$V(r) = {h^2 \over 2r^2} - {GM \over r} \hspace{.25cm}\ldots \parbox{4cm}{``effective potential''}$$ Insert here Lecture 11 Figure 2 \vspace{13cm} A finite angular momentum provides an ``angular momentum barrier'' preventing a particle reaching $r=0$. This is not true in {\bf General Relativity}. In GR our equation sof motion gave \[ \left. \begin{array}{c} \beta(r) \left ({dr \over dp}\right )^2 - {J^2 \over r^2} - {1 \over \alpha(r)} = - E\\ \\ \alpha (r) = {1 \over \beta (r)} = \left (1-{2\mu \over r}\right ) \\ \\ d\tau^2 = E dp^2 \end{array} \right \} \] Hence \[ \left. \begin{array}{c} \left ({dr \over d\tau}\right )^2 + \left (1-{2\mu \over r}\right ) \left (1+{h^2 \over r^2}\right ) = - {1 \over E}\\ \\ h = j/\sqrt{E} \end{array} \right\} \] which I can write as $$\left ({dr \over d\tau}\right )^2 + v(r) = - {1 \over E}$$ where \[ \begin{array}{rcl} V(r) & = & (1-2\mu/r) \ (1+h^2/r^2) \\ & & \\ & = & (1-2/x) \ (1+h^{\prime 2}/x^2 ) \hspace{0.25cm} x = {rc^2 \over GM} \end{array} \] $$V(x) = 1 - 2/x - 2h^{\prime 2}/x^3 + h^{\prime 2} x^2$$ and so to differentiate with respect to $x$ $${dv \over dx} = 2/x^2 + {6h^{\prime 2} \over x^4} - 2h^{\prime 2}/x^3$$ and so extrema of the ``effective potential'' $v(x)$ are located at $$x^4 - h^{\prime 2} x + 3h^{\prime 2} = 0$$ $$x = {h^{\prime 2} \pm \sqrt{h^{\prime 4} - 12h^{\prime 2}} \over 2} = {h^\prime 2 \over 2} \left \{ 1 \pm \sqrt{h^{\prime 2} - 12} \right \}$$ IF $h^\prime = \sqrt{12} = 2\sqrt{3}$ {\bf there is only one extremum}, and there are no turning points in the orbit for lower values of $h^\prime$. See the figure on the next page. The maxima in the potential are {\bf unstable circular orbits}. Insert lecture 11, figure 3 \vspace{7cm} The smallest stable orbit has $$X_{\rm min} = 6,\ \ r_{\rm min} = 6 {GM \over c^2}$$ so, the maximum efficiency of energy generation in an accretion disc is of order, $$\epsilon \approx {1 \over 2} G{M_m \over r_{\rm min}} {1 \over mc^2} \simeq {1 \over 12} \sim 8\%$$ Compare this with the efficiency of nuclear burning of hydrogen to helium $$\triangle E_{\rm nuclear} = 0.007 m_pc^2$$ $$\sim 0.7 \% $$ %now type big table, supermassive black holes \end{document} \begin{center} \begin{tabular}{rlll} \multicolumn{4}{c}{\Large{SUPERMASSIVE BLACK HOLES}}\\ & &$M_n/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 4528 & $4 \times 10^7$ & masing $10_29$ 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\\ \end{tabular} \end{center} \end{document}