General Relativity and Cosmology: Unsolved Questions and Future Directions

Article
General Relativity and Cosmology: Unsolved

Questions and Future Directions
Ivan Debono 1,, and George F. Smoot 1,2,3,
1
arXiv:1609.09781v1 [gr-qc] 28 Sep 2016
Paris Centre for Cosmological Physics, APC, AstroParticule et Cosmologie, Universit Paris Diderot,
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cit, 10, rue Alice Domon et Lonie Duquet,
75205 Paris CEDEX 13, France; [email protected]
Physics Department and Lawrence Berkeley National Laboratory, University of California, Berkeley,
94720 CA, USA
Helmut and Anna Pao Sohmen Professor-at-Large, Hong Kong University of Science and Technology,
Clear Water Bay, Kowloon, 999077 Hong Kong, China
Correspondence: [email protected]; Tel.: +33-1-57276991
These authors contributed equally to this work.
Academic Editors: Lorenzo Iorio and Elias C. Vagenas

Received: 21 August 2016; Accepted: 14 September 2016; Published: 28 September 2016
Abstract: For the last 100 years, General Relativity (GR) has taken over the gravitational theory
mantle held by Newtonian Gravity for the previous 200 years. This article reviews the status of GR
in terms of its self-consistency, completeness, and the evidence provided by observations, which
have allowed GR to remain the champion of gravitational theories against several other classes of
competing theories. We pay particular attention to the role of GR and gravity in cosmology, one of
the areas in which one gravity dominates and new phenomena and effects challenge the orthodoxy.
We also review other areas where there are likely conflicts pointing to the need to replace or revise
GR to represent correctly observations and consistent theoretical framework. Observations have long
been key both to the theoretical liveliness and viability of GR. We conclude with a discussion of the
likely developments over the next 100 years.
Keywords: General Relativity; gravitation; cosmology; Concordance Model; dark energy; dark
matter; inflation; large-scale structure
1. Perspective
Scientists have been fascinated by General Relativity ever since it was developed. It has been
described as poetic, beautiful, elegant, and, at times, as impossible to understand.
General Relativity is often described as a simple theory. It is hard to define simplicity in science.
One can always construct an entire theory encapsulated in one equation. Richard Feynman famously
~ = 0,
demonstrated this in a thought experiment where he rewrote all the laws of physics as U
~
where each element of U contained the hidden structure [1]. His point was that simplicity does not
automatically bring truth.
An examination of the mathematical structure of General Relativity gives us a more sober
definition of simplicity. Under certain assumptions about the structure of physical theories, and of
the properties of the gravitational field, General Relativity is the only theory that describes gravity.
Alternative theories introduce additional interactions and fields.
General Relativity is also unique among theories of fundamental interactions in the Standard
Model. Like electromagnetism, but unlike the strong and weak interactions, its domain of validity
covers the entire range of length scales from zero to infinity. However, unlike the other forces, gravity
as described by General Relativity acts on all particles. This implies that the theory does not fail
Universe 2016, 2, 23; doi:10.3390/universe2040023
www.mdpi.com/journal/universe
Universe 2016, 2, 23
2 of 82
below the Planck scale. All gravitational phenomena, from infinitesimal scales to distances beyond
the observable universe, may be modelled by General Relativity. We may therefore formulate a
mathematically rigorous description of General Relativity: it is the most complete theory of gravity
ever developed.
All gravitational phenomena that have ever been observed can be modelled by General Relativity.
It describes everything from falling apples, to the orbit of planets, the bending of light, the dynamics
of galaxy clusters, and even black holes and gravitational waves. The domain of validity of the theory
covers a wide range of energy levels and scales. That is why it has survived so long, and that is why it
survives today, one hundred years after it was formulated, in an age in which the amount of data and
knowledge increases by orders of magnitude every few years.
Lookback time
Present
z=0
13.8 billion years

geological
record
5 Gy
z=0.5
13.4 Gy
z=12
13.5 Gy
z=15
5 billion years
quasars
1.8 billion years
population III stars
200 million years
neutral hydrogen
21 cm
z=100
z=1100
z=4108
13.8 Gy
z=6109
5 million years
CMB
380,000 years
BBN
3 minutes
neutrinos
1 second
Time since Big Bang
Figure 1. How we observe the universe. The lookback time is the difference between the age of the
universe now, and the age of the universe when photons from an object were emitted. The more
distant an object, the farther in its past we are observing its light. This distance in both space and
time is expressed by the cosmological redshift z. We obtain most of our astrophysical information
from the surface of our past light cone, because it is carried by photons. The only information from
within the cone come from local experiments and observations, such as geological records. The green
dotted line is the world-line of the atoms and nuclei providing the material for our geological data.
Local experiments are carried out along this bundle of world-lines. They provide a useful test of physical
constants. One example is the observation of the Oklo phenomenon [2]. The earliest information
we have collected so far comes from the cosmic microwave background (CMB). Earlier than the
CMB time-like slice is the cosmic neutrino background. We observe Big Bang nucleosynthesis (BBN)
indirectly, through observations of the abundances of chemical elements.
Why, then, are we still testing General Relativity? Why do we still develop, discuss, test and
fine-tune alternative theories? Because there are some very fundamental open questions in physics,
particularly in cosmology. Moreover, the big questions in cosmology happen to be the ones that are not
answered by General Relativity: the accelerated expansion of the universe, the presence of a mysterious
form of matter which cannot be observed directly, and the initial conditions in the early universe.
3 of 82
The theoretical completeness described above is both a necessary and aesthetic feature of a
fundamental theory. However, it creates experimental difficulties, for it compels us to test the theory at
extreme scales, where experimental errors may be large enough to allow several alternative theories.
At extremely small scales below the Planck length, classical mechanics should break down.
This compels us to question whether General Relativity is still accurate at these scales, whether it
needs to be modified, and whether a quantum description of gravity can be formulated. At the other
end of the scale, at cosmological distances, we may question whether General Relativity is valid, given
that the universe cannot be modelled sufficiently accurately by General Relativity without invoking
either a cosmological constant, or some additional, unknown component of the universe. Finally, we
may question the accuracy of our solutions to the equations of General Relativity, which depend on
some approximation scheme. These approximations provided analytical solutions which enabled most
of the early progress in General Relativistic cosmology and astrophysics. However, one century after
the formulation of the theory, we now have a flood of data from increasingly accurate observations
(as shown in Figure 1), coupled with computing power which was hitherto unheard of. Tests of the
higher-order effects predicted by General Relativity and some of its competitors are now within reach.
The purpose of this review is to examine the motivation for the development of alternative
theories throughout the history of GR, to give an overview of the state of the art in General
Relativistic cosmology, and to look ahead. In the next few decades, some of the open questions
in cosmology may well be answered by a new generation of experiments, and GR may be challenged
by alternative theories.
2. A Brief History
Let us start this review by breaking our own rule about unscientific adjectives. General Relativity
is a beautiful theory of gravity. It has not only thrilled us, but has survived 100 years of challenges, both
by experimental tests and by alternative theories. The beauty of the theory was clear at the beginning,
but the initial focus was on whether it was right. When General Relativity provided an explanation for
the 43 seconds of arc per century discrepancy in the advance of the perihelion of Mercury [3], it got
the attention of the scientific community. However, it was the prediction and the observation of the
bending of light by the Sun [4] that confirmed GRs place as the new reigning theory of gravity [5].
The setting at the Royal Society under the portrait of Newton for the report of the eclipse light
bending observations led by Arthur Eddington, and reported by the great writer Aldous Huxley,
was perfect to describe to the world the ascendancy of a new theory replacing Newtons gravity
(see, e.g., [6]). From this point onwards, the scientific community started to take General Relativity
seriously, and theorists worked hard to understand this new theory, beguilingly simple but hard to
apply, and to advance its predictions.
Shortly after its publication, GR quickly became the framework for astrophysics, and for the
Standard or Concordance Model of cosmology. However, it was still challenged by alternative theories.
Initially, the alternatives were motivated by theoretical considerations. This early period led to a fuller
understanding of GR and its predictions. Some of the predictions, such as black holes and gravitational
waves, divided the scientific community. Did they exist as physical objects, or just as mathematical
artifacts of the theory?
By the time GR turned 50, the model of cosmology had been established, GR had been tested, and
things had started to stagnate. However, advances in observations led to new discoveries, which in
turn led to renewed challenges.
First came the missing mass in the universe. Could GR be modified to account for it? Then came
the theories about the very early universe, and the behaviour of the quantum-scale, tiny initial universe.
Finally, twenty years ago, came the confirmation of cosmic acceleration. This had a twofold effect.
On one hand, it spurred the development of a whole range of alternative theories of gravity. On the
other hand, it confirmed GR like never before, for General Relativity, with a cosmological constant,
can account for the observations perfectly.
4 of 82
In 2015, on the 100th birthday of General Relativity, gravitational waves were observed for the
first time. This had been the last major untested prediction of General Relativity. It was a remarkable
achievement, and in many ways it heralds a new age of astrophysical observations. The experimental
capabilities and the computing power have finally caught up with the theory. Cosmology and
astrophysics have now entered the era of Big Data, and much of the theoretical effort is now driven
by data. However, the foundation for almost the entire scientific endeavour is still this theory of
chronogeometrodynamics, developed 100 years ago when todays instruments and computers were
still a distant dream.
From Aristotle to Einstein
General Relativity is the basis for the Standard Model of physical cosmology, and here we shall
discuss the development of General Relativity (GR). The history of cosmology and GR are intertwined.
We shall discuss why the theory has been so successful, and the criteria that must be satisfied by any
alternative theory of physics, and by cosmological models.
Cosmology, in its broadest definition, is the study of the cosmos. It aims to provide an accurate
description of the universe. Throughout much of the history of science, the development of cosmology
was hampered by the lack of a universal physical theory. Observational tools were extremely
limited, and there was no mathematical formulation for physical laws. The cosmos was described in
metaphysical, rather than physical terms.
Discussions on the history of physics often refer to Karl Poppers concept of Falsifizierbarkeit
(falsifiability) [7]. In this formalism, scientific discovery proceeds by successive falsifications of theories.
A falsifiable theory that covers observations, and that has not yet been proven false can be regarded
as provisionally acceptable. Yet we know that in reality it is not quite as straightforward. A theory
that is considered to be correct acquires this status by accumulation of evidence rather than by a
single falsification of a previous theory [8]. This is especially true in cosmology, where the selection of
theoretical models often depends on the outcome of statistical calculations.
The scientific revolution which brought about the development of a precise mathematical language
for physical theories heralded the scientific age of cosmology. Physical laws, tested here on Earth and
later in the Solar System, could be applied to the entire universe, and could thus provide a precise
physical description of the cosmos. Modern cosmology is based upon this epistemological framework.
Cosmology depends upon a fundamental premise. As a science, it must deal strictly with what can
be observed, but the observable universe forms only a fraction of the whole cosmos. One is forced
to make the fundamental but unverifiable assumption that the portion of the universe which can be
observed is representative of the whole, and that the laws of physics are the same throughout the
whole universe [9]. Once we make this assumption, we can construct a model of the universe based on
a description of its observable part.
Any cosmological model which assumes the universality of physical laws must be based upon
some physical theory. Since cosmology aims to describe the universe on the largest possible scales,
it must be based upon an application long-range physical interactions. Since the theory of gravitation
is the physical theory at the basis of standard cosmology, and is also at the centre of the big questions
facing modern cosmology, we shall give an overview of the development of theories of gravitation.
The development of physical theories of gravity was far from smooth, nor did it always conform
to Poppers scheme. Before the logical tools (mathematics) for the phenomenological description
(physics) were invented, progress was rather haphazard.
According to Poppers scheme, this development should be driven by the search for ever more
general principles. Yet Aristotelian theory, to take one example, considered itself to be general
enoughits claimed region of validity was the entire universe, except that rising smoke, floating
feathers, falling apples and orbiting celestial spheres each had their own rules.
5 of 82
The real revolution came when it was realised that the behaviour of all bodies could be described
by a single rulea universal theory of gravitation. This theory is a description of the long range forces
that electrically neutral bodies exert on one another because of their matter content.
Whether they choose to or not, scientists will always stand on the shoulders of giants. No theory is
invented in a scientific vacuum. This goes all the way back to the cosmology of the Euro-Mediterranean
Ancient World, codified in the Aristotelian teachings of the 4th century B.C. This Hellenic natural
philosophy provided qualitative rather than quantitative descriptions for what we would call
today the free parameters of the theory [10]. It stands to reasonthe instruments had not yet
been invented that could test the theory of gravity to within numerical accuracy. Without accurate
timekeeping instruments, processes could at best be described as slower than or faster than.
However, instruments to measure the movement of the celestial bodies, such as sundials, quadrants
and astrolabes, were invented and improved upon, and measurements were carried out [11].
Astronomy flourished.
There is a certain logic to the development of physical theories from the Ancient World, to the
Middle Ages, and right up to the Renaissance [1214]. The basic tenet of the physics of Aristotle
is that actions follow logically from causes. He distinguished between natural and violent motion.
Natural motion implies falling at a speed proportional to the weight of the object and inversely
proportional to the density of the medium. Violent motion happens whenever there is a force acting
on an object, and the speed of the object is proportional to this force. Strato of Lampsacus replaced
Aristotles explanation of unnatural motion with one that is very close to the modern notion of inertia.
He identified natural motion as a form of acceleration, and demonstrated experimentally that falling
bodies accelerate. In the 14th century, Jean Buridan came up with the notion of impetus, where the
initial force imparts motion to the object, which gradually diminishes as gravity and air resistance
act against this initial force. Concurrently, Nicole Oresme was using a crude early form of graph to
describe motion, and unwittingly showing the complicated notions of differentiation and integration
in pictorial form[15,16].
The cosmological observations, limited to the innermost five planets of the Solar System (Mercury,
Venus, Mars, Jupiter, and Saturn) and the sphere of stars, seemed to confirm the Aristotelian-Ptolemaic
theory. Celestial bodies moved in regular patterns made up of repeating circles. Small discrepancies
were explained by circles within circles.
The fact that the theories were based on these regular patterns is no accident. Patterns are the
keyword in all of physics. Human beings are wired to recognise patterns. We can only build theories
because we recognise patterns in the universe. This characteristic of valid theories has been called
sloppiness. The patterns fall within some hyper-ribbon of stability in the theory [17].
The revolution in physics came with the development of mathematical, quantitative, models to
describe physical reality. Starting in the 1580 Galileo carried out a series of observations in which
he subjected kinematics to rigorous experiment, and showed that naturally-falling objects really do
accelerate. Crucially, he showed that the composition of the body has no effect whatsoever on this
acceleration. He also realised that for violent motion, the speed is constant in the absence of friction.
Galileo also took rigorous observations of astronomical objects. In 1610 he made the first observation
of Jupiters satellites, and the first observation of the phases of Venus, which is impossible according to
the Ptolemaic geocentric model. His observations were important in putting to rest the Aristotelian
theory of perfect and unchanging heavens.
By the time Newton came along, telescopes had been invented. Galileo had observed moons
orbiting the Solar System planets, and hundreds of stars invisible to the naked eye. His 1610 treatise,
aptly called Siderus Nuncius (Starry Message, or Astronomical Report in modern language) [18],
was the first scientific work based on observations through a telescope. Mechanical clocks had been
invented. The sphere of observed data had expanded [19]. Calculus provided the tool to make sense
of this new flood of data. Thus, physicists of Newtons generation found a very different scientific
environment than the one in which Galileo had started off.
6 of 82
In 1687, Isaac Newton published in his Mathematical Principles of Natural Philosophy, known
by its abbreviated Latin title as Principia [20]. This was a significant milestone in physics. Newtons
model of gravitation was, in his own words, a universal law. It applied to all bodies in the
universe, whether it was cannonballs on Earth, or planets orbiting the Sun. For more than two
centuries, Newtons theory, was the standard physical description of gravity. There was no other
attempt to find a different theory for the gravitational force, although the intervening years between
Newtonian gravity and Relativity produced some important physical concepts such as de Maupertuiss
Principle of Least Action [21], further developed by Euler [22], Lagrange [23] and Hamilton [24,25].
The path of each particle is assigned a number called an action, which is the integral of the Lagrangian.
In classical mechanics, the action principle is equivalent to Newtons Laws. Lagrangian field theory
is an important cornerstone of modern physics. The Lagrangian of any physical interaction, when
subjected to an action principle, give us field equations and conservation laws for the theory. It is an
expression of the symmetries in physical laws.
Newtonian gravity was the great success story of nineteenth century physics, the golden age of
mathematical astronomy. It allowed astronomers to calculate the position of planets and asteroids
with ever greater precision, and to confirm their calculations by observation. Thus the size of the
known universe grew. Evidence started to accumulate suggesting that there might be other galaxies
in the universe besides our own. In 1845, the planet Neptune was discovered, after Urbain le Verrier
suggested pointing telescopes in a region of the Solar System which he predicted by Newtonian
calculations [26,27]. The search was motivated in the first place by an anomaly in the orbit of Uranus
which could not be otherwise explained using Newtonian theory [28]. The discovery of Neptune
showed that Newtonian theory was valid even in the very farthest limits of the Solar System.
There was another anomaly which could not be explainedthe excessive perihelion precession of
Mercury by 43 arcseconds per century, confirmed by le Verrier himself. Urbain le Verrier thus holds
the distinction of being one of the few experimentalists to have confirmed Newtons theory and then
disproved it. Astronomers attempted to explain this perihelion anomaly using Newtonian mechanics,
which led them to speculate on the existence of Vulcan, a hypothetical planet whose orbit was even
closer to the Sun [29].
The first doubts on Newtonian theory began to take shape just at the time when theorists were
examining the full implications of the theory for complex, multi-body dynamical systems such as the
Solar System. In 1890, Henri Poincar published his magnum opus on the three-body problem [30],
a masterpiece of celestial mechanics. At the time, Poincar was working on another open question in
physics: the aether. This led him to formulate a theory which was very close to Special Relativity [31],
but which did not quite fit with Maxwells electromagnetism [32], and was ultimately flawed.
By the end of the 19th century, the necessary mathematical tools were in place which would enable
the development of Special and then General Relativity. There is an intimate connection between
physics and the language of mathematics which is often overlooked. The former, especially in modern
times, depends on the latter. Could Aristotle have developed General Relativity? No. Because he had
not the mathematical language. Equations and mathematical formulations are relatively recent in the
history of physics. Even Newton, for all his fame as a mathematical genius, never wrote the equation
F = GMm/r2 . He wrote a series of statements implying this law in (Latin) words: Gravitatem, qu
Planetam unumquemque respicit, ese reciproc ut quadratum distanti locorum ab ipsius centro, and so on.
It is hard to imagine how human beings could manipulate tensors and solve the field equations of
Relativity in anything but numbers and symbols. Theories and physics do not happen in a cultural and
scientific vacuum. They are human creations, and they depend intimately on tools for the transmission
and communication of human knowledge.
The physical theory of gravitythe laws that govern gravitational interactionsremained
unchanged until Einsteins time. In 1905, Einstein published his Theory of Special Relativity
(SR) [33]. Soon after, he turned to the problem of including gravitation within four-dimensional
spacetime [3437].
7 of 82
Newtons formulation of the gravitational laws is expressed by the equations:
d2 xi
= i,
dt2
x
4 = 4G ,
(1)
(2)
where is the gravitational potential, G is the universal gravitational constant, is the mass density,
and 4 = 2 is the Laplace operator. These equations cannot be incorporated into Special Relativity
as they stand. The equation of motion (1) for a particle is in three-dimensional form, so it must be
modified into a four-dimensional vector equation for d2 x / d 2 . Similarly, the field Equation (2) is
not Lorentz-invariant, since the three-dimensional Laplacian operator instead of the four-dimensional
dAlembertian = means that the gravitational potential responds instantaneously to changes
in the density at arbitrarily large distances. The conclusion is that Newtonian gravitational fields
propagate with infinite velocity. In other words, instantaneous action in Newtonian theory implies
action at a distance when reconsidered in the light of Special Relativity. This violates one of the
postulates of SR. How do we reconcile gravity and Special Relativity?
3. The Development of General Relativity
3.1. From Special to General Relativity
The simplest relativistic generalisation of Newtonian gravity is obtained by representing the
gravitational field by a scalar . Since matter is described in Relativity by the stress-energy tensor T ,
the only scalar with dimensions of mass density (which corresponds to ) is T . A consistent scalar
relativistic theory of gravity would thus have the field equation
= 4GT .
(3)
However, when the equation of motion from this theory are applied to a static, spherically
symmetric field , such as that of the sun, acting on an orbiting planet, they would result in a
negative precession, or retardation of the perihelion. Experimental evidence since the time of Le Verrier
and his observation of the orbit of Mercury [38] clearly shows that planets experience a prograde
precession of the perihelion. Moreover, in the limit of a zero rest-mass particle, such as a photon,
the equations of motion show that the particle experiences no geodesic deviation. The existence
of an energy-momentum tensor due to an electromagnetic field would also be impossible, since
( Telectromagnetic ) = 0. The theory therefore allows neither gravitational redshift, nor deviation of light
by matter, both of which are clearly observable phenomena [39]. Another route to generalisation could
be to represent the gravitational field by a vector field , analogous to electromagnetism. Following
through with this strategy, the Coulomb law in this theory gives a repulsion between two massive
particles, which clearly contradicts observations. The theory also predicts that gravitational waves
should carry negative energy, and, like the scalar theory, predicts no deviation of light. Like the scalar
theory, then, the vector theory must be discarded.
What about a flat-space tensor theory? The gravitational field in this theory is described by a
symmetric tensor h = h . The choice of the Lagrangian in this theory is dictated by the requirement
that h be a Lorentz-covariant, massless, spin-two field.
In the 1930s, Wolfgang Pauli and Markus Fierz [40] were the first to write down this Lagrangian
and investigate the resulting theory. The predictions of the theory for deviation of light agree with
those of General Relativity, and are consistent with observations. Since the field equations and gauge
properties are identical to those of Einseins linearised theory, the predictions for the properties of
gravitational waves, including positive energy, agree with those obtained using the linearised theory
in General Relativity. However, the theory differs from General Relativity in its predicted value for the
8 of 82
perihelion precession, which is 43 of that given by GR. This disagrees with the value obtained from
observations of Mercurys orbit.
The theory has an even worse deficiency. If two gravitating bodies (that is, not test particles) are
considered, and the field equations are applied to them, then the theory predicts that gravitating bodies
cannot be affected by gravity, since they all move along straight lines in a global Lorentz reference
frame. This holds for bodies made of arbitrary stress-energy, and since all bodies gravitate, then one
must conclude that no body can be accelerated by gravity, which is a obvious self-inconsistency in
the theory.
The only way in which a consistent theory of gravity can be constructed within Special
Relativity is to consider the geometry of spacetime as the gravitational field itself. In other words,
all matter moves in an effective Riemann space of metric g + h , where is the
Minkowski metric. The requirement of consistency leads us to universal coupling, which implies the
Equivalence Principle.
The existence of curved spacetime can be deduced from purely physical arguments. In 1911,
before he had fully developed General Relativity, Einstein [34] showed that a photon must be
affected by a gravitational field, using conservation of energy applied to Newtonian gravitation
theory. Schild [4143] showed by a simple thought experiment, formulated within Special Relativity,
that a consistent theory of gravity cannot be constructed within this framework. His argument is
based upon a gravitational redshift experiment carried out in the field of the Earth, using a global
Lorentz frame tied to the Earths centre. Successive pulses of light rising to the same height should
experience a redshift, and therefore the pulse rate at the top should be slower than that at the bottom.
But light rays are drawn at 45 degrees in Minkowski spacetime diagrams, so that top and bottom
time intervals are equal, which is impossible if redshift occurs. Hence the spacetime must be curved.
One therefore concludes that in the presence of gravity, Special Relativity cannot be valid over any
sufficiently extended region.
General Relativity may be understood as a generalisation of Special Relativity over extended
regions. Since Special Relativity can comfortably be described using tensor calculus, it was only
natural to extend the flat Minkowski spacetime of Special Relativity to the curved spacetime of General
Relativity. This was a physical application of Riemannian geometry [44,45], which had been developed
in the second half of the 19th century. The idea of tensor calculus on curved manifolds was already
mathematically well-established. Einsteins innovation lay in identifying the Einstein tensor, itself
related to the Riemann curvature tensor, as the gravitational field in the theory.
Einstein had been working on the problem for some years, starting in 1907. He arrived at the final,
correct form in 1915 [46,47]. He was well-aware of the significance of his publication, and he gave it
the succinct title of The Field Equations of Gravitation (Feldgleichungen der Gravitation). The correct
field equations for the theory contained in this publication served as the starting point or subsequent
derivations.
3.2. The Formalism of General Relativity
General Relativity is based on two independent but mutually supporting postulates.
The first postulate is sometimes referred to collectively as the Einstein Equivalence Principle:
The Strong Equivalence Principle: The laws of physics take the same form in a freely-falling reference
frame as in Special Relativity
The Weak Equivalence Principle: An observer in freefall should experience no gravitational field.
That is to say, an observer cannot determine from a local experiment whether the his laboratory is
being accelerated by a rocket of static at the surface of a gravitating body. Gravity is erased up to
tidal forces, which are determined by the size of the laboratory and its distance to the centre of the
gravitational attraction.
The Equivalence Principle allows us to construct the metric and the equation of motion by
transforming from a freely-falling to an accelerating frame. It can be mathematically expressed by
9 of 82
the assuming that all matter fields are minimally coupled to a single metric tensor g . The distance
between two points in 4-dimensional spacetime, called events, is:
ds2 = c2 d 2 = g dx dx .
(4)
Throughout the text, we follow the Einstein summation convention for repeated indices, so that
n
ci xi = ci xi for i = 1, . . . , n. Greek indices are used for space and time components, while Latin
i =1
indices are spatial ones only. We use the following metric signature: ( + ++).
The metric defines lengths and times measured by laboratory rods and clocks. This metric implies
that the action for any matter field is of the form
Smatter [, g ] ,
(5)
which gives us three important results. First, it implies the universality of freefall. Second, it implies
that all non-gravitational constants are spacetime independent. Third, it implies that the laws of physics
are isotropic. This equation defines how matter behaves in a given curved geometry, how light rays
propagate, how stars, planets and galaxies move, and gives us verifiable observational consequences.
The second postulate is related to the dynamics of the gravitational interaction. This is assumed
to be governed by the Einstein-Hilbert action:
Sgravity =
c3
16G
d4 x
g R
(6)
is a massless spin-2 field called the Einstein metric. General Relativity identifies the Einstein
where g
. This implements the Strong Equivalence Principle.
metric with the physical metric, that is: g = g
The Einstein-Hilbert action defines the dynamics of gravity itself. Relativity is thus a geometrical
approach to fundamental interactions. These are realised though continuous classical fields which are
inseparably connected to the geometrical structures of spacetime, such as the metric, affine connection,
and curvature.
The General Relativistic equation of motion is simply parallel transport on curved spacetime.
It is given by
d2 x
dx dx
+
(7)
d d = 0 ,
d 2
where x is some set of coordinates for a point in spacetime. are the components of the affine
connection (or metric connection). The fundamental theorem of Riemannian geometry states that the
affine connection can be expressed entirely in terms of the metric:
1
g ( g, + g, g, ) ,
2
(8)
where the comma denotes a derivative, i.e., g, = x

.
We need to construct invariant quantities in GR (quantities that are the same for all observers).
To achieve this, we need to contract the covariant A and contravariant A components of a vector or
tensor A by using the metric to raise or lower indices: A = g A . Thus the equation of motion (7)
can be made covariant by recasting it as the covariant derivative of the 4-velocity U = (c, v):
D U
= 0,
d
(9)
D A = dA + A dx .
(10)
where the covariant derivative is defined as
10 of 82
The quantity is the Lorentz factor:

=
1
1 v2 /c2
(11)
The transformation from SR to GR is then carried out by mapping the Minkowski metric to a
general metric: g and by mapping D.
In GR, freely-falling bodies travel along a geodesic. Geometrically, this is the shortest distance
between two points in spacetime. The path length along a geodesic is given by
S=
( g dx dx )1/2 .
(12)
In cosmology, it is essential for us to be able to describe spacetime which is not empty. In the
presence of a perfect fluid (an inviscid fluid with density and isotropic pressure p), the energy and
momentum of spacetime is described by the energy-momentum tensor (or stress-energy tensor)

p
T = + 2 U U pg .
c
(13)
Classical energy and momentum conservation are generalized in GR as the four conservation laws
D T = 0 .
(14)
In other words, the stress-energy tensor has a vanishing covariant divergence. In the absence of a
component possessing pressure or density, or both, the energy-momentum tensor is zero.
The central notion in General Relativity is that gravitation can be described by a metric.
The Einstein equations give us the relation between the metric and the matter and energy in
the universe:
8G
(15)
G = 4 T .
c
The left-hand side of this equation is a function of the metric: G is the Einstein tensor, defined as:
1
G = R g R ,
2
(16)
where R is the Ricci tensor, which depends on the metric and its derivatives, and the Ricci scalar R is
the contraction of the Ricci tensor (R = g R ). The right-hand side of Equation (15) is a function of
the energy: G is Newtons constant, and T is the energy-momentum tensor.
Einsteins Relativity has three main distinguishing characteristics:
it agrees with experiment

it describes gravity entirely in terms of geometry
it is free of any prior geometry
These characteristics are lacking in most of the other theories [48,49]. Apart from the issue of
agreement with experiment, Einsteins theory is unique in its physical simplicity.
Every other theory introduces auxiliary gravitational fields, or involves prior geometry.
Prior geometry is any aspect of the geometry of spacetime which is fixed immutably, that is, it cannot
be changed by changing the distribution of gravitating sources.
A rigorous mathematical definition of the unique simplicity of General Relativity is given by
Lovelocks theorem [5052]. This is a generalisation of an earlier theorem by lie Cartan [53], and may
be formulated as follows:
In 4 spacetime dimensions, the only divergence-free symmetric rank-2 tensor constructed solely
from the metric g and its derivatives up to second differential order, and preserving diffeomorphism
invariance, is the Einstein tensor plus a cosmological term.
11 of 82
In simple terms, the theorem states that GR emerges as the unique theory of gravity if the
conditions of the theorem are followed. In fact, Lovelocks theorem provides a useful scheme for
classifying alternatives to General Relativity.
Einstein described both the demand for no prior geometry and for a geometric,
coordinate-independent formulation of physics by the single phrase general covariance, but the
two concepts are not quite the same.
While many physical theories can be formulated in a generally covariant way, General Relativity is
actually based on the no prior geometry demand. This distinction was not always made, especially in
the first decades after Einsteins publications [54,55]. Erich Kretschmanns famous objection in 1917 [56]
concerned this point, since he regarded general covariance merely as formal feature that any theory
could have, not as a special feature belonging to GR.
3.3. Newtonian Nostalgia: The First Wave of Alternative Theories
Newtons theories had predicted observations of Solar System objects, comets and asteroids,
with astounding precision. Why should they be tampered with? The first wave of alternative
theories were driven more by theoretical considerations than by observations. Equations (1) and (2)
can be generalized so that they are consistent with the postulates of Special and General Relativity.
Several generalisations of this kind were attempted in the first few decades following the development
of GR, motivated by lingering resistance to any deviation from Newtonian gravity.
One early theory, involving prior geometry, was formulated by Nordstrm in 1913 [57]. In this
theory, the physical metric of spacetime g is generated by a background flat spacetime metric , and by
a scalar gravitational field . Stress-energy generates :
, = 4 T
(17)
g = 2 .
(18)
and g is constructed from and :

Prior geometry cannot be removed by rewriting Nordstrms equations in a form devoid of
and [58]. Mass only influences one degree of freedom in the spacetime geometry, while the
other degrees of freedom are fixed a priori. This prior geometry, if it existed, could be detected by
physical experiments.
In the 1920s, Alfred North Whitehead [59] formulated a two-tensor theory of gravity in which
the prior geometry is quite different from later theories such as Nis [48]. Whiteheads theory is
remarkable in that it agrees with Einsteins in its predictions for the four standard tests (bending of
light, gravitational redshift, perihelion shift, and time delay). It was accepted as a viable alternative
for Einsteins theory until Clifford Martin Will [60] showed that it predicts velocity-independent
anisotropies in the Cavendish constant (the gravitational constant G in Newtonian theory). This would
produce time-dependent Earth tides which are clearly contradicted by everyday observations.
Any valid theory of gravity must not only agree with relativistic experiments, but also with past
experiments in the Newtonian regime.
One theory which disagrees violently with non-relativistic experiments is due to George David
Birkhoff [61]. It was developed in the 1940s, and it predicts the same redshift, perihelion shift, deflection
and time-delay as General Relativity, but it requires that the pressure inside gravitating bodies should
be equal to the total density of mass-energy (p = ). This means that sound waves travel with the
speed of light. This clearly contradicts everyday experiments.
Most of the early alternative theories were abandoned either because they were contradicted by
observations, or because of internal inconsistencies in the theories themselves. One notable exception is
Dicke-Brans-Jordan theory, sometimes called Brans-Dicke, or Jordan-Fierz-Brans-Dicke theory [62,63],
developed in the 1960s by Robert H. Dicke and Carl H. Brans following earlier work by Pascual Jordan
12 of 82
and Markus Fierz. The different names arise from the fact that the theory is a special case of Jordans,
with = 1. An alternative mathematical representation of the theory is given by [64].
This theory introduced auxiliary gravitational fields. Brans and Dicke took the equivalence
principle as the starting point of their theory, and thus they describe gravity in terms of spacetime
curvature, but their gravitational field, unlike Einsteins, is a scalar-tensor combination. In this way it
overcomes the difficulties associated with tensor or scalar-only theories mentioned earlier. The trace of
the energy-momentum tensor ( TM ) (representing matter) and a coupling constant generate the
long-range scalar field via the equation
2 = 4( TM ) .
(19)
The scalar field fixes the value of G, which is therefore not a constant, but simply the
coupling strength of matter to gravity. The gravitational field equations relate the curvature to
the energy-momentum tensors of the scalar field and matter:
R 21 g R =

8
( TM ) + ( T ) ,
4
c
(20)
where ( TM ) is the energy-momentum tensor of matter and ( T ) is the energy-momentum tensor

of the scalar field . For historical reasons, it is usual to write the coupling constant as
=
2
,
3 + 2
(21)
where is the dimensionless Dicke coupling constant. In the limit , we have 0, so is not
affected by the matter distribution, and can be set to a constant = 1/G. Hence Dicke-Brans-Jordan
theory reduces to Einsteins theory in the limit .
The equivalence principle is satisfied in this theory since the special-relativistic laws are valid
in the local Lorentz frames of the metric g of spacetime. The scalar field does not exert any direct
influence on matter. It only enters the field equations that determine the geometry of spacetime. On a
conceptual level, Brans-Dicke theory can be seen as more fully Machian than Einsteins theory. Einstein
himself attempted to incorporate Machs Principle into his theory, but in Einsteins General Relativity,
the inertial mass of an object will always be independent of the mass distribution in the universe.
In Brans-Dicke theory, the long-range scalar field is an indirectly coupling field, so it does not directly
influence matter, but the Einstein tensor is determined partly by the energy-momentum tensor, and
partly by the long-range scalar field.
Dicke-Brans-Jordan theory is self-consistent and complete, but experimental evidence based on
Solar System tests, shows that 600 [65], as a conservative estimate. Some calculations raise this
limit even higher, with & 104 [66]. The Cassini mission set a comparable limit of > 40, 000 [67].
Recent cosmological data from the Planck probe show that 890 [68,69]. This is consistent with
the Solar System bounds. Future cosmological experiments and data from large-scale structure could
provide even better constraints [70].
Brans-Dicke theory is a special case of general scalar-tensor theories with () = constant, where
is a value depending on the cosmological epoch. In these theories, the function () could be such
that the theory is very different from GR in the early universe or in future epochs, but very close to GR
in the present. In fact, it has been shown that GR is a natural attractor for such scalar-tensor theories,
since cosmological evolution naturally drives the fields towards large values of [71,72].
3.4. Self-Consistency, Completeness, and Agreement with Experiment
Any viable theory must satisfy three fundamental criteria: self-consistency, completeness,
and agreement with past experiment.
13 of 82
To be self-consistent, a theory must not contain any internal contradictions. The spin-two field
theory of gravity [40] is equivalent to linearised General Relativity but it is internally inconsistent
since it predicts that gravitating bodies should have their motion unaffected by gravity. When one
tries to remedy this inconsistency, the resulting theory is nothing but General Relativity. Another
self-inconsistent theory is due to Paul Kustaainheimo [73,74]. It predicts zero gravitational redshift
when the wave version of light (Maxwell theory) is used, and nonzero redshift when the particle
version (photon) is used.
To be complete, a theory must be able to analyse the outcome of any experiment. This means
that it must be compatible with other physical theories which describe any other forces that are
present in experiments. This can only be achieved if the theory is derived from first principles, since
the theoretical postulates must be as general as possible if the theory is to cover the widest range
of phenomena.
A viable theory must agree with past experiment, which includes experiments in the Newtonian
regime, and the standard tests of General Relativity. Its results must agree with those obtained from
Newtonian theory in the weak field limit, and with GR in relativistic situations. It also means that the
theory must agree with cosmological observations.
The experimental criterion also works the other way. Any alternative to General Relativity that
claims to have a smaller set of limiting cases must be experimentally distinguishable, perhaps by
future experiment. At some point, the divergence between GR and other theories must manifest itself
physically, in the form of predictions which can be verified by experiment. This is perhaps the greatest
challenge of current alternatives to GR.
3.5. Metric Theories and Quantum Gravity
Most theories of gravity incorporate two principles: spacetime possesses a metric; and that metric
satisfies the equivalence principle. Such theories are called metric theories. There are some exceptions.
Soon after the publication of the theory of General Relativity, it became apparent that its
formulation is incompatible with a Quantum Mechanical description of the gravitational field. It was
Einstein himself who pointed out that quantum effects must lead to a modification of General
Relativity [75]. Back then, the first successful applications of Quantum Mechanics to electromagnetism
were starting to give useful results. These developments led to the question of whether General
Relativity can be quantized.
This is a difficult question. First, Einsteins field equations are much more complicated than
Maxwells equations, and in fact are nonlinear. The physical reason for this is that the gravitational
field is coupled to itselfthe stress-energy tensor acts as the source for spacetime curvature, which
in general contributes to the stress-energy tensor. This means that the equations seem to violate the
superposition principle, which requires the existence of a linear vector space (see, e.g., [76,77]). This is
the mathematical expression of wave-particle dualitya central tenet of Quantum Theory.
Second, to quantize the gravitational field we would have to quantize spacetime itself.
The physical meaning of this is not completely clear.
Finally, there are experimental problems. Maxwells equations predict electromagnetic radiation,
which was first observed by Hertz [78]. Quantization of the field results in being able to observe
individual photons, and these were first seen in the photoelectric effect predicted by Einstein [79].
Similarly, Einsteins equations for the gravitational field predict gravitational radiation [75], so there
should be, in principle, the possibility of observing individual gravitons, which are the quanta of
the field. The direct observation of gravitational waves was finally achieved in September 2015 by
the LIGO instrument [80]. The detection of individual gravitons is more difficult and is beyond the
capability of current experiments.
To develop a quantum theory of General Relativity, the fundamental interactions in GR would
have to follow quantum rules. In Quantum Theory, particle interactions are described by gauge
theories, so GR would have to follow the gauge principle. Although the gauge principle was first
14 of 82
recognized in electromagnetism, modern gauge theory, formulated initially by Chen Ning Yang and
Robert Mills [81,82], emerged entirely within the framework of the quantum field programme. As more
particles were discovered after the 1940s, various possible couplings between those elementary particles
were being proposed. It was therefore necessary to have some principle to choose a unique form out of
the many possibilities suggested. The principle suggested by Yang and Mills in 1954 is based on the
concept of gauge invariance, and is hence called the gauge principle.
3.6. The Gauge Approach and Non-Metric Theories
The idea of gauge invariance, and the term itself, originated earlier, from the following
consideration due to Hermann Weyl in 1918 [83,84]. In addition to the requirement of General
Relativity that coordinate systems have to be defined only locally, so likewise the standard of length,
or scale, should only be defined locally. It is therefore necessary to set up a separate unit of length at
every spacetime point. Weyl called such a system of unit-standards a gauge system (analogous to the
standard width, or gauge, of a railway track).
The gauge principle therefore may be formulated as follows: If a physical system is invariant
with respect to some global (spacetime independent) group of continuous transformations G, then
it remains invariant when that group is considered locally (spacetime dependent), that is G 7 G ( x ).
Partial derivatives are replaced by covariant ones, which depend on some new vector field.
In Weyls view, a gauge system is as necessary for describing physical events as a
coordinate system. Since physical events are independent of our choice of descriptive framework,
Weyl maintained that gauge invariance, just like general covariance, must be satisfied by any
physical theory.
In Euclidean geometry, we know that translation of a vector preserves its length and direction.
In Riemannian geometry, the Christoffel connection [85] (or affine connection) guarantees length
preservation, but a vectors orientation is path dependent. However, the angle between two vectors,
following the same path, is preserved under translation. Weyl wondered why the remnant of planar
geometry, length preservation, persisted in Riemannian geometry. After all, our measuring standards
(rigid rods and clocks), are known only at one point in spacetime. To measure lengths at another
point, we must carry our measuring tools along with us. Weyl maintained that only the relative lengths
of any two vectors at the same point, and the angle between them, are preserved under parallel
transport. The length of any single vector is arbitrary. To encode this mathematically, Weyl made the
following substitution:
g ( x ) 7 ( x ) g ( x ),
(22)
where the conformal factor ( x ) is an arbitrary, positive, smooth function of position. Weyl required
that in addition to GRs coordinate invariance, formulae must remain invariant under the substitution
(Equation (22)). He called this a gauge transformation. The scale therefore becomes a local property of
the metric.
Weyls theory enabled him to unify gravity and electromagnetism, the only two forces known
at the time. However, Weyls original scale invariance was abandoned soon after it was proposed,
since its physical implications seemed to contradict experiments. In particular, if two identical clocks
C1 and C2 are transported on two different paths, which both end at the same point Q, the time-like
vectors l1 and l2 given by C1 and C2 at Q would be different in the presence of an electromagnetic
field. Therefore the two clock rates would differ. As Einstein (probably the only expert who could keep
an eye on Weyls theory at the time) pointed out, this concept meant that spectral lines with definite
frequencies could not exist, since the frequency of the spectral lines of atomic clocks would depend on
the atoms location, both past and present. However, we know the atomic spectral lines to be definite,
and independent of spacetime position [8689].
Despite its initial failures, Weyls idea of a local gauge symmetry survived, and acquired
new meaning with the development of Quantum Mechanics. According to Quantum Mechanics,
interactions are realized through quantum (that is, non-continuous) fields which underlie the local
15 of 82
coupling and propagation of field quanta, but which have nothing to do with the geometry of spacetime.
The question is whether General Relativity can be formulated as a gauge theory. This question has
been discussed by ever since it was first posed [9096].
If features of General Relativity could be recovered from a gauging argument, then that would
show that the two formulations are not inconsistent. The first to succeed in this was Kibble [91], who
elaborated on an earlier, unsuccessful attempt by Utiyama [90]. Kibble arrived at a set of gravitational
field equations, although not the Einstein equations, constructing a slightly more general theory,
known as spin-torsion theory. The inclusion of torsion in Einsteins General Relativity had long been
theorized. In fact the necessary modifications to General Relativity were first suggested by lie Cartan
in the 1920s [97100], who identified torsion as a possible physical field.
The connection between torsion and quantum spin was only made later [91,101,102], once it
became clear that the stress-energy tensor for a massive fermion field must be asymmetric [103,104].
The Einstein-Cartan (1920s) and the Kibble-Sciama (late 1950s) developments occurred independently.
For historical reasons, spin-torsion theories are sometimes referred to as Einstein-Cartan-Kibble-Sciama
(ECKS) theories, but Einstein-Cartan Theory (ECT) is the term more commonly employed.
The Einstein-Cartan Theory of gravity is a modification of GR allowing spacetime to have
torsion in addition to curvature, and, more importantly, relating torsion to the density of intrinsic
angular momentum. This modification was put forward by Cartan before the discovery of quantum
spin, so the physical motivation was anything but quantum theoretic. Cartan was influenced by
the works of the Cosserat brothers [105] who considered a rotation stress tensor in a generalized
continuous medium besides a force tensor.
Cartan assumed the linear connection to be metric and derived, from a variational principle,
a set of gravitational field equations. However, Cartan required, without justification, that the
covariant divergence of the energy-momentum tensor be zero, which led to algebraic constraint
equations, thus severely restricting the geometry. This probably discouraged Cartan from pursuing his
theory. It is now known that the conservation laws in relativistic theories of gravitation follow from
the Bianchi identities and in the presence of torsion, the divergence of the energy-momentum tensor
need not vanish.
In simple mathematical terms, a non-zero torsion tensor means that
T = 6= 0 .
(23)
Geometrically, it means that an infinitesimal geodesic parallelogram forms a non-closed loop.

Torsion is therefore a local property of the metric. The Lagrangian action of Einstein-Cartan theory
takes the usual Einstein-Hilbert form:

Z
p
g R ()
4
+ Lm ,
(24)
S=
d x g
16G
where is a general affine connection and Lm is the matter Lagrangian. The theory differs from GR in
the structure of , leading to a field theory with additional interactions.
Torsion vanishes in the absence of spin and the Einstein-Cartan field equation is then the classical
Einstein field equation. In particular, there is no difference between the Einstein and Einstein-Cartan
theories in empty space. Since practically all tests of relativistic theory are based on free space
experiments, the two theories are, to all effects, indistinguishable via the standard tests of GR.
The inclusion of torsion only results in a slight change in the energy-momentum tensor. Cartans
theory holds the distinction of being complete, self-consistent and in agreement with experiment,
but of being a non-metric theory of gravitation. The link between torsion and quantum spin means
that it could be possible to study the divergence between the GR and ECKS theories at the quantum
level. Such experiments have recently been proposed [106].
16 of 82
Kibbles theory contains some features which were criticized [107]. It is now accepted that torsion
is an inevitable feature of a gauge theory based upon the Poincar group. Classical GR must be
modified by the introduction of a spin-torsion interaction if it is to be viewed as a gauge theory.
The gauge principle alone fails to provide a conceptual framework for GR as a theory of gravity.
In the 1990s, Anthony Lasenby, Chris Doran and Stephen Gull proposed an alternative formulation
of General Relativity which is derived from gauge principles alone [108113]. Their treatment is very
different from earlier ones where only infinitesimal translations are considered [91,107]. There are a
few other theories similar in their approach to that of Lasenby, Doran and Gull (e.g., [114,115]).
4. Why Consider Alternative Theories?
The motivation for considering alternatives to GR comes mainly from theoretical arguments, like
scale invariance of the gravitational theory, additional scalar fields that emerge from string theories,
Dark Matter, dark energy or inflation, or additional degrees of freedom that arise in the framework of
brane-world theories.
In Table 1, we draw up a list of some of the more well-known alternatives to General Relativity.
This list is far from exhaustive, but it serves to highlight the major elements which differentiate these
theories. There are several works containing a more detailed listing and discussion of the various
alternative theories (e.g., [39,116,117]).
Table 1. A comparative morphology of some of the major alternatives to General Relativity,
in approximate chronological order. We have only listed the theories of particular historical significance.
The current landscape, in which cosmologists seek to explain Dark Matter, dark energy, and inflation,
offers far more theories. It is generally easier to incorporate the non-gravitational laws of physics within
metric theories, since other theories would result in greater complexity, rendering calculations difficult.
The only way in which metric theories significantly differ from each other is in their laws for the
generation of the metric. Abbreviations: Tensor (T), V (Vector), S (Scalar), P (Potential), Dy (Dynamic),
Einstein Equivalence Principle (EEP), i.e., uniqueness of freefall, Local Lorentz Invariance (LLI), Local
Position Invariance (LPI), param (Parameter), ftn (Function).
Theory
Metric
Other Fields
Free Elements
Status
Newton 1687 [20]
Nonmetric
None
Nonrelativistic, implicit action at a distance
Nordstrm 1913 [57]
Minkowski
None
Fails to predict observed light detection
General Relativity 1915 [46]
Dy
None
None
Viable
Poincar 1890s1900s [31,118]
Fails; does not mesh with electromagnetism
Whitehead 1922 [59]
Violates LLI; contradiction by everyday observation of tides
Cartan 19221925 [98]
ST
Kaluza-Klein 1920s [119,120]
Still viable; introduces matter spin
Birkhoff 1943 [61]
Fails Newtonian test; demands speed of

sound equal to speed of light
Milne 1948 [121]
Machian
background
Incomplete; no gravitational redshift prediction;

contradicts cosmological observations.
Thiry 1948 [122]
ST
Belifante-Swihart 1957 [123]
Nonmetric
BransDicke 1961 [63]
Generic S
Dy
Viable for > 500
Ni 1972 [48]
Minkowski
T, V, S
1 param, 3 ftns
Violates LPI; predicts preferred-frame effects
Will-Nordtvedt 1972 [124]
Dy T
Barker 1978 [125]
ST
Rosen 1973 [126,127]
Fixed
None
Contradicted by binary pulsar data
Rastall 1976 [128]
Minkowski
S, V
None
Contradicted by gravitational wave data
f ( R) models 1970s [129,130]
n + 1ST
Free ftn
Consistent with Solar System tests;

viable but severely constrained
MOND 1983 [131133]
Nonmetric
Free ftn
DGP 2000 [134]
ST/Quantum
TeVeS 2004 [135]
T,V,S
Extra
dimensions
Violates Equivalence Principle
Unlikely; extremely constrained by results on PPN

T
K param
Violates EEP; contradicted by DickeBraginsky experiments
Viable but can only be significant at high energy regimes

Unlikely; severely constrained.
Nonrelativistic theory
Appears to be contradicted by BAOs, CMB
and Supernovae Ia unless DE added
Dy S
Free ftn
Highly unstable [136]; ruled out by SDSS data [137]
17 of 82
5. From General Relativity to Standard Cosmology

When Einstein published his seminal GR papers it became almost immediately apparent that
the theory could be applied to the whole universe, under certain assumptions, to obtain a relativistic
cosmological description. If the content of the universe is known, then the energy-momentum tensor
can be constructed, and the metric derived using Einsteins equations. Einstein himself was the first
to apply GR to cosmology in 1917 [138]. The first expanding-universe solutions to the relativistic
field equations, describing a universe with positive, zero and negative curvature, were discovered by
Alexander Friedmann [139,140]. This occurred before Edwin Hubbles observations and the empirical
confirmation, in 1929, that the redshift of a galaxy is proportional to its distance. Hubble formulated
the law which bears his name: v = H0 r, where H0 is the constant of proportionality [141]. The problem
of an expanding universe was independently followed up during the 1930s by Georges Lematre [142],
and by Howard P. Robertson [143145] and Arthur Geoffrey Walker [146].
These exact solutions define what came to be known as the Friedmann-Lematre-Robertson-Walker
(FLRW) metric, also referred to as the FRW, RW, or FL metric. This metric starts with the assumption of
spatial homogeneity and isotropy, allowing for time-dependence of the spatial component of the metric.
Indeed, it is the only metric which can exist on homogeneous and isotropic spacetime. The assumption
of homogeneity and isotropy, known as the Cosmological Principle, follows from the Copernican
Principle, which states that we are not privileged observers in the universe. This is no longer true
below a certain observational scale of around 100 Mpc (sometimes called the End of Greatness),
but it does simplify the description of the distribution of mass in the universe.
The FLRW metric describes a homogeneous, isotropic universe, with matter and energy uniformly
distributed as a perfect fluid. Using the definition of the metric in Equation (4), it is written as:
ds2 = c d 2 R2 (t)[ dr2 + Sk0 (r )( d 2 + sin2 d2 )] ,
(25)
where r is a time independent comoving distance, and are the transverse polar coordinates, and t
is the cosmic or physical time. R(t) is the scale factor of the universe. The function Sk0 (r ) is defined as:
Sk0 (r )
sin(r )
( k = +1)
sinh(r )
( k = 0)
(26)
( k = 1)
where k is the geometric curvature of spacetime, the values 0, +1, and 1 indicating flat, positively
curved, and negatively curved spacetime, respectively.
Another common form of the metric defines the comoving distance as Sk0 (r ) r, so that

dr2
2
2
2
2
ds = c dt R (t)
+ r ( d + sin d ) ,
1 kr2
2
(27)
where t is again the physical time, and r, and are the spatial comoving coordinates, which label the
points of the 3-dimensional constant-time hypersurface.
The dimensionless scale factor a(t) is defined as
a(t)
R(t)
,
R0
(28)
where R0 is the present scale factor (i.e., a = 1 at present). The scale factor is therefore a function of
time, so it can be abbreviated to a. The metric can then be written in a dimensionless form:
h
i
ds2 = c2 d 2 = c2 dt2 a2 dr2 + Sk2 (r )( d 2 + sin2 d2 ) ,
(29)
18 of 82
where Sk (r ) can be redefined as
Sk (r ) =
R0 sin(r/R0 )
R sinh(r/R )
0
0
( k = +1)
( k = 0)
(30)
( k = 1) .
Equivalently, using the definition in Equation (27),
ds2 = c2 dt2 a2

dr2
2
2
2
2
+
r
(
d
+
sin
d
)
.
1 k(r/R0 )2
(31)
The comoving distance is distance between two points measured along a path defined at the
present cosmological time. It means that for objects moving with the Hubble flow, the comoving
distance remains constant in time. The proper distance, on the other hand, is dynamic and changes in
time. At the current age of the universe, therefore, the proper and comoving distances are numerically
equal, but they differ in the past and in the future. The comoving distance from an observer to a distant
object such as a galaxy can be computed by the following formula:
=
Z t
te
dt0
a(t0 )
(32)
where a(t0 ) is the scale factor, te is the time of emission of photons from the distant object, and t is the
present time.
The comoving distance defines the comoving horizon, or particle horizon. This is the maximum
distance from which particles could have travelled to the observer since the beginning of the universe.
It represents the boundary between the observable and the unobservable regions of the universe.
If we take the time at the Big Bang as t = 0, we can define a quantity called the conformal time
at a time t as:
Z t
dt0
=
.
(33)
0
0 a(t )
This is useful, because the particle horizon for photons is then simply the conformal time
multiplied by the speed of light c. The conformal time is not the same as the age of the universe. In fact
it is much larger. It is rather the amount of time it would take a photon to travel from the furthest
observable regions of the universe to us. Because the universe is expanding, the conformal time is
continuously increasing.
The concept of particle horizons is important. It defines causal contact. The only objects not in
causal contact are those for which there is no event in the history of the universe that could have sent a
beam of light to both. This is at the origin of some of the big questions about the universe associated
with the Big Bang model, which gave rise to the Inflationary paradigm (see [147]). We shall discuss
this later.
5.1. Cosmological Expansion and Evolution Histories
The FLRW metric relates the spacetime interval ds to the cosmic time t and the comoving
coordinates through the scale factor R(t). The scale factor is the key quantity of any cosmological
model, since it describes the evolution of the universe. The notion of distance is fairly straightforward in
Euclidean geometry. In General Relativity, however, where we work with generally curved spacetime,
the meaning of distance is no longer unique. The separation between events in spacetime depends on
the definition of the distance being used.
19 of 82
By combining the GR field equation (Equation (15)) and the definition of the metric (Equation (31)),
we obtain two independent Einstein equations, known as the Friedmann equations:
and
2
8G
a
kc2
+ 2 =
a
3
a
(34)

8G
a
a
kc2
2
+
+ 2 = 2 p.
a
a
a
c
(35)
The Friedmann equations relate the total density of the universe, including all contributions,
to its global geometry. There exists a critical density c for which k = 0. By rearranging the Friedmann
equation and using the definition of the Hubble parameter we then obtain
c ( t ) =
3H 2 (t)
.
8G
(36)
A universe whose density is above this value will have a positive curvature, that is, it will be
spatially closed (k = +1); one whose density is less than or equal to this value will be spatially open
(k = 0 or k = 1).
A dimensionless density parameter for any fluid component of the universe (i.e., a component
for whose gravitational field is produced entirely by the mass, momentum, and stress density) can be
defined by
8G(t)
(t)
=
(t) =
.
(37)
c ( t )
3H 2 (t)
The current value of the density parameter is denoted 0 .
Subtracting Equation (34) from Equation (35) yields the acceleration equation:

p
a
= 4G
+ 2 .
a
3 c
(38)
The geodesic Equation (12) allows us to compute the evolution in time of the energy and
momentum of the various components particles which make up the universe. From this evolution,
we can construct the fluid equation, or continuity equation, which describes the relation between the
density and pressure:
a
p
+ 2 = 0.
(39)
+ 3
a
c
This is valid for any fluid component of the universe, such as baryonic and nonbaryonic matter,
or radiation.
The foundations of the Concordance Model of cosmology depend on General Relativity.
Any modification to the theory that changes the Einstein equations will have solutions that differ
from the Friedmann equations.
The FLRW universe contains different mass-energy components which are assumed to evolve
independently. This is physically valid at late cosmological times, when the components are decoupled,
so the density evolutions are distinct. In Table 2, we give the equation of state and the evolution of
the density and scale factor for different components of the universe. The quantities in this table are
explained in detail in the following sections.
20 of 82
Table 2. The evolution of the various cosmological components. The quantities are the equation of
state w p/c2 , the density , the pressure p, and the scale factor a(t).
w
= a3(1+ w )
a(t ) = t 2/3(1+w)
1/3
0
1/3
1/3
2/3
1
1
a 4
a 3
a 2 a 4
a 2 a 4
a 1
1 2
2 + V ()
constant
t1/2
t2/3
t
t
t2
e Ht
e Ht
Component
Radiation (photons and relativistic neutrinos)
Dust (includes CDM, baryons and non-relativistic neutrinos)
Curvature
Cosmic strings
Domain walls
Inflation
Vacuum energy
5.2. Matter (Dust)

Matter which is pressureless is referred to as dust. This is a useful approximation for
cosmological structures which do not interact, such as individual galaxies. Substituting pm = 0
in the equation of state for dust shows that the density of this component scales as:
m ( a) =
m,0
,
a3
(40)
where m,0 is the current density. Assuming spatial flatness, the time evolution of the scale factor
is then
2/3
t
a(t) =
,
(41)
t0
which gives us
H (t) =
2
.
3t
(42)
This is known as the Einstein-de Sitter (EdS) solution, and it describes the evolution of H in a
constant-curvature homogeneneous universe with a pressureless fluid as the only component. It was
first described by Einstein and Willem de Sitter in 1932 [148].
5.3. Radiation
In the early universe, the energy content was dominated by photons and relativistic particles
(especially neutrinos). The expansion of the universe dilutes the radiation fluid, and the wavelength is
increased by the expansion so that the energy decreases. From thermodynamics,
Erad = rad c2 = T 4 ,
(43)
where T is the radiation temperature and is the Stefan-Boltzmann constant. The equation of state for
radiation can then be derived from the fluid Equation (39):
rad ( a) =
rad,0
a4
prad =
rad c2
.
3
(44)
Combining this with the Friedmann equations, and assuming flatness (k = 0), we obtain the time
dependence of the scale factor and the Hubble parameter:

a(t) =
t
t0
1/2
;
H (t) =
1
.
2t
(45)
21 of 82
6. The Components and Geometry of the Universe and Cosmic Expansion

How do we relate the expansion of the universe to its contents? The total density of the universe
in terms of its constituent components can be written as the sum of the densities of these components
at any given time or scale factor:
= m + rad + DE ,
(46)
where the subscript DE denotes another component of the universe, called Dark Energy.
10 5
Energy density [kg m 3 ]
10 10
Radiationdominated era
Matter-dominated era
Dark energydominated
era
10 15
10 20
Radiation
Matter
10 25
Dark energy
10
30
Matter-radiation equality
10 4
10 3
Present epoch
10 2
10 1
Scale factor a(t)
Figure 2. The density evolution of the main components of the universe. The early universe was
radiation-dominated, until the temperature dropped enough for matter density to being to dominate.
The energy density of dark energy is constant if its equation of state parameter w = 1. Because the
matter energy density drops as the scale factor increased, dark energy began to dominate in the recent
past. At the present time (a(t) = 1), we live in a universe dominated by dark energy. For dark
energy, the green band represents an equation of state parameter w = 1 0.2, showing how a small
change in the value of this parameter can give very different evolution histories for dark energy. If the
Concordance Model is correct, the universe will be completely dominated by dark energy in future
epochs (shown by the dashed lines). The matter density will keep decreasing as the universe expands.
Our Milky Way will merge with the Andromeda Galaxy, and eventually, the entire Local Group will
coalesce into one galaxy. The luminosities of galaxies will begin to decrease as the stars run out of
fuel and the supply of gas for star formation is exhausted. In the very far future, this galaxy will be in
the only one in our Hubble patch, as all the other galaxies will pass behind the cosmological horizon.
The night sky, save for the stars in the Local Group, will be very dark indeed. Stellar remnants will
either escape galaxies or fall into the central supermassive black hole. Eventually, baryonic matter
may disappear altogether as all nucleons including protons decay, or all matter may decay into iron.
In either scenario, the universe will end up being dominated by black holes, which will evaporate
by Hawking radiation. The end result is a Dark Era with an almost empty universe, and the entire
universe in an extremely low energy state, with a possible heat death as entropy production ceases
(see, e.g., [149,150]) What happens after that is speculative.
22 of 82
The total dimensionless density can then be written:

= m + rad + DE ,
(47)
where we have dropped the subscript for clarity, i.e., m,0 = m , etc. The Friedmann Equation (34)
can now be rewritten using the equations of state for the different components:
H 2 ( a) =
kc2
R1
0
0
8G
m a3 + rad a4 + DE e3 a [1+w(a )] d ln a 2 .
3
a
(48)
This can be rearranged to give:

i
h
R1
0
0
H 2 ( a) = H02 m a3 + rad a4 + DE e3 a [1+w(a )] d ln a + (1 ) a2 ,
(49)
or, in terms of redshift:

i
h
Rz
0
0
0
H 2 (z) = H02 m (1 + z)3 + rad (1 + z) a4 + DE e3 0 [1+w(z )]/(1+z ) d ln z + (1 )(1 + z)2 . (50)
The term 1 is sometimes replaced by k , the density due to the intrinsic geometry of
spacetime. Equation (50) is of central importance since it relates the redshift of an object to the global
density components and geometry of the universe.
The density evolution of the various components of the universe is shown in Figure 2.
7. The Hot Big Bang
In the Standard Model, it is generally accepted that the universe arose from an initial singularity,
often termed the Big Bang, which occurred some 13.8 billion years ago (as measured by Planck [151]).
This is not discussed here, but it should be noted that there are several proposals for the mechanism
of this singularity. During this epoch, we are dealing with Planck scale physics, so most of these
mechanisms involve quantum gravity. Other proposals (such as some superstring and braneworld
theories) do away with the need for an initial singularity altogether.
7.1. The Cosmic Microwave Background
The radiation density rad a4 , so the temperature evolution of the universe from an initial T0 is:
T=
T0
.
a
(51)
In other words, the universe cools down as it expands. Conversely, this means at early times,
when the scale factor was close to zero, the temperature was very high (hence the term Hot Big
Bang). The radiation left from the early hot universe, cooled by expansion, is known as the Cosmic
Microwave Background, or CMB.
The properties of atomic and nuclear processes in an expanding universe provided the first
clue for the existence of a hot Big Bang. This was a remarkable achievement of the Big Bang model,
because it provided an explanation for the observed abundances of chemical elements in terms of
nucleosynthesis. The processes that created nuclei and atoms could only have been possible in an
early universe in thermal equilibrium, with black-body spectrum which cooled down as the universe
expanded. This allowed Ralph Alpher, Robert Herman, Hans Bethe and George Gamow to predict the
existence and temperature of the CMB in 1948 [152155]. The universe therefore has a thermal as well
as an expansion history. Hence the Hot Big Bang.
The first direct evidence for the Hot Big Bang came two decades later, with the observation of the
CMB by Arno Penzias and Robert Woodrow Wilson in 1964 [156].
23 of 82
The confirmation of the thermal history of the universe, together with the discovery of charge
parity violation in 1964 [157], provided clues about baryogenesis and the observed matter-antimatter
imbalance in the universe. This inspired the first proposals for a mechanism for baryogenesis by Andrei
Sakharov in1967 [158], followed by electroweak symmetry breaking by Vadim Kuzmin in 1970 [159].
This is a remarkable demonstration of the success of the Concordance Model. The cosmological
model fits very well with the predictions of particle physics, which in turn can be tested by cosmological
observations. The Concordance Model of the structure and evolution of the universe requires a
mechanism for baryogenesis as well as an explanation for Dark Matter and dark energy. The challenge
for physical theories beyond (or within) the Standard Model is to explain the preference of matter over
antimatter, and to explain the magnitude of this asymmetry. Cosmological observations can be used to
address these challenges [160].
The CMB is an extremely isotropic source of microwave radiation, with a spectrum corresponding
to a perfect blackbody at a temperature T0 = 2.7260 0.0013 K [161]. Using the current temperature
and Erad = rad c2 = T 4 , the radiation density today is given by:
rad = 2.47 105 h2 .
(52)
At some time in the early universe, the ambient radiation temperature corresponded to the
ionisation potential of hydrogen, which is 13.6 eV. During this epoch, the universe was filled with
a sea of highly energetic particles and photonsa hot ionised plasma. The particles were mainly
electrons and protons. Other fundamental particles (quarks) existed earlier when the ambient energy
corresponded to their rest mass. At some point, as the universe expanded and cooled, the energy
of the photons was no longer sufficient to ionise the hydrogen, and within a relatively short time,
all of the electrons and protons combined to form neutral hydrogen. The photons were then free to
move through the universe. This process is known as decoupling and it occurred at a temperature
of 2500 K, when the universe was approximately 380, 000 years old [162]. It is these decoupled
photons which make up the CMB. The surface on the sky from which these photons originate is known
as the surface of last scattering.
7.2. Matter-Radiation Equality
At the present epoch, neglecting dark energy, the universe is dominated by matter. This component
is characterised by the fact that the matter particles can be treated in a non-relativistic regime,
whereas photons and relativistic neutrinos both behave like radiation. The total contribution to
the energy density from non-relativistic components (matter) and relativistic components (radiation
and relativistic neutrinos) can be written as NR and R = rad + , respectively. Using the fact
that m = m,0 a3 , the ratio of the contributions of the components is a function of the scale factor a:
R
+
4.15 105
= rad
=
,
NR
m
m,0 h2 a3
(53)
where we explicitly use the subscript 0 for the present-day values.

Then there must exist a scale factor for which the ratio is unity. This is given by:
aeq =
4.15 105
,
m,0 h2
(54)
or, in terms of redshift,

1 + zeq = 2.4 104 m,0 h2 .
(55)
The epoch at which the matter energy density equals the radiation energy density is called
matter-radiation equality, and it has a special role in large-scale structure formation.
24 of 82
7.3. Neutrinos
Neutrinos have particular properties which give rise to a distinct evolution history. They are
known to exist from the Standard Model of particle physics, and the Hot Big Bang model predicts
the amount of neutrinos in the universe. Neutrinos can be thought of as dark matter because of
their very small reaction cross-section, which implies negligible self-interaction. However, they are
not cold Dark Matter. They are simply extremely light particle that can stream out of high-density
regions. They therefore cause the suppression of perturbations on scales smaller than the free-streaming
scale. Unlike photons and baryons, cosmic neutrinos have not been observed. However, particle
physics allows us to chart the history of this particle during nucleosynthesis, and to relate the neutrino
temperature to the photon temperature today [163165].
The scale on which perturbations are damped by neutrinos is determined by the comoving
distance that a neutrino can travel in one Hubble time at equality. For a neutrino mass 1 eV,
the average velocity, T /m is of order unity at equality. This leads to a suppression of power on
all scales smaller than keq . Note that this phenomenon depends on the individual neutrino mass,
rather than the total neutrino mass. A lighter neutrino can free-stream out of larger scales, so the
suppression begins at lower k for the lighter neutrino species. Heavier neutrinos constitute more of the
total neutrino density, and so suppress small-scale power more than lighter neutrino species, which
means that we need at least two parameters to model massive neutrino phenomenology to sufficient
accuracy: the neutrino mass fraction , or some expression of this quantity in terms of the total
neutrino mass m , and the number of massive neutrino species N .
Neutrinos introduce a redshift and scale dependence in the transfer function. We know that the
perturbation modes of a certain wavelength can grow if they are greater than the Jeans wavelength.
Above the Jeans scale, perturbations grow at the same rate independently of the scale. For the baryonic
and cold Dark Matter components, the time and scale dependence of the power spectrum can therefore
be separated at low redshifts. This is not the case with massive neutrinos, which introduce a new
length scale given by the size of the comoving Jeans length when the neutrinos become non-relativistic.
In terms of the comoving wavenumber knr , this scale is given by:
knr = 0.026
m 1/2
1
1/2
m h Mpc
1 eV
(56)
for three neutrinos of equal mass, each with mass m . The growth of Fourier modes with k > knr
is suppressed because of neutrino free-streaming. From the equation above, it is evident that
the free-streaming scale varies with the cosmological epoch (since there is a dependence on m ),
and therefore the scale and time dependence of the power spectrum cannot be separated.
Neutrinos are fermions, with a Fermi-Dirac distribution with assumed zero chemical potential.
When they decoupled from the plasma, their distribution remained Fermi-Dirac, with their temperature
falling as a1 . This decoupling occurred slightly before the annihilation of electrons and positrons,
which occurred when the cosmic temperature was of the order of the electron mass (T me ). Neutrinos
decoupled when the cosmic plasma had a temperature of around 1 MeV. The energy associated
with this annihilation was therefore not inherited by the neutrinos, and the entropy was completely
transferred to the entropy of the photon background. Thus:
(Se + S )before = (S )after ,
(57)
where Se and S are respectively the entropy of the electron-positron pairs and the photon background,
and before and after refer to the annihilation time.
The entropy per particle species, ignoring constant factors, is S gT 3 , where g is the statistical
weight of the species. For bosons, g = 1 and for fermions, g = 7/8 per spin state. According to
the Standard Model, the neutrino has one spin degree of freedom, each neutrino has an antiparticle,
and there are three generations of neutrinos, also called families or species (, and electron
25 of 82
neutrinos). This means that the degeneracy factor of neutrinos is equal to 6. Before annihilation,
the fermions are electrons (2 spin states), positrons (2), neutrinos and antineutrinos (6 spin states).
The bosons are photons (2 spin states). We therefore have gbefore = 4(7/8) + 2 = 11/2, while after
annihilation g = 2 because only photons remain. Applying entropy conservation and counting
relativistic degrees of freedom, the ratio of neutrino and photon temperatures below me is therefore:
T
=
T
4
11
1/3
,
(58)
TCMB = 1.945 K .
(59)
so that the present neutrino temperature is

T,0 =
4
11
1/3
The number density of neutrinos is then

n =
6 (3) 3
T
,
11 2 CMB
(60)
where (3) 1.202, which gives n 112 cm3 at the present epoch [166]. In the early universe,
neutrinos are relativistic and behave like radiation. So they contribute to the total radiation energy
density rad , which includes the photon energy density :
"
rad
#
4/3
7
4
= 1+
Neff ,
8
11
(61)
where Neff is the effective number of neutrino species. At late times, when massive neutrinos become
non-relativistic, their contribution to the mass density is m n , giving
N
m,i
i
,
c
93.14 eV h2
(62)
where m,i is the mass of individual neutrino species and N is the number of massive neutrino species.
This expression relates the total neutrino mass m to the neutrino fraction .
It can be seen from the above that this equation can be modified through a change in the effective
number of neutrino species by many factors: a non-zero initial chemical potential, or a sizeable
neutrino-antineutrino asymmetry, or even a fourth, sterile neutrino [166168]. The Standard Model
predicts a value of Neff = 3.046 for the effective number of neutrino species. This accounts for
the three neutrino families together with relativistic degrees of freedom, since neutrinos are not
completely decoupled at electron-positron annihilation and are subsequently slightly heated [169].
Any significant deviation from this value could be a signature of hidden physical effects, possibly
requiring a modification of General Relativity [170].
Neutrino oscillation experiments do not, at present, determine absolute neutrino mass scales,
since they only measure the difference in the squares of the masses between neutrino mass eigenstates.
Cosmological observations, on the other hand, can constrain the neutrino mass fraction, and can
distinguish between different mass hierarchies [166].
Observations of neutrino flavour oscillations in atmospheric and solar neutrinos, provide evidence
of a difference between the masses of the different species or flavours, as well as for a non-zero mass.
For three neutrino mass eigenstates m1 , m2 and m3 , the squared mass differences are [171]:
26 of 82
|m221 | = m22 m21 |

= 7.5 105 eV2
|m2 | = |m2 m2 |
= 2.5 103 eV2
31
1
2
|m21 |
|m231 |
(63)
= 0.03 .
The ambiguity in the sign of the mass differences m allows for two possible mass hierarchies:
the normal hierarchy given by the scheme m3 m2 > m1 , or the inverted hierarchy m2 > m1 m3 .
Given Equation (63), constraining the total neutrino mass to a small enough maximum value could
exclude an inverted hierarchy. Conversely, a total neutrino mass m 2 eV is only possible with a
degenerate neutrino mass scheme. Hence the interest in finding cosmological neutrino mass bounds.
The fact that cosmological constraints could be stronger than constraints from particle accelerators
was noticed quite early (see [172]). The closure limit gives us m < 90 eV. This was first derived in the
late 1960s and 1970s [173176]. Since then, cosmological neutrino bounds have improved significantly,
with different methods being used e.g., luminous red galaxies [177], CMB anisotropies [151,178],
or weak lensing [179181].
Joint Planck CMB and BAO observations give us m < 0.23 eV, but various data combinations
can change this figure, and strong priors on the value of the Hubble constant can provide tighter
constraints [151].
8. Inflation: The Second Wave of Alternative Theories
In the late 1970s, General Relativity had been largely accepted by the scientific community.
But a series of cosmological considerations led to renewed interest in alternative theories. These were
not so much attempt to solve problems in the theory itself, but to find explanations for observations
that were not explained by the theory.
General Relativity applied to the universe gave us the Hot Big Bang model: a universe expanding
out of an initial highly energetic, dense state. The Hot Big Bang model was successful in explaining
many interlinked phenomena which were subsequently confirmed by observation: the Hubble Law
and the expansion of the universe, the thermal history of the universe, primordial nucleosynthesis,
the existence of the cosmic microwave background, the relation between the temperature and scale
factor, and finally the blackbody nature of the CMB. The remarkable fact is that these phenomena
occur on extremely different scales, and are observed via different physical processes, and yet they all
fit neatly within one model.
However, there are some observations which the Hot Big Bang model fails to explain.
These cosmological problems are linked to the primordial universe, the most obvious being the
following (for details see [182186], and references therein):
The Horizon Problem

The Flatness Problem
The Monopole Problem
The horizon problem arises from the structure of spacetime. In the standard cosmological model
described by the FLRW equations, different regions of the universe observed today could have not been
in causal contact with because of the great distances between them which are greater than the distance
that could have been traversed by light since the Big Bang. The transfer of information (i.e., any
physical interaction) or energy can occur, at most, at the speed of light, but these regions have the
same temperature and other physical characteristics. In particular, we observe causally-disconnected
regions of the CMB to be in thermal equilibrium. How could this have happened?
The horizon problem was first identified in the late 1960s. This led to to early attempts to model
chaotic solutions to Einsteins field equations near the initial singularity [187,188]. In the late 1970s,
Alexei Starobinsky noted that quantum corrections to General Relativity should be important in the
27 of 82
very early universe. These corrections would lead to a modification of gravity, which induces an
inflationary phase [189]. Starobinskys was the first model of inflation.
The flatness problem is one the so-called coincidence problems of modern cosmology.
By the 1960s, observations had determined that the density of matter in the universe is comparable
to the critical density necessary for a flat universe. So the contribution of curvature had to be of
the same order of magnitude as the contribution of matter throughout the history of the universe.
This represents a fine-tuning problem. Observations of the CMB have confirmed that the universe is
spatially flat to within a few percent. Why is the global geometry of the universe so flat?
The magnetic monopole problem arises from the Hot Big Bang model. Grand Unified Theories
predict the production of a large number of magnetic monopoles [190,191] in the early, extremely hot
universe. Why have none ever been observed? If they exist at all, they are much more rare than the
Big Bang theory predicts. This was noted by Zeldovich and others in the late 1970s [192,193].
Hot Big Bang Plus Inflation
This gave rise to the idea of a model in which the early universe undergoes a period of exponential
accelerated expansion. This theory, called inflation, was first formulated by Alan Guth in the
1980s [194] while he was trying to investigate why no magnetic monopoles are observed. It was
realised that inflation solves the horizon and flatness problems, as well as explaining the absence of
relic monopoles. Better still, it explains the origin of structure in the universe.
In the Standard CDM Model, the initial perturbations from which structure evolved are assumed
to have been seeded by the inflationary potential. Reconstructing the primordial power spectrum
is no easy task, and poses two main problems. Observationally, we want to extract the amplitude
and scale variation from the data. Theoretically, we seek to explain the origin of the perturbations.
At present, the leading theoretical paradigm for the primordial fluctuations is inflation, which provides
initial conditions for both large-scale structure and the cosmic microwave background radiation.
The theory of inflation offers a plethora of models, each of which predicts a certain power spectrum of
primordial fluctuations P (k). Since the inflationary paradigm is linked to the theoretical description
of the primordial power spectrum, it is necessary to briefly explain some of the main concepts here
(for the full details, see [183,185]).
The precise definition of inflation is any period during which the scale factor of the universe is
accelerating, that is, a > 0. This expression is equivalent to other definitions of inflation:
d H 1
H
d2 a
> 0 + 3p < 0 .
< 0 = e 2 < 1
dt a
H
dt2
(64)
The first expression above has a remarkable physical interpretation. It means that the observable
universe becomes smaller during inflation.
The basic theory of inflation states that from the initial Big Bang singularity to approximately
1037 s, there existed a set of highly energetic scalar fields. By definition, is driven towards 1
during inflation. Inflationary theories assume that gravity is described by GR, which means that the
component driving inflation must satisfy + 3P < 0. If for example, the universe was dominated
during the inflationary phase by a scalar field (or set of fields) with a self-interaction potential V ().
It is the form of this potential which differentiates the various inflationary theories. Most theories
assume a Mexican hat potential, with a single field, while chaotic inflation assumes a simple power
law potential with a slowly varying field [195]. The action for this potential is then [196]
S=
"
d4 x
#
m2Pl R 1
g
()2 + V () ,
16
2
(65)
28 of 82
where mPl is the Planck mass. As the universe cooled, the scalar field became trapped in a false vacuum,
so its energy density became constant. The potential energy, however, is nonzero, so the pressure is
negative. The scale factor during inflation has the de Sitter form:
a(t) = e(I /3)
1/2 t
(66)
where I represents the energy density of the inflationary field (sometimes called the inflaton).
Since the energy density of the inflaton field was very high, the associated magnitude of the
negative pressure would have been very large as well. The scale factor is thought to have increased
during inflation by e65 , and any point in the universe which found itself in a false vacuum state
would have undergone inflation. The accelerated expansion lasted until the field rolled down to a
minimum, when it decays into the familiar particles of the Standard Model, and the universe can then
be described by an FLRW model.
The inflationary paradigm provides an explanation for the origin of structure and for the observed
geometry of the universe, in addition to solving the aforementioned cosmological problems.
First, inflation solves the flatness problem. Using Equation (66), the evolution of during
inflation can be written as:
1/2
|(t) 1| e(4I /3) t ,
(67)
so that | 1| is driven very close to 0 as t increases. This explains why the universe is flat.
It also means that this value has not deviated significantly from its initial value right after expansion.
We can therefore safely assume spatial flatness throughout the history of the universe. Given the
observational difficulties, this provides a theoretical motivation for taking the idea of a large
DE seriously.
Second, inflation solves the horizon problem (Figure 3). Regions of the universe which are
causally disconnected today evolved out of the same causally-connected region in the early universe.
The observed uniformity of the CMB is no longer a problem.
Third, inflation explains why we have never observed magnetic monopoles. Due to rapid
expansion of the universe during inflation, they become so rare in any given volume of space that
we would be very unlikely to ever encounter one. Nor would they have sufficient density to alter
the gravity and thereby the normal expansion of the universe following inflation. The problem of
magnetic monopoles motivated the Guths development of his theory in 1981 [194]. The solution of
the monopole problem, and problems related to other relics, was an early success of the inflationary
paradigm, and inspired similar theories [197,198].
Fourth, the inflationary scenario provides a natural explanation for the origin of structure,
providing a link between quantum mechanics and relativistic cosmological paradigm. This was
realised soon after the development of the theory of inflation, and the details were worked out in the
early 1980s [199204].
An initially smooth background needs seed fluctuations around which gravitational collapse
can occur. The inflationary scenario attributes their origin to quantum fluctuations in the inflaton
field potential, so that the universe is not perfectly symmetric. Different points in the universe inflate
from slightly different points on the potential, separated by . Inflation for these two points ends
This induces a density fluctuation = Ht (see [184]).
at different times, separated by t = /.
Since all the points undergoing inflation are part of the same potential field, the initial fluctuations are
nearly scale invariant. This means that the density amplitude on the horizon scale will also be constant:
H = Ht =
H2
= constant .
2
(68)
In summary, inflation solves the three cosmological problems listed above:
The Horizon Problem. Solution: the entire universe evolved out of the same causally-connected region.
29 of 82
The Flatness Problem. Solution: any initial curvature is diluted by the inflationary epoch and
driven to zero.
The Monopole Problem. Solution: the rapid expansion of the universe drastically reduces the
predicted density of magnetic monopoles, if they exist.
Conformal time
st
Pa
ne
co
ht
lig
Last scattering surface
rec
0
Recombination
Reheating
Inflation
Particle horizon
Causal contact
initial =
Big Bang singularity
Figure 3. How inflation solves the horizon problem. The light cones on the causal diagram of an
inflationary FLRW model are at 45 . The worldlines of comoving matter are vertical on this kind of
diagram. The particle horizons are horizontal lines. Here we have shown the particle horizon for the
CMB. Without inflation, conformal time would only go back to 0 , and different regions of the CMB
which we observe today along our past light cone would never have been in causal contact. Because of
inflation, conformal time is extended to the Big Bang singularity, so these regions would have been in
causal contact at some point in our past light cone.
How long did inflation last? The answer is given by looking closely at Equation (67) above.
A convenient measure of expansion is the so-called e-fold number, defined as:

N ln
af
ai
Z t
f
ti
H dt .
(69)
Here, ai and af are the values of the scale factor at the beginning and end of inflation, while ti
and tf are the corresponding proper times. The scale factor a is only physically meaningful up to a
normalisation constant, so the e-fold number is defined with respect to some chosen origin. The reason
30 of 82
is that in cosmology, what is fixed is not the initial condition, but the current expansionwe cannot
measure any H, but we measure H0 then extrapolate backwards.
We can search for the minimum duration of inflation required to solve the horizon problem.
At the very least, we require that the observable universe today fits in the comoving Hubble radius at
the beginning of inflation:
( a0 H0 )1 < ( ai Hi )1 .
(70)
The condition is the same for the horizon and flatness problems.
If we assume that the universe was radiation-dominated since the end of inflation (giving us
H a2 ), and ignore the relatively recent matter- and dark energy-dominated epochs, we obtain
af
a
= 0.
ai
af
(71)
In the general case, the condition becomes:

af
a
0,
ai
af
(72)
Nf Ni N0 Nf .
(73)
or in terms of the number of e-folds,
In other words, there should be as much expansion during inflation as after inflation.
The solution to the horizon problem is the same as the solution to the flatness problem. Taking into
account the present energy density of the universe, we need a minimum of about 50 to 60 e-folds. This
already gives us a useful criterion for realistic inflation models. The most recent Planck results show a
preference for a higher number of e-folds: 78 < N < 157 [205].
Most models of inflation are slow-roll models, in which the Hubble rate varies slowly [185,206,207].
This model was first developed by Andrei Linde in 1982 [208]. It solved a major problem in Guths
early theory. Instead of tunnelling out of a false vacuum state, inflation occurrs by a scalar field rolling
down a potential energy gradient. When the field rolls very slowly compared to the expansion of
the universe, inflation occurs. Hence the name slow-roll inflation. However, when the gradient
becomes steeper, inflation ends and reheating can occur. It is beyond the scope of this review to go
into the detail of the theory, but it is necessary for us to briefly refer to the link between this theory and
the spectral index of primordial fluctuations, which is an important observational parameter in the
Concordance Model of cosmology.
To quantify slow roll, cosmologists typically use two parameters e and which vanish in the limit
that becomes constant. The first parameter is defined as:
d
e
dt
1
H
H
,
aH 2
(74)
which is always positive, since H is always decreasing. The second complementary variable which
defines how slowly the field is rolling is:
h
i
1
(0) + a 2 V 0 ,
3aH
aH (0)
(75)
where (0) is the zero-order field, and V is the potential.

The scalar spectral index can be defined in terms of some function, usually a polynomial, involving
the two slow-roll parameters e and . As an example we shall give two such parameterisations:
n = 1 4e 2 [209] and n = 1 6e + 2 [206]. The rate of change of n can also be expressed in terms
of inflationary parameters: dn/ d ln k = 16e + 24e2 + 2 2 [210], where
31 of 82
2 m4Pl
V 0 ( d3 / d3 )
,
V2
(76)
mPl being the reduced Planck mass (4.342 106 g).

Therefore, by extracting the values of e and from the data, using methods such as weak lensing,
we can directly probe the potential of of the inflaton field. Likewise for the tilt or spectral index
of the primordial power spectrum. Slow-roll inflation predicts that the spectral index of primordial
fluctuations should be slightly less than 1. The reason for this is simple. For inflation to end, the Hubble
parameter H has to change in time. This time-dependence changes the conditions at the time when
each fluctuation mode exits the Hubble horizon and therefore gets translated into a scale-dependence.
Inflation accounts for the observed spatial flatness of the universe, and the absence of
magnetic monopoles. These predictions have been confirmed by various probes, most notably by
precision measurements of CMB anisotropies, starting with the Cosmic Background Explorer (COBE)
in 1992 [211213], then with the Wilkinson Microwave Anisotropy Probe (WMAP) which ran for nine
years from 2001 to 2010. WMAP data placed tight constraints on the predicted burst of growth in the
very early universe, providing compelling evidence that the large-scale fluctuations are slightly more
intense than the small-scale ones, which is a subtle prediction of many inflation models [162,214219].
Significantly, WMAP found evidence that the scalar spectral index is less than 1 (a 2 deviation),
implying a deviation from scale invariance for the primordial power spectrum. As explained above,
this is a major prediction of inflation, and this observation reinforced the evidence in favour of the
theory. Conclusive proof of a scale-dependent primordial power spectrum (a 5 deviation from ns = 1)
was provided by the Planck CMB anisotropy probe in 2013 [220,221] and confirmed in 2015 [151,205].
One current experimental challenge is to observe the B-modes of polarisation of the CMB caused
by primordial gravitational waves produced by inflation. Their detection by the BICEP2 experiment
was announced in early 2014. However, more accurate modelling of the signal over the next few
months, which allowed the observation to be explained by polarised dust emission in our Galaxy,
decreased the statistical confidence of the initial result [222]. This was confirmed by Planck data
in 2016 [223]. Upcoming large-scale structure surveys, such as the Euclid satellite mission, or 21-cm
radiation surveys such as the Square Kilometre Array, may measure the power spectrum with greater
precision than current CMB probes, and could provide further evidence in favour of the inflationary
paradigm [224226].
9. The First Unknown Component: Dark Matter
The first evidence for Dark Matter came from astronomy rather than cosmology.
Newtonian physics and General Relativity both provide very precise rules for the dynamics of
galaxies: the mass determines the rotation velocity. Starting in the 1920s, stronomers noticed that
amount of visible matter in galaxies did not match the observed rotation curves. These curves relate
the tangential velocity of the constituent stars (or gas) about the centre of the galaxy to their radial
distance. Observations of the velocities of globular clusters about galaxies showed that at large radii
the velocities are approximately constant, implying that the amount of mass in the galaxies is much
higher than the visible mass.
The first suggestion of the existence of hidden matter, motivated by stellar velocities, was made by
Jacobus Kapteyn in 1922 [227]. Radio astronomy pioneer Jan Oort also hypothesized the existence of
Dark Matter in 1932 [228]. Oort was studying stellar motions in the local galactic neighbourhood and
found that the mass in the galactic plane must be greater than what was observed. This measurement
was later determined to be erroneous.
In 1933, Fritz Zwicky, who studied galactic clusters while working at the California Institute of
Technology, made a similar inference [229]. Zwicky applied the virial theorem to the Coma galaxy
cluster and obtained evidence of unseen mass that he called dunkle Materie in German, or Dark
Matter. Zwicky estimated its mass based on the motions of galaxies near its edge and compared that
32 of 82
to an estimate based on its brightness and number of galaxies. He estimated that the cluster had about
400 times more mass than was visually observable. The gravity effect of the visible galaxies was far too
small for such fast orbits, thus mass must be hidden from view. Based on these conclusions, Zwicky
inferred that some unseen matter provided the mass and associated gravitation attraction to hold the
cluster together.
In 1937, Zwicky made the bold assertion that galaxies would be unbound without some form of
invisible matter [230]. Zwickys estimates were off by more than an order of magnitude, mainly due to
an obsolete value of the Hubble constant. The same calculation today shows a smaller fraction, using
greater values for luminous mass. However, Zwicky did correctly infer that the bulk of the matter
was dark.
More evidence started to accumulate for the existence of some non-emitting component which
was now being called Dark Matter [231]. In 1959, Kahn and Woltjer [232] pointed out that the
motion of Andromeda towards us implied that there must be Dark Matter in our Local Group of
galaxies. Dynamical evidence for massive Dark Matter halos around individual galaxies came later,
starting in the 1970s, when rotation curve data from multiple galaxies confirmed the Dark Matter halo
hypothesis [233236]. Like baryonic matter, Dark Matter is a fluid with vanishingly small pressure.
Unlike baryonic matter, it has no interaction with photons, making it both dark and transparent. It also
has a vanishingly small self-interaction beside gravity. One result of this is that the Dark Matter halos
surrounding galaxies are rounder than the galaxies themselves [237] .
In the last few decades, cosmology has contributed one important piece of information:
the amount of Dark Matter. The observed value of the matter density in the universe is
m = 0.3089 0.0062. However, the density of baryonic matter is b = 0.0486 0.0010 [151]. The
missing mass is made up of Dark Matter.
The name Dark Matter is an indication of its nonbaryonic nature: it cannot be observed by
emission of photons, so observers need to find a way around this problem. Current evidence for the
existence of Dark Matter comes from a variety of sources besides galactic dynamics [231]. The two most
important ones are CMB anisotropies and gravitational lensing. In addition, Big Bang nucleosynthesis
provides evidence that some of the Dark Matter may be baryonic. The inventory of observed baryons
in the local universe falls short of the total anticipated abundance from Big Bang nucleosynthesis,
implying that most of the baryons in the universe are unseen [238].
Anisotropies in the CMB are related to anisotropies in the baryonic density field by the Sachs-Wolfe
effect [239]. This means that the baryon density field variation at the time of decoupling can be linked
to CMB anisotropies. If all matter were made of baryons, the amplitude of the density fluctuations
should have reached 102 at the present epoch. However, we observe structures with 1 at
the present epoch (e.g., galaxies and galaxy clusters). The discrepancy can only be explained by the
presence of additional matter, which created potential wells for the baryons to fall into after decoupling.
These potential wells would have had to be formed by a weakly interacting fluid that decoupled well
before baryons and began to cluster much earlier. Such a fluid would only interact via the gravitational
and possibly the weak nuclear force. As the baryons accumulated in the potential wells, their pressure
would have built up, leading to oscillations in the baryon fluid, termed baryon acoustic oscillations
(BAO) [240,241]. These oscillations leave an imprint on the CMB power spectrum, which has been
confirmed observationally, and which constrains the mass density, leading to a further confirmation of
the existence of this missing mass.
The phenomenon of gravitational lensing includes cosmic shear, weak lensing, cosmic
magnification. Although the theory of cosmic shear had been worked out from the 1960s to the
early 1990s [242], the first detection had to await the development of instruments sensitive enough to
make the required observations, and image analysis software to accurately correct for unwanted effects
when measuring the shapes of galaxies. In March 2000, four groups independently announced the first
discovery of cosmic shear [243246]. Since then, cosmic shear has established itself as an important
technique in observational cosmology.
33 of 82
Gravitational lensing shows that the amount of lensing of galaxies around galaxy clusters is too
high to be caused by the visible matter. Apart from the stars themselves, a galaxy cluster also has a gas
component, but X-ray observations show that this is still not enough to account for the extra mass.
The cluster must therefore have a non-emitting halo of Dark Matter around it.
Various Dark Matter candidates have been proposed [247251]. However, all these candidates
have one common characteristic: a very small reaction cross-section, making them extremely difficult
to detect directly [252255]. Experiments have, however, placed limits on the mass of Weakly Interacting
Massive Particles (WIMPs), which are the current best candidate for Dark Matter (together with axions).
WIMPs are an entire new class of fundamental particle outside of the Standard Model that result from
supersymmetry [256,257]. These results show that even the lightest Dark Matter particle should have a
mass which is not below 10 MeV. We also know that CDM = 0.2589 0.0057 [151]. The conclusion
is that m , implying that hot Dark Matter (i.e., neutrinos) cannot account for the Dark Matter
density CDM .
An alternative to Dark Matter is to explain the missing mass by means of a modification
of gravity at large distances or more specifically at small accelerations. In 1983, Morderhai
Milgrom proposed a phenomenological modification of Newtons law which fits galaxy rotation
curves [131133]. The theory, known as Modified Newtonian Dynamics (MOND) automatically
recovers the Tully-Fischer law. The theory modifies the acceleration of a particle below a small
acceleration a0 1010 ms2 , which therefore enters the theory as a universal constant. The gravitational
acceleration
at large distances then reads a = GMa0 /r at large distances, instead of the Newtonian
law a = GM/r2 .
There are two main difficulties with MOND. First, it does not explain how galaxy clusters can
be bound without the presence of some hidden mass [258,259]. Second, attempts to derive MOND
from a consistent relativistic field theory have failed. One such attempt is the Tensor-Vector-Scalar
Theory (TeVeS) [135] is more successful and actually relativistic but not apparently necessary since it
still requires dark energy and Dark Matter. Many models are unstable [260], or require actions which
depend on the mass M of the galaxy, thereby giving a different theory for each galaxy. Moreover,
modified gravity theories have serious difficulties reproducing the CMB power spectrum and the
evolution of large-scale structure [261263].
The greatest challenge to modified gravity theories, and also the clearest direct evidence of Dark
Matter, comes from observations of a pair of colliding galaxy clusters known as the Bullet Cluster [264]
in which the stars and Dark Matter separate from the substantial mass of ionised gas. The Dark Matter
follows the less substantial stars and not the more massive gas.
The modifications of gravity proposed as alternative to the Dark Matter paradigm illustrate the
need for tests of GR at large distances and low accelerations. They also illustrate the problems faced by
models which favour goodness of fit over parsimony. Modified gravity theories can give an excellent
phenomenological fit through an adjustment of the values of the extra parameters, but there is no
universal principle to determine these values. This requirement for simplicity and predictivity is met
by General Relativity.
10. The Second Unknown: Dark Energy and the New Wave Alternative Theories
The current motivation for alternative theories seems to be the search for an explanation of the
observed accelerating expansion of the universe. Let us consider the justification for the dark energy
paradigm within the inflationary CDM model, and the process which led to its acceptance by the
scientific community (see [265]).
Round about the time that GR was developed, the universe was thought to be static. There was
no compelling reason to think otherwise. Einstein realised that his equations implied a non-static
universe, so in 1917 he revised his field equations of GR to read [138]:
G g = G 8G g = 8GT
(77)
34 of 82
where = /8G is proportional to the cosmological constant .

It can be seen from this equation that Einstein did not consider the cosmological constant to be
part of the stress-energy term. One could, of course, put g on the right-hand side of the equation
and count it as part of the source term of the stress-energy tensor and simply consider to be the
vacuum energy. This is not just a semantic distinction. When takes part in the dynamics of the
universe, then the field equation is properly written with , or its generalisation, as part of the
stress-energy tensor:
G = 8G ( T + g ) .
(78)
The equation describing gravity is then unchanged from its original formthere is no new
physical theory. Instead, there is a new component in the content of the universe.
This component must satisfy Special Relativity (that is, an observer can choose coordinates so that
the metric tensor has Minkowskian form). An observer moving in spacetime in such a way that the
universe is observed to be homogeneous and isotropic would measure the stress-energy tensor to be
=
0
0
0
p
0
0
0
0
p
0
0
0
.
0
p
(79)
This means that the new component in the stress-energy tensor looks like an ideal fluid with
negative pressure:
p = .
(80)
In modern concordance cosmology, this component is usually termed dark energy. If the
equation of state parameter of dark energy is constant, i.e., w(z) = 1, then its energy density will be
constant regardless of the expansion of the universe.
Einstein inserted the cosmological constant because he felt that the non-static universe predicted
by the formalism of GR was incorrect, given the data available in 1917 [138]. At the time, observations
of the universe were limited primarily to stars in our own galaxy, with observed low velocities, so there
was solid observational evidence justifying the assumption that the universe was static. Einsteins
goal was to obtain a universe that satisfied Machs principle of the relativity of inertia. However,
observational evidence started to accumulate for another paradigm. In 1917, Vesto Slipher [266]
published his measurements of the spectra of spiral nebulae, which showed that most were shifted
towards the red. The breakthrough came when the linear redshift-distance relation was formulated by
Hubble [141], who showed that the universe was expanding. Einstein then dropped his support for
the cosmological constant.
In the FLRW cosmological model, the expansion history of the universe is determined by the
mass density of the different components, whose sum is normalised to unity:
m,0 + rad,0 + X,0 + k,0 = 1 ,
(81)
where the 0 subscript indicates the present epoch. We use the term X to show that this equation does
not assume anything about the nature of the additional energy component (dark energy). In fact we
could have used or DE in the current Concordance Model.
Big bang nucleosynthesis and observations of large scale structure provide a good determination of
the mass content of the universe, allowing m and rad to be fixed. However, observations in the 1980s
and 1990s started to show inconsistencies with the cosmological model at the time, which was a
matter-dominated, expanding universe with a present-epoch Hubble constant of H0 ' 0.7 kms1 Mpc
and = 0 [267,268]. This was the so-called age problem, where the predicted age of the universe
seemed to be younger than the age of the oldest stars. Angular-diameter distances to the last scattering
surface at z = 1100 measured from the CMB are in fact 1.7 times smaller than those predicted by
35 of 82
an isotropic and homogeneous universe containing only pressureless matter (see [269]). Since the
inflationary scenario, which by then was well established, predicts a flat total = 1 universe, there was
a problem with the cosmological model.
It was realised that one of the three assumptions of the cosmological model had to be wrong.
Either the universe contains exotic matter with a negative pressure, or standard General Relativity
is wrong, or the universe is not homogeneous and isotropic. The solution could also lie in some
combination of the three. Most of the research since the late 1990s has followed the first approach, and
the term Concordance Model refers to an FLRW universe, following General Relativistic cosmology,
containing dark energy.
Within the FLRW framework, two main proposals were put forward: one was CDM, in which
there is a contribution to the energy density from a term similar to the cosmological constant (or the
cosmological constant itself), and the other was + CDM, where the missing mass came from massive
neutrinos (m ' 7eV) (e.g., [270,271]).
The first strong evidence of dark energy came in 1998 and 1999, when observations of the
luminosities of type Ia supernovae indicated that the expansion of the universe is accelerating [272,273].
Concurrently, other observations constrained the neutrino mass to m 7eV, thus discounting the
+ CDM model and confirming CDM as the Concordance Model (e.g., [265,274]). It is not clear
when the term dark energy was first used, but it seems to have been around 1998. The term is
analogous to Dark Matter, which had been in use for some time [275].
Since then, numerous observations have confirmed cosmic acceleration, including supernovae,
the cosmic microwave background, large-scale structure and baryon acoustic oscillations
(see, e.g., [276278]). The values of the present epoch matter and radiation components are
well established:
m,0
8Gm,0
0.3
3H02
rad,0
8Grad,0
1 104 ,
3H02
(82)
where H0 is the present value of the Hubble parameter H ( a = 1).

The data also indicate that the universe is currently nearly spatially flat:
|K | 1 .
(83)
This is normally taken to imply that the spatial curvature K = 0, since

k,0 = 0
K
0.
a20 H02
(84)
Thus it also justifies the inflationary paradigm. However, inflation only tells us that K 0,
so that the curvature may have had a nonzero value in the past. In the present universe, however,
the distinction is negligible. In any case, Equation (81) implies that there has to be a nonzero
(a constant term added to the Einstein equation) such that
,0
0.7 .
3H02
(85)
Inserting these values into the Friedmann equation leads to the dramatic conclusion that the
expansion of the universe is accelerating:
a 0 =
H02
1
m rad
2
where a0 is the present value of the scale factor a(t).
> 0,
(86)
36 of 82
Note that this conclusion only holds if the universe is homogeneous and isotropic (i.e.,
a Friedmann-Lematre model). In such a universe, the distance to a given redshift z and the time
elapsed since that redshift are tightly related via the only free function, a(t). If the universe is isotropic
around us, but not homogeneous, that is, a non-Copernican Tolman-Bondi-Lematre model [279],
then this relation would be lost and present data might not imply acceleration. A Copernican model
where this relation again breaks down is the inhomogeneous universe, where the acceleration can be
produced via nonlinear averagingthe backreaction of inhomogeneities.
Dark energy is a fluid component whose equation of state is:
pDE = wc2 DE .
(87)
This is the equation of state in its most general form, since w can be any function of redshift, scale
factor or cosmic time, with the constraint that w 0 (i.e., the fluid has a negative pressure). Assuming
that w = w( a), we have the following density-scale relation:
DE ( a) = DE,0 e3
R1
a
[1+w( a0 )] d(ln a0 )
(88)
It can be seen that in the special case of a constant w = 1, the fluid equation implies that the
density is constant.
If dark energy is a cosmological constant, it still leaves the question of its physical nature.
Its observed value of 3 10122 c3 /h G is so small that it is hard to interpret as the vacuum energy.
One possibility is that the final value of the cosmological constant is zero, and that cosmic
acceleration is due to the potential energy of a scalar field, with some sort of mechanism to dynamically
relax it to a small value. This notion leads to models of dark energy which invoke a slowly-rolling
cosmological scalar field to source accelerated expansion, smilar to cosmological inflation.
How can a scalar field drive cosmic acceleration? The action of a scalar field minimally coupled
to Einstein gravity is
!
Z
p
m2Pl R 1
4
2
S=
d x g
() V () .
(89)
2
2
The stress-energy tensor for the scalar field is given by
= g
1
()2 + V ()
2

.
(90)
For a homogeneous field such that = (t), a cosmological scalar acts like a perfect fluid with
equation of state w = P/ given by
1 2
V ()
w = 21
.
(91)
2 + V ()
2
The observed expansion leads us to a value of w ' 1, which requires a very slowly-rolling
field: 2 V (). There has also been a lot of interest in constructing quintessence models which can
produce an equation of state of the phantom type (w < 1).
The equation of state of dark energy that has the potential to distinguish between dark energy
candidates. The most important distinction that can be made between different dark energy models is
whether the energy density of this component is constant, filling space homogeneously, or whether it is
some form of quintessence field whose energy density can vary in time and space. There is a multitude
of alternative models, such as f ( R) [280282], or Chameleon Models [283,284]. It is therefore useful to
consider the redshift evolution of w, so that w is an arbitrary function of redshift z. There are a number
of different parameterisations of w(z) ([285], and references therein).
37 of 82
The most common parameterisation is sometimes termed the Chevallier-Linder-Polarski (or CPL)
parameterisation [285,286]:
dw
w ( z ) = w0
(1 a ) ,
(92)
da
where the scale factor a = 1/(1 + z). If we define
wa =
dw
,
a d ln a
(93)
then the equation becomes

w ( a ) = w0 + w a ( 1 a ) ,
(94)
which is the most commonly used form.

It has been shown that this parameterisation is stable and robust over large redshift ranges [287].
A wide range of functional forms of w( a) can be parameterised by the w0 wa combination. However,
there are some dark energy models which it cannot reproduce (see [288291]).
The problem with the dark energy paradigm, stated simply, is that the parameters are not
constrained well enough to rule out certain models. We have fairly good bounds on the dark energy
density, but the dark energy equation of state is still poorly constrained. Even for a constant w model,
corresponding to CDM, the bounds are such that a time-varying w( a) could mimic a constant w,
thereby disguising underlying physics (see [292294]).
11. The Evolution of Large-Scale Structure
After the epoch of matter-radiation equality, and before the onset of dark energy domination,
the mass-energy content of the universe became dominated by matter. From an initially smooth
background (as evidenced by CMB observations), structures have evolved to a scale of more than
100 Mpc, with the term large scale structure being used to refer to objects modelled on this scale.
At this scale, the mass-energy inhomogeneities can be modelled as perturbations on a homogeneous
and isotropic unperturbed background spacetime. Below this scale we observe galaxy clusters,
individual galaxies, and stars. The model of structure formation must be accurate enough to provide
an good description of the universe on a wide range of scales.
The standard model for the formation of structure assumes that at some early time there existed
small fluctuations, which grew by gravitational instability. The origins of these fluctuations are unclear,
but they are thought to arise from quantum fluctuations of the primordial universe, uncorrelated
and with Gaussian amplitudes, which were then amplified during a later inflationary phase [185].
The assumption that the amplitudes of the relative density contrasts is much smaller than unity means
that we can think of the primordial fluctuations as small perturbations on a homogeneous and isotropic
background density. This ensures that we can describe them using linear theory.
Heuristically, the mechanism of structure formation can be understood in terms of gravitational
self-collapse. Matter collapses gravitationally around initial mass overdensities. This increases the
relative density of that region, causing further collapse of more matter, and amplifying the effect.
The linear theory of structure formation needs to be relativistic, because the perturbations on any
length scale are comparable or larger than the horizon size at sufficiently early times. The horizon size
is defined as the distance ct which light can travel in time t since the Big Bang. Dissipative effects and
pressure also affect structure formation, as explained below (for details of the theory, see [295,296]).
The relative density is the density at a particular point in space x relative to the mean at some
time parameterised by the scale factor a, and can be expressed as a dimensionless density contrast:
(x, a) =
(x, a) ( a)
.
( a)
(95)
38 of 82
This quantity can be understood as the dimensionless density perturbation of some background
matter distribution.
There are two types of density perturbations that can occur within a matter-radiation fluid. If the
fluid could be compressed adiabatically in space, the perturbations have a constant matter-to-radiation
ratio everywhere. Since the energy density of radiation is proportional to T 4 , and the number density
is proportional to T 3 , the energy densities of radiation and matter are related by:
rad =
4
m .
3
(96)
Isocurvature perturbations occur when the entropy density is perturbed, but not the energy
density. Since the total energy density remains constant, there is no change in the spatial curvature and
rad rad = m m .
(97)
Perturbations can occur at different scales, or modes. The latter term is used when the amount
of perturbation on a particular scale is expressed using Fourier analysis. The Fourier transform pair of
(x) is:
(k) =
(x) =
d3 x(x)eik.x
d3 k
(k)eik.x ,
(2 )3
;
(98)
with each mode assumed to evolve independently. In the Einstein-de Sitter regime, linear adiabatic
perturbations scale with time as follows:
(
a ( t )2
(radiation domination)
a(t)
(matter domination)
(99)
while isocurvature perturbations are initially constant and then decline:

(
constant
(radiation domination)
a ( t ) 1
(matter domination) .
(100)
In both cases, the overall shape of the spectrum of the perturbations over all modes is preserved,
while the amplitude changes with time. The evolution described above is affected on small scales by a
number of processes.
11.1. Evolution on Small Scales
During the radiation-dominated epoch the growth of certain modes is suppressed. This behaviour
can be modelled in terms of the horizon scale H ( a), which is the distance ct that light could have
travelled since the initial singularity (a comoving horizon size). A mode k is said to enter the horizon
when = H ( a), where = (2 )/k. If < H ( aeq ) then a mode enters the horizon during the
radiation-dominated epoch. The time scale for collapse of matter during this epoch is larger than the
typical expansion time scale (t 1/H ( a)) due to the relatively rapid expansion rad a4 . The growth
of these modes is therefore suppressed. After the epoch of matter-radiation equality (a = aeq ), these
perturbations can then start to collapse gravitationally. We can define the suppression factor for a
particular mode as the factor by which the amplitude is reduced had it not entered the horizon:

f sup =
aenter
aeq
2
k0
k
2
(101)
39 of 82
where the mode evolves as a2 until it enters the horizon at aenter and is suppressed until aeq , when
its evolution resumes as a. The second equality in the above equation comes from applying an
Einstein-de Sitter approximation where k0 = 1/ H ( aeq ) (see [297]).
Pressure opposes gravitational collapse for modes with a wavelength less than the Jeans
length [298], sometimes called the free-streaming scale, defined as
J = cs
.
G
(102)
During the radiation-dominated epoch, the sound speed cs = c/ 3 and the Jeans length is always
close to the horizon size. The Jeans length then reaches a maximum at a = aeq and then begins to
decrease as the sound speed declines. This means that on scales larger than the comoving horizon size,
perturbations are only affected by gravity, and the spectrum starts to turn over at this point (where the
effects of pressure begin to dominate). The comoving horizon size at zeq is given by:
R0 r H (zeq )
16.0
Mpc .
m h2
(103)
Another important scale occurs where photon diffusion erases perturbations in the
matter-radiation fluid. This process is termed Silk damping [299]. The scale at which it occurs
is characterised by the distance travelled by the photon in a random walk by the time of last scattering:
S 16.3(1 + z)5/4 (2b m h6 )1/4 Gpc .
(104)
All of the effects mentioned above are particularly important where the behaviour of massive
neutrinos is concerned. Heuristically we can understand the complexity of their behaviour by
considering them as a component whose equation of state changes as the universe evolves. From a
component which behaves like photons (since the particles have a very small mass and relativistic
speeds), massive neutrinos lose energy and start behaving like baryonic matter [300,301].
11.2. Growth oF Perturbations in the Presence of Dark Energy
All of the above effects were described in an Einstein-de Sitter universe. In a universe with a
smooth non-clustering dark energy component below the horizon scale, the matter perturbation fields
evolves according to:
+ 2H (3/2) H 2 m = 0
00 + (2 q) a1 0 (3/2)m a2 = 0,
(105)
where a dot denotes a time derivative and a dash denotes a derivative with respect to a. The term q
is the deceleration parameter. This can be interpreted in the following way: the perturbations grow
according to a source term which involves the amount of matter (m ) but the growth is suppressed by
the friction term due to the expansion of the universe. The latter is also known as the Hubble drag.
If we define the growth as the ratio of the amplitude of a perturbation at a time a to some initial
amplitude, i.e.,
( a)
D ( a) =
,
(106)
( ainitial )
the equation becomes, for a general dark energy scenario where w = w( a) (see [302])
D 00 +
3
2

1
w( a)
1 + X ( a)
D0 ( a) 3
a
2
X ( a)
1 + X ( a)
D
= 0,
a2
(107)
40 of 82
where
X ( a) =
m 3 R 1
a
e
DE
d ln a0 w( a0 )
(108)
is the ratio of the matter density to the dark energy density. For large X (i.e., m 1 where
DE 1 m ) we recover the matter-dominated behaviour (D a). To parameterise deviations
from this behaviour we define the normalised growth as G = D/a. The evolution equation is then:

7 3
G +
2 2
00
w( a)
1 + X ( a)

G 00
3
+
a
2
1 w( a)
1 + X ( a)
G
= 0.
a2
(109)
This equation allows us to physically interpret the effects of dark energy. In the presence of
dark energy, the Hubble drag term is increased, so that growth is suppressed in a universe with an
accelerating expansion. This is similar to the suppression due to radiation dominance.
11.3. The Power Spectrum of Matter
In an FLRW universe, the homogeneity and isotropy assumption means that any statistical
properties must also be homogeneous and isotropic. The implication for the matter perturbation field
is that its Fourier modes must be uncorrelated (due to homogeneity). Usually, we assume that the mode
amplitudes are Gaussian. This assumption is well motivated since the theory for the seed fluctuations
assumes that they have a quantum origin. Due to the central limit theorem, the sum of a sufficiently
large number of mode amplitudes will tend towards a Gaussian distribution (see, e.g., [303,304]).
Such a field, with uncorrelated modes, and a Gaussian distribution of mode amplitudes is called
a Gaussian random field, and can be entirely described by its two-point correlation function:
h(x) (y)i = C (|x y|) .
(110)
The angled brackets denote an ensemble average (an average over a multitude of realisations).
The value of at a given point in the universe will have a different value in each realisation, with a
variance h2 i. Since we can only observe one realisation of our universe (in other words, at most only
a finite region in this one universe), we apply the ergodic principle: The average over a sufficiently
large volume is equal to the ensemble average.
In Fourier space, the correlation function can be written as:
h(k) (k0 )i =
d3 xeik.x
d3 x 0 eik .x h(x) (x0 )i.
(111)
Replacing x0 = x + y, and substituting Equation (110), this can be written as:
h(k) (k0 )i =
d3 xeik.x
d3 yeik .(x+y) C (|y|)
= (2 )3 D (k k0 )
d3 yeik.(y)C (|y|)
= (2 )3 D (k k0 ) P (|k|) .
(112)
(113)
(114)
The power spectrum has been defined as the Fourier transform of the correlation function:
P (|k|) =
d3 yeik.(y) C (|y|) .
(115)
The standard convention in cosmology is to abbreviate P (|k|) to P(k), where k = |k|. The power
spectrum can be expressed in dimensionless form as the variance per ln k, so that:
2 ( k ) =
k3 P ( k )
.
2 2
(116)
41 of 82
11.3.1. Nonlinear Evolution

The power spectrum gives us the evolution of the initial matter density fluctuations. However,
the linear evolution breaks down at small scales, when complex structures begin to form,
and overdensities can no longer treated as perturbations on a smooth background. This is the
nonlinear regime of gravitational evolution. The scale above which nonlinearities cannot be ignored
is approximately set by (kNL ) ' 1, which corresponds to kNL ' 0.2 h Mpc1 in most cosmological
models. The standard way to model nonlinear evolution is by using phenomenological fits based on
N-body simulations.
One strategy is to build a models based on a stable clustering hypothesis [305307], which assumes
that the nonlinear collapsed objects form isolated, virialised systems that are decoupled from the
expansion of the universe.
A different approach is used in the halo model [308]. Here, the density field is decomposed
into individual clumps of matter with some density profile and varying mass. By using this model
to calculate the number of clumps within a given volume, the galaxy halo profile can be calculated.
This is the equivalent of the power spectrum for these matter halos. A functional relation between
the linear power spectrum and this halo profile is then derived and calibrated using large N-body
simulations. This relation is then used to calculate the nonlinear power spectrum, using a fitting
formula, for instance (see, e.g., [309,310]).
Whichever approach is used, it must account for a range of small-scale physical processes, such as
baryonic physics, stellar formation, galactic magnetic fields, and Dark Matter and neutrino properties,
which have become more significant with the ability of future astrophysical experiments to probe the
nonlinear regime with ever increasing precision (see, e.g., [311317]).
The need for realistic models of the universe has come to the fore in recent years, due the
massive improvement in the quality and volume of cosmological measurements. Most of the
N-body simulations rely on a perturbative approach. They use an FLRW cosmological background
(perfectly homogeneous and isotropic)), and assume that any sub-horizon inhomogeneous structure
of the universe will contribute to an average expansion on horizon-sized volumes driven by the
horizon-averaged density. With the next generation of cosmological probes, these may not be accurate
enough to model the universe realistically. There are several ongoing efforts to build fully relativistic,
nonlinear, inhomogeneous and asymmetric models using numerical methods [318320]. This is a
significant contribution to the study of the backreaction effect and the question of the expansion rate of
the universe.
11.3.2. The Primordial Perturbations
In the very early universe, the tiny initial perturbations are thought to have formed a Gaussian
random field whose covariance function is diagonal and nearly scale-invariant. This form, known as
the Harrison-Peebles-Zeldovich spectrum, was assumed in most cases within the Standard Model,
as it corresponds very closely to the observed power spectrum in the universe.
This type of spectrum was first proposed in the 1970s by Edward Robert Harrison [321],
Yakov Zeldovich [322], and Phillip James Edwin Peebles [323], who were working independently,
as the spectrum for initial density fluctuations. This hypothesis was subsequently closely borne out
by observations. The defining characteristic of a Harrison-Peebles-Zeldovich spectrum spectrum is
that it describes a fractal metric, where the degree of perturbation is the same on all scales (hence the
term scale-invariant), so that P(k) k. If we assume scale invariance for the power spectrum on
large scales, and combine this with Equation (101), this implies the following general shape for the
matter power spectrum in the Einstein-de Sitter scenario:
(
P(k)
for
k 3
for
k k0
k k0 .
(117)
42 of 82
The actual form of the spectrum depends in non-trivial ways on the parameters in the cosmological
model, including the slope of the initial power spectrum ns , where P(k ) kns . In a scale-invariant
spectrum in the linear regime, the fiducial value of ns is taken to be 1.
Cosmological observations were consistent with scale invariance up until the early data
releases from the Wilkinson Microwave Anisotropy Probe (WMAP) [214,324], but showed some tension.
Subsequent observations cast further doubts on scale invariance [216]. In 2013, the Planck probe,
mapping the anisotropies in the CMB, led to an important and conclusive result: it ruled out
scale invariance at over 5. The primordial power spectrum was found to be scale-dependent,
with ns = 0.9603 0.0073 [221]. This was confirmed by the 2015 data release, which found that
ns = 0.968 0.006 [205]. This significant result is a powerful demonstration of the importance of
multiple sources of data, joint observations from different probes, and the increasing reliance on
complex statistical techniques for cosmological model selection [325328].
Thus it can be seen that the current Concordance Model, consisting of CDM with an inflationary
epoch in the early universe, was built in stages over the last 100 years. It is the result of a process of
accumulation of evidence and testing of competing models and theories. At each step, theory provides
the basis for adjustments to the model, and observations from different probes provided the evidence
(Figure 4). The current model should in no way be seen as true. It is merely the best model that fits all
the data available so far. Future data may very well require an adjustment to the Concordance Model.
Theory
Cosmological model
General
Relativity
Principle
Big Bang Theory
(FRW)
Observations
Expansion of
the Universe
Copernican
Principle
Standard
Model of
Particle
Physics
Hot Big Bang

Theory
(FRW)
CDM Universe
Element
abundances
Galaxy rotation
curves and bound
clusters
Quantum
Theory
Inflation
CDM Universe
+
Inflation
CDM Universe
+
Inflation
CMB
observations
Accelerating
expansion
Flat Universe
Figure 4. How the Concordance Model of Cosmology was developed. Theories and observations
motivated the development of cosmological models, which were adjusted as new observations
challenged the older models.
12. How Do We Test General Relativity?

The assumption of metric coupling given by Equation (5) has been tested accurately many
times over the last 100 years, at scales from 104 m in laboratories [329,330], up to 1014 m in the
Solar System [331]. These experiment test the implications of the metric coupling: the spacetime
43 of 82
independence of non-gravitational constants, the isotropy of the coupling of all matter field to a unique
metric tensor, the universality of free-fall, and gravitational redshift.
The standard way to test constraints on Sgravity and therefore the dynamics of General Relativity
is by using the parameterised post-Newtonian (PPN) formalism [65]. This assumes that gravity is
described by a metric over all scales. The idea is to write the most general form that g can take in
the presence of matter, when considering correction of order 1/c2 with respect to the Newtonian limit.
This method was first used in 1923 by Arthur Eddington [332]. In its simplest form, it provides
us with two phenomenological parameters PPN and PPN entering the Schwarzchild metric in
isotropic coordinates:
g00
2Gm
= 1 + 2 + 2PPN
rc
2Gm
rc2
2

and
gij =
2Gm
1 + 2PPN 2
rc

ij .
(118)
According to General Relativity, PPN = PPN = 1. The experimental constraints on these

parameters are summarised in Table 3.
Table 3. Some constraints on parameterised post-Newtonian (PPN) parameters from recent tests.
Method
Constraint
Shift of perihelion of Mercury
|2PPN PPN 1| < 3 103
Lunar laser ranging

Very long baseline interferometry
Time-delay variation
Planetary perihelion precessions
|4PPN
|PPN
PPN
PPN
Experiment
PPN
1| =
1 =
3| =
4 104
(4.4 4.5) 104
(1.2 2.3) 105
(2 3) 105
1 =
PPN 1 = (4 6) 105
Data to 1990 [333]

Data to 2004 [334,335]
Data from 1979 to 1999 [336]
Cassini spacecraft [67]
Solar System
ephemerides to 2013 [337]
It should be noted that the PPN formalism assumes that there is no characteristic length
scale for the gravitational interaction, and therefore it does not allow testing of finite-range effects.
These too have been constrained to be very close their General Relativistic value of zero. The deviation
of amplitude from a Newton potential on a characteristic scale is typically < 102 on scales
ranging from a few millimetres to Solar System size [329]. This implies no deviation from GR over
more than 15 orders of magnitude in length scale.
The Solar System tests constrain Eddingtons two PPN parameters to a tiny region very close
to 1. The formalism has been generalized to include eight additional phenomenological parameters [65]
to describe any possible deviation from GR at the first post-Newtonian order. They have all been
constrained to be very close to their General Relativistic values. The latest data use Solar System
ephemerides, which include perihelion measurements for Mercury and the other planets, as well as
lunar ephemerides. Independent teams of astronomers have estimated corrections to the standard
first-order post-Newtonian General Relativistic formalism to be all statistically compatible with
zero [337,338]. This leads us to conclude that General Relativity is the only theory consistent with Solar
System experiments at the post-Newtonian order.
It should be noted that, so far, it is only the first-order post-Newtonian, static, Schwarzschild-like
part of the spacetime metric that has been modelled and tested using Solar System dynamics.
The first-order post-Newtonian gravitomagnetic or Lense-Thirring part of the spacetime metric has
neither been modelled nor tested yet in the Solar System. However, it has recently been pointed out
that this could be possible in the next few years by focussing on particular models (Sun and planets,
planets and spacecraft, or planets and planets) [339]. This would open the field to the possibility of
constraining the PPN parameters using a more complete model of General Relativistic effects.
What about larger or smaller scales? The size of the universe was around 1035 m in the
beginning, and the present size is around 1026 m. GR remains untested at these extreme scales.
44 of 82
But one fundamental assumption of GR is that it describes gravity at all scales (i.e., the theory assumes
the analytical continuity of solutions).
The evolution of the universe itself therefore provides a useful test of General Relativity.
The challenge lies in testing a wide range of potentials (weak and strong field regimes) over
cosmological scales. The strong field regime is tested using compact objects. However, over
cosmological scales, the kind of objects that would produce gravitational potentials approaching
the strong field regime simply do not exist. One solution is to observe the evolution of the
universe and check whether the evolution of large-scale structure corresponds to the predictions
of General Relativity.
13. Cosmological Tests
General Relativity has been submitted to 100 years of Solar System tests, which it has passed with
flying colours. The discoveries of the last two decades: cosmic acceleration, the scale dependence of
the primordial power spectrum, and also the open question of Dark Matter, have made it clear that GR
must also be tested at astrophysical and cosmological scales.
This presents a serious challenge. We derive our knowledge of the universe from measurements
of distances and times, and statistical properties like the distribution of matter. The main difficulty in
extending the Solar System tests to cosmological scales is that these measurements depends strongly
on the construction of cosmological models.
These models depend on four hypotheses:
(1)
(2)
(3)
(4)
A theory of gravitational interactions.

A description of the matter in the universe and the non-gravitational interactions such as
electromagnetic emissions.
A hypothesis on the symmetry.
A hypothesis on the topology, or the global structure of the universe.
Some of the hypotheses are hard to verify, and some have unverifiable implications. Assuming
the symmetry of our solutions means that we also assume the laws of physics are the same throughout
the universe, including its unobservable part outside the cosmological horizon (delimited by a radius
of around 15.7 Gpc). This is a very strong assumption, but one that is unverifiable.
Any cosmological model needs all four hypotheses. The first two are the physical theories, but
their equations cannot be solved without some kind of assumption on the symmetry of the solutions,
given by the third hypothesis. The fourth hypothesis is then an assumption on the global properties of
these solutions.
The simplest CDM Concordance Model assumes that gravity is described by General Relativity
(Hypothesis 1), and that the universe contains the particles and fields of the Standard Model of particle
physics [171], together with Dark Matter and a cosmological constant (Hypothesis 2). The Einstein
equations require an effective stress-energy tensor averaged out on large scales, so the model requires
an extra assumption on the averaging procedure. This averaging and the validity or otherwise of this
assumption is the subject of research on the backreaction effect. The model assumes the Copernican
Principle (Hypothesis 3), and the continuity of the solutions of Einstein equations across all spatial
sections (Hypothesis 4).
The number of alternative cosmological models proposed, especially since the 1990s, is too
numerous to list here. Many of them are impossible to rule out simply by fitting them to the data,
because their predictions are so close to those of CDM. It is easier to test the four hypotheses rather
than trying to test the observables predicted by the models. As with any null test, a significant violation
would indicate the need to modify one or more hypotheses.
As far as tests of GR (the theory of gravity in CDM) are concerned, many experiments test more
than one assumption. We list the main experiment types, and the assumptions they test, in Table 4.
45 of 82
Table 4. What we are actually testing. Experimental tests of General Relativity (GR) often probe more
than one assumption, and isolating the effects is a challenge in itself.
Experiment
Assumption Tested
Solar System tests

Quadrupolar shift of nuclear energy levels
Lunar laser ranging and orbiting gyroscopes
Space-borne clocks
Shift in perihelion of planets
Time invariance of physical constants
Detection of gravitational waves
Metric coupling
Isotropy [340342]
Universality of freefall [343] and structure of metric
Gravitational redshift [344]
Structure of metric [345]
Metric coupling [346]
Lorenz gauge condition and inhomogeneous wave equation
13.1. Testing the Description of Matter and Non-Gravitational Interactions

General Relativity can be tested by measuring the fundamental constants of the theory [346349].
Any variation would require a modification of GR [350,351]. A local measurement of a fundamental
constant, such as a determination of the fine-structure constant from the Oklo phenomenon [352354],
is actually a cosmological-scale measurement along the time dimension. Astrophysical probes such as
21-cm radiation [355] or the cosmic microwave background [356], can be used to test the constancy of
the fine-structure constant , or to constrain simultaneous variations of and Newtons gravitational
constant G [357]. A non-constant G would have serious implications both for Newtonian physics and
for General Relativity. This motivates ongoing efforts to devise new methods to measure this quantity,
and to push the limits of experimental accuracy in order to test the constancy of G [358362].
13.2. Testing the Assumption of Symmetry
The assumption of symmetry can be tested by checking for any deviations from isotropy.
This requires statistical ensembles of data, so the best observables are the CMB, and, on a smaller scale,
large-scale structure [363365].
13.3. Testing the Gravitational Interactions
The gravitational interactions have been tested many times at laboratory scales and all the way
up to Solar System scales. The challenge today is to test them at cosmological scales.
Tests of General Relativity may be placed in three broad categories: laboratory, astrophysical, and
cosmological. Experiments cannot span all length scales, and so they cannot test all theories. We are
forced to design experiments which can test alternative theories, or the effects of GR, at particular
length scales (Table 5).
Laboratory tests probe effects from sub-millimetre scales to a few hundred metres. Galileos
experiments on weights dropped from the Leaning Tower of Pisa [366] are one example of a
laboratory-scale test. The torsion balance experiments of Cavendish and Etvs [367], and their
modern versions [368371], are another.
Astrophysical tests probe gravity from Solar System scales all the way up to galactic cluster scales.
They include experiments such as laser ranging off the Moon (lunar laser ranging) and several of
the Earths artificial satellites [372375], and proposals to extend laser ranging to other planets in
the Solar System (planetary laser ranging) [376380], radar astrometry on near-Earth objects [381],
observations of the precession of the perihelion as well as higher-order effects for Mercury and
other planets [338,382,383], observations of spacecraft in the Solar System [337], actual or proposed
measurements of the acceleration of the Pioneer [384386] and New Horizons probes [387389],
observations of the orbits of compact objects such as neutron stars [390392], and measurements of
galaxy rotation curves [393,394]. In the last decade, observations of extrasolar planetary systems have
been proposed and used in order to test GR [395404]. All of these systems are characterised by their
spherical symmetry. Many can be approximated by a test object orbiting at a distance r around a
46 of 82
central mass M. The gravitational field of the more massive central object is then probed by observing
the orbit of the test body.
Table 5. Bridging the length scales to test the cosmological model. Experimental tests of gravity and
dark sector couplings, at their typical length scales. Massive gravity (MG) screening mechanisms
would show up at short ranges, while smooth dark energy manifests itself at cosmological scales. We
give the experimental accuracy from current and future experiments planned over the next decade.
The comparison between growth and expansion history comes from combined BAO, supernova, weak
lensing, redshift distortion, cosmic microwave background (CMB) lensing, and cluster data. Lensing
effects and dynamical mass comparisons can be carried out over a range of scales: inside galaxies,
using strong lensing and stellar velocities, and on cosmological scales, using cross-correlations. This is
a route to testing screening effects in alternative theories. Laboratory and Solar System tests can also
probe dark sector couplings besides short-range effects, but many of the constraints obtained depend
on the cosmological model.
Test
Length Scale
Theories Probed
Current Status (and Future)
Growth vs. expansion history

Lensing vs. Dynamical mass
Astrophysical tests
100 Mpc 1 Gpc

0.01 100 Mpc
0.01 AU 1 Mpc
GR with smooth dark energy

Test of GR
MG screening mechanisms
Laboratory and Solar System tests
1 mm 1 AU
PPN parameters in MG
10% accuracy (2%4%)

20% accuracy (5%)
10% (Up to 10 times improvement)
Constraints are model dependent
(Up to a tenfold improvement)
For such systems, the deviation of the metric from the Minkowski form is characterised by the
magnitude of the Newtonian gravitational potential
e=
GM
.
rc2
(119)
The strongest gravitational fields accessible to an observer occur in the limit e O(1). In this
limit, the central object is a black hole and the test object is close to the event horizon.
We know that the Riemann curvature tensor is an essential quantity in General Relativity.
The approximate magnitude of this tensor is expressed by the Kretschmann scalar (the fully contracted
Ricci scalar) for the Schwarzchild metric:
GM
3 2.
(120)
r c
Note that this expression is more complicated for rotating objects [405]. For the purposes of our
description, however, we may use this simple expression, which we call the curvature.
The third broad category is the cosmological tests. First, the measurements are made on a
statistical ensemble. The masses must be treated as power spectra. Second, the gravitational field
assigned to each wavenumber k is very weak. The description of the system often depends on the
cosmological model. At this scale, the description must take into account the background perturbations
and the expansion of the universe.
The position of the various systems on a potential-curvature parameter space is shown in Figure 5.
The potential accessible to observers is bounded by the line e ' 1. There is no limit to the maximum
curvature that can be observed except for the Planck limit, where r ' 1.6 1033 cm. This lies many
orders of magnitude above the boundaries of the figure. General Relativity is a complete theory, and
does not fail below the Planck scale. To test or even probe GR below the Planck scale is a different
matter, and a goal that has not yet been achieved. Many alternative theories of gravity that attempt to
explain cosmic acceleration do so through modifications to GR at the Planck scale, so it is important to
explore ways in which GR can be tested at these scales.
47 of 82
10 10
Black holes
Neutron star
Gravity Probe B
10 30
10 40
R
1 AU
10 50
10 60
Moon
Voyager 1
Hulse-Taylor
Mercury
New Horizons
Pioneer 10 and 11
Dark matter
10 kpc
10 Mpc
Supermassive BH
galactic centre
Beyond the horizon
Curvature GM/r 3 c 2 [cm 2 ]
10 20
Dark energy
10 1510 1410 1310 1210 1110 10 10 9 10 8 10 7 10 6 10 5 10 4 10 3 10 2 10 1 1

Gravitational potential GM/rc 2
Figure 5. The parameter space for quantifying the strength of a gravitational field. The horizontal axis
measures the potential. The vertical axis measures the spacetime curvature of the gravitational field at
a radius r away from a central object of mass M. Regions of this parameter space with potential greater
than 1 represent distances from a gravitating object that are smaller than the event horizon radius and
are therefore inaccessible to observers. The red vertical line on the right-hand side of the plot marks
the horizon limit. This is a schematic plot, and in no way do we show an exhaustive list of objects
and systems that have been used or could be used to test GR. The region of Solar System-scale tests is
broadly bounded by the Moon, Gravity probe B [406,407], Mercury, and the Pioneer and New Horizons
spacecraft. We have included the Voyager spacecraft in the diagram, even though, unlike the Pioneer
probe, it was never suitable for tests of GR. The famous Hulse-Taylor binary pulsar [408], although a
neutron star binary, is also roughly in the region of Solar System tests. Black holes and neutron stars
are in the strong field regime, and the former are at at the limit of the event horizon boundary. Adapted
from [409].
What is the minimum curvature? The unperturbed FLRW metric is isotropic, and the unperturbed
Kretschmann scalar is a function of time only. The curvature for the homogeneous universe drops as the
universe expands. The present universe has a curvature which is just above the boundary of the region
marked Dark energy, since dark energy is not yet completely dominant over pressureless matter.
However, the curvature will approach this limit asymptotically. In this paradigm, this represents a
fundamental minimum curvature scale. It is shown in the figure by the region labelled Dark energy.
Galaxies are astrophysical probes of GR. Their innermost regions can be modelled as test particles
orbiting a central supermassive black hole. Galactic rotation curves, which can only be explained by
introducing Dark Matter, describe galaxy velocities up to the outermost regions. Systems below a
constant acceleration scale of approximately 1.2 1010 ms2 cannot be modelled without adding
a contribution to the gravitational filed in the form of Dark Matter. This constant acceleration is the
diagonal boundary of the region labelled Dark matter in the figure.
48 of 82
Note that the regions of parameter space occupied by Dark Matter and dark energy overlapthere
is a degeneracy between these two unknown components of the cosmological model. However, it is
not impossible to distinguish between their effects, since their properties are different. In particular,
Dark Matter forms clumps just like baryonic matter, while dark energy does not.
GR has not been tested in the region between curvatures of 1040 and 1050 . This corresponds
to the region between Solar System scales and cosmological scales probed by galaxy surveys and the
CMB. The challenge lies in finding systems which span these scales. In theory they do exist, in the form
of galaxies, clusters and superclusters. Their rotation curves transition from Schwarzchild orbits in
their innermost regions, to outer regions dominated by Dark Matter. However, our observations are
hampered by the fact that we are limited by the resolution of our telescopes, so we can only observe
the outer regions. In addition, the untested region is situated between a region where GR is extremely
well-constrained (by Solar System tests), and a region where Dark Matter and dark energy have
to be invoked, and where we must take the cosmological model into account, because the effect of
large-scale structure on background dynamics becomes non-negligible. It has been suggested that this
backreaction effect may be at the origin of the observed cosmic acceleration.
Cosmological tests of General Relativity may be broadly classified as follows:
Tests of the consistency between the expansion history and the growth of structure. A discrepancy
in the equation of state parameter of dark energy w, inferred from the two approaches can indicate
a breakdown of the GR-based smooth dark energy cosmological paradigm.
Detailed measurements of the linear growth factor across different scales and redshifts.
Comparison of the cosmological mass distribution inferred from different probes, especially
redshift space distortions and lensing.
14. Possible Modifications of GR and Cosmological Implications

It is useful to identify the regimes in which modifications to GR may appear. This enables us to
get a clearer picture of the capabilities and limitations of current and future experiments to tests these
alternative theories.
14.1. Weak and Strong-Field Regimes
In order to test gravity in the strong-field regime, we need to observe compact objects with a very
high density [410412]. Black holes are good candidates. Their compactness and mass takes GM/rc2
close to the maximum limit of unity. However, they have a serious drawback. Because of the no hair
theorem, they are not characterised by any coupling to a scalar field and therefore cannot be used to
test for this effect, and to discriminate between scalar-tensor theories and GR.
Neutron stars, on the other hand, are still very compact bodies, but they can be strongly
coupled to a hypothetical scalar field. This property has been used to test relativistic parameters by
observing the Hulse-Taylor binary pulsar PSR B1913+16 [413], and the neutron starwhite dwarf binary
PSR J1141-6545 [414]. Gravitational time delay has been tested using other binary pulsars [415419], and
pulsars have also been suggested as ideal probes of the Lense-Thirring effect [420,421]. General
Relativity passes these tests with flying colours.
14.2. Small and Large Distances
Distance-dependent modifications can be induced by a massive degree of freedom, which will
cause a Yukawa-like coupling. General Relativity is very well constrained on the size of the Solar
System, and there are several tests constraining Yukawa interactions at Solar System scales [422424],
but there are no constraints on scales larger than 10 h Mpc1 , at least without assuming some
cosmological model. Some theories put forward to explain cosmic acceleration, such as Chameleon
mechanisms [425], are essentially modifications of GR at cosmological distance scales.
49 of 82
14.3. Low and High Accelerations

Galaxy rotation curves and galaxy dynamics motivated the Dark Matter paradigm. The Tully-Fischer
law [426] tells us that Dark Matter cannot be explained by a modification of General Relativity at a
fixed distance. MOND instead explains it by modifying gravity at low accelerations below the typical
acceleration a0 108 cm s2 .
14.4. Low and High Curvature
Curvature R is important in distinguishing possible extensions of the Einstein-Hilbert action.
A curvature-dependent may become important even if the potential remains small. In the Solar
System, the curvature R 4 1028 cm2 . The curvature of the homogeneous universe according to
the Friedmann equation is

RFLRW (z) = 3H02 m (1 + z)3 + 4 ,
(121)
from which we can see that the curvature of the universe evolved with time, from 1033 cm2
at the time of nucleosynthesis, to 1056 cm2 at z = 1. The curvature scale associated with the
cosmological constant is R = (1/6), so the phenomenology of the cosmological constant occurs in
low curvature regime
R < R 1.2 1030 R .
(122)
This is of particular interest to paradigms which seek to explain cosmic acceleration through the
backreaction effect. For cosmological-scale perturbations, we are always in the weak field regime.
However, the curvature perturbation associated with large-scale structure becomes of the order of the
background curvature at redshift z0, even if we are still in the weak field limit. This means that the
effect of large-scale structure on the background dynamics may be non-negligible [427].
In summary, in order to explain the dark energy or Dark Matter problem by modifying General
Relativity, the modifications have to be either at large scales (typically Hubble scales), low accelerations
(typically below a0 , or small curvatures (typically R ). The regions corresponding to dark energy and
Dark Matter curvature and potentials are shown in Figure 5.
14.5. Cosmological Probes
The idea of testing General Relativity using large-scale structure was first proposed in 2001 [428].
This relies on the ingenious idea that if gravity is well-described by General Relativity, and the universe
well-described by CDM, then on sub-Hubble scales, and considering only scalar perturbations,
the spacetime metric can be written as
ds2 = (1 + 2) dt2 + (1 2) a2 (t)ij dxi dx j
(123)
where and are the two potentials, and ij is the metric of the spatial section [429431]. The Einstein
equations reduce to the Poisson equation
4 = 4Gmatter a2 matter
(124)
= 0,
(125)
and
since the matter anisotropic stress is negligible. The spectrum of the two gravitational potentials has to
be proportional to the matter power spectrum. The scale dependence of the gravitational potential P
and of the matter distribution P(k ) are related by:
k4 P (k, a) =
9
m H 2 a2 P(k, a) .
4
(126)
50 of 82
If the Poisson equation is modified by some modification of gravity, the matter power
spectrum changes shape according to the cosmological model assumed, so the above relation is not
model-independent. However, the fact that the two spectra differ is independent of the cosmological
model. It therefore provides a test of the underlying gravitational theory. Such a test can be carried out
by comparing weak lensing data to galaxy surveys.
Various similar approaches have been proposed [432]. This approach allows us to test
various classes of alternative theories by means of large-scale structure. Some of these theories
include Dvali-Gabadadze-Porrati models [134,433,434], quintessence [435439], and scalar-tensor
theories [440443]. The difficulty lies in finding a parameterisation of the perturbation equations which
is consistent with the one used for the background evolution, since both assume the same theory for
gravitational interactions.
The cosmic microwave background is another potential testing ground for General Relativity.
The amplitude and position of the peaks of the cosmic microwave background allows us to probe
the potential wells present during the recombination era by extracting information on the primordial
power spectrum created during inflation. Planck data, alone or in combination with weak lensing has
been used to test GR [444,445] and modified gravity [446448].
Despite the fact that GR is extremely well tested on laboratory and Solar System scales, cosmology
provides plenty of scope for alternatives to General Relativity. Let us consider cosmic acceleration,
which is now a confirmed observation. It requires an explanation. The cosmological constant paradigm
rests on three assumptions: that the observations are correct, that GR is the correct theory of gravity,
and that our FLRW model of the universe is correct.
More than two decades of precision measurements have removed any doubt that the observed
acceleration may be due to incorrect modelling of experimental errors. So in order to do away with the
cosmological constant, we must seek possible answers in the other two assumptions.
Backreaction is an alternative to the alternative theories, or to the Dark Matter paradigm. It is
firmly within the General Relativity paradigm, and seeks to explain cosmic acceleration by modelling
the universe and the structures within it, and therefore its expansion history, in more detail than is the
case with the FLRW model. It keeps the first and last assumption, but does away with the assumption
of an FLRW universe. However, is it enough to explain the observed acceleration?
If GR is assumed to be the correct theory together with the other two assumptions, cosmological
measurements are usually interpreted as providing evidence for Dark Matter and a nonzero
cosmological constant or dark energy. This poses conceptual problems. Why is the observed value
of the cosmological constant so small in Planck units? It also poses a coincidence problem. Why is
the energy density of the cosmological constant so close to the present matter density? No dynamical
solution of the cosmological constant problem is possible within GRthe cosmological constant is not
the attractor of some dynamical function.
This opens the field to possible modifications of GR. Should GR be modified at low and high
energies? This is a serious challenge for theorists. Einsteins theory is the unique interacting theory of
a Lorentz-invariant massless spin-2 particle. New physics in the gravitational sector must introduce
additional degrees of freedom. These additional degrees of freedom must modify the theory at low
or high energies, or both, while being consistent with GR in the intermediate-energy regime, that is,
at length scales 1 6 ` 6 1011 m, where the theory is extremely well tested.
Figure 6 illustrates some of the difficulties in testing the completeness of General Relativity.
There is an obvious scale desert between Solar System scales, and cosmological scales. At the other
end of the scale, there is a gap which is often overlooked: between sub-millimetres scales, which is the
current limit where GR has been tested, and the Planck scale. A modification of gravity at very small
scales would be apparent in this regime. Even ignoring quantum effects, it is a serious challenge to test
gravity in a regime in which vastly stronger forces come into play.
51 of 82
10 20
10 30
10 40
10 50
Beyond the quantum limit
Potential GM/rc 2 [metres 2 ]
10 10
Limit of short-range tests
Neutron star
Sun Milky Way
Earth
Limit of
homogeneity scale
Last scattering surface
Proton
Beyond the Hubble patch
Black hole
Perturbation regime
43
10
Planck scale
10
30
10
15
10 15
10 27
Observable Universe
Scale size [metres]

Figure 6. The parameter space for experiments. The horizontal axis is the typical length scale of
the object in question. The vertical axis measures the gravitational potential. The red vertical line
on the left-hand side marks the Planck scale. The vertical line on right-hand side of the plot marks
radius of the observable universe, or Hubble radius. Experimental verification of GR is impossible
beyond these limits. The radius of the surface of last scattering is only slightly smaller than our Hubble
radius. Assuming GR implies assumptions far beyond the range that has been experimentally tests.
For instance, if we define the Planck mass as mPlanck = h c/G, we are assuming that the gravitational
constant remains constant down to the Planck length. This extrapolates the inverse square law over
a scale of more than 1030 from what ha s been tested. The green region is where Solar System tests
have been carried out. Beyond 100 Mpc, assumption of a homogeneous and isotropic metric becomes
accurate enough to use in physical models. The blue region shows length scales at which the FLRW
metric is valid. By way of comparison to this parameter space, the nonlinear regime for perturbation
theory, which gives use the matter power spectrum, covers 1022 to 1023 m, the linear, quasi-static
regime covers 1023 to 1025 m, and beyond that is the superhorizon regime. These length scales
ranging from 1 Mpc to above 1 Gpc fit on a tiny part of the horizontal axis above, shown by the thick
blue horizontal line.
15. The Nature of Dark Energy and the Implications for General Relativity
The Concordance Model of cosmology assumes that General Relativity is correct, an assumption
which is justified by the tests which GR has undergone. Within this model, 95% of the content
of the universe is unaccounted for. Dark matter, which makes up around 25% of the mass-energy
of the universe is a matter-like component which is cold (sub-relativistic) and weakly interacting.
The discrepancy between the observed acceleration of the expansion of the universe and the predictions
of GR leads to the conclusion that there must be a cosmological component with a negative equation
of state parameter making up around 70% of the mass-energy content of the universe: dark energy.
But the dark energy paradigm does not fix the nature of this component. There exist many theories
which attempt to explain its nature. In particular, one can ask whether the basic assumptions of the
Concordance Modelhomogeneity and isotropyare correct.
52 of 82
If the universe is not homogeneous and isotropic, the FLRW equations are no longer valid.
Over the last decade, one line of research has attempted to explain the accelerated expansion by
exploring the implications of an inhomogeneous universe on a General Relativistic cosmology. This
effect is the backreaction. This approach does away with dark energy.
There are various ongoing investigations on the effect of the backreaction due to an
inhomogeneous universe [269,427,449455], with different lines of research offering different
interpretations of the Buchert equations, where the Friedmann equations are supplemented by an
additional backreaction term [456]. Whether one can explain all of the observed expansion history of
the universe as a consequence of the growth of inhomogeneities without invoking some additional
fluid component is the subject of ongoing debate [452,457]. The backreaction, even if it turns out to be
incapable of replacing the dark energy paradigm, is still a subject worth investigating, if anything as a
correction to the homogeneity assumption, which obviously breaks down at small scales.
If we assume that the Copernican Principle holds, then the universe is well described by a
Friedmann-Robertson-Walker spacetime. The dynamics of the background expansion are determined
by the content of the universe: the list of fluids (perfect fluids, due to the assumption of the Copernican
Principle), with their equations of state. Within the dark energy paradigm, we can distinguish two
main strategies for formulating hypotheses:
(1)
(2)
there is some new kind of component in the universe, or

there is some new property of gravity.
Let us first recall that General Relativity rests on two assumptions: the gravitational interaction is
described by a massless spin-2 field, and matter is minimally coupled to the metric, which implies the
weak equivalence principle. The Einstein-Hilbert action described by Equation (6) implies that
Sgravity =
c3
16G
g d4 x + Smatter [matter; g ] ,
(127)
where R is the Ricci scalar of the metric tensor g , and Smatter [matter; g ] is the action of the matter
fields. If we only consider field theories, this gives us a useful classification scheme for the different
theories that seek to explain the nature of dark energy.
The first strategy listed above assumes that gravitation is described by General Relativity, and
introduces new forms of gravitating components beyond the Standard Model of particle physics.
This adds a new term SDE [; g ] to the action in Equation (127) while keeping the Einstein-Hilbert
action of all the standard and Dark Matter unchanged.
The second strategy modifies gravity, and therefore extends the action, either by modifying the
Einstein-Hilbert action of the coupling of matter. These theories also involve new forms of matter.
A cosmological constant is the simplest modification which can be made to gravity, and it is
equivalent to dark energy with a constant equation of state. To explain the observed acceleration,
the new form of matter must have an equation of state w < 1/3. Dark energy can also be attributed
to the energy of the vacuum, although the energy predicted by the Standard Model of particle physics
is either 0 (using super-symmetry), or 10120 orders of magnitude larger than the observed cosmological
value [458]. There are ongoing attempts to solve this fine-tuning problem using string theory [459461],
causal sets [462,463], or by using anthropic arguments [464,465] .
The other approach is to attribute dark energy to a scalar field whose potential has evolved in
some way that it currently exerts a negative pressure. Such fields, in theories within the framework
of GR, are termed Quintessence. Their distinguishing feature is that they allow the equation of state
of dark energy to evolve. Alternatives to Quintessence within the same approach include K-essence,
Phantom Fields, or the Chaplygin Gas [466469].
Another strategy is to depart from General Relativity and modify the laws of gravity and posit
dark energy as the manifestation of an effect arising from extra dimensions, or higher-order corrections.
Within this category, the more successful theories have been of two types. Dvali-Gabadadze-Porrati
53 of 82
(DGP) dark energy considers the universe as a 4D brane within a 5D Minkowskian bulk [134].
The weakness of gravity relative to the other forces is explained by gravity leaking into the higher
dimensions as it acts through the bulk [470], whereas the other forces act within the brane. The other
class of theories is f ( R) gravity, where the Ricci scalar R in the Lagrangian is replaced by some function
f ( R). Such theories correspond to scalar-tensor gravity with vanishing Brans-Dicke parameter [471].
We can therefore identify four main classes of theories, as shown in Figure 7. Classes 1 and 2
assume GR and introduce new forms of gravitating matter, the difference being that in class 2, the
distance-duality relationship may be violated due to mechanisms such as photon decay. Classes 3 and
4 modify gravity in some way. In class 3, a new field introduces a long-range force so that gravity
is no longer described by a massless spin-2 graviton. In this class, there may be a variation of the
fundamental constants. Class 4 includes theories in which there may exist an infinite number of
new degrees of freedom, such as brane models or multidimensional models. These models predict
a violation of the Poisson equation on large scales. The constraints which may be placed on these
theories, and the experiments to test them, correspond to the regimes in which gravity is modified, as
summarised in Section 14.
An exhaustive list of dark energy or cosmic acceleration theories is beyond the scope of this
review. We simply note that most of them give different predictions for the equation of state of dark
energy, its evolution, or the expansion history of the universe . To distinguish between these proposals
we need to track the evolution of these parameters throughout the history of the universe.
Gravity: Assume GR
1
Gravity: Modify GR
3
SDE [DE; g ]
New fields dominate

the matter content of
the Universe at low
redshift
e.g. quintessence
2
SField [Field; g ] ! SNew Field [New Field; g ]
e.g. photon-axion mixing,Chameleon models
SMatter [Matter; g ] ! SMatter [Matteri ; A2i ( )g ]
e.g. scalar-tensor theories
Extra dimensions
e.g. brane-induced gravity, multigravity
Figure 7. The four main classes of dark energy theories, within the two broad strategies, classified
as modifications of the General Relativistic action. Classes 1 and 2 assume the gravitational metric
coupling of GR, whereas classes 3 and 4 modify this metric coupling, and are therefore modifications
of gravity. In the upper line of classes (1 and 3), new fields dominate the matter content of the recent
universe. Adapted from [431].
16. The Current Status of General Relativity

General Relativity has been subjected to a multitude of tests in its 100 years of existence. As of 2016,
the main predictions of GR have been tested and confirmed. Whereas it is sufficient for most purposes at
ordinary accelerations and energy scales to use Newtonian calculations, General Relativity has found its
way into daily life, in the Global Positioning System and geodesy. The postulates of General Relativity
have been confirmed with ever-increasing accuracy. Deviations from the Einstein Equivalence Principle
are now constrained to below 1014 , deviations from Local Lorentz Invariance down to 1020 , and
deviations from Local Position Invariance down to 106 .
General Relativity has been probed down to scales of 106 metres in laboratories, and up to
1000 AU in space-based experiments and observations. Astrophysical and cosmological observations
54 of 82
have probed GR at scales from 1 Mpc up to gigaparsec scales. In the Solar System, the dynamics of GR
have been tested with radar and laser-ranging. We can track the ephemerides of the planets in the Solar
System, right up to the minor outlying planets (such as Sedna, with an aphelion just short of 1000 AU).
The farthest objects whose trajectory has been followed from the moment they were thrown are
now outside the Solar System. They are the Pioneer 10 and 11 spacecraft, and are now around 70 AU
distant, and Voyager 1, which is at a distance of 135 AU, making it the only object to have reached
interstellar space. Communication with the Pioneer spacecraft ceased in 2003. The Pioneer spacecraft
exhibited an anomalous constant acceleration towards the Sun which could not be explained using GR.
This prompted a reexamination of all the recorded data. It is now generally accepted that the anomaly
is caused by thermal radiation, and that once this is accounted for, there is no remaining anomalous
acceleration [472].
Gravitational waves, among the last untested predictions of GR, were first detected in late 2015,
ushering a whole new era of observational astrophysics, in which the strong field regime can be probed
and where the predictions of GR can be tested.
The evidence for General Relativity is extremely strong. The theory has passed all tests in the
weak-field limit at Solar System scales, including tests of the assumptions (the Equivalence Principle)
and the predictions specific to GR (frame-dragging, gravitational time dilation), and in the strong field
with the observation of a binary black hole merger and the resulting gravitational waves. We give a
summary of the experimental milestones in Table 6. However, there are still issues that allow room
for speculation.
The first is the question of the completeness of GR. Is it valid at all scales? There is a scale gap in
our tests of GR between laboratory, Solar System and galactic scales, and cosmological ones. This gives
rise to a multitude of domains of validity for different alternative theories, whereas ideally we should
seek a universal theory that can explain phenomena at all scales.
The second issue is the accuracy of our approximations. General Relativity may be conceptually
simple in that it is based on a minimum number of postulates. But the resulting field equations,
when applied to real physical systems, can be very hard to solve. We get around this by making
approximations such as spherical symmetry, or homogeneous and isotropic perfect fluids, which allows
us to obtain analytical solutions. However, the accuracy of these approximations may not always
be good enough. This is evident in the case of cosmological perturbations and large-scale structure.
At which scales is it valid to use a Friedmann-Lematre-Robertson-Walker metric? Can the backreaction
explain some or all of the observed cosmic acceleration? For smaller systems such as aspherical
collapsing bodies, we still need accurate models in order to match theory and observation. Such
questions have spurred the development of numerical methods in General Relativity. Physicists now
have the necessary computing power to go far beyond simple first-order approximations.
The remaining open questions are of a cosmological nature: Dark Matter and dark energy,
and inflation. Dark matter and dark energy account for around 25% and 70% of the mass-energy
content of the universe, respectively. They are not a problem for the theory of gravity itself, but it
does mean that we do not know the nature of 95% of the content of the universe. Inflation solves
a number of cosmological problems, but whatever theory we choose to explain inflation, we still need
to introduce new physics beyond GR.
55 of 82
Table 6. The major predictions of GR, and their experimental status.

Effect
Milestones
Current, Future and Proposed Experiments
Status
Mass equivalence
Galileo [366,473]
Etvs [367]
Et-Wash Group [371]

Lunar Laser Ranging [343]
MICROSCOPE[474]
STEP [475]
Galileo Galilei satellite [476]
Confirmed
Gravitational time
dilation [478]
Eddington solar eclipse [4]

Pound-Rebka experiment [479]
Space-borne hydrogen masers [481]
Quantum interference of atoms [477]

ACES [480]
Galileo 5 and 6 satellites[344]
Einstein Gravity Explorer [482]
Confirmed
Precession of orbits
Orbit of Mercury (Einstein [3])
Binary pulsar observations [483485]

Solar System and extrasolar planets [337,402]
Confirmed
De Sitter precession
Gravity Probe B [486]

Lunar laser ranging [488490]
Binary pulsars[493,494]
Binary pulsars [487]

Improved lunar laser ranging [491,492]
Confirmed
Lense-Thirring precession
Gravity Probe B [486]

LARES [420,421,496]
LARASE [495]
Laboratory tests [497]
Solar System bodies [339]
Binary pulsars [420,421]
Black holes [478,496,498,499]
Confirmed
Gravitational waves
LIGO [80]
Advanced LIGO [500]

eLISA https://www.elisascience.org
Recently confirmed
in two events [80,501]
Strong field effects
PSR J0348 + 0432 [483]
Black holes [411,502]

Binary pulsars [117]
Recently confirmed
Low-orbit satellites [503506]

Stars orbiting black holes [506]
Juno spacecraft around Jupiter [507]
Not yet observed
Orbital precession due to

oblateness of central body
56 of 82
17. Future Developments

The current Concordance Model of cosmology was built in successive (and sometimes concurrent)
steps. General Relativity applied to a spacetime under the Copernican Principle, filled with
pressureless matter, produced the Einstein-de Sitter model. Motivated by the observed Hubble
expansion, it resulted in the Big Bang model. Clues from element abundances, baryon assymmetry, and
knowledge of nucleosynthesis from the Standard Model of particle physics meant that the universe had
to have a thermal history, which resulted in the Hot Big Bang Model. When evidence for missing mass
became incontrovertible, cold Dark Matter had to be added to the inventory of cosmic components.
This model worked well, but not well enough. It could not explain the observed homogeneity of the
universe across regions which were causally disconnected, nor its flatness. So Inflation was introduced.
Observations of an accelerated cosmic expansion motivated a search for explanations within the then
current paradigm, which resulted in various hypotheses: a curved geometry, supermassive neutrinos,
or perhaps a particular cosmological topology. In the end, the paradigm had to be shifted yet again
with the introduction of dark energy.
The Concordance Model can explain the observations with just six parameters: the physical
baryon density parameter b h2 , where h is the Hubble parameter, the physical Dark Matter density
parameter c h2 , the age of the universe t0 , the scalar spectral index ns , the curvature fluctuation
amplitude 2R , and the reionisation optical depth . That such a degree of fit is offered by such a simple
model is remarkable.
The success of the Concordance Model has been its ability to include physical effects at
extremely different scales, from primordial nucleosynthesis to large-scale structure evolution,
in one coherent theory. However, this does not allow us to state that the CDM model is correct.
It merely implies that deviations from CDM are too small compared to the current observational
uncertainties to be inferred from cosmological data alone. This leaves room for some very fundamental
open questions, which we have described in this review.
The science of cosmology finds itself at a critical point where it has to make sense of the vast
quantity of data that has become available. Different probes have allowed us to piece together
interlocking information which, so far, confirms the Concordance Model. The cosmic microwave
background has provided conclusive evidence of a flat geometry, super-horizon features, the correct
harmonic peaks, adiabatic fluctuations, Gaussian random fields, and most recently, a departure from
scale invariance. We have not yet observed primordial inflationary gravitational waves. Large-scale
structure observations, which probe the recent universe, provide firm evidence in favour of the
Concordance Models explanation of the evolution of density perturbations and the growth of
structure, and provide a bridge between the effects of long-range gravitational interactions and
shorter-range forces.
The recent detection of gravitational waves in two events (possibly three), one hundred years
after they were predicted by Einstein, directly validates General Relativity in several ways. It shows
that GR is correct in the strong-field regime, that black holes really exist, and black hole binaries too,
and it proves that gravitational waves are a real physical phenomenon, and not just a mathematical
artefact of GR.
We are fast approaching the point where cosmological observations will be limited only by cosmic
variance, i.e., no more data will be available from our Hubble patch [508]. What does that imply for
the development of new theories? How do we test the predictivity of these theories without new data?
We will have to look for new effects in old data, and for new correlations in future data.
Large-scale structure is extremely useful in testing General Relativity and cosmological models,
but the future may bring other observational windows. In particular, the next fifty to one
hundred years may see the development of a gravitational wave astronomy. This is likely
to be an even more significant development than even CMB astronomy. The observation of
primordial gravitational waves would provide vital information on the inflationary epoch [509,510].
B-mode polarisation of the CMB offers an indirect pathway to the observation of this gravitational wave
57 of 82
background [511,512]. The second major development could be the observation of a cosmic neutrino
background [513], which is the result of neutrino decoupling in the lepton era. This would push the
observations along our past light cone even further back in redshift, providing information on the
universe before recombination and the CMB.
17.1. Plausible Conclusions from Incomplete Information
The statistical questions facing cosmologists pose some particular problems. We observe a finite
region of our universe, which is itself a single realisation of the cosmological theory. We can only
observe whatever is on or inside our past light cone, as shown in Figure 1. Not only are we limited by
cosmic variance, we also have just a single data point for the cosmological model.
There exist several alternatives to GR and to CDM that have not been ruled out by experiment.
Constructing viable physical models is not just a question of fitting the model to the data. It is a question
of model selection, which requires robust statistical techniques that allow us to make sensible decisions
using our incomplete information. Bayesian inference provides a quantitative framework for plausible
conclusions [514,515]. We can identify three levels of Bayesian inference.
(1)
(2)
(3)
Parameter inference (estimation). We assume that a model M is true, and we select a prior for the
parameters , or the Prob(| M ).
Model comparison. There are several possible models Mi . We find the relative plausibility of each
in the light of the data D, that is we calculate the ratio Prob( D | Mi )/Prob( D | M0 ).
Model averaging. There is no clear evidence for a best model. We find the inference on the
parameters which accounts for the model uncertainty.
At the first level of Bayesian inference, we can estimate the allowed parameter values of the theory.
If we assume General Relativity as our theory, we still need to fix the values of the various constants.
This is the rationale behind ongoing efforts to measure quantities like Newtons constant G ever more
accurately, despite the fact that the parameter has been around for three centuries. What are the energy
densities of the various components of a CDM universe? If we include a dark energy equation of
state parameter w, what value does it take?
Next, we can ask which parameters we should include in the theory. Should we include a
cosmological constant in General Relativity? Or should we include dynamical dark energy parameter?
Although current data are consistent with the six-parameter CDM model based on GR, there are
more than twenty candidate parameters which might be required by future data (see [516]). We cannot
simply include all possible parameters to fit the data, since each one will give rise to degeneracies
that weaken constraints on other parameters, including the CDM parameter set (e.g., [517519]).
The landscape of alternative cosmological models is even larger if we relax our assumptions on the
theory of gravity [520].
The goal in data analysis is usually to decide which parameters need to be included in order
to explain the data. For physicists, those extra parameters must be physically motivated. That is,
we need to know the physical effects to which our data are sensitive, so that we can relate these effects
to physics. At the current state of knowledge, we have to acknowledge the possibility of more than
one model. We therefore require a consistent method to discard or include parameters. This is the
second level of Bayesian inferencemodel selection.
Bayesian model selection penalises models which introduce wasted parameter space. We can
always construct a theory that fits the data perfectly, even better than GR, but we would need to
introduce extra free parameters (e.g., extra fields or couplings between matter and the metric). This is
the mathematical equivalent of Occams razor. We seek a balance between goodness of fit (the degree
of complexity) and predictive power (consistency with prior knowledge). General Relativity fits all the
data with the minimum number of parameters. In cosmology, CDM is the best model because it only
involves one new parameter and no new fields.
The problem of model selection in relation to GR is as old as the theory itself. But only
recently have cosmologists have started to use Bayesian methods for cosmological model
58 of 82
selection (e.g., [521524]), when the astrophysical data began to have the necessary statistical power to
enable model testing. Bayesian techniques are starting to be applied to General Relativity itself [525].
With the next generation of astrophysical probes in the pipeline, model selection is likely to grow in
importance [514].
Bayesian model selection cannot be completely free of assumptions. In cosmology, there is some
model structure which depends on a number of unverifiable hypotheses about the nature of the
universe. The Copernican Principle is one such hypothesis [9,526,527].
The third level of Bayesian inference is model averaging. In the current scenario, there is firm
evidence for General Relativity as the best model for gravitational interactions. However, it is still
useful to quantify our degree of certainty (or doubt), for one simple reason: we do not have the final
list of alternative theories. In other words, when we choose GR against any number of alternative
theories, we have no knowledge about other alternatives outside that list (such as theories yet to be
developed, for example). At best, we know which alternatives are ruled out by the data. This level of
Bayesian inference is the application of a principle that has been called Cromwells Rule [528]: even if
all the data show our theory to be correct, we should allow a non-zero probability, even if tiny, that the
theory is false.
The utility of alternative theories becomes evident when we apply Bayesian model selection.
For science to advance by falsification, it is not enough to claim that the present theory is false. We need
to know which alternative theory is favoured instead. Newtonian gravity would have likely have
survived the 1919 eclipse if General Relativity had not been formulated.
17.2. Experimental Progress
There has been rapid progress in constraining cosmological parameters and models over
the last two decades, with a multitude of experiments observing the CMB, large-scale structure,
galaxies and supernovae. We will just provide a summary of the most recent data sets.
The first are anisotropies in the CMB, where the main statistic is the angular power spectrum
of fluctuations C` , and polarisation of the CMB. The most recent and current are: WMAP
(http://map.gsfc.nasa.gov), Planck (http://www.cosmos.esa.int/web/planck), Atacama Cosmology
Telescope (ACT) (http://act.princeton.edu), South Pole Telescope (SPT) (https://pole.uchicago.
edu), Atacama Cosmology Telescope polarisation-sensitive receiver (ACTPol) [529], SPTPol, Spider
(http://spider.princeton.edu), Polarbear (http://bolo.berkeley.edu/polarbear), Background Imaging
of Cosmic Extragalactic Polarization (BICEP2) (http://bicepkeck.org,https://www.cfa.harvard.edu/
CMB/bicep2), Keck Array (http://bicepkeck.org,https://www.cfa.harvard.edu/CMB/keckarray).
The second source of data are surveys cataloguing the angular positions and redshifts of
individual galaxies, leading to the power spectrum of fluctuations P(k, z), or the two-point correlation
function (r ). The recent and current experiments are: Baryon Oscillation Spectroscopic Survey
(BOSS) using Sloan Digital Sky Survey (SDSS) data (http://www.sdss3.org/surveys/boss.php), Dark
Energy Survey (DES) (https://www.darkenergysurvey.org), Weave (http://www.ing.iac.es/weave/
science.html), Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) (http://www.hetdex.org),
Extended Baryon Oscillation Spectroscopic Survey (eBOSS) (https://www.sdss3.org/future/eboss.
php), Mid-Scale Dark Energy Spectroscopic Instrument (MS-DESI) (https://www.skatelescope.org),
Canadian Hydrogen Intensity Mapping Experiment (CHIME) (http://chime.phas.ubc.ca), Baobab,
MeerKAT (http://www.ska.ac.za/science-engineering/meerkat), and ASKAP (http://www.atnf.csiro.
au/projects/askap/index.html).
The third source of data are weak lensing, which use the fact that images of distant galaxies are
distorted and correlated by intervening gravitational potential wells to produce statistics such as the
convergence power spectrum C` [530,531]. Some current experiments are: Dark Energy Survey (DES)
(https://www.darkenergysurvey.org), Red Cluster Sequence Lensing Survey (RCSLens) (http://www.
rcslens.org), Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) (http://www.cfhtlens.org),
59 of 82
New Instrument of Kids Arrays (NIKA2) (http://ipag.osug.fr/nika2), and Hyper Suprime-Cam (HSC)
(http://hsc.mtk.nao.ac.jp/ssp).
The final source of data are catalogues of peculiar velocities. By measuring redshifts and
radial distances of galaxies and clusters it is possible to reconstruct a radially projected map of
large-scale motions. Progress in this field will come from all three data sets above.
The main science goal of the next generation of cosmological probes is to test the Concordance
Model of cosmology. Some major experiments will be operational in the next decade.
Planck, decommissioned in 2013, marked a major milestone in CMB experiments. Its proposed
successor is the ground-based programme CMB-S4 [532], which should reach sensitivities below
103 K whose main aim is to achieve higher resolutions, probe larger scales, and measure new
observables such as polarisation.
There are various future experiments to map the mass-energy content of the universe
(including baryons and Dark Matter), either in the planning phase, or close to completion. The Large
Synoptic Space Telescope (LSST) (https://www.lsst.org), which should achieve first light in 2019,
is a ground-based telescope which will map the entire sky.
The Euclid space telescope (http://www.euclid-ec.org), due for launch in 2020, will map galaxies
and large-scale structure over the whole sky at visible and near-infrared wavelengths, providing
a catalogue of 12 billion sources at 50 million . Euclid will probe the recent universe, when galaxies
have formed and dark energy starts to dominate. Its main scientific objective is to understand the origin
of the accelerated expansion of the universe by probing the nature of dark energy using weak-lensing
observables (which include cosmic shear, higher-order distortions, and cosmic magnification), and
galaxy-clustering observations [225].
The Square Kilometre Array (SKA) (https://www.skatelescope.org), which will begin operations
in 2020, is a multi-wavelength radio telescope, built across multiple sites to achieve the largest collecting
area ever [70,533]. It should provide the highest resolution images of the radio sky, thus providing
maps of large-scale structure, and observing pulsars which should provide direct tests of General
Relativity. The SKA will observe the epoch between the emission of the CMB and the formation of
the first galaxies. Neutral hydrogen surveys (or 21 cm intensity mapping) [534], offer yet another
promising probe, as do galaxy redshift surveys [535].
Joint observables will be key in the next generation of experiments. These include the Sachs-Wolfe
effect, and the Sunyaev-Zeldovich effect [536]. We will also need to extract more observables from the
CMB. We need information on in order to constrain inflationary models.
As experiments probe larger scales at better resolutions (low ` and high `), the data analysis
and the cosmological tests (dark energy, Dark Matter, the properties of cosmic neutrinos) will require
accurate calculations on the growth of large-scale structure, which can only be achieved using N-body
simulations [537544].
At the other end of the length scale, there is particle physics. The Very Large Hadron Collider,
the successor to the Large Hadron Collider (https://home.cern/topics/large-hadron-collider)
(which can achieve energies around 14 TeV), is still in its conceptual phase. It could, if built,
probe energies around 100 TeV, allowing it to test physics beyond the Standard Model, possibly
including Supersymmetry and Grand Unified Theories, and thus provide clues on the nature of Dark
Matter and dark energy (see, e.g., [545]).
Tremendous progress is also being made in Milky Way astrophysics. With the launch on 2013 of
the Gaia mission (http://www.cosmos.esa.int/web/gaia), we will soon have improved data on Solar
System ephemerides, and on the orbits and tidal streams in the Milky Way. This will allow precise
tests of GR at Galactic scales. In particular, tidal streams provide an opportunity to probe GR at late
cosmological times and to close the gap between astronomical and cosmological scales.
Now that gravitational waves have been detected, gravitational wave physics is set to become
an established branch of astrophysics. The Laser Interferometer Gravitational-Wave Observatory
(LIGO) (https://www.ligo.caltech.edu), Advanced LIGO (aLigo) (https://www.advancedligo.mit.edu)
60 of 82
and Virgo (http://www.virgo-gw.eu) will subject GR to a battery of test in the strong regime.
Also targeting the regime of strong curvatures and potentials, the planned Event Horizon Telescope
(http://www.eventhorizontelescope.org) is a network of millimetre and sub-millimetres telescopes
being used for very-long baseline interferometry to directly image supermassive black holes in
galactic centres.
17.3. Theoretical and Computational Progress
On the theoretical front, there are four main lines of development in theories of gravity. They are
all motivated by current open questions in physics. The ongoing attempt to find a Grand Unified
Theory continues to motivate the development of string theory [546] as a framework for gravity.
Aside from the theoretical difficulties of a mathematically complex theory, the challenge for string
theory is to produce physical predictions which can be experimentally tested. The second approach
is brane theories or supergravity, in which spacetime has more than four dimensions. The resulting
field theory combines Supersymmetry (from particle physics) and General Relativity [547]. The third
approach is quantum gravity [548], in which spacetime, as a dynamical field, is a quantum object.
This implies a violation of Lorentz Invariance near the Planck scale, which in turn means that some
particle decays forbidden by Special Relativity are allowed, and possibly charge-parity-time violations
too. This motivates the search for signatures of quantum gravity in particle physics experiments.
The final approach, which is closer to the theoretical framework of the Concordance Model of
cosmology, is a phenomenological use of General Relativity. This includes the various alternative
theories that seek to explain cosmic acceleration and the missing mass, and also the various theories
which provide dark energy and Dark Matter candidates.
The early alternatives to General Relativity were motivated by theoretical considerations.
The current alternatives are mainly motivated by the open questions in cosmology. Cosmology is now
in the age of Big Data. In the last decade, the data finally caught up with the theory, and we are now in
a position to test many of the current alternative theories using statistical techniques.
The principles that underpin statistical techniques are simple enough. Under minimal
assumptions about signal and noise, it is simply a question of maximising the Gaussian likelihoods.
However, in practice they are extremely complex. The foreground contamination, the signal, the
possibly non-Gaussian noise and the systematics all have to be modelled. This requires data simulation,
so the number of maps that is generated is far larger than the number that is actually observed by the
instrument (see, e.g., [549554]).
A low signal-to-noise ratio, and higher-resolution observations of fainter signals over a larger
frequency range, result in massively big data sets. For statistical probes such as the CMB and large-scale
structure, these data sets have to be analysed as a whole, in order to correlate data. In addition to
the volume of data collected, this generates a huge amount of computational data, and requires the
appropriate computing power to carry out multiple complex calculations. We are witnessing a Moores
Law in cosmology, where the data volume of experiments increases by a factor of around 1000 every
ten years, as shown in Figure 8.
In addition to the data volume challenge, we have an algorithmic challenge. The number
of computations for cosmological data analysis depends on four main factors: the number of
observations Nt , the number of pixels Np , the number of multipoles (a function of the resolution
and frequency range) N` , and the number of iterations Ni . The first three quantities determine the
data volume. The data simulation scales at least as O( Nt ). The map-making scales as O( Ni Nt log Nt ).
The maximum-likelihood power spectrum estimation scales as O( Ni Nl Np3 ).
61 of 82
10 18
LSST 2019
Relative data volume
10 15
SDSS 2000
10 12
COBE
1989
10 9
10 6
SKA 2020
CMB-S4 2022 CMBPol

after 2020
Euclid 2020
PolarBear 2010
Planck 2008
WMAP 2001
BOOMERanG 2000
U2 and balloon
CMB experiments
1978
10 3
Bells Labs 1965
1
1960 1970 1980
1990
2000
2010
2020
2030
2040
Figure 8. The growing data volume of experiments. The data volume of each experiment is shown as an
order-of-magnitude multiple of the data volume of the 1965 Bells Labs experiment which detected the
CMB. The labels show the year of first light for each experiment. CMB surveys are marked by blue dots,
while red dots show large-scale structure surveys. A first light date for CMBPol (light blue triangle)
has not yet been fixed. Note that the vertical axis is logarithmic: the date volume increases about a
thousandfold every ten years (grey dashed line). Note too that we plot CMB, large-scale structure,
space and ground-based probes on the same graph. Ground-based probes will always tend to have a
larger data volume than space-borne probes, due to the bandwidth limit on data transmission from
spacecraft. Longer-running experiments will also have a larger data volume.
Planck marked a milestone in data science. It was the first CMB experiment in which the whole
data treatment process was parallelised, and where Monte Carlo methods were used in order to cut
down on the number of data realisation iterations that were carried out.
There are two main considerations in data analysis: the amount of data, and the complexity of
the theory. Future experiments will require sophisticated techniques, and considerable computing
power to process the vast amounts of data. Gaia is already collecting data [555], while Euclid, the SKA,
Enhanced LIGO, and CMB-S4 will soon be operational. Farther ahead, the Evolved Laser Interferometer
Space Antenna (eLISA) (https://www.elisascience.org), due for launch in 2034, will also generate huge
volumes of data, and will require the necessary computing power to test GR directly, by comparing
the data to simulated black hole collapse, inspiralling binaries.
What direction will experimental cosmology take over the coming decade? Extracting fainter
signals, such as CMB polarisation, or going to very high resolutions requires larger data volumes to
provide a higher signal-to-noise ratio, and it requires more complex models to control fainter systematic
effects. Even the Solar System tests of General Relativity will depend upon vastly greater data volumes
and computational complexity to get the full relativistic ephemerides, given the ever-increasing number
of objects being tracked in the Solar System and the significantly greater precision of the data for
each object.
62 of 82
This enables us to make a sensible prediction on future developments. The science we are able
to extract from present and future data sets will be determined by the limits on our computational
capability, and our ability and willingness to exploit it.
17.4. Conclusions
GR may well survive for another 100 years. After all, Newtonian gravity was around for 200 years.
GR has just reached its peak, when data and computing power have caught up with the theory. We are
at a pivotal moment in the history of GR. We are on the point of confirming beyond reasonable doubt
all its predictions throughout its entire domain of validity.
We have seen how modern cosmology is faced with big questions which touch the very
foundations of physics. What is this form of matter which interacts only with gravity and apparently
with nothing else? Why is the expansion of the universe accelerating? What caused the universe
to undergo a period of rapid expansion soon after the Big Bang? These questions, motivated by
cosmological observations, lead to questions about fundamental physics. Are there forces and
interactions besides the four we know of, that is, gravity, electromagnetism, and the strong and weak
nuclear forces? Are there particles beyond the Standard Model? What determines the value of the
fundamental constants of nature? What is the real structure of spacetime? Are there extra dimensions?
Science needs data, so each of these questions must be addressed through careful experiment.
The challenge of modern experimental physics is to probe nature at extreme distances and energies,
well outside the capabilities of the instruments that were available to Einstein. It has certainly come a
long way, as shown by the detection of gravitational waves in 2015, a feat which was thought to be
impossible by many of Einsteins contemporaries.
General Relativity is not the final theory of gravity, for there is no such thing. As General Relativity
turns 100, we would do well to celebrate it with a healthy does of scientific scepticism. Long live
General Relativity, and a big welcome to its eventual replacement, whether in our lifetime or not.
Acknowledgments: George F. Smoot acknowledges support through his Chaire dExcellence Universit Sorbonne
Paris Cit and the financial support of the UnivEarthS Labex programme at Universit Sorbonne Paris Cit
(ANR-10-LABX-0023 and ANR-11-IDEX-0005-02). I.D. acknowledges that the research work disclosed in this
publication is partially funded by the REACH HIGH Scholars ProgrammePost-Doctoral Grants. The grant is
part-financed by the European Union, Operational Programme IICohesion Policy 20142020.
Author Contributions: All the authors conceived the idea and contributed equally.
Conflicts of Interest: The authors declare no conflict of interest.
Abbreviations
The following abbreviations are used in this work:
AU
CDM
CMB
EEP
ECKS
eV
FLRW
GR
Gy
CDM
LLI
LLR
LPI
Mpc
PPN
SR
TeV
Astronomical Unit
Cold Dark Matter
Cosmic Microwave Background
Einstein Equivalence Principle
Einstein-Cartan-Kibble-Sciama
electronvolt
Friedmann-Lematre-Robertson-Walker
General Relativity
Gigayear (109 years)
Cold Dark Matter
Local Lorentz invariance
Lunar laser ranging
Local position invariance
Megaparsec
Parameterised post-Newtonian
Special Relativity
teraelectronvolt (1012 electronvolts)
63 of 82
References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
Feynman, R.; Leighton, R.B.; Sands, M.L. The Feynman Lectures on Physics; Addison-Wesley: Redwood City,
CA, USA, 1989; Volume 2.
Gauthier-Lafaye, F.; Holliger, P.; Blanc, P.L. Natural fission reactors in the Franceville basin, Gabon:
A review of the conditions and results of a critical event in a geologic system. Geochim. Cosmochim. Acta
1996, 60, 48314852.
Einstein, A. Erklrung der Perihelbewegung des Merkur aus der allgemeinen Relativittstheorie.
Sitzungsber. Preuss. Akad. Wiss. 1915, 1915, 831839. (In German)
Dyson, F.W.; Eddington, A.S.; Davidson, C. A Determination of the deflection of light by the suns
gravitational field, from observations made at the total eclipse of 29 May 1919. Philos. Trans. R. Soc. Lond.
Ser. A 1920, 220, 291333.
Will, C.M. The 1919 measurement of the deflection of light. Class. Quantum Gravity 2015, 32, 124001.
Coles, P. Einstein, Eddington and the 1919 Eclipse. In Historical Development of Modern Cosmology;
Martnez, V.J., Trimble, V., Pons-Bordera, M.J., Eds.; Astronomical Society of the Pacific: San Francisco,
CA, USA, 2001; Volume 252, p. 21.
Popper, K.R. Logik der Forschung; Mohr Siebeck: Tbingen, Germany, 1934.
Lahav, O.; Massimi, M. Dark energy, paradigm shifts, and the role of evidence. Astron. Geophys. 2014, 55,
3.123.15.
Ellis, G.F.R. Cosmology and verifiability. Q. J. R. Astron. Soc. 1975, 16, 245264.
Crombie, A.C. Quantification in medieval physics. Isis 1961, 52, 143160.
Bennett, J. Early modern mathematical instruments. Isis 2011, 102, 697705.
Franklin, A. Principle of inertia in the middle ages. Am. J. Phys. 1976, 44, 529543.
Sorabji, R. Matter, Space and Motion: Theories in Antiquity and Their Sequel; Duckworth: London, UK, 1988.
Cushing, J.T. Philosophical Concepts in Physics; Cambridge University Press: Cambridge, UK, 1998.
Grant, E. Scientific Thought in Fourteenth-Century Paris: Jean Buridan and Nicole Oresme. Ann. N. Y.
Acad. Sci. 1978, 314, 105126.
Babb, J. Mathematical concepts and proofs from nicole oresme: Using the History of calculus to
teach mathematics. Sci. Educ. 2005, 14, 443456.
Transtrum, M.K.; Machta, B.B.; Brown, K.S.; Daniels, B.C.; Myers, C.R.; Sethna, J.P. Perspective: Sloppiness
and emergent theories in physics, biology, and beyond. J. Chem. Phys. 2015, 143, 010901.
Galilei, G. Sidereus Nuncius Magna, Longeque Admirabilia Spectacula Pandens; Tommaso Baglioni: Venice,
Italy, 1610.
Perryman, M. The history of astrometry. Eur. Phys. J. H 2012, 37, 745792.
Newton, I. Philosophiae Naturalis Principia Mathematica, 1st ed.; Joseph Streater: London, UK, 1687.
De Maupertuis, P.L.M. Les Loix du mouvement et du repos dduites dun principe metaphysique.
In Histoire de lAcadmie Royale des Sciences et des Belles Lettres; Acadmie Royale des Sciences et des
Belles Lettres de Berlin: Berlin, Germany, 1746.
Euler, L. Rflxions sur quelques loix gnrales de la nature qui sobservent dans les effets des
forces quelconques. Acad. R. Sci. Berl. 1750, 4, 189218.
Lagrange, L. Mchanique Analytique; Chez La Veuve Desaint: Paris, France, 1788.
Hamilton, W.R. On a general method in dynamics. Philos. Trans. R. Soc. Lond. 1834, 124, 247308.
Hamilton, W.R. Second essay on a general method in dynamics. Philos. Trans. R. Soc. Lond. 1835, 125,
95144.
Galle, J.G. Account of the discovery of Le Verriers planet Neptune, at Berlin, Sept. 23, 1846. Mon. Not. R.
Astron. Soc. 1846, 7, 153.
Danjon, A. Le centenaire de la dcouverte de Neptune. Ciel Terre 1946, 62, 369.
Lexell, A.J. Recherches sur la nouvelle planete, decouverte par M. Herschel & nominee Georgium Sidus.
Acta Acad. Sci. Imp. Petropolitanae 1783, 1, 303329.
Sheehan, W. News from Front (of the Solar System): The problem with Mercury, the Vulcan hypothesis,
and General Relativitys first astronomical triumph. In Proceedings of the 227th AAS Meeting, Kissimmee,
FL, USA, 48 January 2016; Volume 227.
Poincar, H. Sur le problme des trois corps et les quations de la dynamique. Acta Math. 1890, 13, 1270.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
58.
59.
60.
64 of 82
Poincar, H. Thorie Mathmatique de la Lumire; Carr & C. Naud: Paris, France, 1889.
Maxwell, J.C. A dynamical theory of the electromagnetic field. Philos. Trans. R. Soc. Lond. 1865, 155, 489512.
Einstein, A. Zur Elektrodynamik bewegter Krper. Ann. Phys. (Germany) 1905, 17, 891921.
Einstein, A. ber den Einfluss der Schwerkraft auf die Ausbreitung des Lichtes. Ann. Phys. (Germany)
1911, 35, 898908.
Einstein, A. Lichtgeschwindigkeit und Statik des Gravitationsfeldes. Ann. Phys. (Germany) 1912, 38,
355369.
Einstein, A. Theorie des statischen Gravitationsfeldes. Ann. Phys. (Germany) 1912, 38, 443458.
Einstein, A.; Grossmann, M. Entwurf einer verallgemeinerten Relativittstheorie und eine Theorie der
Gravitation. I. Physikalischer Teil von A. Einstein II. Mathematischer Teil von M. Grossmann. Z. Math. Phys.
1913, 62, 225244 and 245261.
Le Verrier, U. Lettre de M. Le Verrier M. Faye sur la thorie de Mercure et sur le mouvement du prihlie
de cette plante. C. R. Hebd. Sances LAcad. Sci. 1859, 49, 379383.
Will, C.M. The confrontation between general relativity and experiment. Living Rev. Relativ. 2014, 17, 4.
Fierz, M.; Pauli, W. Relativistic wave equations for particles of arbitrary spin in an electromagnetic field.
Proc. R. Soc. Lond. 1939, 173, 211232.
Schild, A. Time. Tex. Q. 1960, 3, 4262.
Schild, A. Gravitational Theories of the Whitehead type and the principle of equivalence. In Evidence for
Gravitational Theories; Mller, C., Ed.; Academic Press: New York, NY, USA, 1962.
Schild, A. Lectures on General Relativity Theory. In Relativity Theory and Astrophysics: I, Relativity and
Cosmology: II, Galactic Structure: III, Stellar Structure; Ehlers, J., Ed.; American Mathematical Society:
Providence, RI, USA, 1967.
Riemann, B. Bernhard Riemanns Gesammelte mathematische Werke und wissenschaftlicher Nachlass, 1st ed;
Weber, H., Ed.; Teubner: Leipzig, Germany, 1876.
Riemann, B. Ueber die Hypothesen, welche der Geometrie zu Grunde liegen. Abh. Kn. Ges. Wiss. Gtt.
1868, 13, 133150.
Einstein, A. Die Feldgleichungen der Gravitation. Sitzungsber. Preuss. Akad. Wiss. 1915, 1915, 844847.
(In German)
Einstein, A. Zur Allgemeinen Relativittstheorie. Sitzungsber. Preuss. Akad. Wiss. 1915, 1915, 799801.
(In German)
Ni, W.T. Theoretical Frameworks for Testing Relativistic Gravity. IV. A Compendium of Metric Theories of
Gravity and Their Post-Newtonian Limits. Astrophys. J. 1972, 176, 769796.
Thorne, K.S.; Ni, W.T.; Will, C.M. Theoretical frameworks for testing relativistic gravity: A review.
In Proceedings of the Conference on Experimental Tests of Gravitational Theories, Pasadena, CA, USA,
1113 November 1970; pp. 1031.
Lovelock, D. The uniqueness of the Einstein field equations in a four-dimensional space. Arch. Ration.
Mech. Anal. 1969, 33, 5470.
Lovelock, D. The Einstein Tensor and Its Generalizations. J. Math. Phys. 1971, 12, 498501.
Lovelock, D. The Four-Dimensionality of Space and the Einstein Tensor. J. Math. Phys. 1972, 13, 874876.
Cartan, E. Sur les quations de la gravitation dEinstein. J. Math. Pures Appl. 1922, 1, 141204.
Anderson, J.L. Principles of Relativity Physics; Academic Press: New York, NY, USA, 1967.
Norton, J.D. General covariance, gauge theories and the Kretschmann objection. In Symmetries in Physics:
Philosophical Reflections; Brading, K.A., Castellani, E., Eds.; Cambridge University Press: Cambridge,
UK, 2003.
Kretschmann, E. ber den physikalischen Sinn der Relativittspostulate, A. Einsteins neue und seine
ursprungliche Relativittstheorie. Ann. Phys. (Germany) 1917, 53, 575614.
Nordstrm, G. Zur Theorie der Gravitation vom Standpunkt des Relativitsprinzip. Ann. Phys. (Germany)
1913, 42, 533554.
Einstein, A.; Fokker, A.D. The North current gravitation theory from the viewpoint of absolute
differential calculus. Ann. Phys. (Germany) 1914, 44, 321328.
Whitehead, A.N. The Principle of Relativity; Cambridge University Press: Cambridge, UK, 1922.
Will, C.M. Relativistic gravity in the solar system, II: Anisotropy in the Newtonian gravitational constant.
Astrophys. J. 1971, 169, 141156.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
80.
81.
82.
83.
84.
85.
86.
87.
88.
89.
90.
91.
65 of 82
Birkhoff, G.D. Matter, electricity and gravitation in flat spacetime. Proc. Natl. Acad. Sci. USA 1943, 29,
231239.
Jordan, P. Zum gegenwrtigen Stand der Diracschen kosmologischen Hypothesen. Z. Phys. 1959, 157,
112121.
Brans, C.; Dicke, R.H. Machs principle and a relativistic theory of gravitation. Phys. Rev. 1961, 124, 925935.
Dicke, R.H. Machs principle and invariance under transformation of units. Phys. Rev. 1962, 125, 21632167.
Will, C.M. Theory and Experiment in Gravitational Physics, 2nd ed.; Cambridge University Press: Cambridge,
UK, 1993.
Psaltis, D. Constraining Brans-Dicke Gravity with Accreting Millisecond Pulsars in Ultracompact Binaries.
Astrophys. J. 2008, 688, 12821287.
Bertotti, B.; Iess, L.; Tortora, P. A test of general relativity using radio links with the Cassini spacecraft.
Nature 2003, 425, 374376.
Avilez, A.; Skordis, C. Cosmological Constraints on Brans-Dicke Theory. Phys. Rev. Lett. 2014, 113, 011101.
Ooba, J.; Ichiki, K.; Chiba, T.; Sugiyama, N. Planck constraints on scalar-tensor cosmology and the variation
of the gravitational constant. Phys. Rev. D 2016, 93, 122002.
Bull, P. Extending Cosmological Tests of General Relativity with the Square Kilometre Array. Astrophys. J.
2016, 817, 26.
Damour, T.; Nordtvedt, K. General relativity as a cosmological attractor of tensor-scalar theories.
Phys. Rev. Lett. 1993, 70, 22172219.
Damour, T.; Nordtvedt, K. Tensor-scalar cosmological models and their relaxation toward general relativity.
Phys. Rev. D 1993, 48, 34363450.
Kustaanheimo, P. Route dependence of the gravitational redshift. Phys. Lett. 1966, 23, 7577.
Kustaanheimo, P.E.; Nuotio, V.S. Relativistic Theories of Gravitation I: The One-Body Problem; University of
Helsinki: Helsinki, Finland, 1967; unpublished.
Einstein, A. Nherungsweise Integration der Feldgleichungen der Gravitation. Sitzungsber. Preuss.
Akad. Wiss. 1916, 1916 Pt 1, 688696. (In German)
Kumar, D.; Soni, V. Single particle Schrdinger equation with gravitational self-interaction. Phys. Lett. A
2000, 271, 157166.
Deser, S. Gravity from self-interaction redux. Gen. Relativ. Gravit. 2010, 42, 641646.
Hertz, H. Untersuchungen ber die Ausbreitung der Elektrischen Kraft; J.A. Barth: Leipzig, Germany, 1894.
Einstein, A. ber einen die Erzeugung und Verwandlung des Lichtes betreffenden heuristischen
Gesichtspunkt. Ann. Phys. (Germany) 1905, 17, 132148.
LIGO Collaboration. Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett.
2016, 116, 061102.
Yang, C.L.; Mills, R.L. Conservation of isotopic spin and isotopic gauge invariance. Phys. Rev. 1954, 96,
191195.
Yang, C.L.; Mills, R.L. Isotopic Spin Conservation and a Generalized Gauge Invariance. Phys. Rev. 1954,
95, 631.
Weyl, H. Gravitation und Elektrizitt. In Das Relativittsprinzip; J. A. Barth: Leipzig, Germany, 1918.
Weyl, H. Raum-Zeit-Materie; Springer: Berlin, Germany, 1918.
Christoffel, E.B. ber die Transformation der homogenen Differentialausdrcke zweiten Grades. J. Reine
Angew. Math. 1869, 70, 4670.
Goenner, H.F.M. On the History of Unified Field Theories. Living Rev. Relativ. 2004, 7, 2.
Smiciklas, M.; Brown, J.M.; Cheuk, L.W.; Smullin, S.J.; Romalis, M.V. New test of local lorentz invariance
using a 21 Ne-Rb-K comagnetometer. Phys. Rev. Lett. 2011, 107, 171604.
Peck, S.K.; Kim, D.K.; Stein, D.; Orbaker, D.; Foss, A.; Hummon, M.T.; Hunter, L.R. Limits on local Lorentz
invariance in mercury and cesium. Phys. Rev. A 2012, 86, 012109.
Hohensee, M.A.; Leefer, N.; Budker, D.; Harabati, C.; Dzuba, V.A.; Flambaum, V.V. Limits on violations of
lorentz symmetry and the einstein equivalence principle using radio-frequency spectroscopy of atomic
dysprosium. Phys. Rev. Lett. 2013, 111, 050401.
Utiyama, R. Invariant theoretical interpretation of interaction. Phys. Rev. 1956, 101, 15971607.
Kibble, T.W.B. Lorentz invariance and the gravitational field. J. Math. Phys. 1961, 2, 212221.
92.
93.
94.
95.
96.
97.
98.
99.
100.
101.
102.
103.
104.
105.
106.
107.
108.
109.
110.
111.
112.
113.
114.
115.
116.
117.
118.
119.
120.
66 of 82
Carmeli, M.; Malin, S. Reformulation of general relativity as a gauge theory. Ann. Phys. (Germany) 1977,
103, 208232.
Ivanenko, D.; Sardanashvily, G. The gauge treatment of gravity. Phys. Rep. 1983, 94, 145.
Sardanashvily, G.; Zakharov, O. Gauge Gravitation Theory; World Scientific: Singapore, 1992.
Jackson, J.D.; Okun, L.B. Historical roots of gauge invariance. Rev. Mod. Phys. 2001, 73, 663680.
Sardanashvily, G. Classical gauge theory of gravity. Theor. Math. Phys. 2002, 132, 11631171.
Cartan, . Sur une gnralisation de la notion de courboure de Riemann et les espaces torsion. C. R. Acad.
Sci. (Paris) 1922, 174, 593595.
Cartan, . Sur les varits connexion affine et la thorie de la relativit gnralise (premire partie).
Ann. Sci. cole Norm. Super. 1923, 40, 325412.
Cartan, . Sur les varits connexion affine et la thorie de la relativit gnralise (suite). Ann. Sci. cole
Norm. Super. 1924, 41, 125.
Cartan, . Sur les varits connexion affine et la thorie de la relativit gnralise (deuxime partie).
Ann. Sci. cole Norm. Super. 1925, 42, 1788.
Weyl, H. A remark on the coupling of gravitation and electron. Phys. Rev. 1950, 77, 699701.
Sciama, D. The physical structure of general relativity. Rev. Mod. Phys. 1964, 36, 1103.
Weyssenhoff, J.; Raabe, A. Relativistic Dynamics of spin-fluids and spin-particles. Acta Phys. Pol. 1947,
9, 718.
Costa de Beauregard, O. Translational inertial spin effect. Phys. Rev. 1963, 129, 466471.
Cosserat, E.; Cosserat, F. Thorie des Corps dformables; Hermann: Paris, France, 1909.
Puetzfeld, D.; Obukhov, Y.N. Prospects of detecting spacetime torsion. Int. J. Mod. Phys. D 2014, 23,
1442004.
Hehl, F.W.; von der Heyde, P.; Kerlick, G.D.; Nester, J.M. General relativity with spin and torsion:
Foundations and prospects. Rev. Mod. Phys. 1976, 48, 393416.
Lasenby, A.N.; Doran, C.J.; Gull, S.F. Cosmological consequences of a flat-space theory of gravity. In Clifford
Algebras and Their Applications in Mathematical Physics; Brackx, F., Delanghe, R., Serras, H., Eds.; Kluwer:
Dordrecht, The Netherlands, 1993; p. 387.
Lasenby, A.N.; Doran, C.J.; Gull, S.F. Grassmann calculus, pseudoclassical mechanics and geometric algebra.
J. Math. Phys. 1993, 34, 36833712.
Lasenby, A.N.; Doran, C.J.L.; Gull, S.F. A multivector derivative approach to Lagrangian field theory.
Found. Phys. 1993, 23, 12951327.
Lasenby, A.N.; Doran, C.J.L.; Gull, S.F. 2-spinors, twistors and supersymmetry in the spacetime algebra.
In Spinors, Twistors, Clifford Algebras and Quantum Deformations; Oziewicz, Z., Jancewicz, B.,
Borowiec, A., Eds.; Kluwer: Dordrecht, The Netherlands, 1993; p. 233.
Lasenby, A.N.; Doran, C.J.L.; Gull, S.F. Astrophysical and cosmological consequences of a gauge theory
of gravity. In Current Topics in Astrofundamental Physics; Snchez, N., Zichichi, A., Eds.; World Scientific:
Singapore, 1995; p. 359.
Lasenby, A.N.; Doran, C.J.; Gull, S.F. Gravity, gauge theories and geometric algebra. Philos. Trans. R.
Soc. Lond. 1998, A356, 487582.
Carmeli, M.; Leibowitz, E.; Nissani, N. Gravitation: SL(2, C) Gauge Theory and Conservation Laws;
World Scientific: Singapore, 1990.
Mukunda, N. An elementary introduction to the gauge theory approach to gravity. In Gravitation, Gauge
Theories and the Early Universe; Iyer, B.R., Mukunda, N., Vishvershwara, C.V., Eds.; Kluwer: Dordrecht,
The Netherlands, 1989.
Misner, C.W.; Thorne, K.S.; Wheeler, J.A. Gravitation; Freeman: San Francisco, CA, USA, 1973.
Berti, E.; Barausse, E.; Cardoso, V.; Gualtieri, L.; Pani, P.; Sperhake, U.; Stein, L.C.; Wex, N.; Yagi, K.; Baker, T.;
et al. Testing general relativity with present and future astrophysical observations. Class. Quantum Gravity
2015, 32, 243001.
Poincar, H. Ltat actuel et lavenir de la physique mathmatique. Bull. Sci. Math. 1904, 2, 302324.
Kaluza, T. Zum Unittsproblem in der Physik; Sitzungsberichte der Kniglich Preuischen Akademie der
Wissenschaften (Berlin): Berlin, Germany, 1921; pp. 966972.
Klein, O. Quantentheorie und fnfdimensionale Relativittstheorie. Z. Phys. 1926, 37, 895906.
121.
122.
123.
124.
125.
126.
127.
128.
129.
130.
131.
132.
133.
134.
135.
136.
137.
138.
139.
140.
141.
142.
143.
144.
145.
146.
147.
148.
149.
150.
151.
152.
153.
67 of 82
Milne, E.A. Kinematic Relativity; A Sequel to Relativity, Gravitation and World Structure; Clarendon Press:
Oxford, UK, 1948.
Thiry, Y. Sur la rgularit des champs gravitationels et lectromagntiques dans les thories unitaires. C. R.
Hebd. Sances Acad. Sci. 1948, 226, 18811882.
Belinfante, F.; Swihart, J. Phenomenological linear theory of gravitation: Part I. Classical mechanics.
Ann. Phys. (Germany) 1957, 1, 168195
Will, C.M.; Nordtvedt, K., Jr. Conservation laws and preferred frames in relativistic gravity. I.
Preferred-frame theories and an extended PPN formalism. Astrophys. J. 1972, 177, 757774.
Barker, B.M. General scalar-tensor theory of gravity with constant G. Astrophys. J. 1978, 219, 511.
Rosen, N. A bi-metric theory of gravitation. Gen. Relativ. Gravitat. 1973, 4, 435447.
Rosen, N. A bi-metric theory of gravitation. II. Gen. Relativ. Gravitat. 1975, 6, 259268.
Rastall, P. A theory of gravity. Can. J. Phys. 1976, 54, 6675.
Buchdahl, H.A. Non-linear lagrangians and cosmological theory. Mon. Not. R. Astron. Soc. 1970, 150, 18.
Starobinsky, A. A new type of isotropic cosmological models without singularity. Phys. Lett. B 1980, 91,
99102.
Milgrom, M. A modification of the Newtonian dynamics as a possible alternative to the hidden
mass hypothesis. Astrophys. J. 1983, 270, 365370.
Milgrom, M. A modification of the newtonian dynamicsImplications for galaxies. Astrophys. J. 1983,
270, 371389.
Milgrom, M. A modification of the newtonian dynamicsImplications for galaxy systems. Astrophys. J.
1983, 270, 384415.
Dvali, G.; Gabadadze, G.; Porrati, M. 4D gravity on a brane in 5D Minkowski space. Phys. Lett. B 2000,
485, 208214.
Bekenstein, J.D. Relativistic gravitation theory for the modified Newtonian dynamics paradigm.
Phys. Rev. D 2004, 70, 083509.
Seifert, M.D. Stability of spherically symmetric solutions in modified theories of gravity. Phys. Rev. D 2007,
76, 064002.
Reyes, R.; Mandelbaum, R.; Seljak, U.; Baldauf, T.; Gunn, J.E.; Lombriser, L.; Smith, R.E. Confirmation of
general relativity on large scales from weak lensing and galaxy velocities. Nature 2010, 464, 256258.
Einstein, A. Kosmologische Betrachtungen zur allgemeinen Relativittstheorie. Sitzungsber. Preuss.
Akad. Wiss. 1917, 1917, 142152.
Friedmann, A. ber die Krmmung des Raumes. Z. Phys. 1922, 10, 377386.
Friedmann, A. ber die Mglichkeit einer Welt mit konstanter negativer Krmmung des Raumes. Z. Phys.
1924, 21, 326332.
Hubble, E. A relation between distance and radial velocity among extra-galactic nebulae. Proc. Natl. Acad.
Sci. USA 1929, 15, 168173.
Lematre, G. Un univers homogne de masse constante et de rayon croissant rendant compte de la vitesse
radiale des nbuleuses extra-galactiques. Ann. Soc. Sci. Brux. 1927, A47, 4959.
Robertson, H.P. Kinematics and world structure. Astrophys. J. 1935, 82, 248301.
Robertson, H.P. Kinematics and world structure II. Astrophys. J. 1936, 83, 187201.
Robertson, H.P. Kinematics and world structure III. Astrophys. J. 1936, 83, 257271.
Walker, A.G. On Milnes theory of world-structure. Proc. Lond. Math. Soc. 1937, 42, 90127.
Ellis, G.F.R.; Stoeger, W. Horizons in inflationary universes. Class. Quantum Gravity 1988, 5, 207220.
Einstein, A.; de Sitter, W. On the relation between the expansion and the mean density of the universe.
Proc. Natl. Acad. Sci. USA 1932, 18, 213214.
Adams, F.C.; Laughlin, G. A dying universe: The long-term fate and evolution of astrophysical objects.
Rev. Mod. Phys. 1997, 69, 337372.
Caldwell, R.R.; Kamionkowski, M.; Weinberg, N.N. Phantom energy: Dark energy with w < 1 causes
a cosmic doomsday. Phys. Rev. Lett. 2003, 91, 071301.
Planck Collaboration. Planck 2015 Results. XIII. Cosmological Parameters. Astron. Astrophys. 2016,
594, A13.
Alpher, R.A.; Bethe, H.; Gamow, G. The origin of chemical elements. Phys. Rev. 1948, 73, 803804.
Alpher, R.A.; Herman, R. Evolution of the universe. Nature 1948, 162, 774775.
154.
155.
156.
157.
158.
159.
160.
161.
162.
163.
164.
165.
166.
167.
168.
169.
170.
171.
172.
173.
174.
175.
176.
177.
178.
179.
180.
181.
182.
183.
184.
68 of 82
Gamow, G. The evolution of the universe. Nature 1948, 162, 680682.

Gamow, G. The origin of elements and the separation of galaxies. Phys. Rev. 1948, 74, 505506.
Penzias, A.A.; Wilson, R.W. A measurement of excess antenna temperature at 4080 Mc/s. Astrophys. J.
1965, 142, 419421.
Christenson, J.H.; Cronin, J.W.; Fitch, V.L.; Turlay, R. Evidence for the 2 decay of the K20 meson.
Phys. Rev. Lett. 1964, 13, 138140.
Sakharov, A.D. Violation of CP invariance, c asymmetry, and baryon asymmetry of the universe. JETP Lett.
1967, 5, 2427.
Kuzmin, V.A. CP-noninvariance and baryon asymmetry of the universe. J. Exp. Theor. Phys. Lett. 1970,
12, 228230.
Khlopov, M. Cosmological reflection of particle symmetry. Symmetry 2016, 8, 81.
Fixsen, D.J. The temperature of the cosmic microwave background. Astrophys. J. 2009, 707, 916920.
Dunkley, J.; Komatsu, E.; Nolta, M.R.; Spergel, D.N.; Larson, D.; Hinshaw, G.; Page, L.; Bennett, C.L.;
Gold, B.; Jarosik, N.; et al. Five-year Wilkinson Microwave Anisotropy Probe observations: Likelihoods
and parameters from the WMAP data. Astrophys. J. Suppl. 2009, 180, 306329.
Hannestad, S. Primordial neutrinos. Annu. Rev. Nucl. Part. Sci. 2006, 56, 137161.
Lesgourgues, J.; Pastor, S. Massive neutrinos and cosmology. Phys. Rep. 2006, 429, 307379.
Lesgourgues, J.; Pastor, S. Neutrino cosmology and Planck. New J. Phys. 2014, 16, 065002.
Elgary, .; Lahav, O. Neutrino masses from cosmological probes. New J. Phys. 2005, 7, 61.
Abazajian, K.; Fuller, G.M.; Patel, M. Sterile neutrino hot, warm, and cold Dark Matter. Phys. Rev. D 2001,
64, 023501.
Boyarsky, A.; Ruchayskiy, O.; Shaposhnikov, M. The role of sterile neutrinos in cosmology and astrophysics.
Annu. Rev. Nucl. Part. Sci. 2009, 59, 191214.
Mangano, G.; Miele, G.; Pastor, S.; Peloso, M. A precision calculation of the effective number of cosmological
neutrinos. Phys. Lett. B 2002, 534, 816.
Feeney, S.M.; Peiris, H.V.; Verde, L. Is there evidence for additional neutrino species from cosmology?
J. Cosmol. Astropart. Phys. 2013, 4, 036.
Olive, K.; Group, P.D. Review of particle physics. Chin. Phys. C 2014, 38, 090001.
Primack, J.R.; Gross, M.A.K. Hot Dark Matter in cosmology. In Current Aspects of Neutrino Physics;
Caldwell, D.O., Ed.; Springer: Berlin, Germany, 2001; pp. 287308.
Gerstein, G.; Zeldovich, Y.B. Rest mass of the muonic neutrino and cosmology (English translation).
Lett. J. Exp. Theor. Phys. 1966, 4, 120122.
Marx, G.; Szalay, A.S. Cosmological limit on the neutretto rest mass. In Proceedings of the Neutrino 72,
Balatonfured, Hungary, 1117 June 1972; Volume 1, p. 123.
Cowsik, R.; McClelland, J. An upper limit on the neutrino rest mass. Phys. Rev. Lett. 1972, 29, 669670.
Szalay, A.S.; Marx, G. Neutrino rest mass from cosmology. Astron. Astrophys. 1976, 49, 437441.
Tegmark, M.; Eisenstein, D.J.; Strauss, M.A.; Weinberg, D.H.; Blanton, M.R.; Frieman, J.A.; Fukugita, M.;
Gunn, J.E.; Hamilton, A.J.S.; Knapp, G.R.; et al. Cosmological constraints from the SDSS luminous
red galaxies. Phys. Rev. D 2006, 74, 123507.
Shiraishi, M.; Ichikawa, K.; Ichiki, K.; Sugiyama, N.; Yamaguchi, M. Constraints on neutrino masses from
WMAP5 and BBN in the lepton asymmetric universe. J. Cosmol. Astropart. Phys. 2009, 7, 5.
Tereno, I.; Schimd, C.; Uzan, J.P.; Kilbinger, M.; Vincent, F.H.; Fu, L. CFHTLS weak-lensing constraints on
the neutrino masses. Astron. Astrophys. 2009, 500, 657665.
Ichiki, K.; Takada, M.; Takahashi, T. Constraints on neutrino masses from weak lensing. Phys. Rev. D 2009,
79, 023520.
Xia, J.Q.; Granett, B.R.; Viel, M.; Bird, S.; Guzzo, L.; Haehnelt, M.G.; Coupon, J.; McCracken, H.J.; Mellier, Y.
Constraints on massive neutrinos from the CFHTLS angular power spectrum. J. Cosmol. Astropart. Phys.
2012, 6, 010.
Olive, K.A. Inflation. Phys. Rep. 1990, 190, 307403.
Guth, A.H. (Ed.) The Inflationary Universe. The Quest for a New Theory of Cosmic Origins; Addison-Wesley:
Reading, MA, USA, 1997.
Peacock, J.A. Cosmological Physics; Cambridge University Press: Cambridge, UK, 1999.
185.
186.
187.
188.
189.
190.
191.
192.
193.
194.
195.
196.
197.
198.
199.
200.
201.
202.
203.
204.
205.
206.
207.
208.
209.
210.
211.
212.
213.
69 of 82
Liddle, A.R.; Lyth, D.H. Cosmological Inflation and Large-Scale Structure; Cambridge University Press:
Cambridge, UK, 2000.
Uzan, J.P. Inflation in the standard cosmological model. C. R. Phys. 2015, 16, 875890.
Misner, C.W. Mixmaster universe. Phys. Rev. Lett. 1969, 22, 10711074.
Belinskij, V.A.; Khalatnikov, I.M.; Lifshits, E.M. Oscillatory approach to a singular point in the relativistic
cosmology. Adv. Phys. 1970, 19, 525573.
Starobinsky, A.A. Spectrum of relict gravitational radiation and the early state of the universe. Sov. J. Exp.
Theor. Phys. Lett. 1979, 30, 682685.
t Hooft, G. Magnetic monopoles in unified gauge theories. Nucl. Phys. B 1974, 79, 276284.
Preskill, J. Magnetic monopoles. Annu. Rev. Nucl. Part. Sci. 1984, 34, 461530.
Zeldovich, Y.B.; Khlopov, M.Y. On the concentration of relic monopoles in the universe. Phys. Lett. B 1979,
79, 239241.
Preskill, J.P. Cosmological production of superheavy magnetic monopoles. Phys. Rev. Lett. 1979, 43,
13651368.
Guth, A.H. Inflationary universe: A possible solution to the horizon and flatness problems. Phys. Rev. D
1981, 23, 347356.
Linde, A. Inflation and Quantum Cosmology; Academic Press: New York, NY, USA, 1989.
Lidsey, J.E.; Liddle, A.R.; Kolb, E.W.; Copeland, E.J.; Barreiro, T.; Abney, M. Reconstructing the inflaton
potentialAn overview. Rev. Mod. Phys. 1997, 69, 373410.
Sato, K. Cosmological baryon-number domain structure and the first order phase transition of a vacuum.
Phys. Lett. B 1981, 99, 6670.
Einhorn, M.B.; Sato, K. Monopole production in the very early universe in a first-order phase transition.
Nucl. Phys. B 1981, 180, 385404.
Mukhanov, V.F.; Chibisov, G.V. Energy of vacuum and the large-scale structure of the universe.
Zhurnal Eksper. Teor. Fiz. 1982, 83, 475487.
Mukhanov, V.F.; Chibisov, G.V. Quantum fluctuations and a nonsingular universe. Sov. J. Exp. Theor.
Phys. Lett. 1981, 33, 532535.
Guth, A.H.; Pi, S.Y. Fluctuations in the new inflationary universe. Phys. Rev. Lett. 1982, 49, 11101113.
Hawking, S.W. The development of irregularities in a single bubble inflationary universe. Phys. Lett. B
1982, 115, 295297.
Starobinsky, A.A. Dynamics of phase transition in the new inflationary universe scenario and generation
of perturbations. Phys. Lett. B 1982, 117, 175178.
Bardeen, J.M.; Steinhardt, P.J.; Turner, M.S. Spontaneous creation of almost scale-free density perturbations
in an inflationary universe. Phys. Rev. D 1983, 28, 679693.
Planck Collaboration. Planck 2015 Results. XX. Constraints on Inflation Astron. Astrophys. 2016, 594, A20.
Liddle, A.R.; Lyth, D.H. COBE, gravitational waves, inflation and extended inflation. Phys. Lett. B 1992,
291, 391398.
Liddle, A.R.; Lyth, D.H. The cold Dark Matter density perturbation. Phys. Rep. 1993, 231, 1105.
Linde, A.D. A new inflationary universe scenario: A possible solution of the horizon, flatness, homogeneity,
isotropy and primordial monopole problems. Phys. Lett. B 1982, 108, 389393.
Dodelson, S. Modern Cosmology; Academic Press: New York, NY, USA, 2003.
Kosowsky, A.; Turner, M.S. CBR anisotropy and the running of the scalar spectral index. Phys. Rev. D 1995,
52, 17391743.
Smoot, G.F.; Bennett, C.L.; Kogut, A.; Wright, E.L.; Aymon, J.; Boggess, N.W.; Cheng, E.S.; de Amici, G.;
Gulkis, S.; Hauser, M.G.; et al.; Structure in the COBE differential microwave radiometer first-year maps.
Astrophys. J. 1992, 396, L1L5.
Fixsen, D.J.; Cheng, E.S.; Cottingham, D.A.; Eplee, R.E., Jr.; Isaacman, R.B.; Mather, J.C.; Meyer, S.S.;
Noerdlinger, P.D.; Shafer, R.A.; Weiss, R.; et al. Cosmic microwave background dipole spectrum measured
by the COBE FIRAS instrument. Astrophys. J. 1994, 420, 445449.
Dwek, E.; Arendt, R.G.; Hauser, M.G.; Fixsen, D.; Kelsall, T.; Leisawitz, D.; Pei, Y.C.; Wright, E.L.;
Mather, J.C.; Moseley, S.H.; et al. The COBE diffuse infrared background experiment search for the cosmic
infrared background. IV. Cosmological implications. Astrophys. J. 1998, 508, 106122.
214.
215.
216.
217.
218.
219.
220.
221.
222.
223.
224.
225.
226.
227.
228.
229.
230.
231.
232.
233.
234.
235.
236.
237.
238.
239.
70 of 82
Spergel, D.N.; Verde, L.; Peiris, H.V.; Komatsu, E.; Nolta, M.R.; Bennett, C.L.; Halpern, M.; Hinshaw, G.;
Jarosik, N.; Kogut, A.; et al. First-year Wilkinson Microwave Anisotropy Probe (WMAP) observations:
Determination of cosmological parameters. Astrophys. J. Suppl. 2003, 148, 175194.
Peiris, H.V.; Komatsu, E.; Verde, L.; Spergel, D.N.; Bennett, C.L.; Halpern, M.; Hinshaw, G.; Jarosik, N.;
Kogut, A.; Limon, M.; et al. First-year Wilkinson Microwave Anisotropy Probe (WMAP) observations:
Implications for inflation. Astrophys. J. Suppl. 2003, 148, 213231.
Spergel, D.N.; Bean, R.; Dor, O.; Nolta, M.R.; Bennett, C.L.; Dunkley, J.; Hinshaw, G.; Jarosik, N.;
Komatsu, E.; Page, L.; et al. Three-year Wilkinson Microwave Anisotropy Probe (WMAP) observations:
Implications for cosmology. Astrophys. J. Suppl. 2007, 170, 377408.
Komatsu, E.; Dunkley, J.; Nolta, M.R.; Bennett, C.L.; Gold, B.; Hinshaw, G.; Jarosik, N.; Larson, D.;
Limon, M.; Page, L.; et al. Five-year Wilkinson Microwave Anisotropy Probe observations: Cosmological
interpretation. Astrophys. J. Suppl. 2009, 180, 330376.
Komatsu, E.; Smith, K.M.; Dunkley, J.; Bennett, C.L.; Gold, B.; Hinshaw, G.; Jarosik, N.; Larson, D.;
Nolta, M.R.; Page, L.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) observations:
Cosmological interpretation. Astrophys. J. Suppl. 2011, 192, 18.
Hinshaw, G.; Larson, D.; Komatsu, E.; Spergel, D.N.; Bennett, C.L.; Dunkley, J.; Nolta, M.R.; Halpern, M.;
Hill, R.S.; Odegard, N.; et al. Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) observations:
Cosmological parameter results. Astrophys. J. Suppl. 2013, 208, 19.
Planck Collaboration. Planck 2013 results. XVI. Cosmological parameters. Astron. Astrophys. 2014,
571, A16.
Planck Collaboration. Planck 2013 results. XXII. Constraints on inflation. Astron. Astrophys. 2014, 571, A22.
BICEP2 Collaboration. Detection of B-mode polarization at degree angular scales by BICEP2. Phys. Rev. Lett.
2014, 112, 241101.
Planck Collaboration. Planck intermediate results XXX. The angular power spectrum of polarized dust
emission at intermediate and high Galactic latitudes. Astron. Astrophys. 2016, 586, A133.
Huang, Z.; Verde, L.; Vernizzi, F. Constraining inflation with future galaxy redshift surveys.
Amendola, L.; Appleby, S.; Bacon, D.; Baker, T.; Baldi, M.; Bartolo, N.; Blanchard, A.; Bonvin, C.; Borgani, S.;
Branchini, E.; et al. Cosmology and fundamental physics with the Euclid satellite. Living Rev. Relativ. 2013,
16, 6.
Namikawa, T.; Yamauchi, D.; Sherwin, B.; Nagata, R. Delensing cosmic microwave background B modes
with the Square Kilometre Array Radio Continuum Survey. Phys. Rev. D 2016, 93, 043527.
Kapteyn, J.C. First attempt at a theory of the arrangement and motion of the sidereal system. Astrophys. J.
1922, 55, 302328.
Oort, J.H. The force exerted by the stellar system in the direction perpendicular to the galactic plane and
some related problems. Bull. Astron. Inst. Neth. 1932, 6, 249287.
Zwicky, F. Die Rotverschiebung von extragalaktischen Nebeln. Helv. Phys. Acta 1933, 6, 110127.
Zwicky, F. On the Masses of Nebulae and of Clusters of Nebulae. Astrophys. J. 1937, 86, 217246.
Freese, K. Review of observational evidence for Dark Matter in the universe and in upcoming searches for
dark stars. In Proceedings of the Dark Energy and Dark Matter: Observations, Experiments and Theories,
Lyons, France, 711 July 2008; Volume 36, pp. 113126.
Kahn, F.D.; Woltjer, L. Intergalactic matter and the galaxy. Astrophys. J. 1959, 130, 705717.
Roberts, M.S.; Rots, A.H. Comparison of rotation curves of different galaxy types. Astron. Astrophys. 1973,
26, 483485.
Einasto, J.; Kaasik, A.; Saar, E. Dynamic evidence on massive coronas of galaxies. Nature 1974, 250, 309310.
Ostriker, J.P.; Peebles, P.J.E.; Yahil, A. The size and mass of galaxies, and the mass of the universe.
Astrophys. J. 1974, 193, L1L4.
Rubin, V.C.; Thonnard, N.; Ford, W.K., Jr. Extended rotation curves of high-luminosity spiral galaxies.
IVSystematic dynamical properties, SA through SC. Astrophys. J. 1978, 225, L107L111.
Navarro, J.F.; Frenk, C.S.; White, S.D.M. The structure of cold Dark Matter halos. Astrophys. J. 1996, 462, 563.
Fukugita, M.; Hogan, C.J.; Peebles, P.J.E. The cosmic baryon budget. Astrophys. J. 1998, 503, 518530.
Sachs, R.K.; Wolfe, A.M. Perturbations of a cosmological model and angular variations of the microwave
background. Astrophys. J. 1967, 147, 7390.
240.
241.
242.
243.
244.
245.
246.
247.
248.
249.
250.
251.
252.
253.
254.
255.
256.
257.
258.
259.
260.
261.
262.
263.
264.
265.
266.
267.
268.
71 of 82
Eisenstein, D.J. Dark energy and cosmic sound. New Astron. Rev. 2005, 49, 360365.
Bassett, B.; Hlozek, R. Baryon acoustic oscillations. In Dark Energy: Observational and Theoretical Approaches;
Ruiz-Lapuente, P., Ed.; Cambridge University Press: Cambridge, UK, 2010; p. 246.
Schneider, P. Weak gravitational lensing. In Gravitational Lensing: Strong, Weak and Micro, Saas-Fee Advanced
Courses 33; Schneider, P., Kochanek, C.S., Wambsganss, J., Eds.; Springer: Berlin/Heidelberg, Germany,
2006; p. 269.
Bacon, D.J.; Rfrgier, A.R.; Ellis, R.S. Detection of weak gravitational lensing by large-scale structure.
Mon. Not. R. Astron. Soc. 2000, 318, 625640.
Kaiser, N. A new shear estimator for weak-lensing observations. Astrophys. J. 2000, 537, 555577.
Van Waerbeke, L.; Mellier, Y.; Erben, T.; Cuillandre, J.C.; Bernardeau, F.; Maoli, R.; Bertin, E.;
Mc Cracken, H.J.; Le Fvre, O.; Fort, B.; et al. Detection of correlated galaxy ellipticities from CFHT data:
First evidence for gravitational lensing by large-scale structures. Astron. Astrophys. 2000, 358, 3044.
Wittman, D.M.; Tyson, J.A.; Kirkman, D.; DellAntonio, I.; Bernstein, G. Detection of weak gravitational
lensing distortions of distant galaxies by cosmic Dark Matter at large scales. Nature 2000, 405, 143148.
Bertone, G.; Hooper, D.; Silk, J. Particle Dark Matter: Evidence, candidates and constraints. Phys. Rep. 2005,
405, 279390.
Salati, P. The bestiary of Dark Matter species. In Proceedings of the Dark Energy and Dark Matter:
Observations, Experiments and Theories, Lyons, France, 711 July 2008; Volume 36, pp. 175186.
Garrett, K.; Duda,

G. Dark matter: A primer. Adv. Astron. 2011, 2011, 968283.
Kappl, R.; Winkler, M.W. New limits on Dark Matter from Super-Kamiokande. Nucl. Phys. B 2011, 850,
505521.
Arina, C.; Hamann, J.; Wong, Y.Y.Y. A Bayesian view of the current status of Dark Matter direct searches.
Porter, T.A.; Johnson, R.P.; Graham, P.W. Dark matter searches with astroparticle data. Annu. Rev.
Astron. Astrophys. 2011, 49, 155194.
Calore, F.; de Romeri, V.; Donato, F. Conservative upper limits on WIMP annihilation cross section from
Fermi-LAT rays. Phys. Rev. D 2012, 85, 023004.
Klasen, M.; Pohl, M.; Sigl, G. Indirect and direct search for Dark Matter. Prog. Part. Nucl. Phys. 2015, 85,
132.
Mayet, F.; Green, A.M.; Battat, J.B.R.; Billard, J.; Bozorgnia, N.; Gelmini, G.B.; Gondolo, P.; Kavanagh, B.J.;
Lee, S.K.; Loomba, D.; et al. A review of the discovery reach of directional Dark Matter detection. Phys. Rep.
2016, 627, 149.
Jungman, G.; Kamionkowski, M.; Griest, K. Supersymmetric Dark Matter. Phys. Rep. 1996, 267, 195373.
IceCube Collaboration. Search for Dark Matter annihilations in the Sun with the 79-string IceCube detector.
Phys. Rev. Lett. 2013, 110, 131302.
Aguirre, A.; Schaye, J.; Quataert, E. Problems for modified Newtonian dynamics in clusters and the Ly
forest? Astrophys. J. 2001, 561, 550558.
Pointecouteau, E.; Silk, J. New constraints on modified Newtonian dynamics from galaxy clusters. Mon. Not.
R. Astron. Soc. 2005, 364, 654658.
Contaldi, C.R.; Wiseman, T.; Withers, B. TeVeS gets caught on caustics. Phys. Rev. D 2008, 78, 044034.
Lue, A.; Starkman, G.D. Squeezing MOND into a cosmological scenario. Phys. Rev. Lett. 2004, 92, 131102.
Dodelson, S. The real problem with MOND. Int. J. Mod. Phys. D 2011, 20, 27492753.
Ferreras, I.; Mavromatos, N.E.; Sakellariadou, M.; Yusaf, M.F. Confronting MOND and TeVeS with strong
gravitational lensing over galactic scales: An extended survey. Phys. Rev. D 2012, 86, 083507.
Clowe, D.; Bradac, M.; Gonzalez, A.H.; Markevitch, M.; Randall, S.W.; Jones, C.; Zaritsky, D. A direct
empirical proof of the existence of Dark Matter. Astrophys. J. 2006, 648, L109L113.
Peebles, P.J.; Ratra, B. The cosmological constant and dark energy. Rev. Mod. Phys. 2003, 75, 559606.
Slipher, V.M. Nebulae. Proc. Am. Philos. Soc. 1917, 56, 403409.
Paal, G.; Horvath, I.; Lukacs, B. Inflation and compactification from galaxy redshifts? Astrophys. Space Sci.
1992, 191, 107124.
Krauss, L.M. The end of the age problem, and the case for a cosmological constant revisited. Astrophys. J.
1998, 501, 461466.
269.
270.
271.
272.
273.
274.
275.
276.
277.
278.
279.
280.
281.
282.
283.
284.
285.
286.
287.
288.
289.
290.
291.
292.
293.
294.
295.
72 of 82
Rsnen, S. Structure formation as an alternative to dark energy and modified gravity. In Proceedings of
the Dark Energy and Dark Matter: Observations, Experiments and Theories, Lyons, France, 711 July 2008;
Volume 36, pp. 6374.
Blasone, M.; Capolupo, A.; Capozziello, S.; Carloni, S.; Vitiello, G. Neutrino mixing contribution to the
cosmological constant. Phys. Lett. A 2004, 323, 182189.
Capolupo, A.; Capozziello, S.; Vitiello, G. Neutrino mixing as a source of dark energy. Phys. Lett. A 2007,
363, 5356.
Riess, A.G.; Filippenko, A.V.; Challis, P.; Clocchiatti, A.; Diercks, A.; Garnavich, P.M.; Gilliland, R.L.;
Hogan, C.J.; Jha, S.; Kirshner, R.P.; et al. Observational evidence from supernovae for an accelerating
universe and a cosmological constant. Astron. J. 1998, 116, 10091038.
Perlmutter, S.; Aldering, G.; Goldhaber, G.; Knop, R.A.; Nugent, P.; Castro, P.G.; Deustua, S.; Fabbro, S.;
Goobar, A.; Groom, D.E.; et al. The supernova cosmology project. Measurements of Omega and Lambda
from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565586.
Sahni, V.; Starobinsky, A. The case for a positive cosmological -term. Int. J. Mod. Phys. D 2000, 9, 373443.
Turner, M.S. Dark Matter and Dark Energy in the universe. In The Third Stromlo Symposium: The Galactic Halo;
Gibson, B.K., Axelrod, R.S., Putman, M.E., Eds.; Astronomical Society of the Pacific Conference Series;
The Astronomical Society of the Pacific: San Francisco, CA, USA 1999; Volume 165, p. 431.
Eisenstein, D.J.; Zehavi, I.; Hogg, D.W.; Scoccimarro, R.; Blanton, M.R.; Nichol, R.C.; Scranton, R.; Seo, H.J.;
Tegmark, M.; Zheng, Z.; et al. Detection of the baryon acoustic peak in the large-scale correlation function
of SDSS luminous red galaxies. Astrophys. J. 2005, 633, 560574.
Frieman, J.A.; Turner, M.S.; Huterer, D. Dark Energy and the Accelerating Universe. Annu. Rev.
Astier, P.; Pain, R. Observational evidence of the accelerated expansion of the universe. C. R. Phys. 2012, 13,
521538.
Enqvist, K. Lemaitre Tolman Bondi model and accelerating expansion. Gen. Relativ. Gravit. 2008, 40,
451466.
Sotiriou, T.P.; Faraoni, V. f ( R) theories of gravity. Rev. Mod. Phys. 2010, 82, 451497.
Felice, A.D.; Tsujikawa, S. f ( R) theories. Living Rev. Relativ. 2010, 13, 3.
Capozziello, S.; de Laurentis, M. Extended theories of gravity. Phys. Rep. 2011, 509, 167321.
Khoury, J.; Weltman, A. Chameleon cosmology. Phys. Rev. D 2004, 69, 044026.
Durrer, R.; Maartens, R. Dark energy and dark gravity: Theory overview. Gen. Relativ. Gravit. 2008, 40,
301328..
Linder, E.V. Exploring the expansion history of the universe. Phys. Rev. Lett. 2003, 90, 091301.
Chevallier, M.; Polarski, D. Accelerating Universes with scaling Dark Matter. Int. J. Mod. Phys. D 2001, 10,
213223.
Linden, S.; Virey, J.M. Test of the Chevallier-Polarski-Linder parametrization for rapid dark energy equation
of state transitions. Phys. Rev. D 2008, 78, 023526.
Wang, Y.; Freese, K. Probing dark energy using its density instead of its equation of state. Phys. Lett. B
2006, 632, 449452.
Johri, V.B.; Rath, P.K. Parametrization of the dark energy equation of state. Int. J. Mod. Phys. D 2007, 16,
15811591.
Avelino, P.P.; Martins, C.J.A.P.; Nunes, N.J.; Olive, K.A. Reconstructing the dark energy equation of state
with varying couplings. Phys. Rev. D 2006, 74, 083508.
Scherrer, R.J. Mapping the Chevallier-Polarski-Linder parametrization onto physical dark energy models.
Phys. Rev. D 2015, 92, 043001.
Chongchitnan, S.; King, L. Imprints of dynamical dark energy on weak-lensing measurements. Mon. Not.
R. Astron. Soc. 2010, 407, 19891997.
Astashenok, A.V.; Odintsov, S.D. Confronting dark energy models mimicking CDM epoch with
observational constraints: Future cosmological perturbations decay or future Rip? Phys. Lett. B 2013, 718,
11941202.
Debono, I. Bayesian model selection for dark energy using weak lensing forecasts. Mon. Not. R. Astron. Soc.
2014, 437, 887897.
Bardeen, J.M. Gauge-invariant cosmological perturbations. Phys. Rev. D 1980, 22, 18821905.
296.
297.
298.
299.
300.
301.
302.
303.
304.
305.
306.
307.
308.
309.
310.
311.
312.
313.
314.
315.
316.
317.
318.
319.
320.
321.
322.
323.
73 of 82
Lifshitz, E.M. On the gravitational stability of the expanding universe. J. Phys. (USSR) 1946, 46, 587602.
Bartelmann, M.; Schneider, P. Weak gravitational lensing. Phys. Rep. 2001, 340, 291472.
Jeans, J.H. The stability of a spherical nebula. Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci. 1902, 199,
153.
Silk, J. Cosmic black-body radiation and galaxy formation. Astrophys. J. 1968, 151, 459471.
Bond, J.R.; Efstathiou, G.; Silk, J. Massive neutrinos and the large-scale structure of the universe.
Phys. Rev. Lett. 1980, 45, 19801984.
Silk, J.; Mamon, G.A. The current status of galaxy formation. Res. Astron. Astrophys. 2012, 12, 917946.
Linder, E.V.; Jenkins, A. Cosmic structure growth and dark energy. Mon. Not. R. Astron. Soc. 2003,
346, 573583.
Bernardeau, F.; Colombi, S.; Gaztaaga, E.; Scoccimarro, R. Large-scale structure of the Universe and
cosmological perturbation theory. Phys. Rep. 2002, 367, 1248.
Percival, W.J. Cosmological structure formation in a homogeneous dark energy background.
Hamilton, A.J.S.; Kumar, P.; Lu, E.; Matthews, A. Reconstructing the primordial spectrum of fluctuations
of the universe from the observed nonlinear clustering of galaxies. Astrophys. J. 1991, 374, L1L4.
Jain, B.; Mo, H.J.; White, S.D.M. The evolution of correlation functions and power spectra in gravitational
clustering. Mon. Not. R. Astron. Soc. 1995, 276, L25L29.
Peacock, J.A.; Dodds, S.J. Non-linear evolution of cosmological power spectra. Mon. Not. R. Astron. Soc.
1996, 280, L19L26.
Smith, R.E.; Peacock, J.A.; Jenkins, A.; White, S.D.M.; Frenk, C.S.; Pearce, F.R.; Thomas, P.A.; Efstathiou, G.;
Couchman, H.M.P. Stable clustering, the halo model and non-linear cosmological power spectra. Mon. Not.
R. Astron. Soc. 2003, 341, 13111332.
Seljak, U. Analytic model for galaxy and Dark Matter clustering. Mon. Not. R. Astron. Soc. 2000, 318,
203213.
Peacock, J.A.; Smith, R.E. Halo occupation numbers and galaxy bias. Mon. Not. R. Astron. Soc. 2000, 318,
11441156.
Huterer, D.; Takada, M. Calibrating the nonlinear matter power spectrum: Requirements for future weak
lensing surveys. Astropart. Phys. 2005, 23, 369376.
Van Daalen, M.P.; Schaye, J.; Booth, C.M.; Dalla Vecchia, C. The effects of galaxy formation on the matter
power spectrum: A challenge for precision cosmology. Mon. Not. R. Astron. Soc. 2011, 415, 36493665.
Semboloni, E.; Hoekstra, H.; Schaye, J.; van Daalen, M.P.; McCarthy, I.G. Quantifying the effect of baryon
physics on weak lensing tomography. Mon. Not. R. Astron. Soc. 2011, 417, 20202035.
Bird, S.; Viel, M.; Haehnelt, M.G. Massive neutrinos and the non-linear matter power spectrum. Mon. Not.
R. Astron. Soc. 2012, 420, 25512561.
Hearin, A.P.; Zentner, A.R.; Ma, Z. General requirements on matter power spectrum predictions for
cosmology with weak lensing tomography. J. Cosmol. Astropart. Phys. 2012, 4, 034.
Takahashi, R.; Sato, M.; Nishimichi, T.; Taruya, A.; Oguri, M. Revising the halofit model for the nonlinear
matter power spectrum. Astrophys. J. 2012, 761, 152.
Villaescusa-Navarro, F.; Viel, M.; Datta, K.K.; Choudhury, T.R. Modeling the neutral hydrogen distribution
in the post-reionization Universe: Intensity mapping. J. Cosmol. Astropart. Phys. 2014, 9, 050.
Bentivegna, E.; Bruni, M. Effects of nonlinear inhomogeneity on the cosmic expansion with numerical
relativity. Phys. Rev. Lett. 2016, 116, 251302.
Mertens, J.B.; Giblin, J.T.; Starkman, G.D. Integration of inhomogeneous cosmological spacetimes in the
BSSN formalism. Phys. Rev. D 2016, 93, 124059.
Giblin, J.T.; Mertens, J.B.; Starkman, G.D. Departures from the Friedmann-Lemaitre-Robertston-Walker
cosmological model in an inhomogeneous universe: A numerical examination. Phys. Rev. Lett. 2016, 116,
251301.
Harrison, E.R. Fluctuations at the threshold of classical cosmology. Phys. Rev. D 1970, 1, 27262730.
Zeldovich, Y.B. A hypothesis, unifying the structure and the entropy of the Universe. Mon. Not. R.
Astron. Soc. 1972, 160, 13.
Peebles, P.J.E.; Yu, J.T. Primeval adiabatic perturbation in an expanding universe. Astrophys. J. 1970, 162,
815836.
324.
325.
326.
327.
328.
329.
330.
331.
332.
333.
334.
335.
336.
337.
338.
339.
340.
341.
342.
343.
344.
345.
346.
347.
348.
349.
74 of 82
Mukherjee, P.; Wang, Y. Model-independent reconstruction of the primordial power spectrum from
Wilkinson Microwave Anistropy Probe data. Astrophys. J. 2003, 599, 16.
Trotta, R. Forecasting the Bayes factor of a future observation. Mon. Not. R. Astron. Soc. 2007, 378, 819824.
Trotta, R. Applications of Bayesian model selection to cosmological parameters. Mon. Not. R. Astron. Soc.
2007, 378, 7282.
Bridges, M.; Feroz, F.; Hobson, M.P.; Lasenby, A.N. Bayesian optimal reconstruction of the primordial
power spectrum. Mon. Not. R. Astron. Soc. 2009, 400, 10751084.
Vzquez, J.A.; Bridges, M.; Hobson, M.P.; Lasenby, A.N. Model selection applied to reconstruction of the
Primordial Power Spectrum. J. Cosmol. Astropart. Phys. 2012, 6, 006.
Hoyle, C.D.; Kapner, D.J.; Heckel, B.R.; Adelberger, E.G.; Gundlach, J.H.; Schmidt, U.; Swanson, H.E.
Submillimeter tests of the gravitational inverse-square law. Phys. Rev. D 2004, 70, 042004.
Kapner, D.J.; Cook, T.S.; Adelberger, E.G.; Gundlach, J.H.; Heckel, B.R.; Hoyle, C.D.; Swanson, H.E. Tests of
the gravitational inverse-square law below the dark-energy length scale. Phys. Rev. Lett. 2007, 98, 021101.
Baker, T.; Psaltis, D.; Skordis, C. Linking tests of gravity on all scales: From the strong-field regime to
cosmology. Astrophys. J. 2015, 802, 63.
Eddington, A.S. The Mathematical Theory of Relativity; Cambridge University Press: Cambridge, UK, 1923.
Shapiro, I.I. Solar system tests of general relativity: Recent results and present plans. In General Relativity
and Gravitation, 1989; Ashby, N., Bartlett, D.F., Wyss, W., Eds.; Cambridge University Press: Cambridge,
UK, 2005; p. 313.
Williams, J.G.; Turyshev, S.G.; Boggs, D.H. Progress in lunar laser ranging tests of relativistic gravity.
Phys. Rev. Lett. 2004, 93, 261101.
Williams, J.G.; Turyshev, S.G.; Boggs, D.H. Lunar laser ranging tests of the equivalence principle with the
Earth and Moon. Int. J. Mod. Phys. D 2009, 18, 11291175.
Shapiro, S.S.; Davis, J.L.; Lebach, D.E.; Gregory, J.S. Measurement of the solar gravitational deflection of
radio waves using geodetic very-long-baseline interferometry data, 19791999. Phys. Rev. Lett. 2004, 92,
121101.
Pitjeva, E.V.; Pitjev, N.P. Relativistic effects and Dark Matter in the Solar system from observations of
planets and spacecraft. Mon. Not. R. Astron. Soc. 2013, 432, 34313437.
Fienga, A.; Laskar, J.; Kuchynka, P.; Manche, H.; Desvignes, G.; Gastineau, M.; Cognard, I.; Theureau, G.
The INPOP10a planetary ephemeris and its applications in fundamental physics. Celest. Mech. Dyn. Astron.
2011, 111, 363385.
Iorio, L.; Lichtenegger, H.I.M.; Ruggiero, M.L.; Corda, C. Phenomenology of the Lense-Thirring effect in
the solar system. Astrophys. Space Sci. 2011, 331, 351395.
Hughes, V.W.; Robinson, H.G.; Beltran-Lopez, V. Upper Limit for the anisotropy of inertial mass from
nuclear resonance experiments. Phys. Rev. Lett. 1960, 4, 342344.
Drever, R.W.P. A search for anisotropy of inertial mass using a free precession technique. Philos. Mag. 1961,
6, 683687.
Allmendinger, F.; Heil, W.; Karpuk, S.; Kilian, W.; Scharth, A.; Schmidt, U.; Schnabel, A.; Sobolev, Y.;
Tullney, K. New limit on Lorentz-invariance- and CPT-violating neutron spin interactions using a
free-spin-precession He3-Xe129 comagnetometer. Phys. Rev. Lett. 2014, 112, 110801.
Williams, J.G.; Turyshev, S.G.; Boggs, D.H. Lunar laser ranging tests of the equivalence principle.
Class. Quantum Gravity 2012, 29, 184004.
Delva, P.; Hees, A.; Bertone, S.; Richard, E.; Wolf, P. Test of the gravitational redshift with stable clocks in
eccentric orbits: Application to Galileo satellites 5 and 6. Class. Quantum Gravity 2015, 32, 232003.
Iorio, L. Gravitational anomalies in the solar system? Int. J. Mod. Phys. D 2015, 24, 1530015.
Uzan, J.P. Varying constants, gravitation and cosmology. Living Rev. Relativ. 2011, 14, 2.
Flambaum, V.V. Variation of the fundamental constants: Theory and observations. Int. J. Mod. Phys. A
2007, 22, 49374950.
Lea, S.N. Limits to time variation of fundamental constants from comparisons of atomic frequency
standards. Rep. Prog. Phys. 2007, 70, 14731523.
Rich, J. Which fundamental constants for cosmic microwave background and baryon-acoustic oscillation?
Astron. Astrophys. 2015, 584, A69.
350.
351.
352.
353.
354.
355.
356.
357.
358.
359.
360.
361.
362.
363.
364.
365.
366.
367.
368.
369.
370.
371.
372.
373.
374.
375.
75 of 82
Barrow, J.D.; Sandvik, H.B.; Magueijo, J.A. Behavior of varying-alpha cosmologies. Phys. Rev. D 2002, 65,
063504.
Barrow, J.D.; Graham, A.A.H. General dynamics of varying-alpha universes. Phys. Rev. D 2013, 88, 103513.
Fujii, Y.; Iwamoto, A.; Fukahori, T.; Ohnuki, T.; Nakagawa, M.; Hidaka, H.; Oura, Y.; Mller, P. The nuclear
interaction at Oklo 2 billion years ago. Nucl. Phys. B 2000, 573, 377401.
Uzan, J.P. The fundamental constants and their variation: Observational and theoretical status.
Rev. Mod. Phys. 2003, 75, 403455.
Lamoreaux, S.K.; Torgerson, J.R. Neutron moderation in the Oklo natural reactor and the time variation
of . Phys. Rev. D 2004, 69, 121701.
Khatri, R.; Wandelt, B.D. 21-cm radiation: A new probe of variation in the fine-structure constant.
Phys. Rev. Lett. 2007, 98, 111301.
Nakashima, M.; Ichikawa, K.; Nagata, R.; Yokoyama, J. Constraining the time variation of the coupling
constants from cosmic microwave background: Effect of QCD. J. Cosmol. Astropart. Phys. 2010, 2010, 030.
Martins, C.J.A.P.; Menegoni, E.; Galli, S.; Mangano, G.; Melchiorri, A. Varying couplings in the early
universe: Correlated variations of and G. Phys. Rev. D 2010, 82, 023532.
Anderson, J.D.; Schubert, G.; Trimble, V.; Feldman, M.R. Measurements of Newtons gravitational constant
and the length of day. Europhys. Lett. 2015, 110, 10002.
Pitkin, M. Comment on Measurements of Newtons gravitational constant and the length of day.
Europhys. Lett. 2015, 111, 30002.
Anderson, J.D.; Schubert, G.; Trimble, V.; Feldman, M.R. Reply to the comment by M. Pitkin. Europhys. Lett.
2015, 111, 30003.
Iorio, L. Does Newtons gravitational constant vary sinusoidally with time? Orbital motions say no.
Class. Quantum Gravity 2016, 33, 045004.
Feldman, M.R.; Anderson, J.D.; Schubert, G.; Trimble, V.; Kopeikin, S.M.; Lmmerzahl, C. Deep space
experiment to measure G. Class. Quantum Gravity 2016, 33, 125013.
Lahav, O. Observational tests of FRW world models. Class. Quantum Gravity 2002, 19, 35173526.
Hansen, F.K.; Banday, A.J.; Grski, K.M. Testing the cosmological principle of isotropy: Local
power-spectrum estimates of the WMAP data. Mon. Not. R. Astron. Soc. 2004, 354, 641665.
Schwarz, D.J.; Bacon, D.; Chen, S.; Clarkson, C.; Huterer, D.; Kunz, M.; Maartens, R.; Raccanelli, A.;
Rubart, M.; Starck, J.L. Testing foundations of modern cosmology with SKA all-sky surveys. In Proceedings
of the Advancing Astrophysics with the Square Kilometre Array (AASKA14), Giardini Naxos, Italy,
913 June 2014; p. 32.
Segre, M. Galileo, Viviani and the tower of Pisa. Stud. Hist. Philos. Sci. Part A 1989, 20, 435451.
Etvos, R. ber die Anziehung der Erde auf Verschiedene Substanzen. Math. Naturwissenschaft. Ber. Ung.
1890, 8, 6568.
Gundlach, J.H.; Smith, G.L.; Adelberger, E.G.; Heckel, B.R.; Swanson, H.E. Short-range test of the
equivalence principle. Phys. Rev. Lett. 1997, 78, 25232526.
Schlamminger, S.; Choi, K.Y.; Wagner, T.A.; Gundlach, J.H.; Adelberger, E.G. Test of the equivalence
principle using a rotating torsion balance. Phys. Rev. Lett. 2008, 100, 041101.
Adelberger, E.; Gundlach, J.; Heckel, B.; Hoedl, S.; Schlamminger, S. Torsion balance experiments:
A low-energy frontier of particle physics. Prog. Part. Nucl. Phys. 2009, 62, 102134.
Wagner, T.A.; Schlamminger, S.; Gundlach, J.H.; Adelberger, E.G. Torsion-balance tests of the weak
equivalence principle. Class. Quantum Gravity 2012, 29, 184002.
Dickey, J.O.; Bender, P.L.; Faller, J.E.; Newhall, X.X.; Ricklefs, R.L.; Ries, J.G.; Shelus, P.J.; Veillet, C.;
Whipple, A.L.; Wiant, J.R.; et al. Lunar Laser ranging: A continuing legacy of the Apollo program. Science
1994, 265, 482490.
Murphy, T.W. Lunar laser ranging: The millimeter challenge. Rep. Prog. Phys. 2013, 76, 076901.
Pearlman, M.R.; Degnan, J.J.; Bosworth, J.M. The international laser ranging service. Adv. Space Res. 2002,
30, 135143.
Appleby, G.; Rodrguez, J.; Altamimi, Z. Assessment of the accuracy of global geodetic satellite laser
ranging observations and estimated impact on ITRF scale: Estimation of systematic errors in LAGEOS
observations 19932014. J. Geod. 2016, doi:10.1007/s00190-016-0929-2.
376.
377.
378.
379.
380.
381.
382.
383.
384.
385.
386.
387.
388.
389.
390.
391.
392.
393.
394.
395.
396.
397.
398.
399.
400.
401.
76 of 82
Degnan, J.J. Laser transponders for high-accuracy interplanetary laser ranging and time transfer. In Lasers,
Clocks and Drag-Free Control: Exploration of Relativistic Gravity in Space; Dittus, H., Lammerzahl, C.,
Turyshev, S.G., Eds.; Astrophysics and Space Science Library; Springer: Berlin/Heidelberg, Germany, 2008;
Volume 349, p. 231.
Iorio, L. Effects of standard and modified gravity on interplanetary ranges. Int. J. Mod. Phys. D 2011, 20,
181232.
Dirkx, D.; Noomen, R.; Visser, P.N.A.M.; Bauer, S.; Vermeersen, L.L.A. Comparative analysis of one- and
two-way planetary laser ranging concepts. Planet. Space Sci. 2015, 117, 159176.
Smith, D.; Zuber, M.; Sun, X.; Neumann, G.; Cavanaugh, J.; McGarry, J.; Zagwodzki, T. Two-way laser link
over interplanetary distance. Science 2006, 311, 53.
Chen, Y.; Birnbaum, K.M.; Hemmati, H. Active laser ranging over planetary distances with millimeter
accuracy. Appl. Phys. Lett. 2013, 102, 241107.
Margot, J.L.; Giorgini, J.D. Probing general relativity with radar astrometry in the inner solar system.
Proc. Int. Astron. Union 2010, 261, 183188.
Fienga, A.; Laskar, J.; Manche, H.; Gastineau, M. Tests of GR with INPOP15a Planetary Ephemerides:
Estimations of Possible Supplementary Advances of Perihelia for Mercury and Saturn. 2016,
arXiv:1601.00947.
Iorio, L. The Solar Lense-Thirring Effect: Perspectives for a Future Measurement. 2016, arXiv:1601.01382.
Anderson, J.D.; Laing, P.A.; Lau, E.L.; Liu, A.S.; Nieto, M.M.; Turyshev, S.G. Study of the anomalous
acceleration of Pioneer 10 and 11. Phys. Rev. D 2002, 65, 082004.
Nieto, M.M.; Anderson, J.D. Using early data to illuminate the Pioneer anomaly. Class. Quantum Gravity
2005, 22, 53435354.
Turyshev, S.G.; Toth, V.T. The pioneer anomaly. Living Rev. Relativ. 2010, 13, 4.
Nieto, M.M. New Horizons and the onset of the Pioneer anomaly. Phys. Lett. B 2008, 659, 483485.
Iorio, L. Perspectives on effectively constraining the location of a massive trans-Plutonian object with the
New Horizons spacecraft: A sensitivity analysis. Celest. Mech. Dyn. Astron. 2013, 116, 357366.
Iorio, L.; Ruggiero, M.L.; Radicella, N.; Saridakis, E.N. Constraining the Schwarzschild-de Sitter solution in
models of modified gravity. Phys. Dark Universe 2016, 13, 111120.
Damour, T.; Taylor, J.H. Strong-field tests of relativistic gravity and binary pulsars. Phys. Rev. D 1992, 45,
18401868.
Lyne, A.G.; Burgay, M.; Kramer, M.; Possenti, A.; Manchester, R.N.; Camilo, F.; McLaughlin, M.A.;
Lorimer, D.R.; DAmico, N.; Joshi, B.C.; et al. A double-pulsar system: A rare laboratory for relativistic
gravity and plasma physics. Science 2004, 303, 11531157.
Kramer, M.; Wex, N. The double pulsar system: A unique laboratory for gravity. Class. Quantum Gravity
2009, 26, 073001.
McGaugh, S.S.; Schombert, J.M.; Bothun, G.D.; de Blok, W.J.G. The baryonic tully-fisher relation.
Astrophys. J. 2000, 533, L99L102.
Famaey, B.; McGaugh, S.S. Modified Newtonian Dynamics (MOND): Observational phenomenology and
relativistic extensions. Living Rev. Relativ. 2012, 15, 10.
Adams, F.C.; Laughlin, G. Relativistic effects in extrasolar planetary systems. Int. J. Mod. Phys. D 2006, 15,
21332140.
Adams, F.C.; Laughlin, G. Effects of secular interactions in extrasolar planetary systems. Astrophys. J. 2006,
649, 9921003.
Adams, F.C.; Laughlin, G. Long-term evolution of close planets including the effects of secular interactions.
Astrophys. J. 2006, 649, 10041009.
Iorio, L. Are we far from testing general relativity with the transitting extrasolar planet HD 209458b Osiris?
New Astron. 2006, 11, 490494.
Jordn, A.; Bakos, G.. Observability of the general relativistic precession of periastra in exoplanets.
Astrophys. J. 2008, 685, 543552.
Pl, A.; Kocsis, B. Periastron precession measurements in transiting extrasolar planetary systems at the
level of general relativity. Mon. Not. R. Astron. Soc. 2008, 389, 191198.
Iorio, L. Classical and relativistic node precessional effects in WASP-33b and perspectives for detecting
them. Astrophys. Space Sci. 2011, 331, 485496.
402.
403.
404.
405.
406.
407.
408.
409.
410.
411.
412.
413.
414.
415.
416.
417.
418.
419.
420.
421.
422.
423.
424.
425.
426.
427.
77 of 82
Iorio, L. Classical and relativistic long-term time variations of some observables for transiting exoplanets.
Iorio, L. Accurate characterization of the stellar and orbital parameters of the exoplanetary system WASP-33
b from orbital dynamics. Mon. Not. R. Astron. Soc. 2016, 455, 207213.
Iorio, L. Post-Keplerian corrections to the orbital periods of a two-body system and their measurability.
Henry, R.C. Kretschmann scalar for a kerr-newman black hole. Astrophys. J. 2000, 535, 350353.
Will, C.M. Focus Issue: Gravity Probe B. Class. Quantum Gravity 2015, 32, 220301.
Everitt, C.W.F.; Muhlfelder, B.; DeBra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Silbergleit, A.S.; Acworth, E.B.;
Adams, M.; Adler, R.; Bencze, W.J.; et al. The Gravity Probe B test of general relativity. Class. Quantum
Gravity 2015, 32, 224001.
Hulse, R.A.; Taylor, J.H. Discovery of a pulsar in a binary system. Astrophys. J. 1975, 195, L51L53.
Psaltis, D. Probes and tests of strong-field gravity with observations in the electromagnetic spectrum.
Living Rev. Relativ. 2008, 11, 9.
Psaltis, D. Two approaches to testing general relativity in the strong-field regime. J. Phys. Conf. Ser. 2009,
189, 012033.
Johannsen, T. Testing general relativity in the strong-field regime with observations of black holes in the
electromagnetic spectrum. Publ. Astron. Soc. Pac. 2012, 124, 11331134.
Kramer, M. Precision tests of theories of gravity using pulsars. In Thirteenth Marcel Grossmann Meeting:
On Recent Developments in Theoretical and Experimental General Relativity, Astrophysics and Relativistic Field
Theories; Rosquist, K., Ed.; World Scientific: Singapore, 2015; pp. 315332.
Weisberg, J.M.; Taylor, J.H. The Relativistic Binary Pulsar B1913+16: Thirty Years of Observations and
Analysis. Astron. Soc. Pac. Conf. Ser. 2005, 328, 2532.
Bailes, M.; Ord, S.M.; Knight, H.S.; Hotan, A.W. Self-consistency of relativistic observables with general
relativity in the white dwarf-neutron star binary PSR J1141-6545. Astrophys. J. 2003, 595, L49L52.
Stairs, I.H.; Thorsett, S.E.; Taylor, J.H.; Wolszczan, A. Studies of the relativistic binary pulsar PSR B1534+12.
I. Timing analysis. Astrophys. J. 2002, 581, 501508.
Kramer, M.; Lorimer, D.R.; Lyne, A.G.; McLaughlin, M.; Burgay, M.; DAmico, N.; Possenti, A.; Camilo, F.;
Freire, P.C.C.; Joshi, B.C.; et al. Testing GR with the double pulsar: Recent results. In Proceedings of
the 22nd Texas Symposium on Relativistic Astrophysics, Stanford, CA, USA , 1317 December 2004;
pp. 142148.
Kramer, M.; Stairs, I.H.; Manchester, R.N.; McLaughlin, M.A.; Lyne, A.G.; Ferdman, R.D.; Burgay, M.;
Lorimer, D.R.; Possenti, A.; DAmico, N.; et al. Tests of general relativity from timing the double pulsar.
Science 2006, 314, 97102.
Desai, S.; Kahya, E.O. Galactic one-way Shapiro delay to PSR B1937+21. Mod. Phys. Lett. A 2016, 31,
1650083.
Huang, Y.; Weisberg, J.M. Timing the relativistic binary pulsar PSR B1913+16. In Proceedings of the 228th
American Astronomical Society Meeting, San Diego, CA, USA, 1216 June 2016.
Iorio, L. Prospects for measuring the moment of inertia of pulsar J0737-3039A. New Astron. 2009, 14, 4043.
Kehl, M.S.; Wex, N.; Kramer, M.; Liu, K. Future Measurements of the Lense-Thirring Effect in the Double
Pulsar. 2016, arXiv:1605.00408.
Iorio, L. Constraints from orbital motions around the Earth of the environmental fifth-force hypothesis for
the OPERA superluminal neutrino phenomenology. J. High Energy Phys. 2012, 5, 73.
Deng, X.M.; Xie, Y. Yukawa effects on the clock onboard a drag-free satellite. Mon. Not. R. Astron. Soc. 2013,
431, 32363239.
Li, Z.W.; Yuan, S.F.; Lu, C.; Xie, Y. New upper limits on deviation from the inverse-square law of gravity in
the solar system: A Yukawa parameterization. Res. Astron. Astrophys. 2014, 14, 139143.
Khoury, J.; Weltman, A. Chameleon fields: Awaiting surprises for tests of gravity in space. Phys. Rev. Lett.
2004, 93, 171104.
Tully, R.B.; Fisher, J.R. A new method of determining distances to galaxies. Astron. Astrophys. 1977, 54,
661673.
Dunsby, P.; Goheer, N.; Osano, B.; Uzan, J.P. How close can an inhomogeneous universe mimic the
concordance model? J. Cosmol. Astropart. Phys. 2010, 6, 017.
428.
429.
430.
431.
432.
433.
434.
435.
436.
437.
438.
439.
440.
441.
442.
443.
444.
445.
446.
447.
448.
449.
450.
451.
452.
453.
454.
455.
456.
78 of 82
Uzan, J.P.; Bernardeau, F. Lensing at cosmological scales: A test of higher dimensional gravity. Phys. Rev. D
2001, 64, 083004.
Peebles, P.J.E. The Large-Scale Structure of the Universe; Princeton Series in Physics: Princeton, NJ, USA, 1980.
Rich, J. Fundamentals of Cosmology; Springer: Berlin, Germany, 2001.
Uzan, J.P. The acceleration of the universe and the physics behind it. Gen. Relativ. Gravit. 2007, 39, 307342.
Uzan, J.P. Tests of general relativity on astrophysical scales. Gen. Relativ. Gravit. 2010, 42, 22192246.
Lue, A.; Scoccimarro, R.; Starkman, G.D. Probing Newtons constant on vast scales:
Dvali-Gabadadze-Porrati gravity, cosmic acceleration, and large scale structure. Phys. Rev. D
2004, 69, 124015.
Song, Y.S. Large scale structure formation of the normal branch in the DGP brane world model. Phys. Rev. D
2008, 77, 124031.
Benabed, K.; Bernardeau, F. Testing quintessence models with large-scale structure growth. Phys. Rev. D
2001, 64, 083501.
Koivisto, T.; Mota, D.F. Gauss-Bonnet quintessence: Background evolution, large scale structure,
and cosmological constraints. Phys. Rev. D 2007, 75, 023518.
Tsujikawa, S. Quintessence: A review. Class. Quantum Gravity 2013, 30, 214003.
Baldi, M.; Pettorino, V.; Amendola, L.; Wetterich, C. Oscillating non-linear large-scale structures in growing
neutrino quintessence. Mon. Not. R. Astron. Soc. 2011, 418, 214229.
Koivisto, T.S.; Saridakis, E.N.; Tamanini, N. Scalar-fluid theories: Cosmological perturbations and
large-scale structure. J. Cosmol. Astropart. Phys. 2015, 2015, 047.
Schimd, C.; Uzan, J.P.; Riazuelo, A. Weak lensing in scalar-tensor theories of gravity. Phys. Rev. D 2005,
71, 083512.
Rodrguez-Meza, M.A. Power spectrum of large-scale structure cosmological models in the framework of
scalar-tensor theories. J. Phys. Conf. Ser. 2010, 229, 012063.
Goenner, H. Some remarks on the genesis of scalar-tensor theories. Gen. Relativ. Gravit. 2012, 44, 20772097.
Takushima, Y.; Terukina, A.; Yamamoto, K. Bispectrum of cosmological density perturbations in the most
general second-order scalar-tensor theory. Phys. Rev. D 2014, 89, 104007.
Dossett, J.N.; Ishak, M.; Parkinson, D.; Davis, T.M. Constraints and tensions in testing general relativity
from Planck and CFHTLenS data including intrinsic alignment systematics. Phys. Rev. D 2015, 92, 023003.
Di Valentino, E.; Melchiorri, A.; Silk, J. Cosmological hints of modified gravity? Phys. Rev. D 2016, 93,
023513.
Hu, B.; Liguori, M.; Bartolo, N.; Matarrese, S. Parametrized modified gravity constraints after Planck.
Phys. Rev. D 2013, 88, 123514.
Pettorino, V. Testing modified gravity with Planck: The case of coupled dark energy. Phys. Rev. D 2013,
88, 063519.
Planck Collaboration. Planck 2015 Results. XIV. Dark Energy and Modified Gravity. Astron. Astrophys.
2016, 594, A14.
Rsnen, S. Dark energy from back-reaction. J. Cosmol. Astropart. Phys. 2004, 2, 003.
Wiltshire, D.L. Gravitational energy as dark energy: Towards concordance cosmology without Lambda.
In Proceedings of the Dark Energy and Dark Matter: Observations, Experiments and Theories, Lyons,
France, 711 July 2008; Volume 36, pp. 9198.
Rsnen, S. Backreaction: Directions of progress. Class. Quantum Gravity 2011, 28, 164008.
Buchert, T.; Carfora, M.; Ellis, G.F.R.; Kolb, E.W.; MacCallum, M.A.H.; Ostrowski, J.J.; Rsnen, S.;
Roukema, B.F.; Andersson, L.; Coley, A.A.; et al. Is there proof that backreaction of inhomogeneities is
irrelevant in cosmology? Class. Quantum Gravity 2015, 32, 215021.
Rsnen, S.; Bolejko, K.; Finoguenov, A. New test of the Friedmann-Lematre-Robertson-Walker metric
using the distance sum rule. Phys. Rev. Lett. 2015, 115, 101301.
Rsnen, S.; Vliviita, J.; Kosonen, V. Testing distance duality with CMB anisotropies. J. Cosmol.
Astropart. Phys. 2016, 4, 050.
Boehm, C.; Rsnen, S. Violation of the FRW consistency condition as a signature of backreaction. J. Cosmol.
Buchert, T. On average properties of inhomogeneous fluids in general relativity: Dust cosmologies.
Gen. Relativ. Gravit. 2000, 32, 105126.
457.
458.
459.
460.
461.
462.
463.
464.
465.
466.
467.
468.
469.
470.
471.
472.
473.
474.
475.
476.
477.
478.
479.
480.
481.
482.
483.
484.
485.
79 of 82
Buchert, T. Dark Energy from structure: A status report. Gen. Relativ. Gravit. 2008, 40, 467527.
Banks, T. Supersymmetry breaking and the cosmological constant. Int. J. Mod. Phy. A 2014, 29, 1430010.
Arkani-Hamed, N.; Dimopoulos, S. Supersymmetric unification without low energy supersymmetry and
signatures for fine-tuning at the LHC. J. High Energy Phys. 2005, 6, 073.
Susskind, L. The anthropic landscape of string theory. In Universe or Multiverse? Carr, B., Ed.; Cambridge
University Press: Cambridge, UK, 2007.
Bauer, F.; Sol, J.; tefancic, H. Dynamically avoiding fine-tuning the cosmological constant: The Relaxed
Universe. J. Cosmol. Astropart. Phys. 2010, 12, 029.
Bombelli, L.; Lee, J.; Meyer, D.; Sorkin, R.D. Space-time as a causal set. Phys. Rev. Lett. 1987, 59, 521524.
Ahmed, M.; Dodelson, S.; Greene, P.B.; Sorkin, R. Everpresent . Phys. Rev. D 2004, 69, 103523.
Weinberg, S. Anthropic bound on the cosmological constant. Phys. Rev. Lett. 1987, 59, 26072610.
Garriga, J.; Linde, A.; Vilenkin, A. Dark energy equation of state and anthropic selection. Phys. Rev. D 2004,
69, 063521.
Armendariz-Picon, C.; Mukhanov, V.; Steinhardt, P.J. Essentials of k-essence. Phys. Rev. D 2001, 63, 103510.
Kamenshchik, A.; Moschella, U.; Pasquier, V. An alternative to quintessence. Phys. Lett. B 2001, 511,
265268.
de Putter, R.; Linder, E.V. Kinetic k-essence and quintessence. Astropart. Phys. 2007, 28, 263272.
Durrer, R.; Maartens, R. Dark Energy and Modified Gravity. In Dark Energy: Observational and Theoretical
Approaches; Ruiz-Lapuente, P., Ed.; Cambridge University Press: Cambridge, UK, 2008; pp. 4891.
Deffayet, C.; Dvali, G.; Gabadadze, G. Accelerated universe from gravity leaking to extra dimensions.
Phys. Rev. D 2002, 65, 044023.
Amendola, L.; Polarski, D.; Tsujikawa, S. Are f(R) Dark Energy Models Cosmologically Viable?
Phys. Rev. Lett. 2007, 98, 131302.
Turyshev, S.G.; Toth, V.T.; Kinsella, G.; Lee, S.C.; Lok, S.M.; Ellis, J. Support for the thermal origin of the
pioneer anomaly. Phys. Rev. Lett. 2012, 108, 241101.
Adler, C.G.; Coulter, B.L. Galileo and the Tower of Pisa experiment. Am. J. Phys. 1978, 46, 199201.
Berg, J.; Touboul, P.; Rodrigues, M.; The MICROSCOPE Team. Status of MICROSCOPE, a mission to test
the Equivalence Principle in space. J. Phys. Conf. Ser. 2015, 610, 012009.
Overduin, J.; Everitt, F.; Worden, P.; Mester, J. STEP and fundamental physics. Class. Quantum Gravity 2012,
29, 184012.
Nobili, A.M.; Comandi, G.L.; Doravari, S.; Bramanti, D.; Kumar, R.; Maccarrone, F.; Polacco, E.;
Turyshev, S.G.; Shao, M.; Lipa, J.; et al. Galileo Galilei (GG) a small satellite to test the equivalence
principle of Galileo, Newton and Einstein. Exp. Astron. 2009, 23, 689710.
Muller, H.; Peters, A.; Chu, S. A precision measurement of the gravitational redshift by the interference of
matter waves. Nature 2010, 463, 926929.
Bardeen, J.M.; Petterson, J.A. The lense-thirring effect and accretion disks around kerr black holes.
Astrophys. J. 1975, 195, L65.
Pound, R.V.; Rebka, G.A. Apparent weight of photons. Phys. Rev. Lett. 1960, 4, 337341.
Turyshev, S.G.; Yu, N.; Toth, V.T. General relativistic observables for the ACES experiment. Phys. Rev. D
2016, 93, 045027.
Vessot, R.F.C.; Levine, M.W.; Mattison, E.M.; Blomberg, E.L.; Hoffman, T.E.; Nystrom, G.U.; Farrel, B.F.;
Decher, R.; Eby, P.B.; Baugher, C.R. Test of relativistic gravitation with a space-borne hydrogen maser.
Phys. Rev. Lett. 1980, 45, 20812084.
Schiller, S.; Tino, G.M.; Gill, P.; Salomon, C.; Sterr, U.; Peik, E.; Nevsky, A.; Grlitz, A.; Svehla, D.;
Ferrari, G.; et al. Einstein Gravity ExplorerA medium-class fundamental physics mission. Exp. Astron.
2009, 23, 573610.
Antoniadis, J.; Freire, P.C.C.; Wex, N.; Tauris, T.M.; Lynch, R.S.; van Kerkwijk, M.H.; Kramer, M.; Bassa, C.;
Dhillon, V.S.; Driebe, T.; et al. A massive pulsar in a compact relativistic binary. Science 2013, 340, 448.
Iorio, L. Constraining the preferred-frame 1 , 2 parameters from solar system planetary precessions. Int. J.
Mod. Phys. D 2014, 23, 1450006.
Fienga, A.; Laskar, J.; Exertier, P.; Manche, H.; Gastineau, M. Numerical estimation of the sensitivity of
INPOP planetary ephemerides to general relativity parameters. Celest. Mech. Dyn. Astron. 2015, 123,
325349.
486.
487.
488.
489.
490.
491.
492.
493.
494.
495.
496.
497.
498.
499.
500.
501.
502.
503.
504.
505.
506.
507.
508.
509.
510.
80 of 82
Everitt, C.W.F.; Debra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Conklin, J.W.; Heifetz, M.I.; Keiser, G.M.;
Silbergleit, A.S.; Holmes, T.; Kolodziejczak, J.; et al. Gravity Probe B: Final results of a space experiment to
test general relativity. Phys. Rev. Lett. 2011, 106, 221101.
Kramer, M. Pulsars as probes of gravity and fundamental physics. Int. J. Mod. Phys. D 2016, 25, 1630029.
Bertotti, B.; Ciufolini, I.; Bender, P.L. New test of general relativityMeasurement of de Sitter geodetic
precession rate for lunar perigee. Phys. Rev. Lett. 1987, 58, 10621065.
Shapiro, I.I.; Reasenberg, R.D.; Chandler, J.F.; Babcock, R.W. Measurement of the de Sitter precession of the
moonA relativistic three-body effect. Phys. Rev. Lett. 1988, 61, 26432646.
Williams, J.G.; Newhall, X.X.; Dickey, J.O. Relativity parameters determined from lunar laser ranging.
Phys. Rev. D 1996, 53, 67306739.
Merkowitz, S.M. Tests of gravity using lunar laser ranging. Living Rev. Relativ. 2010, 13, 7.
Martini, M.; DellAgnello, S.; Currie, D.; Delle Monache, G.O.; Vittori, R.; Berardi, S.; Boni, A.; Cantone, C.;
Ciocci, E.; Lops, C.; et al. Moonlight: A new lunar laser ranging retroreflector and the lunar geodetic
precession. Acta Polytech. 2013, 53, 745750.
Kramer, M. Determination of the geometry of the PSR B1913 + 16 System by geodetic precession.
Astrophys. J. 1998, 509, 856860.
Breton, R.P.; Kaspi, V.M.; Kramer, M.; McLaughlin, M.A.; Lyutikov, M.; Ransom, S.M.; Stairs, I.H.;
Ferdman, R.D.; Camilo, F.; Possenti, A. Relativistic spin precession in the double pulsar. Science 2008,
321, 104107.
Lucchesi, D.M.; Anselmo, L.; Bassan, M.; Pardini, C.; Peron, R.; Pucacco, G.; Visco, M. Testing the
gravitational interaction in the field of the Earth via satellite laser ranging and the Laser Ranged Satellites
Experiment (LARASE). Class. Quantum Gravity 2015, 32, 155012.
Ingram, A.; van der Klis, M.; Middleton, M.; Done, C.; Altamirano, D.; Heil, L.; Uttley, P.; Axelsson, M.
A quasi-periodic modulation of the iron line centroid energy in the black hole binary H 1743-322. Mon. Not.
R. Astron. Soc. 2016, 461, 19671980.
Bosi, F.; Cella, G.; di Virgilio, A.; Ortolan, A.; Porzio, A.; Solimeno, S.; Cerdonio, M.; Zendri, J.P.;
Allegrini, M.; Belfi, J.; et al. Measuring gravitomagnetic effects by a multi-ring-laser gyroscope. Phys. Rev. D
2011, 84, 122002.
Stone, N.; Loeb, A. Observing lense-thirring precession in tidal disruption flares. Phys. Rev. Lett. 2012,
108, 061302.
Franchini, A.; Lodato, G.; Facchini, S. Lense-Thirring precession around supermassive black holes during
tidal disruption events. Mon. Not. R. Astron. Soc. 2016, 455, 19461956.
Martynov, D.V.; Hall, E.D.; Abbott, B.P.; Abbott, R.; Abbott, T.D.; Adams, C.; Adhikari, R.X.;
Anderson, R.A.; Anderson, S.B.; Arai, K.; et al. Sensitivity of the advanced LIGO detectors at the
beginning of gravitational wave astronomy. Phys. Rev. D 2016, 93, 112004.
LIGO Scientific Collaboration and Virgo Collaboration. GW151226: Observation of gravitational waves
from a 22-solar-mass binary black hole coalescence. Phys. Rev. Lett. 2016, 116, 241103.
Johannsen, T. Sgr A* and general relativity. Class. Quantum Gravity 2016, 33, 113001.
Soffel, M.; Wirrer, R.; Schastok, J.; Ruder, H.; Schneider, M. Relativistic effects in the motion of artificial
satellites: The oblateness of the central body I. Celest. Mech. 1987, 42, 8189.
Heimberger, J.; Soffel, M.; Ruder, H. Relativistic effects in the motion of artificial satellites: The oblateness
of the central body II. Celest. Mech. Dyn. Astron. 1989, 47, 205217.
Soffel, M. Relativity in Astrometry, Celestial Mechanics and Geodesy; Springer-Verlag: Berlin, Germany, 1989.
Iorio, L. Post-Newtonian direct and mixed orbital effects due to the oblateness of the central body. Int. J.
Mod. Phys. D 2015, 24, 1550067.
Iorio, L. A possible new test of general relativity with Juno. Class. Quantum Gravity 2013, 30, 195011.
Starkman, G.D.; Trotta, R.; Vaudrevange, P.M. The virtues of frugalityWhy cosmological observers
should release their data slowly. Mon. Not. R. Astron. Soc. 2010, 401, L15L18.
Domcke, V.; Pieroni, M.; Bintruy, P. Primordial gravitational waves for universality classes of pseudoscalar
inflation. J. Cosmol. Astropart. Phys. 2016, 6, 031.
Ito, A.; Soda, J. MHz gravitational waves from short-term anisotropic inflation. J. Cosmol. Astropart. Phys.
2016, 4, 035.
511.
512.
513.
514.
515.
516.
517.
518.
519.
520.
521.
522.
523.
524.
525.
526.
527.
528.
529.
530.
531.
532.
533.
534.
535.
536.
537.
81 of 82
Fidler, C.; Pettinari, G.W.; Beneke, M.; Crittenden, R.; Koyama, K.; Wands, D. The intrinsic B-mode
polarisation of the Cosmic Microwave Background. J. Cosmol. Astropart. Phys. 2014, 7, 011.
Namikawa, T.; Nagata, R. Non-Gaussian structure of B-mode polarization after delensing. J. Cosmol.
Abazajian, K.N.; Arnold, K.; Austermann, J.; Benson, B.A.; Bischoff, C.; Bock, J.; Bond, J.R.; Borrill, J.;
Calabrese, E.; Carlstrom, J.E.; et al. Neutrino physics from the cosmic microwave background and large
scale structure. Astropart. Phys. 2015, 63, 6680.
Trotta, R. Bayes in the sky: Bayesian inference and model selection in cosmology. Contemp. Phys. 2008, 49,
71104.
Hobson, M.P.; Jaffe, A.H.; Liddle, A.R.; Mukeherjee, P.; Parkinson, D. Bayesian Methods in Cosmology;
Cambridge University Press: Cambridge, UK, 2010.
Liddle, A.R. How many cosmological parameters? Mon. Not. R. Astron. Soc. 2004, 351, L49L53.
Metcalf, R.B.; Silk, J. On breaking cosmic degeneracy. Astrophys. J. 1998, 492, L1L4.
Debono, I.; Rassat, A.; Rfrgier, A.; Amara, A.; Kitching, T.D. Weak lensing forecasts for dark energy,
neutrinos and initial conditions. Mon. Not. R. Astron. Soc. 2010, 404, 110119.
Howlett, C.; Lewis, A.; Hall, A.; Challinor, A. CMB power spectrum parameter degeneracies in the era of
precision cosmology. J. Cosmol. Astropart. Phys. 2012, 4, 027.
Joyce, A.; Jain, B.; Khoury, J.; Trodden, M. Beyond the cosmological standard model. Phys. Rep. 2015,
568, 198.
Jaffe, A. H0 and odds on cosmology. Astrophys. J. 1996, 471, 24.
Mukherjee, P.; Parkinson, D.; Liddle, A.R. A nested sampling algorithm for cosmological model selection.
Astrophys. J. 2006, 638, L51L54.
Kilbinger, M.; Wraith, D.; Robert, C.P.; Benabed, K.; Capp, O.; Cardoso, J.F.; Fort, G.; Prunet, S.;
Bouchet, F.R. Bayesian model comparison in cosmology with Population Monte Carlo. Mon. Not. R.
Astron. Soc. 2010, 405, 23812390.
Wraith, D.; Kilbinger, M.; Benabed, K.; Capp, O.; Cardoso, J.; Fort, G.; Prunet, S.; Robert, C.P. Estimation
of cosmological parameters using adaptive importance sampling. Phys. Rev. D 2009, 80, 023507.
Del Pozzo, W.; Veitch, J.; Vecchio, A. Testing general relativity using Bayesian model selection: Applications
to observations of gravitational waves from compact binary systems. Phys. Rev. D 2011, 83, 082002.
Ellis, G.F.R. Dark Matter and Dark Energy Proposals: Maintaining Cosmology as a True Science; EAS Publications
Series; European Astronomical Society: Versoix, Switzerland, 2009; Volume 36, pp. 325336.
Ellis, G.F.R. Issues in the Philosophy of Cosmology. In Philosophy of Physics; Butterfield, J., Earman, J., Eds.;
Handbook of the Philosophy of Science; Elsevier: Amsterdam, The Netherlands, 2006; pp. 11831285.
Jackman, S. Bayesian Analysis for the Social Sciences; Wiley: Hoboken, NJ, USA, 2009.
Niemack, M.D.; Ade, P.A.R.; Aguirre, J.; Barrientos, F.; Beall, J.A.; Bond, J.R.; Britton, J.; Cho, H.M.; Das, S.;
Devlin, M.J.; et al. ACTPol: A polarization-sensitive receiver for the Atacama Cosmology Telescope.
Proc. SPIE 2010, 7741, 77411S.
Beynon, E.; Bacon, D.J.; Koyama, K. Weak lensing predictions for modified gravities at non-linear scales.
Harnois-Draps, J.; Munshi, D.; Valageas, P.; van Waerbeke, L.; Brax, P.; Coles, P.; Rizzo, L. Testing modified
gravity with cosmic shear. Mon. Not. R. Astron. Soc. 2015, 454, 27222735.
Carlstrom, J.E.; CMB-S4 Collaboration. The next generation ground-based CMB experiment, CMB-S4.
In Proceedings of the 228th AAS Meeting, San Diego, CA, USA, 1216 June 2016.
Santos, M.; Bull, P.; Alonso, D.; Camera, S.; Ferreira, P.; Bernardi, G.; Maartens, R.; Viel, M.;
Villaescusa-Navarro, F.; Abdalla, F.B.; et al. Cosmology from a SKA HI intensity mapping
survey. In Proceedings of the Advancing Astrophysics with the Square Kilometre Array (AASKA14),
Giardini Naxos, Italy, 913 June 2014.
Hall, A.; Bonvin, C.; Challinor, A. Testing general relativity with 21-cm intensity mapping. Phys. Rev. D
2013, 87, 064026.
Lahav, O.; Suto, Y. Measuring our Universe from Galaxy Redshift Surveys. Living Rev. Relativ. 2004, 7, 8.
Dup, F.X.; Rassat, A.; Starck, J.L.; Fadili, M.J. Measuring the integrated Sachs-Wolfe effect.
Astron. Astrophys. 2011, 534, A51.
Stabenau, H.F.; Jain, B. N-body simulations of alternative gravity models. Phys. Rev. D 2006, 74, 084007.
538.
539.
540.
541.
542.
543.
544.
545.
546.
547.
548.
549.
550.
551.
552.
553.
554.
555.
82 of 82
Adamek, J.; Daverio, D.; Durrer, R.; Kunz, M. General relativistic N-body simulations in the weak field
limit. Phys. Rev. D 2013, 88, 103527.
He, J.H.; Li, B.; Hawken, A.J. Effective Dark Matter power spectra in f ( R) gravity. Phys. Rev. D 2015,
92, 103508.
Winther, H.A.; Schmidt, F.; Barreira, A.; Arnold, C.; Bose, S.; Llinares, C.; Baldi, M.; Falck, B.;
Hellwing, W.A.; Koyama, K.; et al. Modified gravity N-body code comparison project. Mon. Not.
R. Astron. Soc. 2015, 454, 42084234.
Eingorn, M. First-order cosmological perturbations engendered by point-like masses. Astrophys. J. 2016,
825, 84.
Hahn, O.; Paranjape, A. General Relativistic Screening in Cosmological Simulations. 2016,
arXiv:1602.07699.
Kamdar, H.M.; Turk, M.J.; Brunner, R.J. Machine learning and cosmological simulationsI. Semi-analytical
models. Mon. Not. R. Astron. Soc. 2016, 455, 642658.
Kamdar, H.M.; Turk, M.J.; Brunner, R.J. Machine learning and cosmological simulationsII.
Hydrodynamical simulations. Mon. Not. R. Astron. Soc. 2016, 457, 11621179.
Smoot, G.F. See Saw Inflation/Dark Matter/Dark Energy/Baryogenesis. 2014, arXiv:1405.2776.
Rickles, D. A Second Superstring Revolution and the Future of String Theory. In A Brief History of String
Theory: From Dual Models to M-Theory; Springer: Berlin, Germany, 2014; pp. 207242.
Maartens, R.; Koyama, K. Brane-world gravity. Living Rev. Relativ. 2010, 13, 5.
Woodard, R.P. How far are we from the quantum theory of gravity? Rep. Prog. Phys. 2009, 72, 126002.
Norris, R.P. Data challenges for next-generation radio telescopes. In Proceedings of the Sixth IEEE
International Conference on eScience, Brisbane, Australia, 710 December 2010; pp. 2124.
Planck Collaboration. Planck 2015 Results. II. Low Frequency Instrument Data Processing. Astron. Astrophys.
2016, 594, A2.
Planck Collaboration. Planck 2015 Results. VIII. High Frequency Instrument Data Processing: Calibration
and Maps. Astron. Astrophys. 2016, 594, A8.
Juric, M.; Kantor, J.; Lim, K.; Lupton, R.H.; Dubois-Felsmann, G.; Jenness, T.; Axelrod, T.S.; Aleksic, J.;
Allsman, R.A.; AlSayyad, Y.; et al. The LSST Data Management System. 2015, arXiv:1512.07914.
Desai, S.; Mohr, J.J.; Henderson, R.; Kmmel, M.; Paech, K.; Wetzstein, M. CosmoDM and its application to
Pan-STARRS data. J. Instrum. 2015, 10, C06014.
Dodson, R.; Vinsen, K.; Wu, C.; Popping, A.; Meyer, M.; Wicenec, A.; Quinn, P.; van Gorkom, J.; Momjian, E.
Imaging SKA-scale data in three different computing environments. Astron. Comput. 2016, 14, 822.
Arviset, C.; Gonzlez, J.; Gutirrez, R.; Hernndez, J.; Salgado, J.; Segovia, J.C. The Gaia archive: VO
in action in the big data era. In Proceedings of the International Symposium Dedicated to the 50th
Anniversary of the Markarian Survey and the 10th Anniversary of the Armenian Virtual Observatory,
Byurakan, Armenia, 58 October 2015; Volume 505, p. 218.
c 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access

article distributed under the terms and conditions of the Creative Commons Attribution
(CC-BY) license (http://creativecommons.org/licenses/by/4.0/).

General Relativity and Cosmology: Unsolved Questions and Future Directions

Uploaded by

Copyright:

Available Formats

General Relativity and Cosmology: Unsolved Questions and Future Directions

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

General Relativity and Cosmology: Unsolved Questions and Future Directions

Uploaded by

Copyright:

Available Formats

Article

General Relativity and Cosmology: Unsolved

arXiv:1609.09781v1 [gr-qc] 28 Sep 2016

Academic Editors: Lorenzo Iorio and Elias C. Vagenas

Universe 2016, 2, 23; doi:10.3390/universe2040023

13.8 billion years

1.8 billion years

population III stars

200 million years

Time since Big Bang

Newtons formulation of the gravitational laws is expressed by the equations:

where the comma denotes a derivative, i.e., g, = x

where the covariant derivative is defined as

The quantity is the Lorentz factor:

it agrees with experiment

and g is constructed from and :

where ( TM ) is the energy-momentum tensor of matter and ( T ) is the energy-momentum tensor

Geometrically, it means that an infinitesimal geodesic parallelogram forms a non-closed loop.

Newton 1687 [20]

Nonrelativistic, implicit action at a distance

Nordstrm 1913 [57]

Fails to predict observed light detection

General Relativity 1915 [46]

Poincar 1890s1900s [31,118]

Fails; does not mesh with electromagnetism

Whitehead 1922 [59]

Violates LLI; contradiction by everyday observation of tides

Cartan 19221925 [98]

Kaluza-Klein 1920s [119,120]

Still viable; introduces matter spin

Birkhoff 1943 [61]

Fails Newtonian test; demands speed of

Milne 1948 [121]

Incomplete; no gravitational redshift prediction;

Thiry 1948 [122]

Belifante-Swihart 1957 [123]

BransDicke 1961 [63]

Viable for > 500

Violates LPI; predicts preferred-frame effects

Will-Nordtvedt 1972 [124]

Barker 1978 [125]

Rosen 1973 [126,127]

Contradicted by binary pulsar data

Rastall 1976 [128]

Contradicted by gravitational wave data

f ( R) models 1970s [129,130]

Consistent with Solar System tests;

MOND 1983 [131133]

DGP 2000 [134]

TeVeS 2004 [135]

Violates Equivalence Principle

Unlikely; extremely constrained by results on PPN

Violates EEP; contradicted by DickeBraginsky experiments

Viable but can only be significant at high energy regimes

Highly unstable [136]; ruled out by SDSS data [137]

5. From General Relativity to Standard Cosmology

ds2 = c d 2 R2 (t)[ dr2 + Sk0 (r )( d 2 + sin2 d2 )] ,

where Sk (r ) can be redefined as