Arxiv:0706.1988V1 [Physics.Ed-Ph] 14 Jun 2007 Ae H Ih Otae Ie Itne O Xml,1 It 3 Example, Time for the S Distance

Home , Arthur Geoffrey Walker

arXiv:0706.1988v1 [physics.ed-ph] 14 Jun 2007 ae h ih otae ie itne o xml,1 it 3 example, time For the s distance. 1 of given = terms a second in travel light to them light specify the we wetakes that involved and distances large center The so the larger. are vastly at as ourselves universe place the dra- view longer universe no the observa (see we of telescopic first tions: planet view on Galileo’s Our here after last changed we, center. matically the and its contained at planets), beyond self were was (the Earth, sphere universe move The a to 1). Fig. on appeared that fixed few where a for except beyond motionless. Far appear stars system. the solar planets, our the within phenomena motions local merely by are stars. caused these of that the know background and we the course shrinks against Of on and move and planets grows star, and itself polar Moon Moon the the across around times, rotates drift longer sky clouds the that Moon, know We the universe. changeless a codn oteacetcsooia eif h stars, the belief, cosmological ancient the to According of impression strong a provides sky night the at look A I.1 eeta pee facetcosmology. ancient of spheres Celestial 1: FIG. eateto hsc,Uiest fWisconsin-Milwaukee of University Physics, of Department .SASADGALAXIES AND STARS I. 07Srn eetr tU-iwue.Te r o td pu enrol study students for are to They aid UW-Milwaukee. an at as semester, prepared Spring were 2007 notes lecture These × 10 etrso srnm,Atohsc,adCosmology and Astrophysics, Astronomy, on Lectures 8 m / 3 = s × 10 8 300 = m usA Anchordoqui A. Luis , Dtd pig2007) Spring (Dated: 0 km, 000 - h nvra a fGravitation: of Law Universal the iclrobtaon h etro h aay(ftotal (of Galaxy the of center the mass around orbit circular a a ih iue=18 = minute light 1 h u n h oa ytm(fttlmass that total (of and mass system Galaxy the solar the of the of most and center that Sun the the assume near so, center concentrated the do esti- To about is be Sun can Galaxy. the Galaxy the of the data of in orbital stars the the using all mated of mass total The hmi em flgt h at-ondsac is is distance distance Earth-Sun Earth-Moon The The ls. 1.28 specify 1 is light. could which we of we km, Moon, 384,000 but terms the kilometers, in and or them Sun meters the use to usually distances specifying For iue.Froti h oa ytm lt sabout is Pluto system, solar the in out 6 Far minutes. ilo er rs,s t pe is speed its so 250 so, every or once years approximately million center Galactic the orbits oteeg,sm 26 center some “nu- Galactic just edge, the be central from the to way bulging seems to half a which located is Sun, has star, Our It another arms. spiral ly. and 2,000 cleus” roughly for of (Greek ness Galaxy the call we dis- which great way”). to plane, “milky extending the a stars that of in disc suggested tances flat Wright a century was Thomas Way half Milky A 1750) stars. (about individual ob- later of first comprised numbers is Galileo Way countless Milky Way. long of the that the Milky telescope as his the well with as served as known seen, be strip can cloudy brightness of degrees ing the that system. us solar from the farther of times reach 10,000 Therefore, outer is away. ly star 4.3 nearest about the is Centauri, Proxima us, to . × 5 = u aayhsadaee f1000l n thick- a and ly 100,000 of diameter a has Galaxy Our vary- with stars of thousands nights, moonless clear On × 10 ..Bx43 iwue,W 30,USA 53201, WI Milwaukee, 413, Box P.O. , v 2 10 9 M /r y=2 = ly 1 mfo h u,o 6 or Sun, the from km 11 .Te,apyNwo’ Law, Newton’s apply Then, ). en h etiea ceeainand acceleration centripetal the being r100000k;ti seult . light 8.3 to equal is this km; 150,000,000 or m e n(G srnm 2,i the in 320, Astronomy (UG) in led v pssonly. rposes ≈ 9 = 0 km 200 = = 10 . . 2 46 998 13 . 5 × 2 × km × GMm × π , 10 10 0 yfo h etr h Sun The center. the from ly 000 r 10 10 . 2 / 15 8 26 s 6 8 . r3 yr m n ih year light 1 and km, m , m = 000 / × m 3 s . 156 10 × v r 2 − 10 . × 156 4 , 13 y h ers star nearest The ly. 10 × km 7 F s 10 / yr 7 = m s / oein move ) ma yr F with , being (3) (2) (1) 2 where G =6.67 10−11 N m2 kg−2 is Newton’s constant. method employs simple geometry to measure the parallax All in all, × of a star. By parallax we mean the apparent motion of a star against the background of more distant stars, due r v2 M = 2 1041 kg . (4) to Earth’s motion around the Sun. The sighting angle G ≈ × of a star relative to the plane of Earth’s orbit (usually indicated by θ) can be determined at different times of Assuming all the stars in the Galaxy are similar to our 30 the year. Since we know the distance d from the Earth Sun (M⊙ = 2 10 kg), we conclude that there are to the Sun, we can determine the distance D to the star. roughly 1011 stars× in the Galaxy. For example, if the angle θ of a given star is measured In addition to stars both within and outside the Milky to be 89.99994◦, the parallax angle is p φ =0.00006◦. Way, we can see with a telescope many faint cloudy From trigonometry, tan φ = d/D, and since≡ the distance patches in the sky which were once all referred to as “neb- to the Sun is d =1.5 108 km the distance to the star is ulae” (Latin for clouds). A few of these, such as those × in the constellations of Andromeda and Orion, can actu- d d 1.5 108 km D = = × =1.5 1014 km , (5) ally be discerned with the naked eye on a clear night. In tan φ ≈ φ 1 10−6 × the XVII and XVIII centuries, astronomers found that × these objects were getting in the way of the search for or about 15 ly. comets. In 1781, in order to provide a convenient list of Distances to stars are often specified in terms of paral- objects not to look at while hunting for comets, Charles lax angles given in seconds of arc: 1 second (1”) is 1/60 Messier published a celebrated catalogue. Nowadays as- of a minute (1’) of arc, which is 1/60 of a degree, so 1” tronomers still refer to the 103 objects in this catalog by = 1/3600 of a degree. The distance is then specified in their Messier numbers, e.g., the Andromeda Nebula is parsecs (meaning parallax angle in seconds of arc), where M31. the parsec is defined as 1/φ with φ in seconds. For ex- − ◦ Even in Messier’s time it was clear that these extended ample, if φ =6 10 5 , we would say the the star is at × objects are not all the same. Some are star clusters, a distance D =4.5 pc. It is easily seen that groups of stars which are so numerous that they ap- 1pc=3.26 ly=3.08 1016 m . (6) peared to be a cloud. Others are glowing clouds of gas × or dust and it is for these that we now mainly reserve Parallax can be used to determine the distance to stars as the word nebula. Most fascinating are those that be- far away as about 100 ly ( 30 pc) from Earth. Beyond long to a third category: they often have fairly regular that distance, parallax angles≈ are two small to measure elliptical shapes and seem to be a great distance beyond and more subtle techniques must be employed. the Galaxy. Immanuel Kant (about 1755) seems to have been the first to suggest that these latter might be circular discs, but appear elliptical because we see them at III. LUMINOSITY AND BRIGHTNESS an angle, and are faint because they are so distant. At first it was not universally accepted that these objects A useful parameter for a star or galaxy is its luminosity were extragalactic (i.e. outside our Galaxy). The very (or “absolute luminosity”), L, by which we mean the to- large telescopes constructed in the XX century revealed tal power radiated in watts. Also important is the appar- that individual stars could be resolved within these ex- ent brightness, l, defined as the power crossing unit area tragalactic objects and that many contain spiral arms. at the Earth perpendicular to the path of light. Given Edwin Hubble (1889-1953) did much of this observational that energy is conserved and ignoring any absorption in work in the 1920’s using the 2.5 m telescope on Mt. Wil- space, the total emitted power L when it reaches a dis- son near Los Angeles, California. Hubble demostrated tance D from the star will be spread over a sphere of that these objects were indeed extragalactic because of surface area 4πD2. If D is the distance from the star to their great distances [1]. The distance to our nearest spi- Earth, then ral galaxy, Andromeda, is over 2 million ly, a distance 20 times greater than the diameter of our Galaxy. It L =4πD2 l. (7) seemed logical that these nebulae must be galaxies similar to ours. Today it is thought that there are roughly Careful analyses of nearby stars have shown that the ab- 4 1010 galaxies in the observable universe – that is, as solute luminosity for most of the stars depends on the many× galaxies as there are stars in the Galaxy [2]. mass: the more massive the star, the greater the luminosity.

II. DISTANCE MEASUREMENTS BY PARALLAX IV. SURFACE TEMPERATURE

Last class we have been talking about the vast distance Another important parameter of a star is its surface of the objects in the universe. Today, we will discuss temperature, which can be determined from the spec- diﬀerent methods to estimate these distances. One basic trum of electromagnetic frequencies it emits. 3

The rate at which an object radiates energy has been found to be proportional to the fourth power of the Kelvin temperature T and to the area A of the emitting object, i.e., L AT 4. At normal temperatures ( 300 K) we are not aware∝ of this electromagnetic radiation≈ because of its low intensity. At higher temperatures, there is sufficient infrared radiation that we can feel heat if we are close to the object. At still higher temperatures (on the order of 1000 K), objects actually glow, such as a red-hot electric stove burner. At temperatures above 2000 K, objects glow with a yellow or whitish color, such as the filament of a lightbulb. Planck’s law of blackbody radiation predicts the spectral intensity of electromagnetic radiation at all wavelengths from a blackbody at temperature T . The spectral radiance or brightness (i.e., the energy per unit time per unit surface area per unit solid angle per unit frequency ν) is FIG. 2: Measured spectra of wavelengths emitted by a black- 2hν3 1 body at different temperatures. I(ν,T )= , (8) c2 ehν/kT 1 − where h =6.626 10−34 J s and k =1.38 10−23 JK−1 spectrum of the other peaks at about 350 nm (bluish). are the Planck and× Boltzmann constants× and c is the Using Wien’s law, the temperature of the reddish star is speed of light. The law is sometimes written in terms of 2.90 10−3 m K the spectral energy density T = × = 4140 K . (12) r 700 10−9 m × 4π u(ν,T )= I(ν,T ) , (9) The temperature of the bluish star will be double because c its peak wavelength is half; just to check which has units of energy per unit volume per unit fre- 2.90 10−3 m K quency. Integrated over frequency, this expression yields T = × = 8280 K . (13) b 350 10−9 m the total energy density. Using the relation λ = c/ν, × the spectral energy density can also be expressed as a For a blackbody the Steffan-Boltzmann equation reads function of the wavelength, L = σAT 4 , (14) 8πhc 1 u(λ, T )= . (10) −8 −2 −4 λ5 ehc/λkT 1 where σ = 5.67 10 W m K . Thus, the power − radiated per unit× of area from a star is proportional to Stars are fairly good approximations of blackbodies. the fourth power of the Kelvin temperature. Now the As one can see in Fig. 2, the 5500 K curve, correspond- temperature of the bluish star is double that of the redish ing to the temperature of the Sun, peaks in the visible star, so the bluish must radiate 16 times as much energy part of the spectrum (400 nm <λ< 750 nm). For lower per unit area. But we are given that they have the same temperatures the total radiation drops considerably and luminosity, so the surface area of the blue star must be the peak occurs at longer wavelengths. It is found exper- 1/16 that of the red one. Since the surface area is 4πr2, imentally that the wavelength at the peak of the spec- we conclude that the radius of the redish star is 4 times trum, λp, is related to the Kelvin temperature by larger than the radius of the bluish star (and its volume 64 times larger) [3]. λ T =2.9 10−3 mK . (11) p × This is known as Wien’s law. V. HR DIAGRAM We can now use Wien’s law and the Steffan-Boltzmann equation (power output or luminosity AT 4) to deter- An important astronomical discovery, made around mine the temperature and the relative size∝ of a star. Sup- 1900, was that for most of the stars, the color is related pose that the distance from Earth to two nearby stars to the absolute luminosity and therefore to the mass. can be reasonably estimated, and that their measured A useful way to present this relationship is by the so- apparent brightnesses suggest the two stars have about called Hertzsprung-Russell (HR) diagram [4]. On the the same absolute luminosity, L. The spectrum of one of HR diagram, the horizontal axis shows the temperature the stars peaks at about 700 nm (so it is reddish). The T , whereas the vertical axis the luminosity L, each star 4

A star on the main sequence of the HR diagram at this temperature has absolute luminosity of about L 1026 W. Then, using Eq. (7) we can estimate its distance≈ from us,

L D = r4πl 1026 W −12 −2 ≈ r4π 10 W m 3 1018 m , (16) ≈ × or equivalently 300 ly.

VII. STELLAR EVOLUTION

The stars appear unchanging. Night after night the heavens reveal no significant variations. Indeed, on human time scales, the vast majority of stars change very little. Consequently, we cannot follow any but the tini- est part of the life cycle of any given star since they live FIG. 3: Hertzsprung-Russell diagram. The vertical axis de- for ages vastly greater than ours. Nonetheless, in today’s picts the inherent brightness of a star, and the horizontal axis class we will follow the process of stellar evolution from the surface temperature increasing from right to left. The the birth to the death of a star, as we have theoretically spectral class is also indicated. You can see that the Sun is a reconstructed it. G-star. There is a general consensus that stars are born when gaseous clouds (mostly hydrogen) contract due to the pull of gravity. A huge gas cloud might fragment into is represented by a point on the diagram shown in Fig. 3. numerous contracting masses, each mass centered in an Most of the stars fall along the diagonal band termed the area where the density was only slightly greater than main sequence. Starting at the lowest right, we find the at nearby points. Once such “globules” formed, grav- coolest stars, redish in color; they are the least luminous ity would cause each to contract in towards its center- and therefore low in mass. Further up towards the left we of-mass. As the particles of such protostar accelerate in- find hotter and more luminous stars that are whitish like ward, their kinetic energy increases. When the kinetic en- our Sun. Still farther up we find more massive and more ergy is sufficiently high, the Coulomb repulsion between luminous stars, bluish in color. Stars that fall on this the positive charges is not strong enough to keep hydro- diagonal band are called main-sequence stars. There are gen nuclei appart, and nuclear fussion can take place. In also stars that fall above the main sequence. Above and a star like our Sun, the “burning” of hydrogen occurs to the right we find extremely large stars, with high lu- when four protons fuse to form a helium nucleus, with minosity but with low (redish) color temperature: these the release of γ rays, positrons and neutrinos.1 are called red giants. At the lower left, there are a few The energy output of our Sun is believed to be due stars of low luminosity but with high temperature: these principally to the following sequence of fusion reactions: are white dwarfs. 1H+1H 2H+2 e+ +2 ν (0.42 MeV) , (17) 1 1 →1 e

VI. DISTANCE TO A STAR USING HR 1H+2H 3He + γ (5.49 MeV) , (18) 1 1 →2 Suppose that a detailed study of a certain star suggests and that it most likely fits on the main sequence of the HR 3He +3He 4He +1H+1H (12.86 MeV) , (19) diagram. Its measured apparent brightness is l = 1 2 2 →2 1 1 10−12 W m−2, and the peak wavelength of its spectrum× is λ 600 nm. We can first find the temperature using p ≈ Wien’s law and then estimate the absolute luminosity 1 using the HR diagram; namely, The word “burn” is put in quotation marks because these high- temperature fusion reactions occur via a nuclear process, and must not be confused with ordinary burning in air, which is a −3 2.9 10 m K chemical reaction, occurring at the atomic level (and at a much T × 4800 K . (15) ≈ 600 10−9 m ≈ lower temperature). × 5 where the energy released for each reaction (given in evidence that stars are actually being born at this mo- parentheses) equals the difference in mass (times c2) be- ment in the Eagle Nebula. tween the initial and final states. Such a released energy As hydrogen fuses to form helium, the helium that is is carried off by the outgoing particles. The net effect formed is denser and tends to accumulate in the central of this sequence, which is called the pp-cycle, is for four core where it was formed. As the core of helium grows, 4 protons to combine to form one 2He nucleus, plus two hydrogen continues to fuse in a shell around it. When positrons, two neutrinos, and two gamma rays: much of the hydrogen within the core has been consumed, the production of energy decreases at the center and is no 1 4 + 4 1H 2He+2e +2νe +2γ . (20) longer sufficient to prevent the huge gravitational forces → from once again causing the core to contract and heat up. Note that it takes two of each of the first two reactions to 3 The hydrogen in the shell around the core then fuses even produce the two 2He for the third reaction. So the total more fiercely because of the rise in temperature, causing energy released for the net reaction is 24.7 MeV. How- the outer envelope of the star to expand and to cool. The ever, each of the two e+ quickly annihilates with an elec- 2 surface temperature thus reduced, produces a spectrum tron to produce 2mec = 1.02 MeV; so the total energy of light that peaks at longer wavelength (reddish). By released is 26.7 MeV. The first reaction, the formation this time the star has left the main sequence. It has of deuterium from two protons, has very low probability, become redder, and as it has grown in size, it has become and the infrequency of that reaction serves to limit the more luminous. Therefore, it will have moved to the right rate at which the Sun produces energy. These reactions 7 and upward on the HR diagram. As it moves upward, it requiere a temperature of about 10 K, corresponding to enters the red giant stage. This model then explains the an average kinetic energy (kT ) of 1 keV. origin of red giants as a natural step in stellar evolution. In more massive stars, it is more likely that the energy Our Sun, for example, has been on the main sequence output comes principally from the carbon (or CNO) cy- for about four and a half billion years. It will probably cle, which comprises the following sequence of reactions: remain there another 4 or 5 billion years. When our Sun leaves the main sequence, it is expected to grow in size 12C+1H 13N+ γ , (21) 6 1 → 7 (as it becomes a red giant) until it occupies all the volume out to roughly the present orbit of the planet Mercury. 13 13 + 7N 6C+ e + ν , (22) If the star is like our Sun, or larger, further fusion can → occur. As the star’s outer envelope expands, its core is shrinking and heating up. When the temperature reaches 13C+1H 14N+ γ , (23) 8 6 1 → 7 about 10 K, even helium nuclei, in spite of their greater charge and hence greater electrical repulsion, can then reach each other and undergo fusion. The reactions are 14N+1H 15O+ γ , (24) 7 1 → 8 4He +4He 8Be , (27) 2 2 →4 15O 15N+ e+ + ν , (25) 8 → 7 and 4 8 12 15 1 12 4 2He +4Be 6C , (28) 7N+1H 6C+2He . (26) → → with the emission of two γ-rays. The two reactions must It is easily seen that no carbon is consumed in this cycle occur in quick succession (because 8Be is very unstable), (see first and last equations) and that the net effect is 4 and the net effect is the same as the pp cycle. The theory of the pp cycle and the carbon cycle as the source of energy for the Sun and 4 12 3 2He 6C (7.3 MeV) . (29) the stars was first worked out by Hans Bethe in 1939 [5]. → The fusion reactions take place primarily in the core This fusion of helium causes a change in the star which of the star, where T is sufficiently high. (The surface moves rapidly to the “horizontal branch” of the HR di- 4 temperature is of course much lower, on the order of a agram. Further fusion reactions are possible, with 2He 12 16 few thousand K.) The tremendous release of energy in fusing with 6C to form 8O. In very massive stars, higher 20 24 these fusion reactions produces an outward pressure suf- Z elements like 10Ne or 12Mg can be made. This process ficient to halt the inward gravitational contraction; and of creating heavier nuclei from lighter ones (or by absorp- our protostar, now really a young star, stabilizes in the tion of neutrons at higher Z) is called nucleosynthesis. main sequence. Exactly where the star falls along the The final fate of the star depends on its mass. Stars can main sequence depends on its mass. The more massive lose mass as parts of their envelope drift off into space. the star, the further up (and to the left) it falls in the Stars born with a mass less than about 8 solar masses HR diagram. To reach the main sequence requires per- eventually end up with a residual mass less than about haps 30 million years, if it is a star like our Sun, and 1.4 solar masses, which is known as the Chandrasekhar it is expected to remain there 10 billion years (1010 yr). limit [6]. For them, no further fusion energy can be ob- Although most of stars are billions of years old, there is tained. The core of such low mass star (original mass 6

. 8M⊙) contracts under gravity; the outer envelope ex- degeneracy pressure can hold it up. As the core collapses pands again and the star becomes an even larger red under the hughe gravitational forces, protons and elec- giant. Eventually the outer layers escape into space, the trons are pushed together to form neutrons and neutrinos core shrinks, the star cools, descending downward in the (even though neutrinos don’t interact easily with matter, HR diagram, becoming a white dwarf. A white dwarf at these extremely high densities, they exert a tremen- contracts to the point at which the electron clouds start dous outward pressure), and the star begins to contract to overlap, but collapses no further because of the Pauli rapidly towards forming an enormously dense neutron exclusion principle (no two electrons can be in the same star [7]. The outer layers fall inward when the iron core quantum state). Arriving at this point is called electron collapses. When the core stops collapsing (this happens degeneracy. A white dwarf continues to lose internal en- when the neutrons start getting packed too tightly – neu- ergy by radiation, decreasing in temperature and becom- tron degeneracy), the outer layers crash into the core and ing dimmer until its light goes out. It has then become rebound, sending shock waves outward. These two effects a cold dark chunk of ash. – neutrino outburst and rebound shock wave – cause the Stars whose residual mass is greater than the Chan- entire star outside the core to be blow apart in a huge drasekhar limit are thought to follow a quite different explosion: a type II supernova! It has been suggested scenario. A star with this great mass can contract un- that the energy released in such a catastrophic explosion der gravity and heat up even further. In the range could form virtually all elements in the periodic table.2 9 56 56 T = 2.5 5 10 K, nuclei as heavy as 26Fe and 28Ni The presence of heavy elements on Earth and in our solar can be made.− × As massive red supergiants age, they pro- system suggests that our solar system formed from the duce “onion layers” of heavier and heavier elements in debris of a supernova. their interiors. However, the average binding energy per The core of a neutron star contracts to the point at nucleon begins to decrease beyond the iron group of iso- which all neutrons are as close together as they are in a topes. Thus, the formation of heavy nuclei from lighter nucleus. That is, the density of a neutron star is on the ones by fusion ends at the iron group. Further fusion order of 1014 times greater than normal solids and liquids would require energy, rather than release it. As a conse- on Earth. A neutron star that has a mass 1.5M⊙ would quence, a core of iron builds up in the centers of massive have a diameter of only about 20 km. ∼ supergiants. If the final mass of a neutron star is less than about Elements heavier than Ni are thought to form mainly 3 M⊙, its subsequent evolution is thought to be similar by neutron capture, particularly in supernova explosions. to that of a white dwarf. If the mass is greater than Large number of free neutrons, resulting from nuclear this, the star collapses under gravity, overcoming even reactions, are present inside highly evolved stars and they the neutron exclusion principle [8]. The star eventually 56 can readily combined with, say, a 26Fe nucleus to form (if collapses to the point of zero volume and infinite density, 59 59 59 three are captured) 26Fe which decays to 27Co. The 27Co creating what is known as a “singularity” [9]. As the can capture neutrons, also becoming neutron rich and density increases, the paths of light rays emitted from the decaying by β− to next higher Z element and so on. The star are bent and eventually wrapped irrevocably around highest Z elements are thought to form by such neutron the star. Any emitted photons are trapped into an orbit capture during supernova (SN) explosions when hordes by the intense gravitational field; they will never leave of neutrons are available. it. Because no light escapes after the star reaches this Yet at these extremely high temperatures, well above infinite density, it is called a black hole. 109 K, the kinetic energy of the nuclei is so high that Binary X-ray sources are places to find strong black fusion of elements heavier than iron is still possible even hole candidates [10]. A companion star is a perfect source though the reactions require energy input. But the high of infalling material for a black hole. As the matter falls energy collisions can also cause the breaking apart of iron or is pulled towards the black hole, it gains kinetic energy, and nickel nuclei into helium nuclei and eventually into heats up and is squeezed by tidal forces. The heating ion- protons and neutrons izes the atoms, and when the atoms reach a few million degrees Kelvin, they emit X-rays. The X-rays are sent 56Fe 13 4He+4 n (30) 26 → 2 off into space before the matter crosses the event horizon and (black hole boundary) and crashes into the singularity. Thus we can see this X-ray emission. Another sign of 4He 2 p +2 n . (31) 2 → the presence of a black hole is random variation of emit- These are energy-requiring reactions, but at such ex- ted X-rays. The infalling matter that emits X-rays does tremely high temperature and pressure there is plenty not fall into the black hole at a steady rate, but rather of energy available, enough even to force electrons and protons together to form neutrons in inverse β decay − e + p n + ν . (32) 2 Type Ia supernovae are different. They are believed to be binary → stars, one of which is a white dwarf that pulls mass from its com- Eventually, the iron core reaches the Chandrasekhar panion. When the total mass reaches the Chandrasekhar limit, mass. When something is this massive, not even electron the star begins to collapse and then explodes as a supernova. 7 more sporadically, which causes an observable variation In order to avoid this paradox, both Loys de Chéseaux in X-ray intensity. Additionally, if the X-ray source is in (1744) [11] and Heinrich Olbers (1826) [12] postulated the a binary system, the X-rays will be periodically cut off existence of an interstellar medium that absorbs the light as the source is eclipsed by the companion star. When from very distant stars responsible for the divergence of looking for black hole candidates, all these things are the integral in Eq. (33). However, this resolution of the taken into account. Many X-ray satellites have scanned paradox is unsatisfactory, because in an eternal universe the skies for X-ray sources that might be possible black the temperature of the interstellar medium would have hole candidates. For a hole with the mass of the sun, to rise until the medium was in thermal equilibrium with the density is about 1019 kg/m3. Such a density is about the starlight, in which case it would be emitting as much the highest that can be created through gravitational col- energy as it absorbs, and hence could not reduce the av- lapse. Observationally, the lightest black hole candidates erage radiant energy density. The stars themselves are 30 are about 6 M⊙ 6 10 kg [3]. of course opaque, and totally block out the light from ∼ × sufficiently distant sources, but if this is the resolution of the so-called “Olbers paradox” then every line of segment VIII. THE OLBERS PARADOX must terminate at the surface of a star, so the whole sky should have a temperature equal to that at the surface The XVI century finally saw what came to be a water- of a typical star. shed in the development of Cosmology. In 1543 Nicolas Copernicus published his treatise “De Revolutionibus Or- bium Celestium” (The Revolution of Celestial Spheres) where a new view of the world is presented: the helio- centric model. It is hard to underestimate the importance of this IX. THE EXPANSION OF THE UNIVERSE work: it challenged the age long views of the way the universe worked and the preponderance of the Earth and, by extension, of human beings. The realization that we, Even though at first glance the stars seem motion- our planet, and indeed our solar system (and even our less we are going to see that this impression is illusory. Galaxy) are quite common in the heavens and reproduced There is observational evidence that stars move at speeds by myriads of planetary systems, provided a sobering ranging up to a few hundred kilometers per second, so (though unsettling) view of the universe. All the reas- in a year a fast moving star might travel 1010 km. ∼ surances of the cosmology of the Middle Ages were gone, This is 103 times less than the distance to the clos- and a new view of the world, less secure and comfort- est star, so their apparent position in the sky changes able, came into being. Despite these “problems” and the very slowly. For example, the relatively fast moving many critics the model attracted, the system was soon star known as Barnard’s star is at a distance of about accepted by the best minds of the time such as Galileo. 56 1012 km; it moves across the line of sight at about × The simplest and most ancient of all astronomical ob- 89 km/s, and in consequence its apparent position shifts servations is that the sky grows dark when the Sun goes (so-called “proper motion”) in one year by an angle of down. This fact was first noted by Johannes Kepler, who, 0.0029 degrees. The apparent position in the sky of the in the XVII century, used it as evidence for a finite uni- more distant stars changes so slowly that their proper verse. In the XIX century, when the idea of an unending, motion cannot be detected with even the most patient unchanging space filled with stars like the Sun was wid- observation. spread in consequence of the Copernican revolution, the The observations that we will discuss in this class re- question of the dark night sky became a problem. To veal that the universe is in a state of violent explosion, clearly ascertain this problem, note that if absorption is in which the galaxies are rushing appart at speeds ap- neglected, the aparent luminosity of a star of absolute proaching the speed of light. Moreover, we can extrapo- luminosity L at a distance r will be l = L/4πr2. If the late this explosion backwards in time and conclude that number density of such stars is a constant n, then the all the galaxies must have been much closer at the same number of stars at distances r between r and r + dr is time in the past – so close, in fact, that neither galaxies dN = 4πnr2dr, so the total radiant energy density due nor stars nor even atoms or atomic nuclei could have had to all stars is a separate existence. ∞ L ρ = l dN = 4π n r2dr Our knoweledge of the expansion of the universe rests s 4πr2 Z Z0 entirely on the fact that astronomers are able to measure ∞ the motion of a luminous body in a direction directly = Ln dr . (33) 0 along the line of sight much more accurately than they Z can measure its motion at right angles to the line of sight. The integral diverges, leading to an infinite energy den- The technique makes use of a familiar property of any sity of starlight! sort of wave motion, known as Doppler effect. 8

A. Doppler Effect example, the galaxies of the Virgo cluster are moving away from our Galaxy at a speed of about 103 km/s. ′ When we observe a sound or light wave from a source Therefore, the wavelength λ of any spectral line from at rest, the time between the arrival wave crests at our the Virgo cluster is larger than its normal value λ by a instruments is the same as the time between crests as ratio they leave the source. However, if the source is mov- λ′ T ′ 103 km/s = 1+ 1.0033 . (38) ing away from us, the time between arrivals of successive λ T ≈ 3 105 km/s ≈ wave crests is increased over the time between their de- × partures from the source, because each crest has a little This effect was apparently first pointed out for both farther to go on its journey to us than the crest before. light and sound waves by Johann Christian Doppler in The time between crests is just the wavelength divided 1842. Doppler thought that this effect might explain the by the speed of the wave, so a wave sent out by a source different colors of stars. The light from the stars that moving away from us will appear to have a longer wave- happen to be moving away from the Earth would be length than if the source were at rest. Likewise, if the shifted toward longer wavelengths, and since red light source is moving toward us, the time between arrivals has a wavelength longer than the average wavelength of of the wave crests is decreased because each successive visible light, such stars might appear more red than aver- crest has a shorter distance to go, and the waves appear age. Similarly, light from stars that happen to be moving to have a shorter wavelength. A nice analogy was put for- toward the Earth would be shifted toward shorter wave- ward by Steven Weinberg in “The First Three Minutes.” lengths, so that the star might appear unusually blue. It He compared the situation with a travelling man that was soon pointed out by Buys-Ballot and others that the has to send a letter home regularly once a week during Doppler effect has essentially nothing to do with the color his travels: while he is travelling away from home, each of a star – it is true that the blue light from a receding successive letter will have a little farther to go than the star is shifted toward the red, but at the same time some one before, so his letters will arrive a little more than of the star’s normally invisible ultraviolet light is shifted a week apart; on the homeward leg of his journey, each into the blue part of the visible spectrum, so the over- succesive letter will have a shorter distance to travel, so all color hardly changes. As we discussed in the previous they will arrive more frequently than once a week. class, stars have different colors chiefly because they have Suppose that wave crests leave a light source at regular different surface temperatures. intervals separated by a period T . If the source is moving at velocity V c away from the observer, then during the time between≪ successive crests the source moves a B. Hubble Law distance VT . This increases the time required for the wave crest to get from the source to the observer by an The Doppler effect began to be of enormous impor- amount V T/c, where c is the speed of light. Therefore, tance to astronomy in 1968, when it was applied to the time between arrival of successive wave crests at the the study of individual spectral lines. In 1815, Joseph observer is Frauenhofer first realized that when light from the Sun is allowed to pass through a slit and then through a glass VT T ′ T + . (34) prism, the resulting spectrum of colors is crossed with ≈ c hundreds of dark lines, each one an image of the slit. The The wavelength of the light upon emission is dark lines were always found at the same colors, each corresponding to a definite wavelength of light. The same λ = cT (35) dark spectral lines were also found in the same position in the spectrum of the Moon and brighter stars. It was soon and the wavelength when the light arrives is realized that these dark lines are produced by the selec- tive absorption of light of certain definite wavelengths, λ′ = cT ′ . (36) as light passes from the hot surface of a star through its cooler outer atmosphere. Each line is due to absorption Hence the ratio of these wavelengths is of light by a specific chemical element, so it became possible to determine that the elements on the Sun, such as λ′ T ′ V = 1+ . (37) sodium, iron, magnesium, calcium, and chromium, are λ T ≈ c the same as those found on Earth. The same reasoning applies if the source is moving toward the observer, except that V is replaced with V .3 For − p(1 + V/c)/(1 − V/c), where λ is the emitted wavelength as seen in a reference frame at rest with respect to the source and λ′ is the wavelength measured in a frame moving with 3 For source and observer moving away from each other, velocity V away from the source along the line of sight; for Special Relativity’s [13] correction leads to λ′/λ = relative motion toward each other, V < 0 in this formula. 9

In 1868, Sir William Huggins was able to show that the dark lines in the spectra of some of the brighter stars are shifted slightly to the red or the blue from their normal position in the spectrum of the Sun. He correctly inter- preted this as a Doppler shift, due to the motion of the star away from or toward the Earth. For example, the wavelength of every dark line in the spectrum of the star Capella is longer than the wavelength of the corresponding dark line in the spectrum of the Sun by 0.01%, this FIG. 4: Homogeneity and the Hubble Law: A string of equally shift to the red indicates that Capella is receding from spaced galaxies Z, A, B, C, ... are shown with velocities as us at 0.01% c (i.e., 30 km/s). measured from A or B or C indicated by the lengths and di- In the late 1920’s, Hubble discovered that the spec- rections of the attached arrows. The principle of homogeneity tral lines of galaxies were shifted towards the red by an requires that the velocity of C as seen by B is equal to the amount proportional to their distances [14]. If the red- velocity of B as seen by A, adding these two velocities gives shift is due to the Doppler effect, this means that the the velocity of C as seen by A, indicated by an arrow twice as galaxies move away from each other with velocities pro- long. Proceeding in this way we can fill out the whole pattern portional to their separations. The importance of this of velocities shown in the figure. As can be seen the velocities observation is that it is just what we should predict ac- obey the Hubble law: the velocity of any galaxy, as seen by others is proportional to the distance between them. This is cording to the simplest possible picture of the flow of the only pattern of velocities consistent with the principle of matter in an expanding universe. In terms of the red- homogeneity. shift z (λ′ λ)/λ, the “linear” Hubble law can be written≡ as − pattern of velocities, whatever typical galaxy the ob- z (H /c) r , (39) ≈ 0 server happens to be riding in. It is a direct mathemat- ical consequence of this principle that the relative speed where c is the speed of light, H0 is the present value of 4 of any two galaxies must be proportional to the distance the Hubble constant and r the distance to the galaxy. between them, just as found by Hubble. To see this con- For small velocities (V c), from Eq. (37) the Doppler ≪ sider three typical galaxies A, B, and C, strung out in a redshift is z V/c. Therefore, V H0r, which is the straight line, as shown in Fig. 4. Suppose that the dis- most commonly≈ used form of Hubble≈ law. The present −1 −1 tance between A and B is the same as the distance be- day Hubble expansion rate is H0 = 100 h kms Mpc , +0.04 tween B and C. Whatever the speed of B as seen from A, where h =0.71−0.03. the Cosmological principle requires that C should have the same speed relative to B. But note that C, which is twice away from A as is B, is also moving twice as fast C. The Cosmological Principle relative to A as is B. We can add more galaxies in our chain, always with the result that the speed of recession We would expect intuitively that at any given time the of any galaxy relative to any other is proportional to the universe ought to look the same to observers in all typi- distance between them. cal galaxies, and in whatever direction they look. (Here- As often happens in science, this argument can be used after we will use the label “typical” to indicate galax- both forward and backward. Hubble, in observing a pro- ies that do not have any large peculiar motion of their portionality between the distances of galaxies and their own, but are simply carried along with the general cos- speeds of recession, was indirectly verifying the Cosmo- mic flow of galaxies.) This hypothesis is so natural (at logical Principle. Contrariwise, we can take the Cosmo- least since Copernicus) that it has been called the Cos- logical Principle for granted on a priori grounds, and mological Principle by the English astrophysicist Edward deduce the relation of proportionality between distance Arthur Milne. and velocity. In this way, through the relatively easy As applied to the galaxies themselves, the Cosmologi- measurement of Doppler shifts, we are able to judge the cal Principle requires that an observer in a typical galaxy distance of very remote objects from their velocities. should see all the other galaxies moving with the same Before proceeding any further, two qualifications have to be attached to the Cosmological Principle. First, it is obviously not true on small scales – we are in a Galaxy which belongs to a small local group of other galaxies 4 To avoid confusion, it should be kept in mind that λ denotes the (including M31 and M33), which in turn lies near the wavelength of the light if observed near the place and time of enormous cluster of galaxies in Virgo. In fact, of the emission, and thus presumably take the values measured when 33 galaxies in Messier’s catalogue, almost half are in one the same atomic transition occurs in terrestrial laboratories, while λ′ is the wavelength of the light observed after its long small part of the sky, the constellation of Virgo. The Cos- journey to us. If z > 0 then λ′ > λ and we speak of a redshift; mological Principle, if at all valid, comes into play only if z < 0 then λ′ < λ , and we speak of a blueshift. when we view the universe on a scale at least as large 10 as the distance between clusters of galaxies, or about intensity of this radiation was found initially not to vary 100 million light years. Second, in using the Cosmologi- by day or night or time of the year, nor to depend on cal Principle to derive the relation of proportionality be- the direction. It came from all directions in the universe tween galactic velocities and distances, we suppose the with equal intensity, to a precision of better than 1%. It usual rule for adding V c. This, of course, was not could only be concluded that this radiation came from a problem for Hubble in≪ 1929, as none of the galaxies the universe as a whole. he studied then had a speed anywhere near the speed of The intensity of this CMB as measured at λ =7.35 cm light. Nevertheless, it is important to stress that when corresponds to a blackbody radiation at a temperature one thinks about really large distances characteristic of of about 3 K. When radiation at other wavelengths was the universe, as a whole, one must work in a theoretical measured, the intensities were found to fall on a black- framework capable of dealing with velocities approaching body curve, corresponding to a temperature of 2.725 K. the speed of light. The CMB provides strong evidence in support of the Big Bang, and gives us information about conditions in the very early universe. In fact, in the late 1940s, George D. The Cosmic Microwave Background Gamow calculated that the Big Bang origin of the universe should have generated just such a CMB [16]. The Cosmological Principle has observational support To understand why, let us look at what a Big Bang of another sort, apart from the measurements of the might have been like. The temperature must have been Doppler shifts. After making due allowances for the dis- extremely high at the start, so high that there could not tortions due to our own Galaxy and the rich nearby clus- have been any atoms in the very early stages of the uni- ter of galaxies in the constellation of Virgo, the universe verse. Instead the universe would have consisted solely seems remarkably isotropic; that is, it looks the same in of radiation (photons) and a plasma of charged electrons all directions. Now, if the universe is isotropic around and other elementary particles. The universe would have us, it must also be isotropic about every typical galaxy. been opaque - the photons in a sense “trapped,” trav- However, any point in the universe can be carried into elling very short distances before being scattered again, another point by a series of rotations around fixed cen- primarily by electrons. Indeed, the details of the CMB ters, so if the universe is isotropic around every point, provide strong evidence that matter and radiation were it is necessary also homogeneous. In what follows we once in thermal equilibrium at very high temperature. will discuss how the observation of the cosmic microwave As the universe expanded, the energy spread out over an background (CMB) provides convincing evidence for an increasingly larger volume and the temperature dropped. isotropic universe. Only when the temperature had fallen to about 3,000 K The expansion of the universe seems to suggest that was the universe cool enough to allow the combination typical objects in the universe were once much closer to- of nuclei and electrons into atoms. (In the astrophys- gether than they are right now. This is the idea for the ical literature this is usually called “recombination,” a basis that the universe began about 13.7 billion years singularly inappropriate term, for at the time we were ago as an expansion from a state of very high density considering the nuclei and electrons had never in the pre- and temperature known affectionately as the Big Bang. vious history of the universe been combined into atoms!) The Big Bang was not an explosion, because an explo- The sudden disappearence of electrons broke the thermal sion blows pieces out into the surrounding space. Instead, contact between radiation and matter, and the radiation the Big Bang was the start of an expansion of space itself. continued thereafter to expand freely. The volume of the observable universe was very small at At the moment this happened, the energy in the radi- the start and has been expanding ever since. The ini- ation field at various wavelengths was governed by the tial tiny volume of extremely dense matter is not to be conditions of the thermal equilibrium, and was there- thought of as a concentrated mass in the midst of a much fore given by the Planck blackbody formula, Eq. (10), larger space around it. The initial tiny but dense volume for a temperature equal to that of the matter 3, 000 K. ∼ was the universe – the entire universe. There would not In particular, the typical photon wavelength would have have been anything else. When we say that the universe been about one micron, and the average distance between was once smaller than it is now, we mean that the av- photons would have been roughly equal to this typical erage separation between galaxies (or other objects) was wavelength. less. Therefore, it is the size of the universe itself that What has happened to the photons since then? Indi- has increased since the Big Bang. vidual photons would not be created or destroyed, so the In 1964, Arno Penzias and Robert Wilson were ex- average distance between photons would simply increase periecing difficulty with what they assumed to be back- in proportion to the size of the universe, i.e., in pro- ground noise, or “static,” in their radio telescope [15]. portion to the average distance between typical galaxies. Eventually, they became convinced that it was real and But we saw that the effect of the cosmological redshift that it was coming from outside the Galaxy. They made is to pull out the wavelength of any ray of light as the precise measurements at wavelength λ =7.35 cm, in the universe expands; thus the wavelength of any individual microwave region of the electromagnetic spectrum. The photon would also simply increase in proportion to the 11 size of the universe. The photons would therefore remain about one typical length apart, just as for blackbody radiation. Before proceeding we will pursue this line of argument quantitavively. From Eq. (10) we can obtain the Planck distribution that gives the energy du of a blackbody radiation per unit volume, in a narrow range of wavelengths from λ to λ + dλ, 8πhc 1 du = dλ . (40) λ5 ehc/λkT 1 − For long wavelengths, the denominator in the Planck distribution may be approximated by FIG. 5: The CMB over the entire sky, color-coded to represent differences in temperature from the average 2.725 K: hc/λkT e 1 hc/λkT , (41) the color scale ranges from +300 µK (red) to −200 µK (dark − ≃ blue), representing slightly hotter and colder spots (and also Hence, in this wavelength region, variations in density.) Results are from the WMAP satellite. ◦ 8πkT The angular resolution is 0.2 [17]. du = dλ . (42) λ4 This is the Rayleigh-Jeans formula. If this formula held Therefore, the number photon density is down to arbitrarily small wavelengths, du/dλ would become infinite for λ 0, and the total energy density in kT 3 → n = 60.42198 the blackbody radiation would be infinite. Fortunately, hc as we saw before, the Planck formula for du reaches a T 3 maximum at a wavelength λ =0.2014052 hc/kT and then = 20.28 photons cm−3 , (48) falls steeply off for decreasing wavelengths. The total en- K ergy density in the blackbody radiation is and the average photon energy is ∞ 8πhc 1 u = dλ . (43) −16 λ5 ehc/λkT 1 Eγ = u/n =3.73 10 (T/K) erg . (49) Z0 − h i × Integrals of this sort can be looked up in standard ta- Now, let’s consider what happens to blackbody radi- bles of definite integrals; the result gives the Stefan- ation in an expanding universe. Suppose the size of the Boltzmann law universe changes by a factor f, for example, if it doubles 5 4 in size, then f = 2. As predicted by the Doppler effect, 8π (kT ) − u = =7.56464 10 15 (T/K)4erg/cm3 . (44) the wavelengths will change in proportion to the size of 15(hc)3 × the universe to a new value λ′ = fλ. After the expansion, (Recall that 1 J 107 erg = 6.24 1018 eV.) the energy density du′ in the new wavelength range λ′ to We can easily≡ interpret the Planck× distribution in λ′ + dλ′ is less than the original energy density du in the terms of quanta of light or photons. Each photon has old wavelength range λ + dλ, for two different reasons: an energy E = hc/λ. Hence the number dn of photons (i) Since the volume of the universe has increased by a per unit volume in blackbody radiation in a narrow range factor of f 3, as long as no photons have been created or of wavelengths from λ to λ + dλ is destroyed, the numbers of photons per unit volume has decreased by a factor of 1/f 3. du 8π 1 dn = = dλ . (45) (ii) The energy of each photon is inversely proportional hc/λ λ4 ehc/λkT 1 − to its wavelength, and therefore is decreased by a factor Then the total number of photons per unit volume is of 1/f. It follows that the energy density is decreased by an overall factor 1/f 3 1/f =1/f 4: ∞ × n = dn 1 8πhc 1 0 ′ Z du = 4 du = 5 4 dλ hc/λkT . (50) 3 ∞ 2 f λ f e 1 kT x dx − = 8π x , (46) hc 0 e 1 If we rewrite Eq. (50) in terms of the new wavelengths Z − λ′, it becomes where x = hc/(λkT ). The integral cannot be expressed in terms of elementary functions, but it can be expressed ′ 8πhc ′ 1 du = ′ dλ ′ , (51) as an infinite series λ 5 ehcf/λ kT 1 ∞ ∞ − x2 dx 1 =2 2.4 . (47) which is exactly the same as the old formula for du in ex 1 j3 ≈ 0 j=1 terms of λ and dλ, except that T has been replaced by a Z − X 12 new temperature T ′ = T/f. Therefore, we conclude that freely expanding blackbody radiation remains described by the Planck formula, but with a temperature that drops in inverse proportion to the scale of expansion. The existence of the thermal CMB gives strong support to the idea that the universe was extremely hot in its early stages. As can be seen in Fig. 5, the background is very nearly isotropic, supporting the isotropic and homogeneous models of the universe. Of course one would expect some small inhomogeneities in the CMB that would provide “seeds” around which galaxy formation could have started. These tiny inhomogeneities were first detected by the COBE (Cosmic Background Ex- plorer) [18] and by subsequent experiments with greater detail, culminating with the WMAP (Wilkinson Mi- crowave Anisotropy Probe) results [19]. FIG. 6: Spherical region of galaxies with a larger radius than the distance between clusters of galaxies, but smaller radius than any distance characterizing the universe as a whole. X. HOMOGENEOUS AND ISOTROPIC FRIEDMANN-ROBERTSON-WALKER UNIVERSES density, 4 π r3 M = ρ . (52) In 1917 Albert Einstein presented a model of the uni- 3 verse based on his theory of General Relativity [20]. It described a geometrically symmetric (spherical) space with We can now consider the motion of a galaxy of mass m at finite volume but no boundary. In accordance with the the edge of our spherical region. According to Hubble’s Cosmological Principle, the model was homogeneous and law, its velocity will be V = Hr and the corresponding isotropic. It was also static: the volume of the space did kinetic energy not change. 1 In order to obtain a static model, Einstein had to in- K = mV 2 . (53) troduce a new repulsive force in his equations [21]. The 2 size of this cosmological term is given by the cosmological In a spherical distribution of matter, the gravitational constant Λ. Einstein presented his model before the red- force on a given spherical shell depends only on the mass shifts of the galaxies were known, and taking the universe inside the shell. The potential energy at the edge of the to be static was then reasonable. When the expansion of sphere is the universe was discovered, this argument in favor of GMm 4πmr2ρG a cosmological constant vanished. Einstein himself later U = = , (54) called it the biggest blunder of his life. Nevertheless the − r − 3 most recent observations seem to indicate that a non-zero −8 3 −1 −2 cosmological constant has to be present. where G =6.67 10 cm g s is Newton’s constant of gravitation. Hence,× the total energy is The St. Petersburg physicist Alexander Friedmann [22] studied the cosmological solutions of Einstein equations. 1 GMm E = K + U = m V 2 , (55) If Λ = 0, only evolving, expanding or contracting mod- 2 − r els of the universe are possible. The general relativistic derivation of the law of expansion for the Friedmann which has to remain constant as the universe expands. models will not be given here. It is interesting that the The value of ρ corresponding to E = 0 is called the crit- existence of three types of models and their law of ex- ical density ρc. We have, pansion can be derived from purely Newtonian consider- 1 GMm ations, with results in complete agreement with the rel- E = mH2r2 ativistic treatment. Moreover, the essential character of 2 − r 1 2 2 4π 2 the motion can be obtained from a simple energy argu- = mH r Gm r ρc ment, which we discuss next. 2 − 3 1 4π Consider a spherical region of galaxies of radius r. (For = mr2 H2 Gρ = 0 (56) 2 − 3 c the purposes of this calculation we must take r to be larger than the distance between clusters of galaxies, but smaller than any distance characterizing the universe as whence a whole, as shown in Fig. 6. We also assume Λ = 0.) The 3H2 ρ = . (57) mass of this sphere is its volume times the cosmic mass c 8πG 13

The density parameter Ω is deﬁned as Ω = ρ/ρc. Now consider two points at separation r, such that their relative velocity is V . Let R(t) be a time dependent quantity representing the scale of the universe. If R increases with time, all distances, including those between galaxies, will grow. Then

R(t) r = r0 , (58) R(t0) and FIG. 7: The two dimensional analogues of the Friedmann models. A spherical surface, a plane, and a pseudo-sphere. R˙ (t) Note that the global geometry of the universe affects the sum V =r ˙ = r0 , (59) of angles of a triangle. R(t0) where dots denote derivative with respect to t. Therefore, the Hubble constant is the same as in the relativistic context. The global geometry of Friedmann models can only be understood within V R˙ (t) the general theory of relativity. H = = . (60) r R(t) What is meant by a curve space? To answer this question, recall that our normal method of viewing the 3 From the conservation of mass it follows that ρ0R0 = world is via Euclidean plane geometry. In Euclidean ge- ρR3. Using Eq. (57) for the critical density one obtains ometry there are many axioms and theorems we take for granted. Non-Euclidean geometries which involve 8πG ρ R3 Ω= 0 0 . (61) curved space have been independently imagined by Carl 3 R3H2 Friedrich Gauss (1777-1855), Janos Bólyai (1802 - 1860), and Nikolai Ivanovich Lobachevski (1793-1856). Let us The deceleration expansion is described by the decelera- try to understand the idea of a curve space by using two tion parameter q defined as dimensional surfaces. RR¨ Consider for example, the two-dimensional surface of q = . (62) a sphere. It is clearly curved, at least to us who view it − R˙ 2 from outside – from our three dimensional world. But The deceleration parameter describes the change of ex- how do the hypothetical two-dimensional creatures de- pansion R˙ . The additional factors have been included in termine whether their two-dimensional space is flat (a order to make it dimensionless, i.e., independent of the plane) or curved? One way would be to measure the choice of units of length and time. sum of the angles of a triangle. If the surface is a plane, The expansion of the universe can be compared to the the sum of the angles is 180◦. But if the space is curved, motion of a mass launched vertically from the surface of a and a sufficiently large triangle is constructed, the sum celestial body. The form of the orbit depends on the ini- of the angles would not be 180◦. To construct a triangle tial energy. In order to compute the complete orbit, the on a curve surface, say the sphere of Fig. 7, we must use mass M of the main body and the initial velocity have to the equivalent of a straight line: that is the shortest dis- be known. In cosmology, the corresponding parameters tance between two points, which is called a geodesic. On are the mean density and the Hubble constant. a sphere, a geodesic is an arc of great circle (an arc in a The E = 0 model corresponds to the “flat” Friedmann plane passing through the center of the sphere) such as model, so-called Einstein-de Sitter model. If the density Earth’s equator and longitude lines. Consider, for exam- exceeds the critial density, the expansion of any spherical ple, the triangle whose sides are two longitude lines pass- region will turn to a contraction and it will collapse to a ing from the north pole to the equator, and the third side point. This corresponds to the closed Friedmann model. is a section of the equator. The two longitude lines form ◦ Finally, if ρ<ρc, the ever-expanding hyperbolic model 90 angles with the equator. Thus, if they make an angle is obtained. with each other at the north pole of 90◦, the sum of the These three models of the universe are called the stan- angles is 270◦. This is clearly not Euclidean space. Note, dard models. They are the simplest relativistic cosmolog- however, that if the triangle is small in comparison to ical models for Λ = 0. Models with Λ = 0 are mathemat- the radius of the sphere, the angles will add up to nearly ically more complicated, but show the6 same behaviour. 180◦, and the triangle and space will seem flat. On the The simple Newtonian treatment of the expansion saddlelike surface, the sum of the angles of a triangle is problem is possible because Newtonian mechanics is ap- less than 180◦. Such a surface is said to have negative proximately valid in small regions of the universe. How- curvature. ever, although the resulting equations are formally sim- Now, what about our universe? On a large scale what ilar, the interpretation of the quantities involved is not is the overall curvature of the universe? Doe it have 14 positive curvature, negative curvature or is it flat? By solving Einstein equations, Howard Percy Robertson [23] and Arthur Geoffrey Walker [24], showed that the three 3 hypersurfaces of constant curvature (the hyper-sphere, No Big Bang the hyper-plane, and the hyper-pseudosphere) are indeed possible geometries for a homegeneous and isotropic universe undergoing expansion. If the universe had a positive curvature, the universe 2 would be closed, or finite in volume. This would not mean that the stars and galaxies extended out to a cer- Supernovae tain boundary, beyond which there is empty space. There is no boundary or edge in such a universe. If a particle were to move in a straight line in a particular direction, it 1 SNAP would eventually return to the starting point – perhaps Target Statistical Uncertainty eons of time later. On the other hand, if the curvature CMB of the space was zero or negative, the universe would be Boomerang open. It could just go on forever. expands forever In an expanding universe, the galaxies were once much 0 y (cosmological constant) ll nearer to each other. If the rate of expansion had Maxima lapses eventua vacuum energy density recol been unchanging, the inverse of the Hubble constant, −1 closed tage = H0 , would represent the age of the universe. In Friedmann-Robertson-Walker models, however, the ex- Clusters flat pansion is gradually slowing down (i.e., q< 0), and thus open the Hubble constant gives an upper limit on the age of -1 the universe, tage 14 Gyr. Of course, if Λ = 0 this upper limit for the≈ age of the universe no longer6 holds. In an expanding universe the wavelength of radiation 0 1 2 3 is proportional to R, like all other lengths. If the wavelength at the time of emission, corresponding to the scale mass density factor R, is λ, then it will be λ0 when the scale factor has increased to R0: λ0/λ = R0/R. The redshift is FIG. 8: Shown are three independent measurements of the cosmological parameters (ΩΛ, Ωm). The high-redshift super- λ λ R novae [27], galaxy cluster abundance [28] and the CMB [29] z = 0 − = 0 1 ; (63) λ R − converge nicely near ΩΛ = 0.7 and Ωm = 0.3. The upper- left shaded region, labeled “no Big Bang,” indicates bounc- i.e., the redshift of a galaxy expresses how much the scale ing cosmologies for which the universe has a turning point in factor has changed since the light was emitted. its past [30]. The lower right shaded region corresponds to a universe which is younger than long-lived radioactive iso- −1 −1 topes [31], for any value of H0 ≥ 50 km s Mpc . Also XI. THE DARK SIDE OF THE UNIVERSE shown is the expected confidence region allowed by the future SuperNova / Acceleration Probe (SNAP) mission [32]. According to the standard Big Bang model, the universe is evolving and changing. Individual stars are being created, evolving and dying as white dwarfs, neutron galaxies are moving away from us relatively faster than stars, and black holes. At the same time, the universe as those very distant ones, meaning the expansion of the a whole is expanding. One important question is whether universe in more recent epochs has sped up. This accel- the universe will continue to expand forever. Just before eration in the expansion of the universe seems to have the year 2000, cosmologists received a surprise. As we begun roughly 5 billion years ago, 8 to 9 Gyr after the discussed in the previous class, gravity was assumed to Big Bang. What could be causing the universe to acceler- be the predominant force on a large scale in the uni- ate in its expansion? One idea is a sort of quantum field, verse, and it was thought that the expansion of the uni- so-called “quintessence” [26]; another possibility suggests verse ought to be slowing down in time because gravity an energy latent in space itself (vacuum energy) related acts as an attractive force between objects. But mea- to the cosmological constant. Whatever it is, it has been surements on type Ia supernovae unexpectedly showed given the name of “dark energy.” that very distant (high z) supernovae were dimmer than expected [25]. That is, given their great distance d as We will now re-examine the problem of the galaxy at determined from their low brightness, their speed V de- the edge of the massive sphere undergoing expansion to termined from the measured z was less than expected include the dark energy effect. The galaxy will be af- according to Hubble’s law. This suggests that nearer fected by a central force due to gravity and the cosmo- 15 logical force,

4πG r3ρm 1 mr¨ = + mΛr , (64) − 3r2 3 or 4πG rρ 1 r¨ = + Λr . (65) − 3 3 In this equations, the radius r and the density ρ are changing with time. They may be expressed in terms of the scale factor R: FIG. 9: The three possibilities for expansion in Friedmann- r = (R/R0) r0 , (66) Robertson-Walker models, plus a forth possibility that in- cludes dark energy. where R is deﬁned to be R0 when the radius r = r0;

3 ρ = (R0/R) ρ0 , (67) neutrons – by weight). WMAP data indicate that the amount of normal baryonic matter in the universe is only where the density ρ = ρ0 when R = R0. Substituting 4% of the critical density. What is the other 96%? Eqs. (66) and (67) into Eq. (65), one obtains There is a strong evidence for a significant amount of a 1 nonluminous matter in the universe referred to as dark R¨ = + ΛR , (68) matter. For example, observations of the rotation of −R2 3 galaxies suggest that they rotate as they had consider- 3 where a =4πGR0ρ0/3. If Eq. (68) is multiplied on both ably more mass than we can see. From Eq. (3) it is easily sides by R˙ , the left hand side yields seen that the speed of stars revolving around the Galactic center, 1 d(R˙ 2) R˙ R¨ = , (69) 2 dt v = GM/r , (74) and thus Eq. (68) takes the form depends on the galactic massp M and their distance from the center r. Observations show that stars farther from 2a 2 the Galactic center revolve much faster than expected d(R˙ 2)= dR + Λ R dR . (70) −R2 3 from visible matter, suggesting a great deal of invisible matter. Similarly, observations of the motions of galaxies We now define R0 = R(t0). Integrating Eq. (70) from t0 within clusters also suggest they have considerably more to t gives mass than can be seen. What might this nonluminous 1 1 1 matter in the universe be? We do not know yet. It cannot R˙ 2 R˙ 2 =2a + Λ(R2 R )2 . (71) be made of ordinary (baryonic) matter, so it must consist − 0 R − R 3 − 0 0 of some other sort of elementary particle. The evolution of the scale factor for an accelerating universe as derived The constants R˙ and a can be eliminated in favor of the 0 from Eq. (73) assuming 70% of dark energy, 26% of dark Hubble constant H0 and the density parameter Ω0 Ωm. 2 ≡ matter, and 4% of baryons, is shown in Fig. 9. Because ρc =3H0 /8πG, 2a =8πGR3 ρ /3= H2R3 ρ /ρ = H2R3Ω , (72) 0 0 0 0 0 c 0 0 0 XII. GRAVITATIONAL REDSHIFT where Ω0 = ρ0/ρc. Using Eq.(72) and R˙ 0 = H0R0 in Eq. (71), and defining the vacuum energy density ΩΛ = We have seen that the global geometry of spacetime 2 Λ/(3H0 ), one obtains may (in principle) be non-Euclidean; but how does the universe look locally (i.e., on a small scale)? According 2 2 R˙ R0 R to Einstein’s theory, spacetime is curved near massive 2 2 =Ω0 +ΩΛ +1 Ω0 ΩΛ (73) bodies. We might think of space as being like a thin H0 R0 R R0 − − rubber sheet: if a heavy weight is hung from it, it curves. as the basic differential equation governing R(t). As The weight corresponds to a huge mass that causes space shown in Fig. 8, the data from the WMAP survey and (space itself!) to curve. Thus, in Einstein’s theory we other recent experiments agree well with theories when do not speak of the “force” of gravity acting on bodies. they input dark energy as providing 70% of the energy in Instead we say that bodies and light rays move as they do the universe, and when the total energy density ρ equals because spacetime is curved. A body at rest or moving the critical density ρc. Interestingly, this ρ cannot be only slowly near a great mass would follow a geodesic towards baryonic matter (atoms are 99.9% baryons – protons and that body. 16

The extreme curvature of spacetime could be produced The coordinates r1 and t1 are related via r1 = r(t1), by a black hole. A black hole, as we previously learned, where r(t) is defined by is so dense that even light cannot escape from it. To become a black hole, a body of mass M must undergo t0 dt′ r(t) ′ dr . (79) gravitational collapse, contracting by gravitational self t R(t ) ≡ 0 attraction to within a radius called the Schwarzschild ra- Z Z dius [33], Differentiation of Eq. (79) leads to dr1 = dt1/R(t1) and so Eq. (78) can be re-written as − 2GM 2 2 rs = 2 , (75) dN =4π R (t ) n(t ,L) r dt dL . (80) c 1 1 1 | 1| where G is the gravitational constant and c the speed of The total energy density of starlight is therefore light. The Schwarzschild radius also represents the event t0 horizon of the black hole. By event horizon we mean the ρ = l dN = (t ) [R(t )/R(t )]4 dt , (81) s L 1 1 0 1 surface beyond which no signals can ever reach us, and Z Z−∞ thus inform us of events that happen. As a star collapses where is the proper luminosity density toward a black hole, the light it emits is pulled harder L and harder by gravity, but we can still see it. Once the (t ) n(t ,L) dL . (82) matter passes within the event horizon the emitted light L 1 ≡ 1 cannot escape, but is pulled back in by gravity. Thus, Z a light ray propagating in a curve spacetime would be In the hot Big Bang model there is obviously no paradox, redshifted by the gravitational field. For example, the because the integral in Eq. (81) is effectively cut off at a redshift of radiation from the surface of a star of radius lower limit t1 = 0, and the integrand vanished at t1 =0, R and mass M is roughly like R(t1). In other words, the modern explanation of the paradox is that the stars have only existed for 1 zg = 1 . (76) a finite time, so the light from very distant stars has not 1 rs/R − yet reached us. Rather than providing the world to be − finite in space, as was suggested by Kepler in 1610, the p The gravitational redshift of the radiation from the galax- Olbers paradox has shown it to be of a finite age. ies is normally insignificant [34].

XIV. LOOKBACK TIME XIII. ENERGY DENSITY OF STARLIGHT When we look out from the Earth, we look back in To see how modern cosmological models avoid the Ol- time. Any other observer in the universe would see the bers paradox, we discuss next how the apparent lumi- same view. The farther an object is from us, the earlier nosity of a star l and the number of sources per unit in time the light we see left it. No matter which direction volume become corrected due to the expansion of the we look, our view of the very early universe is blocked by universe. First, note that each photon emitted with en- the “early plasma” - we can only see as far as its surface, ergy hν1 will be red-shifted to energy hν1R(t1)/R(t0) and called “the surface of last scattering,” but not into it. photons emitted at time intervals δt1 will arrive at time Wavelengths from there are red-shifted by z 1000. ≈ intervals δt1R(t0)/R(t1), where t1 is the time the light We now discuss how the parameters of the universe leaves the source, and t0 is the time the light arrives at change when looking back toward the Big Bang. As we Earth. Second, it is easily seen that the total emitted have seen, in the matter-dominated era, the total mass power by a source at a distance r1 at the time of emis- within a co-moving sphere of radius R(t) is proportional 2 2 sion will spread over a sphere of surface area 4πR (t0)r1. to the number of nuclear particles within the sphere, and Thus, the apparent luminosity is it must remain constant (i.e., 4π ρ(t) R3(t)/3 = constant). Consequently, ρ(t) R−3(t). Nevertheless, if the 2 ∝ R(t1) 1 mass density is dominated by the mass equivalent to the l = L 2 2 . (77) R(t0) 4π R (t0) r1 energy of radiation (the radiation-dominated era), then ρ(t) is proportional to the fourth power of the temper- −1 Now, we assume that at time t1 there are n(t1,L) sources ature; and because the temperature varies like R (t), per unit volume with absolute luminosity between L ρ(t) R−4(t). and L + dL. The proper volume element of space is Let∝ us assume that at a time t a typical galaxy of 3 1 dV = R (t1) dr sin θ1 dθ1 dφ1, so the number of sources mass m is at a distance R(t) from some arbitrarily chosen between r1 and r1 +dr1 with absolute luminosity between galaxy, say the Milky Way. From Eq. (55) we see that L and L + dL is 2 1 2 4 3 2 E = mR (t) H (t) πρ(t)G . (83) dN =4π R (t1) r n(t1,L) dr1 dL . (78) 2 − 3 1 17

Now, ρ(t) increases as R(t) 0 at least as fast as R−3(t), can only be affected by events occurring close enough so and thus ρ(t)R2(t) grows at→ least as fast as R−1(t) for that a ray of light would have had time to reach us since R(t) going to zero. Therefore, in order to keep the energy the beginning of the universe. Any event that occurred E constant, the two terms in brackets must nearly cancel beyond this distance could as yet have no effect on us – it each other, so that for R(t) 0, H2(t) 8πρ(t)G/3 . is beyond the horizon. Note that as we look back towards → → The characteristic expansion time texp is the reciprocal the beginning, the distance to the horizon shrinks faster of the Hubble constant, hence than the size of the universe. As can be read in Eq. (89) the size of the universe is proportional to the one-half or 1 3 two-thirds power of the time, whereas the distance to the texp(t) = . (84) horizon is simply proportional to the time. Therefore, for ≡ H(t) s8πρ(t)G earlier and earlier times, the horizon encloses a smaller However, the Hubble constant is proportional to √ρ, and and smaller portion of the universe. This implies that therefore we conclude that H R−n/2(t). The velocity at a sufficiently early time any given “typical” particle is beyond the horizon. of a typical galaxy is then ∝ Though we can “see” only as far as the surface of last V (t)= H(t) R(t)= R˙ R1−n/2(t) . (85) scattering, in recent decades a convincing theory of the ∝ origin and evolution of the “early universe” has been de- Integration of Eq. (85) leads to the relation veloped. Most of this theory is based on recent theoretical and experimental advances in elementary particle 2 R(t1) R(t2) physics. Hence, before continuing our look back through t1 t2 = , (86) − n V (t1) − V (t2) time, we make a detour to discuss the current state of the art in High Energy Physics. or equivalently,

2 1 1 t t = . (87) XV. ELEMENTARY PARTICLES 1 − 2 n H(t ) − H(t ) 1 2 We can express H(t) in terms of ρ(t) and ﬁnd that A. The Four Forces in Nature

2 3 1 1 Since the years after World War II, particle acceler- t1 t2 = . (88) ators have been a principal means of investigating the − n 8πG ρ(t ) − ρ(t ) r " 1 2 # structure of nuclei. The accelerated particles are pro- Hence, whatever the valuep of n, the timep elapsed is projectiles that probe the interior of the nuclei they strike portional to the change in the inverse square root of the and their constituents. An important factor is that faster density. This general result can also be expressed more moving projectiles can reveal more detail about the nu- simply by saying that the time required for the density clei. The wavelength of the incoming particles is given by to drop to a value ρ from some value very much greater de Broglie’s wavelength formula λ = h/p, showing that than ρ is the greater the momentum p of the bombarding particle, the shorter the wavelength and the finer the detail that can be obtained. (Here h is Planck’s constant.) 2 3 1/2 texp radiation dominated t = = − , The accepted model for elementary particle physics to- n 8πGρ 2/3 texp matter dominated r − day views quarks and leptons as the basic constituents of where we have neglected the second term in Eq. (88) ordinary matter. To understand our present-day view because ρ(t2) ρ(t1). we need to begin with the ideas up to its formulation. In summary,≫ we have shown that the time required for While the classic discoveries of Thomson [35] (the elec- the density of the universe to drop to a value ρ from much tron) and Rutherford [36] (the proton) opened successive higher earlier values is proportional to 1/√ρ, while the doors to subatomic and nuclear physics, particle physics density ρ R−n. The time is therefore proportional to may be said to have started with the discovery of the Rn/2, or equivalently∝ positron in cosmic rays by Carl Anderson at Pasadena in 1932 [37], verifying Paul Dirac’s almost simultaneous t1/2 radiation dominated era prediction of its existence [38]. By the mid 1930s, it was R t2/n = − . (89) ∝ t2/3 matter dominated era recognized that all atoms can be considered to be made − up of neutrons, protons and electrons. The structure of This remains valid until the kinetic and potential energies matter seemed fairly simple in 1935 (with a total of six have both decreased so much that they are beginning to elementary particles: the proton, the neutron, the elec- be comparable to their sum, the total energy. tron, the positron, the neutrino and the photon), but in The universe then has a sort of horizon, which shrinks the decades that followed, hundreds of other elementary rapidly as we look back toward the beginning. No signal particles were discovered. The properties and interac- can travel faster than the speed of light, so at any time we tions of these particles, and which ones should be consid- 18 ered as basic or elementary, is the essence of elementary particle physics. TABLE I: Relative strength for two protons in a nucleus of the four forces in nature. In 1935 Hideki Yukawa predicted the existence of a new particle that would in some way mediate the strong Type Relative Strength Field Particle nuclear force [39]. In analogy to photon exchange that Strong Nuclear 1 gluons mediates the electromagnetic force, Yukawa argued that Electromagnetic 10−2 photon there ought to be a particle that mediates the strong nu- Weak Nuclear 10−6 W ± Z0 clear force – the force that holds nucleons together in the Gravitational 10−38 graviton (?) nucleus. Just as the photon is called the quantum of the electromagnetic force, the Yukawa particle would represent the quantum of the strong nuclear force. Yukawa matter. Soon after their discovery in cosmic rays, pions predicted that this new particle would have a mass in- were produced in the laboratory using a particle accel- termediate between that of the electron and the proton. erator [42]. Reactions observed included pp p p π0 Hence it was called a meson (meaning “in the middle”). and pp pnπ+. The incident proton from the→ accelera- We can make a rough estimate of the mass of the meson tor must→ have sufficient energy to produce the additional as follows. For a nucleon at rest to emit a meson would mass of the pion. require energy (to make the mass) that would have to Yukawa’s theory of pion exchange as the carrier of the come from nowhere, and so such a process would violate strong force is now out of date, and has been replaced by conservation of energy. But the uncertainty principle al- quantum chromodynamics in which the basic entities are lows non-conservation of energy of an amount ∆E if it quarks, and the carriers of the strong force are gluons. occurs only for a time ∆t given by However, the basic idea of Yukawa’s theory, i.e., that forces can be understood as the exchange of particles ∆E ∆t ℏ = h/2π . (90) ≥ remains valid. We set ∆E = mc2, the energy needed to create the mass There are four known types of force – or interaction m of the meson. Now conservation of energy is violated – in nature. The electromagnetic force is carried by the only as long as the meson exists, which is the time ∆t photon, the strong force by gluons. What about the other required for the meson to pass from one nucleon to the two: the weak nuclear force and gravity? These two are also mediated by particles. The particles that transmit other. If we assume that the meson travels at relativistic + − 0 speeds, close to the speed of light c, then ∆t would be the weak force are referred to as the W , W , and Z , at most about ∆t = d/c, where d 1.5 10−15 cm and were detected in 1983 [43]. The quantum (or carrier) is the maximum distance that can separate≈ × interacting of the gravitational force is called the graviton, and has nucleons. Thus, replacing in Eq. (90) we have not yet been observed. A comparison of the four forces is given in Table I, mc2 2.2 10−11 J = 130 MeV . (91) where they are listed according to their (approximate) ≈ × relative strengths. Note that although gravity may be The mass of the predicted meson is thus very roughly the most obvious force in daily life (because of the huge 2 130 MeV/c or about 250 times the electron mass, me mass of the Earth), on a nuclear scale it is the weakest 0.51 MeV/c2. (Note, incidentally, that the electromag-≈ of the four forces, and its effect at the particle level can netic force that has an infinite range, d = , requires a nearly always be ignored. massless mediator, which is indeed the case∞ of the photon.) Just as photons can be observed as free particles, as B. Particle’s Zoo well as acting in an exchange, so it was expected that mesons might be observed directly. Such a meson was In the decades following the discovery of the pion, a searched on the cosmic radiation that enters the Earth’s great many other subnuclear particles were discovered. atmosphere. In 1937 a new particle was discovered whose One way of arranging the particles in categories is ac- mass is 106 MeV (207 times the electron mass) [40]. This cording to their interactions, since not all particles inter- is quite close to the mass predicted, but this new parti- act by means of all four of the forces known in nature cle called the muon, did not interact strongly with mat- (though all interact via gravity). ter. Then it cannot mediate the strong nuclear force if The gauge bosons include the gluons, the photon and it does not interact by means of the strong nuclear force. the W and Z particles; these are the particles that medi- Thus, the muon, which can have either a + or charge ate the strong, electromagnetic, and weak interactions re- and seems to be nothing more than a very massive− elec- spectively. The leptons are particles that do not interact tron, is not the Yukawa particle. The particle predicted via the strong force but do interact via the weak nuclear by Yukawa was finally found in 1947 [41]. It is called force. The leptons include the electron e−, the electron − the pion (π). It comes in three charged states, +, , neutrino νe, the muon µ , the muon neutrino νµ, the + − − − or 0. The π and π have a mass of 140 MeV and the tau τ , the tau neutrino ντ and their corresponding an- 0 + + + π a mass of 135 MeV. All three interact strongly with tiparticles (e , νe, µ , νµ, τ , ντ ). The third category 19 of particles is the hadron. Hadrons are those particles and Λ0-baryon where found to behave strangely in two that interact via the strong nuclear force. Therefore they ways. (a) They were always produced in pairs; for exam- are said to be strongly interacting particles. They also ple the reaction π−p K0Λ0 occurred with high prob- interact via the other forces, but the strong force pre- ability, but the similar→ reaction π−p K0n was never dominates at short distances. The hadrons include the observed to occur. (b) These strange→ particles, as they proton, the neutron, the pion, the kaon K, the Λ, and a came to be called, were produced via the strong interac- large number of other particles. Hadrons can be divided tion (i.e., at high rate), but did not decay at a fast rate into two subgroups: baryons, which are those particles characteristic of the strong interaction (even though they that have baryon number B = +1 (or B = 1 in case of decay into strongly interacting particles). their antiparticles), and mesons, which have−B =0. To explain these observations, a new quantum num- Conservation laws are indispensable in ordering this ber, strangeness, and a new conservation law, conserva- subnuclear world. The laws of conservation of energy, tion of strangeness were introduced. By assigning the of momentum, of angular momentum, of electric charge, strangeness numbers S = +1 for the kaon and S = 1 − and baryon number are found to hold precisely in all par- for the Λ, the production of strange particles in pairs ticle interactions. Also useful are the conservation laws was explained. Antiparticles were assigned opposite for the three lepton numbers (for a given lepton l, l− strangeness from their particles. Note that for the reac- − and ν are assigned L = +1 and l+ and ν are assigned tion π p K0Λ0, both the initial and final states have l l l → Ll = 1, whereas all other particles have Ll = 0), asso- S = 0 and hence strangeness is conserved. However, for − ciated− with weak interactions including decays. Never- π p K0n, the initial state has S = 0 and the final → theless, an important result has come to the fore in our state has S = 1, so strangeness would not be conserved; young XXI century: neutrinos can occasionally change and this reaction is not observed. into one another in certain circumstances, a phenomenon To explain the decay of strange particles, it is postu- called neutrino flavor oscillation [44]. This result suggests lated that strangeness is conserved in the strong interac- that neutrinos are not massless particles and that the lep- tion, but is not conserved in the weak interaction. There- ton numbers Le, Lµ, and Lτ are not perfectly conserved. fore, strange particles are forbidden by strangeness con- The sum Le + Lµ + Lτ , however, is believed to be always servation to decay to non-strange particles of lower mass conserved. via the strong interaction, but could decay by means of The lifetime of an unstable particle depends on which the weak interaction at the observed longer lifetimes of −10 −8 force is more active in causing the decay. When a 10 to 10 s. The conservation of strangeness was the stronger force influences the decay, that decay occurs first example of a partially conserved quantity. more quickly. Decays caused by the weak force typically have lifetimes of 10−13 s or longer (W and Z are exceptions). Decays via the electromagnetic force have much C. The “Standard Model” shorter lifetimes, typically about 10−16 to 10−19 s and normally involve a photon. Hundreds of particles have All particles, except the gauge bosons, are either lep- − been found with very short lifetimes, typically 10 23 s, tons or hadrons. The principal difference between these that decay via the strong interaction. Their lifetimes are two groups is that hadrons interact via the strong inter- so short that they do not travel far enough to be detected action, whereas the leptons do not. before decaying. The existence of such short-lived parti- The six leptons are considered to be truly elementary cles is inferred from their decay products. In 1952, using particles because they do not show any internal structure, a beam of π+ with varying amounts of energy directed and have no measurable size. (Attempts to determine through a hydrogen target (protons), Enrico Fermi [45] the size of the leptons have put an upper limit of about found that the number of interactions (π+ scattered) 10−18 m.) when plotted versus the pion kinetic energy had a promi- On the other hand, there are hundreds of hadrons nent peak around 200 MeV. This led Fermi to conclude and experiments indicate they do have internal struc- + that the π and proton combined momentarily to form ture. In the early 1960s, Murray Gell-Mann [46] and a short-lived particle, before coming apart again, or at George Zweig [47] proposed that none of the hadrons are least that they resonated together for a short time. This truly elementary, but instead are made up of combina- ++ “new particle” – so called the ∆ – is referred to as a tions of three pointlike entities called, somewhat whim- resonance. Since then many other resonances have been sically, quarks. Today quarks are considered the truly found. The width of these excited states is an interesting elementary particles like leptons. The three quarks origi- application of the uncertainty principle. If the particle nally proposed were labeled, u, d, s and have the names −23 only lives 10 s, then its mass will be uncertain by an up, down, and strange; respectively. At present the the- −11 amount ∆E ℏ/∆t 10 J 100 MeV, which is ory has six quarks, just as there are six leptons – based what is observed.≈ Indeed≈ the lifetimes≈ of 10−23 s for ∼ on a presumed symmetry in nature. The other three such resonances are inferred by the reverse process, i.e., quarks are called charmed, bottom, and top (labeled c, from the measured width being 100 MeV. ∼ b, and t). The names apply also to new properties of In the early 1950s, the newly found particles K0-meson each that distinguish the new quarks from the old quarks 20

quarks is referred to as the color force. This theory of the Q, TABLE II: Quark quantum numbers: charge baryon num- strong force is called quantum chromodynamics (QCD). ber B, strangeness S, charm c, “beauty” or bottomness b, and “truth” or topness t The particles that transmit the color force are called the gluons. There are eight gluons, all massless and all have namesymbol Q B S c b t color charge. Thus gluons have replaced mesons as par- 2 1 up u 3 e 3 0 0 0 0 ticles responsible for the strong (color) force. 1 1 down d − 3 e 3 0 0 0 0 One may wonder what would happen if we try to see a 1 1 strange s − 3 e 3 −10 0 0 single quark with color by reaching deep inside a hadron. 2 1 charm c 3 e 3 0 1 0 0 Quarks are so tightly bound to other quarks that extract- 1 1 ing one would require a tremendous amount of energy, so bottom b − 3 e 3 0 0 −1 0 top t − 1 e 1 0 0 0 1 much that it would be sufficient to create more quarks. 3 3 Indeed, such experiments are done at modern particle colliders and all we get is not an isolated quark, but more hadrons (quark-antiquark pairs or triplets). This prop- and which (like strangeness) are conserved in strong, but erty of quarks, that they are always bound in groups that not in weak interactions. The properties of the quarks are colorless, is called confinement. Moreover, the color are given in Table II. Antiquarks have opposite signs of force has the interesting property that, as two quarks ap- electric charge, baryon number, strangeness, charm, bot- proach each other very closely (or equivalently have high tomness, and topness. energy), the force between them becomes small. This All hadrons are considered to be made up of com- aspect is referred to asymptotic freedom [48]. binations of quarks, and their properties are described The weak force, as we have seen, is thought to be me- by looking at their quark content. Mesons consist of a diated by the W +, W −, and Z0 particles. It acts be- + quark-antiquark pair. For example, a π meson is a ud tween the “weak charges” that each particle has. Each combination: note that for ud pair Q = (2/3+1/3)e = e, elementary particle can thus have electric charge, weak + B =1/3 1/3=0, S =0+0, as they must for the π ; charge, color charge, and gravitational mass, although +− and a K = us, with Q =1, B = 0, S = 1. Baryons on one or more of these could be zero. the other hand, consist of three quarks. For example, a The fundamental particles can be classified into spin- neutron is n = ddu, whereas an antiproton is p = uud. 1/2 fermions (6 leptons and 6 quarks), and spin-1 gauge The truly elementary particles are the leptons, the bosons (γ, W ±,Z0, and g). The leptons have 18 degrees quarks, and the gauge bosons. Matter is made up mainly of freedom: each of the 3 charged leptons has 2 pos- of fermions which obey the Pauli exclusion principle, but sible chiralities and its associated anti-particle, whereas the carriers of the forces are all bosons (which do not). the 3 neutrinos and antineutrinos have only one chirality Quarks are fermions with spin 1/2 and therefore should (neutrinos are Left-handed and antineutrinos are Right- obey the exclusion principle. Yet for three particular handed).5 The quarks have 72 degrees of freedom: each − − baryons (∆++ = uuu, ∆ = ddd, and Ω = sss), all of the 6 quarks, has the associated antiparticle, three dif- three quarks would have the same quantum numbers, and ferent color states, and 2 chiralities. The gauge bosons at least two quarks have their spin in the same direction have 27 degrees of freedom: a photon has two possible because there are only two choices, spin up (ms = +1/2) polarization states, each massive gauge boson has 3, and or spin down (ms = 1/2). This would seem to violate each of the eight independent types of gluon in QCD has − the exclusion principle! 2. Not long after the quark theory was proposed, it was One important aspect of on-going research is the at- suggested that quarks have another property called color, tempt to find a unified basis for the different forces in or “color charge” (analogous to the electric charge). The nature. A so-called gauge theory that unifies the weak distinction between the six quarks (u,d,s,c,b,t) was re- and electromagnetic interactions was put forward in the ferred to as flavor. Each of the flavors of quark can have 1960s by Sheldon Lee Glashow [49], Abdus Salam [50], three colors usually designated red, green, and blue. The and Steven Weinberg [51]. In this electroweak theory, the antiquarks are colored antired, antigreen, and antiblue. weak and electromagnetic forces are seen as two differ- Baryons are made up of three quarks, one with each color. ent manifestations of a single, more fundamental elec- Mesons consist of a quark-antiquark pair of a particular troweak interaction. The electroweak theory has had color and its anticolor. Both baryons and mesons are thus colorless or white. Because the color is different for each quark, it serves to distinguish them and allows the exclusion principle to hold. Even though quark color was 5 A phenomenon is said to be chiral if it is not identical to its originally an ad hoc idea, it soon became the central fea- mirror image. The spin of a particle may be used to define a ture of the theory determining the force binding quarks handedness for that particle. The chirality of a particle is Right- together in a hadron. handed if the direction of its spin is the same as the direction of its motion. It is Left-handed if the directions of spin and motion Each quark is assumed to carry a color charge, analo- are opposite. By convention for rotation, a standard clock, tossed gous to the electric charge, and the strong force between with its face directed forwards, has Left-handed chirality. 21 many successes, including the prediction of the W ± particles as carriers of the weak force, with masses of 81 2 GeV/c2 in excellent agreement with the measured± values of 80.482 0.091 GeV/c2 [52] (and similar accuracy for the Z0),± and the prediction of the mass +9.7 2 of the top quark mt = 174.1−7.6 GeV/c in remarkable agreement with the top mass measured at the Tevatron, 2 mt = 178.0 4.3 GeV/c [53]. The combination of the electroweak± theory plus QCD for the strong interaction is referred today as the Standard Model (SM). Another intriguing aspect of particle physics is to explain why the W and Z have large masses rather than being massless like the photon. Electroweak theory suggests an explanation by means of the Higgs boson, which interacts with the Z and the W to “slow them down.” In being forced to go slower than the speed of light, they must acquire mass. The search for the Higgs boson has been a priority in particle physics. So far, searches have FIG. 10: Symmetry around a table. excluded a Higgs lighter than 115 GeV/c2 [55]. Yet it is expected to have a mass no larger than 200 GeV/c2. We are narrowing in on it. teractions exist that would change fermions into bosons Summing up, the SM is conceptually simple and con- and vice versa [59], and that all known fermions have a tains a description of the elementary particles and forces. supersymmetric boson partner. Thus, for each quark we The SM particles are 12 spin-1/2 fermions (6 quarks and know (a fermion), there would be a squark (a boson) or 6 leptons), 4 spin-1 gauge bosons and a spin-0 Higgs bo- SUSY quark. For every lepton there would be a slep- son. Seven of the 16 particles (charm, bottom, top, tau ton. Likewise for every gauge boson (photons and gluons neutrino, W, Z, gluon) were predicted by the SM before for example), there would be a SUSY fermion. But why they were observed experimentally! There is only one haven’t all these particles been detected. The best guess particle predicted by the SM which has not yet been ob- is that SUSY particles may be heavier than their conserved [54]. ventional counterparts, perhaps too heavy to have been produced in today’s accelerators. Until a supersymmetric particle is found, and this may be possible in the now D. Beyond the Standard Model coming-on-line Large Hadron Collider (LHC), SUSY is just an elegant guess. With the success of the unified electroweak theory, at- What happens between the unification distance of tempts are being made to incorporate it and QCD for the 10−32 m and more normal (larger) distances is referred to strong (color) force into a so-called grand unified theory as symmetry breaking. As an analogy, consider a table (GUT). One type of GUT has been worked out, in which that has four identical place settings as shown in Fig. 10. leptons and quarks belong to the same family (and are The table has several kinds of symmetries. For example, able to change freely from one type to the other) and it is symmetric to rotations of 90◦ degrees; that is, the the three forces are different aspects of the underlying table will look the same if everyone moved one chair to force [56]. The unity is predicted to occur, however, only the left or to the right. It is also north-south symmetric on a scale less than about 10−32 m, corresponding to an and east-west symmetric, so that swaps across the table extremely high energy of about 1016 GeV. If two elemen- do not affect the way the table looks. It also does not tary particles (leptons or quarks) approach each other to matter whether any person picks up the spoon to the left whithin this unification scale, the apparently fundamen- of the plate or the spoon to the right. However, once the tal distinction between them would not exist at this level, first person picks up either spoon, the choice is set for and a quark could readily change to a lepton, or vice all the rest of the table as well. The symmetry has been versa. Baryon and lepton numbers would not be con- broken. The underlying symmetry is still there – the blue served. The weak, electromagnetic, and strong (color) glasses could still be chosen either way – but some choice force would blend to a force of a single strength. must be made, and at that moment, the symmetry of the The unification of the three gauge coupling constants diners is broken. of the SM is sensitive to the particle content of the the- Since unification occurs at such tiny distances and huge ory. Indeed if the minimal supersymmetric extension of energies, the theory is difficult to test experimentally. the SM [57] is used, the match of the gauge coupling But it is not completely impossible. One testable pre- constants becomes much more accurate [58]. It is com- diction is the idea that the proton might decay (via, for monly believed that this matching is unlikely to be a example p π0e+) and violate conservation of baryon coincidence. Supersymmetry (SUSY) predicts that in- number. This→ could happen if two quarks approached to 22 within 10−31 m of each other. This is very unlikely at nor- where we have ignored the factor 2/3 in our order of mag- mal temperature and energy; consequently the decay of nitude calculation. At t = 10−44 s, a kind of “phase tran- a proton can only be an unlikely process. In the simplest sition” is believed to have occured during which the grav- form of GUT, the theoretical estimate of the proton life- itational force, in effect, “condensed out” as a separate time for the decay mode p π0e+ is about 1031 yr, and force. The symmetry of the four forces was broken, but this has just come within the→ realm of testability. Proton the strong, weak, and electromagnetic forces were still decays have still not been seen, and experiments put a unified, and the universe entered the grand unified era. lower limit on the proton lifetime for the above mode to There were no distinctions between quarks and leptons; be about 1033 yr [60], somewhat greater than the theo- baryon and lepton numbers were not conserved. Very retical prediction. This may seem a disappointment, but shortly thereafter, as the universe expanded considerably on the other hand, it presents a challenge. Indeed more and the temperature had dropped to about 1027 K, there complex GUTs are not affected by this result. was another phase transition and the strong force con- Even more ambitious than GUTs are the attempts to densed out at about 10−35 s after the Big Bang. Now the also incorporate gravity, and thus unify all forces of na- universe was filled with a “soup” of leptons and quarks. ture into a single theory. Superstring theory [61], in The quarks were initially free, but soon began to “con- which elementary particles are imagined not as points, dense” into mesons and baryons. With this confinement but as one-dimensional strings (perhaps 10−35 m long), of quarks, the universe entered the hadron era. is at present the best hope for unification of all forces.6 About this time, the universe underwent an incredi- ble exponential expansion, increasing in size by a factor of 1040 or 1050 in a tiny fraction of a second, perhaps −32 XVI. THE EARLY UNIVERSE 10 s. The usefulness of this inflationary scenario is that it solved major problems with earlier Big Bang models, such as explaining why the universe is flat, as well In today’s class, we will go back to the earliest of times as the thermal equilibrium to provide the nearly uniform – as close as possible to the Big Bang – and follow the CMB [65]. Inflation is now a generally accepted aspect evolution of the universe. It may be helpful to consult of the Big Bang theory. Fig. 11 as we go along. After the very brief inflationary period, the universe We begin at a time only a minuscule fraction of a would have settled back into its more regular expansion. second after the Big Bang: it is thought that prior to −35 −44 The universe was now a “soup” of leptons and hadrons. 10 s, perhaps as early as 10 s, the four forces of na- We can think of this soup as a grand mixture of parti- ture were unified – the realm Superstring theory. This is cles and antiparticles, as well as photons – all in roughly an unimaginably short time, and predictions can be only 7 equal numbers – colliding with one another frequently speculative. The temperature would have been about and exchanging energy. 1032 K, corresponding to “particles” moving about every which way with an average kinetic energy of By the time the universe was only about a microsec- ond (10−6 s) old, it had cooled to about 1013 K, corre- − sponding to an average kinetic energy of 1 GeV, and the 1.4 10 23 J/K 1023 K K kT × 1019 GeV , (92) vast majority of hadrons disappeared. To see why, let us ≈ ≈ 1.6 10−10 J/GeV ≈ × focus on the most familiar hadrons: nucleons and their antiparticles. When the average kinetic energy of particles was somewhat higher than 1 GeV, protons, neutrons, and their antiparticles were continually being created out 6 A point worth noting at this juncture: Very recently a new framework with a diametrically opposite viewpoint was put for- of the energies of collisions involving photons and other ward. The new premise suggests that the weakness of gravity particles. But just as quickly, particle and antiparticles may be evidence from large extra spatial compactified dimen- would annihilate. Hence the process of creation and anni- sions [62]. This is possible because SM particles are confined to hilation of nucleons was in equilibrium. The numbers of a 4-dimensional world (corresponding to our apparent universe) nucleons and antinucleons were high – roughly as many as and only gravity spills into the higher dimensional spacetime. Hence, if this picture is correct, gravity is not intrinsically weak, there were electrons, positrons, or photons. But as the but of course appears weak at relatively large distances of com- universe expanded and cooled, and the average kinetic mon experience because its effects are diluted by propagation energy of particles dropped below about 1 GeV, which in the extra dimensions. The distance at which the gravita- is the minimum energy needed in a typical collision to tional and electromagnetic forces might have equal strength is unknown, but a particularly interesting possibility is that it is create nucleons and antinucleons (940 MeV each), the around 10−19 m, the distance at which electromagnetic and weak process of nucleon creation could not continue. However, forces are known to unify to form the electroweak force. This the process of annihilation could continue with antinu- would imply a fundamental Planck mass, M∗ ∼ MW ∼ 1 TeV, cleons annihilating nucleons, until there were almost no at the reach of experiment [63]. nucleons left; but not quite zero! To explain our present 7 Though for the moment this is a matter of conjecture, it is ap- pealing that Superstring theory predicts a concrete candidate world, which consists mainly of matter with very little for the quintessence field that can drive the acceleration of the antimatter in sight, we must suppose that earlier in the − universe in the present epoch [64]. universe, perhaps around 10 35 s after the Big Bang, a 23

FIG. 11: Compressed graphical representation of events as the universe expanded and cooled after the Big Bang. The Laser Interferometer Space Antenna (LISA) mission, whose launch is envisaged for 2013, is also shown. LISA will be a trio of spacecrafts orbiting the Sun trying to observe gravitational waves by using laser interferometry over astronomical distances. slight excess of quarks over antiquarks was formed. This other. would have resulted in a slight excess of nucleons over When the temperature of the universe is 3 1010 K, antinucleons, and it is these “leftover” nucleons that we × are made of today. The excess of nucleons over antin- 0.11 s have elapsed. Nothing has changed qualitatively ucleons was about one part in 109. Earlier, during the – the contents of the universe are still dominated by hadron era there should have been about as many nucle- electrons, positrons, neutrinos, antineutrinos, and pho- ons as photons. After it ended, the “leftover” nucleons tons, all in thermal equilibrium, and all high above their thus numbered only about one nucleon per 109 photons, threshold temperatures. By the time the first second has and this ratio has persisted to this day. Protons, neu- past (certainly the most eventful second in history!), the 10 trons, and all other heavier particles were thus tremen- universe has cooled to about 10 K. The average ki- dously reduced in number by about 10−6 s after the Big netic energy is 1 MeV. This is still sufficient energy to Bang. The lightest hadrons, the pions, were the last ones create electrons and positrons and balance their annihi- to go, about 10−4 s after the Big Bang. Lighter particles lations reactions, since their masses correspond to about (those whose threshold temperature are below 1011 K), 0.5 MeV. However, about this time the decreasing den- including electrons, positrons, and neutrinos, were the sity and temperature have increased the mean free path dominant forms of matter, and the universe entered the of neutrinos and antineutrinos so much that they are be- lepton era. The universe was so dense that even the ginning to behave like free particles, no longer in ther- weakly interacting (anti)neutrinos (that can travel for mal equilibrium with the electrons, positrons, or photons. years through lead bricks without being scattered) were From now on neutrinos will cease to play an active role kept in thermal equilibrium with the electrons, positrons in our discussion, except that their energy will continue and photons by rapid collision with them and with each to provide part of the source of the gravitational field of the universe. 24

After 13.8 s, the temperature 3 109 K had dropped value before annihilation, sufficiently so that e+e− could no∼ longer× be formed. An- + − nihilation (e e photons) continued. Electrons and (Tν R) after = (Tν R) before = (T R) before . (95) positrons have now→ also disappeared from the universe, | | | except for a slight excess of electrons over positrons (later We conclude therefore that after the annihilation pro- to join with nuclei to form atoms). Thus, about 14 sec- cess is over, the photon temperature is higher than the onds after the Big Bang, the universe entered the radia- neutrino temperature by a factor tion era. Its major constituents were photons and neutrinos, but neutrinos partaking only in the weak force, T (T R) = |after 1.4 . (96) rarely interacted. T (T R) ≃ ν after ν after As long as thermal equilibrium was preserved, the total | entropy remained fixed. The entropy per unit volume is Therefore, even though out of thermal equilibrium, the found to be: S N T 3, where N is the effective number ∝ T T neutrinos and antineutrinos make an important contri- of species (total number of degrees of freedom) of parti- bution to the energy density. The effective number of cles in thermal equilibrium whose threshold temperatures neutrinos and antineutrinos is 6, or 3 times the effective are below T. In order to keep the total entropy constant, number of species of photons. On the other hand, the S must be proportional to the inverse cube of the size fourth power of the neutrinos temperature is less than of the universe. That is, if R is the separation between the fourth power of the photon temperature by a factor any pair of typical particles, then SR3 N T 3R3 = −4/3 ∝ T of 3 . Thus the ratio of the energy density of neutri- constant. nos and antineutrinos to that of photons is Just before the annihilation of electrons and positrons (at about 5 109 K) the neutrinos and antineutrinos had −4/3 uν /uγ =3 3=0.7 . (97) already gone× out of thermal equilibrium with the rest of the universe, so the only abundant particles in equi- From Eq. (44) we obtain the photon energy density at librium were the electrons, positrons and photons. The photon temperature T ; hence the total energy density effective total number of particle species before annhila- after electron positron annihilation is tion were Nbefore = 6. On the other hand, after the annihilation of electrons and positrons, the only remaining −14 4 3 uν + uγ =1.7uγ 1.3 10 (T/K) erg/cm . (98) abundant particles in equilibrium were photons. Hence ≃ × the effective number of particle species was Nafter =2. It We can convert this to an equivalent mass density, follows from the conservation of entropy that 2 −35 4 3 3 3 ρ = u/c 1.22 10 (T/K) g/cm . (99) 6 (T R) before = 2 (T R) after . (93) ≃ × We have seen that the baryonic matter is about 4% That is, the heat produced by the annihilation of electrons and positrons increases the quantity TR by a factor of the total energy in the universe ρb 0.04 ρc = 0.04 10−26 kg/m3. Using the proton as≈ typical bary- × −27 (T R) after 1/3 onic matter (mp = 1.67 10 kg), this implies that | =3 1.4 . (94) × 3 (T R) before ≃ the number density of baryons is nb =0.24 nucleons/m . | From Eqs. (48) and (97) we see that the neutrino num- 9 Before the annihilation of electrons and positrons, the ber density is about 10 that of nucleons, i.e., nν = neutrino temperature T was the same as the photon 2.4 108 neutrinos/m3. By assuming that neutrinos sat- ν × temperature T . But from then on, Tν simply dropped urate the dark matter density we can set an upper bound −1 like R , so for all subsequent times, Tν R equals the on the neutrino mass

−27 3 8 2 0.26 ρc 2.6 10 kg/m −35 kg 9.315 10 eV/c 2 mν < × 3 10 × −27 6 eV/c . (100) nν ≈ 2.4 108 neutrino/m ≈ neutrino × 1.67 10 kg ∼ × ×

Meanwhile, during the next few minutes, crucial events can strike each other and be able to fuse, but now cool were taking place. Beginning about 2 or 3 minutes after enough so that newly formed nuclei would not be broken the Big Bang, nuclear fusion began to occur. The tem- apart by subsequent collisions. Deuterium, helium, and perature had dropped to about 109 K, corresponding to very tiny amounts of litium nuclei were probably made. an average kinetic energy K 100 keV, where nucleons Because the universe was cooling too quickly, larger nu- h i≈ 25 clei were not made. After only a few minutes, probably not a quarter of an hour after the Big Bang, the temperature dropped far enough that nucleosynthesis stopped, not to start again for millions of years (in stars). Thus, after the first hour or so of the universe, matter consisted mainly of bare nuclei of hydrogen (about 75%) and helium (about 25%) and electrons. Nevertheless radiation continues to dominate. The evolution of the early universe is almost complete. The next important event is presumed to have occurred 380,000 years later. The universe had expanded to about 1/1000 of its present size, and the temperature had cooled to about 3000 K. The average kinetic energy of nuclei, electrons, and photons was less than 1 eV. Since ionization energies of atoms are (eV), as the temperature dropped below this point, electronsO could orbit the bare nuclei and remain there (without being ejected by collisions), thus forming atoms. With the birth of atoms, the photons – which had been continually scattering from free electrons – now became free to spread unhindered throughout the universe, i.e., the photons became decou- FIG. 12: Simulations of structure formation with varying pled from matter. The total energy contained in radia- amounts of matter in the neutrino component, i.e., vary- tion had been decreasing (lengthening in wavelength as ing neutrino mass: (top left) massless neutrinos, (top right) the universe expanded), and the total energy contained mν = 1 eV, (bottom left) mν = 7 eV, (bottom right) in matter became dominant. The universe was said to mν = 4 eV. have become matter-dominated. As the universe continued to expand the photons cooled further, to 2.7 K today, forming the CMB we detect from everywhere in the universe. WIMPs - weakly interacting massive particles. Generic After the birth of atoms, stars and galaxies begin to WIMPs were once in thermal equilibrium, but decou- form – presumably by self-gravitation around mass con- pled while strongly non-relativistic. The most promising centrations (inhomogeneities). Stars began to form about WIMP candidate is probably the lightest SUSY particle 200 million years after the Big Bang, galaxies after almost (hereafter denoted by χ) [67]. The relic abundance of χ’s 109 yr. The universe continued to evolve until today, is determined by its thermally average annihilation cross some 13.7 billion years later [66]. section at freezeout. At high temperatures the number density of χ’s is roughly the same as the number density of photons, but as the temperature drops below the mass XVII. WIMPs Mχ, the number density of WIMPs drops exponentially. This continues until the total annihilation cross section is no longer large enough to mantain equilibrium and the At present, the only dark matter particle which is WIMP number density “freezes out.” The largest an- known to exist from experiment is the neutrino. As we nihilation cross section in the early universe is expected learned in the last class, neutrinos decoupled from ther- −2 to be roughly Mχ . This implies that very massive mal equilibrium while still relativistic, constituting a hot WIMPs have such∼ a small annihilation cross section that dark matter candidate. Their contribution to the uni- their present abundance would be too large. Indeed, to verse’s matter balance is similar to that of light, but neu- account for the observed non-baryonic matter the ther- trinos play a secondary role. The role is however identifi- mally average annihilation cross section at freezeout has able: because of their large mean free path, they prevent to be larger than about 1 pb, yielding M < 10 TeV. the smaller structures in the cold dark matter from fully χ developing. This is visible in the observed distribution Many approaches have been developed to attempt to of galaxies shown in Fig. 12. Simulations of structure detect dark matter. Such endeavors include direct detec- formation with varying amounts of matter in the neu- tion experiments which hope to observe the scattering trino component can match to a variety of observations of dark matter particles with the target material of the of today’s sky, including measurement of galaxy-galaxy detector and indirect detection experiments which are correlations and temperature fluctuations on the surface designed to search for the products of WIMP annihila- of last scattering. These analyses suggest a more restric- tion into gamma-rays, anti-matter and neutrinos. The tive limit for neutrino masses than the one derived in detection sensitivity of current experiments has been im- Eq. (100), namely mν . 1 eV. proving at a steady rate, and new dark matter hunters The simplest model for cold dark matter consists of have been proposed for the coming years [68]. 26

XVIII. MULTI-MESSENGER ASTRONOMY Below 1014 eV the flux of particles is sufficiently large that individual nuclei can be studied by detectors carried A. The Photon Window aloft in balloons or satellites. From such direct experiments we know the relative abundances and the energy Conventional astronomy expands 18 decades in photon spectra of a variety of atomic nuclei, protons, electrons wavelengths, from 104 cm radio waves to 10−14 cm γ rays and positrons as well as the intensity, energy and spa- of GeV energy. The images of the Galactic Plane shown tial distribution of X-rays and γ-rays. Measurements in Fig. 13 summarize the different wave-bands. Our of energy and isotropy showed conclusively that one ob- rudimentary understanding of the GeV γ-ray sky was vious source, the Sun, is not the main source. Only greatly advanced in 1991 with the launch of the Energetic below 100 MeV kinetic energy or so, where the solar Gamma Ray Experiment Telescope (EGRET) on board wind shields protons coming from outside the solar sys- of the Compton Gamma Ray Observatory (CGRO). The tem, does the Sun dominate the observed proton flux. science returns from EGRET observations exceed pre- Spacecraft missions far out into the solar system, well launch expectations. In particular, the number of pre- away from the confusing effects of the Earth’s atmosphere viously known GeV γ ray sources increased from 1–2 and magnetosphere, confirm that the abundances around dozen to the 271 listed in the 3rd EGRET Catalog [69]. 1 GeV are strikingly similar to those found in the ordi- However, of this multitude of sources, only about 100 nary material of the solar system. Exceptions are the have been definitively associated with known astrophys- overabundance of elements like lithium, beryllium, and ical objects. Therefore, most of the γ ray sky, as we boron, originating from the spallation of heavier nuclei currently understand it, consists of unidentified objects. in the interstellar medium. 14 One of the reasons that such a small fraction of the Above 10 eV, the flux becomes so low that only sources were identified is the size of the typical γ-ray er- ground-based experiments with large apertures and long ror box of EGRET, that was about 1◦ 1◦, an area that exposure times can hope to acquire a significant number contains several candidate sources preventing× straight- of events. Such experiments exploit the atmosphere as forward identification. This leaves intringuing puzzles a giant calorimeter. The incident cosmic radiation inter- for the next generation of GeV γ ray instruments to un- acts with the atomic nuclei of air molecules and produces cover [70]. What happens at higher energies? extensive air showers which spread out over large areas. Above a few 100 GeV the universe becomes opaque to Already in 1938, Pierre Auger concluded from the size the propagation of γ rays, because of e+e− production on of extensive air showers that the spectrum extends up to 15 the radiation fields permeating the universe (see Fig. 14). and perhaps beyond 10 eV [73]. Nowadays substantial The pairs synchrotron radiate on the extragalactic mag- progress has been made in measuring the extraordinarily low flux ( 1 event km−2 yr−1) above 1019 eV. Contin- netic field before annihilation and so the photon flux is ∼ significantly depleted. Moreover, the charged particles uously running experiments using both arrays of parti- also suffer deflections on the B~ -field comouflaging the cle detectors on the ground and fluorescence detectors exact location of the sources. In other words, the injec- which track the cascade through the atmosphere, have tion photon spectrum is significantly modified en route detected events with primary particle energies somewhat 20 to Earth. This modification becomes dramatic at around above 10 eV [74]. 106 GeV where interaction with the CMB dominates and The mechanism(s) responsible for imparting an energy the photon mean free path is smaller than the Galactic of more than one Joule to a single elementary parti- radius. cle continues to present a major enigma to high energy Therefore, to study the high energy behavior of dis- physics. It is reasonable to assume that, in order to ac- tance sources we need new messengers. Nowadays the celerate a proton to energy E in a magnetic field B, the best candidates to probe the high energy universe are cos- size R of the accelerator must encompass the gyro radius mic rays, neutrinos and gravitational waves. Of course of the particle: R > Rgyro = E/B, i.e. the accelerat- in doing multi-messenger astronomy one has to face new ing magnetic field must contain the particle’s orbit. By challenges that we discuss next. dimensional analysis, this condition yields a maximum energy E = ΓBR. The Γ-factor has been included to allow for the possibility that we may not be at rest in B. Cosmic Rays the frame of the cosmic accelerator, resulting in the observation of boosted particle energies. Opportunity for In 1912 Victor Hess carried out a series of pioneer- particle acceleration to the highest energies is limited to ing balloon flights during which he measured the levels dense regions where exceptional gravitational forces cre- of ionizing radiation as high as 5 km above the Earth’s ate relativistic particle flows. All speculations involve surface [72]. His discovery of increased radiation at high collapsed objects and we can therefore replace R by the Schwarzschild radius R GM/c2 to obtain E < ΓBM. altitude revealed that we are bombarded by ionizing par- ∼ ticles from above. These cosmic ray (CR) particles are At this point a reality check is in order. Such a di- now known to consist primarily of protons, helium, car- mensional analysis applies to the Fermilab accelerator: bon, nitrogen and other heavy ions up to iron. 10 kilogauss fields over several kilometers (covered with a 27

FIG. 13: The Milky Way as it appears at (a) radio, (b) infrared, (c) visible, (d) X-ray, and (e) γ-ray wavelenghts. Each frame is a panoramic, view covering the entire sky. The center of the Galaxy, which lies in the direction of the constellation Sagittarius is at the center of each map. (NRAO; NASA; Lund Observatory; K. Dennerl and W. Voges; NASA.)

16 ﬁciencies matching the dimensional limit. It is highly

15 3K questionable that nature can achieve this feat. Theorists can imagine acceleration in shocks with an eﬃciency of 14 IR perhaps 1 10%. − 13 Given the microgauss magnetic ﬁeld of our galaxy, no + − γγ e e structures are large or massive enough to reach the en- [E/eV] 12 10 VIS ergies of the highest energy cosmic rays. Dimensional Eew 11 analysis therefore limits their sources to extragalactic

UV objects. A common speculation is that there may be 10

energy log relatively nearby active galactic nuclei powered by a bil- γp e+e−p 9 lion solar mass black holes. With kilo-Gauss fields we 11 Galactic Center Mrk501 maximum of star formation first objects form reach 10 GeV. The jets (blazars) emitted by the cen- 8 tral black hole could reach similar energies in accelerating γe γe sub-structures boosted in our direction by a Γ-factor of 7 10−5 10−4 10−3 10−2 10−1 100 101 102 103 10, possibly higher. redshift z In contrast to the irregular shape of the isotropic electromagnetic background spectrum from, say, 108 FIG. 14: Mean interaction length for photons on the ultravi- 1020 Hz, the CR energy spectrum above 109 eV can− olet (UV), visible (VIS), infrared (IR), and microwave back- be described by a series of power laws, with the flux grounds. The electroweak scale is indicated by a dashed line. falling about 3 orders of magnitude for each decade in- The redshifts of the star formation epoch and the famous γ- crease in energy (see Fig. 15). In the decade centered at ray source Markarian 501 are also indicated [71]. 1015.5 eV (the knee) the spectrum steepens from E−2.7 to∼ E−3.0. This feature, discovered around 40 years ago, is still not consistently explained. The spectrum steepens repetition rate of 105 revolutions per second) yield 1 TeV. further to E−3.3 above 1017.7 eV (the dip) and then The argument holds because, with optimized design and flattens to E−2.7 at 10∼18.5 eV (the ankle). Within the perfect alignment of magnets, the accelerator reaches ef- statistical uncertainty∼ of the data, which is large around 28

relevant to the problem,

θ (d/Mpc) (B/nG) . (102) 0.1◦ ≃ E/1020.5 eV

Speculations on the strength for the inter-galactic magnetic field range from 10−7 to 10−9 G. For the distance to a nearby galaxy at 100 Mpc, the resolution may therefore be anywhere from sub-degree to nonex- istent. Moreover, neutrons with energy & 1018 eV have a boosted cτn sufficiently large to serve as Galactic messengers.8 The decay mean free path of a neutron is 9 c Γn τ n = 9.15 (En/10 GeV) kpc, the lifetime being boosted from its rest-frame value, τ n = 886 s, to its lab value by Γn = En/mn. It is therefore reasonable to expect that the arrival directions of the very highest energy cosmic rays may provide information on the location of their sources. At very high energies, however, the universe becomes opaque to the propagation of cosmic rays. Shortly after the discovery of the CMB, Greisen, Zatsepin, and Kuzmin (GZK) [77] pointed out that this photonic mo- lasses makes the universe opaque to protons of sufficiently high energy, i.e., protons with energies beyond the photopion production threshold,

th mπ (mp + mπ/2) EpγCMB = EγCMB FIG. 15: Compilation of measurements of the differential en- −1 −3 E ergy spectrum of cosmic rays. The dotted line shows an E 1020 γCMB eV , (103) ≈ 10−3 eV power-law for comparison. Approximate integral fluxes (per steradian) are also shown. where mp (mπ) denotes the proton (pion) mass and −3 EγCMB 10 eV is a typical CMB photon energy. Af- ter pion∼ production, the proton (or perhaps, instead, a 20 10 eV, the tail of the spectrum is consistent with a neutron) emerges with at least 50% of the incoming en- simple extrapolation at that slope to the highest ener- ergy. This implies that the nucleon energy changes by an gies, possibly with a hint of a slight accumulation around −1 e-folding after a propagation distance . (σpγ nγ y) 1019.5 eV. A very widely held interpretation of the ankle −3 ∼ 15 Mpc. Here, nγ 400 cm is the number density of is that above 1018.5 eV a new population of CRs with ≈ the CMB photons, σpγ > 0.1 mb is the photopion produc- extragalactic origin begins to dominate the more steeply tion cross section, and y is the average energy fraction (in falling Galactic population. The extragalactic compo- the laboratory system) lost by a nucleon per interaction. nent seems to be dominated by protons [75]. Therefore, if ultra-high energy cosmic rays originate at The main reason why this impressive set of data fails cosmological distances, the net effect of their interactions to reveal the origin of the particles is undoubtedly that would yield a pile-up of particles around 1019.6 eV with their directions have been scrambled by the microgauss the spectrum dropping sharply thereafter. This so-called galactic magnetic fields. However, above 1019 eV pro- GZK cutoff has been recently observed [78], suggesting ton astronomy could still be possible because the ar- that cosmic ray astronomy (if possible) would be limited rival directions of electrically charged cosmic rays are no to 1019 . E/eV . 1020 [79]. longer scrambled by the ambient magnetic field of our Cosmic accelerators are also cosmic beam dumps pro- own Galaxy. Protons point back to their sources with ducing secondary photons and neutrino beams. Particles an accuracy determined by their gyroradius in the inter- accelerated near black holes pass through intense radia- galactic magnetic field B, tion fields or dense clouds of gas leading to production of secondary photons and neutrinos that accompany the d dB θ = , (101) ≃ Rgyro E

8 Neutron astronomy from nearby extragalactic sources may also where d is the distance to the source. Scaled to units be possible [76]. 29 primary cosmic ray beam. The target material, whether Kamiokande experiment [84]) recorded a total of 19 neu- a gas or photons, is likely to be sufficiently tenuous so trino interactions over a span of 13 seconds. Two and that the primary beam and the photon beam are only a half hours later (but reported sooner) astronomers in partially attenuated. Neutrinos propagate in straight the Southern Hemisphere saw the first Supernova to be lines and they can reach the Earth from far distance visible with the unaided eye since the time of Kepler, sources, but as we will discuss in the next class, these 250 years ago, and this occurred in the Large Magel- messengers need large detectors with sensitivity to weak lanic Clouds at a distance of some 50 kiloparsecs (roughly interactions. 150,000 light years). As we have seen, supernovae of the gravitational collapse type, occur when elderly stars run out of nuclear fusion energy and can no longer re- C. Cosmic Neutrinos sist the force of gravity. The neutrinos wind up carry- ing off most of the in-fall energy, some 10% of the total For a deep, sharply focused examination of the uni- mass-energy of the inner part of star of about 1.4M⊙ masses. Approximately 3 1053 ergs are released with verse a telescope is needed which can observe a particle 58 × that is not much affected by the gas, dust, and swirling about 10 neutrinos over a few seconds. This is a stag- magnetic fields it passes on its journey. The neutrino is gering thousand times the solar energy released over its the best candidate. As we have seen, neutrinos consti- whole lifetime! The awesome visible fireworks consist of tute much of the total number of elementary particles in a mere one thousandth of the energy release in neutrinos. the universe, and these neutral, weakly-interacting par- From this spectacular beginning to neutrino astronomy ticles come to us almost without any disruption straight followed many deductions about the nature of neutrinos, from their sources, traveling at very close to the speed of such as limits on mass, charge, gravitational attraction, light. A (low energy) neutrino in flight would not notice and magnetic moment. a barrier of lead fifty light years thick. When we are able Moving up in energy, neutrinos would also be in- to see outwards in neutrino light we will no doubt receive evitably produced in many of the most luminous and en- a wondrous new view of the universe. ergetic objects in the universe. Whatever the source, the machinery which accelerates cosmic rays will inevitably Neutrinos were made in staggering numbers at the time also produce neutrinos, guaranteeing that high energy of the Big Bang. Like the CMB photons, the relic neu- neutrinos surely arrive to us from the cosmos. The burn- trinos now posess little kinetic energy due to expansion ing question for would-be neutrino astronomers is, how- of the universe. There are expected to be at least 114 ever, are there enough neutrinos to detect? neutrinos per cubic centimeter, averaged over all space. Neutrino detectors must be generally placed deep un- There could be many more on Earth because of conden- derground, or in water, in order to escape the back- sation of neutrinos, now moving slowly under the grav- grounds caused by the inescapable rain of cosmic rays itational pull of our galaxy. As of now, we only have a upon the atmosphere. These cosmic rays produce many lower limit on the total mass in this free-floating, ghostly muons which penetrate deeply into the earth, in even the gas of neutrinos, but even so it is roughly equivalent to the total mass of all the visible stars in the universe. deepest mines, but of course with ever-decreasing numbers with depth. Hence the first attempts at high energy These relic neutrinos would be wonderful to observe, neutrino astronomy have been initiated underwater and and much thought has gone into seeking ways to do under ice. The lead project, called DUMAND was can- so [80]. The problem is that the probability of neutrinos celed in 1995 on account of slow progress and budget interacting within a detector decreases with the square of difficulties, but managed to make great headway in pio- the neutrino’s energy, for low energies. And even in the neering techniques, studying backgrounds, and exploring rare case when the neutrino does react in the detector detector designs. the resulting signal is frustratingly miniscule. Nobody High energy neutrinos are detected by observing the has been able to detect these lowest-energy neutrinos as Cherenkov radiation from secondary particles produced of yet. Prospects are not good for a low-energy neutrino by neutrinos interacting inside large volumes of highly telescope, at least in the near future. transparent ice or water, instrumented with a lattice of Next best are neutrinos from the nuclear burning of photomultiplier tubes. To visualize this technique, con- stars. Here we are more fortunate, as we have the Sun sider an instrumented cubic volume of side L. Assume, close-by, producing a huge flux of neutrinos, which has for simplicity, that the neutrino direction is perpendic- been detected by several experiments [81]. A forty year ular to a side of the cube (for a realistic detector the mystery persists in the deficit of about one half of the exact geometry has to be taken into account as well as numbers of neutrinos detected compared to expectations, the arrival directions of the neutrinos). To a first ap- the so-called “Solar Neutrino Problem” [82]. This deficit proximation, a neutrino incident on a side of area L2 is now thought probably to be due to neutrino oscilla- will be detected provided it interacts within the detector tions. volume, i.e. within the instrumented distance L. That A marvelous event occurred at 07:35:41 GMT on probability is 23 February 1987, when two detectors in deep mines in the US (the IMB experiment [83]) and Japan (the P =1 exp( L/λ ) L/λ , (104) − − ν ≃ ν 30

x y x FIG. 16: Initial configuration of test particles on a circle of radius L before a gravitational wave hits them. FIG. 17: The effect of a plus-polarized gravitational wave on a ring of particles. The amplitude shown in the figure is roughly h = 0.5. Gravitational waves passing through the Earth are −1 with λν = (ρNAσν ) . Here ρ is the density of the ice many billion billion times weaker than this. or water, NA Avogadro’s number and σν the neutrino- nucleon cross section. A neutrino flux (neutrinos per 2 F y cm per second) crossing a detector with cross sectional 2 area A ( L ) facing the incident beam, will produce ≃ x = AT P (105) N F events after a time T. In practice, the quantities A, P and depend on the neutrino energy and is obtained

F N by a convolution over neutrino energy above the detector 9 threshold. The Antartic Muon And Neutrino Detector Array (AMANDA), using a 1-mile-deep Antartic ice area as a Cherenkov detector, has operated for more than 5 years. The AMANDA group has reported detection of upcom- ing atmospheric neutrinos produced in cosmic ray showers, a demonstration of feasibility [85]. IceCube, the suc- FIG. 18: The eﬀect of cross-polarized gravitational waves on cesor experiment to AMANDA, would reach the sensitiv- a ring of particles. ity close to the neutrino ﬂux anticipated to accompany the highest energy cosmic rays, dubbed the Waxman- Bahcall bound [86]. This telescope, which is currently from the signatures of dark matter to the evidence for being deployed near the Amundsen-Scott station, com- superstring theory. prises a cubic-kilometer of ultra-clear ice about a mile Summing up, it seems likely that real high energy neu- below the South Pole surface, instrumented with long trino astronomy with kilometer scale projects is only a strings of sensitive photon detectors which record light few years away. Meanwhile, underground detectors wait produced when neutrinos interact in the ice [87]. Com- patiently for the next galactic supernova. In the very panion experiments in the deep Mediterranean are mov- long run, as has been the case with every venture into ing into construction phase. new parts of the electromagnetic spectrum, one can be With the Sun and supernova neutrino observations as sure that neutrino astronomy will teach us many new and proofs of concepts, next generation neutrino experiments unexpected wonders as we open a new window upon the will also scrutinize their data for new particle physics, universe [88].

D. Gravitational Waves 9 The ”eﬀective” telescope area A is not strictly equal to the geometric cross section of the instrumented volume facing the in- Ever since Isaac Newton in the XVII century, we have coming neutrino because even neutrinos interacting outside the learned that gravity is a force that acts immediately on instrumented volume may produce a suﬃcient amount of light inside the detector to be detected. In practice, A is therefore an object. In Einstein theory of General Relativity, how- determined as a function of the incident neutrino direction by ever, gravity is not a ”force” at all, but a curvature in simulation of the full detector, including the trigger. space [20]. In other words, the presence of a very mas- 31

y y mically as the wave passes. To visualize this effect, consider a perfectly flat region of spacetime with a group of motionless test particles lying in a plane, as shown in Fig. 16. When a weak gravitational wave arrives, pass- x x ing through the particles along a line perpendicular to the ring of radius L, the test particles will oscillate in a ”cruciform” manner, as indicated in Figs. 17 and 18. The area enclosed by the test particles does not change, and there is no motion along the direction of propaga- h h + + tion. The principal axes of the ellipse become L + ∆L and L ∆L. The amplitude of the wave, which measures − FIG. 19: Two linearly independent polarizations of a grav- the fraction of stretching or squeezing, is h = ∆L/L. Of itational wave are illustrated by displaying their effect on a course the size of this effect will go down the farther the ring of free particles arrayed in a plane perpendicular to the observer is from the source. Namely, h d−1, where d direction of the wave. The figure shows the distortions in the is the source distance. Any gravitational∝ waves expected original circle that the wave produces if it carries the plus- to be seen on Earth will be quite small, h 10−20. polarization or the cross-polarization. In general relativity ∼ there are only 2 independent polarizations. The ones shown The frequency, wavelength, and speed of a gravita- here are orthogonal to each other and the polarizations are tional wave are related through λ = cν. The polariza- transverse to the direction of the wave. tion of a gravitational wave is just like polarization of a light wave, except that the polarizations of a gravitational wave are at 45◦, as opposed to 90◦. In other words, sive body does not affect probed objects directly; it warps the effect of a “cross”-polarized gravitational wave (h×) the space around it first and then the objects move in the on test particles would be basically the same as a wave ◦ curved space. Inherit from such a redefinition of gravity with plus-polarization (h+), but rotated by 45 . The dif- is the concept of gravitational waves: as massive bodies ferent polarizations are summarized in Fig. 19. move around, disturbances in the curvature of spacetime In general terms, gravitational waves are radiated by can spread outward, much like a pebble tossed into a very massive objects whose motion involves acceleration, pond will cause waves to ripple outward from the source. provided that the motion is not perfectly spherically sym- Propagating at (or near) the speed of light, these distur- metric (like a spinning, expanding or contracting sphere) bances do not travel “through” spacetime as such – the or cylindrically symmetric (like a spinning disk). For fabric of spacetime itself is oscillating! example, two objects orbiting each other in a quasi- The simplest example of a strong source of gravita- Keplerian planar orbit will radiate. The power given off tional waves is a spinning neutron star with a small by a binary system of masses M1 and M2 separated a mountain on its surface. The mountain’s mass will cause distance R is curvature of the spacetime. Its movement will “stir up” spacetime, much like a paddle stirring up water. The 32 G4 (M M )2(M + M ) waves will spread out through the universe at the speed P = 1 2 1 2 . (106) of light, never stopping or slowing down. − π c5 R5 As these waves pass a distant observer, that observer will find spacetime distorted in a very particular way: dis- For the Earth-Sun system R is very large and M1 and tances between objects will increase and decrease rhyth- M2 are relatively very small, yielding

3 −11 m 4 24 30 2 24 30 32 (6.7 10 2 ) (6 10 kg 2 10 kg) (6 10 kg+2 10 kg) P = × kg s × × × × = 313 W . (107) − π (3 108 m/s)5 (1.5 1011m)5 × ×

Thus, the total power radiated by the Earth-Sun system kinetic energy of the Earth orbiting the Sun is about in the form of gravitational waves is truly tiny compared 2.7 1033 J. As the gravitational radiation is given off, to the total electromagnetic radiation given off by the it takes× about 300 J/s away from the orbit. At this rate, Sun, which is about 3.86 1026 W. The energy of the it would take many billion times more than the current gravitational waves comes× out of the kinetic energy of age of the Universe for the Earth to fall into the Sun. the Earth’s orbit. This slow radiation from the Earth- Although the power radiated by the Earth-Sun system Sun system could, in principle, steal enough energy to is minuscule, we can point to other sources for which the drop the Earth into the Sun. Note however that the radiation should be substantial. One important example 32 is the pair of stars (one of which is a pulsar) discovered by age arms which are 2 to 4 km in length. These are at 90◦ Russell Hulse and Joe Taylor [89]. The characteristics of angles to each other, and consist of large vacuum tubes the orbit of this binary system can be deduced from the running the entire 4 kilometers. A passing gravitational Doppler shifting of radio signals given off by the pulsar. wave will then slightly stretch one arm as it shortens the Each of the stars has a mass about 1.4 M⊙. Also, their other. This is precisely the motion to which an interfer- orbit is about 75 times smaller than the distance between ometer is most sensitive. the Earth and Sun, which means the distance between Even with such long arms, a gravitational wave will the two stars is just a few times larger than the diameter only change the distance between the ends of the arms by of our own Sun. This combination of greater masses and about 10−17 m at most. This is still only a fraction of the smaller separation means that the energy given off by the width of a proton. Nonetheless, LIGO’s interferometers Hulse-Taylor binary will be far greater than the energy are now running routinely at an even better sensitivity given off by the Earth-Sun system, roughly 1022 times as level. As of November 2005, sensitivity had reached the much. primary design specification of a detectable strain of one The information about the orbit can be used to predict part in 1021 over a 100 Hz bandwidth. just how much energy (and angular momentum) should All detectors are limited at high frequencies by shot be given off in the form of gravitational waves. As the noise, which occurs because the laser beams are made up energy is carried off, the orbit will change; the stars will of photons. If there are not enough photons arriving in draw closer to each other. This effect of drawing closer a given time interval (that is, if the laser is not intense is called an inspiral, and it can be observed in the pul- enough), it will be impossible to tell whether a measure- sar’s signals. The measurements on this system were car- ment is due to real data, or just random fluctuations in ried out over several decades, and it was shown that the the number of photons. changes predicted by gravitational radiation in General All ground-based detectors are also limited at low fre- Relativity matched the observations very well, providing quencies by seismic noise, and must be very well iso- the first experimental evidence for gravitational waves. lated from seismic disturbances. Passing cars and trains, Inspirals are very important sources of gravitational falling logs, and even waves crashing on the shore hun- waves. Any time two compact objects (white dwarfs, dreds of miles away are all very significant sources of noise neutron stars, or black holes) come close to each other, in real interferometers. they send out intense gravitational waves. As the ob- Space-based interferometers, such as LISA, are also be- jects come closer and closer to each other (that is, as R ing developed. LISA’s design calls for test masses to be becomes smaller and smaller), the gravitational waves be- placed five million kilometers apart, in separate space- come more and more intense. At some point these waves crafts with lasers running between them, as shown in should become so intense that they can be directly de- Fig. 11. The mission will use the three spacecrafts ar- tected by their effect on objects on the Earth. This direct ranged in an equilateral triangle to form the arms of a detection is the goal of several large experiments around giant Michelson interferometer with arms about 5 million the world. kilometers long. As gravitational waves pass through the The great challenge of this type of detection, though, array, they slowly squeeze and stretch the space between is the extraordinarily small effect the waves would pro- the spacecrafts. Although LISA will not be affected by duce on a detector. The amplitude of any wave will seismic noise, it will be affected by other noise sources, fall off as the inverse of the distance from the source. including noise from cosmic rays and solar wind and (of Thus, even waves from extreme systems like merging bi- course) shot noise. nary black holes die out to very small amplitude by the As we have seen in our look back through time, after time they reach the Earth. For example, the amplitude the GUT era, the universe underwent a period of fan- of waves given off by the Hulse-Taylor binary as seen tastic growth. This inflationary phase concluded with a on Earth would be roughly h 10−26. However, some violent conversion of energy into hot matter and radia- gravitational waves passing the≈ Earth could have some- tion. This reheating process also resulted in a flood of what larger amplitudes, h 10−20. For an object 1 m in gravitational waves. Because the universe is transparent length, this means that its≈ ends would move by 10−20 m to the propagation of gravitational waves, LISA will be relative to each other. This distance is about a billionth able to search for these relic ripples probing energy scales of the width of a typical atom. far beyond the eV range, associated to the surface of last A simple device to detect this motion is the laser inter- scattering [91]. ferometer, with separate masses placed many hundreds of meters to several kilometers apart acting as two ends of a bar. Ground-based interferometers are now operating, XIX. QUANTUM BLACK HOLES and taking data. The most sensitive is the Laser Inter- ferometer Gravitational Wave Observatory (LIGO) [90]. As we have seen, black holes are the evolutionary end- This is actually a set of three devices: one in Livingston, points of massive stars that undergo a supernova explo- Louisiana; the other two (essentially on top of each other) sion leaving behind a fairly massive burned out stellar in Hanford, Washington. Each consists of two light stor- remnant. With no outward forces to oppose gravitational 33

forces, the remnant will collapse in on itself. The evaporation is generally regarded as being thermal 11 The density to which the matter must be squeezed in character, with a temperature inversely proportional scales as the inverse square of the mass. For example, to its mass MBH, the Sun would have to be compressed to a radius of 1 1 3 km (about four millionths its present size) to become TBH = = , (108) a black hole. For the Earth to meet the same fate, one 8πGMBH 4 π rs would need to squeeze it into a radius of 9 mm, about and an entropy S = 2 πMBH rs, where rs is the a billionth its present size. Actually, the density of a Schwarzschild radius and we have set c = 1. Note that solar mass black hole ( 1019 kg/m3) is about the high- ∼ for a solar mass black hole, the temperature is around est that can be created through gravitational collapse. 10−6 K, which is completely negligible in today’s uni- A body lighter than the Sun resists collapse because it verse. But for black holes of 1012 kg the temperature is becomes stabilized by repulsive quantum forces between about 1012 K hot enough to emit both massless particles, subatomic particles. such as γ-rays, and massive ones, such as electrons and However, stellar collapse is not the only way to form positrons. black holes. The known laws of physics allow matter den- The black hole, however, produces an effective po- sities up to the so-called Planck value 1097 kg/m3, the tential barrier in the neighborhood of the horizon that density at which the force of gravity becomes so strong backscatters part of the outgoing radiation, modifing the that quantum mechanical fluctuations can break down blackbody spectrum. The black hole absorption cross the fabric of spacetime, creating a black hole with a ra- section, σs (a.k.a. the greybody factor), depends upon dius 10−35 m and a mass of 10−8 kg. This is the the spin of the emitted particles s, their energy Q, and ∼ lightest black hole that can be produced according to the the mass of the black hole [95]. At high frequencies conventional description of gravity. It is more massive (Qrs 1) the greybody factor for each kind of particle but much smaller in size than a proton. must≫ approach the geometrical optics limit. The inte- The high densities of the early universe were a pre- grated power emission is reasonably well approximated requisite for the formation of primordial black holes but taking such a high energy limit. Thus, for illustrative did not guarantee it. For a region to stop expanding simplicity, in what follows we adopt the geometric optics and collapse to a black hole, it must have been denser approximation, where the black hole acts as a perfect than average, so the density fluctuations were also neces- absorber of a slightly larger radius, with emitting area sary. As we have seen, such fluctuations existed, at least given by [95]

on large scales, or else structures such as galaxies and 2 clusters of galaxies would never have coalesced. For pri- A = 27πrs . (109) mordial black holes to form, these fluctuations must have been stronger on smaller scales than on large ones, which Within this framework, we can conveniently write the greybody factor as a dimensionless constant normalized is possible though not inevitable. Even in the absence of fluctuations, holes might have formed spontaneously to the black hole surface area seen by the SM fields Γs = σ /A , such that Γ = 1, Γ 2/3, and Γ at various cosmological phase transitions – for example, s 4 s=0 s=1/2 ≈ s=1 ≈ when the universe ended its early period of accelerated 1/4. All in all, a black hole emits particles with initial total expansion, known as inflation, or at the nuclear density epoch, when particles such as protons condensed out of energy between (Q,Q + dQ) at a rate the soup of their constituent quarks. −1 dN˙ i σs Q The realization that black holes could be so small = Q2 exp ( 1)2s (110) dQ 8 π2 T − − prompted Stephen Hawking to consider quantum effects, BH and in 1974 his studies lead to the famous conclusion that per degree of particle freedom i. The change of variables black holes not only swallow particles but also spit them u = Q/T, brings Eq. (110) into a more familar form, out [92]. The strong gravitational fields around the black hole induce spontaneous creation of pairs near the event 27Γ T u2 N˙ = s BH du. (111) horizon. While the particle with positive energy can es- i 128 π3 eu ( 1)2s cape to infinity, the one with negative energy has to tun- Z − − nel through the horizon into the black hole where there This expression can be easily integrated using are particle states with negative energy with respect to infinity.10 As the black holes radiate, they lose mass and ∞ zn−1 dz = Γ(n) ζ(n) (112) so will eventually evaporate completely and disappear. ez 1 Z0 −

10 One can alternatively think of the emitted particles as coming 11 Indeed both the average number [92] and the probability dis- from the singularity inside the black hole, tunneling out through tribution of the number [94] of outgoing particles in each mode the event horizon to inﬁnity [93]. obey a thermal spectrum. 34

and Using ℏ = 6.58 10−25 GeVs, Eq. (121) can then be × ∞ n−1 re-written as z 1 n dz = (2 2)Γ(n) ζ(n) , (113) −99 3 ez +1 2n − τBH 2.6 10 (MBH/GeV) s Z0 ≃ × 1.6 10−26 (M /kg)3 yr . (122) yielding ≃ × BH This implies that for a solar mass black hole, the lifetime 27Γ ˙ s is unobservably long 1064 yr, but for a 1012 kg one, it is Ni = f 3 Γ(3) ζ(3) TBH , (114) 128 π 1.5 1010 yr, about the present age of the universe. ∼ × where Γ(x) (ζ(x)) is the Gamma (Riemann zeta) function Therefore, any primordial black hole of this mass would and f = 1 (f =3/4) for bosons (fermions).12 Therefore, be completing its evaporation and exploding right now. the black hole emission rate is found to be The questions raised by primordial black holes moti- vate an empirical search for them. Most of the mass of TBH − these black holes would go into gamma rays (quarks and N˙ 7.8 1020 s 1 , (115) i ≈ × GeV gluons would hadronize mostly into pions which in turn would decay to γ-rays and neutrinos), with an energy spectrum that peaks around 100 MeV. In 1976, Hawking 20 TBH −1 and Don Page realized that γ-ray background observa- N˙ i 3.8 10 s , (116) ≈ × GeV tions place strong upper limits on the number of such black holes [96]. Specifically, by looking at the observed γ-ray spectrum, they set an upper limit of 104/pc3 on the T 11 N˙ 1.9 1020 BH s−1 , (117) density of these black holes with masses near 5 10 kg. i ≈ × GeV × Even if primordial black holes never actually formed, thinking about them has led to remarkable physical in- for particles with s =0, 1/2, 1, respectively. sights because they linked three previously disparate ar- At any given time, the rate of decrease in the black eas of physics: general relativity, quantum theory, and hole mass is just the total power radiated thermodynamics [97]. 3 M˙ BH σs Q = ci , (118) dQ − 8π2 eQ/TBH ( 1)2s i XX. HOMEWORKS X − − where ci is the number of internal degrees of freedom of 1. Using the definitions of the parsec and light year, particle species i. A straightforward calculation yields show that 1 pc = 3.26 ly. 27Γ 2. About 1350 J of energy strikes the atmosphere of M˙ = c f˜ s Γ(4) ζ(4) T 2 , (119) the Earth from the Sun per second per square meter of BH − i 128 π3 BH i area at right angle to the Sun’s rays. What is (a) the X apparent brightness l of the Sun, and (b) the absolute where f˜ = 1 (f˜ = 7/8) for bosons (fermions). Assum- luminosity L of the Sun. ing that the effective high energy theory contains ap- 3. Estimate the angular width that our Galaxy would proximately the same number of modes as the SM (i.e., subtend if observed from Andromeda. Compare to the cs=1/2 = 90, and cs=1 = 27), we find angular width of the Moon from the Earth. 4. (a) In a forest there are n trees per hectare, evenly dMBH 73 4 1 spaced. The thickness of each trunk is D. What is the =8.3 10 GeV 2 . (120) dt × MBH distance of the wood not seen for the trees? (b) How is this related to the Olbers paradox? Ignoring thresholds, i.e., assuming that the mass of the 5. According to special relativity, the Doppler shift is black hole evolves according to Eq. (120) during the en- given by tire process of evaporation, we can obtain an estimate for the lifetime of the black hole, λ′ = λ (1 + V/c)/(1 V/c) , (123) − −74 −4 2 where λ is the emittedp wavelength as seen in a reference τBH =1.2 10 GeV M dMBH . (121) × BH frame at rest with respect to the source, and λ′ is the Z wavelength measured in a frame moving with velocity V away from the source along the line of sight. Show that

12 the Doppler shift in wavelength is z V/c for V c. The Gamma function is an extension of the factorial function for 6. Through some coincidence, the≈ Balmer lines≪ from non-integer and complex numbers. If s is a positive integer, then Γ(s) = (s − 1)!. The Riemann zeta function of a real variable single ionized helium in a distant star happen to overlap ∞ s s, defined by the infinite series ζ(s) = Pn=1 1/n , converges with the Balmer lines from hydrogen in the Sun. How ∀s> 1. Using these two definitions one can now verify Eq. (47). fast is that star receding from us? [Hint: the wavelengths 35 from single-electron energy level transitions are inversely and the neutrino-nucleon cross section rises with energy proportional to the square of the atomic number of the as [98] nucleus.] σ(E ) 6.04 (E /GeV)0.358 pb . (125) 7. Astronomers have recently measured the rotation of ν ≃ ν gas around what might be a supermassive black hole of about 2 billion solar masses at the center of the Galaxy. Assume that for this energy range the Earth becomes completely opaque to the propagation of neutrinos. If the radius of the Galactic center to the gas cloud is −28 2 60 light-years, what Doppler shift ∆λ/λ do you estimate [Hint: 1b = 10 m ]. they saw? 20. (a) Estimate the power radiated in gravitational 8. Find the photon density of the 2.7 K background waves by a neutron star of M⋆ =1.4M⊙ orbiting a black hole of MBH = 20M⊙, assuming the orbital radius is radiation. 2 9. What is the maximum sum-of-the-angles for a tri- R =6GMBH/c . (b) If the kinetic energy of the neutron star orbiting the black hole is about 7 1047 J, how much angle on a sphere? × 10. Estimate the age of the Universe using the Eisntein- time will it take the neutron star to fall into the black de Sitter approximation. hole? 11. To explore the distribution of charge within nuclei, 21. Hawking has calculated quantum mechanically very-high-energy electrons are used. Electrons are used that a black hole will emit particles as if it were a hot rather than protons because they do not partake in the body with temperature T proportional to its surface strong nuclear force, so only the electric charge is inves- gravity. Since the surface gravity is inversely proportional to the black hole mass M, and the emitting area is tigated. Experiments at the Stanford linear accelerator 2 were using electrons with 1.3 GeV of kinetic energy to proportional to M , the luminosity or total power emitted is proportional to AT 4, or equivalently L M −2. obtain the charge distribution of the bismuth nucleus. ∝ What is the wavelength, and hence the expected resolu- Show that as M decreases at this rate, the black hole lifetime will be proportional to M 3. Estimate the life- tion for such a beam of electrons. 12 + + time of a black hole with a temperature T = 10 K 12. In the rare decay π e +νe, what is the kinetic 12 energy of the positron? Assume→ that the π+ decays from corresponding to M = 10 kg, i.e., about the mass of a mountain. rest and that neutrinos are massless. −3 13. Show, by conserving momentum and energy, that 22. The spectrum of black holes scales like 1/E , it is impossible for an isolated electron to radiate only a just as the ultra high energy cosmic ray background is single photon. observed to do above the so-called “knee.” Using the up- 14. Are any of the following reactions/decays possible? per bound on the primordial black hole density, estimate For those forbidden, explain what conservation law is vio- an upper limit on the black hole contribution to the high + + − 0 + 0 energy cosmic ray spectrum. Can primordial black holes lated. (a) µ e +νµ, (b) π p K n (c) π p nπ , +→ + +→ + → be the sources of the observed cosmic rays? (d) p ne νe, (e) π p pe , (f) p e νe, and (g) π−→p K0Λ0. → → Solutions: 15. (a)→ Symmetry breaking occurs in the electroweak theory at about 10−18 m. Show that this corresponds to ± 1. The parsec is the distance D when the angle φ an energy that is on the order of the mass of the W . (b) is 1 second of arc. Hence φ = 4.848 10−6, and since −32 Show that the so-called unification distance of 10 m tan φ = d/D, we obtain D = 3.086 ×1016 m = 3.26 ly. in a grand unified theory is equivalent to an energy of (d =1.5 108 km is the distance to× the Sun.) about 1016 GeV. × 16. An experiment uses 3300 tons of water waiting to 2.(a) The apparent brightness is l =1.3 103 W/m2. 0 + see a proton decay of the type p π e . If the experi- 2.(b) The absolute luminosity is L =3.7× 1026 W. ment is run for 4 yr without detecting→ decay, estimate a × lower limit on the proton lifetime. 3. The angular width is the inverse tangent of the 17. What is the temperature that corresponds to diameter of our Galaxy divided by the distance to An- 1.8 TeV collisions at the Fermilab collider? To what era dromeda, in cosmological history does this correspond? 18. Estimate the minimum energy energy of a cosmic Galaxy diameter φ = arctan 2.4◦ . (126) ray proton to excite the ∆ resonance scattering off CMB Distance to Andromeda ≈ photons. 19. Evaluate the prospects for neutrino detection at For the Moon we obtain, the km3 IceCube facility. Estimate the total number of 8 11 Moon diameter ◦ ultra-high energy events (10 GeV < Eν < 10 GeV) φ = arctan 0.52 . (127) Distance to Moon ≈ expected to be detected at IceCube during its lifetime of 15 years, if the cosmic neutrino flux is [86] Therefore, the galaxy width is about 4.5 times the Moon −8 −2 −1 −2 −1 −1 ν (Eν ) 6.0 10 (Eν /GeV) GeV cm s sr , width. F ≃ × (124) 36

4.(a) Imagine a circle with radius x around the ob- 6. The wavelengths from single electron energy level server. A fraction s(x), 0 s(x) 1, is covered by transitions are inversely proportional to the square of the trees. Then we’ll move a distance≤ dx≤outward, and draw atomic number of the nucleus. Therefore, the lines from another circle. There are 2πnxdx trees growing in the singly-ionized helium are usually one fourth the wave- annulus limited by these two circles. They hide a dis- length of the corresponding hydrogen lines. Because of tance 2πxnDdx, or a fraction nDdx of the perimeter of their redshift, the lines have 4 times their usual wave- the circle. Since a fraction s(x) was already hidden, the length (i.e., λ′ =4λ) and so contribution is only [1 s(x)]nDdx. We get − 4λ = λ (1 + V/c)/(1 V/c) V =0.88 c . (135) − ⇒ s(x + dx)= s(x)+[1 s(x)] n D dx , (128) − p which gives a diﬀerential equation for s: 7. We assume that gravity causes a centripetal accel- ds(x) eration on the gas. Using Newton’s law, = [1 s(x)] n D . (129) 2 dx − m M mgasV F = G gas BH = gas , (136) This is a separable equation which can be integrated: gravity r2 r

s x we solve for the speed of the rotating gas, ds 5 = n D dx . (130) Vgas = 6.84 10 m/s. The Doppler shift as com- 1 s × Z0 − Z0 pared to the light coming from the center of the Galaxy is then ∆λ/λ V /c 2 10−3. This yields the solution ≈ gas ≈ × s(x)=1 e−nDx . (131) 8. Substituting the CMB temperature in Eq. (48) we − obtain n 400 photons cm−3. This is the probability that in a random direction we ≈ can see at most to a distance x. This function x is a 9. The limiting value for the angles in a triangle on a cumulative probability distribution. The corresponding sphere is 540◦. Imagine drawing an equilateral triangle probability density is its derivative ds/dx. The mean free near the north pole, enclosing the north pole. If that path λ is the expectation of this distribution triangle were small, the surface would be approximate flat, and each angle on the triangle would be 60◦. Then ∞ ds(x) 1 λ = x dx = . (132) imagine “stretching” each side of the triangle down dx nD Z0 towards the equator, while keeping sure that the north For example, if there are 2000 trees per hectare, and each pole stayed inside the triangle. The angle at each vertex trunk is 10 cm thick, we can see to a distance of 50 m, of the triangle would expand, with a limiting value of 180◦. The three 180◦ angles in the triangle would add on average. ◦ 4.(b) The result can be easily generalized into 3 dimen- up to 540 . sions. Assume there are n stars per unit volume, and each has a diameter D and a surface A = πD2 perpendicular 10. Setting Ω0 =1 and ΩΛ = 0 in Eq. (73) we obtain to the line of sight. Then we have a differential equation, R˙ 2 = H2R3/R , (137) s(x)=1 e−nAx , (133) 0 0 − that can be easily integrated where λ = (nA)−1. For example, if there were one R = (3H t/2)2/3 R . (138) sun per cubic parsec, the mean free path would be 0 0 1.6 104 pc. If the universe were infinite old and infinite Then, for R = R we obtain × 0 in size, the line of sight would eventually meet a stellar 2 1 surface in any direction, although we could see very far t0 = 10 Gyr . (139) 3 H ∼ indeed. 0

5. To show that the Doppler shift in wavelength is 11. The momentum of the electron is given by the ∆λ/λ V/c for v c we use the binomial expansion: ≈ ≪ relativistic identity E2 = p2c2 + m2c4. Each electron has ′ 1/2 −1/2 an energy of 1300 MeV, which is about 2500 times the λ = λ (1 + V/c) (1 V/c) 2 − mass of the electron me 0.51 MeV/c , thus we can V V 2 V V 2 ≈ λ 1+ + 1 + ignore the second term and solve for p = E/c. Therefore ≈ 2c O c2 − −2c O c2 the de Broglie wavelength is λ = h/p = hc/E, 2 2 λ[1 + V/c + (V /c )] . (134) 6.63 10−34 Js 3 108 m/s ≈ O λ = × × 1.3 109 eV 1.6 10−19 J/eV × × 0.96 10−15 m . (140) ≃ × 37

This resolution, of approximately 1 fm, is on the order 15.(a) We use Eq. (90) to estimate the mass of the of the size of a nucleus. particle based on the given distance,

hc 12. We use conservation of energy, mc2 =1.98 1011 eV 200 GeV . (146) ≈ 2πd × ≈ 2 m + c = E + + E , (141) π e ν This value is of the same order of magnitude as the mass of the W ±. and momentum. Since the pion decays from rest, the 15.(b) We take the de Broglie wavelength, λ = h/p, momentum before the decay is zero. Thus, the momen- as the unification distance. The energy is so high that tum after the decay is also zero, and so the magnitude 15 we assume E = pc, and thus EGUT =3.1 10 GeV. A of momenta of the positron and the neutrino are equal, 16 × 2 2 2 2 similar result (EGUT = 2 10 GeV) can be obtained p + = p . Therefore, p + c = p c . The positron and e ν e ν using Eq. (90). Therefore,× this energy is the amount neutrino momenta are obtained from the relativistic iden- 2 2 4 2 that could be violated in conservation of energy if the tity, yielding E + m + c = E . (Recall neutrinos are e e ν universe were the size of the unification distance. taken as massless particles).− From Eq. (141) we first find

2 2 4 2 2 16. As with radioactive decay, the number of decays E + m + c = (m + c E + ) (142) e − e π − e is proportional to the number of parent species (N), the and then solve for time interval ∆t, and the decay constant λ,

2 2 ∆N = λN∆t, (147) 1 2 me+ c − Ee+ = mπ+ c + . (143) 2 2 mπ+ where the minus sign means that N is decreasing. The rate of decay is then ∆N/∆t = λN, or equivalently, Finally the total kinetic energy of the positron is −

2 2 2 dN N −λt mπ+ c 2 me+ c = λN ln = e , (148) Ke+ = me+ c + dt − ⇒ N0 2 − 2 mπ+ 139.6 MeV 0.511 MeV where N is the number of protons present at time t =0 0.511 MeV + 0 ≃ 2 − 2 139.6 MeV and N is the number remaining after a time t. Therefore, × −λt 69.3 MeV . (144) N = N0e . The half-lifetime is deﬁned as the time for ≃ which 1/2 of the original amount of protons in a given sample have decayed,

ln 2 ln 2 13. We work in the rest frame of the isolated electron, T1/2 = λ = . (149) so that it is initially at rest. Energy conservation leads λ ⇒ T1/2 to Hence, m c2 = K + m c2 + E . (145) e e e γ ln 2 ∆N = N ∆t; (150) −T Hence, Ke = Eγ (i.e., Ke = Eγ = 0). Because the 1/2 photon has no− energy, it does not exist. Consequently, thus, for ∆N < 1 over the four-year trial, it has not been emitted. T >N ∆t ln 2 . (151) 14.(a) The reaction is forbidden because lepton number 1/2 is not conserved. To determine N, we note that each molecule of water 14.(b) The reaction is forbidden because strangeness is H2O contains 10 protons. So one mole of water (18 g, not conserved. 6 1023 molecules) contains 6 1024 protons in 18 g 14.(c) The reaction is forbidden because charge is not of× water, or about 3 1026 protons× per kg. One ton is conserved. 103 kg, so the chamber× contains 14.(d) The reaction is forbidden because energy is not conserved. N = 3.3 106 kg 3 1026 protons/kg × × 14.(e) The reaction is forbidden because lepton number 1033 protons . (152) is not conserved. ≈ 14.(f) The reaction is forbidden because baryon num- Then a very rough estimate for the lower limit on the 33 33 ber is not conserved. proton half-life is T1/2 > 10 4 yr0.7 3 10 yr. 14.(g) The reaction is possible via the strong interac- It should be stressed that this na¨ıve calculation≈ × does tion. not take into account either the required study of background events or a statistical analysis of small 38 signals [99]. volume. Assuming the Earth is completely opaque to the propagtion of neutrinos (i.e., an eﬀective aperture 17. We can approximate the temperature from the for detection at IceCube of 2π sr) the total number of kinetic energy relationship given in Eq. (92), events expected in T =15 yr is

2 K 1011 GeV T = 3 k = 2πNT T (Eν ) σ(Eν ) dEν N 8 F 2 1.8 1012 eV 1.6 10−19 J/eV Z10 GeV = × × 7 . (156) 3 1.38 10−23 J/K ≈ × 1.4 1016 K . (153) ≈ × From Sec. XVI, we see that this temperature corresponds 20. Substituting in Eq. (106) we obtain to the hadron era. P = 2.4 1047 W. At this emission rate the neutron star will× fall into the black hole in t 2.9 s. 18. The average square center-of-mass energy in a ≈ pγCMB collision can be expressed in the Lorentz invari- 21. Substituting in Eq. (122) we obtain ant form (note that the cosine of the angle between the 10 τBH 1.5 10 yr, about the present age of the particles average away) universe.∼ × 2 s = m +2 Ep Eγ , (154) p CMB 22. The upper bound on the primordial black hole den- −3 sity yields an upper limit on the black hole contribution where Eγ 10 eV is the typical CMB photon en- CMB ∼ to the high energy cosmic ray background as observed on ergy and Ep is the energy of the incoming proton. The ∆+ decays according to Earth given by

2 2 dJ/dE [1016 eV2/m2 s sr]/(E eV)3 . (157) s = mp + mπ +2mπmp . (155) The details of≃ this tedious calculation have been pub- Thus, from Eqs. (154) and (155), we obtain lished in [100]. Comparing to the observed background Eth 1020 eV. shown in Fig. 15, pγCMB ≈ 19. The total number of target nucleons present in one dJ/dE [1024.8 eV2/m2 s sr]/(E eV)3 , (158) 3 38 ≃ km of ice is: NT = ρice NA V 6 10 , where ρice 3 ∼ × 23 −∼1 1 g/cm is the density of ice, NA = 6.02 10 mol we ﬁnd that the black signal is down by about 9 orders is Avogadro’s number, and V = 1015 cm×3 is the total of magnitude.

[1] E. Hubble, The Realm of Nebulae, (Yale University (1939). Press, New Haven, 1936; reprinted by Dover Publica- [9] R. Penrose, Phys. Rev. Lett. 14, 57 (1965); S. Hawk- tions, Inc., New York, 1958). ing, Phys. Rev. Lett. 15, 689 (1965); S. Hawking, Proc. [2] A riguorus discussion of all the topics discussed in Roy. Soc. Lond. A 294, 511 (1966); S. Hawking, Proc. this course can be found in S. Weinberg, Gravtitation Roy. Soc. Lond. A 295, 490 (1966); S. Hawking, Proc. and Cosmology (Wiley, New York, 1972). For a com- Roy. Soc. Lond. A 300, 187 (1967); S. W. Hawking and prehensive book on introductory physics, which covers R. Penrose, Proc. Roy. Soc. Lond. A 314, 529 (1970). the essential prerequisite reading material to follow the [10] A. G. W. Cameron, Nature 229, 178 (1971); R. E. Wil- course, see e.g., D. C. Giancoli, Physics (6th Edition, son, Astrophys. J. 170, 529 (1971). Pearson Prentice Hall, NJ, 2005). [11] J. P. L. de Cheseaux, Traitéde la Comète (Lausanne, [3] For further details see e.g., H. Karttunen, P. Kröger, H. 1774), pp. 223 ff; reprinted in The Bowl of Night, by F. Oja, M. Poutanen, K. J. Donner, Fundamental Astron- P. Dickson (MIT Press, Cambridge, 1968) Appendix II. omy, (4th Edition, Springer-Verlag Berlin Heidelberg [12] H. W. M. Olbers, Bode’s Jahrbuch, 111 (1826); New York, 2003). reprinted by Dickson, op. cit., Appendix I. [4] E. Hertzsprung, Astron. Nachr. 196, 201 (1913); H. N. [13] A. Einstein, Annalen Phys. 17, 891 (1905) [Annalen Russell, Science 37, 651 (1913). Phys. 14, 194 (2005)]. [5] H. A. Bethe, Phys. Rev. 55, 434 (1939). [14] E. Hubble, Proc. Nat. Acad. Sci. 15, 168 (1927). [6] S. Chandrasekhar, Mon. Not. Roy. Astron. Soc. 95, 207 [15] A. A. Penzias and R. W. Wilson, Astrophys. J. 142, (1935). 419 (1965). [7] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, [16] G. Gamow, Phys. Rev. 70, 572 (1946); R. A. Alpher, 374 (1939). H. Bethe and G. Gamow, Phys. Rev. 73, 803 (1948); G. [8] J. R. Oppenheimer and H. Snyder, Phys. Rev. 56, 455 Gamow, Rev. Mod. Phys. 21, 367 (1949). 39

[17] G. Hinshaw et al. [WMAP Collaboration], 453 (1947); C. M. G. Lattes, G. P. S. Occhialini and arXiv:astro-ph/0603451. C. F. Powell, Nature 160, 486 (1947). [18] G. F. Smoot et al., Astrophys. J. 396, L1 (1992); [42] E. Gardner and C. Lattes, Science 107, 270 (1948); W. C. L. Bennett et al., Astrophys. J. 464, L1 (1996) Barkas, E. Gardner, and C Lattes, Phys. Rev. 74 1558 [arXiv:astro-ph/9601067]. (1948); J. Richardson, Phys. Rev. 74, 1720 (1948); H. [19] D. N. Spergel et al. [WMAP Collaboration], Astrophys. York, B. Moyer and R. Bjorklund, Phys. Rev. 76, 187 J. Suppl. 148, 175 (2003) [arXiv:astro-ph/0302209]; (1949); J. Steinberger, W. Panofsky and J. Steller, Phys. D. N. Spergel et al. [WMAP Collaboration], Rev. 78, 802 (1950). arXiv:astro-ph/0603449. [43] G. Arnison et al. [UA1 Collaboration], Phys. Lett. B [20] A. Einstein, Annalen Phys. 49, 769 (1916) [Annalen 122, 103 (1983). M. Banner et al [UA2 Collaboration], Phys. 14, 517 (2005)]. Phys. Lett. B 122, 476 (1983); G. Arnison et al. [UA1 [21] A. Einstein, Sitzungsber. Preuss. Akad. Wiss. Berlin Collaboration], Phys. Lett. B 126, 398 (1983); P. Bag- (Math. Phys. ) 1917, 142 (1917). naia et al. [UA2 Collaboration], Phys. Lett. B 129, 130 [22] A. Friedmann, Z. Phys. 10, 377 (1922); A. Friedmann, (1983); G. Arnison et al. [UA1 Collaboration], Phys. Z. Phys. 21, 326 (1924). Lett. B 129, 273 (1983). [23] H. P. Robertson, Astrophys. J. 82, 284 (1935); H. P. [44] See e.g., Y. Fukuda et al. [Super-Kamiokande Col- Robertson, Astrophys. J. 83, 187, 257 (1936). laboration], Phys. Rev. Lett. 81, 1562 (1998) [24] A. G. Walker, Proc. Lond. Math. Soc. (2), 42 90 (1936). [arXiv:hep-ex/9807003]; Y. Fukuda et al. [Super- [25] A. G. Riess et al. [Supernova Search Team Kamiokande Collaboration], Phys. Rev. Lett. 82, 1810 Collaboration], Astron. J. 116, 1009 (1998) (1999) [arXiv:hep-ex/9812009]; S. Fukuda et al. [Super- [arXiv:astro-ph/9805201]; S. Perlmutter et al. [Su- Kamiokande Collaboration], Phys. Lett. B 539, 179 pernova Cosmology Project Collaboration], Astrophys. (2002) [arXiv:hep-ex/0205075]; M. H. Ahn et al. [K2K J. 517, 565 (1999) [arXiv:astro-ph/9812133]. Collaboration], Phys. Rev. Lett. 90, 041801 (2003) [26] B. Ratra and P. J. E. Peebles, Phys. Rev. D [arXiv:hep-ex/0212007]; T. Araki et al. [KamLAND 37, 3406 (1988); I. Zlatev, L. M. Wang and Collaboration], Phys. Rev. Lett. 94, 081801 (2005) P. J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999) [arXiv:hep-ex/0406035]. [arXiv:astro-ph/9807002]. [45] H. L. Anderson, E. Fermi, E. A. Long and D. E. Nagle, [27] R. A. Knop et al. [Supernova Cosmology Project Phys. Rev. 85, 936 (1952). Collaboration], Astrophys. J. 598, 102 (2003) [46] M. Gell-Mann, Phys. Lett. 8, 214 (1964). [arXiv:astro-ph/0309368]. [47] G. Zweig, CERN-TH/1964-401 and CERN-TH/1964- [28] S. W. Allen, R. W. Schmidt and A. C. Fabian, 412, unpublished. Mon. Not. Roy. Astron. Soc. 334, L11 (2002) [48] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 [arXiv:astro-ph/0205007]. (1973); H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). [29] A. E. Lange et al. [Boomerang Collaboration], Phys. [49] S. L. Glashow, Nucl. Phys. 22, 579 (1961). Rev. D 63, 042001 (2001) [arXiv:astro-ph/0005004]; [50] A. Salam, in Elementary Particle Physics, N. A. Balbi et al., Astrophys. J. 545, L1 (2000) [Erratum- Svartholm, ed. (Nobel Symposium No. 8, Almqvist & ibid. 558, L145 (2001)] [arXiv:astro-ph/0005124]. Wiksell, Stockholm, 1968), p.367. [30] S. M. Carroll, W. H. Press and E. L. Turner, Ann. Rev. [51] S. Weinberg, Phys. Rev. Lett. 19, 1264 (1967). Astron. Astrophys. 30, 499 (1992). [52] B. Abbott et al. [D0 Collaboration], Phys. Rev. Lett. [31] B. S. Meyer and D. N. Schramm, Astrophys. J. 311, 84, 222 (2000) [arXiv:hep-ex/9909030]. 406 (1986). [53] V. M. Abazov et al. [D0 Collaboration], Nature 429, [32] G. Aldering et al. [SNAP Collaboration], 638 (2004) [arXiv:hep-ex/0406031]. arXiv:astro-ph/0209550. [54] For further details see e.g., F. Halzen and A. D. Martin, [33] K. Schwarzschild, Sitzungsber. Preuss. Akad. Quarks And Leptons: An Introductory Course In Mod- Wiss. Berlin (Math. Phys. ) 1916, 424 (1916) ern Particle Physics, (New York, Usa: Wiley, 1984). [arXiv:physics/9912033]. [55] See e.g., V. M. Abazov et al. [D0 Collaboration], Phys. [34] For further details see e.g., C. W. Misner, K. S. Thorne Rev. Lett. 97, 151804 (2006) [arXiv:hep-ex/0607032]; and J. A. Wheeler, Gravitation, (W. H. Freeman & Co., A. Abulencia et al. [CDF Collaboration], Phys. San Francisco, 1973). Rev. Lett. 97, 081802 (2006) [arXiv:hep-ex/0605124]; [35] J. J. Thomson, Philosophical Magazine Fifth Series 44, V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 293 (1897). 97, 161803 (2006) [arXiv:hep-ex/0607022]; S. Schael et [36] E. Rutherford, Philosoﬁcal Magazine Sixth Series 37, al. [ALEPH Collaboration], Eur. Phys. J. C 47, 547 537 (1919). (2006) [arXiv:hep-ex/0602042]. [37] C. D. Anderson, Science 76, 238 (1932); C. D. Ander- [56] H. Georgi and S. L. Glashow, Phys. Rev. Lett. 32, 438 son, Phys. Rev. 43, 491 (1933). (1974). [38] P. A. M. Dirac, Proc. Roy. Soc. Lond. A 126, 360 [57] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 (1930); P. A. M. Dirac, Proc. Roy. Soc. Lond. A 133, (1981). 60 (1931). [58] S. Dimopoulos, S. Raby and F. Wilczek, Phys. Rev. D [39] H. Yukawa, Proc. Phys. Math. Soc. Jap. 17, 48 (1935). 24, 1681 (1981). [40] J. C. Street and E. C. Stevenson, Phys. Rev. 52, 1003 [59] J. Wess and B. Zumino, Nucl. Phys. B 70, 39 (1974); (1937). J. Wess and B. Zumino, Phys. Lett. B 49, 52 (1974). [41] C. M. G. Lattes, H. Muirhead, G. P. S. Occhialini and [60] M. Shiozawa et al. [Super-Kamiokande Col- C. F. Powell, Nature 159, 694 (1947); C. M. G. Lat- laboration], Phys. Rev. Lett. 81, 3319 (1998) tes, G. P. S. Occhialini and C. F. Powell, Nature 160, [arXiv:hep-ex/9806014]. 40

[61] M. B. Green, J. H. Schwarz and E. Witten, Super- [78] R. Abbasi et al. [HiRes Collaboration], string Theory Vol. 1: Introduction, (Cambridge Uni- arXiv:astro-ph/0703099. versity Press, 1987); M. B. Green, J. H. Schwarz and [79] For further details see e.g., L. Anchordoqui, T. Paul, E. Witten, Superstring Theory Vol. 2: Loop Amplitudes, S. Reucroft and J. Swain, Int. J. Mod. Phys. A 18, 2229 Anomalies And Phenomenology, (Cambridge University (2003) [arXiv:hep-ph/0206072]. Press 1987). [80] T. J. Weiler, Phys. Rev. Lett. 49, 234 (1982). [62] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, [81] P. Anselmann et al. [GALLEX Collaboration], Phys. Phys. Lett. B 429, 263 (1998) [arXiv:hep-ph/9803315]; Lett. B 285, 376 (1992); P. Anselmann et al. I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos [GALLEX Collaboration.], Phys. Lett. B 342, 440 and G. R. Dvali, Phys. Lett. B 436, 257 (1998) (1995); Y. Fukuda et al. [Kamiokande Collabora- [arXiv:hep-ph/9804398]. A similar framework was intro- tion], Phys. Rev. Lett. 77, 1683 (1996); Y. Fukuda duced in L. Randall and R. Sundrum, Phys. Rev. Lett. et al. [Super-Kamiokande Collaboration], Phys. Rev. 83, 3370 (1999) [arXiv:hep-ph/9905221]. Lett. 81, 1158 (1998) [Erratum-ibid. 81, 4279 (1998)] [63] M. Acciarri et al. [L3 Collaboration], Phys. Lett. B 470, [arXiv:hep-ex/9805021]; Y. Fukuda et al. [Super- 281 (1999) [arXiv:hep-ex/9910056]; B. Abbott et al. Kamiokande Collaboration], Phys. Rev. Lett. 82, 2430 [D0 Collaboration], Phys. Rev. Lett. 86, 1156 (2001) (1999) [arXiv:hep-ex/9812011]; S. Fukuda et al. [Super- [arXiv:hep-ex/0008065]; C. D. Hoyle, U. Schmidt, Kamiokande Collaboration], Phys. Rev. Lett. 86, 5651 B. R. Heckel, E. G. Adelberger, J. H. Gundlach, (2001) [arXiv:hep-ex/0103032]. D. J. Kapner and H. E. Swanson, Phys. Rev. Lett. 86, [82] R. J. Davis, D. S. Harmer and K. C. Hoffman, Phys. 1418 (2001) [arXiv:hep-ph/0011014]; L. A. Anchordo- Rev. Lett. 20, 1205 (1968); J. N. Bahcall, N. A. Bahcall qui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. and G. Shaviv, Phys. Rev. Lett. 20, 1209 (1968). Rev. D 65, 124027 (2002) [arXiv:hep-ph/0112247]; [83] R. M. Bionta et al., Phys. Rev. Lett. 58, 1494 (1987). A. Abulencia et al. [CDF Collaboration], Phys. Rev. [84] K. Hirata et al. [KAMIOKANDE-II Collaboration], Lett. 97, 171802 (2006) [arXiv:hep-ex/0605101]. Phys. Rev. Lett. 58, 1490 (1987); [64] L. Anchordoqui, H. Goldberg, S. Nawata and C. Nunez, [85] E. Andres et al., Nature 410, 441 (2001). arXiv:0704.0928 [hep-ph]. [86] E. Waxman and J. N. Bahcall, Phys. Rev. D 59, 023002 [65] A. H. Guth, Phys. Rev. D 23, 347 (1981). (1999) [arXiv:hep-ph/9807282]. [66] For further details see e.g., E. W. Kolb and [87] J. Ahrens et al. [IceCube Collaboration], Nucl. Phys. M. S. Turner, The Early Universe, Front. Phys. 69, 1 Proc. Suppl. 118, 388 (2003) [arXiv:astro-ph/0209556]. (1990). [88] For further details see e.g., F. Halzen, Phys. Rept. [67] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983); J. R. El- 333, 349 (2000); L. Anchordoqui and F. Halzen, An- lis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive and nals Phys. 321, 2660 (2006) [arXiv:hep-ph/0510389]. M. Srednicki, Nucl. Phys. B 238, 453 (1984). [89] R. A. Hulse and J. H. Taylor, Astrophys. J. 195, L51 [68] For further details see e.g., G. Bertone, D. Hooper (1975). and J. Silk, Phys. Rept. 405, 279 (2005) [90] A. Abramovici et al., Science 256, 325 (1992). [arXiv:hep-ph/0404175]. [91] For further details see e.g., K. S. Thorne, Rev. Mod. [69] R. C. Hartman et al. [EGRET Collaboration], Astro- Phys. 52, 285 (1980). K. S. Thorne, Rev. Mod. Phys. phys. J. Suppl. 123 (1999) 79. 52, 299 (1980). [70] For further details see e.g., F. A. Aharonian, Very high [92] S. W. Hawking, Nature 248, 30 (1974); S. W. Hawking, energy cosmic gamma radiation: A critical window on Commun. Math. Phys. 43 (1975) 199. the extreme universe, (Singapore: World Scientific Pub- [93] J. B. Hartle and S. W. Hawking, Phys. Rev. D 13 (1976) lishing, 2004). 2188. [71] J. G. Learned and K. Mannheim, Ann. Rev. Nucl. Part. [94] L. Parker, Phys. Rev. D 12, 1519 (1975); R. M. Wald, Sci. 50, 679 (2000). Commun. Math. Phys. 45, 9 (1975); S. W. Hawking, [72] V. F. Hess, Phys. Z. 13, 1804 (1912). Phys. Rev. D 14, 2460 (1976). [73] P. Auger, R. Maze, T. Grivet-Meyer, Comptes Rendus [95] D. N. Page, Phys. Rev. D 13 (1976) 198. 206, 1721 (1938); P. Auger, P. Ehrenfest, R. Maze, J. [96] D. N. Page and S. W. Hawking, Astrophys. J. 206, 1 Daudin, Robley, and A. Fréon, Rev. Mod. Phys. 11, 288 (1976). (1939). [97] S. W. Hawking, Phys. Rev. D 13, 191 (1976). [74] D. J. Bird et al., Astrophys. J. 441, 144 (1995). [98] L. A. Anchordoqui, A. M. Cooper-Sarkar, D. Hooper [75] D. J. Bird et al. [HIRES Collaboration], Phys. Rev. Lett. and S. Sarkar, Phys. Rev. D 74, 043008 (2006) 71, 3401 (1993). [arXiv:hep-ph/0605086]. [76] L. A. Anchordoqui, H. Goldberg and T. J. Weiler, Phys. [99] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, Rev. Lett. 87, 081101 (2001) [arXiv:astro-ph/0103043]. 3873 (1998) [arXiv:physics/9711021]. [77] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zat- [100] J. Oliensis and C. T. Hill, Phys. Lett. B 143, 92 (1984). sepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma Zh. Eksp. Teor. Fiz. 4, 114 (1966)].