Next Article in Journal
Balancing Asymmetric Dark Matter with Baryon Asymmetry and Dilution of Frozen Dark Matter by Sphaleron Transition
Previous Article in Journal
The Swampland Conjectures: A Bridge from Quantum Gravity to Particle Physics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Relativistic Cosmology with an Introduction to Inflation

by
Muhammad Zahid Mughal
1,
Iftikhar Ahmad
1 and
Juan Luis García Guirao
2,3,*
1
Department of Mathematics, University of Gujrat, Gujrat 50700, Pakistan
2
Department of Applied Mathematics and Statistics, Technical University of Cartagena, 30310 Cartagena, Spain
3
Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 23 June 2021 / Revised: 20 July 2021 / Accepted: 26 July 2021 / Published: 30 July 2021

Abstract

:
In this review article, the study of the development of relativistic cosmology and the introduction of inflation in it as an exponentially expanding early phase of the universe is carried out. We study the properties of the standard cosmological model developed in the framework of relativistic cosmology and the geometric structure of spacetime connected coherently with it. The geometric properties of space and spacetime ingrained into the standard model of cosmology are investigated in addition. The big bang model of the beginning of the universe is based on the standard model which succumbed to failure in explaining the flatness and the large-scale homogeneity of the universe as demonstrated by observational evidence. These cosmological problems were resolved by introducing a brief acceleratedly expanding phase in the very early universe known as inflation. The cosmic inflation by setting the initial conditions of the standard big bang model resolves these problems of the theory. We discuss how the inflationary paradigm solves these problems by proposing the fast expansion period in the early universe. Further inflation and dark energy in f R modified gravity are also reviewed.

1. Introduction

With the advent of general relativity in 1916, spacetime transformed itself into one of the four fundamental interactions of the universe and the geometrical structure attached to it was taken to demonstrate gravity in a dynamical way [1]. The force of gravity was replaced by the curvature of spacetime that is mirrored through the geometric structure of metric tensor g μ ν . The spacetime became an integral part of the universe and a dynamical medium where the whole phenomenal universe exists. Any solution of the field equations of general relativity entails a certain structural geometry of spacetime or just a spacetime that represents a universe itself, therefore determining a solution of the field equations is similar to coming across a specific model of the universe.
Cosmology studies the universe as a whole [2] encompassing its beginning in spacetime or as spacetime itself, its evolution, and its eventual ultimate fate. The history of cosmology dates back to ancient Greeks, Indians, and Iranians with its roots at that time in philosophy and religion. Before modern scientific cosmology emerges, it has been nurtured in the womb of Ibrahamic religions especially Judaism, Christianity, and Islam. Cosmology as modern science begins with the surfacing of general relativity when Einstein first himself put it to use to formulate a cosmological model of the universe mathematically. The model brought about a dynamic universe but was rendered to be static as there was no cosmological evidence of its contraction or expansion at that time [3]. Einstein’s static model was afterward proved to be inconsistent with cosmological observations and was discarded; however, its formulation as the first mathematical model based on the field equations of general relativity laid the foundational stone for the inception of modern relativistic cosmology as science.
Cosmology takes into account the largest scale of spacetime that is the causally connected maximal patch of the cosmos from the perspective of its origin, evolution, and futuristic eventual fate. It gives the universe a mathematical description as large as the cosmological observational parameters reveal and allow consequently. The modern relativistic cosmology was established on general relativity which brought forth the big bang model of the universe. The big bang model was marred with some inward problems related to it, which were removed by introducing an exponentially expanding phase in the early universe known as inflation. de Sitter presented a model of the universe devoid of ordinary matter, however with the cosmological constant term retained. The geometry of the model was proved to be accelerating [4]. The de Sitter universe corresponds to the specific case related to one of the very early solutions of Einstein’s Field Equations (EFE). The importance of the de Sitter model was not recognized until the introduction of inflation in the late 20th century, as the actual universe must be considered as a local set of perturbations in the geometry of de Sitter having validity at large. de Sitter geometry represents Euclidean space with a metric that depends on time. It was found that the inflation could be the de Sitter in general or quasi-de Sitter geometry which has an innate impact on the evolution of the geometry of FLRW spacetimes. It further bears its relation with the late-time accelerated expansion of the universe and to the dynamic geometry of the spacetime intrinsically cohered with it. The paradigm of inflation, as it is propounded, has a profound impact on the evolution of the universe as the geometry of spacetime. The de Sitter universe represents the inflationary phase of the universe with slightly broken time translational symmetry.
Alexander Friedmann predicted theoretically the universe to be dynamic, the one which can expand, contract, or even be born out of a singularity [5]. George Lemaitre, unaware of Friedmann’s work at that time, independently reached the same conclusion. In 1931, he also proposed a theory of the primeval atom which came to be known later on as the big bang theory by Fred Hoyle accidentally [6]. Edwin Hubble first proved the existence of other galaxies besides the Milky Way and afterward in 1929 discovered based on observational evidence that the universe is actually expanding [7]. This was actually discovering what Friedmann already had predicted theoretically in 1922. In the late 1940s, George Gamow (1904–1968) and his collaborators, Ralph Alpher (1921–2007) and Robert Herman (1914–1997), independently worked on Lemaître’s hypothesis and transformed it into a model of the early universe. They made supposition about the initial state of the universe comprising of a very hot, compressed mixture of nucleons and photons, thereby introducing the big bang model on the basis of comparatively strong evidences. They did not associate a particular name with the early state of the universe. Based on this model they were successful in calculating the amount of helium in the universe, but unfortunately, there was no authentic observational evidence through which their calculations could be compared [8].
The standard relativistic model of cosmology underpinning big bang theory could not explain the global structure of the universe and the origin of matter in it. The distribution of matter in it homogeneously on large scales and the spatial flatness also remained enigmatic. The big bang model just made an assumption about these but could not solve them. In the framework of effective field theory, the aspects of nonsingular cosmology were explored by Yong Cai et al. It is shown that the effective field theory assists in having the clarification about the origin of no-go-theorem and helps to resolve this theorem [9].
The inflationary era was proposed in the standard model of cosmology which propounds the big bang theory of the creation of the universe. Inflation solves the problems encountered in the big bang cosmology. Gliner, in 1965, hypothesized an era of exponential expansion for the universe earlier than any significant inflationary model surfaced [10]. It was found that the scalar fields are dynamic in nature, and in 1972 it was proposed that during phase transitions the energy density of the universe as a scalar field changes [11]. Andrei Linde, in 1974, realized that scalar fields can play an important role in describing the phases of the very early universe. He speculated that the energy density of a scalar field can play the role of vacuum energy dubbed as a cosmological constant [12].
In 1978, Englert, Brout, and Gunzig [13] forwarded a proposal of “fireball” hypothesis attempting to resolve the primordial singularity problem. They based their investigations on the entropy contained in the universe and approached the issue of early evolution of the universe by introducing particle production in it. They inquired deep down into it and on the basis of their hypothesis inferred that a universe undergoing a quantum mechanical effect would itself appear in a state of negative pressure and would be subject to a phase of exponential expansion. A work was mentioned by Linde in his review article [14] where he sought, in collaboration with Chibisov, to develop a cosmological model based upon the facts known in connection with the scalar fields. Considering the supercooled vacuum as a self-contained source for entropy, they tried to bring about the exponential expansion of the universe to be concerned with it. They, however discovered instantly that the universe becomes very inhomogeneous after the bubble wall collisions take place.
Slightly before Alan Guth’s original proposal of inflation surfaced, Alexei Starobinsky in 1980 proposed a model of inflation on the base of a conformal anomaly in quantum gravity. His proposal was presented in the framework of general relativity where slight modification of the equations of general relativity in matter sector was proposed and quantum corrections were employed to it in order to have a phase of the early universe. Starobinsky’s model can be considered as the first model of inflation which is of semi-realistic nature and evades from the graceful exit problem [15]. It was hardly concerned with the problems of homogeneity and isotropy which occur in the relativistic cosmological model of the big bang. His model, as he himself accentuated, can be considered the extreme opposite of chaos in Misner’s model. The model is found to agree with cosmological observations with slight deviations form recent measurements. Tensor perturbations that represent gravitational waves have also been predicted in Starobinsky’s model with a spectrum that is flat.
Alan Guth employed the dynamics of a scalar field and with a clear physical motivation presented an inflationary model [16] in 1981 on the base of supercooling theory during the cosmic phase transitions where the universe expands in a supercooled false vacuum state. A false vacuum is a metastable state containing a huge energy density without any field or particle so that when the universe expands from this heavy nothingness state its energy density does not change and empty space remains empty so that the inflation occurs in false vacuum [17]. The duration of inflationary phase in Guth’s original scenario is too short to resolve any problem, although was supposed to solve these problems and consequently the universe becomes very inhomogeneous which leads to the graceful exit problem [18,19]. The problem prevents the universe from evolving to later stages and is inherently existing in the originally proposed version of Guth.
The graceful exit problem was addressed independently by Linde, Steinhard, and Albrecht [20,21,22,23,24,25], where they introduced a phase of slow roll inflation at the end of the normal inflationary phase inclusively known as new inflation. The resolution of the problem was sought by constructing a new inflationary paradigm where the inflation can have its inception either in an unstable state at the top of the effective potential or in the state of false vacuum. In this scenario, the dynamics of the scalar field is such that it rolls gradually down to the lowest of its effective potential. It is of great importance to note that the shifting away of the scalar field from the false vacuum state to other later states has remarkable consequences. When the scalar field rolls slowly towards its lowest so-called slow roll inflation, the density perturbations are generated which seed the structure formation of the universe [26,27,28]. The production of density perturbations during the phase of slow roll inflation is inversely proportional to the motion of the scalar field [29,30]. The basic difference between the new inflationary scenario and that of the old one is that the advantageous portion of the inflation in the new scenario, which is responsible for the large scale homogeneity of the universe, does not take place in the false vacuum state, where the scalar field vanishes. This means that the new inflation could explain why our universe was so large only if it was very large initially and contained many particles from the very beginning.
The course of 20th century has presented many challenges to the standard cosmology. In the framework of the standard model, in addition to inflation, another breakthrough came forth in 1998 when the observation-based accelerated expansion of the universe was discovered [31,32,33]. Before this discovery, however, it was thought that in the perspective of all known forms of matter and energy that obey the strong energy condition ρ + 3 p > 0 , the expansion of the universe would slow down with the passage of time. This was a natural consequence of Friedmann equations that play a central role in the evolution of the universe. From the acceleration equation a ¨ a = 4 π G 3 ρ + 3 p , the universe must be undergoing deceleration characterized by deceleration parameter q 0 = a a ¨ a ˙ 2 ; however, astoundingly the value of q 0 < 0 was observationally determined, meaning that the expansion of the universe is accelerating rather to be decelerating. The discovery of accelerated expansion has won Noble prize in 2011. To explain the cause of accelerated expansion an exotic form of energy density was introduced hypothetically known usually as dark energy. The present budget of the universe from the observational data is contributed by dark energy 70 % , dark matter 25 % and 5 % ordinary baryon matter [34,35]. Dark energy is effective on the largest scales of intra galaxies and does not affect gravitationally bound systems. To explain the origin of dark energy, there is a large number of proposed models. Many independent observations lend support to the existence of dark energy such as CMB, SN Ia, BAO, etc. Today dark energy constitutes a very significant subject of relativistic cosmology with observational data by providing information about its basic nature, for reviews see in [36,37,38]. In this article, we study the standard model of cosmology by investigating the geometric structure of spacetime related with it in the framework of general relativity. Beginning with Euclidean space, we study spacetime in special and general theory of relativity. We discuss problems encountered in the standard big bang cosmology, and the inflationary solutions introduced into it, by proposing a phase of accelerated expansion in the early universe. A discussion of f R modified gravity is also presented with discussion of how inflation and dark energy can be described in its framework.
The layout of the paper is as follows. In Section 2, we discuss the structure of Euclidean space beginning with the axioms of Euclidean geometry and the significant role played by the Pythagoras theorem in its development. It has four subsections discussing space, time, and spacetime in relativity and pre-relativity physics. Section 3 begins with relativistic cosmology with a discussion on its underlying principles. The standard model of cosmology is discussed in Section 4 with nine subsections about its geometric structure. In Section 5, the derivation of Friedmann equations is carried out. Section 6 describes different aspects of embedding a geometrical object in a space of higher dimensions. It is has four subsections. Section 7 presents the very first relativistic model developed by Einstein himself. It has two subsections that discuss the instability of Einstein’s universe and de Sitter’s empty universe model respectively. In Section 8, a discussion on conformal FLRW line elements is presented in addition to the vacuum, radiation, and matter-dominated eras. It has 12 subsections covering related topics. Section 9 with its four subsections is devoted to the discussion of cosmological problems faced by the standard model. In Section 10, we embark on inflation and discuss its dynamics. Section 11 describes how the proposal of exponential expansion in the early universe solves the cosmological problems. Λ CDM and f R are discussed in this section. In the last Section, we provide a summary of the paper. Four indices are added in the end.

2. Euclidean Space

Euclidean geometry is established on a set of simple axioms and the definitions derived from these axioms. These axioms were first stated by Euclid in about 300 B.C. [39]. A space at the level of mathematical abstraction is the set of points where each point represents a specific position in it. When an abstract space is mapped onto a physical space, each point of it represents a physical location in it. Euclidean space is what entails on the base of axioms of Euclidean geometry. Geometrically, a space can be described by reducing it to a certain specification of the distance between each pair of its neighboring points. In order to reduce all of the geometry of a space to a certain specification of the distance between each pair of neighboring points we use the metric or line element which measures the space and describes its nature. A line element specifies a certain geometry and its form varies corresponding to different coordinate systems. Five basic postulate lie at the core of Euclidean space and are the basis of standard laws of geometry:
1.
Any two points can be joined by a straight line, i.e., the shortest distance between two points is a straight line.
2.
A straight line can be extended to any length.
3.
A circle can be drawn with a given a line segment as radius and one end as center of the circle.
4.
All right angles are congruent.
5.
Given a line and a point not on the line, it is possible to draw exactly one line through the given point parallel to the line, i.e., parallel lines remain a constant distance apart. Pythagoras theorem was known before Euclid and can also be derived from the five postulates and is used to find distance between any two points in Euclidean space. A mathematical space is an abstraction used to model the physical space of the universe. The Euclidean space consists of geometric points and has three dimensions. Now the Pythagorean theorem for a right triangle describes how to calculate the length of hypotenuse when the lengths of other two sides namely base and altitude are given. The length of hypotenuse gives distance between two points. Figure 1 below shows Pythagoras theorem:
d 2 = x 2 + y 2
Now, as the space can be described everywhere consisting of geometric points, we can define mutual relation for every infinitesimally close three points of the space forming a right triangle so that we can determine the element of distance between any two points with the help of Pythagoras theorem. Using rectangular Cartesian coordinate system we can express distance between two points in differential form as
d l 2 = d x 2 + d y 2
The distance-measure by Pythagoras theorem in Equation (2) will be known as metric or line element in two dimensions and defines Euclidean metric for two dimensional space. The distance measured between two points by the metric in Equation (2) does not change on rotating the coordinate system in which these two points are specified as Figure 2 manifests it:
x y = cos θ sin θ sin θ cos θ x y
or
x = x cos θ + y sin θ y = x sin θ + y cos θ
The distance between two points remains invariant which means that
d x 2 + d y 2 = d x 2 + d y 2
The Pythagorean theorem in three dimensions can be described as
d 2 = x 2 + y 2 + z 2
Three mutually perpendicular planes along three dimensions of the Cartesian coordinate system divide it in 3-planes as is shown in Figure 3:
Now, in reference to a coordinate system each point of this space will have three coordinates x , y , z if we approach its structure through Cartesian scheme, i.e., in Cartesian coordinates each point of it is represented by three coordinates which are the distances measured starting from the origin of the coordinate axes along the corresponding axes, i.e., x-axis, y-axis, and z-axis, respectively. These three axes stand for three dimensions of space. We find the distance between two points with Cartesian coordinate for three points separated infinitesimally
d l 2 = d x 2 + d y 2 + d z 2
which gives the metric of three-dimensional space The distance between two points with Cartesian coordinates x , y , z and p , q , r will be
d s 2 = x p 2 + y q 2 + z r 2
The infinitesimal distance between any two points x , y , z and x + d x , y + d y , z + d z can be had using the metric written above in Equation (7) in three dimensional Euclidean space.
d s 2 = x x + d x 2 + y y + d y 2 + z z + d z 2
d s 2 = d x 2 + d y 2 + d z 2 = d x 2 + d y 2 + d z 2
or in tensor form
d s 2 = δ μ ν d x μ d x ν
where δ μ ν is the Kronecker delta function representing a symmetric tensor of rank two and can be expressed as a 3 × 3 matrix form
δ μ ν = δ 11 δ 12 δ 13 δ 21 δ 22 δ 23 δ 31 δ 32 δ 33 = 1 0 0 0 1 0 0 0 1
diagonal of
δ μ ν = d i a g δ μ ν = + 1 , + 1 , + 1
and trace of
δ μ ν = μ = ν δ μ ν = 1 + 1 + 1 = 3
therefore, δ μ ν = d i a g + 1 , + 1 + 1 in Equation (13) defines an n-dimensional Euclidean space.
Now Equation (11) can be expanded using Einstein summation convention
d s 2 = δ 1 ν d x 1 d x ν + δ 2 ν d x 2 d x ν + δ 3 ν d x 3 d x ν
d s 2 = δ 11 d x 1 d x 1 + δ 12 d x 1 d x 2 + δ 13 d x 1 d x 3 + δ 21 d x 2 d x 1 + δ 22 d x 2 d x 2 + δ 23 d x 2 d x 3 + δ 31 d x 3 d x 1 + δ 32 d x 3 d x 2 + δ 33 d x 3 d x 3
d s 2 = 1 d x 1 d x 1 + 0 d x 1 d x 2 + 0 d x 1 d x 3 + 0 d x 2 d x 1 + 1 d x 2 d x 2 + 0 d x 2 d x 3 + 0 d x 3 d x 1 + 0 d x 3 d x 2 + 1 d x 3 d x 3
d s 2 = d x 1 d x 1 + 0 + 0 + 0 + d x 2 d x 2 + 0 + 0 + 0 + d x 3 d x 3
d s 2 = d x 1 d x 1 + d x 2 d x 2 + d x 3 d x 3
d s 2 = d x 1 2 + d x 2 2 + d x 3 2
d s 2 = d x 2 + d y 2 + d z 2
Equation (21) can also be written in the form
d s · d s = d x · d x + d y · d y + d z · d z
From Equation (22), we can see that the inner product in three dimensional Euclidean space can be perfectly described, that is why three dimensional Euclidean space is an example of a complete inner product space. An explanatory discussion of maximally symmetric 3 space can be consulted in Appendix B.

2.1. Newtonian Mechanics: The Structure of Space and Time

Space and time are absolute structures in classical physics and can be distinguished from one another in an independent way. Newton’s Mechanics is based specifically on three laws of motion, a law of gravitation and Galilean principle of relativity which are inherently related with the properties of space and time. Newtonian space is a three-dimensional extension around us which constitutes absolute space. Absolute space in Newton’s own words is described as “Absolute space, in its own nature, without relation to anything external remains always similar and immovable”, therefore space is rigid, motionless, and can be viewed as colossally empty three-dimensional cubic or cuboidal box where material objects reside and all physical phenomena take place. Newtonian space has the properties of Euclidean pace where infinitesimal distance between any two points is a straight line and if three points constitute a right angled triangle, then three sides are related by Pythagoras theorem which ascribes to it the properties of a flat space. Sum of angles in a triangle in such space is 180 . Newtonian space is homogeneous and isotropic which entails Newtonian Mechanics. Homogeneity implies translational invariance of the properties of space which means that it has similar properties at every point contained in it. The property of being homogeneous is called homogeneity that leads to the invariance of physical laws performed in two or more coordinate systems. Newton’s 3rd law, law of conservation of momentum and energy, etc. come out as a consequence of homogeneity of space. It is also an isotropic that implies rotational invariance of the properties of space. It means that it has similar properties in all directions and is therefore direction-independent. Thus, isotropy implies homogeneity but the converse is not true. The absolute time has been enunciated as follows “Absolute time, and mathematical time of itself and from its own nature flows equably without relation to anything external, and is otherwise called duration” such time exists independent of space and whatever dynamically happens in it and flows uniformly in one direction. An interval of time possesses always unchanging meaning for all times. This is presented figuratively in Figure 4 below:
According to Newtonian Mechanics, gravitation and relative motion do not affect the rate at which time flows. From Newton’s 2nd law F = m a , the isotropy of time can be viewed in case of a dynamic system that does not change from perpetrating transition from + t to t . This is because it does not incorporate the element of time explicitly which implies that past and future are indistinguishable but this is paradoxical because time is unidirectional and flows always from past to future. Two observers in two inertial frames of reference in relative motion and equipped with standard measuring clocks record the spacetime coordinates of an event written as t , x , y , z and t , x , y , z , respectively. According to Galilean principle of relativity, the coordinate transformations are
x = x v t y = y z = z t = t
We can calculate the addition of velocities according to these transformations by differentiating the spatial parts of Equation (23) with respect to time t, we have
d x d t = d x d t v d y d t = d y d t d z d t = d z d t
As t = t , we infer that d x d t = d x d t . Likewise, acceleration can also be differentiated once again from Equation (24), which gives
d 2 x d t 2 = d 2 x d t 2 d 2 y d t 2 = d 2 y d t 2 d 2 z d t 2 = d 2 z d t 2
We can observe from Equation (25) that the accelerations in both frames are same. The time-coordinate t of one inertial frame remains unaffected during transformation to another inertial frame of reference in classical physics and does not depend on spatial coordinates x, y, and z. The set of equations in Equation (23) are known as Galilean transformations. The motion along y and z spatial dimensions remains unaffected and the time coordinates in the two frames are equivalent which implies that time is absolute as Newton believed meaning that for all the inertial observers the time interval between any two events would be invariant. We notice that the two events having coordinates t , x , y , z and t , x , y , z , respectively, with differential of the distance as Euclidean spatial interval described in Equation (21) as d s 2 = d x 2 + d y 2 + d z 2 and the time interval Δ t = t t both remain separately invariant under the Galilean transformations in Equation (23). This fact makes us consider the nature of space and time as absolute entities in Newtonian Mechanics. We identify the quantity d s 2 as square of the distance between points of three-dimensional Euclidean space and invariance of this differential of distance alludes to the fact that it is geometrical structural property of the space itself in its own right. This describes the geometry of space and time according to Newton’s views.

2.2. Special Theory of Relativity: The Structure of Spacetime

Special relativity is a theory of the structure of spacetime and in this way constitutes a geometric theory [40]. The fields and particles grow over this spacetime structure and relativistic mechanics is developed according to this structure which corresponds to the postulates of special relativity. According to the Lorentz transformations implied by it, space and time are not distinguishable quantities but constitute innately a single continuum to be known as spacetime. One of the Einstein’s 1905 papers brought forward this theory founded upon two postulates [41].
1.
The principle of special covariance.
2.
The principle of invariance of the velocity of light (c). As the laws of physics remain form-invariant, i.e., covariant according to a privileged class of observers known as inertial frames. This is also called principle of relativity. These two principles overthrew the pre-relativity notions of absolute space and absolute time proposing instead relative concepts. In classical physics as we saw earlier the coordinates of two observers are related by Galilean transformations, whereas according to the special relativity, the coordinates in two frames are related using Lorentz transformations.
x = x v t 1 v 2 c 2 y = y z = z t = t v x c 2 1 v 2 c 2
Lorentz transformations contain all the geometric information about space and time, and describe the structure of spacetime. Further, we can see that space and time coordinates are absolute according to the Galilean transformations for two inertial observers which move relative to each other and are connected through space and time coordinates. Time coordinate has the same magnitude in pre-relativity physics; however, according to the special relativity, which obeys Lorentz transformations, the time coordinate in one coordinate system is connected to the time coordinate of the second coordinate system through both time and space coordinates, which alludes to the fact that space and time coordinates are now to be dealt on equal footings. It is obvious from the Lorentz transformations that the time coordinates are not equivalent in two frames, i.e., t t rather t is innately cohered with both of the coordinates of time and space t and x respectively. It means that time t of one coordinate frame converts partially in space and partially in time coordinates. Therefore, t does not remain independent but has partially coalesced with space coordinates losing its absolute nature and the principle of relativity forbade us to locate a preferred frame of reference ensuing that absolute notion of time disappears logically. This fact was first perceived by Minkowski when he was recasting the special relativity in the language of geometry. He has presented a very profound and significant geometrical structure underlying special relativity. While delivering a lecture at the meeting of the Göttingen Mathematical Society on 5 November 1907, he introduced the concept of spacetime continuum whereby he asserted that independent space and time have to doom away into mere shadows and only a union of the two can preserve an independent reality. Minkowski viewed that the principle of special relativity can be described by the metric d t 2 + d x 2 + d y 2 + d z 2 on the four-dimensional space R 4 which familiarized the concept of spacetime continuum and paved the way for the formulation of general relativity. A Minkowski metric g on the linear space R 4 is a symmetric non-degenerate bilinear form with signature , + , + , + . It means that there exists a basis e 0 , e 1 , e 2 , e 3 such that g e μ , e ν = g μ ν where μ , ν 0 , 1 , 2 , 3 and g μ ν is expressed in the form
g μ ν = g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
so that an we have orthonormal basis and can construct a system of coordinates of R 4 as x 0 , x 1 , x 2 , x 3 such that at each point we can have e 0 = t and x j where j = 1 , 2 , 3 . Now, with respect to this coordinate system, we can write the metric tensor 0 , 2 in the form g = g μ ν d x μ d x ν = d t 2 + 1 3 d x j or d s 2 = d t 2 + d x 2 + d y 2 + d z 2 The negative sign with one time component term in the metric indicates that it is not Euclidean space but represents a pseudo-Euclidean known as Minkowski space and also guarantees that the speed of light is same in all inertial frames. An expanding Minkowskian spacetime can be described in the form as written below which represents the simplest of all dynamic spacetimes d s 2 = d t 2 + a 2 ( t ) d x 2 + d y 2 + d z 2 . It was thought convenient on the dimensional grounds to introduce the coordinates in the form x 0 , x 1 , x 2 , x 3 = c t , x , y , z . Pythagoras theorem applied in Euclidean space R 3 of three spatial dimensions gives the distance of two points as an invariant as we observed in previous section.
d s 2 = + d x 2 + d y 2 + d z 2
here d s the length element is a scalar quantity which means that in certain frame of references all the observers will agree upon the length of the measured object. In 1905, Einstein speculated that the measurement of the spacetime interval
d s 2 = d x 2 d y 2 d z 2 + c d t 2 = η μ ν d x μ d x ν
where
η μ ν = g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
would not result in identical either in space or in time [42] for the observers in relative uniform motion. However, Minkowski noted that the four-dimensional entity in Equation (29) would remain invariant for all such observers. The basic significant idea which Minkowski took notice of was that the spacetime interval remains invariant for all the observers in uniform relative motion meaning that it is also a scalar upon which they all will agree. The metric of Minkowski space which is homogenous and isotropic is given by
g μ ν = η μ ν = d i a g 1 , + 1 , + 1 , + 1
thus the geometry of spacetime is flat in special relativity. It is notable here that it is spacetime that is flat, however in classical mechanics, it is space rather than spacetime. If the Minkowskian geometry of spacetime is required to be expanding, it can be made so. However, in the framework of special relativity, it does not need to expand. Figure 5 gives the structure of Minkowskian spacetime as a null cone structure. In the Figure 6, it is shown how the dimension time is converted in as space.
g μ ν = η μ ν = d i a g 1 , a 2 ( t ) , a 2 ( t ) , a 2 ( t )

2.3. General Theory of Relativity: The Structure of Spacetime

The essence of general relativity is that geometry is gravity which comes from Equivalence principle. It models gravity into the dynamic structure of spacetime. In general relativity, the structure of spacetime is described by a fundamental quantity called the spacetime metric g μ ν or line element which gives the nature of the geometry of spacetime by finding the distance between two infinitesimally neighboring points in it. The geometrical structure of spacetime is incarnated [43] in two basic principles.
1.
Principle of general covariance.
2.
The spacetime continuum has, at each of its points, a quadratic structure of coordinate differentials d s 2 = g = g μ ν d x μ d x ν known as “square of the interval” between the two points under consideration.
We consider a four-dimensional continuum every point of which is distinct from the other with four coordinates-a quadruplet x 1 , x 2 , x 3 , x 4 assigned consecutively to each of them
x 1 = x 1 x 1 , x 2 , x 3 , x 4 x 2 = x 2 x 1 , x 2 , x 3 , x 4 x 3 = x 3 x 1 , x 2 , x 3 , x 4 x 4 = x 4 x 1 , x 2 , x 3 , x 4
It is denoted by g μ ν . In matrix form with components, it is written as
g μ ν = g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 a 31 g 32 g 33
The properties of spacetime that are intrinsically related to it, are completely determined by the spacetime metric. An example of the local curved spacetime around the Sun in two dimensions is displayed in Figure 7:
A detailed discussion of space, time and spacetime is presented in Appendix A.

2.4. The Basics of General Relativity

It would be convenient to have a retrospective look into the basics of general relativity whose role has been very fundamental to the modern cosmology. We briefly review the structure of the theory specifically in connection with the geometrical structure of spacetime in it. General relativity in its core describes that gravity is the geometry of four-dimensional spacetime manifested through its curvature. It is a theory of spacetime and gravitation that are the very basic components of the universe. Einstein’s journey towards general relativity in order to introduce gravity in his previous theory sought the fascinating geometry of the structure of spacetime, such that gravity as a field force disappeared and was assimilated in the very geometric structure of spacetime. In constructing the framework of new theory, Einstein was influenced and governed by Mach’s principle, which states that it is a priori existence and distribution of matter which determines the geometry of spacetime, and in the absence of it, there shall be no geometric structure of a spacetime in the universe. Therefore, there will be no inertial properties in an, otherwise, empty universe. In general, relativity gravitation and inertia are essentially indistinguishable. The metric tensor g μ ν describes the effect of both combinedly, and it is arbitrary to ask which one contributes its effect more and which less, therefore to call it with a single name is suitable either inertia or gravitation [4]. In general relativity gravitation, inertia and the geometry of spacetime are coalesced into a single entity represented by a symmetric tensor of second rank g μ ν which owes its existence due to presence and distribution of matter which is represented by an other symmetric tensor T μ ν known as energy-momentum tensor. The metric tensor g μ ν is the fundamental object of study in general relativity and takes into consideration all the causal and geometrical structure of spacetime. General relativity underlies five fundamental principles connotated by it implicitly or explicitly manner:
  • Mach’s principle
  • principle of equivalence
  • principle of covariance
  • principle of minimal gravitational coupling
  • correspondence principle
In the light of the principle of general covariance, the theory requires that the laws of physics might be formulated in a coordinate-independent style. The coordinate independence requires the replacement of partial derivatives by covariant derivatives which introduces connection coefficients Γ μ ν λ as the 2nd kind of Christoffel symbols. All the geometric structure of spacetime is based on the existence of these connection coefficients. The field equations of general relativity read as G μ ν = 8 π T μ ν , where G μ ν = R μ ν 1 2 g μ ν R is the Einstein tensor and is expressed in terms of Ricci tensor, metric tensor, and Ricci scalar, and T μ ν is energy momentum tensor. The spacetime continuum of general relativity is postulated as a 4-dimensional Lorentzian manifold (M, g), where M denotes the Manifold and g is metric defined over it. The geometry of a spacetime is encoded in its metric which has a geodesic structure, though complex and frequently solved numerically for a specific bunch of geodesics. These geodesics specify the physical properties of the geometry of spacetime which are interpreted by drawing graphically in a certain spacelike volume. Gravity is the geometry of spacetime itself which is described through its dynamic structure in the framework of general relativity. The interaction between spacetime and the content it contains which mutually form and the universe is the pith and marrow of general relativity. Matter tells spacetime how to curve and spacetime tells the matter how to move. General relativity thus transforms gravitation from being a force to being it a property of spacetime, so that gravity does remain a force but curvature of the geometric structure of spacetime. Einstein worked out a relation between matter–energy content of the universe and its gravitating effects in the form of geometry of spacetime. He employed the language of tensors to describe it. The invariant interval between two events separated infinitesimally with coordinates t , x , y , z and t + d t , x + d x , y + d y , z + d z has been defined according to special relativity
d s 2 = η μ ν d x μ d x ν
Which defines a Lorentz invariant Minkowski flat spacetime whose geometry of spacetime is encoded in η μ ν . Under the change of coordinates d s 2 remains invariant and is spacelike for d s 2 > 0 , timelike for d s 2 < 0 and light-like for d s 2 = 0 . Photon path is described by d s = 0 and baryonic matter follows a path between two events for which
d s = 0
i.e., it generates stationary values and conforms to the shortest distance between two points to be straight line which means that there are no external forces to set their path deviated. General relativity was based on five principles incorporated in it explicitly or implicitly, namely, equivalence principle, relativity principle, Mach’s principle, and Correspondence principle. Tensors are geometric objects defined on a manifold M, which remain invariant under the change of coordinates. It is composed of a set of quantities which are called its components, therefore a it is the generalization of a vector which means that it has more than three components. They represent mathematical entities which conform to certain laws of transformations. The properties of components of a tensor do not depend on a coordinate system which is used to describe the tensors. Transformation laws of a tensor relate its components in two different coordinate systems. The mathematical representation of a tensor is displayed through considering usually a bold face alphabetical letter like A, B, T, P, etc. with an index or a set of indices in the form of superscripts or subscripts or both in mixed form. These superscripts and subscripts in case of a tensor are called contravariant and covariant indices. Contravariant indices of a tensor are used to give the meaning of a contravariant components of it like A μ , A μ ν , A μ ν ξ . . . Covariant indices of a tensor are used to signify the meaning of a contravariant components of it like A μ , A μ ν , A μ ν ξ . . . The indices of both types namely contravariant and covariant are used to specify the components of a mixed tensor like A μ ν , A μ ξ ν , A μ ν σ , A μ ξ υ ν σ . A mixed tensor is a tensor which has contravariant as well as covariant components. The number of indices appearing in the symbol representing certain type of a tensor is known as its rank. The appearing indices in the symbol representing a tensor can be contravariant or covariant or both type of indices in it. The order of a tensor is the same thing as rank, only the name differs. The number of components of a tensor is related with its rank or order and the dimensions of the space in which the is being described. In an n-dimensional space, a tensor of rank, say, k will have number of components equal to number of components of a tensor in n-dimensional space is equivalent to n k = n u m b e r o f d i m e n s i o n s o f s p a c e r a n k . However, the spacetime of general relativity is pseudo-Riemannian having four dimensions, three spatial and one temporal. Coordinate patches are necessarily considered to map whole of the spacetime. Each point-event of a coordinate patch in the four-dimensional pseudo-Riemannian spacetime is labeled by a general coordinate system, which conventionally runs over 0, 1, 2, and 3, where 0 stands for time and the rest for space coordinates. An inertial or otherwise frame of reference characterized by a coordinate system can be attached to every point event of the spacetime and coordinate transformations between any two coordinate systems can be found. These can be written as
A μ = x ν x μ A ν B μ = x μ x ν B ν A ν μ = x μ x ζ x σ x ν A σ ζ
while switching to Riemannian geometry for non-Euclidean spaces ordinary partial differentiation is generalized to covariant differentiation and is defined using a semi-colon; as
B ν ; μ = , μ B ν Γ ν μ σ B σ B ν ; μ = , μ B ν + Γ μ σ ν B σ
where comma, denotes an ordinary partial differentiation with respect to the corresponding variable and; signifies covariant differentiation. In the covariant differentiation, indices can also be raised or lowered with metric tensor, however the covariant differentiation of it vanishes, i.e., g μ ν ; α = 0 . The interval between infinitesimally separated events x μ and x μ + d x μ is given by
d s 2 = g μ ν d x μ d x ν
The corresponding contravariant tensor of g μ ν is given by g μ ν and they result in Kronecker delta. Moreover, indices can be lowered or raised using the metric tensor in either form as
g μ ν g μ ζ = δ ν ζ g μ ν B ν = B μ g μ ν B ν = B μ
In general relativity, all the geometry of curved spacetime is contained in the second-rank symmetric tensor g μ ν known as fundamental or metric tensor and is the function of four coordinates g μ ν = g μ ν x 0 , x 1 , x 2 , x 3 and g μ ν encodes all the information about gravitational field induced by presence of matter. It governs the other matter as a response mimicking the role of gravitational potential similar to that of Newtonian gravity so that the paths remain no more straight, and the action in Equation (36) determines the path of a free particle known as geodesic
d 2 x μ d s 2 + Γ ν ζ μ d x ν d s d x ζ d s = 0
where
Γ ν ζ μ = g μ λ Γ ν ζ λ = 1 2 g μ λ g ν λ x ζ + g ζ λ x ν + g ν ζ x λ
are the Christoffel symbols which through the geodesic equation specify the world lines of free particles. The “acceleration due to gravity” in Newtonian gravitation law is described by these symbols in Einstein’s picture of gravity as the geometric properties of spacetime encoding the similar information. Locally these symbols vanish in the inertial frame of reference in free fall and under coordinate transformation from x μ and x μ do not constitute components of a tensor and therefore do not represent a tensor.
Γ μ ν σ = x σ x λ x ζ x μ x ρ x ν Γ ζ ρ λ + 2 x ζ x μ x ν x σ x ζ
The Riemann tensor is defined as
R μ ν λ σ = x ν Γ μ λ σ x λ Γ μ ν σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ
It has symmetry properties and satisfies the following Bianchi identity:
R μ ν λ ; ζ σ + R μ ζ ν ; λ σ + R μ λ ζ ; ν σ = 0
The Ricci tensor is obtained from Riemann tensor contracting
R μ ν = g λ σ R μ ν λ σ = x ν Γ μ λ σ x λ Γ μ ν σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ
Another expression of Ricci tensor is written in the form given below when determinant of the metric tensor g μ ν is envisaged as a matrix and denoted by g
R μ ν = Γ μ ν , λ λ ln g , μ ν + ln g , λ Γ μ ν λ Γ π μ λ Γ λ ν π
The Ricci scalar or scalar curvature is described as
R = g μ ν R μ ν
Contraction of the Bianchi Identity in Equation (45) gives
R μ ν 1 2 g μ ν R
which is the Einstein tensor. Now we can write basic equations of general relativity
R μ ν 1 2 g μ ν R = 8 π T μ ν
or
G μ ν = 8 π T μ ν
G μ ν T μ ν
These are written with cosmological constant also. From Equation (52)
G μ ν + Λ g μ ν = 8 π T μ ν
Energy-momentum tensor T μ ν is the source term for the metric tensor g μ ν which for a most general matter-energy fluid that is consistent with the assumption of homogeneity and isotropy represents a perfect fluid and has the form
T μ ν = ρ + p u μ u ν p g μ ν
where u μ = 1 , 0 , 0 , 0 is the four velocity in a comoving frame of reference and
T μ ν = ρ 0 0 0 0 p 0 0 0 0 p 0 0 0 0 p

3. Relativistic Cosmology

Relativistic cosmology was founded on three fundamental principles
  • Cosmological principle;
  • Weyl’s principle;
  • General relativity.
These are illustrated in the following subsections.

3.1. Cosmological Principle

The cosmological principle states that on sufficiently large scale, the universe is homogenous and isotropic at any time. Therefore, it is the same for all observers and has similar properties on larger scales. The principle is the generalization of Copernican principle and almost all the standard cosmological models of the spacetime underpin it. It has two forms:
(1)
Cosmological principle with respect to spatial invariance
(2)
Cosmological principle with respect to temporal invariance
In special invariance, we suppose the invariance of space with respect to translational and rotational properties known as homogeneity and isotropy, respectively, and the principle may be regarded as cosmological principle. Under both the invariant properties the space remains isomorphic. A perfect cosmological principle incorporates temporal homogeneity and isotropy which was employed by the steady state theory of the eternal universe and was not supported by the observation and was disfavored. For a local observer the principle might not be satisfied as the Earth and the solar system are not homogeneous and isotropic since the matter clumps together to form objects like planets, stars, galaxies with voids of vacuum-like in between them but on the larger scales of about MP > 1000 Pc the universe obeys the cosmological principle. The uniformity of CMBR in all directions (homogeneity and isotropy) provides the confirmatory proof of the cosmological principle. It is the generalization of Copernican Principle which incorporates homogeneity and isotropy. Homogeneity means location independence, i.e., all places in the universe at galactic scales are indistinguishable. Isotropy gives direction independence, i.e., in whatever direction we look in the universe it appears same. Certainly Isotropy connotes homogeneity, but this not true vice versa. To better understand its geometric properties, we begin with 1-dimensional spaces and revise to the four-dimensional spaces and then observe how the four-dimensional spacetime geometrical properties can be understood in this perspective. It is necessary to understand what we mean by embedding of a geometric object in an n-dimensional space because of the reason FLRW metric incorporates example of embedding three dimensional spaces in four dimensional spacetime. Figure 8 delineates homogeneity and isotropy properties of space:

3.2. Weyl’s Principle

Weyl’s principle helps us consider the universal stuff as consisting of a fluid, the particles of which are constituted by galaxies. Therefore, what we name “the universe” is just cosmic fluid. In the cosmological spacetime, the world lines of the fundamental observers form a smooth bundle of time-like geodesics which would never meet except in the past singularity from where the universe emerged or at the future singularity if it would happen. The fundamental observers are those who comove with the cosmic fluid. The world lines of galaxies as fluid particles are always and everywhere orthogonal to the family of spatial hypersurfaces. The postulate was presented by Hermann Weyl (1885–1955) in 1923 which is essentially about the nature of matter in the universe [44]. He regarded the material content of the universe in the form of fluid whose constituent particles make a substratum in the cosmic fluid.
It means that in the substratum of spacetime it allows us to consider the structure of the universe as fluid. The Weyl principle introduces further symmetry in the structure of spacetime described by the metric tensor by considering the galaxies as test particles and postulates that the geodesics on which these galaxies move do not intersect. It states that the world lines of galaxies considered as “test particles” form a 3-bundle of non-intersecting geodesics orthogonal to a series of spacelike hypersurfaces. A simple illustration of the Weyl’s principle is given in the Figure 9:

3.3. General Relativity

General relativity provides the best existing theory of gravitation on cosmological scales and models it structured into the geometric structure of spacetime. In Section 3, we discussed its basic ingredients.

4. The Standard Model of Cosmology

The standard model in cosmology has been established on the most general homogeneous and isotropic spacetime. The standard model that propounds the hot big bang model of the universe is known as Friedmann–Lemaitre–Robertson–Walker (FLRW) line element which reads as in the Cartesian coordinates
d s 2 = d t 2 + d x 2 + d y 2 + d z 2
and in the spherical coordinates, we have
d s 2 = g = g μ ν d x μ d x ν = d t 2 + a ( t ) 2 1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
Or equivalently
d s 2 = g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 = 1 0 0 0 0 a 2 ( t ) 1 k r 2 0 0 0 0 a 2 ( t ) r 2 0 0 0 0 a 2 ( t ) r 2 sin 2 θ
The predictions for the quantitative behavior of the expanding universe is enunciated suitably by the metric tensor and the scale factor as a function of time, i.e., a ( t ) describes the scale of coordinate grid interrelating the coordinate distance with physical distance, i.e., in a smooth and homogeneously expanding universe.

4.1. Geometric Properties of the FLRW Line Element

From the line element in Equation (57)
d s 2 = d t 2 + a ( t ) 2 1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
As time flows only in one direction and the space obeys cosmological principle, therefore we are allowed to separate the metric in temporal and spatial parts. To understand the four dimensional spacetime geometry of FLRW universe we begin with the geometry of spatial part of the line element that is
a ( t ) 2 1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
This is the spatial part of the metric in Equation (59) and is characterized by the scale factor a ( t ) , which is the function of time and 2nd curvature of the space k. These are obviously determined by the self-gravitating properties of the matter–energy content in the universe. The spatial part of the metric incorporates cosmological principle implying homogeneity and isotropy which provides the kinematics for the geometry of spacetime while we will observe afterwards that Einstein equations provide the dynamics into it through the scale factor a t .

4.2. Comoving Coordinates and Peculiar Velocities

The coordinates r , θ , ϕ form the cosmological rest frame and are known as comoving coordinates. They can be considered constant because the particles remain at rest in these coordinates. Peculiar velocity is the motion of the particles with respect to comoving coordinates. Peculiar velocities of the galaxies and supernovae are ignored in cosmology in the expanding spacetime. As p a 1 a t , therefore momentum in expanding spacetime is red-shifted and freely moving particles come to rest in comoving coordinates. Physical distance between two points is calculated as thee scale factor a t times the coordinate distance. The expression without scale factor inside the bracket is the pure kinematical statement of the geometry of spacetime
1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
and represents the line element of the three-dimensional space with hidden symmetry of being homogeneous and isotropic. It represents three geometries for three values of k.

4.3. The Geometry of Spherical World

For k = + 1 , the hypersurface is
1 1 r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
and represents a three dimensional sphere embedded in a four dimensional Euclidean space. This space is finite and closed.

4.4. The Geometry of Hyperbolic World

For k = 1 , the hypersurface is
1 1 + r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
and represents a three-dimensional hypersphere or hyperbola embedded in a four-dimensional pseud-Euclidean space. This space is infinite and open.

4.5. The Geometry of Euclidean World

For k = 0 , the hypersurface is
d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
and represents a three-dimensional Euclidean flat space. This space is also infinite and open. Now to determine Friedmann equations, we write first the components of the metric tensor, since the metric is diagonal due to homogeneity and isotropy therefore we have these diagonal components
g 00 = g t t = 1 g 11 = g r r = a 2 ( t ) 1 k r 2 g 22 = g θ θ = a 2 ( t ) r 2 g 33 = g ϕ ϕ = a 2 ( t ) r 2 sin 2 θ
Now, we turn to solve the FLRW metric and begins with finding Christoffel symbols of 2nd kind or the affine connections which are given by
Γ μ ν σ = g σ λ Γ μ ν λ = 1 2 g σ λ g μ λ x ν + g ν λ x μ + g μ ν x λ
In four dimensions these will have 4 3 = 64 components. The four generalized cases emerge in four dimensions for μ , ν , λ , and σ .
Case I: μ = ν = λ
Γ μ μ μ = 1 2 x μ log g μ μ
In four dimensions
Γ 00 0 Γ 11 1 Γ 22 2 Γ 33 3
will emerge.
Case II: σ = μ , μ ν
Γ μ λ μ = 1 2 x λ log g μ μ
In four dimensions the following twelve cases
Γ 01 0 Γ 02 0 Γ 03 0 Γ 10 1 Γ 12 1 Γ 13 1 Γ 20 2 Γ 21 2 Γ 23 2 Γ 30 3 Γ 31 3 Γ 32 3
will emerge.
Case III: σ = μ , μ = λ
Γ λ λ μ = 1 2 g μ μ g λ λ x μ
In four dimensions the following twelve cases
Γ 11 0 Γ 22 0 Γ 33 0 Γ 00 1 Γ 22 1 Γ 33 1 Γ 00 2 Γ 11 2 Γ 33 2 Γ 00 3 Γ 11 3 Γ 22 3
will emerge.
Case IV: σ μ ν
Γ μ ν σ = g σ λ Γ μ ν λ = 1 2 g σ λ g μ λ x ν + g ν λ x μ + g μ ν x λ = 0
In four dimensions the following twenty four cases:
Γ 12 0 Γ 21 0 Γ 13 0 Γ 31 0 Γ 23 0 Γ 32 0 Γ 02 1 Γ 20 1 Γ 03 1 Γ 30 1 Γ 23 1 Γ 32 1 Γ 01 2 Γ 10 2 Γ 03 2 Γ 30 2 Γ 13 2 Γ 31 2 Γ 01 3 Γ 10 3 Γ 12 3 Γ 21 3 Γ 20 3 Γ 02 3
emerge and vanish.

4.6. Non-Vanishing Christoffel Symbols

We determine the following non-vanishing Christoffel symbols of the 2nd kind with the help of formula given in Equation (66) for the metric in Equation (57), which is the metric of universe in standard cosmology.
Γ 11 1 = k r 1 k r 2 Γ 10 1 = Γ 01 1 = Γ 20 2 = Γ 02 2 = Γ 30 3 = Γ 03 3 = a ˙ ( t ) a ( t ) Γ 21 2 = Γ 12 2 = Γ 31 3 = Γ 13 3 = 1 r Γ 32 3 = Γ 23 3 = sin θ cos θ Γ 11 0 = a ( t ) a ˙ ( t ) 1 k r 2 Γ 22 0 = a ( t ) a ˙ ( t ) r 2 Γ 33 0 = a ( t ) a ˙ ( t ) r 2 sin 2 θ Γ 22 1 = r 1 k r 2 Γ 33 1 = r sin 2 θ 1 k r 2 Γ 33 2 = sin θ cos θ

4.7. Riemann Curvature Tensor

The Riemann curvature tensor R μ ν λ σ has 4 4 = 256 components in four dimensions from which only twenty components can possibly be non-vanishing. The Riemann tensor is given by
R μ ν λ σ = x ν Γ μ λ σ x λ Γ μ ν σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ
The possibly non-vanishing twenty components are given by
R 110 0 R 220 0 R 330 0 R 221 0 R 331 0 R 001 1 R 221 1 R 331 1 R 332 1 R 002 2 R 112 2 R 332 2 R 021 2 R 003 3 R 113 3 R 223 3 R 031 3
The non-vanishing components are
R 110 0 = a ( t ) a ¨ ( t ) 1 k r 2 R 101 0 = a ( t ) a ¨ ( t ) 1 k r 2 R 010 1 = a ¨ ( t ) a ( t ) R 001 1 = a ¨ ( t ) a ( t ) R 220 0 = a ( t ) a ¨ ( t ) r 2 R 202 0 = a ( t ) a ¨ ( t ) r 2 R 020 2 = R 030 3 = a ¨ ( t ) a ( t ) R 002 2 = R 003 3 = a ¨ ( t ) a ( t ) R 330 0 = a ( t ) a ¨ ( t ) r 2 sin 2 θ R 303 0 = a ( t ) a ¨ ( t ) r 2 sin 2 θ R 221 1 = R 223 3 = r 2 k + a ˙ 2 ( t ) R 212 1 = R 232 3 = r 2 k + a ˙ 2 ( t ) R 121 2 = R 131 3 = k + a ˙ 2 ( t ) 1 k r 2 R 112 2 = R 113 3 = k + a ˙ 2 ( t ) 1 k r 2 R 331 1 = R 332 2 = r 2 sin 2 θ k + a ˙ 2 ( t ) R 313 1 = R 323 2 = r 2 sin 2 θ k + a ˙ 2 ( t )

4.8. Ricci Curvature Tensor and Ricci Scalar

Ricci tensor R μ ν is obtained by contracting Riemann tensor R μ ν λ σ . We contract it by placing λ = σ , so that R μ ν λ σ = R μ ν λ λ = R μ ν In four dimensions it has 4 2 = 16 components:
R 00 R 11 R 22 R 33 R 01 R 10 R 02 R 20 R 03 R 30 R 12 R 21 R 31 R 13 R 23 R 32
The non-vanishing components are
R 00 = 3 a ¨ a R 11 = a ( t ) a ¨ ( t ) + 2 k + 2 a ˙ 2 1 k r 2 R 22 = r 2 a ( t ) a ¨ ( t ) + 2 k + 2 a ˙ 2 R 33 = r 2 sin 2 θ a ( t ) a ¨ ( t ) + 2 k + 2 a ˙ 2
Ricci scalar R is obtained by contracting Ricci tensor
R = g μ ν R μ ν
Using double sums and simplifying in four dimensions, we have
R = g 00 R 00 + g 11 R 11 + g 22 R 22 + g 33 R 33
R = 6 a ¨ ( t ) a ( t ) + a ˙ ( t ) a ( t ) 2 + k a 2 ( t )

4.9. Einstein Tensor G μ ν

Einstein tensor is defined in terms of Ricci tensor R μ ν , Ricci scalar R, and the metric tensor g μ ν . It is expressed as
G μ ν = R μ ν 1 2 g μ ν R
In four dimensions it has 4 2 = 16 , components. These are
G 00 G 11 G 22 G 33 G 01 G 10 G 02 G 20 G 03 G 30 G 12 G 21 G 31 G 13 G 23 G 32
The non-vanishing components are
G 00 = 3 a ˙ a 2 + k a 2 G 11 = g 11 2 a ¨ a + a ˙ a 2 + k a 2 G 22 = g 22 2 a ¨ a + a ˙ a 2 + k a 2 G 33 = g 33 2 a ¨ a + a ˙ a 2 + k a 2
In Equation (86), the spatial components of Einstein tensor can be written in a single equation of tensorial nature.
G μ ν = g μ ν 2 a ¨ a + a ˙ a 2 + k a 2
where μ = ν = 1 , 2 , 3 , and mixed Einstein tensor can be found by g ζ ν G μ ν = G μ ζ
G 0 0 = 3 a ˙ a 2 + k a 2 G 1 1 = g 11 2 a ¨ a + a ˙ a 2 + k a 2 G 2 2 = g 22 2 a ¨ a + a ˙ a 2 + k a 2 G 3 3 = g 33 2 a ¨ a + a ˙ a 2 + k a 2
Now, we calculate the energy-momentum tensor of a perfect fluid in mixed form. Cosmological principle and Weyl’s postulate imply the material content of the universe to be regarded as perfect fluid [1,2,3].
g ζ ν T μ ν = T μ ζ = T 00 T 01 T 02 T 03 T 10 T 11 T 12 T 13 T 20 T 21 T 22 T 23 T 30 T 31 T 32 T 33 = ρ 0 0 0 0 p 0 0 0 0 p 0 0 0 0 p
The non-vanishing components of energy-momentum tensor are
T 0 0 = ρ T 1 1 = p T 2 2 = p T 3 3 = p
Putting the values of Einstein tensor G μ ν and energy-momentum tensor T μ ν from Equations (88) and (89), respectively, in Einstein field equations
G μ ζ = 8 π T μ ζ
3 a ˙ a 2 + k a 2 = 8 π G ρ g 11 2 a ¨ a + a ˙ a 2 + k a 2 = 8 π G p g 22 2 a ¨ a + a ˙ a 2 + k a 2 = 8 π G p g 33 2 a ¨ a + a ˙ a 2 + k a 2 = 8 π G p

5. Derivation of Friedmann’s Equations

Now, using the Einstein field equations, we set to derive the Friedmann’s Equations that describe the evolution of the universe by relating the large-scale geometrical characteristics of spacetime to the large-scale distribution of matter–energy and momentum. From Equation (92), we can write
3 a ˙ a 2 + k a 2 = 8 π G ρ
2 a ¨ a + a ˙ a 2 + k a 2 = 8 π G p
For other two components listed in Equation (92) the 2nd and 3rd components repeat, therefore we will write only one time from the three components. From Equations (93) and (94) we derive the Friedmann’s Equations and an equation for the conservation of matter. Substituting Equation (93) in Equation (94) and performing simplification we get
a ¨ a = 4 π G 3 ρ + 3 p
and from Equation (93) which is the time-time component of the Einstein Equations.
a ˙ a 2 + k a 2 = 8 π G 3 ρ
for
a ˙ a = H
which is Hubble parameter and gives expansion rate. The above Equation (96) can be written as
H 2 + k a 2 = 8 π G 3 ρ
differentiating Equation (97) with respect to time ‘t
t H = t a ˙ a
we obtain
H ˙ = a ¨ a H 2
which gives
H ˙ + H 2 = a ¨ a
Therefore, that Equation (95) takes the form in terms of Hubble parameter.
H ˙ + H 2 = 4 π G 3 ρ + 3 p
we can also find
H ˙ = 4 π G 3 ρ + 3 p H 2
From Equation (98) H 2 = 8 π G 3 ρ with k = 0 , for flat universe substituting it in Equation (103) above
H ˙ = 4 π G 3 ρ + 3 p 8 π G 3 ρ
which results in
t H = 4 π G ρ + p
Now, differentiating Equation (93) with respect to time after shifting the factor 3 on the right side, we have
a ˙ a 2 a ¨ a 2 a ˙ a 2 2 k a 2 = 8 π G 3 ρ ˙
subtracting now Equation (93) from Equation (94), we obtain
2 a ¨ a 2 a ˙ a 2 2 k a 2 = 8 π G ρ + p
substituting Equation (107) in Equation (106), after simplification we have
ρ ˙ + 3 a ˙ a ρ + p = 0
Cosmological principle compels us to consider a fluid in which inhomogeneities will be considered smoothed out and evolution of the universe shall be considered in the form of perfect fluid characterized by energy density ρ and isotropic pressure p. Further we consider that the pressure of the fluid depends only on the density neglecting its impact on the volume and the temperature, i.e., p = p ρ which defines a barotropic fluid. In addition, pressure and density bear a linear relationship
p ρ p = w ρ
where w = p ρ is a dimensionless constant known as equation of state parameter. Substituting Equation (109) in Equation (108), we have another form of energy conservation for the equation of state parameter w,
ρ ˙ ρ + 3 a ˙ a 1 + w = 0
Now, Equations (95), (96), and (108) represent two Friedmann’s Equations, namely, acceleration and evolution equations, and the equation of conservation, respectively. According to this equation, the evolution of all kinds of matter is determined by the conservation of energy and momentum.

Friedmann Equations with Cosmological Constant Λ

We have to incorporate dark matter and dark energy in the matter–energy content due to the significance of their role in current accelerated expansion and the present Minkowskian flat geometry of the universe. Therefore, their role is however unavoidable in the evolution of the universe. The solution of FLRW line element gives the Friedman’s equations using Einstein field equations with cosmological constant Λ written usually in the form
G μ ν + Λ g μ ν = 8 π T μ ν
and Friedmann’s equations with cosmological constant Λ can be worked out
a ˙ 2 a 2 + k a 2 = 8 π G 3 ρ + Λ 3
a ¨ a = 4 π G 3 ρ + 3 p + Λ 3
The equation of energy conservation can also be calculated from these Friedman equations in the presence of cosmological constant Λ . Multiplying Equation (112) with 3 a 2 , differentiating it with respect to time and then dividing by a ˙ , we have
6 a ¨ = 8 π G a 2 ρ + a a ˙ ρ ˙ + 2 Λ a
dividing Equation (114) by a.
6 a ¨ a = 8 π G 2 ρ + a a ˙ ρ ˙ + 2 Λ
Substituting now the 2nd Friedman Equation from Equation (113) in it, we have
6 4 π G 3 ρ + 3 p + Λ 3 = 8 π G 2 ρ + a a ˙ ρ ˙ + 2 Λ
after simplification, we obtain
ρ ˙ + 3 a ˙ a ρ + p = 0
where ρ and p are contributed by all whatever exists and constitutes the universe.

6. A Geometric Object Embedded in an n-Dimensional Space

An object cannot be placed in a space whose dimensions are equal or less than the object to be placed, rather the space must have larger number of dimensions in order to let the object allow rest in it. The presence of an object in a space having larger dimensions than the object is called embedding of it in that space.

6.1. Intrinsic Geometry

The properties of the geometry that we have access to, based on visualization of the two dimensional beings are called intrinsic because two dimensional beings cannot observe how surfaces are shaped in three or higher dimensional spaces.

6.2. Extrinsic Geometry

The properties of the geometry that we have access to, based on visualization of higher dimensional creature are called extrinsic because higher-dimensional creature can observe how surfaces are shaped in three- or higher-dimensional spaces. The geometrical properties related to an object describing how it has been embedded in some higher dimensional space. Extrinsic geometrical properties depend on how the bodies are placed in the space and how they affect it. The geometry which comes into existence due to interaction between space and the body placed in it describes the extrinsic properties. General relativity considers the geometry of spacetime as the extrinsic property of an object and owes its existence due to the body being present in it.

6.3. The Geometry of 2-Sphere Embedded in Three-Dimensional Space

We consider a three-dimensional Euclidean space where three dimensions namely length, width, and height are represented by three coordinate axes, respectively, as we know this space consists of points separate from time, and therefore we do not call its points as events. We assign the triplet of three Cartesian coordinates x , y , z to each point of it, where x, y, and z are measured along the three axes of it. The sketch of embedding the geometry of 2-sphere in three dimensional space is drawn in Figure 10:
The line element in this space is given by
d s 2 = d x 2 + d y 2 + d z 2
Considering now a sphere with its center at the origin of this coordinate system and envisaging its radius to be a, the surface in Cartesian coordinates x , y , z where x, y, and z are along the three axes of three-dimensional Euclidean space. The equation of sphere of this sphere is
x 2 + y 2 + z 2 = a 2
Differentiating Equation (119) with respect to time
2 x d x d t + 2 x d x d t + 2 x d x d t = 0
Moreover, in differential form
2 x d x + 2 y d y + 2 z d z = 0
Solving Equation (121) for d z , we have
d z = x d x + y d y z
Finding the value of z from Equation (119)
z = a 2 ( x 2 + y 2 )
Substituting in Equation (122)
d z = x d x + y d y a 2 ( x 2 + y 2 )
The value of d z = x d x + y d y a 2 ( x 2 + y 2 ) 1 2 comes up with a sort of constraint on d z which despite of being displaced by infinitesimally small amounts d x and d y from an arbitrary point on the surface of the sphere holds us on the surface of the sphere. Squaring d z in Equation (124)
d z 2 = ( x d x + y d y ) 2 a 2 ( x 2 + y 2 )
Putting in Equation (118), the line element takes the form by substituting for d z 2
d s 2 = d x 2 + d y 2 + ( x d x + y d y ) 2 a 2 ( x 2 + y 2 )
The value of the line element in Equation (126) represents the line element for a sphere in terms of Cartesian coordinates x , y , z . We further observe that the line element in Equation (126) has a coordinate singularity at a 2 = x 2 + y 2 in correspondence with the equator of the sphere and in relation to the point A, otherwise at the equator in the intrinsic geometry of 2-sphere there exists no such physical situation. The embedding scenario manifests how the coordinates ( x , y ) cover the whole surface of the sphere uniquely up to this point A. The geometry of 2-sphere in these coordinates becomes geometrically meaningful in three-dimensional Euclidean space. We can transform the line element in Equation (126) above into spherical polar coordinates by taking
x = r sin θ cos ϕ y = r sin θ sin ϕ
where we differentiate each of x and y with respect to θ and ϕ alternately to find
d x = sin θ cos ϕ d r + r cos ϕ cos θ d θ r sin θ sin ϕ d ϕ d y = sin θ sin ϕ d r + r sin ϕ cos θ d θ + r sin θ cos ϕ d ϕ
adding the values of x, y given in Equation (127) after taking square of both equations in it, we get
x 2 + y 2 = r 2 sin 2 θ
adding d x and d y in Equation (128) after taking square of both equations in it, we possess
d x 2 + d y 2 = sin θ d r + r cos θ d θ 2 + r 2 sin 2 θ d ϕ 2
we find the expression
x d x + y d y = sin θ d r + r cos θ d θ x cos ϕ + y sin ϕ r sin θ d ϕ x sin ϕ y cos ϕ
squaring Equation (131), we have
x d x + y d y 2 = sin θ d r + r cos θ d θ x cos ϕ + y sin ϕ r sin θ d ϕ x sin ϕ y cos ϕ 2
Now substituting Equations (129), (130), and (132) in Equation (126) and simplifying to have the following form:
d s 2 = d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
The value of the line element in Equation (133) gives the line element for a sphere in terms of Spherical polar coordinates r , θ , ϕ . The line element in Equation (126) results in an alternative form for
x = ξ cos ϕ y = ξ sin ϕ
d s 2 = a 2 a 2 ξ 2 ) d ξ 2 + ξ 2 d ϕ
The line element in Equation (135) above gives us, in addition, freedom to choose an arbitrary point on the surface of the sphere by ξ = 0 as the origin of the coordinate system. This freedom connotes in it as a hidden symmetry. We can develop ξ and ϕ coordinate curves on the surface of the sphere by generating a standard coordinate system ( ξ , ϕ ) on the tangent plane at the point A that projects vertically downward onto the surface of the sphere. We further observe that the line element in Equation (135) has a coordinate singularity at a = ξ in correspondence with the equator of the sphere in relation to the point A, otherwise at the equator in the intrinsic geometry of 2-sphere there exists no shade of occurrence of such situation. The embedding picture manifests how the coordinates ( ξ , ϕ ) cover the whole surface of the sphere uniquely up to this point A. The geometry of 2-sphere in these coordinates becomes geometrically meaningful in three dimensional Euclidean space.

6.4. The Geometry of 3-Sphere Embedded in Four Dimensional Euclidean Space

Spaces with dimensions higher than three are now significant in mathematical sciences to have proper description of the physical universe. We consider a four-dimensional Euclidean space which can be considered mathematical extension of three-dimensional Euclidean space. Minkowski used a four-dimensional spacetime to explain the phenomena of the physical world as required by special relativity. The structure of Euclidean four-dimensional space is simple as compared to the Minkowskian structure of spacetime. Minkowskian four dimensional spacetime is pseudo-Euclidean space. In four-dimensional Euclidean space, we assign the quadruplet of four Cartesian coordinates x , y , z , w to each point of it, where x, y, z, and w are along the four axes of it. The line element in this space is given by
d s 2 = d x 2 + d y 2 + d z 2 + d w 2
Considering now a sphere with its center at the origin of this coordinate system with radius a, the surface in Cartesian coordinates x , y , z , w where x, y, z and w are along the four axes of four dimensional Euclidean space. The equation of the sphere reads as
x 2 + y 2 + z 2 + w 2 = a 2
Differentiating Equation (137) with respect to time,
2 x d x d t + 2 y d y d t + 2 z d z d t + 2 w d w d t = 0
And in differential form
2 x d x + 2 y d y + 2 z d z + 2 w d w = 0
Finding out the value of d w from Equation (138), we get
d w = x d x + y d y + z d z w
Now finding the value of w from Equation (137)
w = a 2 ( x 2 + y 2 + z 2 )
Substituting in Equation (140), we obtain
d w = x d x + y d y + z d z a 2 ( x 2 + y 2 + z 2 )
The value of d w = x d x + y d y + z d z a 2 ( x 2 + y 2 + z 2 ) 1 2 provides a sort of constraint on d w which, though displaced by infinitesimally small amounts d x , d y , d z from an arbitrary point on the surface of the sphere holds us stuck on the surface of the sphere. squaring Equation (142)
d w 2 = x d x + y d y + z d z 2 a 2 ( x 2 + y 2 + z 2 )
substituting now in Equation (136), the line element takes the form for the value of d w 2
d s 2 = d x 2 + d y 2 + d z 2 + x d x + y d y + z d z 2 a 2 ( x 2 + y 2 + z 2 )
The value of the line element in Equation (144) gives the line element for a sphere in terms of Cartesian coordinates x , y , z , w , We further observe that the line elements in Equation (144) has a coordinate singularity at a 2 = x 2 + y 2 + z 2 in correspondence with the equator of the sphere relative to the point A; otherwise, at the equator in the intrinsic geometry of the 3-sphere, there does not exist any situation like this. The embedding picture manifests how the coordinates ( x , y , z ) cover the whole surface of the sphere uniquely up to this point A. The geometry of 3-sphere in these coordinates becomes geometrically meaningful in four dimensional Euclidean space. We transform the line element in Equation (144) into spherical polar coordinates which are given below,
x = r sin θ cos ϕ y = r sin θ sin ϕ z = r cos θ
where we differentiate x, and y with respect to θ and with respect to ϕ each and differentiate z with respect to θ only to find
d x = sin θ cos ϕ d r + r cos ϕ cos θ d θ r sin θ sin ϕ d ϕ d y = sin θ sin ϕ d r + r sin ϕ cos θ d θ + r sin θ cos ϕ d ϕ d z = r cos θ
Adding x, y, and z in Equation (145) after taking square of all three equations in it to obtain
x 2 + y 2 + z 2 = r 2
Adding now d x , d y , and d z in Equation (146) after taking square of all three equations in it to have
d x 2 + d y 2 + d z 2 = d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
and the expression we calculate
x d x + y d y + z d z = sin θ d r + r cos θ d θ x cos ϕ + y sin ϕ r sin θ d ϕ x sin ϕ y cos ϕ + r cos θ cos θ d r r sin θ d θ
Squaring Equation (149), we have
x d x + y d y + z d z 2 = sin θ d r + r cos θ d θ x cos ϕ + y sin ϕ r sin θ d ϕ x sin ϕ y cos ϕ + r cos θ cos θ d r r sin θ d θ 2
Now, substituting Equations (147), (148), and (150) in Equation (144), after simplification we get
d s 2 = a 2 a 2 r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
It can further be expressed in the form
d s 2 = 1 1 r 2 a 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
It is important to note here that for a in the Equation (152) above. It reduces to the metric of ordinary three-dimensional Euclidean space
d s 2 = d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
which we calculated in Equation (133). The metric in Equation (151) has a singularity at r = a , which is just a coordinate singularity and has nothing to do with physical reality of the sphere as we can observe. The line element in Equation (152) results in an alternative form for
x = ξ cos ϕ y = ξ sin ϕ
d s 2 = a 2 ( a 2 ξ 2 ) d ξ 2 + ξ 2 d ϕ
The line element in Equation (155) gives us, in addition, freedom to choose an arbitrary point on the surface of the sphere by ξ = 0 as the origin of the coordinate system. This freedom is implied by it as a hidden symmetry in it. We can develop ξ and ϕ coordinate curves on the surface of the sphere by generating a standard coordinate system ( ξ , ϕ ) on the tangent plane at the point A that projects vertically downward onto the surface of the sphere. We further observe that the line element in Equation (155) has a coordinate singularity at a = ξ with respect to the equator of the sphere in relation with the point A, otherwise at the equator in the intrinsic geometry of 2-sphere there exists no such situation. The embedding picture manifests how the coordinates ( ξ , ϕ ) cover the whole surface of the sphere uniquely up to the point A. The geometry of the 2-sphere in these coordinates becomes geometrically meaningful in three dimensional Euclidean space.

7. Einstein’s Static Universe

Albert Einstein himself applied general relativity to the largest scale of spacetime [3] and presented the very first relativistic model of the universe laying the foundations of modern theoretical cosmology. The model was later on called as Einstein world or universe. For this purpose, Einstein modified his field equations by proposing an inbuilt energy density known as cosmological constant Λ in the geometrical structure of spacetime itself that provides repulsive gravity to keep the universe from expanding
G μ ν + Λ g μ ν = 8 π T μ ν
Equation (156) when solved for the most homogeneous and isotropic geometry of FLRW spacetime produces Friedmann equations as we derived earlier
a ˙ a 2 + k a 2 = 8 π G 3 ρ + Λ 3
a ¨ a = 4 π G 3 ρ + 3 p + Λ 3
As for a static universe H = 0 , which implies that a ¨ a = 0 . Now a static universe possesses cold matter which means it does not has pressure, i.e., p = 0 , so Equations (157) and (158) reduce to the form, respectively,
8 π G 3 ρ + Λ 3 k a 2 = 0
and
4 π G 3 ρ + Λ 3 = 0
From above Equation (160), we have
Λ = 4 π G ρ
Substituting this value of Λ in Equation (159), and having the equation simplified, we again get the value of Λ in terms of the curvature term k and the scalar factor a t , that is
Λ = k a 2
The line element for the static Einstein universe can be written now using FLRW metric. From above Equation (162) for k = + 1 , we have a 2 t = Λ 1 , substituting in Equation (59), the static solution for closed universe becomes
d s 2 = d t 2 + Λ 1 1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
Using the Schwarzschild coordinates with the re-scale of radial coordinate and by defining R = r a , we have
d s 2 = d t 2 + 1 1 Λ R 2 d R 2 + R 2 d θ 2 + R 2 sin 2 θ d ϕ 2
Case-I (empty universe)
substituting Λ = 0 in Equation (161) gives 4 π G ρ = 0 ρ = 0 which implies that k = 0 a Euclidean flat universe. It does not belong to Einstein static universe because it is empty.
Case II (non-empty universe)
Einstein universe belongs to Λ 0 and ρ 0 implying that k > 0 which represents a universe with hypersurface of Riemannian geometry. In Einstein’s universe ρ > 0 , therefore a positive cosmological constant Λ > 0 would be allowed which also implies k > 0 .

7.1. Instability of Einstein’s Universe

Equation of energy conservation can be had from Equations (157) and (158) by multiplying Equation (158) with 3 a 2 , differentiating it with respect to time and then dividing by a ˙ , we have
6 a ¨ = 8 π G a 2 ρ + a a ˙ ρ ˙ + 2 Λ a
Dividing Equation (163) by a and substituting the 2nd Friedman Equation from Equation (159) in it, we have
6 a ¨ a = 8 π G 2 ρ + a a ˙ ρ ˙ + 2 Λ
6 4 π G 3 ρ + 3 p + Λ 3 = 8 π G 2 ρ + a a ˙ ρ ˙ + 2 Λ
after simplification, we get
ρ ˙ + 3 a ˙ a ρ + p = 0
For the cold matter universe p = 0 , with this the resulting equation is a separable universe
ρ ˙ + 3 a ˙ a ρ = 0
ρ ˙ ρ d t = 3 a ˙ a d t
ln ρ = 3 ln a + ln Z
ln ρ = ln a 3 + ln Z
ln ρ ln a 3 = ln Z
ln ρ a 3 = ln Z
which gives
ρ a 3 = Z
where Z is some positive constant of integration Z > 0
ρ = Z a 3
Further, as the universe does not expand so that a ( t ) = a ( t 0 ) = a 0 , therefore replacing a ( t ) with a 0 in Equation (174)
ρ = Z a 0 3
Substituting the value of ρ from Equation (176) in Equation (158), i.e., 2nd Friedmann Equation with p = 0 , we obtain
a ¨ a = 4 π G Z 3 a 0 3 + Λ 3
and substituting in Equation (161) gives
Λ = 4 π G ρ = 4 π G Z a 0 3
where 4 π G Z a 0 3 > 0 since Z > 0 , Now, we perturb the solution slightly with the following perturbation
ε < < 1
a ( t ) = a ( t 0 ) + ε ( t ) a ( t 0 ) = a 0 ( 1 + ε ( t ) )
substituting this in Equation (175), we have
d 2 d t 2 a 0 ( 1 + ε ( t ) ) a 0 ( 1 + ε ( t ) ) = 4 π G Z 3 a 0 ( 1 + ε ( t ) ) 3 + Λ 3
Or
d 2 d t 2 a 0 ( 1 + ε ( t ) ) = 4 π G Z 3 a 0 2 ( 1 + ε ( t ) ) 2 + Λ 3 a 0 ( 1 + ε ( t ) )
Using the Maclaurin series expansion as ε < < 1 , and ( 1 + ε ( t ) ) 2 = 1 2 ε + O ( ε 2 ) , Now Equation (181) becomes by neglecting O ( ε 2 ) as ε < < 1 , so that
a 0 d 2 ε d t 2 = 4 π G Z 3 a 0 2 ( 1 2 ε ( t ) ) + Λ 3 a 0 ( 1 + ε ( t ) ) + O ( ε 2 )
ε ¨ = 4 π G Z 3 a 0 3 + 8 π G Z 3 a 0 3 ε + Λ 3 ε + Λ 3 = 8 π G Z 3 a 0 3 + Λ 3 ε + Λ 3 4 π G Z 3 a 0 3
Using the value of Λ = 4 π G Z a 0 3 from Equation (177) in Equation (183), it can be expressed in the form
d 2 ε d t 2 Λ ε = 0
As the cosmological constant is Λ > 0 , the solution of above equation will read as
ε = P exp ( Λ t ) + Q exp ( Λ t )
Due to existence of the 1st term in the above solution as positive and in the case of an arbitrary perturbation considered initially, both of the constants P 0 , Q 0 will help the perturbation grow and it will not remain small which will imply that the static solution is unstable, although P = 0 can be possible only for specialized initial conditions such as singular one.

7.2. De Sitter Universe

In Einstein’s static model with positive cosmological constant when energy density of the matter is removed de Sitter model results. The de Sitter model of the universe presented in 1917 was proposed just after Einstein presented his static closed model of the universe. Einstein resorting to the Mach’s principle was of the view that it is merely matter density in universe that is the cause of inertia and gravitation. For checking the status of this Einstein’s belief de Sitter posed the 2nd model of the universe devoid of matter density T μ ν = 0 , however retaining the cosmological constant that is G μ ν = g μ ν Λ . The de Sitter model is the maximally symmetric solution of Einstein’s field equations with vanishing matter density. The geometric theoretic structure of spacetime of the de Sitter model is comparatively more complicated than that of Einstein’s model of the universe. The characteristic of the de Sitter model is that it predicts redshift despite it contains neither matter density nor radiation. We review de Sitter model using Fiedmann’s equations, however it is important to note that these equations were worked out after the development of de Sitter model. We derived Friedmann equations above in the presence of cosmological constant term Λ , which are
a ˙ a 2 + k a 2 = 8 π G 3 ρ + Λ 3
a ¨ a = 4 π G 3 ρ + 3 p + Λ 3
de Sitter universe corresponds to ρ = 0 , so that k ρ = 0 , Equation (186) takes the form
a ˙ a = Λ 3
Integrating with respect to time
a ˙ ( t ) a ( t ) d t = Λ 3 d t
a ( t ) = e Λ 3 t
From Equation (188), H = a ˙ a , so the Equation (190) can be expressed as
a ( t ) = e H t
Which corresponds to the modified Einstein field equations
G μ ν = g μ ν Λ

8. The Conformal FLRW Line Element

The metric in Equation (57) can be conformally recast by defining conformal time as
d τ = d t a ( t )
so that
d t = a ( t ) d τ
After substituting Equation (194) in Equation (57) and simplifying, we get the line element in the form
d s 2 = a ( τ ) 2 d τ 2 1 k r 2 1 d r 2 r 2 Ω 2
Due to conformal time the scale factor a τ becomes a factor of spatial as well as temporal components in the metric. Now, a function f t which depends upon time can be differentiated as
f ˙ t = f τ a τ f ¨ t = f τ a 2 τ f τ a 2 τ H
where dot “.” and “,” represent derivatives with respect to cosmic and conformal times, respectively, and H = a τ a τ . Now, replacing f t and its derivatives with a t both in correspondence with cosmic time ‘t’ and conformal time ‘ τ
a ˙ t = a τ a τ
a ¨ t = a τ a 2 τ H 2 a τ
and
H = a ˙ a = H a τ
H ˙ = H a 2 τ H 2 a 2 τ
Similarly
H 2 = 8 π G a 2 3 ρ k
and
ρ + 3 H ρ + p = 0
Now we solve the energy conservation equation From Equation (108)
ρ ˙ + 3 a ˙ a ρ + p = 0
in order to get the relation between energy density ρ , scale factor a and equation of state parameter w = p ρ we solve
ρ ˙ = 3 a ˙ a ρ + p = 3 a ˙ a ρ 1 + p ρ
ρ ˙ ρ = 3 a ˙ a 1 + w
where p ρ = w . Integrating Equation (205)
1 ρ d ρ = 3 1 + w 1 a d a
which gives
ρ = a 3 1 + w
Now from 1st Friedmann equation, after simplification and doing integration, we find
a = t 2 3 1 + w
For w = 1 , 0, 1 3 , we find pressure, energy density and scale factor characterizing the expansion of the universe which depicts three phases of the universe namely vacuum dominated, radiation dominated and matter dominated, respectively.

8.1. Vacuum Domination ( Λ -Dominated Era)

For w = 1
ρ = a 3 1 + w = a 0
and
a = t 2 3 1 1 = t

8.2. Radiation Domination

For w = 1 3
ρ = a 3 1 + 1 3 = a 4
and
a = t 2 3 1 + 1 3 = t 1 2

8.3. Matter Domination

For w = 0
ρ = a 3 1 + 0 = a 3
and
a = t 2 3 1 + 0 = t 2 3

8.4. Critical Density ( ρ c ) and Density Parameter ( Ω )

Now from 1st Friedman Equation (112) with Λ = 0 and H = t ln a , we relate the curvature of spacetime k and the expansion characterized by the scale factor a t to the energy density ρ t of the universe and find the expression for the critical density required to keep the current rate of the expansion.
H 2 = 8 π G 3 ρ k a 2
For critical density ρ c the curvature of spacetime geometry k must vanish, so that Equation (215) reduces to
H 2 = 8 π G 3 ρ
where we obtain the expression for critical density
ρ = ρ c = 3 H 2 8 π G
From Equation (215) dividing both sides by H 2 and rearranging
1 = 8 π G 3 H 2 ρ k a 2 H 2 = ρ 3 H 2 8 π G k a 2 H 2
where 3 H 2 8 π G = ρ c , therefore Equation (218) becomes
1 = ρ ρ c k a 2 H 2 = Ω k a 2 H 2
Ω 1 = k a 2 H 2
where Ω = ρ ρ c is the density parameter and we can predict in terms of it about the geometry of universe. The local geometry of the universe is investigated by this parameter by observing whether the relative density is smaller than unity, greater than or equal to it. In the Figure 11 all three geometries are represented as the density parameter would allow:
Equation (220) can also be derived from Equation (215) in an alternative style. Writing Equation (215) by multiplying and dividing the 1st term on the right side with ρ c
H 2 = 8 π G ρ 3 ρ c ρ c k a 2
Using the density parameter Ω = ρ ρ c , in Equation (221) we can write
H 2 = 8 π G 3 ρ c Ω k a 2
Now, from the critical density expression in Equation (217),
3 8 π G = ρ c H 2 8 π G 3 = H 2 ρ c
Substituting the 2nd part in Equation (223) in Equation (222) and using the density parameter, we get
H 2 = H 2 Ω k a 2
which gives the following form similar to Equation (220)
Ω 1 = k a 2 H 2
Now
Ω = ρ ρ c
is considered decisive in describing the evolution of the universe. The present value of it is denoted by Ω 0 and it gives following three geometries of the universe
Ω 0 > 1
a closed universe implying the universe with spherical geometry
Ω 0 < 1
an open universe implying the universe with hyperbolic geometry and
Ω 0 = 1
a flat universe implying the universe with Euclidean or Minkowskian geometry. The present value of critical density can be calculated with present value of Hubble constant H 0 , gravitational constant G and π .
ρ c , 0 = 3 H 0 2 8 π G = 3 73.8 2 8 3.14 6.67 × 10 11 = 1.1 × 10 5 h 2
where the scaled Hubble parameter h is defined by H = 100 hkm s 1 Mp c 1 and H 1 = 9.778 h 1 Gyr H 1 = 2998   h 1 Mpc.

8.5. Particle Horizon

When the scale factor a t is multiplied with the co-moving coordinates we get the proper distance. In cosmology causality is one directional since we just receive photons from the outer world that serves to be self-sufficient approach. The horizon or horizon distance of the universe is defined as the maximum distance that light could have traveled to our reference Earth as the time after the beginning of the universe when for the first time it became exposed to electromagnetic radiation [45], thus horizon represents the causal distance in the universe.
d H ( t ) = a ( t ) 0 t d t a ( t )
Such that d H ( t ) H 1 ( t ) Particle horizon is defined to be the distance traveled by a photon from the time of big bang up to a certain later time, t. Particle horizon puts limits on communication from the deep inward past.

8.6. Event Horizon

An event horizon defines such a set of points from which light signals sent at some given time will never be received by an observer in the future. It sets limits on the horizon distance and on communication to the future so that as long as it exists, the size of the causal patch of the universe will be finite.

8.7. Deceleration Parameter ( q 0 )

A Taylor series is a series expansion of a function about a given point. We require here a one dimensional Taylor series which is the expansion of a real function f x about a point x = a and is given by
f t = f x x = a = f a + f a x a + f a 2 ! x a 2 + f a 3 ! x a 3 + · · · + f n a n ! x a n + · · ·
We take the function f x = a t which is scale factor and find its Taylor series expansion about the present time t = t 0
a t = a t t = t 0 = a t 0 + a ˙ t 0 t t 0 + a ¨ t 0 2 ! t t 0 2 + t 0 3 ! t t 0 3 + · · · + a n t 0 n ! t t 0 n + · · ·
dividing Equation (233) by a t 0 throughout, we have
a t a t 0 = a t 0 a t 0 + a ˙ t 0 a t 0 t t 0 + 1 2 a ¨ t 0 a t 0 t t 0 2 + 1 6 t 0 a t 0 t t 0 3 + · · · + 1 n ! a n t 0 a t 0 t t 0 n + · · ·
ignoring the higher terms we have the following remaining expression
a t a t 0 = 1 + a ˙ t 0 a t 0 t t 0 + 1 2 a ¨ t 0 a t 0 t t 0 2
multiplying and dividing now by a ˙ t with 3rd term of Equation (235) on the right hand side:
a t a t 0 = 1 + a ˙ t 0 a t 0 t t 0 + 1 2 a ˙ t 0 a t a ¨ t 0 a ˙ t 0 t t 0 2
Multiplying again the 3rd term on the right hand side of Equation (236) with a ˙ t 0 a t 0 and its reciprocal a t 0 a ˙ t 0 , we have
a t a t 0 = 1 + a ˙ t 0 a t 0 t t 0 + 1 2 a ˙ t 0 a t 0 × a t 0 a ˙ t 0 a ˙ t 0 a t a ¨ t 0 a ˙ t 0 t t 0 2
Putting for a ˙ t 0 a t 0 = H 0 , the present value of Hubble parameter and a t 0 a ¨ t 0 a ˙ t 0 2 = q 0 , Equation (237) reduces to the following:
a t a t 0 = 1 + H 0 t t 0 + 1 2 H 0 2 q 0 t t 0 2
where
q 0 = a t 0 a ¨ t 0 a ˙ t 0 2 = a ¨ t 0 a ˙ t 0 H 0 1 = a ¨ t 0 a t 0 H 0 2
is called the deceleration parameter. It tells us that greater the value of q 0 , the faster will be speed of deceleration. It can be further related with the acceleration equation
a ¨ t a t = 4 π G 3 ρ + 3 p
Putting Equation (240) in Equation (239)
q 0 = 4 π G 3 ρ + 3 p H 0 2
With p = 0 for a universe having matter domination and present energy density ρ = ρ 0 with dividing and multiplying by 2, we possess
q 0 = 1 2 8 π G 3 H 0 2 ρ 0
Now, as the critical density is given by ρ c = 3 H 0 2 8 π G from the 1st Friedmann equation. Therefore Equation (242) takes the form
q 0 = 1 2 1 ρ c ρ 0 = 1 2 ρ 0 ρ c = 1 2 Ω 0
The measurement of deceleration parameter q 0 determines how much bigger the universe was in earlier times. The explorations of redshift measures of supernovae of Type S N I a to measure the value of q 0 has shown astoundingly that q 0 < 0 at the present which means that the expansion of the universe is accelerating rather than to be decelerating which affirms that the concept of dark energy must be acknowledged. Accelerated expansion of the universe corresponds to q 0 < 0 , whereas q 0 > 0 corresponds decelerated expansion. It is interesting to notice that for all of these components we have H > 0 , i.e., an increasing scale factor which gives the expansion rate of the universe. Moreover, to get a better understanding of the properties of each species, it is useful to introduce the deceleration parameter q 0 as
q 0 = a ¨ a ˙ a ˙ 2 = a ¨ a ˙ a a ˙ = a ¨ a ˙ H 1
such that for both matter-dominated or radiation-dominated universe the expansion is decelerating. It is also interesting to notice that components with w < 1 3 give an accelerated expansion.

8.8. Friedmann Equations in Terms of Density Parameter

We found earlier Friedmann equations
a ˙ a 2 + k a 2 = 8 π G 3 ρ
a ¨ a = 4 π G 3 ρ + 3 p
In Equation (245) in order to incorporate vacuum energy we can write energy density as the sum of all energy components ρ = ρ m + ρ r + ρ Λ such that the equation can be written as
H 2 = 8 π G 3 ρ m + ρ r + ρ Λ k a 2
where a ˙ a = H is the Hubble parameter, writing ρ Λ as ρ Λ = Λ = Λ 8 π G and ρ = ρ m + ρ r which further can be written as the contributing ingredients ρ m = ρ b + ρ C D M and ρ r = ρ γ + ρ ν , also we found earlier the critical density to be 3 H 2 8 π G = ρ c d , which for the present value can be expressed as ρ c , 0 = 3 H 0 2 8 π G we find from it the value of 8 π G = 3 H 0 2 ρ c d , 0 and substitute in Equation (247) which comes to be
H 2 = H 0 2 ρ m ρ c , 0 + ρ r ρ c , 0 + ρ Λ ρ c , 0 k a 2
or
H 2 = H 0 2 Ω m , 0 + Ω r , 0 + Ω Λ , 0 k a 2
where
Ω m , 0 = ρ m ρ c , 0 , Ω r , 0 = ρ r ρ c , 0 , Ω Λ , 0 = ρ Λ ρ c , 0
It might be suitable to write the curvature term k in terms of density parameter k = Ω k , 0 = ρ k ρ c , 0 , further the present value of the scale factor a t = 1 so that Equation (249) takes the form
H 2 = H 0 2 Ω m , 0 + Ω r , 0 + Ω Λ , 0 Ω k , 0
Now, for the present value of Hubble parameter, i.e., H = H 0 , Equation (251) can be written for the curvature density parameter
Ω k , 0 = H 0 2 Ω m , 0 + Ω r , 0 + Ω Λ , 0 1
Equation (249) can be written in general form, i.e., H H 0 and a a 0 = 1
H 2 = H 0 2 Ω m + Ω r + Ω Λ Ω k a 2
Equation (253) can also be written for the present values of all the energy density parameters
H 2 = H 0 2 Ω m , 0 + Ω r , 0 + Ω Λ , 0 Ω k , 0 a 2
We know that energy density ρ for matter, radiation and vacuum domination eras changes with the scale factor that characterizes the expansion of the universe according to
ρ a 3 ρ a 4 ρ a 0
respectively. Thus, Equation (254) takes the following form using Equation (255):
H 2 = H 0 2 Ω m , 0 a 3 + Ω r , 0 a 4 + Ω Λ , 0 a 0 Ω k , 0 a 2
The Equation (256) represents Friedmann equation in terms of density parameters.
For a = a t a t 0 = 1 1 + z , Equation (256) can be expressed in terms of redshift as follows
H 2 = H 0 2 Ω m , 0 1 + z 3 + Ω r , 0 1 + z 4 + Ω Λ , 0 1 + z 0 Ω k , 0 1 + z 2
We can discuss various models for the universe using Equation (256) for matter, radiation, Λ and curvature-dominated eras.
For matter domination, Equation (256) with Ω m , 0 = 1 and with the rest of terms vanishing gives
a = 3 2 H 0 t 2 3 t = 2 a 3 2 3 H 0
which gives an expanding universe with expansion rate inversely proportional to time, i.e., H = 2 3 t 1 and age of the universe would be t 0 = 2 3 H 0 1 . Such model must be subject to deceleration as the time goes on.
For radiation domination, Equation (256) with Ω r , 0 = 1 and with the rest of terms vanishing, gives
a = 2 H 0 t t = a 2 2 H 0
which gives an expanding universe with expansion rate inversely proportional to time, i.e., H = 1 2 t 1 and age of the universe would be t 0 = 1 2 H 0 1 . The expansion is subject to deceleration in this radiation-dominated era.
For Λ domination, Equation (256) with Ω Λ , 0 = 1 and the rest of terms vanishing gives
a = e H 0 t t =
which gives an exponentially expanding universe with expansion rate inversely proportional to time, i.e., H = ln a t 1 and infinite age.
For k domination or the otherwise empty universe, Equation (256) with Ω k , 0 = 1 and with the rest of terms vanishing gives
a = H 0 t t =
which gives an expanding universe with expansion rate inversely proportional to the time, i.e., “t”.

8.9. Cosmological Redshift

we considering the FLRW geometry
d s 2 = d t 2 + a ( t ) 2 1 k r 2 1 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
Note here that the coordinates r , θ , ϕ in the metric Equation (262) are comoving spatial coordinates; therefore, galaxies which are considered as point particles constituting the particles of cosmological fluid in cosmology remain at fixed coordinates and it is the geometry of the spacetime that expands itself and is characterized by the scale factor a t completely. Three intervals—spacelike, timelike and lightlike, or null—expressed in the form d s 2 > 0 , d s 2 < 0 , or d s 2 = 0 , respectively. In the spacetime geometry light propagates following the interval d s 2 = 0 or d s = 0 which means that it does not travel at all any distance through the spacetime. We consider a ray of light propagating along the radius as all the points in space are equivalent at a given time from some zero value radius to some certain value of it in later times. As the light ray travels radially therefore only one spatial dimension is retained and the vanishing of time dimension follows from d s = 0 and other two spatial dimensions vanish due to radial propagation of light therefore d t = d θ = d ϕ = 0 , then Equation (262) gives
0 = d t 2 + a ( t ) 2 1 k r 2 1 d r 2
or
d t a ( t ) = 1 1 k r 2 d r
In order to calculate the total time elapsed from r = 0 to some certain later time value r = r 0 , we shall integrate Equation (264) between emission and reception times t e and t r respectively.
t = t e m i t = t r e c 1 a ( t ) d t = r = 0 r = r 0 1 1 k r 2 d r
A ray of light, now, given off after a short interval of time d t e m i so that time of emission of light ray becomes t e m i + d t e m i and accordingly we can have the time of reception to be t r e c + d t r e c from an integral of the same nature given in Equation (265) because of comoving coordinates the galaxies remain at the same coordinates, so that
t = t e m i t = t r e c 1 a ( t ) d t = t = t e m i + d t e m i t = t r e c + d t r e c 1 a ( t ) d t
now
t = t e m i t = t e m i + d t e m i 1 a ( t ) d t = t = t r e c t = t r e c + d t r e c 1 a ( t ) d t
For deriving the redshift relation as the universe expands we use Figure 12 that is given below:
The slices are very narrow, so the area is just the area of a rectangle, i.e., width times height, i.e.,
d t r e c a t r e c = d t e m i a t e m i
For an expanding universe
a t r e c > a t e m i ,
it implies from Equation (268) d t r e c > d t e m i that as the universe expands the time interval between two rays increases. We consider now successive crests or troughs of a single ray instead of two rays as we did earlier so that wave length λ is directly proportional to the time interval between two successive crests or troughs λ d t and d t a ( t ) and we have
λ r e c λ e m i = a t r e c a t e m i
We define now the redshift
1 + z a t r e c a t e m i

8.10. Luminosity ( L ), Brightness, Luminosity Distance ( d L ) and Angular Diameter Distance  d A

We can deduce relations from the properties of electromagnetic radiation and the quantities contained in FLRW line element. The velocity of electromagnetic waves is constant and finite. Light and electromagnetic radiation acts as cosmological messenger and all the distances measured cosmologically are extracted from the properties of it. The velocity of light being finite has to take time to reach us and universe might have expanded significantly during this time.

8.11. Luminosity L

Luminosity is defined as the absolute measure of the electromagnetic power or energy radiated per unit time by an astronomical object like star, galaxy or cluster of galaxies. It is denoted by L and is measure in Joule per second J s 1 which is also known as watts. Usually luminosity is measured in terms of the luminosity of the sun denoted by L .

8.12. Brightness

It refers to how bright an object appears to an observer and depends upon luminosity, distance between the observer and the object and absorption of light along the path between observer and the object.

8.13. Luminosity Distance ( d L )

We consider a point source S radiating electromagnetic light equally in all directions spherically; the amount of light passing through elements of surface areas varies with the distance of it from the light source.
Given below in Figure 13, the light of luminosity L is being radiated. We consider a spherical hollow centered on the point source S as shown in the Figure 13 below The interior of a hollow sphere gets illuminated thoroughly. As the radius of the sphere increases, the surface area of the imagined hollow sphere also increases, such that a constant or absolute measure of luminosity has to spread in expanding sphere illuminating it, i.e., as the radius increases the constant luminosity has more and more surface area to illuminate which leads to decrease in the observed brightness. If an observer at a distance equivalent to the radius of sphere receives the electromagnetic radiation L per unit time and F be the energy flux per unit time per unit area from the source or the point source, say O, then in Euclidean geometry we will have
F = L A = L 4 π r 2
where F = Flux density of the illuminated sphere, L = luminosity, and A = area of the illuminated sphere From Equation (272) for r = d L
F = L 4 π d L 2
which gives
d L = L 4 π F
We next look how the luminosity distance is related with expansion of the universe. In expanding sphere we might have its radius as the product of scale factor and the radius, i.e., a ( t ) r = a ( t ) d L , so that the energy emitted gets diluted
4 π r 2 4 π a t r 2
and a photon loses energy as F a ( t e ) a ( t 0 ) and redshift relation we have 1 + z = λ ( t 0 ) λ ( t e ) = a ( t 0 ) a ( t e ) which implies F a ( t e ) a ( t 0 ) 1 1 + z Equation (273) becomes
F = L 4 π a ( t ) r 2
further
a ( t e ) a ( t 0 ) = L 4 π a ( t ) d L 2
If L is known for a source, it is known as standard candle. Supernovae type I a were used as standard candles for larger cosmic redshifts which led to accelerated expansion.

8.14. Angular Diameter Distance d A

It is the ratio of the proper distance measured when the light left the surface of an object to the later measured distance by redshifting of light in some later time. Certainly the redshift of light measured would be smaller measured at the time when the light left the surface of the object to be measured in later times. The schematic diagram is shown in Figure 14 for angular diameter distance:
It is defined in terms of objects physical distance known as proper distance and the angular size of the object seen from the surface of earth. If size of the source be S and angular size θ , then
θ = S D A
where D A is the angular diameter distance of the source. From FLRW line element for photons d r 2 d ϕ 2 0 , we have
d s 2 = a 2 ( t ) r 2 d θ 2
d s = D A = a ( t ) r d θ
d θ = D A r a ( t ) = d A r a ( t )

9. Problems Faced by the Standard Model of Cosmology

From 1st Friedmann equation
a ˙ a 2 + k a 2 = 8 π G 3 ρ
we see that curvature k is negligible depending on observation and Ω 1 which means it would have been created tuned finely in the very early universe. From 2nd Friedmann equation
a ¨ a = 4 π G 3 ρ + 3 p
we see that if ρ + 3 p remains positive, the acceleration is negative which means that the expansion of the universe will go on slowing down. Further, far flung parts of the universe display the same properties as observation evidence despite the fact that they have not been in causal contact with each other.

9.1. Monopole Problem

The problem is about the question of why do we not observe magnetic monopoles in the universe today. It results from combining the big bang model with GUT in particle physics, thus it is related to particle cosmology where during symmetry breaking phase transitions are considered. In the very early universe, when the phase transitions are considered to occur, it is expected that these phase transitions will create magnetic monopoles with enormous energy density which might dominate the total energy density of the universe. During symmetry breaking when phase transitions take place, these give rise to flaws known to be as topological defects. GUT predict that during GUT phase transitions these point-like topological defects are created which act as magnetic monopoles. It is considered that the radiation and matter dominated eras could not take place as these monopoles do not get diluted as they are supposed to be non-relativistic and their energy density would decay like a 3 [46], but as we observe the universe evolved to the later eras so question arises how this occurred which is at the heart of this problem.

9.2. Horizon Problem

On the basis of the standard big bang model it is difficult to understand the uniform distribution of the temperature of CMB to 1 part in 10 5 . The horizon problem is related with the issue of the causal contact as it has been revealed by the uniform distribution of temperature of the cosmic background radiation (CMB) across all parts of the universe. In order to understand the problem we have to understand the horizon size and causal contact or communication. At any instant of time horizon size is defined as the largest distance, i.e., maximal distance over which two events could be in causal with each other. Therefore it is the maximum distance a photon could have traveled since the birth of the universe or since the time when the universe became transparent. It can be found from the FLRW metric to be d s 2 = R H = c 0 t d t a ( t ) which reveals the fact that size of the horizon depends upon the history of the universe as it evolves through time. It is also called comoving horizon as causal contact develops between two events and the universe is expanding so that they are getting separated apart mutually. In the standard big bang theory the universe was matter dominated at the time of last scattering ( t l s ) so that the horizon distance at that time can be approximated by the value d H t l s = 2 c H 1 t l s . Now, the Hubble distance at the time of last scattering was c H 1 t l s 0.2 Mpc and the horizon distance at last scattering was d H t l s 0.4 Mpc . Therefore, the points which were separated more than 0.4 Mpc distance apart at the time of last scattering ( t l s ) were not connected causally in the big bang scenario. Further, angular diameter distance d A to the last scattering surface is 13 Mpc; therefore. points on the last scattering surface that were separated by a horizon distance shall have angular separation θ H = d H t l s d A 0.4 Mpc 13 Mpc 0.03 rad 2 0 as viewed today from the Earth. It means that the points separated by an angle as small as ∼ 2 0 on the last scattering surface were not in causal contact with each other when CMB emitted with temperature fluctuations. However, we come to know that δ T T is as small as 10 5 on the scales with angular separation θ H > 2 0 . So here we state the problem that the regions which were not connected through causal contact with each other at the time of last scattering have similar properties homogeneously.

9.3. Flatness Problem

When we consider Friedmann’s equations evolve in a universe where only radiation and baryonic matter exist without vacuum energy density present there, then flatness problem arises in such a universe [47]. From the 1st Friedmann Equation
H 2 = 8 π G 3 ρ k a 2
1 = 8 π G 3 H 2 ρ k a 2 H 2 = ρ 3 H 2 8 π G k a 2 H 2
where 3 H 2 8 π G = ρ c , therefore Equation (285) becomes
1 = ρ ρ c k a 2 H 2 = Ω k a 2 H 2
Ω 1 = k a 2 H 2
so that the spatial curvature of the universe is related to the density parameter Ω through Friedmann’s equation. Observational evidence shows that the universe is nearly flat today, i.e., ρ = ρ c , Ω = ρ ρ c 1 . This means that the value of Ω would have to be very close to 1 at Planck era t p l . This means that the initial conditions of the universe were tuned finely. Because of this, the flatness problem is also known as the fine-tuning problem and the flatness problem arises because in a comoving volume the entropy remains conserved. Further, from Equation (284) above, the energy density of the universe without considerations of vacuum energy as is the case of big bang model is ρ = ρ R + ρ M and we can also write
H 2 = 8 π G 3 ρ R + ρ M k a 2
The term k a 2 is clearly proportional to a 2 , while the energy density terms ρ R and ρ M fall off faster than scale factor a ( t ) , i.e., ρ R 1 a 3 ( t ) and ρ M 1 a 4 ( t ) . This ratio k a 2 ( t ) 8 π G 3 ( ρ R + ρ M ) = k a 2 ( t ) ρ 3 M p l 2 then is much smaller than unity when the scale factor a ( t ) has increased by a factor of 10 30 since the Planck era.

9.4. Entropy Problem

The adiabatic expansion of the universe following the first law of thermodynamics is related to the flatness problem [48] discussed above. Temperature plays a significant role in the early universe because at early epochs the age and expansion rate H = t ln a are described in terms of it with the number of relativistic degrees of freedom. From the 1st Friedmann’s equation we have the expression for density parameter Ω 1 = k a 2 ( t ) H 2 and expansion rate in radiation-dominated era in terms of temperature is H ρ R 2 8 π G T 4 = T 4 M p l 2 , so that the density parameter expression becomes Ω 1 = k M p l 2 a 2 ( t ) T 4 . Now the entropy density is s T 3 and the entropy per commoving volume S a 3 ( t ) s a 3 ( t ) T 3 and we have Ω 1 = k M p l 2 S 2 3 T 2 . The entropy per co-moving volume S remains constant throughout the evolution of the universe as the hypothesis of adiabaticity requires so that we Ω 1 t = t p l = ( 1 ) M p l 2 S U 2 3 T p l 2 10 60 . It comes clear that at early epochs the value of Ω 1 is very close to zero as the total entropy of the universe is very large.

10. Introduction to Inflation

Inflation is the period of superluminally accelerated expansion of the universe taking place sometime in the very early history of the universe. It is now a widely accepted paradigm which is described as the monumental outgrowth gushing out during the tiniest fraction of the first second between ( 10 36 10 32 ) s. Inflation maintains that just after the occurrence of the big bang, exponential stretching of spacetime geometry took place, i.e., becoming twice in size again and again at least about (60–70) times over before slowing down. Alexei Strobinsky approached the exponentially expanding phase in the early universe by modifying Einstein Field Equations whereas Alan Guth approached the scenario in the realm of particle physics proposing a new picture of the time elapsed in the very small fraction of the first second in the 1980. He suggested that the universe spent its earliest moments growing exponentially faster than it does today. There is a large number of inflation models in hand today but every model has its own limitations to draw the true picture of what happened actually in the early universe.
As the theory of inflation is recognized today, it has myriad models describing inflationary phase in the early universe. From among the heap of these competing models which differ slightly from one to the other, no model can claim a complete and all-embracing prospectus of what happened actually in the universe so that the fast expansion of or in spacetime takes place. All the energy density that can be adhered to the early exponentially expanding phase of the universe was in the very fabric of spacetime itself despite ti being in the form of radiation or particles. The early accelerating phase can be now best described with de Sitter model with slightly broken time symmetry. With the creation of spacetime that purports to be the earliest patch of the universe that comes to being would be stretched apart in an incredibly small time span of the order of a tiniest fraction of first second to such a colossally larger size that its geometry and topology would be hardly indiscernible from Euclidean geometry. It will logically ensue similar initial conditions for the energy density to be dispersed at every point in the fabric of spacetime and the same will be the condition of temperature in this early phase. That’s why the quantum fluctuations which seed in later times the structure formation in the universe impart the uniform temperature to all parts of the universe thereby resolving the homogeneity problem of the universe. This is because all the quantum fluctuations which cause the observable universe were once causally connected in the deep past of the universe. It might have attained a highest temperature which was within or lesser than the limits of Planck scale ( 10 19 GeV). The energy scale mentioned earlier when the inflation comes to an end and transforms into the uniform, very hot, largely dense that is a cooling and expanding state we ascribe to the hot big bang. This will take place for a universe inflating from a lower entropy state to an entropy state at higher level in the panorama of the hot big bang, where the entropy would carry on to get larger as it happens in our observed universe. The point of time in the earliest where the universe can be viewed approximately and hardly as classical is known as the Planck Era. It is thought that prior to this era the universe might be described as the hitherto unsuspected theory of certain quantum nature like quantum gravity etc. This era corresponds to E P 10 19 GeV > E > E G U T 10 15 GeV and the energies, temperature and times of particles are E P 10 19 GeV , T P 10 32 K , t P 10 43 s , respectively. Grand unified theories describe that at high energies as described above the Electroweak and strong force are unified into a single force and that these interactions bring the particles present into thermal equilibrium Electroweak Era corresponds to phase transitions that occur through spontaneous symmetry breaking (SSB) which can be characterized by the acquisition of certain non-zero values by scalar parameters known as Higgs fields. Until the Higgs field has zero values, symmetry remains observable and spontaneously breaks at the moment the Higgs field becomes non-zero. The idea of phase transitions in the very early universe suggests the existence of the scalar fields and provides the motivation for considering their effect on the expansion of the universe.
The power spectrum of CMB is calculated by measuring the magnitude of temperature variations versus the angular size of hot and cold spots. To understand the nature of CMB the spectrum of a perfect blackbody is given in Appendix C. During these measurements, a series of peaks with different strengths and frequencies are determined which conforms to the predictions of inflation theory which confirms that all sound waves were indeed produced at the same moment by inflation. It is believed that inflation might have given rise to sound waves-the waves traveling in the primordial vacuum-like medium with different frequencies after the big bang at 10 35 s starting in phase and would have been oscillating in radiation era for 380,000 years. Now, in the acoustic oscillations of the early universe, these must be measurable as power spectrum similar to that of measuring the sound spectrum of a musical instrument. The history of evolution of the universe to the present epoch is sketched in Figure 15, Figure 16 and Figure 17 show how the inflationary period is driven by the inflaton field:

10.1. Starobinsky R 2 -Inflation

Alexei Starobinsky proposed a cosmological inflationary phase of the universe shortly before Alan Guth in 1980 working in the framework of general relativity. The model is founded on the semiclassical Einstein field equations which provide a self-consistent solution for an exponentially accelerating era [49]. Starobinsky modified the general relativity to describe the behavior of very early universe undergoing an exponential period by suggesting quantum corrections to the energy momentum tensor T μ ν . The quantum corrections are calculated by taking the expectation value of the energy momentum tensor. Beginning with Einstein equations
G μ ν = 8 π T μ ν
R μ ν 1 2 g μ ν R = 8 π < T μ ν >
where < T μ ν > represents the expectation value of the energy momentum tensor. The expectation value of energy momentum tensor is the probabilistic value of a result or measurement which is fundamentally rooted in all quantum mechanical systems. Intuitively, it is the arithmetic mean of a large number of independent values of a variable under consideration. The energy momentum tensor T μ ν usually takes care of classical components of the universe in the form of matter and radiation in the context of flat spacetime as the parametric observations evidence in the recent data. In the case of curved spacetimes, nonetheless T μ ν might be vanishing gradually and < T μ ν > must be imparted contributions from quantum regime non-trivially. In the absence of classical components of the universe in the form of matter and radiation, the curvature of spacetime from quantum fluctuations of matter fields contribute to < T μ ν > non-trivially which Starobinsky utilized. These are known as quantum corrections to the energy-momentum tensor T μ ν . The quantum fluctuations of matter fields give non-trivial contributions to the expectation value of the energy momentum tensor < T μ ν > in the presence of cosmologically curved spacetime, regardless, matter and radian do not exist in classical style. In the background we consider FLRW spacetime
d s 2 = d t 2 + a ( t ) 2 1 k r 2 1 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
The spatial part 1 1 k r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2 of the metric represents the three geometries depending on the values of k. For k = + 1 , it represents a spherical geometry of 3-sphere which is finite, closed, and without boundary. For k = 0 , it represents a flat Euclidean geometry of 3-planes which is, in principle, infinite in extent, open, and without boundary. k = 1 , it represents a hyperbolic geometry of 3-hyperboloids which is infinite, open and without boundary. In the presence of conformally-invariant, free and massless fields, the quantum corrections adapt a simple form such that we can describe the expectation value of energy momentum tensor as
< T μ ν > = k 1 H μ ν 1 + k 2 H μ ν 2
where k 1 and k 2 are numerical coefficients in standard notation. In order to find < T μ ν > we have to compute these constants k 1 and k 2 and H μ ν 1 and H μ ν 2 . The coefficient k 1 is determined experimentally and can assume any value. The H μ ν 1 is a tensor and is conserved identically when expressed as the action given below and varied with respect to metric tensor g , i.e.,
H μ ν 1 = 2 g δ δ g μ ν d 4 x g R 2 = 0
and
H μ ν 2 = 2 R , μ , ν 2 g μ ν R , λ , λ + 2 R R μ ν 1 2 g μ ν R 2
The coefficient k 2 of H μ ν is defined uniquely in the following form
k 2 = N 0 + 11 2 N 1 / 2 + 31 N 1 1 1440 π 2
where N 0 , N 1 / 2 and N 1 denote the number of quantum fields with the subscripts of all three N’s 0, 1 / 2 and 1 representing spins of zero, half, and one respectively. In certain GUT theories due to larger multiplier factor of N 1 , the value of k 2 is majorly contributed by vector fields. Now, H μ ν is also a tensor and it does not conserve generally but conserves only in those spacetimes which are conformally flat like FLRW spacetimes in particular and cannot be obtained by varying a local action as in the case of H μ ν . The Equation (292) multiplying with 8 π G to both sides can be written as
8 π G < T μ ν > = 8 π G k 1 H μ ν 1 + 8 π G k 2 H μ ν 2
or
8 π G < T μ ν > = 48 π G 6 k 1 H μ ν 1 + 8 π G k 2 H μ ν 2
Now, we introduce the following parameters for convenience
M = 1 48 π G k 1 H 0 = 1 8 π G k 2
where both the parameters are positive i.e., H 0 > 0 and M > 0 . Now Equation (297) takes the form
8 π G < T μ ν > = 1 6 M 2 H μ ν 1 + H 0 2 H μ ν 2
Equation (299) can serve as the reasonable approximation in case of certain GUT models for the limit R > μ 2 , where μ represents the unified energy scale. Conformally invariant field equations usually describe the spinor and massless vector fields and contribute to < T μ ν > in the form of Equation (299). Further, if the number of matter fields is sufficiently bigger, then the corrections to Einstein’s field equations due to gravitons can also be ignored.

10.2. Trace Anomaly

The trace of expectation value of energy-momentum tensor < T μ ν > does not vanish rather it has a non-zero anomalous trace and this is what we call as trace anomaly. It is, however, interesting to note here that the trace of energy-momentum tensor without expectation value, i.e., T μ ν , vanishes for all those classical fields which are conformally invariant. Therefore, the trace of < T μ ν > is given by
< T ν ν > = M p l 2 H 0 2 1 3 R 2 R ν σ R ν σ M 2 R ; ν ; ν
The masses of the fields can be looked over in the limit of higher curvature, i.e., when R m 2 and in the same limit it remains true for the case of asymptotically free gauge theories where interactions between the fields become negligible. In de sitter space we can have
R μ ν = 1 4 R g μ ν
where R is constant curvature term i.e., Ricci scalar. Substituting now Equations (300) and (301) in Equation (290), we have R = 12 H 0 2 for non-trivial solution and the corresponding de Sitter solutions come about for k = 0 , + 1 , 1 , respectively,
a ( t ) = a 0 e t H 0
a ( t ) = 1 H 0 cosh t H 0
a ( t ) = 1 H 0 sinh t H 0
Equation (302) corresponds to k = 0 and gives a flat universe, Equation (303) gives a closed solution for k = + 1 , and the 3rd Equation (304) for k = 1 propounds the open de Sitter model of the universe. These solutions are impelled completely by the quantum corrections rendered to classical EFE and serve the purpose of inflationary epoch in the very early universe. Starobinsky inflation corresponds to a potential parameterized in terms of scalar field ϕ is V ϕ = 3 4 1 e 2 3 ϕ 2 .

10.3. Inflation and de Sitter Universe

In a very shorter period of time about 10 35 after spacetime came into being, the inflationary era of accelerating superluminal expansion known to be de Sitter phase took place. de Sitter phase removed all the wrinkles of curvature and warpage of spacetime so that the universe is to be observed flat. It further smoothed out all energy density stuff for the distribution of radiation and matter. One significant remnant as the traces of this fast expansion remains there known later on to be cosmic background radiation. In de Sitter universe there exists no ordinary matter, however, de Sitter retained cosmological constant which represents vacuum energy smeared out into the structure of spacetime. We can define the energy density of this non-relativistic matter
ρ Λ = Λ 8 π G
As p Λ = ρ Λ which gives an exotic form of matter with negative pressure, that is where the scale factor a ( t ) goes on increasing but a ˙ ( t ) is decreasing. We write
a ˙ a 2 = 8 π G 3 ρ Λ k a 2
a ¨ a = 4 π G 3 ρ Λ + 3 p Λ
and
ρ ˙ Λ + 3 H ρ Λ + p Λ = 0
From Equation (308) with ρ ˙ Λ = d ρ Λ d t = 0 and p Λ ρ Λ = w
3 H ρ Λ + p Λ = 0
or
3 H ρ Λ 1 + w = 0
Now from Equation (307) for p Λ ρ Λ = w
a ¨ a = 4 π G 3 ρ Λ 1 + 3 w
For w = 1 , Equation (311) becomes
a ¨ a = 4 π G 3 ρ Λ 1 3 = 8 π G 3 ρ Λ
or
d 2 a d t 2 8 π G 3 ρ Λ a = 0
Equation (313) is the equation of a harmonic oscillator. From Equation (306) for vanishing curvature, i.e., k = 0 where Λ dominates and a ˙ a = H
H Λ 2 = 8 π G 3 ρ Λ
Finding the value of ρ Λ
ρ Λ = 3 H Λ 2 8 π G
Substituting Equation (315) in Equation (313), and simplifying we have
d 2 a d t 2 8 π G 3 3 H Λ 2 8 π G a = 0
or
d 2 a d t 2 3 H Λ 2 a = 0
We can write the solution of above Equation (317) as
a ( t ) = C 1 exp H Λ t + C 2 exp H Λ t
differentiating Equation (318) twice with respect to time t
a ˙ ( t ) = C 1 H Λ exp H Λ t C 2 H Λ exp H Λ t
again
a ¨ ( t ) = C 1 H Λ 2 exp H Λ t + C 2 H Λ 2 exp H Λ t
using Equation (318) in Equation (320), we can write
a ¨ ( t ) = H Λ 2 C 1 exp H Λ t + C 2 exp H Λ t = H Λ 2 a ( t )
substituting the value of ρ Λ from Equation (315) in Equation (306), we have
a ˙ a 2 = 8 π G 3 3 H Λ 2 8 π G k a 2 = H Λ 2 k a 2
simplifying Equation (322) gives
k = H Λ 2 a 2 a ˙ 2
substituting the values of a ( t ) and a ˙ ( t ) from Equations (318) and (319) in above Equation (323)
k = H Λ 2 C 1 exp H Λ t + C 2 exp H Λ t 2 C 1 H Λ exp H Λ t C 2 H Λ exp H Λ t 2
simplification gives
k = 4 H Λ 2 C 1 C 2
Equation (325) means that the curvature term k depends upon the constants of integration C 1 and C 2 . For flat universe either C 1 = 0 or C 2 = 0 . The solution in Equation (318) becomes accordingly
a ( t ) = C 2 e H Λ t
and
a ( t ) = C 1 e H Λ t
Further Einstein equations are given by
G μ ν + Λ g μ ν = 8 π T μ ν
where
G μ ν = R μ ν 1 2 g μ ν R
and the form of solution of these equations upon which big bang standard cosmology is based, as worked out by Alexander Friedman (1922), George Lemaitre (1927), and afterwards by Robertson and Walker (1935) independently on the base of cosmological principle which put to use the homogeneity and isotropy, is
d s 2 = d t 2 + a ( t ) 2 1 k r 2 1 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2
where Ω 2 = d θ 2 + sin 2 θ d ϕ 2 . The metric in Equation (330) is characterized by scale factor a ( t ) and the curvature of spacetime k which are obviously determined by the self-gravitation of all the matter–energy content in the universe. We have incorporated dark matter and dark energy in the matter–energy content because their role is not avoidable at all in accelerated expansion and the present Minkowskian flat geometry of the universe. The solution of this line element gives Friedman equations using Einstein field equations that govern the time evolution of the universe and are given as
a ˙ 2 a 2 + k a 2 = 8 π G 3 ρ + Λ 3
and
a ¨ a = 4 π G 3 ρ + 3 p + Λ 3
The presence of cosmological term Λ in the above equations would be equivalent to that of a fluid having an equation of state p = ρ which is satisfied by
ρ + 3 p > 0
Looking at the things classically, we may approach the classical period of exponential expansion by using the first Friedmann equation by vanishing density ρ of radiation and baryons and the entailing curvature k in Λ -dominated Era which corresponds to equivalently having a fluid with p = ρ , thus Equation (331) becomes
a ˙ 2 a 2 + ( 0 ) a 2 = 8 π G 3 ( 0 ) + Λ 3 a ˙ 2 a 2 = Λ 3 a ˙ = d a d t = Λ 3 a d a a = Λ 3 d t
after integrating and simplifying, we get
a = e Λ 3 t
Equation (335) gives the exponential expansion of the scale factor. It describes the fact that when the universe was dominated by cosmological constant Λ , the rate expansion was much faster than the present day scenario. From Equation (332)
a ¨ = 4 π G 3 ρ + 3 p a + Λ 3 a
Considering a closed volume with energy U = ρ V = ρ 4 π 3 a 3 and now we see how inflationary period is obtained in the perspective of particle physics where a negative pressure is achieved for it to take place. Friedmann solved EFE with Λ = 0 , so
a ˙ 2 a 2 + k a 2 = 8 π G 3 ρ
a ¨ a = 4 π G 3 ρ + 3 p
Equation (338) is known as acceleration equation. The inflationary period, as its definition implies, is the acceleratingly expanding phase of the universe in a very small fraction of first second, as the expansion is characterized by the scale factor a; therefore, we have such an era as
a ¨ > 0
thus Inflationary era
a ¨ > 0
dividing both sides of Equation (317) by scale factor a
a ¨ a > 0
which is LHS of Equations (338) and (341) imposes the condition on RHS of Equation (338)
4 π G 3 ρ + 3 p > 0 ρ + 3 p < 0 p < 1 3 ρ ρ > 3 p
For the inflation to occur and set the universe in an accelerating phase, we require the matter to possess an equation of state with negative pressure. The possibility of this negative pressure p which is less than negative of one-third of density is in perspective of symmetry breaking in modern models of particle physics. From
a ˙ 2 a 2 + k a 2 = 8 π G 3 ρ
a ˙ 2 = 8 π G 3 ρ a 2 k
For a ¨ > o , the scale factor shall increase faster than a t t and the term 8 π G 3 ρ a 2 shall increase during this accelerated era such that the curvature term k will become negligibly small and shall vanish. Inflationary era is also defined by considering the shrinking of Hubble Sphere [43] due to its direct linkage to the horizon problem and because it provides a fundamental role in producing of quantum fluctuations. Shrinking Hubble Sphere is defined as
d a H 1 d t < 0
d a H 1 d t = d a a ˙ a 1 d t = d a ˙ 1 d t = a ¨ a 2
a ¨ a 2 < 0
which will imply accelerated expansion
a ¨ > 0
At t = 0 , the scale factor a characterizing expansion of the universe comes out to be of a specific value. In Equation (337), when ρ = ρ ϕ is of very larger value and the scale factor a dominates over the curvature term k, then we have
a ˙ a 2 = H 2 = 8 π G 3 ρ ϕ
a = a 0 e H t
de Sitter line element is given by
d s 2 = d t 2 + e 2 H t d x 2 + d y 2 + d z 2
inflation has to terminate and H is constant, meaning that the de Sitter phase cannot give perfect inflationary era, however for H ˙ H 2 , it would compensate. It would be interesting here to note that Z. G. Lie and Y.S. Piao have shown that the universe we observe today may have emerged from a de Sitter background without having the requirement of a large tunneling in potential and with low energy scale [50].

10.4. The Conditions under Which the Inflation Occurs

Shrinking Hubble sphere has been considered as basic definition of inflationary era due to its direct connection to the horizon problem and with mechanism of quantum fluctuation generations [51]. differentiating the comoving Hubble radius a H 1 with respect to time we find the acceleratedly expanding Hubble sphere
t a H 1 = a ¨ a ˙ 2
We see that a ¨ a ˙ 2 < 0 , multiplying the inequality by 1 and simplifying, we have
a ¨ > 0
which means that shrinking comoving Hubble sphere a H 1 points toward accelerated expansion a ¨ > 0 . As Hubble sphere H remains nearly constant, in order to understand the meaning of nearly constant we see how its slow roll variation takes place, so taking H as variable
t 1 a H = a ˙ H + a H ˙ a H 2 = 1 a 1 + H ˙ H 2
where H ˙ H 2 = ε known as slow roll parameter. It can be inferred that H ˙ H 2 < 0 implies shrinking Hubble sphere.

10.5. Slow Roll Inflation—The Dynamics of Scalar Field

Elementary particles in modern physics are represented by quantum fields and oscillations of these fields are translated as particles. Scalar fields represent spin zero particles in field theories and look like vacuum states because they have same quantum numbers as vacuum. The matter with negative pressure ρ = p represents physical vacuum-like state where the quantum fluctuations of all types of physical fields exist. These fluctuations can be considered as waves of all possible wavelengths related with physical fields, i.e., wavy physical fields moving freely in all directions. The negative pressure violates the strong energy condition which is necessary for the inflation to occur. To keep things simpler a single scalar field namely inflaton ϕ = ϕ x , t is considered present in the very early universe, as the value of the scalar field depends upon position x in space which assigns potential energy to each field value. It is also dynamical due to being function of time t and has kinetic energy as well, i.e., energy density ρ ϕ associated with the inflaton ϕ is ρ ϕ = ρ p + ρ k . The ratio of the potential and kinetic energy terms of ϕ = ϕ x , t , decides the evolution of the universe. The Lagrangian of the scalar inflaton field ϕ is expressed as the energy difference between its kinetic and potential terms, that is
L = 1 2 g μ ν 2 ϕ V ( ϕ )
It is assumed that the background of FLRW universe has been sourced by energy-momentum associated with the inflaton that dominates the universe in the beginning. We shall observe under what conditions this causes accelerated expansion of the FLRW universe.
S = d 4 x g L = d 4 x g 1 2 g μ ν 2 ϕ V ( ϕ )
The energy-momentum tensor of the inflaton field is given as
T μ ν = μ ϕ ν ϕ g μ ν L
T μ ν = μ ϕ ν ϕ g μ ν L
which for μ = 0 , ν = 0 results as
T 00 = 1 2 ϕ ˙ 2 + 1 2 a 2 2 ϕ + V ( ϕ )
and for μ = ν = j
T j j = 1 2 ϕ ˙ 2 1 6 a 2 2 ϕ V ( ϕ )
The gradient term vanishes, in the otherwise condition, the pressure gained is much less than the required value to impart impetus for inflation to take place, therefore we obtain the following values for energy density and pressure
ρ ϕ = T 00 = 1 2 ϕ ˙ 2 + V ( ϕ )
and
p ϕ = T j j = 1 2 ϕ ˙ 2 V ( ϕ )
The condition V ( ϕ ) > > ϕ ˙ 2 corresponds to the negative pressure condition ρ ϕ = p ϕ which means that the potential (vacuum) energy of the inflaton derives inflation. Now using Euler–Lagrange equations
μ δ g L δ μ ϕ δ g L δ ϕ = 0
we can find equation for inflaton field that comes to be
ϕ ¨ + 3 a ˙ a ϕ ˙ 1 a 2 ( t ) 2 ϕ + V , ϕ ( ϕ ) = 0
It can also be computed from the energy density and the pressure terms given in Equations (361) and (362) respectively by substituting in equation of energy conservation
d ρ d t + 3 H ρ + p = 0
Equation (365) in terms of inflation field ϕ
d ρ ϕ d t + 3 H ρ ϕ + p ϕ = 0
By substituting Equations (361) and (362) in Equation (366), we have
d 1 2 ϕ ˙ 2 + V ( ϕ ) d t + 3 H 1 2 ϕ ˙ 2 + V ( ϕ ) + 1 2 ϕ ˙ 2 V ( ϕ ) = 0
ϕ ¨ + V ( ϕ ) + 3 H ϕ ˙ ϕ ˙ = 0
ϕ ¨ + V ( ϕ ) + 3 H ϕ ˙ = 0
where V ( ϕ ) = d V ( ϕ ) d ϕ and the term 3 H ϕ ˙ is known as friction term and offers friction to the inflaton field when it rolls down ( ϕ ˙ ) its potential during expansion of the universe H = a ˙ a . Figure 18 and Figure 19 manifest how the scalar field drives the universe evolution in the beginning and how does it slow roll afterward respectively:

10.6. Conditions of the Slow Roll Inflation

According to the big bang model, that is, the currently accepted model, the universe is about 14 billion years old. At the point of existence the curvature of spacetime was very large or equivalently can be described in other words that space was largely warped and curved where only quantum effects can prevail and the question of time to exist is likely to become absurd. From this state how the very brief era of exponential expansion can be had is fulfilled by assumption of scalar field which take the responsibility of such state mentioned. We know from the 2nd Friedmann’s equation which is acceleration equation
a ¨ a = 4 π G 3 ρ ϕ + 3 p ϕ
For a ¨ > 0
ρ ϕ + 3 p ϕ < 0 p ϕ < 1 3 ρ ϕ
From Equations (361) and (362), substituting for p ϕ and ρ ϕ in Equation (371)
1 2 ϕ ˙ 2 V ( ϕ ) < 1 3 1 2 ϕ ˙ 2 + V ( ϕ )
solving the inequality and keeping in mind that ϕ ˙ is a squared term, we have
ϕ ˙ 2 < < V ( ϕ )
which means that the inflaton field is slowly rolling down its potential. Differentiating Equation (373) with respect to time, we have
ϕ ¨ < 1 2 V ( ϕ )
Now from Equation (369), we obtain
ϕ ¨ + V ( ϕ ) = 3 H ϕ ˙
We neglect the acceleration providing term ϕ ¨ = d 2 ϕ d t 2 as the inflaton field has to roll now slowly to escape from graceful exit problem in inflation i.e., it is decelerating, so we write
V ( ϕ ) = 3 H ϕ ˙
plugging Equation (376) in Equation (374)
ϕ ¨ < 1 2 ( 3 H ϕ ˙ )
On neglecting the constant factor, it gives
ϕ ¨ < < 3 H ϕ ˙
differentiating now Equation (376) with respect to time,
3 H ˙ ϕ ˙ + H ϕ ¨ = V ( ϕ ) ϕ ˙
As H remains constant during inflation, therefore H ˙ vanishes and we have
ϕ ¨ = V ( ϕ ) ϕ ˙ 3 H
Putting Equation (380) in Equation (378), we have
V ( ϕ ) ϕ ˙ 3 H < < 3 H ϕ ˙
It gives
V ( ϕ ) < < H 2

10.7. Parameters for the Slow Roll Inflation

Two slow roll parameters ε and η are defined in terms of Hubble parameter H as well as potential V which quantify slow roll inflation.
ε H = H ˙ H 2
Using the relation a ( t ) e N N = ln a , it can also be expressed in the form
ε H = d ln H d N
where N is the number of e-folds and 2nd is defined as
η H = 1 2 H ¨ H ˙ H
From 1st Friedmann equation
a ˙ a 2 k a 2 = 8 π G 3 ρ
For ρ = ρ ϕ and from Equation (361) ρ ϕ = 1 2 ϕ ˙ 2 + V ( ϕ ) , as during inflation V ϕ > > ϕ ˙ 2 , so that ρ ϕ = V ( ϕ ) also curvature term k is negligibly small, so that Equation (386) becomes
H 2 = 8 π G 3 V ϕ
differentiating Equation (387) with respect to time and simplifying
H ˙ = 4 π G 3 H V ϕ ϕ ˙
And from Equation (376) substituting in Equation (388), we have
H ˙ = 4 π G ϕ ˙ 2
Substituting above in Equation (383), we have
ε H = H ˙ H 2 = 4 π G ϕ ˙ 2 H 2 = 4 π G H 2 ϕ ˙ 2
Again from Equation (376), we have
ϕ ˙ = V ( ϕ ) 3 H
squaring
ϕ ˙ 2 = V 2 ( ϕ ) 9 H 2
substituting in Equation (390)
ε H = 4 π G H 2 V 2 ( ϕ ) 9 H 2 = 4 π G V 2 ( ϕ ) 9 H 2 2
From Equation (387) putting for H 2
ε V = 4 π G V 2 ( ϕ ) 9 8 π G 3 V ( ϕ ) 2 = 1 16 π G V ( ϕ ) V ( ϕ ) 2 = M p l 2 2 V ( ϕ ) V ( ϕ ) 2
η H can also be expressed as
η H = ϕ ¨ H ϕ ˙
η V = 1 8 π G V ( ϕ ) V ( ϕ ) = M p l 2 V ( ϕ ) V ( ϕ )
From Equation (387) H 2 = 8 π G 3 V ( ϕ ) , which gives 8 π G V ( ϕ ) = 3 H 2 substituting above in Equation (396), we have
η V = V ( ϕ ) 3 H 2

10.8. Number of e-Folds

It is usual practice to have the inflation quantified and the quantity which does this is called number of e-fold denoted by N before the inflation ends. As the time goes by N goes on decreasing and becomes zero when inflation ends. It is counted or measured backwards in time from the end of inflation which means that N = 0 at the end of inflation grows to maximal value towards the beginning of inflation. It measures the number of times the space grows during inflationary period. The amount of e-folds necessarily required to resolve the big bang problems of Horizon, Flatness, Monopole, Entropy, etc. is N 60–75 depending upon the different models and on the reasonable estimation of the observational parameters. To find the number of e-folds between beginning and end of inflation we know that during inflation the scale factor evolves as
a t = a t 0 e H t
or
a t = a t 0 e H t t i
The factor H t constitute the number of e-folds denoted by N, i.e.,
N = H t
differentiating Equation (400) with respect to time
d N d t = H = t ln a
N = t i t f H d t = t i t f a ˙ a d t = ln a t f a t i
Further, the relation between Hubble parameter H and the number of e-folds N can be written. we have derived earlier the evolution equation for inflaton field that reads
ϕ ¨ + 3 H ϕ ˙ + V , ϕ = 0
During slow roll inflation ϕ ¨ = 0 , so that Equation (403) becomes
3 H ϕ ˙ + V , ϕ = 0
3 H ϕ ˙ = V , ϕ
Moreover, during slow roll the Friedmann’s 1st equation evolves as with k = 0 and ρ = V ( ϕ ) + 1 2 ϕ ˙ 2
H 2 = 8 π G 3 V ( ϕ ) + 1 2 ϕ ˙ 2
During slow roll ϕ ˙ 2 < < V ϕ and only ϕ ˙ works, thus Equation (406) becomes
H 2 = 8 π G 3 V ( ϕ )
Dividing Equation (405) by Equation (407)
ϕ ˙ H = V , ϕ 8 π G V ϕ
Now from Equation (400), we can write because t = t f t i , so t = t i t f d t and with dividing and multiplying by d ϕ
N = H t = t i t f H d t = t i t f H d t d ϕ d ϕ
where ϕ ˙ = d ϕ d t , Equation (409) takes the form
N = ϕ i ϕ f H ϕ ˙ d ϕ
Substituting from Equation (408) after inverting
N = ϕ i ϕ f 8 π G V ϕ V , ϕ d ϕ = 8 π G ϕ i ϕ f V ϕ V , ϕ d ϕ
or
N = 8 π G ϕ f ϕ i V ϕ V , ϕ d ϕ
Thus number of e-folds can be found in terms of potential of the inflaton field. Further slow roll parameter ε H can be described in terms of number of e-fold N, we know
ε H = H ˙ H 2 = 1 H 2 d H d t = 1 H 2 d H d N d N d t
ε H = 1 H 2 d ln N d t

11. Inflationary Solutions to the Big Bang Problems

Horizon, flatness, entropy, and monopole problems are initial value problems which inflation solves in one go. Inflation explains why the observable universe is spatially flat, isotropically homogeneous and so large in size.

11.1. Inflation and Horizon Problem

We consider that the inflation begins at a time t i and comes to an end at some time t f and the expansion rate H = t ln a , curvature term k and energy density of matter and radiation ρ = ρ M + ρ R during inflation vanishes, we know that
a ( t ) = a ( t ) e H t = a ( t ) e H ( t t i )
and
a ( t ) = a ( t ) e H t = a ( t ) e H ( t f t )
We will find how long the inflation must sustain to resolve the horizon problem. We can find the corresponding e-folding number N that is
N = H t = t i t f H d t
As H = a ˙ a ,
N = t i t f a ˙ a d t = t i t f d a a
N = ln a t i t f = ln a f a i
or
a f a i = e N
or
a i a f = e N
Now the horizon scale observed today H 0 1 was reduced during inflation to a value of λ H 0 ( t i ) which is smaller than the horizon length during inflation.
λ H 0 ( t i ) = R H 0 a t i a t 0
Dividing and multiplying Equation (422) by a t f
λ H 0 ( t i ) = R H 0 a t i a t 0 × a t f a t f
Now, from Equation (421) a i a f = e N and using the relation between scale factor and temperature during this phase a 1 T a i 1 T i and a f 1 T f so that we have
λ H 0 ( t i ) = H 0 1 T 0 T f e N
where R H 0 = H 0 1 . Now λ H 0 ( t i ) < H I 1 where H I 1 is the horizon length during inflation. Therefore, Equation (424) can be expressed as
H 0 1 T 0 T f e N H I 1
H 0 1 H I 1 T 0 T f e N
or
e N T 0 H 0 T f H I N ln T 0 H 0 T f H I
N ln T 0 H 0 ln T f H I
or
N 67 + ln H I T f
N 70

11.2. Inflation and Flatness Problem

From 1st Friedmann equation,
H 2 + k a 2 = 8 π G 3 ρ
We found the density parameter expression
Ω 1 = k H 2 a 2 ( t )
during the inflationary period Hubble parameter giving expansion rate remains almost constant so that Equation (432) is
Ω 1 = k H 2 a 2 ( t ) 1 a 2 ( t )
We observed earlier that
Ω 1 t = t p l 10 60
which means that to have the value of the density parameter as observed today, i.e., Ω 0 to be of the order of unity, the initial value of Ω at the beginning of the radiation-dominated era must be same as given in Equation (434) above, and from Equation (432), we can write for the time at the beginning of inflationary era
Ω 1 t = t i = k H I 2 a i 2 ( t )
and for the time when inflationary period comes to end
Ω 1 t = t f = k H I 2 a f 2 ( t )
Further, the beginning of the radiation-dominated era can be recognized with the beginning of inflationary phase such that it is required
Ω 1 t = t f = 10 60
dividing now Equation (436) by Equation (435)
Ω 1 t = t f Ω 1 t = t i = k H I 2 a f 2 ( t ) k H I 2 a i 2 ( t ) = a i 2 ( t ) a f 2 ( t )
We calculated a i a f = e N , the above Equation (438) takes the form
Ω 1 t = t f Ω 1 t = t i = e N 2
e 2 N = Ω 1 t = t f Ω 1 t = t i
With taking Ω 1 t = t i 1 ,
N 1 2 ln Ω 1 t = t f
N 70

11.3. Inflation and Entropy Problem

The entropy problem can be resolved if a large amount of entropy is created in the very early universe non-adiabatically [48,51] which is accomplished by inflationary era in a finite time in the early history of the universe. Let the entropy at the end of inflation is S f and in the beginning it was S i such that S f S i , then
S f = M 3 S i
where M is the numerical factor with value M 3 = 10 10 M = 10 30 . Now S f = S U . We know that S a T 3 , so that we can write for
S i a i T i 3
and for
S f a f T f 3
where T i and T f are the measures of temperature at the beginning and end of the inflationary period. dividing Equation (445) by Equation (444) we have
S f S i a f a i 3 T f T i 3
or
a f a i 3 S f S i T i T f 3
a f a i S f S i 3 T i T f
Now a f a i = e N , and considering that at the beginning of inflationary phase the total entropy of the universe was of the order 1, i.e., S i 1 and S f = S U , thus Equation (448) takes the form
e N S U 3 T i T f
N ln S U 3 T i T f
N 70
therefore, entropy problem is resolved by inflationary period.

11.4. Inflation and Monopole Problem

In grand unified theories (GUT), the standard model S U 3 × S U 2 × U 1 in particle physics emerges out of a simple symmetry group breaking. In these theories, heavy particles of very high density are predicted to be created which are known as magnetic monopoles. The cosmological monopoles prior to the period of inflation are considered to take place and are allowed supposedly to exist. It means that inflation allows the existence of magnetic monopoles, that is to say that as they are created earlier to era of inflation, these magnetic monopoles are supposed to form in symmetry breaking during phase transitions and where the inflationary era is considered to take place just after it. Inflation dilutes the density of these magnetic monopoles n m p N m p a 3 0 to the negligibly small size such that these become so small to be detected today [52]. During inflation monopoles collapse in an exponential way and their abundant presence falls to the level of being hardly detectable.

11.5. Inflation and Observations

Cosmological perturbations are an important relic of the inflation used to describe the anisotropies of the cosmic microwave background (CMB) and structure evolution and formation of the universe. The seeds of inhomogeneities, which represent all the structure in the universe were produced during inflationary phase and were stretched to the astronomical scales with the exponential expansion. These inhomogeneities are what we see today as stars, galaxies, etc. in the form of baryonic matter. From the theory of linear perturbations and from the relation δ ϕ δ g μ ν we know how to categorize FLRW metric perturbations at the first order in the form of scalar, vector, and tensor perturbations of spins 0, 1, and 2, respectively. A very important parameter which determines the properties of the perturbations of the scalar field is power spectrum p ϕ k . The scalar field power spectrum at the time of horizon crossing comes out to be
p ϕ k = H 2 π k = a H 2
and the curvature power spectrum is calculated to be
p R k = H ϕ ˙ H 2 π k = a H 2
and the power spectrum for tensor perturbations is given by
p T k = 64 π G H 2 π k = a H 2
The tensor-to-scalar ratio of the spectrum is defined in the following expression and is found in terms of the first slow roll parameter:
r = p T k p R k = 16 ε V
We can now define scalar and tensor spectral indices as
n s 1 = d ln p R k d ln k
and
n T = d ln p T k d ln k
scale invariance of the scalar power spectrum is characterized by n s 1 = 0 n s = 1 . Deviations from the scale invariance in an inflationary model gives it specific features. In this case, the spectral indices to their lowest order can be described in terms of potential slow roll parameters ε V and η V . In case of large field model, a general polynomial potential is
V ϕ = Λ 4 ϕ μ p
and we find
r = 4 p N
n s 1 = 1 N 1 + p 2
For a particular inflationary model p must be assigned a value greater than unity [53].

11.6. Λ CDM

The standard model of cosmology describes a universe that evolves from a singularity at t = 0 . This singularity is known as the big bang and marks the instant when the universe begins its evolution in time. A detailed discussion of big bang theory of creation is presented in Appendix D. The kinematics of it is described by FLRW spacetime, and its dynamics is governed by the Friedmann equations in the framework of general relativity. The standard model is usually known as big bang model [54] due to extrapolation of redshift towards big bang singularity. The observational parameters are not fixed by the standard model or big bang which means that the big bang is parameterizable. Λ CDM model constitutes one of such parameterizations and shows remarkable consistency with the recent observations. That is why it has gathered support by majority of cosmologists. It incorporates the ingredients, namely, cosmological constant Λ and cold matter (CDM), in addition to ordinary matter. It is interesting to note that the nature of both ingredients contained in it are colossally unknown for which theoretical and observational developments are underway. It is also called the Concordance model for being it in agreement with the recent measurements of parameters. Obviously the dynamics of Λ CDM is governed by general relativity. Λ was introduced by Einstein himself [3] to balance the gravitational effect of the ordinary matter in order to show a static model of the universe. Therefore, the energy density of Λ is contained in the structure of spacetime itself or in other words it is vacuum energy of space. However, it was dropped by Einstein after the expansion was confirmed in 1929 calling it as the biggest blunder of his life ever made. After the accelerated expansion was discovered in 1998, Λ is coming back once again to accommodate the effect of accelerated expansion, but this time it is expected to have a dynamic nature. According to the recent observations [35,36], dark energy is ~ 70 % , dark matter constitutes 25 % , and the ordinary matter (baryonic) is 5 % only. In the framework of the Λ CDM model, the nature of dark energy presents one of the most challenging issues to the present day cosmology. In Λ CDM, space is spatially flat and the radius of curvature, therefore becomes infinitely large. The Λ CDM adapts a minimum number of parameters to describe the universe that is six. The recent measurements of these parameters from different sources are given [35,55] below in the Table 1.
From these six free parameters, we can deduce other parameters like Hubble constant with some assumption about the cosmological model. The detail can be found in [35]. In the context of Λ CDM we can find how the energy densities is related with the parameters like scale factor a t , time t, Hubble parameter, etc. in the corresponding sectors namely radiation, matter (cold matter), and dark energy. Using Equations (96) and (108), we find for the radiation sector for which ρ = 3 p w = 1 3
ρ r a = ρ r 0 a 4 t = 3 M p l 2 H 0 2 Ω r 0 a 4 t
Ω r a = H 0 H 2 Ω r 0 a 4 t
and
a t t 1 2
for the matter sector incorporating baryonic and cold dark matter for which ρ = 3 p w = 1 3 we have
ρ m a = ρ m 0 a 3 t = 3 M p l 2 H 0 2 Ω m 0 a 3 t
Ω m a = H 0 H 2 Ω m 0 a 3 t
and
a t t 2 3
for the dark energy sector for which ρ Λ = p w Λ = 1 , we have
ρ Λ = c o n s t a n t
and
a t exp Λ 3 t
Using Friedmann equation and total energy density of all forms, we can determine the background dynamics in the form of two very significant parameters
H 2 a = H 0 2 Ω m 0 a 3 + Ω r 0 a 4 + Ω Λ 0 a 0
q z = 1 2 Ω m z + 2 Ω r z Ω Λ z
The observations of Planck collaboration on cosmic microwave background radiation (CMB) give [35] the values of these parameters to be Ω m 0 = 0.3089 , Ω r 0 = 5.38916 × 10 5 and Ω Λ 0 = 0.691046 and for a flat Λ CDM model, the observations from Type SNe Ia supernovae give [3] Ω m 0 = 0.295 and for other we can estimate from Ω Λ = 1 Ω D M Ω b . The Figure 20 gives the dark energy and cold matter densities [56] in terms of density parameters Ω Λ and Ω m .

11.7. Inflation and Dark Energy in f R Modified Gravity

Einstein Field Equation (EFE) of general relativity is well known
G μ ν = R μ ν 1 2 g μ ν R = 8 π T μ ν
Equation (471) corresponds to Einstein-Hilbert action
S E H = 1 16 π G d 4 x g R + d 4 x g L m
In scalar field models we usually modify RHS, i.e., energy-momentum tensor (matter sector) and accordingly add some terms for a scalar field. If the RHS is kept unaltered and LHS is modified that stands for the geometry of spacetime mimicking the role of gravity. Due to this reason, it is called the model of modified gravity. The LHS of EFE is derived merely from the curvature term, i.e., Ricci Scalar R, however in the modified gravity we replace it by a general function of it [57,58,59,60]. Replacing the Ricci scalar R in the Einstein–Hilbert action given in Equation (471) by a general function of R, that is, f R , i.e., R f R , we have
S E H f R = 1 16 π G d 4 x g f R + d 4 x g L M
The variation of Equation (473) would be
δ S E H f R = d 4 x δ g f R + d 4 x δ L
Equation (474) yields through tedious calculations the following modified gravity equation
F R R μ ν 1 2 f R g μ ν μ ν F R + g μ ν F R = k T μ ν
where = ν ν and F R = f , R R In Equation (474), the LHS is the modified form of R μ ν 1 2 g μ ν R = G μ ν . We contract Equation (475) with g μ ν to find the trace of modified EFE
F R R + 3 F R 2 f R = k T
For a vacuum solution T = 0 , the de Sitter space with curvature term R to be constant
F R R + 3 F R 2 f R = 0
Equation (477) represents an inflationary solution with the term 3 F R = 0 .
F R R 2 f R = 0
If the condition in Equation (478) is fulfilled, the late time de Sitter solution can be obtained in a f R -based dark energy model. The Friedmann Equations for the modified gravity can be determined in the following way by a spatially flat expanding FLRW universe
d s 2 = d t 2 + a 2 t η μ ν d x μ d x ν
now going through the lengthy calculations, we find
R 00 = 3 H ˙ + H 2
R j j = 2 a ˙ 2 + a a ¨
and
R = 6 H ˙ + 2 H 2
Equation (480) can be further re-expressed using Equation (482)
R 00 = 1 2 R + 3 H 2
for
T μ ν = d i a g ρ M , p M , p M , p M
where the trace of T μ ν is ρ M + 3 p M . Now μ = ν = 0 in Equation (475) gives
F R R 00 1 2 f R g 00 0 0 F R + g 00 F R = k T 00
solving through tedious calculations by making use of Equation (480) and Equation (484), we reach at the result
3 H 2 F = 1 2 F R f 3 H F ˙ + k ρ M
again for μ = ν = j in Equation (475), we have
F R R j j 1 2 f R g j j j j F R + g j j F R = k 2 T j j
we find
g j j F R = 2 a 2 H F ˙ a 2 F ¨ + j j F
Using Equations (481), (484) and (488) in Equation (487), we obtain
2 H ˙ F = F ¨ + H F ˙ k ρ M + p M
Equations (486) and (489) together determine the background dynamics of a flat FLRW universe governed by f R . From Equation (486), dividing by 3 H 2 F and for ρ = ρ R + ρ M , we can construct a dynamical system in the framework of f R , that is,
F ˙ H F + f 6 H 2 F + R 6 H 2 + k ρ R 3 H 2 F + k ρ M 3 H 2 F = 1
here we can define the following parameters
x 1 = F ˙ H F
x 2 = f 6 H 2 F
x 3 = R 6 H 2
x 4 = Ω R = k ρ R 3 H 2 F
x 5 = Ω M = k ρ M 3 H 2 F
Equation (493) can be recast by using Equation (482)
x 3 = H ˙ H 2 + 2 = R 6 H 2
we have
x 1 + x 2 + x 3 + x 4 + x 5 = 1
Now for p M = 4 3 ρ R , N = ln a and H ˙ = H H , we can determine the following dynamical system:
d x 1 d N = x 1 2 x 1 x 3 3 x 2 x 3 + x 4 1
d x 2 d N = x 1 x 2 2 x 2 x 3 2 + x 1 m
where
m = R f , R R f , R
d x 3 d N = 2 x 3 x 3 2 x 1 x 3 m
d x 4 d N = x 4 x 1 2 x 3
d x 5 d N = x 5 x 1 2 x 3
The effective equation of state can also be written form Equations (486) and (489) by division
w e f f = 1 2 H ˙ F 3 H 2 F
Or using Equation (496)
w e f f = 1 3 1 + 2 x 3
Another form can also be written as
w e f f = 1 + F ¨ H F ˙ + k ρ + p 1 2 F R f 3 H F ˙ + k ρ M
Now we find f R inflation by considering first a general form of f R and determine its dynamics. Afterwards, Starobinsky inflation in f R will be discussed. Let us consider
f R = R + b R n
we have
F = f , R R = 1 + n b R n 1
and
F ˙ = F t = n n 1 b R n 2 R t
and
F R = R + n b R n
substituting Equations (507)–(510) in Equation (486) where the k ρ M vanishes during inflationary phase, we have
3 H 2 1 + n b R n 1 = 1 2 n b 1 R n 3 n n 1 H b R n 2 R t
The cosmological acceleration could be realized in the regime F > > 1 1 + n b R n 1 > > 1 and n b R n 1 > > 0 , which implies 1 + n b R n 1 n b R n 1 . Dividing Equation (511) by 3 n b R n 1 , we obtain after simplification
H 2 = 1 6 n 1 n R 6 H n R R t
also taking the time derivative of Equation (482)
R ˙ = 6 H ¨ + 4 H H ˙
finding the ratio
1 R R t = H H ˙ H 2 4 + H ¨ H H ˙ 2 + H ˙ H 2
In Equation (514), the following approximations are validated during the inflation era
H ˙ H 2 < < 1 H ˙ H 2 0
and
H ¨ H H ˙ < < 1 H ¨ H H ˙ 0
Equation (514) takes the form
1 R R t = 2 H ˙ H
Using Equation (482) and Equation (517) in Equation (512), we get
H 2 = n 1 n H 2 2 + 1 2 n H ˙ H 2
Or
H ˙ H 2 = 2 n n 1 2 n 1
Let
2 n n 1 2 n 1 = ε
then Equation (519) becomes
H ˙ H 2 = ε
Equation (519) gives the slow roll parameter. Integrating Equation (521), we get
H = 1 ε t
Equation (522) gives the expansion rate. We can determine the evolution of the scale factor a t
a t = t 1 ε
Equation (523) implies that the inflationary period is the one during which cosmic expansion is very fast in a very short span of time. Now to find out the condition on n, we use the condition ε < 1 . From Equation (520) we have
2 n n 1 2 n 1 < 1
which gives
n > 1 2 1 ± 3
We observe here that for n = 2 , the slow roll parameter ε = 0 , which means the Hubble parameter H remains nearly constant in the limit F > > 1 . For n > 2 it implies super-inflation with H ˙ > 0 . It means that the viable condition to achieve a phase of standard inflation must be
1 2 1 + 3 < n < 2
For the Starobinsky R 2 -inflation, we put n = 2 and b = 1 6 M 2 in Equation (507)
f R = R + 1 6 M 2 R 2
using Equations (482) and (513), after substituting for n and b in Equation (511), we get after simplification
H ¨ 1 2 H ˙ 2 H + 1 2 M 2 H + 3 H H ˙ = 0
as the Hubble parameter remains almost constant during inflation or varies very slowly, therefore the first two terms H ¨ and 1 2 H ˙ 2 H can be neglected, and we get
H ˙ = 1 6 M 2
Integration of above Equation (529) between the limits H i , H f and t i , t f gives
H f H i = 1 6 M 2 t f t i
For H f = H = t a , we have on integrating again
a t = a i t exp H i t t i 1 12 t t i 2 M 2
Equation (482) becomes on making use of Equation (529)
R = M 2 + 12 H 2
The slow roll parameter on using Equation (529) takes the form
ε = M 2 6 H 2 < 1
for ε 1 , we have from Equation (533)
H 2 1 6 M 2 6 H 2 M 2 H 2 M 2
for t = t f ε 1 , from Equations (521) and (529) for H = H f we obtain
H f = 1 6 M
substituting Equation (535) in Equation (530) we have
t f t i = 6 H i M 2
where the term 1 6 M 2 is ignored. In addition, Equation (532) reduces to R M 2 when the inflationary phase approaches towards its ending. The number of e-folds for the inflationary phase turns out to be
N = t i t f H d t = 1 2 ε

12. The Line Element of the Perturbed Universe

When the perturbations of inflaton field ϕ are considered the energy momentum tensor T μ ν is also perturbed. The perturbations of the inflaton field are thus consequently reflected through the metric tensor g μ ν such that δ ϕ δ g μ ν . We discuss here only the scalar perturbations for which the metric takes the form as
d s 2 = a 2 ( t ) 1 2 A d τ 2 + 2 i B d τ d x i + 1 2 ψ δ i j + D i j E ) d x i d x j
δ g μ ν = a 2 ( t ) δ g 00 δ g 0 i δ g i 0 δ g i j = a 2 ( t ) 1 2 A i B i B 1 2 ψ δ i j + D i j E

12.1. Inverse of δ g μ ν

Let the inverse of δ g μ ν be δ g μ ν and we write
δ g μ ν = 1 a 2 ( t ) δ g 00 δ g 0 i δ g i 0 δ g i j = 1 a 2 ( t ) 1 + X 1 a 2 ( t ) i Y 1 a 2 ( t ) i Y 1 a 2 ( t ) 1 + 2 Z δ i j + D i j K
We find the inverse of Equation (454) that is g μ ν , so that we have
g μ ζ g ζ ν = g 0 μ ζ g ζ ν 0 = δ ν μ g 0 μ ζ + δ g μ ζ g ζ ν 0 + δ g ζ ν = δ ν μ g 0 μ ζ g ζ ν 0 + g 0 μ ζ δ g ζ ν + δ g μ ζ g ζ ν 0 + δ g μ ζ δ g ζ ν = δ ν μ
where δ ν μ is the Kronecker delta function and is defined as
δ ν μ = 1 , if μ = ν 0 , if μ ν
and g ζ ν is the simply unperturbed FLRW line element described as
d s 2 = a 2 ( t ) d t 2 + δ i j d x i d x j
g μ ν 0 = g μ ν = a 2 ( t ) g 00 g 0 i g i 0 g i j = a 2 ( t ) 1 0 0 δ i j
The inverse of g μ ν 0 = g μ ν is simply g 0 μ ν = g μ ν since for diagonal unperturbed metric g μ ν = 1 g μ ν , so that we can write
g 0 μ ν = g μ ν = 1 a 2 ( t ) g 00 g 0 i g i 0 g i j = 1 a 2 ( t ) 0 0 1 a 2 ( t ) δ i j
For μ = 0 , ν = 0 in Equation (541), we have
g 0 00 g 00 0 + g 0 0 i g i 0 0 + g 0 00 δ g 00 + g 0 0 i δ g i 0 + δ g 00 g 00 0 + δ g 0 i g i 0 0 + δ g 00 δ g 00 + δ g 0 i δ g 0 i = δ 0 0
which give
δ g 00 δ g 00 + δ g 0 i δ g i 0 = 1
substituting the values, we have
1 a 2 ( t ) 1 + X a 2 ( t ) 1 2 A + 1 a 2 ( t ) i Y a 2 ( t ) i B = 1
Simplifying and neglecting second order product terms 2 A X and i Y · i B , we get
X = 2 A
Again from Equation (541) for μ = 0 , ν = i , we have after simplification
g 0 00 g 0 i 0 + g 0 0 i g j i 0 + g 0 00 δ g 0 i + g 0 0 i δ g j i + δ g 00 g 0 i 0 + δ g 0 i g j i 0 + δ g 00 δ g 0 i + δ g 0 i δ g j i = δ i 0
δ g 00 δ g 0 i + δ g 0 j δ g j i = 0
Substituting the values
1 a 2 ( t ) 1 + X a 2 ( t ) i B + 1 a 2 ( t ) j Y a 2 ( t ) 1 2 ψ δ j i + D j i E = 0
neglecting the higher product terms 2 A · i B , i Y · 2 ψ and j D j i Y · E , we have
i B + i Y = 0
On integrating, we get
Y = B
Now from Equation (541) for μ = i , ν = j , we have
g 0 i 0 g 0 j 0 + g 0 i j g i j 0 + g 0 i 0 δ g 0 j 0 + g 0 i j δ g i j 0 + δ g i 0 g 0 j 0 + δ g i j g i j 0 + δ g i 0 δ g 0 j + δ g i j δ g i j = δ j i
The non-vanishing terms are
δ g i 0 δ g 0 j + δ g i k δ g k j = δ j i
substituting values suitable change of indices
1 a 2 ( t ) i Y a 2 ( t ) i B + 1 a 2 ( t ) 1 + 2 z δ i k + D i k E · a 2 ( t ) 1 2 ψ δ k j + D k j E = δ j i
using properties δ i k δ k j = δ j i , δ i k D k j = D j i , δ k j D i k = D k i and neglecting the higher order product terms i B · j B , 4 Z ψ , 2 Z D j i E , 2 ψ · D j i K and D j i K E , we have
1 2 ψ + 2 Z δ j i + E + K D j i = δ j i + 0 D j i
Comparing the coefficients of δ j i and D j i , we get Z = ψ and K = E , so that inverse metric of the perturbed line element becomes
δ g μ ν = 1 a 2 ( t ) δ g 00 δ g 0 i δ g i 0 δ g i j = 1 a 2 ( t ) 1 + 2 A 1 a 2 ( t ) i B 1 a 2 ( t ) i B 1 a 2 ( t ) 1 + 2 ψ δ i j D i j E

12.2. The Unperturbed Line Element

The unperturbed line element is given by
d s 2 = a 2 ( t ) d t 2 + δ i j d x i d x j
or
g μ ν 0 = g μ ν = a 2 ( t ) g 00 g 0 i g i 0 g i j = a 2 ( t ) 1 0 0 δ i j
The inverse of g μ ν 0 = g μ ν is simply g 0 μ ν = g μ ν as for diagonal unperturbed metric g μ ν = 1 g μ ν , so that we can write
g 0 μ ν = g μ ν = 1 a 2 ( t ) g 00 g 0 i g i 0 g i j = 1 a 2 ( t ) 0 0 1 a 2 ( t ) δ i j
We calculate now affine connections-the 2nd kind of Christoffel symbols
Γ μ ν σ = g σ λ Γ μ ν λ = 1 2 g σ λ g μ λ , ν + g ν λ , μ + g μ ν , λ
We can compute the following possible components
Γ 00 0 Γ i i i Γ j j j Γ 00 i Γ 00 j Γ 0 j i Γ i 0 j Γ i j 0 Γ 0 i 0 Γ 0 j 0 Γ j k i Γ i k j Γ i j k
For σ = μ = ν = 0 , we have
Γ 00 0 = 1 2 a 2 ( t ) , 0 g 00 = a ˙ a
For σ = i , μ = 0 , ν = j , we have
Γ 0 j i = 1 2 a 2 ( t ) δ i k , 0 a 2 ( t ) δ j k = a ˙ a δ j i
For σ = 0 , μ = i , ν = j , we have
Γ i j 0 = 1 2 a 2 ( t ) , 0 a 2 ( t ) δ i j = a ˙ a δ i j
For σ = i , μ = 0 , ν = 0 , we have
Γ 00 i = 0
and
For σ = 0 , μ = 0 , ν = i , we have
Γ 0 i 0 = 0
and for σ = i , μ = j , ν = k , we have
Γ j k i = 0
Now, we can calculate perturbed element when we have necessary components of perturbed and unperturbed metric and its inverse metric tensor components,
δ Γ μ ν σ = δ g σ λ Γ μ ν λ = 1 2 δ g σ λ g μ λ , ν + g ν λ , μ + g μ ν , λ = 1 2 δ g σ λ g μ λ , ν + g ν λ , μ g μ ν , λ + 1 2 g σ λ δ g μ λ , ν + δ g ν λ , μ + δ g μ ν , λ
We are on our stake now to calculate the following components:
δ Γ 00 0 δ Γ i i i δ Γ j j j δ Γ 00 i δ Γ 00 j δ Γ 0 j i δ Γ i 0 j δ Γ i j 0 δ Γ 0 i 0 δ Γ 0 j 0 δ Γ j k i δ Γ i k j δ Γ i j k
The non-vanishing components are For σ = μ = ν = 0 , we have
δ Γ 00 0 = A ˙
For σ = i , μ = 0 , ν = j , we have
δ Γ 0 j i = ψ ˙ δ j i + 1 2 D i j E ˙
For σ = 0 , μ = i , ν = j , we have
δ Γ i j 0 = 2 a ˙ a A δ i j i j B 2 a ˙ a ψ δ i j ψ δ i j a ˙ a D i j E + 1 2 D i j E ˙
For σ = i , μ = 0 , ν = 0 , we have
δ Γ 00 i = a ˙ a i B + i B ˙ + i A
and For σ = 0 , μ = 0 , ν = i , we have
δ Γ 0 i 0 = i A + a ˙ a i B
and for σ = i , μ = j , ν = k , we have
δ Γ j k i = j ψ δ k i k ψ δ j i + i ψ δ j k a ˙ a i B δ j k + 1 2 j D k i E + 1 2 k D j i E + 1 2 i D j k E
Now the unperturbed Ricci tensor is given
R μ ν = g λ σ R μ ν λ σ = , ν Γ μ λ σ , λ Γ μ ν σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ
For μ = 0 , ν = 0 , we have
R 00 = , σ Γ 00 σ , 0 Γ 0 σ σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ = 0
For μ = 0 , ν = i , we have
R 0 i = , σ Γ 00 σ , 0 Γ 0 σ σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ = 0
For μ = i , ν = j , we have
R i j = , σ Γ 00 σ , 0 Γ 0 σ σ + Γ μ λ n Γ n ν σ Γ μ ν n Γ n λ σ = 0
Now. we calculate the perturbed Ricci tensor components:
δ R μ ν = , ν δ Γ μ λ σ , λ δ Γ μ ν σ + δ Γ μ λ n Γ n ν σ + Γ μ λ n δ Γ n ν σ δ Γ μ ν n Γ n λ σ Γ μ ν n δ Γ n λ σ
For μ = 0 , ν = 0 , we have
δ R 00 = , ν δ Γ μ λ σ , λ δ Γ μ ν σ + δ Γ μ λ n Γ n ν σ + Γ μ λ n δ Γ n ν σ δ Γ μ ν n Γ n λ σ Γ μ ν n δ Γ n λ σ = a ˙ a i i B + i i B + i i A + 3 ψ + 3 a ˙ a ψ + 3 a ˙ a A
For μ = 0 , ν = i , we have
δ R 0 i = , ν δ Γ μ λ σ , λ δ Γ μ ν σ + δ Γ μ λ n Γ n ν σ + Γ μ λ n δ Γ n ν σ δ Γ μ ν n Γ n λ σ Γ μ ν n δ Γ n λ σ = a ¨ a i B + a ˙ a 2 i B + 2 i ψ + 2 a ˙ a i A + 1 2 k D i k E
For μ = i , ν = j , we have
δ R i j = , ν δ Γ μ λ σ , λ δ Γ μ ν σ + δ Γ μ λ n Γ n ν σ + Γ μ λ n δ Γ n ν σ δ Γ μ ν n Γ n λ σ Γ μ ν n δ Γ n λ σ = a ˙ a A ˙ 5 a ˙ a ψ ˙ 2 a ¨ a A 2 a ˙ a 2 A 2 a ¨ a ψ 2 a ˙ a 2 ψ ψ ¨ + k k ψ a ˙ a k k B δ i j i j B ˙ + a ˙ a D i j E ˙ + a ¨ a D i j E + a ˙ a 2 D i j E + 1 2 D i j E ¨ + i j ψ i j A 2 a ˙ a i j B + 1 2 k i D j k E + 1 2 k j D i k E 1 2 k k D i j E
Unperturbed Ricci scalar is obtained by contracting the unperturbed Ricci tensor
R = g μ ν R μ ν
Using double sum and simplifying, we have
R = g 00 R 00 + g 11 R 11 + g 22 R 22 + g 33 R 33 = 6 a ¨ a 3
and perturbed Ricci scalar, using double sum and simplifying, we have
δ R = δ g μ ν R μ ν + g μ ν δ R μ ν = δ g 00 R 00 + δ g 11 R 11 + δ g 22 R 22 + δ g 33 R 33 + g 00 δ R 00 + g 11 δ R 11 + g 22 δ R 22 + g 33 δ R 33 = 6 a ˙ a 3 i i B 2 a 2 i i B 2 a 2 i i A 6 a 2 ψ ¨ 6 a ˙ a 3 A ˙ 18 a ˙ a 3 ψ ˙ 12 a ¨ a 3 A + 4 a 2 i i ψ + 1 a 2 k i D i k E
now unperturbed Einstein tensor is
G μ ν = R μ ν + 1 2 g μ ν R
For μ = 0 , ν = 0 , we have
G 00 = R 00 + 1 2 g 00 R = 3 a ˙ a 2
For μ = 0 , ν = i , we have
G 0 i = R 0 i + 1 2 g 0 i R = 0
For μ = i , ν = j , we have
G i j = R i j + 1 2 g i j R = 2 a ¨ a + a ˙ a 2 δ i j
Now, the perturbed Einstein tensor at first order of perturbation is
δ G μ ν = δ R μ ν 1 2 δ g μ ν R + g μ ν δ R
For μ = 0 , ν = 0 , we have
δ G 00 = δ R 00 1 2 δ g 00 R + g 00 δ R = 2 a ˙ a i i B 6 a ˙ a ψ ˙ + 2 i i ψ + 1 2 k i D i k E
For μ = 0 , ν = i , we have
δ G 0 i = δ R 0 i 1 2 δ g 0 i R + g 0 i δ R = 2 a ¨ a i B + a ˙ a 2 i B + 2 i ψ ˙ + 1 2 k D i k E ˙ + 2 a ˙ a i A
For μ = i , ν = j , we have
δ G i j = δ R i j 1 2 δ g i j R + g i j δ R = 2 a ˙ a A ˙ + 4 a ˙ a ψ ˙ + 4 a ¨ a A 2 a ˙ a 2 A + 4 a ¨ a ψ 2 a ˙ a 2 ψ + 2 ψ ¨ k k ψ + 2 a ˙ a k k B + k k B ˙ + k k A + 1 2 k ρ D ρ k E δ i j i j B ˙ + a ˙ a D i j E ˙ 2 a ¨ a D i j E i j A + i j ψ + a ˙ a 2 D i j E + 1 2 D i j E ¨ 2 a ˙ a i j B + 1 2 k i D j k E + 1 2 k j D i k E 1 2 k k D i j E
Now, the unperturbed stress energy tensor is
T μ ν = μ ϕ ν ϕ g μ ν 1 2 g σ ρ σ ϕ ρ ϕ V ϕ
For μ = 0 , ν = 0 , we have
T 00 = 0 ϕ 0 ϕ g 00 1 2 g σ ρ σ ϕ ρ ϕ V ϕ = 1 2 ϕ ˙ 2 + a 2 ( t ) V ϕ
For μ = 0 , ν = i , we have
T 0 i = 0 ϕ i ϕ g 0 i 1 2 g σ ρ σ ϕ ρ ϕ V ϕ = 0
For μ = i , ν = j , we have
T i j = i ϕ j ϕ g i j 1 2 g σ ρ σ ϕ ρ ϕ V ϕ = 1 2 ϕ ˙ 2 a 2 ( t ) V ϕ δ i j
Now, the perturbed stress energy tensor at first order of perturbation is
δ T μ ν = μ δ ϕ ν ϕ + μ ϕ ν δ ϕ δ g μ ν 1 2 g σ ρ σ ϕ ρ ϕ + V ϕ g μ ν 1 2 δ g σ ρ σ ϕ ρ ϕ + g σ ρ σ δ ϕ ρ ϕ + g σ ρ σ ϕ ρ δ ϕ + ϕ δ V ϕ + ϕ V ϕ δ ϕ
For μ = 0 , ν = 0 , we have
δ T 00 = ϕ ˙ δ ϕ ˙ + 2 a 2 ( t ) A V ϕ + δ ϕ a 2 ( t ) V ϕ ϕ
For μ = 0 , ν = i , we have
δ T 0 i = ϕ ˙ i δ ϕ + ϕ ˙ 2 2 i B a 2 ( t ) V ϕ i B
For μ = i , ν = j , we have
δ T i j = ϕ ˙ δ ϕ ˙ ϕ ˙ 2 A a 2 ( t ) V ϕ ( ϕ ) δ ϕ ϕ ˙ 2 ψ + 2 a 2 ( t ) V ( ϕ ) ψ δ i j a 2 ( t ) V ( ϕ ) D i j E + 1 2 ϕ ˙ 2 D i j E
Therefore, the perturbed Einstein field equations are
δ G μ ν = 8 π δ T μ ν
comparing the corresponding components, we have
2 a ˙ a i i B 6 a ˙ a ψ ˙ + 2 i i ψ + 1 2 k i D i k E = 8 π ϕ ˙ δ ϕ ˙ + 2 a 2 ( t ) A V ϕ + δ ϕ a 2 ( t ) V ϕ ϕ
2 a ¨ a i B + a ˙ a 2 i B + 2 i ψ ˙ + 1 2 k D i k E ˙ + 2 a ˙ a i A = 8 π ϕ ˙ i δ ϕ + ϕ ˙ 2 2 i B a 2 ( t ) V ϕ i B
2 a ˙ a A ˙ + 4 a ˙ a ψ ˙ + 4 a ¨ a A 2 a ˙ a 2 A + 4 a ¨ a ψ 2 a ˙ a 2 ψ + 2 ψ ¨ k k ψ + 2 a ˙ a k k B + k k B ˙ + k k A + 1 2 k ρ D ρ k E δ i j i j B ˙ + a ˙ a D i j E ˙ 2 a ¨ a D i j E i j A + i j ψ + a ˙ a 2 D i j E + 1 2 D i j E ¨ 2 a ˙ a i j B + 1 2 k i D j k E + 1 2 k j D i k E 1 2 k k D i j E = 8 π ϕ ˙ δ ϕ ˙ ϕ ˙ 2 A a 2 ( t ) V ϕ ( ϕ ) δ ϕ ϕ ˙ 2 ψ + 2 a 2 ( t ) V ( ϕ ) ψ δ i j a 2 ( t ) V ( ϕ ) D i j E + 1 2 ϕ ˙ 2 D i j E
perturbed equations can also be determined in mixed tensor form. Expressing Equation (606) in mixed form and working out on the same lines as we did earlier
δ G μ ν = 8 π δ T μ ν

13. Summary

Relativistic cosmology was founded on the general theory of relativity with the introduction of the cosmological principle and Weyl’s principle implicitly implied. In the beginning, Einstein’s and de Sitter’s cosmological models were presented, though now of historical interest, yet they both are very significant as the first initiates the modern cosmology relativistically and scientifically and the latter, later on, was used to provide the initial conditions of the big bang model with a slight change. The first theoretical models for the possibility of dynamic universe evolved beginning with Friedmann, Lemaitre and were observationally determined by E. Hubble. In 1929, E. Hubble found exactly the same expanding universe that Friedmann did theoretically in 1922. Therefore, it was Friedmann who championed the cause of dynamical universes; however, his work was recognized later when he was no more in the world. The theory of big bang based on the standard cosmological model faces Horizon, Flatness, Entropy problems etc. To resolve these problems a phase of exponentially expanding universe was introduced in its very early history which occurred in a very small fraction of time (about 1 10 43 s of the very 1st second after time creation) known as inflation. The inflation is identified as the initial conditions under which the big bang might have taken place. The introduction of inflation caused the name inflationary cosmology and it is about forty years since its birth to date. The inflationary paradigm stands now on firm observational footing and is accepted irrevocably in cosmology as the viable description for the early universe. Starobinsky, Guth, and Linde are credited with setting the foundations of inflationary cosmology. The inflationary cosmology is being hailed as successful in explaining the origin of structure formation through cosmological quantum fluctuations as relicts of cosmic inflation. The observations conducted on microwave background radian and the recent discoveries of gravitational waves and black holes lend the confirmatory support to the underlying principles of the inflationary cosmology. Dark energy is the one of most challenging issues of the standard cosmology both on theoretical and observational grounds. In the framework of Λ CDM it has equation of state (EoS) w = 1 , however δ is facing fine-tuning problem. An alternative remedy to tackle the problems of δ are the model consisting of canonical and non-canonical scalar fields. The scalar field models modify matter sector of EFE on the right hand side, nonetheless in f R geometry is modified as curvature of spacetime. Λ CDM model is accepted for its being in good agreement with the recent observations. Note that there exists well-elaborated scenario to unify inflation with Dark Energy in modified gravity which was first proposed in S. Nojiri and S.D. Odintsov [61].

Author Contributions

All of the authors shared their contribution to the design and structure of the study concept of the manuscript. The 2nd author I.A. presented the idea of the article and the 1st author M.Z.M. did search for the literature, analyzed the data, and drafted with revision of the data. The 3rd author J.L.G.G. contributed in the conclusion section with financial support. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by Ministerio de Ciencia, Innovación y Universidades, grant number PGC2018-097198-B-I00, and by Fundación Séneca of Región de Murcia, grant number 20783/PI/18.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Mirza Faisal Shafique (SSE-Physics) Govt. High School Chawinda, Distt. Sialkot Punjab for having useful discussions and for the help rendered by him in composing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Space, Time and Spacetime

A background arena of space and time is necessarily required for all the physical phenomena to play over it and the compatibility of the known physical laws is made with structure of space and time. Space, time, and motion are concomitant ingredients cohered to matter and can never be disengaged from each other. The universe exists in space and evolves in time so that universe, space and time are in separable from each other and are coherently related to each other. Space is understood as possessing three dimensions and time is speculated to have only one dimension. Therefore, Newtonian Mechanics has been formulated in such a way to consider the spatial dimensions existing independently from the only one dimension of time. The Euclidean geometry provides necessary mechanism in dealing with such notions of space and time. In this regard Euclidean space becomes important which proposes three independent perpendicular dimensions of space and the dimension of time does not get affected by it. Space and time are envisaged as independent absolute entities which are not affected by each other. The Euclidean structure of space is flat and distances are measured by using the standard Pythagoras theorem for three dimensions as
d s 2 = x 2 + y 2 + z 2 ,
or in differential of the distances
d s 2 = d x 2 + d y 2 + d z 2 ,
where d s = x , y , z or d s = d x , d y , d z , respectively. The time coordinate does appear anywhere in this distance-measuring formula which means in the geometry of space, the dimension of time will be dealt separately. Newton’s notions of space and time as described in Principia Mathematica are given as “Absolute space, in its own nature, without regard to anything external, remains always similar and immovable. Relative space is some movable dimension or measure of the absolute spaces which our senses determine by its position to bodies: and which is vulgarly taken for immovable space. Absolute motion is the translation of a body from one absolute place into another: and relative motion, the translation from one relative place into another” and absolute time is defined in these words “Absolute, true and mathematical time, of itself, and from its own nature flows equably without regard to anything external, and by another name is called duration. Relative, apparent and common time, is some sensible and external (whether accurate or inequable) measure of duration by the means of motion, which is commonly used instead of true time”.
In 1905, Einstein’s paper entitled “On the electrodynamics of moving bodies” put forth on the base of two postulates that time might be dealt on equal footing with space as one of the dimensions of space. Minkowski (1864–1909) translated the mixing of space and time coordinates as requiring a four-dimensional scenario where physical phenomena take place and the geometry of such four dimensional spacetime, where time is one dimension, is described by spacetime interval which is the generalized form of Pythagoras theorem
d s 2 = d t 2 + d x 2 + d y 2 + d z 2
or
d s 2 = η μ ν d x μ d x ν
where
η μ ν = η 00 η 01 η 02 η 03 η 10 η 11 η 12 η 13 η 20 η 21 η 22 η 23 η 30 η 31 η 32 η 33 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
Minkowski first understood that the spacetime interval given in Equation (A3) remains invariant for all the observers and carries the similar meaning for all the observers in uniform relative motion, however, Einstein considered either with respect time or space the interval does not remain identical for all relative observers in uniform motion. Minkowski avowedly said in a conference addressing to the German scientists that “Ladies and gentlemen! the views of space and time which i wish to lay before you have sprung from the soil of experimental physics, therein lies their strength, they are radical. Henceforth space by itself and time by itself are doomed to fade away into mere shadows and only a union of the two will preserve an independent reality” [42]. General relativity was formulated on the base of four dimensional spacetime as Minkowski has laid it but in order to incorporate the gravity into it Einstein utilized the power of tensors and modeled the curved geometry of spacetime describing its curvature as gravity. The geometry of curved spacetime is encoded into a two rank symmetric tensor known as fundamental tensor and given as the spacetime metric or line element as
d s 2 = g μ ν d x μ d x ν
where g μ ν is given by
g μ ν = g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 = g 00 0 0 0 0 g 11 0 0 0 0 g 22 0 0 0 0 g 33
In the absence of matter, the curvature of spacetime vanishes and the geometry of spacetime becomes flat, i.e., g μ ν = μ μ ν , yet non-Euclidean that is required by the special relativity.

Appendix B. Maximally Symmetric 3-Space (Spherically Symmetric Space)

In order to have a space more symmetrical we require comparatively lesser number of functions as much as possible to determine its properties. It is the curvature of a space and its nature that determines how much the space is symmetric maximally. If the curvature K of a space under consideration does not depend upon the coordinates of the points constituting it and has a constant value, then the space shall be maximally symmetric and the spaces possessing the curvature of this kind logically entail cosmological principle, i.e., homogeneity and isotropy. Spacelike coordinates x 1 , x 2 , x 3 obviously span 3-space which we require to be maximally symmetric. The Riemann curvature tensor R μ ν ρ σ in three dimensional space has 3 4 = 81 components which depend on the coordinates. From these, only six components are independent and are the functions of coordinates and require six functions to be described in order to specify intrinsically the geometric properties of the three dimensional space. The Riemann curvature tensor R μ ν ρ σ depends on curvature K and the metric tensor g μ ν for the maximally symmetric spaces which is the simplest form for it to adopt. It is given by
R μ ν ζ π = K g μ ζ g ν π g μ π g ν ζ
g μ π R μ ν ζ π = K g μ π g μ ζ g ν π g μ π g ν ζ = K g μ π g μ ζ g ν π g μ π g μ π g ν ζ
R ν ζ = K δ ζ π g ν π δ ζ ζ g ν ζ = K g ν ζ δ 1 1 + δ 2 2 + δ 3 3 g ν ζ
R ν ζ = K g ν ζ 3 g ν ζ = K 2 g ν ζ
Then, Ricci scalar or curvature scalar from above Equation (A11) can be had by contraction with inverse metric tensor g ν ζ
g ν ζ R ν ζ = 2 g ν ζ g ν ζ K
R = 2 δ ζ ζ K = 2 δ 1 1 + δ 2 2 + δ 3 3 K = 2 1 + 1 + 1 K = 6 K
The metric of an isotropic 3-space must depend on rotational invariants given by
x · x = r 2 d x · d x , x · d x
and in spherical polar coordinates r , θ , ϕ , it should take the form
d σ 2 = C r x d x 2 + D r d x d x 2
d σ 2 = C r r 2 d r 2 + D r d r 2 + r 2 d θ 2 + r 2 sin 2 θ d θ 2
Redefining the radial coordinate r ¯ 2 = r 2 D r and dropping the bars on the variables, we can write the above Equation (A16) in the form
d σ 2 = B r d r 2 + r 2 d θ 2 + r 2 sin 2 θ d θ 2
where B(r) is an arbitrary function of r. Solving the metric in Equation (A17), the components are
g 11 = B ( r ) g 22 = r 2 g 33 = r 2 sin 2 θ
The non-vanishing Christoffel symbols we find, are
Γ 11 1 = Γ r r r = 1 2 B ( r ) d B ( r ) d r Γ 22 1 = Γ θ θ r = r B ( r ) Γ 33 1 = Γ ϕ ϕ r = r sin 2 θ B ( r ) Γ 13 2 = Γ r ϕ θ = 1 r Γ 13 3 = Γ r ϕ ϕ = 1 r Γ 33 2 = Γ ϕ ϕ θ = sin θ cos θ Γ 32 2 = Γ ϕ θ ϕ = cot θ
Now, from the Ricci tensor
R μ ν = ν Γ μ ρ ρ ρ Γ μ ν ρ + Γ μ ρ σ Γ σ ν ρ + Γ μ ν σ Γ σ ρ ρ
We calculate non-vanishing components
R 11 = R r r = 1 r B d B d r R 22 = R θ θ = 1 + 1 B r 2 B 2 d B d r R 33 = R ϕ ϕ = R θ θ sin 2 θ = 1 + 1 B r 2 B 2 d B d r sin 2 θ
and the Ricci scalar is
R = 2 δ ζ ζ K
1 r B d B d r = 2 K B ( r ) 1 + r r B 2 d B d r 1 B = 2 K r 2
Integrating 1st part of Equation (A23), we obtain
B ( r ) = 1 A K r 2
where A being a constant of integration can be found by substituting Equation (A24) into 2nd part of Equation (A23), we get
1 A + K r 2 = K r 2 A = 1
so we obtain the metric
d σ 2 = d r 2 1 K r 2 d r 2 + r 2 d θ 2 + r 2 sin 2 θ d θ 2
Equation (A26) incorporates a hidden symmetry characterized by homogeneity and isotropy, and represents the line element of a maximally symmetric 3-space. Due to arbitrary origin of radial coordinate system, we considered and due to symmetry of space we can take all the points of this space equivalent and the origin of this coordinate system can be chosen arbitrarily at any point which means that there exists no center in this space. Therefore the maximally symmetric space is infinite and open. Further the line element is equivalent perfectly to the metric of a 3-sphere embedded in a four dimensional Euclidean space which has spherical symmetry as well.

Appendix C. Spectrum of the Black Body

A blackbody can absorb hypothetically radiation of all wavelengths falling on it and reflecting nothing at all. How at the different wavelengths distribution of radiation occurs in a blackbody is given below in Figure A1:
Figure A1. Radiation distribution of blackbody at different wavelengths.
Figure A1. Radiation distribution of blackbody at different wavelengths.
Universe 07 00276 g0a1
In the early universe when matter and radiation decoupled from each other, the so-called decoupling, the primordial radiation given off gives a snapshot of the universe at that time and is known as cosmic microwave background radiation (CMBR) observed accidentally in the 60 s. The recent observations conducted on cosmic microwave background radiation reveals the fact that this is the perfect black body radiation with a temperature of 2.7255 Kelvin on average. We know that the wavelength distribution of a black body is given by
u λ , T d λ = 8 π h c λ 5 1 e h c λ k B T 1 d λ
where u λ , T d λ is the energy per unit volume of the radiation with wavelength between λ and λ + d λ emitted by a blackbody at temperature T. We consider now a black body radiation from the big bang when the universe first became transparent to photons after 400,000 years after big bang to this time about 4,000,000,000 years. The wavelength of the primordial photons λ is Doppler shifted to λ due to expansion of universe, certainly λ > λ . Let f λ , T d λ be the current per unit volume of the residual big bang radiation as measured from the earth. As the shell of charged particles that emitted the radiation is moving away from the Earth at extremely relativistic speed so we should use the relativistic Doppler shift for light from a receding source to relate λ to λ that is
λ = 1 + v / c 1 v / c λ = B λ
where we put B = 1 + v / c 1 v / c , and v is the speed of recession of the charged shell. As v < c , clearly λ > λ by a factor
1 + v / c 1 v / c
Equation (A29) can be interpreted by generalization that all the distances have grown since first radiation emitted. In order to have a relation between currently observed spectrum f λ , T d λ and original black body radiation distribution
u λ , T d λ
we put from Equation (A28) λ = λ B into Equation (A27)
u λ , T d λ = 8 π h c λ B 5 1 e h c λ B k B T 1 d λ B
u λ , T d λ B 4 = 8 π h c λ 5 1 e h c λ k B T 1 d λ
where T = T B and RHS of Equation (A32) can be identified with current black body spectrum f λ , T d λ which has standard functional form of a blackbody spectrum with wavelength λ and temperature T . Equation (A30) becomes
u λ , T d λ B 4 = f λ , T d λ
Equation (A33) says that the radiation from a receding blackbody has same spectral distribution as yet but its temperature T and energy
u λ , T d λ
dropped by factors of B and B 4 respectively.

Appendix D. Big Bang Theory of Creation

Historically the name of this theory as big bang is due to Fred Hoyle (1915–2001), one of the inventors and staunch proponents of steady state theory who coined the term accidentally with showing abhorrence towards it. The steady-state theory, once a rival theory of the big bang, lends support to an eternally evolving universe without a beginning and an end. The big bang theory explains the evolutionary phases of the universe beginning with a very small span of a fraction of very first second to the present age. The warp and woof of the theory is woven by equation of general relativity and the developments made in its context. The theory traces its theoretical origin back to Friedmann equations and the discovery of expansion of the universe by Edwin Hubble in thirties. Further it rests upon the relative abundance of light elements by George Gamow forties and CMB accidental discovery in seventies by Penzias and Wilson. The theory comes forth on the base of standard cosmological model and describes that our universe had had a beginning and had erupted from an extremely dense, point-like singularity about 14 billion years ago. At the singularity state, all basic interactions of nature had coalesced symmetrically where all the matter-energy melted down into an indistinguishable quark-gluon primordial soup. Einstein has expressed his views on the nature of this singularity in his later years: “The theory is based on a separation of the concepts of the gravitational field and matter. While this may be a valid approximation for weak fields, it may presumably be quite inadequate for very high densities of matter. One may not therefore assume the validity of the equations for very high densities and it is just possible that in a unified theory there would be no such singularity” [62]. It is speculated that during the Planck time of the order of 10 43 s all the forces of nature, namely, electroweak nuclear, strong nuclear, electromagnetic, and gravitational, were so merged into one another such that they were indistinguishable bearing perfect symmetry. From the beginning of time, t = 0 s to Planck time t p 10 43 s within the time span of very first second is known as the Trans-Planckian era whose physics is yet incomplete and is open hitherto to investigation. It is being conjectured that during the time ranging from 10 43 s to 10 35 s, the gravitational force freed itself from the rest of interactions, and during this period there exist the particles that supersymmetry predicts and are known as quarks, leptons, their antiparticles, and some certain massive particles. After the time interval that begins with 10 35 s to some shortly later time 10 32 s, the universe expanded exponentially and gradually cooled down where the strong and electroweak forces get separated from the rest. As the universe continues to cool after the big bang, around the time 10 10 s, the electroweak force splits into weak force and electromagnetic force and within few minutes after it, protons and neutrons start to condense out of the cooling quark–gluon plasmic soup. During the first half of creation, the universe can be viewed as a thermonuclear bomb fusing protons and neutrons into deuterium and then helium producing most of the helium nuclei that exist now. After the big bang until about 400,000 years radiation-dominated era prevailed. Vibrant photonic radiation halted itself to become a clumped matter rather even forming single atom hydrogen or helium due to photon–atom collisions which would result in ionization instantly in the case if any atom happened to form, therefore no chance occurs for the formation of atoms and the universe remains opaque to electromagnetic radiation due to incessant Compton scattering experienced by photons with free electrons that abound in. On further cooling electrons could bind to protons forming helium nuclei with the reduction in the number of charged particles, absorption or scattering of photons consequently the universe suddenly became transparent to photons and radiation dominated era diminishes and neutral matter domination begins in the form of atoms, molecules, gas clouds, stars and in the end galaxies-the universe today. This is the whole saga of the big bang theory of creation.

References

  1. Einstein, A. The Foundation of General Theory of Relativity. In The Berlin Years: Writings; Springer: Dordrecht, The Netherlands, 1916; Volume 6, pp. 1914–1917. [Google Scholar]
  2. Mukhanov, V.F. Physical Foundations of Cosmology; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  3. Einstein, A. Cosmological considerations in the General Theory of Relativity. Volume 6: The Berlin Years: Writings, 1914–1917, Doc 43. In Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys.); Pinceton Univeristy Press: Pinceton, NJ, USA, 1917; pp. 421–450. [Google Scholar]
  4. de Sitter, W. On Einstein’s Theory of Gravitation and its Astronomical Consequences. Third Paper. Mon. Not. R. Astron. Soc. 1917, 78, 3–28. [Google Scholar] [CrossRef] [Green Version]
  5. Friedmann, A. Über die Möglichkeit einer Welt mit konstanter negativer Krümmung des Raumes. Z. Phys. 1924, 21, 326–332, English translation: “On the Possibility of a World with Constant Negative Curvature of Space”. Gen. Rel. Grav. 1999, 31, 2001–2008.. [Google Scholar] [CrossRef]
  6. Lemaître, G. Un univers homogène de masse constante et de rayon croissant, rendant compte de la vitesse radiale des nébuleuses extra-galactiques. Ann. Soc. Sci. Brux. 1927, 47A, 49–59, English translation: “A homogeneous universe with constant mass and increasing radius accounting for the radial velocities of the radial velocities of the extra-galactic nebulae”. Gen. Rel. Grav. 2013, 45, 1635–1646.. [Google Scholar]
  7. Hubble, E. A relation between distance and radial velocity among extra-galactic nebulae. Proc. Natl. Acad. Sci. USA 1929, 15, 168–173. [Google Scholar] [CrossRef] [Green Version]
  8. Alpher, R.; Herman, R. Reflections on early work on ‘big bang’ cosmology. Phys. Today 1988, 41, 24–34. [Google Scholar] [CrossRef]
  9. Cai, Y.; Li, H.G.; Qiu, T.; Piao, Y.S. The effective field theory of nonsingular cosmology: II. Eur. Phys. J. C 2017, 77, 369. [Google Scholar] [CrossRef]
  10. Gliner, E. Algebraic properties of the energy-momentum tensor and vacuumlike states of matter. Sov. J. Exp. Theor. Phys. 1966, 22, 378. [Google Scholar]
  11. Kirzhnits, D.A.; Linde, A.D. Macroscopic Consequences of the Weinberg Model. Phys. Lett. B 1972, 42, 471. [Google Scholar] [CrossRef]
  12. Linde, A.D. Is the Cosmological Constant Really a Constant? JETP Lett. 1974, 19, 183. [Google Scholar]
  13. Brout, R.; Englert, F.; Gunzig, E. The creation of the Universe as a Quantum phenomenon. Ann. Phys. 1978, 115, 78–106. [Google Scholar] [CrossRef]
  14. Linde, A.D. Phase Transitions In Gauge Theories and Cosmology. Rept. Prog. Phys. 1979, 42, 389. [Google Scholar] [CrossRef]
  15. Starobinsky, A.A. A new type of isotropic cosmological models without singularity. Phys. Lett. B 1980, 91, 99–102. [Google Scholar] [CrossRef]
  16. Guth, A.H. The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems. Phys. Rev. D 1981, 23, 347. [Google Scholar] [CrossRef] [Green Version]
  17. Linde, A.D. The inflationary universe. Rep. Prog. Phys. 1984, 47, 925. [Google Scholar] [CrossRef] [Green Version]
  18. Hawking, S.W.; Moss, I.G.; Stewart, J.M. Bubble Collisions in the Very Early Universe. Phys. Rev. D 1982, 26, 2681. [Google Scholar] [CrossRef]
  19. Guth, A.H.; Weinberg, E.J. Could The Universe Have Recovered from a Slow First Order Phase Transition? Nucl. Phys. B 1983, 212, 321. [Google Scholar] [CrossRef]
  20. 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, 389–393. [Google Scholar] [CrossRef]
  21. Albrecht, A.; Steinhardt, P.J.; Turner, M.S.; Wilczek, F. Reheating an inflationary universe. Phys. Rev. Lett. 1982, 48, 1437. [Google Scholar] [CrossRef]
  22. Linde, A.D. Coleman-Weinberg Theory and A New Inflationary Universe Scenario. Phys. Lett. B 1982, 114, 431. [Google Scholar] [CrossRef]
  23. Linde, A.D. Temperature Dependence of Coupling Constants and the Phase Transition in the Coleman-Weinberg Theory. Phys. Lett. B 1982, 116, 340–342. [Google Scholar] [CrossRef] [Green Version]
  24. Linde, A.D. Scalar Field Fluctuations in Expanding Universe and the New Inflationary Universe Scenario. Phys. Lett. B 1982, 116, 335. [Google Scholar] [CrossRef]
  25. Albrecht, A.; Steinhardt, P.J. Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking. Phy. Rev. Lett. 1982, 48, 1220. [Google Scholar] [CrossRef]
  26. Hawking, S.W. The Development of Irregularities In A Single Bubble Inflationary Universe. Phys. Lett. B 1982, 115, 295. [Google Scholar] [CrossRef]
  27. Starobinsky, A.A. Dynamics of Phase Transition in the New Inflationary Universe Scenario and Generation of Perturbations. Phys. Lett. B 1982, 117, 175. [Google Scholar] [CrossRef]
  28. Guth, A.H.; Pi, S.Y. Fluctuations In The New Inflationary Universe. Phys. Rev. Lett. 1982, 49, 1110. [Google Scholar] [CrossRef]
  29. Linde, A.D. Towards inflation in string theory. Int. J. Phys. Conf. Ser. 2004, 24, 151–160. [Google Scholar] [CrossRef]
  30. Mukhanov, V.F.; Chibisov, G.V. Quantum Fluctuation And ‘Nonsingular’ Universe. JETP Lett. 1981, 33, 532. [Google Scholar]
  31. 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. Measurements of Ω and Λ from 42 high-redshift supernovae. Astrophys. J. 1999, 517, 565. [Google Scholar] [CrossRef]
  32. Perlmutter, S.; Turner, M.S.; White, M. Constraining dark energy with type Ia supernovae and large-scale structure. Phys. Rev. Lett. 1999, 83, 670. [Google Scholar] [CrossRef] [Green Version]
  33. 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, 1009. [Google Scholar] [CrossRef] [Green Version]
  34. Jedamzik, K.; Pogosian, L.; Zhao, G.B. Why reducing the cosmic sound horizon can not fully resolve the Hubble tension. arXiv 2020, arXiv:2010.04158. [Google Scholar]
  35. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results-VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar]
  36. Bamba, K.; Capozziello, S.; Nojiri, S.I.; Odintsov, S.D. Dark energy cosmology: The equivalent description via different theoretical models and cosmography tests. Astrophys. Space Sci. 2012, 342, 155–228. [Google Scholar] [CrossRef] [Green Version]
  37. Nojiri, S.I.; Odintsov, S.D. Introduction to modified gravity and gravitational alternative for dark energy. Int. J. Geom. Methods Mod. Phys. 2007, 4, 115–145. [Google Scholar] [CrossRef] [Green Version]
  38. Nojiri, S.; Odintsov, S.D.; Oikonomou, V.K. Modified gravity theories on a nutshell: Inflation, bounce and late-time evolution. Phys. Rep. 2017, 692, 1–104. [Google Scholar] [CrossRef] [Green Version]
  39. Simson, R. (Ed.) The Elements of Euclid; Desilver, Thomas: Philadelphia, PA, USA, 1838. [Google Scholar]
  40. Dray, T. The Geometry of Special Relativity; CRC Press: Boca Raton, FL, USA, 2012. [Google Scholar]
  41. Einstein, A. On the electrodynamics of moving bodies. Ann. Phys. 1905, 17, 891–921. [Google Scholar] [CrossRef]
  42. Minkowski, H. Space and Time: Minkowski’s Papers on Relativity; Minkowski Institute Press: Montreal, QC, Canada, 2013. [Google Scholar]
  43. Schrödinger, E. Space-Time Structure; Cambridge University Press: Cambridge, UK, 1985. [Google Scholar]
  44. Goenner, H. Weyl’s contributions to cosmology. In Hermann Weyl’s Raum—Zeit—Materie and a General Introduction to His Scientific Work; Birkhäuser: Basel, Switzerland, 2001; pp. 105–137. [Google Scholar]
  45. Kolb, E.W.; Turner, M.S. The early universe. Nature 1981, 294, 521–526. [Google Scholar] [CrossRef]
  46. Lesgourgues, J. Inflationary Cosmology. Lecture Notes of a Course Presented in the Framework of the “3i‘eme Cycle de Physique de Suisse Romande” Troisième cycle de la Physique en Suisse Romande. EPFL, March 2006. Available online: https://lesgourg.github.io/courses/Inflation_EPFL.pdf (accessed on 1 June 2021).
  47. Carroll, S.M. Spacetime and Geometry/Sean Carroll; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  48. Riotto, A. Inflation and the theory of cosmological perturbations. arXiv 2002, arXiv:hep-ph/0210162. [Google Scholar]
  49. Vilenkin, A. Classical and quantum cosmology of the Starobinsky inflationary model. Phys. Rev. D 1985, 32, 2511. [Google Scholar] [CrossRef]
  50. Liu, Z.G.; Piao, Y.S. Galilean islands in eternally inflating background. Phys. Rev. D 2013, 88, 043520. [Google Scholar] [CrossRef] [Green Version]
  51. Baumann, D. Cosmology: Part III mathematical tripos. Univ. Lect. Notes 2014, 56, 34. [Google Scholar]
  52. Pradhan, A. Inflationary Cosmology: Scalar Field Models and Structure Formation. Ph.D. Dissertation, Marieke Postma, Inflation, NIKHEF, Science Park 105 1098 XG, Amsterdam, The Netherlands, 2015. [Google Scholar]
  53. Vázquez, J.A.; Padilla, L.E.; Matos, T. Inflationary cosmology: From theory to observations. arXiv 2018, arXiv:1810.09934. [Google Scholar]
  54. Coles, P.; Lucchin, F. Cosmology: The Origin and Evolution of Cosmic Structure; John Wiley and Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  55. Chudaykin, A.; Gorbunov, D.; Nedelko, N. Exploring an early dark energy solution to the Hubble tension with Planck and SPTPol data. Phys. Rev. D 2021, 103, 043529. [Google Scholar] [CrossRef]
  56. Kowalski, M.; Rubin, D.; Aldering, G.; Agostinho, R.J.; Amadon, A.; Amanullah, R.; Balland, C.; Barbary, K.; Blanc, G.; Challis, P.J.; et al. Improved cosmological constraints from new, old, and combined supernova data sets. Astrophys. J. 2008, 686, 749. [Google Scholar] [CrossRef] [Green Version]
  57. Capozziello, S.; De Laurentis, M. Extended theories of gravity. Phys. Rep. 2011, 509, 167–321. [Google Scholar] [CrossRef] [Green Version]
  58. Odintsov, S.D.; Oikonomou, V.K.; Banerjee, S. Dynamics of inflation and dark energy from F (R, G) gravity. Nucl. Phys. B 2019, 938, 935–956. [Google Scholar] [CrossRef]
  59. Samart, D. Cosmological Dynamics from Modified f(R) Gravity in Einstein Frame. arXiv 2006, arXiv:astro-ph/0606612. [Google Scholar]
  60. Capozziello, S.; Nojiri, S.I.; Odintsov, S.D.; Troisi, A. Cosmological viability of f (R)-gravity as an ideal fluid and its compatibility with a matter dominated phase. Phys. Lett. B 2006, 639, 135–143. [Google Scholar] [CrossRef] [Green Version]
  61. Nojiri, S.I.; Odintsov, S.D. Modified gravity with negative and positive powers of curvature: Unification of inflation and cosmic acceleration. Phys. Rev. D 2003, 68, 123512. [Google Scholar] [CrossRef] [Green Version]
  62. Einstein, A. On the generalized theory of gravitation. Sci. Am. 1950, 182, 13–17. [Google Scholar] [CrossRef]
Figure 1. Pythagoras theorem: c 2 = a 2 + b 2 .
Figure 1. Pythagoras theorem: c 2 = a 2 + b 2 .
Universe 07 00276 g001
Figure 2. The rotation of two-dimensional rectangular coordinate system through angle θ .
Figure 2. The rotation of two-dimensional rectangular coordinate system through angle θ .
Universe 07 00276 g002
Figure 3. Three-dimensional rectangular Cartesian plane representing Euclidean space-three mutual perpendicular planes.
Figure 3. Three-dimensional rectangular Cartesian plane representing Euclidean space-three mutual perpendicular planes.
Universe 07 00276 g003
Figure 4. Newtonian space.
Figure 4. Newtonian space.
Universe 07 00276 g004
Figure 5. A spacetime frame as null cone structure.
Figure 5. A spacetime frame as null cone structure.
Universe 07 00276 g005
Figure 6. Structure of spacetime where one second of time along time axis equals 300,000 km along the space axis.
Figure 6. Structure of spacetime where one second of time along time axis equals 300,000 km along the space axis.
Universe 07 00276 g006
Figure 7. Curved spacetime around the Sun-spacetime in general relativity.
Figure 7. Curved spacetime around the Sun-spacetime in general relativity.
Universe 07 00276 g007
Figure 8. An illustration of the cosmological principle. The figure describes how homogeneity and isotropy of space are related.
Figure 8. An illustration of the cosmological principle. The figure describes how homogeneity and isotropy of space are related.
Universe 07 00276 g008
Figure 9. Illustration of the Weyl postulate.
Figure 9. Illustration of the Weyl postulate.
Universe 07 00276 g009
Figure 10. The geometry of 2-sphere embedded in three dimensional Euclidean space.
Figure 10. The geometry of 2-sphere embedded in three dimensional Euclidean space.
Universe 07 00276 g010
Figure 11. The spherical geometry Ω 0 > 1 and for hyperbolic geometry Ω 0 < 1 and Ω 0 = 1 represents flat geometry.
Figure 11. The spherical geometry Ω 0 > 1 and for hyperbolic geometry Ω 0 < 1 and Ω 0 = 1 represents flat geometry.
Universe 07 00276 g011
Figure 12. Diagrammatic scheme for deriving redshift.
Figure 12. Diagrammatic scheme for deriving redshift.
Universe 07 00276 g012
Figure 13. A source S radiating electromagnetic energy.
Figure 13. A source S radiating electromagnetic energy.
Universe 07 00276 g013
Figure 14. Angular diameter distance.
Figure 14. Angular diameter distance.
Universe 07 00276 g014
Figure 15. Inflationary universe.
Figure 15. Inflationary universe.
Universe 07 00276 g015
Figure 16. History of the universe beginning with big bang and expanding with inflation.
Figure 16. History of the universe beginning with big bang and expanding with inflation.
Universe 07 00276 g016
Figure 17. How the scalar field drive inflationary era.
Figure 17. How the scalar field drive inflationary era.
Universe 07 00276 g017
Figure 18. How the universe springs into being through a scalar field. Old Inflation: (a) The scalar field is in stable false vacuum (b) It shows that the scalar field causes inflation through quantum tunneling which comes to an end suddenly (c) Due to abrupt ending of inflation, energy is dissipated without evolution of the universe or an empty universe results. New Inflation: (d) The scalar field begins in right false vacuum (e) Despite quantum tunneling, the scalar field decays by slowly rolling down towards its minimum hence the name slow roll inflation (f) The energy does not dissipate, instead reheat occurs and the universe evolves to radiation and other phases.
Figure 18. How the universe springs into being through a scalar field. Old Inflation: (a) The scalar field is in stable false vacuum (b) It shows that the scalar field causes inflation through quantum tunneling which comes to an end suddenly (c) Due to abrupt ending of inflation, energy is dissipated without evolution of the universe or an empty universe results. New Inflation: (d) The scalar field begins in right false vacuum (e) Despite quantum tunneling, the scalar field decays by slowly rolling down towards its minimum hence the name slow roll inflation (f) The energy does not dissipate, instead reheat occurs and the universe evolves to radiation and other phases.
Universe 07 00276 g018
Figure 19. How inflation ends-slow roll inflation.
Figure 19. How inflation ends-slow roll inflation.
Universe 07 00276 g019
Figure 20. Plot of data for Λ CDM form SNe Ia, CMB, and BAO between Ω Λ and Ω m .
Figure 20. Plot of data for Λ CDM form SNe Ia, CMB, and BAO between Ω Λ and Ω m .
Universe 07 00276 g020
Table 1. Values of six parameters of Λ CDM from different sources.
Table 1. Values of six parameters of Λ CDM from different sources.
Sr. NoParameters Λ CDM
Planck-2018Planck-Low + SPTPlanck-Low + SPT + BOSSPlanck-Low + SPT + BOSS + SHOES
1 Ω b h 2 0.02233 ± 0.000150.02269 ± 0.000250.02250 ± 0.000210.02255 ± 0.00020
2 Ω c h 2 0.1198 ± 0.00120.1143 ± 0.00200.1166 ± 0.00120.1159 ± 0.0010
3 n s 0.9652 ± 0.00420.9785 ± 0.00740.9716 ± 0.00560.9735 ± 0.0054
4 ln 10 10 A S 3.043 ± 0.0143.021 ± 0.0173.014 ± 0.0173.008 ± 0.017
5 τ 0.0540 ± 0.00740.0510 ± 0.00860.0456 ± 0.00820.0437 ± 0.0087
6 100 θ M C 1.04089 ± 0.00031---
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mughal, M.Z.; Ahmad, I.; García Guirao, J.L. Relativistic Cosmology with an Introduction to Inflation. Universe 2021, 7, 276. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7080276

AMA Style

Mughal MZ, Ahmad I, García Guirao JL. Relativistic Cosmology with an Introduction to Inflation. Universe. 2021; 7(8):276. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7080276

Chicago/Turabian Style

Mughal, Muhammad Zahid, Iftikhar Ahmad, and Juan Luis García Guirao. 2021. "Relativistic Cosmology with an Introduction to Inflation" Universe 7, no. 8: 276. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7080276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop