Peebles 2011 Seeing Cosmology Grow
Peebles 2011 Seeing Cosmology Grow
Peebles 2011 Seeing Cosmology Grow
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org ARI 19 July 2012 14:56
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Jim Peebles ( far left), Alison Peebles (center left), Eunice Wilkinson (center right)
and Dave Wilkinson ( far right) on the occasion of Dave’s award of the National Academy
of Sciences Watson Medal.
AA50CH01-Peebles ARI 16 July 2012 11:15
ANNUAL
REVIEWS Further Seeing Cosmology Grow
Click here for quick links to
Annual Reviews content online,
including:
P.J.E. Peebles
• Other articles in this volume Joseph Henry Laboratories, Princeton University, Princeton, New Jersey 08544;
• Top cited articles email: pjep@princeton.edu
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
1
AA50CH01-Peebles ARI 16 July 2012 11:15
1. INTRODUCTION
These recollections1 of how I came to study the large-scale structure of the observable Universe—
the science of cosmology—may afford some understanding of what formed the approach I have
taken to this subject, the lessons I draw from what happened, and my thinking about where the
science may be headed.
I remember admiring the starry skies over the plains of Southern Manitoba, and aurora, but
took no particular interest in the names of planets or stars or constellations. I still admire the sky,
and my career has fed on what astronomers are doing, but I still cannot get interested in where
the planets are.
I was fascinated by how things around me work, particularly things I was allowed to take apart.
I remember reading, in an older sister’s schoolbook, an account of compound pulleys. I thought
that was neat and still do: To me, physics is compound pulleys, all the way down. The faculty at
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
the University of Manitoba (U of M) showed me that I like doing physics and taught me a lot
about how to do it. Fellow students added to my education in physics, instructed me on rudiments
of social behavior, introduced me to Alison, who was studying microbiology at the U of M, and
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
1
More recollections of how cosmology grew are in Finding the Big Bang (FTBB; Peebles, Page & Partridge 2009). I have
taken a single story line in this account, and hope friends and colleagues who would have been mentioned in a more complete
recollection of my life in science understand. I apologize also to colleagues whose contributions to this story are not mentioned
through failed judgment or memory.
2 Peebles
AA50CH01-Peebles
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org ARI 16 July 2012 11:15
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Figure 1
Bob and Annie Dicke with their first child, Nancy, in Cambridge, Massachusetts, in 1945. This is about the
time Bob was moving from the MIT Radiation Laboratory to the Department of Physics at Princeton
University.
a globular star cluster (Dicke 1939). As part of war research at the Radiation Laboratory at MIT
Bob invented a microwave radiometer (Dicke 1946); Dicke & Beringer (1946) used it to measure
microwave radiation from the Sun and Moon, and Dicke et al. (1946) used it to show that “there is
very little (<20◦ K) radiation from cosmic matter at the radiometer wave-lengths” of 1 to 1.5 cm
(Dicke et al. 1946, p. 340). Figure 1 shows the Dickes at about the time of this very interesting
measurement. Figure 2 shows Bob 35 years later with two colleagues who also had large effects
on my career.
When Dicke joined the Princeton faculty he decided he ought to do laboratory physics. His
remarkable contributions include (Dicke 1953) collisional line narrowing and (Dicke 1954) su-
perradiance that describes many particles radiating as a pure quantum state. Appreciation of the
latter took time. His paper on superradiance received 15 ADS citations in the first decade after
publication and 1042 citations in the past decade.
Dicke (1957, p. 52) explained his change of direction in the proceedings of a conference at
Chapel Hill, North Carolina:
. . . the situation with respect to the experimental checks of general relativity theory is not much better
than it was a few years after the theory was discovered, [despite] the tremendous improvement in
experimental techniques now available.
Figure 2
From left to right: Dave Wilkinson, Bob Dicke, Ed Groth, and Jim Peebles in February 1981 on the occasions of a celebration of the
Nobel Prize award to Val Fitch and Jim Cronin and the Ph.D. defense of Bob’s student Jeff Kuhn. At the time, Bob was working on
tests of gravity physics, paying particular attention to the possibility that a rapidly rotating solar interior produces a gravitational
quadrupole moment large enough to affect the orbit of Mercury and the test of general relativity. Dave was occupied by ground-,
balloon-, and satellite-based measurements of cosmic radiation backgrounds. Ed and I with other colleagues had nearly completed a
program of statistical measures of how galaxies are distributed, oriented, and moving. Some of this is mentioned in Section 4. I had just
published The Large-Scale Structure of the Universe (Peebles 1980) and was exploring ideas of how structure formed. The photograph
was taken by Jane Groth.
Dicke mentioned work in his group to improve the experimental situation: a pendulum for preci-
sion measurement of the gravitational acceleration to “detect possible annual variations” (Dicke
1957, p. 59), a repetition of Eötvös’s demonstration that the gravitational acceleration of a free
test particle depends very little on its composition, and the development of a rubidium atomic
clock to better measure the motion of the Moon. He recalled Dirac’s idea that the strength of the
gravitational interaction might be small because it has been decreasing for a long time, and Bob
discussed how that might be detected. At the conference John Wheeler discussed the argument
by Ernst Mach and others in the nineteenth century that inertial motion has relative rather than
absolute meaning. Bob did not mention this, but we heard a lot about it in his gravity group
meetings.
4 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
The gravity research group discussed a wonderfully free flow of ideas to be explored by experi-
ments, mostly benchtop, consultation of the literature on geology and astronomy, and exploratory
calculations of the sort I enjoy. I published a few, some with Bob. He decided I was better suited
FTBB: Finding the Big
to his new research interest rather than particle physics and invited me to look into constraints on Bang (Peebles, Page &
possible variations of the strength of the electromagnetic interaction, which usually is represented Partridge 2009)
by the fine-structure constant α. A change of α over geological times would have different effects
on the Coulomb energies of different atomic nuclei, depending on their charges and structures,
which would change decay rates, resulting in inconsistent radioactive decay ages derived from
different isotopes and different alpha, beta, and fission decays. I noticed that the β-decay energy of
187
Re to 187 Os is particularly small, so the decay rate is particularly sensitive to the value of α. My
source, Strominger, Hollander & Seaborg (1958), lists only rough estimates of the decay energy,
so I drew up a plan to measure it. Bob suggested I stay with theory. My dissertation included a
covariant gravity theory with variable α that does not violate the Eötvös experiment, and I obtained
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
a limit from terrestrial and meteorite radioactive decay ages on the change of the fine-structure
constant, |δα/α| < 1 × 10−3 in the past 4.5 Gyr. Current bounds are two orders of magnitude
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
The path was not straight. In Relativity, Thermodynamics and Cosmology, Tolman (1934, p. 363)
wrote that homogeneity
In my edition of The Classical Theory of Fields (Landau & Lifshitz 1951, the English translation of
the 1948 Russian first edition), a footnote on page 332 cautions that although observations offer
some basis for the assumption of homogeneity, “it remains an open question whether this situation
will not be changed even qualitatively as new data are obtained.” This sensible remark is the more
striking because the magisterial series of books on theoretical physics by Landau & Lifshitz is
almost entirely free of phenomenology. In Cosmology, Bondi (1952, p. 168) warned that it is “quite
possible that looking into more distant spaces will completely upset present concepts.” There
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
was some evidence of homogeneity: Hubble’s (1936) galaxy counts scale with apparent magnitude
about as expected if the galaxy distribution were close to homogeneous, and the linear relation
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
between galaxy redshifts and distances (Hubble & Humason 1931) agrees with homogeneous
and isotropic expansion. These measurements probe to recession velocities ∼10% of the velocity
of light, an impressively large distance but far from the edge of the observable Universe. In an
authoritative assessment of the observational situation, Oort (1958) stressed that there are large
fluctuations in the counts of galaxies across the sky, but he was willing to conclude his review with
an estimate of the cosmic mean mass density, assuming there is one. Tolman (1934), Landau &
Lifshitz (1951), and Bondi (1952) went further, working through detailed analyses of properties of
homogeneous world models. That was fair enough: It was a working hypothesis, but little more.
Gérard de Vaucouleurs (1970) was right to emphasize this, but tests were emerging. For exam-
ple, in Klein’s (1966) and Alfvén’s (1966) world picture, a bounded cloud of galaxies is expanding
into empty flat space-time. Velocity sorting would put the faster-moving galaxies further away,
producing Hubble’s linear velocity–distance relation. The highest velocities would have to exceed
a tenth of the velocity of light to fit Hubble’s counts and Hubble & Humason’s (1931) redshift
measurements. That would be quite an explosion, but you might imagine it. The problem for this
picture is that we are in a sea of radiation that is close to uniformly bright across the sky. The
X-ray background was known to be isotropic to a few percent (Schwartz 1970), and the microwave
background was known to be isotropic to one part in 103 (Partridge & Wilkinson 1967, Conklin
& Bracewell 1967). Because space was known to be transparent to microwave radiation—radio
sources were observed at high redshifts—radiation produced in the Klein-Alfvén explosion at the
start of the expansion would have long since left the cloud of galaxies. If the radiation we observe
came from the galaxies, then we would have to be at the center of the cloud to a part in 1,000,
and the cloud would have to be isotropic to like accuracy, which seems absurd. If instead the sea
of radiation uniformly filled the space into which the galaxies are moving then our Galaxy would
have to be moving relative to the sea at less than about 0.001 times the velocity of light while
the vast majority of galaxies—apparently equally suitable homes—are moving relative to the sea
much faster, up to a tenth of the velocity of light. That also seems absurd. The more reasonable
interpretation is that, contrary to Klein and Alfvén, the radiation is flowing with the expansion
of the homogeneous Universe of galaxies (Peebles 1971, p. 40). I used to hear the complaint that
this draws a large conclusion from limited and indirect data, but that is how science operates. The
case for homogeneity is still tighter when there is independent evidence for the same conclusion,
of course. We have that, as is discussed below.
Another alternative to Einstein’s homogeneity picture is an unbounded clustering hierarchy, a
fractal galaxy distribution. The idea has a long history (Bondi 1952, p. 14), but Benoit Mandelbrot
6 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
(1975), in Les objects fractals, made the community pay attention. He argued for a scale-invariant
galaxy distribution with fractal dimension D = 1, meaning the typical number of galaxies within
distance r of a galaxy scales as N (<r) ∝ r D ∝ r, because that makes the gravitational potential
Cosmic microwave
GM/r independent of length scale: There is no relativistic divergence of space-time curvature on background (CMB):
large or small scales. At the time, the statistical analyses recalled in Section 4 showed that on scales 3K cosmic thermal
of less than about 10 Mpc the low-order galaxy position correlation functions agree with D = microwave
1.23. More importantly, Mandelbrot would have it that the fractal distribution continues to the background radiation
largest observable scales, while I interpreted the evidence to be that on scales larger than 30 Mpc
the fractal dimension is close to D = 3, which is statistical homogeneity. Consider that in a scale-
invariant fractal Universe the fluctuations δN/N of galaxy counts across the sky are independent of
the depth to which they are counted (assuming a sensible nonsingular distribution of distances of
galaxies selected by apparent magnitude). This is what scale invariance means. It predicts that the
large fluctuations across the sky in the Shapley & Ames (1932) catalog of bright nearby galaxies
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
are repeated in the Lick counts of fainter more distant galaxies (Shane & Wirtanen 1967, Seldner
et al. 1977). That is quite contrary to the galaxy maps I used to show fractal advocates (usually to
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
which is the idea that led Einstein to argue for a homogeneous Universe, might the Universe have
other effects on local physics? Perhaps motion relative to the rest frame defined by the Universe
has an effect on laboratory physics in some more subtle way than the velocity of light, maybe local
gravity physics. Perhaps cosmic evolution affects laboratory physics, maybe changing the strength
of gravity or the electromagnetic interaction. Bob was attracted to this line of thought because
it suggests new and maybe informative measurements, which is the point of natural science, and
because he shared with Dennis Sciama the philosophy of the Unity of the Universe (Sciama 1959).
Donald Lynden-Bell (2010, p. 1) began his contribution to this series with “Is space-time only
brought into being by its energy content? The jury is still out.” Examples of variants of Mach’s
philosophy are in Are the Fundamental Constants Varying in Spacetime? (Molaro & Vangioni 2009).
Philosophy this durable can’t be all bad. (d ) An even more durable idea is that the world operates
by rules we can discover, in successive approximations to a reality thus defined. We are surrounded
by results of this spectacularly productive philosophy. I like the example of homogeneity, a simple
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
result established in a roundabout way. I continue with this line of thought in Section 6.
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
8 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
temperature scales at T ∝ λ−1 ∝ a(t)−1 . Because this temperature scaling applies to all modes,
the spectrum remains thermal as the radiation cools. Tolman (1934, p. 427) showed this another
way. Interaction with matter can change the CMB spectrum, but the CMB energy density is so
large that a significant perturbation to the spectrum is difficult to arrange.
George Gamow found an application. In papers leading up to Alpher, Bethe & Gamow (1948),
he developed the idea that the chemical elements were built up by successive neutron captures
during the early rapid expansion of a Friedmann-Lemaı̂tre Universe. In his Ph.D. dissertation,
Gamow’s student Ralph Alpher (1948) saw an inconsistency. The condition for appreciable but
not excessive production of heavier elements in expansion time t at baryon number density n is
σ vnt ∼ 1, where σ is the radiative neutron capture cross section and v is the relative velocity.
In relativistic cosmology, the expansion time is t ∼ (Gρ)−1/2 . It was natural to take the mass
density to be ρ = mn, where m is the baryon mass. These conditions, with measured σ v, make t
much larger than the neutron half-life. One could imagine neutrons were created at the Big Bang,
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
along with everything else, but where would neutrons have come from so late in the career of the
Universe? Gamow (1948) did not acknowledge the problem but proposed a brilliant solution: that
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
the early Universe was so hot that the mass density was dominated by thermal radiation, which
would speed the expansion. [Alpher, Bethe & Gamow (1948) had noted another advantage of a
hot Universe: Hot neutrons avoid low-energy resonances that would spoil the Gamow-Alpher
anticorrelation of σ v with element abundance.] Gamow analyzed the first step in the proposed
buildup of the elements, which is radiative capture of neutrons by protons to form deuterons,
n + p → d + γ . At high temperature, the reverse process, d + γ → n + p, suppresses accumulation
of deuterons. When the temperature falls to T 109 K, the photon energies become too low for
photodissociation, so deuterium can accumulate and burn to heavier elements by faster reactions
like d + d → 3 He + n. The temperature T 109 K fixes the thermal energy density, which
determines the expansion time, t ∼ (Gρ)−1/2 ∼ 100 s, and the condition σ vnt ∼ 1 with Gamow’s
estimate σ v ∼ 10−20 cm3 s−1 gives baryon number density n ∼ 1018 cm−3 at T ∼ 109 K. This
with refinements has become the well-tested theory of the origin of helium and deuterium.
Gamow certainly understood that the thermal radiation would remain after element formation.
He saw (I think by intuition, not calculation) that the epoch of equal mass densities in matter and
radiation marks the onset of the gravitational growth of structure, and he found the temperature
then, ∼103 K (Gamow 1948). Gamow did not compute the present temperature at his usual
estimate of the present mass density, mn = 10−30 g cm−3 , but the easy exercise gives T ∼ 10 K.
Alpher & Herman (1948) took this bold step. The Dicke et al. (1946) experiment showed T <
20 K. Bernie Burke’s judgment in FTBB is that this experiment could have been improved to
detect the CMB at T = 2.7 K in the 1950s. But it was not attempted until 1964.
The first indirect detection of the CMB is the observation of absorption of starlight from the
first rotationally excited level of the interstellar molecule CN. McKellar (1941) converted the ratio
of populations in the excited and ground levels to an effective thermal excitation at temperature
T = 2.3 K. Herzberg (1945, p. 496) wrote that this has “a very restricted meaning.” I expect
he meant the excitation could be caused by particle collisions rather than radiation. In FTBB,
George Field recalls concluding that particle excitation is unlikely because radiative relaxation of
CN to the ground level is fast, but he did not publish. Nick Woolf recalls mentioning the CN spin
temperature to Dicke when he was considering a search for the CMB, but Bob didn’t reply; maybe
he didn’t understand. Fred Hoyle (1950) understood. In a review of Theory of Atomic Nucleus and
Nuclear Energy-Sources (Gamow & Critchfield 1949), Hoyle (1950, p. 195) wrote that Gamow’s hot
Big Bang theory, which is presented in an appendix, would lead to a “temperature of the radiation
at present maintained throughout the whole of space much greater than McKellar’s determination
for some regions within the galaxy” (Hoyle 1950, p. 195). Hoyle recalled Gamow (in about 1956)
“explaining his conviction that the universe must have a microwave background” (Hoyle 1981,
p. 522) and his telling Gamow that McKellar “set an upper limit of 3 K for any such background”
Hoyle also recalled that he “missed it again in exactly the same way in a discussion with Bob Dicke”
(Hoyle 1981, p. 522) in 1961. I have seen no evidence that Bob remembered this discussion, or
that Hoyle influenced Bob’s thinking in some more subtle way, or that anyone noticed that the
CN spin temperature might be a signature of a Hot Big Bang until after direct detection of the
CMB [at T = 2.7 K, which is consistent with McKellar (1941) within the uncertainty].
A second signature is the high cosmic abundance of helium. At a summer school at the Uni-
versity of Michigan, Gamow (1953a) explained that, while element buildup in a Hot Big Bang has
problems producing heavier elements, it readily yields helium abundance ∼30% by mass, “in good
agreement with the observed relative amount of Hydrogen and Helium in the universe” (Gamow
1953b, p. 12). He proved to be right, but I am not sure how convincing those observations were,
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
and Gamow tended not to bother with documentation. Geoff Burbidge and Don Osterbrock were
at the summer school. (They are pictured with Gamow in the conference photograph taken by
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Ed Spiegel and published in figure 1 in Gingerich 1994). Burbidge (1958) later concluded that the
galaxy contains at least 10% helium by mass, and that this large amount is difficult to understand at
the present rate of nuclear burning. He did not mention a Hot Big Bang. Osterbrock & Rogerson
(1961, p. 133) put the helium mass fraction at about 32%, noting that this
. . . could be at least in part the original abundance of helium from the time the universe formed, for the
build-up of elements to helium can be understood without difficulty on the explosive universe picture.
Their reference for the explosive Universe picture was Gamow (1949). Hoyle & Tayler (1964)
knew the cosmic helium abundance is high and naturally explained by a Hot Big Bang, or, they
pointed out, by many little bangs. The spectrum of the radiation accompanying the process would
distinguish them. Tayler recalled that early drafts of Hoyle & Tayler’s paper mentioned the
thermal radiation in a Hot Big Bang, but “with such a low temperature that it was not surprising
that it had not been discovered” (Tayler 1990, p. 372).
In FTBB Yuri Smirnov recalls that Yakov Borisovich Zel’dovich (who is pictured in
Figure 3) thought that the cosmic helium abundance is low, perhaps 2.5% by mass. Zel’dovich
(1963, p. 372) wrote that “We cannot make any estimate of the reliability of these results,” but
he went on to show how a Cold Big Bang could avoid producing helium: In a degenerate sea
of equal number densities of protons, electrons, and neutrinos, the neutrino degeneracy energy
would prevent formation of neutrons by inverse beta decay.
Key information may be overlooked, as happened to Zel’dovich, or forgotten. Hoyle & Tayler
seem to have forgotten the CN temperature Hoyle mentioned to Dicke, who forgot the conver-
sation. We had to remind Dicke that he already had placed a bound on the CMB temperature,
T < 20 K (Dicke et al. 1946). In an analysis of the structure of Jupiter, I used the Osterbrock &
Rogerson helium abundance in an estimate of the equation of state (Peebles 1964). The explo-
sive Universe picture meant nothing to me then, and I had forgotten about the Osterbrock &
Rogerson helium abundance when Dicke suggested I look into theoretical implications of a
radiometer search for the CMB.
Thoughts of a Hot Big Bang made me think of an exploding pressure cooker, which made me
think of out-of-equilibrium particle reactions and light element formation. I didn’t know about
Zel’dovich, but I, too, wanted to prepare for a low helium abundance or a null result from the
radiometer measurement, which led me to think of effects of degenerate neutrinos. I talked about
Hot and Cold Big Bangs in a colloquium at Wesleyan University on December 2, 1964. There
is no evidence in my notes that I was aware of Gamow’s work. A letter dated February 1, 1965
10 Peebles
AA50CH01-Peebles
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org ARI 16 July 2012 11:15
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Figure 3
A photograph I took of Zel’dovich and his wife in Hungary in 1987 on the occasion of IAU Symposium 130.
I met Zel’dovich only twice, at this meeting and in 1977 at IAU Symposium 79 in Estonia. I was not eager to
travel to the Soviet Union, and Zel’dovich could not travel out of it, so we tended to work in parallel, which
can be productive. At the 1987 symposium, the CDM model was generally accepted as a promising
framework for analyses of structure formation, and simulations presented at the meeting suggested galaxies
were assembled by mergers and accretion at low redshifts. I don’t remember what Zel’dovich thought. I was
skeptical and still see problems with the successor, CDM, on this score.
to Hoyle and Tayler shows I had learned I was reinventing the wheel. I set out to do a more
careful computation, as did Bob Wagoner with Fred Hoyle and Willie Fowler. Bob Wagoner and
I compared notes at a December 1965 conference at the Fontainebleau Hotel in Miami Beach,
Florida. Not atypical for the time was the fancy setting—we shared conference facilities with the
Jewish Funeral Directors of America—and the participation of a good fraction of the cosmology
community, apart from the usual absence of people from the USSR. Computation of element
buildup in a Hot Big Bang is complicated by the extreme changes of the various nuclear reaction
times relative to the cosmic expansion time as the model Universe expands and cools. I learned
from Wagoner that Hoyle and Fowler knew about this from analyses of stellar evolution. I relied
on carefully checked numerical fixes. I was ready to submit a paper (Peebles 1966) that takes
account of the relatively few reactions relevant for helium and deuterium. Wagoner, Fowler &
Hoyle (1967) looked at more reactions relevant for trace production of elements heavier than
helium.
The evidence reviewed in FTBB is that the first serendipitous detection of the CMB was by
communication systems at the Bell Telephone Laboratories. Dave Hogg describes in FTBB the
several generations of receiving systems with low-noise traveling wave maser preamplifiers, and
he and Arno Penzias and Bob Wilson recall the increasingly awkward problem of reconciling
detected and expected system noise. Doroshkevich & Novikov (1964, p. 113) were the first to
point out the relevance of the Bell work reported by Ohm (1961) “for experimental checking of
the Gamow theory.” They were not aware that Bell had an excess noise problem, and the Bell
people were not thinking about a sea of microwave radiation.
In 1964, at Bob Dicke’s suggestion, Peter Roll and Dave Wilkinson built a Dicke radiometer
to look for the CMB. I was the unwitting medium for recognition that Bell had the signal Roll
and Wilkinson were seeking. I had asked Dave if I might mention the Roll-Wilkinson experiment
in colloquia. He saw no problem: “no one could catch up with us now.” At the Applied Physics
Laboratory at Johns Hopkins University (where, I later realized, Alpher and Herman were based in
1948), on February 19, 1965, I mentioned the experiment. Ken Turner from Princeton graduate
student days attended the colloquium. He recalls in FTBB telling Bernie Burke, who advised
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
Penzias to talk to Dicke. That led to recognition of a new type of cosmic radiation, the CMB
(Penzias & Wilson 1965), maybe a fossil from a Hot Big Bang (Dicke et al. 1965). The reaction I
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
recall at Princeton was excitement that there actually is something to detect, as Roll & Wilkinson
(1966) did a few months later, maybe from the Hot Big Bang, but at least something to be analyzed.
The Nobel Prize Committee was right to recognize Penzias & Wilson: They demonstrated that
Bell had a problem, and they complained about it until the community responded. But Dicke
should have been included. He invented key technology, and he initiated the experiment that led
to the recognition of this most informative fossil from the Big Bang.
Despite Gamow’s elegant Hot Big Bang theory, the high spin temperature of interstellar CN,
the high cosmic abundance of helium, and the unexplained radiation in the Bell communications
systems, the catalyst for recognition of the CMB came from yet another direction: Dicke’s idea
of entropy from a bouncing Universe. I hope we cosmologists are not exceptionally obtuse. The
extenuating circumstance I see is that there were too few actors to encounter each other often
enough to be forced to bear in mind many points of view. In this major advance in slow motion, it
is easy to see the missteps. I would like to think that if you examined major advances in any other
natural science in comparable detail you would see similar confusion.
12 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
function may be horribly complicated. Tukey (1966, p. 27) put it that “the tragic accident that
killed H. R. Seiwell . . . destroyed the only man who could usefully look at a plot of autocorrelation
against lag.” I had been hearing similar comments in Dicke’s group. So I started with power spectra.
We might pause to note that measured galaxy distances are not accurate enough for small-
scale mapping of the spatial distribution (except in our immediate neighborhood), so one instead
measures angular distributions in samples selected by some measure of distance, such as apparent
magnitude, and translates to spatial statistics under the assumption we have a fair sample of
a statistically homogeneous and isotropic point process in an imagined statistical ensemble of
universes. The two-point spatial and angular correlation functions, ξ (r) and w(θ), are defined by
the probabilities of finding a point in the volume element δV at distance r from a randomly chosen
point, and of finding a point in the element of solid angle δ at angular distance θ from a point,
The mean number densities of points are n per unit volume and N per steradian. The analog of
a Fourier representation for point positions θi , φi on the sky is
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
1
alm = Y lm (θi , φi ), ul = |alm |2 /N − 1 = 2π N w(θ)Pl (cos θ)d cos θ. (2)
i −1
The second expression is the ensemble average relation between the angular power spectrum
ul and the angular correlation function w(θ), where Pl is a Legendre polynomial. The distances
between zeros of the real and imaginary parts of the spherical harmonic Y lm are close to π/l
(except near the poles where Y lm → 0). This corresponds to the spacing π/k of the zeros of the
real and imaginary parts of e i k·r . If a homogeneous and isotropic process is sampled at different
effective depths, d, with distributions of distances r that scale as dP/dr ∝ f (r/d ), as in selection
by apparent magnitude, and clustering is measured on scales much less than d, then the angular
functions scale with depth as
with N ∝ d 3 . These scaling relations were of prime importance for a check of statistical ho-
mogeneity, and, more important, a test for degradation of the data by systematic errors such as
irregular local obscuration or inhomogeneous sampling. The limited part of the sky in an astro-
nomical sample requires a correction to the estimates of ul and introduces a correlation among
the alm . This and more is in a paper (Peebles 1973) with the pretentious subtitle “Paper I,” though
I don’t think I had any idea how many more there might be. (The last is Paper XII, when work in
statistics was in directions quite different from Paper I.)
John Tukey was at Princeton, and in November 1972 I presented a seminar for his group.
He, Geoff Watson, and the others were encouraging but we agreed that there was no need to
collaborate. The complications to be feared were in the astronomical data. As Tukey might have
expected, power spectra are the better statistic for the large-scale galaxy distribution and the
CMB anisotropy. But I found that the two-point and higher order correlation functions are better
suited to the strongly nonlinear clustering of galaxies on scales 10 Mpc. I have to mention Nelson
Limber’s (1954, p. 656) account of estimating the correlation of Lick galaxy counts:
The number of nebulae per square degree along [a parallel of latitude] was recorded separately on each
of two strips . . . one strip was displaced φ degrees with respect to the other, and then the values on the
two strips which were adjacent after the displacement were multiplied together, and the mean value of
these products were obtained.
Jer Yu and I had the great advantage of programmable computers when we introduced the use of
power spectra and spherical harmonic expansions.
satellite). Paper II (Hauser & Peebles 1973) showed that the spectrum estimates for the Abell
cluster positions at fixed l and different m have the exponential distribution P (ulm > x) ∼ e −x/β(l)
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
that Jer Yu pointed out follows from the central limit theorem. We better understood the scaling
relations (Equation 3), and we found the predicted scaling of ul and w(θ) with limiting Abell
distance class. This evidence that our statistics were reliable measures seems trivial now, but it was
really encouraging then.
At the Lick Observatory in California, Donald Shane led a major survey of the large-scale
galaxy distribution (well described by Vasilevskis & Osterbrock 1989). Shane and Carl Wirtanen,
with help from others, counted 1.2 million galaxies (some more than once where fields overlap;
there are 800,000 different galaxies) brighter than apparent magnitude mpg 19 and north of
declination δ = −23◦ . Galaxies were identified by eye using a traveling microscope to scan 17-
inch square photographic plates, each covering a 6◦ × 6◦ field. The counts in each of the 1,296
10 arcmin × 10 arcmin cells in each field were recorded by hand. A sample data sheet is shown
Supplemental Material in Supplemental Figure 1 (follow the Supplemental Material link from the Annual Reviews
home page at http://www.annualreviews.org). The counts commenced in late 1947 and were
completed in 1959. Donald is on my list of heroic figures.
Shane & Wirtanen (1967) published the counts summed in 1◦ × 1◦ cells, and Totsuji & Kihara
(1969) used these data to show that the galaxy spatial two-point function is close to a power law,
14 Peebles
AA50CH01-Peebles
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org ARI 16 July 2012 11:15
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Figure 4
A photograph I took of Donald Shane when he visited Princeton in about 1976. Some of the people who
worked on the Lick data reduction are Jim Fry (center left), Ray Soneira (center right), and Mike Seldner ( far
right).
the data took four years). Despite astronomers’ tendencies to enter the occasional inscrutable dot
or colon, we were happy with the data. Our computer was modest compared to what is on my
desk now, but far superior to when the counts commenced. Donald could have had no reason
to expect the 10-arcmin data ever would be fully analyzed, but he had the courage to take great
care, including the valuable allowance of generous plate overlaps, 1◦ strips on the four sides of
each 6◦ × 6◦ field. The independent counts of the same parts of the sky by different observers,
at different times, and maybe different plate emulsions enabled detailed checks for systematics [as
discussed by Shane & Wirtanen (1967); in more detail by Seldner et al. (1977); and still further
in Groth & Peebles 1986a,b)]. The offer by Seldner et al. (1977, p. 253) of the raw and corrected
data “on receipt at Princeton of an IBM compatible 2,400-ft magnetic tape” pleased Donald: He
considered the Lick data published at last (though now the several generations of storage of these
data are obsolete).
I earlier made maps from Zwicky’s catalog of nearby galaxies and sent him copies. I received
a polite reply from a great but not always happy scientist. I used a mechanical plotter driven
by the university’s computer (a big operation then) to represent the Lick 1◦ counts as squares
with sizes proportional to the counts, and Alison and I sat at our dining room table to black in
the squares. This map is the frontispiece of The Large-Scale Structure of the Universe (Peebles
1980). Supplemental Figure 2 is a better copy. The 10-arcmin Lick maps were made by a film
scanner at the Princeton Observatory that was acquired and modified by John Lowrance to read
digital data as well as plot it. This was part of Lyman Spitzer’s goal of a space telescope. Heiles &
Jenkins (1976) describe the scanner and show its application to astronomical data. Supplemental
Figure 3 shows improved digital versions of the 10-arcmin Lick maps.
It was a pleasure to bring Donald the 10-arcmin Lick maps in 1976 and ask, Is this what you
saw? I believe I have an accurate memory of his answer: “I was looking at this one galaxy at a
time.” Stewart Brand, editor of the Whole Earth Catalog in its many editions, contacted me with
the proposal to publish the Lick map as a poster. I never thought to ask how he knew about the
map. The poster was published in 1978. We received no royalties, but enjoyed seeing this elegantly
printed poster taped to office walls.
It was good also to find a positive scaling check of consistency of the two-point functions from
CDM: CDM
adjusted to include a the 10-arcmin Lick counts, from the 10,000 galaxies cataloged in a single, much deeper 6◦ ×
cosmological constant 6◦ field by Rudnicki et al. (1973), and from the shallower Zwicky sample (figures 13 and 14 of
that allows lower mean Groth & Peebles 1977). We had moved on to the three-point function (Peebles & Groth 1975;
mass density while while I enjoyed the hospitality of UC Berkeley), which was harder to measure, and the scaling
keeping flat space
check was less detailed but still meaningful. Jim Fry worked through the four-point function (Fry
sections
& Peebles 1978). No one volunteered for the five-point function. At the range of galaxy separations
we could measure, 0.03 r 10 Mpc, the galaxy two-point function ξ (r) is close to the power
law in Equation 4, and the three- and four-point functions are well approximated by products of
two-point functions in the form of a fractal with dimension D = 3 − γ = 1.23.
Ray Soneira and I checked this picture by making model galaxy distributions to compare to the
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
Lick map. We placed galaxies in randomly placed fractals of dimension D = 1.23 with bounded
size and a spread of richnesses (constrained by the two-, three-, and four-point functions), and
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
used a galaxy luminosity function to project the model distributions onto the sky at the depth
of the Lick catalog. We at first used fractals far broader than 10 Mpc, meaning the power-law
forms of the correlation functions extend to a large separation, which I supposed wouldn’t matter
because the correlation functions at large separations are small. The maps showed I was wrong:
They were, in Ray’s term, “blotchy.” This is because the mean number of model galaxies in excess
r
of random within distance r of a randomly chosen galaxy scales as 0 ξ (r)d 3 r ∝ r D = r 1.23 . If
the power law continues to large r it makes this excess large, placing the galaxies in disgusting
blotches. I checked that by cutting a model map into squares of size ∼5◦ (linear size ∼20 Mpc) and
gluing the squares onto a sheet of paper at random. That removed the large-scale correlation and
made a more pleasing map. A two-point function ξ (r) ∝ r −γ at r < 20 Mpc truncated to ξ (r) = 0
at r ≥ 20 Mpc produces an angular function that slowly bends away from a power law, contrary
to the measurements. To remedy that we made ξ (r) rise above the power law at r ∼ 10 Mpc and
then fall below it at r ∼ 20 Mpc, producing a shoulder. The functional form in figure 6 of Soneira
& Peebles (1978) shares these features with modern measurements (Zehavi et al. 2004).
Our argument for the existence of the shoulder in ξ (r) was sound, I believe, but too qualitative
to have much chance of being influential. I nevertheless count it as an example of how qualitative
impressions can complement objective measures. A more ambiguous example is my impression
that linear features in the Lick Map may only illustrate the tendency to see patterns in noise.
Indeed, our model map presents linear features that are in the eye of the beholder; they were not
built in. I don’t know whether some statistical measure could have deduced the existence of walls
of galaxies from the Lick data: recall the danger of designing a statistic to respond to what you
have seen. The frothy character of the galaxy distribution seen in dense redshift surveys (Davis
et al. 1982; de Lapparent, Geller & Huchra 1986) is a qualitative impression, and real.
Qualitative impressions can inspire useful quantitative models. Jerzy Neyman and Elisabeth
Scott proposed to use the galaxy N-point correlation functions to constrain a model that in its
original version (Neyman & Scott 1952) places all galaxies in clusters (now termed halos). The
number of halo members is a random variable, the position of each galaxy within its halo is a
universal random variable, and the halos are placed at random. The modern halo occupation
distribution (HOD) framework generalizes this: The halo mass is a random variable matched to
the dark matter halo mass function in the CDM cosmology (discussed in Section 5); the number
of galaxies in a halo and their luminosities, morphologies, radial distributions, and motions are
random variables that depend on the halo mass; and the halo space distribution is modeled on
the CDM dark matter power spectrum (Cooray & Sheth 2002, Berlind & Weinberg 2002).
16 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
The successful HOD interpretation of the galaxy two-point correlation function ξ (r), from pairs
of galaxies in the same halo and in different halos, and the variation of the two-point function
with galaxy luminosity and redshift, establish HOD as a useful representation (Zehavi et al. 2004;
Einstein-de Sitter
Conroy, Wechsler & Kravtsov 2006). There is no arguing with success, but qualitative impressions (EdeS) model:
tell me there is more to the story. Neyman, Scott & Shane (1956, p. 96) wrote that, relativistic
Friedman-Lemaı̂tre
Because of the observations, originally made by Shapley . . . and more fully confirmed by Shane and cosmological model
with no space
Wirtanen . . .that the clusters of galaxies tend to congregate in small groups, it appears interesting to
curvature and zero
deduce formulas allowing for the possibility of “superclustering” . . . . cosmological constant
The Nearby Galaxies Atlas (Tully & Fisher 1987) beautifully illustrates the local situation. We are
in a distinctly nonlinear flat concentration of galaxies, the supergalactic plane. In addition to the
galaxies around the Milky Way and M31 in the Local Group, the supergalactic plane contains the
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
M81 group, the galaxies around Centaurus A, and the more massive concentrations around Maffei
and IC342. Just above the supergalactic plane is a very empty part of the Local Void (illustrated
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
in figure 1 of Peebles & Nusser 2010). And, with Peter Shaver (1991), I expect it is significant
that the nearest six Abell clusters are close to the supergalactic plane (figure 3.7 of Peebles 1993).
This hierarchy of nonlinear structures includes the supergalactic plane with diameter ∼100 Mpc
that contains groups and clusters, some of which contain substructures, that contain the structures
observed on scales ∼10 kpc as galaxies. I doubt this range of nonlinear structures can be accom-
modated within HOD; it speaks of some version of the clustering hierarchies Neyman, Scott,
Shane, Soneira, and I modeled. The issue is recognized in a modified HOD that places galaxies
in subclusters in clusters (Sheth & Jain 2003, Giocoli et al. 2010). I expect HOD will continue to
evolve into some kind of clustering hierarchy over the considerable range of scales of structures
seen in our neighborhood, maybe to be constrained by the higher order correlation functions
Neyman, Scott, and I so valued.
In the early days, before establishment of the gravitational instability picture for cosmic struc-
ture formation that inspired HOD, I liked to point out that if structure formed by explosions that
pushed matter into piles, moving matter by distances ∼R, it would make ξ (r) negative at r ∼ R.
The inverse square law of gravity allows gravitational growth of mass concentrations that nowhere
lowers the density (in the ensemble average), so it need not drive ξ (r) negative anywhere, in line
with what we were finding (Peebles 1974a). This is no surprise now, but it was interesting then.
I also liked the evidence that the small-scale, nonlinear clustering of galaxies has N-point
correlation functions that are close to scale-invariant power laws. Classical gravity has no scale.
Among the Friedman-Lemaı̂tre models, Einstein-de Sitter (hereafter EdeS) has no scale. What
could be more natural than scale-invariant clustering of matter in a scale-invariant Universe
(Peebles 1974a,b; Davis, Groth & Peebles 1977)? Scales are built in by the interaction of matter
and radiation in the early Universe (Section 5) and by the transition from linear fluctuations on
large scales to small-scale nonlinear clustering, but the thought was that fully developed nonlinear
clustering may forget initial conditions, in the manner of fully developed turbulence. Nature had
other ideas.
Recall that the expansion rate in the relativistic Friedman-Lemaı̂tre models satisfies
2 a 3 a 2
ȧ o o
= H o2 m + k + , m + k + = 1. (5)
a a a
Hubble’s constant is H o , the mean separation between galaxies scales as a(t), with present value a o ,
the constant m is a measure of the mean mass density in matter (radiation may be ignored here),
k represents the effect of space curvature, and represents Einstein’s cosmological constant, .
18 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
relative line-of-sight peculiar velocities is exponential, and it indicated m = 0.60 ± 0.25 (Peebles
1981a).
These early analyses used galaxies selected by apparent magnitude. It is better to select to fixed
CfA sample: the
minimum luminosity within a volume large enough to approach a fair sample. Kirshner, Oemler pioneering redshift
& Schechter (1978) took a first step in this direction in samples of eight deeper fields, from which catalog obtained at the
I found m = 0.4 ± 0.2 (Peebles 1979). On p. 734, I entered an enthusiastic conclusion: CfA
It is remarkable that a sample of just 166 galaxies gives . . . such a clear presentation of peculiar redshifts.
This suggests that a considerable improvement of the situation should come fairly soon, with the larger
“random” redshift surveys now in progress.
I was thinking of Marc Davis, who was well aware of the value of a better redshift survey from
life as a graduate student in the Princeton gravity group. He moved to Harvard where he saw
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
that the Harvard-Smithsonian Center for Astrophysics (CfA) had the needed facilities. The CfA
sample (Davis et al. 1982) gave m = 0.2e ±0.4 (Davis & Peebles 1983). There are much better data
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
now—the two-point correlation function in redshift space in our figure 4 is a crude approximation
to the beautiful result found by Peacock et al. (2001)—but the situation was clear: Either m is
well below unity or mass is less strongly clustered than galaxies.
Reactions to this result were conditioned by the inflation picture (Guth 1981), which offers an
elegant explanation of large-scale homogeneity: Early rapid expansion driven by a temporarily very
large (as ideas in particle physics suggest is not unreasonable) could stretch wrinkles to scales too
large to see, and stretch away space curvature in the process. Cosmological space curvature thus
fell out of vogue. Because Einstein’s cosmological constant in the present Universe never seemed
reasonable to many, EdeS ( = 0, m = 1) was promoted to a general community standard. But
if m = 1, then dynamical mass estimates must be biased low. The generally accepted idea was
that this is because galaxy formation was more efficient where the mass density was higher, leaving
low-density regions with few galaxies and a lot of mass that would not be counted in the mass
concentrated around galaxies. This biasing idea had to be explored, but it had a problem from the
start. Galaxies that did manage to form in supposedly inhospitable low-density regions would be
predicted to bear the stigmata of a troubled youth, tending to be small and irregular. That is not
what the CfA sample maps show [see figure 2a and 2d of Davis et al. (1982)]. It should be noted
that the rare most luminous cD galaxies are more strongly clustered, along with the rich clusters
of galaxies where they are found, but recent data confirm that normal galaxies, such as the Milky
Way, and fainter ones, have quite similar distributions (Peebles 2001) and are useful mass tracers.
Related to the old biasing picture is the idea that galaxies are affected by environment. Numerical
simulations indicate this and, indeed, cDs prefer dense environments, whereas gas-rich galaxies
avoid them. But this line of thought does not so readily fit the observations that properties of
elliptical galaxies are remarkably insensitive to environment (Bernardi et al. 2006; Park et al. 2007;
Nair, van den Bergh & Abraham 2011) and that pure disk galaxies seem not to have interacted
with their environment at all (Kormendy et al. 2010)—fascinating puzzles for future work.
My reading of the evidence from the CfA sample maps along with other arguments for m < 1
is presented in Peebles (1986). I heard complaints that I only did it to annoy. I was serious, but
I have to admit that I enjoyed commenting on the struggle to save EdeS from all the contrary
evidence.
Trust in ideas, as in the compelling logic of EdeS, has a solid historical basis—consider that
general relativity passes a tight network of tests on cosmological length scales some 14 orders of
magnitude larger than the scale of Einstein’s test from the orbit of Mercury—as well as a history of
unintended consequences. I offer personal examples. Mike Seldner and I compared the measured
cluster-galaxy cross correlation function ξc g to the cluster-mass function ξcρ modeled as a limiting
isothermal sphere normalized to cluster galaxy velocity dispersions. We assumed ξcρ ξc g , and
matched the two at radius ∼3 Mpc, where the mix of morphological types is fairly close to the
field. The result, m = 0.69 ± 0.11 (Seldner & Peebles 1977), might be expected to be biased high
because of the high ratio of mass to starlight in clusters. Indeed, it is on the high side of our other
early estimates and below unity. But before the CfA sample changed my mind, I liked the EdeS
scale-invariance and told myself the uncertainties in these measurements are hard to assess. I had
earlier invited Marc Davis to join in a dynamical analysis of galaxy clustering using theoretical
methods for nonideal gases [following Bill Saslaw (1972)]. The thought was that galaxy clustering
may be simple enough to analyze, and the forms of the galaxy N-point functions suggested an
elegant closure ansatz for the moments of the distribution of particle positions and velocities (Davis
& Peebles 1977). Marc spent a lot of time on this, but the approach did not enter mainstream
cosmology. People still are working on such lines, which is good: We ought to piece together
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
the physical considerations that account for phenomena. I am reminded that Marc’s doctoral
dissertation, with Dave Wilkinson and Mike Hauser, was a search for the “primeval galaxies”
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Bruce Partridge and I proposed (Partridge & Peebles 1967, Davis & Wilkinson 1974). This was
good science too, though again maybe a quarter century ahead of its time. The consequences of
these premature ideas may have slowed Marc’s career. But maybe they helped too: He went on to
do many things to shape cosmology.
20 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
evolution of what is encoded in initial conditions (Peebles 1967). Third, the Universe is gravi-
tationally unstable: An initially quite smooth mass distribution grows clumpy. Lifshitz made this
point, as did Lemaı̂tre (1933). All this is commonplace now, but there was a time when serious
people did not agree that “the remarkable uniformity of the Universe must await a deeper theory,
perhaps a quantized gravity theory” (Peebles 1972, p. 58). I made my first numerical simulations
to illustrate this gravitational instability at the Los Alamos National Laboratory in 1969. Because
I was an alien (Canadian), I was not allowed to be alone with computers I suppose were used for
classified research, but I was allowed to use them, which may seem as archaic now as the 300 par-
ticles in the simulations. Supplemental Figure 4 is a somewhat later illustration of the evolution Supplemental Material
of our Universe from order to chaos. I used to enjoy making this point in popular and technical
lectures.
Recognition of the CMB (Section 3) suggested that matter in the nearly homogeneous early
Universe was thermally ionized and recombined (or combined) to largely neutral atomic hydrogen
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
and helium when the temperature fell to about 3,000 K. Before that, scattering of the radiation
by free electrons causes plasma and radiation to behave as a fluid, and the radiation pressure
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
prevents formation of bound nonrelativistic concentrations of matter (or, as we now have to say,
concentrations of baryons) until matter recombines and decouples from the CMB (Peebles 1965).
Residual free electrons after recombination interact with the CMB to keep the baryons warm, and
the matter temperature fixes the Jeans mass of the first generation of bound baryonic systems. Free
electrons also catalyze formation of molecular hydrogen (e + H → H− + γ , H− + H → H2 + e),
and radiation by H2 can cause the Jeans mass to shrink, maybe to a stellar mass. The situation has
been analyzed (e.g., Wise et al. 2012) in far more detail than in our early calculation (Peebles &
Dicke 1968), but we knew enough then for clear motivation for computing recombination. (I still
wonder whether structure formation really failed to leave a signature of the primeval Jeans mass,
so tantalizingly in the range of masses of globular star clusters.)
A hydrogen atom forms when a proton captures an electron, in the process releasing an ionizing
photon or a Lyman-α or other resonance photon. An ionizing photon almost inevitably ionizes
another atom. Resonance photons place atoms in excited levels that are more easily thermally
ionized, but some Lyman-α photons are redshifted out of resonance, and others are destroyed by
two-photon decay from the 2s level. Rashid Sunyaev recalls in FTBB a conversation in September
1966 that was the first step by him and colleagues in the Soviet Union toward the computation of
cosmological recombination (Zel’dovich, Kurt & Sunyaev 1968). My first step was a conversation
with Bruce Partridge in the spring of 1966, followed by a summer project for Russell Hensel, who
had just graduated from Princeton University with a major in physics. I completed the computation
in the fall of 1967 (Peebles 1968). This close timing of independent work in the US and USSR is
not surprising. Both groups knew the relevance for cosmic evolution.
There were even more independent analyses of the evolution of fluctuations in the distributions
of plasma and radiation approximated as a viscous fluid. Sunyaev in FTBB recalls the approach
taken in the Soviet Union. Richard Michie had worked out main elements of the physics by 1967
(Michie 1969), but illness prevented his work from getting past the preprint stage. Joe Silk recalls in
FTBB how he came to this problem. He published (Silk 1967). I worked it out too, and presented
it at the January 1967 Texas Symposium on Relativistic Astrophysics in New York City. Turmoil
at the publishing house prevented publication (but inclusion here as Supplemental Text 1 is
publication of sorts).
I took the next step, the radiative transfer analysis needed for a more complete computation
of the residual patterns in the distributions of matter and radiation, while on sabbatical leave at
Caltech in 1968–1969. Jer Yu, whose earlier work with me is recalled in Section 4.1, joined in the
numerical solutions (Peebles & Yu 1970). We derived what is now termed the baryon acoustic
oscillations in the galaxy power spectrum. The corresponding peak in the galaxy correlation
function is in Peebles (1981b, figure 5). Both are quite different from what is measured (Percival
et al. 2001, Eisenstein et al. 2005), of course, because we did not take account of nonbaryonic dark
CDM: EdeS with
mass dominated by matter. Jer and I gave an awkward representation of the relic CMB anisotropy, and our numerical
nonbaryonic cold dark methods were archaic, but we were close to the standard physics of analysis of large-scale structure
matter and scale- in the distributions of matter and the CMB temperature. I did not continue with this, in part
invariant adiabatic because I had trouble imagining that such tiny disturbances to the CMB could be detected and
initial conditions
in part because of recurring indications that radiation energy released during structure formation
may have caused much larger disturbances to the distributions of matter and radiation. As Dick
Bond (1990) put it (on page 51),
On fairly broad theoretical grounds we expect backgrounds such as the submillimetre excess found by
the Berkeley-Nagoya team (Matsumoto et al. 1988) to arise. Although current observations allow large
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
energy releases, the surprising part of the BN data is its apparent magnitude.
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
Bond went on to discuss how energy releases could disturb the spectrum of the radiation and its
spatial distribution, a messy situation. I preferred to think about something else and saw ready
work in measures of the galaxy distribution and motions (Section 4). It was a better subject for
me; I like the instant gratification of data.
I returned to the CMB anisotropy when two groups announced detection of a quadrupole
δT /T ∼ 10−4 (Fabbri et al. 1980; Boughn, Cheng & Wilkinson 1981). I dutifully produced a
theory to fit the measurements (Peebles 1981c). The Princeton group soon withdrew their result
and lowered the possible anisotropy. I responded with a theory that predicts anisotropy below
their new upper limit (Peebles 1982). This model assumes the mass of the Universe is dominated
by a gas of particles that interact weakly, if at all, with baryonic matter and radiation. These dark
matter particles have small primeval velocities, hence the eventual name, cold dark matter, and
the model name, CDM. The CDM model grew out of the proposal that the mass of the Universe
is dominated by a massive neutrino family (e.g., Doroshkevich et al. 1981). That soon led to many
dark matter candidates offered by ideas in particle physics (e.g., Bond, Szalay & Turner 1982;
Blumenthal, Pagels & Primack 1982; Abbott & Sikivie 1983; Preskill, Wise & Wilczek 1983).
Because cold dark matter would not be dragged by radiation, its distribution could have been
growing clumpy while the baryons were still ionized and coupled to the CMB. This is a big help in
reconciling the small CMB anisotropy with the present very clumpy matter distribution. I assumed
the primeval mass fluctuations produce curvature fluctuations (in the sense of Equation 6) that
diverge only as the logarithm of the length scale. It was known then that an approximation to this
spectrum is naturally produced by inflation, but I was more taken by scale invariance (Harrison
1970, Peebles & Yu 1970, Zel’dovich 1972). I first used EdeS, but the CfA sample led me to work
out the result of going to lower mass density with a cosmological constant to allow flat space
sections (Peebles 1984). I computed only the quadrupole CMB anisotropy in the CDM model, to
compare to the measured upper bound, and estimated the anisotropy in the low-density version
by scaling from CDM. Kofman & Starobinskii (1985) did that better. This low-density version
came to be known as the CDM model.
I liked CDM, but distrusted its early use in analyses of structure formation. I had just made
up CDM to save the gravitational instability picture, which didn’t make it right, and I set out to
show that by inventing other viable models. I gave up (a last example is Peebles 1999) because
the rapidly advancing observations were ruling out my models as fast as I could put them on
astro-ph, and they were agreeing wonderfully well with CDM. Earlier than that, the beautiful
measurements by Mather et al. (1990) and Gush, Halpern & Wishnow (1990) had shown that the
22 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
CMB spectrum is very close to thermal, removing the worry about serious perturbations to the
CMB by explosions. Still earlier, Bond & Efstathiou (1987) were willing to put aside my worries and
compute. Their CMB anisotropy spectrum looks much like modern measurements. Sometimes
it pays to prefer ideas to observations (though observations trumped ideas in the addition of the
previously unpopular cosmological constant).
The CMB temperature anisotropy is part of the network of evidence that now establishes
CDM as a very successful approximation. In the 1990s, there were other models of structure
formation to consider (Peebles & Silk 1990), there were arguments for EdeS and for a lower
density Friedman-Lemaı̂tre solution, and there was the less widely discussed but very relevant
issue of reliability of the extrapolation of general relativity to the enormous scales of cosmology.
The conference Cosmic Velocity Fields gives a fair sample of the arguments for the EdeS mass density
that led to Sandy Faber’s summary remark that “. . . might actually be close to 1” (Bouchet &
Lachièze-Rey 1993, p. 491). But this was work in progress, and it is worth listing the growth of a
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
considerable variety of evidence for lower m and CDM. Sandage (1961) suspected the expansion
time in EdeS may be shorter than stellar evolution ages of the oldest stars, a problem that could
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
be avoided if > 0 with m < 1. This gained weight (Gunn & Tinsley 1975; Turner, Steigman
& Krauss 1984; Chaboyer et al. 1996) with advances in stellar evolution theory and measures of
the cosmic distance scale (leading to Freedman et al. 2001). Sandage (1961) concluded that the
relation between redshifts and apparent magnitudes of distant objects is a promising cosmological
test. He had in mind galaxies. The successful application to supernovae indicates m < 1 and
> 0 (Riess et al. 1998, Perlmutter et al. 1999). Other early evidence for m < 1 includes
the evolution of cluster galaxy velocity dispersions (Bahcall, Fan & Cen 1997) and intracluster
plasma temperatures (Henry 1997), the cluster plasma mass fraction (White et al. 1993), and the
relative velocities of galaxies (Davis & Peebles 1983; Shaya, Peebles & Tully 1995). Assuming
some variant of the CDM model, measures of the maximum extent of correlated structure could
be interpreted as m < 1 with open space curvature (Blumenthal, Dekel & Primack 1988), perhaps
with a power-law tilt of the initial conditions (Vittorio, Matarrese & Lucchin 1988) or CDM
(Silk & Vittorio 1987; Efstathiou, Sutherland & Maddox 1990). This constraint also allows EdeS
with a feature in the initial conditions that happens to mimic scale-invariant initial conditions in a
low-density Universe (Efstathiou, Bond & White 1992), but we are saved from that ugly prospect
by all the evidence that the mass density is less than EdeS. By the end of the decade, the first
peak of the CMB temperature anisotropy spectrum was well identified (figure 2 of Miller et al.
1999). Fitting the peak to CDM with EdeS requires an absurdly long distance scale, whereas the
fit to CDM with low m agrees with the astronomers’ distance scale. These last two constraints
assume variants of the CDM model and help establish this structure formation model and the
relativistic cosmology. There are more tests now, and they are tighter (as reviewed in FTBB), but
these examples make the key point. Observations of a broad variety of phenomena allow cross-
checks of reliability of the measurements and searching tests of the cosmology by looking at the
Universe in many independent ways. This means that if there is a better cosmology than CDM,
it will predict a Universe that looks much like CDM.
6. FUTURE DIRECTIONS
Cosmology has grown far beyond anything I, or I suppose anyone, dreamed a half century ago,
but it remains incomplete. Recall Lifshitz (1946, p. 116): “. . . gravitational instability is not the
source of condensation of matter into separate nebulae.” If early Universe studies prosper to the
point of substantiating inflation or some other theory, it may inform us of Lifshitz’s source of
condensation, but I expect the theory to remain incomplete in Lifshitz’s sense: It will require
hypotheses about still deeper physics. Maybe there is a final theory in Weinberg’s (1992) sense,
and from it an effective theory for cosmology; but if so I imagine deriving that effective theory will
require approximations to be debated by theorists and tested by observers. Thus I suspect it will
never be known whether cosmology is being played with a full deck, with all the physics relevant
for all that could be measured. But this lack of completeness is normal in all natural science.
These days it usually is taken as well demonstrated that we have a complete theory for cosmology
after baryogenesis: CDM with allowance for modest adjustments of initial conditions, maybe
a variable , maybe dark matter that does something a little more interesting than CDM. My
thinking about this has to be colored by my experience: The many wrong turns taken on the path
to CDM, and the many right turns overlooked, do not inspire confidence that there will be no
more unexpected developments. Some may be quite modest: I expect the great range of levels of
nonlinear structures will drive generalizations of HOD (Section 4.1). Some may be more serious:
I see the tendency of galaxies to act as island universes (e.g., Peebles & Nusser 2010; Nair, van den
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
Bergh & Abraham 2011; Peebles 2011) as a challenge to established ideas about galaxy formation,
maybe to be remedied by better methods of analysis of CDM, maybe by adjustments of this
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
theory. But CDM is a good approximation, and it is good tactics for the community to act as if
we had the complete theory. It concentrates research, and if this leads to general agreement that
an adjustment to CDM is required, then the evidence very likely will be right and point to a
better theory.
While awaiting unexpected developments it can be good tactics for some to seek issues that
seem seriously problematic for conventional lines of thought. I have a few favorites (Peebles 2001,
Peebles & Nusser 2010); others are thinking in other directions. When something demonstrably
challenging to established wisdom turns up, the community should hear about it. I like Fred
Hoyle’s advice: Publish and be damned—but keep it short.
DISCLOSURE STATEMENT
The author is not aware of any affiliations, memberships, funding, or financial holdings that might
be perceived as affecting the objectivity of this review.
ACKNOWLEDGMENTS
I am grateful to Marc Davis, Sandy Faber, Jim Fry, Ed Groth, Mike Hauser, Ed Jenkins, Adi
Nusser, Alison Peebles, Bob Pollock, Ravi Sheth, Ken Standing, and David Weinberg for their
help in framing this essay.
LITERATURE CITED
Abbott LF, Sikivie P. 1983. Phys. Lett. 120B:133–36
Abell GO. 1958. Ap. J. Suppl. 3:211–88
Alfvén H. 1966. Worlds-Antiworlds. San Francisco: WH Freeman
Alpher RA. 1948. On the origin and relative abundance of the elements. PhD thesis. The George Washington
Univ., Washington, DC
Alpher RA, Bethe H, Gamow G. 1948. Phys. Rev. 73:803–4
Alpher RA, Herman R. 1948. Nature 162:774–75
Bahcall NA, Fan X, Cen R. 1997. Ap. J. 485:L53–56
Berlind AA, Weinberg DH. 2002. Ap. J. 575:587–616
Bernardi M, Nichol RC, Sheth RK, Miller CJ, Brinkmann J. 2006. Astron. J. 131:1288–317
Blackman RB, Tukey JW. 1959. The Measurement of Power Spectra. New York: Dover
24 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
Gamow G, Critchfield CL. 1949. Theory of Atomic Nucleus and Nuclear Energy-Sources. London, UK: Oxford
Univ. Press
Geller MJ, Peebles PJE. 1973. Ap. J. 184:329–41
Giocoli C, Bartelmann M, Sheth RK, Cacciato M. 2010. MNRAS 408:300–13
Gingerich O. 1994. Phys. Today 47:34–40
Gott JR III, Schramm DN, Tinsley BM, Gunn JE. 1974. Ap. J. 194:543–53
Groth EJ, Peebles PJE. 1977. Ap. J. 217:385–405
Groth EJ, Peebles PJE. 1986a. Ap. J. 310:499–506
Groth EJ, Peebles PJE. 1986b. Ap. J. 310:507–17
Gunn JE, Tinsley BM. 1975. Nature 257:454–57
Gush HP, Halpern M, Wishnow EH. 1990. Phys. Rev. Lett. 65:537–40
Guth A. 1981. Phys. Rev. D 23:347–56
Harrison ER. 1970. Phys. Rev. D 1:2726–30
Hauser MG, Peebles PJE. 1973. Ap. J. 185:757–85
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
Herzberg G. 1945. Molecular Spectra and Molecular Structure. II. Infrared and Raman Spectra of Polyatomic
Molecules. New York: Van Nostrand Reinhold
Hogg DW, Eisenstein DJ, Blanton MR, Bahcall NA, Brinkmann J, et al. 2005. Ap. J. 624:54–58
Hoyle F. 1950. Observatory 70:194–97
Hoyle F. 1981. New Sci. 92:521–24
Hoyle F, Tayler RJ. 1964. Nature 203:1108–10
Hubble EP. 1936. Realm of the Nebulae. New Haven, CT: Yale Univ. Press
Hubble EP, Humason ML. 1931. Ap. J. 74:43–80
Irvine WM. 1961. Local irregularities in a universe satisfying the cosmological principle. PhD thesis. Harvard Univ.,
Cambridge, MA
Jenkins A, Frenk CS, Pearce FR, Thomas PA, Colberg JM, et al. 1998. Ap. J. 499:20–40
Katz L, Mulders GFW. 1942. Ap. J. 95:565–68
Kiang T. 1967. MNRAS 135:1–22
Kiang T, Saslaw WC. 1969. MNRAS 143:129–38
Kirshner RP, Oemler A Jr, Schechter PL. 1978. Astron. J. 83:1549–63
Klein O. 1966. Nature 211:1337–41
Kofman LA, Starobinskii AA. 1985. Sov. Astron. Lett. 11:271–74
Komatsu E, Smith KM, Dunkley J, Bennett CL, Gold B, et al. 2011. Ap. J. Suppl. 192(2):18
Kormendy J, Drory N, Bender R, Cornell ME. 2010. Ap. J. 723:54–80
Landau L, Lifshitz E. 1951. The Classical Theory of Fields. Reading, MA: Addison-Wesley
Larson D, Dunkley J, Hinshaw G, Komatsu E, Nolta MR, et al. 2011. Ap. J. Suppl. 192(2):16
Layzer D. 1956. Astron. J. 61:383–85
Lemaı̂tre G. 1933. Comptes Rendus 196:1085–87
Lifshitz E. 1946. J. Phys. 10:116–19
Limber DN. 1954. Ap. J. 119:655–81
Lynden-Bell D. 2010. Annu. Rev. Astron. Astrophys. 48:1–19
Mach E. 1893. The Science of Mechanics. Transl. TJ McCormack, 1960. La Salle, IL.: The Open Court. 6th ed.
Mandelbrot BB. 1975. Les Objects Fractals. Paris: Flamarrion
Mather JC, Cheng ES, Eplee RE Jr, Isaacman RB, Meyer SS, et al. 1990. Ap. J. 354:L37–40
Matsumoto T, Hayakawa S, Matsuo H, Murikami H, Sato S, et al. 1988. Ap. J. 329:567–71
McKellar A. 1941. Publ. Dom. Astrophys. Obs. 7:251–72
Michie RW. 1969. Kitt Peak Natl. Obs. Preprint
Miller AD, Caldwell R, Devlin MJ, Dorwart WB, Herbig T, et al. 1999. Ap. J. 524:L1–4
Molaro P, Vangioni E. 2009. Mem. Soc. Astron. Ital. 80:749:932
Nair P, van den Bergh S, Abraham RG. 2011. Ap. J. 734:L31:1–5
Neyman J, Scott EL. 1952. Ap. J. 116:144–63
Neyman J, Scott EL, Shane CD. 1956. Berkeley Symp. Math. Stat. Probab. 3:75–111
26 Peebles
AA50CH01-Peebles ARI 16 July 2012 11:15
Sciama DW. 1959. The Unity of the Universe. Garden City, NY: Doubleday
Seldner M, Peebles PJE. 1977. Ap. J. 214:L1–4
Seldner M, Siebers B, Groth EJ, Peebles PJE. 1977. Astron. J. 82:249–314
Shakeshaft JR, Ryle M, Baldwin JE, Elsmore B, Thomson JH. 1955. MmRAS 67:106–54
Shane CD, Wirtanen CA. 1967. Publ. Lick Obs. 22:1–60
Shapley H, Ames A. 1932. Ann. Harv. Coll. Obs. 88:41–75
Shaver PA. 1991. Aust. J. Phys. 44:759–69
Shaya EJ, Peebles PJE, Tully RB. 1995. Ap. J. 454:15–31
Sheth RK, Jain B. 2003. MNRAS 345:529–38
Silk J. 1967. Nature 215:1155–56
Silk J, Vittorio N. 1987. Ap. J. 317:564–75
Soneira RM, Peebles PJE. 1978. Astron. J. 83:845–60
Strominger D, Hollander JM, Seaborg GT. 1958. Rev. Mod. Phys. 30:585–904
Sylos Labini F, Montuori M, Pietronero L. 1998. Phys. Rep. 293:61–226
Annu. Rev. Astron. Astrophys. 2012.50:1-28. Downloaded from www.annualreviews.org
28 Peebles
AA50-FrontMatter ARI 19 July 2012 14:56
Annual Review of
Astronomy and
Astrophysics
P.J.E. Peebles p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 1
Magnetic Fields in Molecular Clouds
Richard M. Crutcher p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p29
The Formation and Early Evolution of Low-Mass Stars
and Brown Dwarfs
Kevin L. Luhman p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p65
Presupernova Evolution of Massive Single and Binary Stars
N. Langer p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 107
Critical Reactions in Contemporary Nuclear Astrophysics
M. Wiescher, F. Käppeler, and K. Langanke p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 165
Planet-Disk Interaction and Orbital Evolution
W. Kley and R.P. Nelson p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 211
Galactic Stellar Populations in the Era of the Sloan Digital Sky Survey
and Other Large Surveys
Željko Ivezić, Timothy C. Beers, and Mario Jurić p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 251
Adaptive Optics for Astronomy
R. Davies and M. Kasper p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 305
Formation of Galaxy Clusters
Andrey V. Kravtsov and Stefano Borgani p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 353
Microlensing Surveys for Exoplanets
B. Scott Gaudi p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 411
Observational Evidence of Active Galactic Nuclei Feedback
A.C. Fabian p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 455
Gaseous Galaxy Halos
M.E. Putman, J.E.G. Peek, and M.R. Joung p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p p 491
vii
AA50-FrontMatter ARI 19 July 2012 14:56
Indexes
Errata
Access provided by 103.216.207.232 on 11/29/23. For personal use only.
An online log of corrections to Annual Review of Astronomy and Astrophysics articles may
be found at http://astro.annualreviews.org/errata.shtml
viii Contents