In this chapter we consider the issues of the origin and evolution of relic gravitational waves (GW), which appear as a result of quantum fluctuations of the scalar field and the corresponding perturbations of the space-time metric at the early inflationary stage of the evolution of the universe. The main provisions of the inflationary paradigm and the methods of the construction of current cosmological models on its basis are considered. The influence of relic gravitational waves on the anisotropy and polarization of the relic radiation and the importance of estimating such an effect on the verification of cosmological models are discussed as well.
- scalar field
- cosmological perturbations
- gravitational waves
The general relativity (GR) theory proposed by A. Einstein more than a 100 years ago currently finds new confirmations. The possible existence of gravitational waves was predicted by A. Einstein on the basis of solving the equations of general relativity when calculating the power of gravitational radiation [1, 2].
Gravitational waves (GW) are space-time curvature disturbances, which propagate at the speed of light. They occur at any movements of material bodies, leading to inhomogeneous gravity force variation in the environment. Gravitational radiation was predicted by A. Einstein in the general relativity (GR) theory, but so far not detected by direct measurements.
According to general relativity, space-time is curved around the bodies due to the action of gravity and is represented by a symmetric tensor with 10 independent components. However, far from the masses (the case of weak gravitational fields), the tensor can be divided into two terms where the first term, i.e., tensor , corresponds to the flat space-time of the special theory of relativity and has only four components. The second tensor contains information about the curvature caused by the gravitational field and makes small corrections. In the case of gravitational disturbances propagating far from their sources, the components of the tensor can be calculated by the method proposed by Einstein , similar to that used in electrodynamics for delayed potentials.
The first evidence was received by experimental studies of Joseph H. Taylor and Joel M. Weisberg et al., who studied the effect of slowing down the period of the binary star system PSR 1913 + 16 due to energy losses on gravitational radiation . Until recently, however, there has remained the main task: the direct recording of gravitational waves from space radiation sources by means of ground-based or space gravitational antennas.
Over the years, several methods have been proposed for recording gravitational radiation. Experimental work began in the 1960s of the twentieth century, but before the beginning of the twenty-first century, there was no reliable experimental proof of the ground-based recording of gravitational radiation .
This is due to the fact that gravitational waves have small amplitude; in addition, the proposed detection methods have insufficient sensitivity and are rather complicated in technical implementation.
These are broadband gravitational antennas, which offer a lot of opportunities as to the methods of recording gravitational waves and extracting signals, as well as the use of quantum non-perturbative measurements and the inclusion of gravitational antennas in the network.
The creation of new-generation gravitational antennas designed to reliably receive gravitational waves from remote space sources involves the use of high-power lasers, complex computer systems for processing large data arrays, the use of complex seismic protection systems, and the solution of other complex engineering and physical problems.
At present large international experience has been gained in the field of creating laser gravitational antennas, which ensured the ground-based recording of gravitational waves from black hole collision [5, 6] and neutron star merger . Furthermore, the gravitational wave propagation velocity was estimated, which appeared to be equal to the speed of light in vacuum with an accuracy of based on almost simultaneous recording of gravitational waves and a short gamma-ray burst from neutron star merger .
The modern theory of the early universe is based on the inclusion of the inflationary stage which precedes the stage of the hot universe. The theory of cosmological inflation  explains the origin of a large-scale structure and corresponds to observational data . Inflationary expansion of the universe during very early times, once the universe emerged from the quantum gravity (Planck) era, has been proposed in the late 1970s and mainly in the beginning of the 1980s and is becoming more accepted as a necessary stage which modifies the standard Big Bang theory model. According to the theory of inflation, the primordial perturbations appear from quantum fluctuations. These fluctuations had essential amplitudes in scales of Planck length, and during inflation they generate the primordial perturbations which then lead nearer to scales of galaxies with almost the same amplitudes .
Thus, the theory of cosmological inflation connects large-scale structure of the universe with microscopic scales. The resultant range of inhomogeneities practically doesn’t depend on scenarios of inflation and has a universal form. It leads to unambiguous predictions for a range of anisotropy of the background radiation.
The small quantum perturbations of the scalar field and the corresponding perturbations of the metric generate the relic gravitational waves. This type of gravitational waves was not directly observed; however, the possibility of such observations plays a key role to understand the physical processes in the early universe.
2. Inflationary stage of the early universe
The models of inflationary (accelerated) expansion of the universe at the early stage of its evolution, that is, at times close to the Planck time, are becoming increasingly convincing as a necessary step modifying the standard Big Bang theory, which is based on solutions of Einstein’s equations for the universe filled with ordinary baryonic matter with a positive energy density obtained by Friedmann. However, the extrapolation of Friedmann solutions to early times leads to many insoluble problems when constructing on their basis the evolution scenarios of the universe .
The exponential (de Sitter) expansion, suggesting , or a close expansion of the early universe based on the evolution of a certain substance with the equation of state , i.e., with a negative pressure, is a feature of inflation models which allow to solve the problems of the standard model of the Big Bang theory, namely, the problems of the horizon, flatness, homogeneity, isotropy, low concentration of exotic states of matter (domain walls, monopoles, etc.), anisotropy of the background radiation, the initial singularity, and some other problems .
Thus, the cosmological models containing a combination of Friedmann solutions and (quasi) de Sitter solutions provide the basis for a current description of the evolution of the universe. In the context of the inflationary paradigm, the early universe expands for some time accelerated and, further, goes into a power-law expansion mode without acceleration corresponding to Friedmann solutions.
In most cosmological models, the geometric description of the universe is based on the Friedmann-Robertson-Walker (FRW) homogeneous isotropic space (space-time) model, which is associated with a high degree of isotropy of space, measured on the basis of the cosmic microwave background (CMB) radiation research. This identification also relies on a formal result known as the Ehlers-Geren-Sachs theorem, which refers to the universe filled with any ideal barotropic fluid .
The metric of Friedmann-Robertson-Walker (FRW) space-time is written as follows:
where is a scale factor which characterizes the size of the universe, are the spherical coordinates, and the values , , and correspond to a spatially flat, closed, and open model of the universe.
The source of the accelerated expansion of the early universe with the equation of state is a vacuum; the equation of the state corresponds to a scalar (bosonic) field. The Bose-Einstein statistics for an ensemble of bosons, in contrast to an ensemble of fermions obeying the Pauli exclusion principle, implies that there can be several particles in one quantum state, which leads to the formation of boson condensate in which the increase in the concentration of massless bosons is associated with a decrease in the effective pressure corresponding to the equation of state . The initial (quasi) exponential expansion associated with a negative pressure, due to the exotic equation of state, is unstable, which leads to a phase transition, the termination of accelerated expansion, and the fragmentation of the original volume into many areas in which further evolution corresponds to the Friedmann solutions.
Also, the presence of a scalar field violates the symmetry of the system, which leads to the appearance of a mass of initially massless particles, for example, in the Higgs field .
Thus, the inclusion of the scalar field into cosmological models makes it possible to move from (quasi) de Sitter solutions to the Friedmann ones.
To prevent the rapid decay of the state of , it is necessary to assume the existence of some potential barrier, that is, the minimum potential energy of a scalar field. Consequently, in realistic inflation models, the scalar field evolves from a state of “false vacuum” with a non-zero potential energy to a state of “true vacuum,” corresponding to the minimum of a potential . In other words, the scalar field rolls down (or tunnels) from some initial state to the minimum of , and the nature of this process is determined by the shape of the potential.
At the moment, there are many models of cosmological inflation with different potentials of a scalar field and different specifics of its evolution. A large number of current models of the early universe on the basis of the inflationary paradigm are considered in the review .
The physical justification for the inclusion of scalar fields in cosmological models is based on the experimental detection of the Higgs boson in the experiment at the Large Hadron Collider . Thus, the scalar field corresponding to the Higgs bosons can be considered as the source of the gravitational field of the early universe. Moreover, the Higgs field can be considered as “inflation,” leading to early accelerated expansion of the universe.
Now, in the system of units , we write the action that determines the dynamics of a scalar field based on Einstein’s theory of gravity:
where is the Ricci scalar and is the potential of a scalar field.
where is the Hubble parameter, is the kinetic energy of a canonical scalar field , and the dot denotes the derivative with respect to cosmic time .
Also, the state parameter of a scalar field can be calculated as
To build a consistent model of cosmological inflation, the following conditions must be met:
The presence of the stage of accelerated expansion, which implies
The completion of the stage of accelerated expansion
The reheating of the scalar field with the subsequent formation of photons, i.e., the transition to the stage of predominance of radiation with
Currently, along with other models, several types of cosmological inflationary models are considered, which differ in both by the type of potential and the initial conditions under which an inflationary stage occurs: namely, a scalar field can be located at one of its potential minima, or accelerated expansion occurs for any conditions permitting the onset of inflation for scalar field energy density values comparable to the Planck mass .
The form of the scalar field’s potential is determined from the physics of elementary particles and theories of the unification of fundamental interactions, such as supersymmetric theories and string theories in the context of the inflationary paradigm. Physical mechanisms corresponding to a large number of inflationary potentials were discussed in the review . Due to the fact that the potential of the scalar field has a great importance for determining the physical processes at the stage of cosmological inflation, the potential is given to build models of the early universe.
However, the finding of exact solutions to the system of Eqs. (3)–(5) for a given potential is impossible in most cases due to their nonlinearity. For this reason, a convenient tool for analyzing inflationary models based on a given scalar field potential is the “slow-roll approximation” which implies that and and, therefore, simplifies the initial dynamic equations .
The dynamics of the expansion of the universe which determined by the scale factor is no less important when analyzing cosmological models. By setting the expansion law , it is often possible to find the exact solutions of the system of Eqs. (3)–(5) and restore the evolution of the scalar field and the potential . The different methods for constructing exact and approximate solutions of the equations of cosmological dynamics (3)–(5) can be found, for example, in the papers [15, 16]. We also note that the system of Eqs. (3)–(5) has many solutions that satisfy all the conditions for the inflationary stage that were outlined earlier.
Now, we consider the parameters that are necessary for the analysis of inflationary stage, namely, the e-fold number and the slow-roll parameters.
The e-fold number is usually noted as the natural logarithm of the ratio of the scale factor at the end of inflation to the scale factor at the beginning of inflation :
where and are the times of the beginning and ending of the inflation. The value of the number of e-folds at the end of the inflationary stage is estimated as .
When analyzing inflationary models, the slow-roll parameters are important, and these parameters are defined as follows :
3. Cosmological perturbations
Cosmological perturbations are the source of the evolution of large-scale structure of the universe. An explanation of the distribution of galaxies and clusters of galaxies at large distances in the observable part of the universe on the basis of cosmological perturbations was originally proposed in the works of Harrison  and Zeldovich . In the context of inflationary paradigm, the source of cosmological perturbations is quantum fluctuations of a scalar field and the corresponding fluctuations of the metric, which, in a linear order, correspond to three modes that evolve independently.
It is known from the classical theory of cosmological perturbations that the analysis of metric inhomogeneities can be simplified to the study of one perturbed quantity . Thus, the quantum theory of cosmological perturbations can be reduced to the quantum description of the fluctuations of a certain scalar field.
Since the background in which the scalar field evolves depends on time, the field mass will also depend on time. This dependence of the field’s mass on time will lead to the appearance of particles if the evolution begins with a certain vacuum state. Quantum particle production corresponds to the development and growth of cosmological perturbations.
In inflation models with one scalar field, at the crossing of the Hubble radius, cosmological perturbations “freeze,” and their quantum state begins to change in such a way that the condition of constant amplitude is satisfied. The freezing of the vacuum state leads to the appearance of the classical properties . Thus, the theory of cosmological perturbations provides a consistent approach for considering the generation and evolution of cosmological perturbations.
The influence of cosmological perturbations on the anisotropy and polarization of the background radiation is determined on the basis of spectral parameters and observational restrictions on the values of which form the basis of the experimental verification of theoretical models of the early universe. Also, within the framework of the cosmological perturbation theory, it is possible to calculate the spectra of initial density perturbations and relic gravitational waves depending on the values of the parameters of theoretical models .
After the end of the inflationary stage, the scalar field reheating and the formation of the first light particles of baryon matter begin. In the hot dense plasma, due to scattering on electrons, photons propagate much slower than the speed of light. When the universe expands so much that the plasma cools down to the recombination temperature, the electrons begin to connect with the protons, forming neutral hydrogen, and the photons begin to spread freely.
The points from which the photons reach the observer form the last scattering surface, whose temperature at the time of recombination is ∼3000 K and rapidly decreases with the expansion of the universe. The background radiation temperature is isotropic with an accuracy of . The low anisotropy manifests itself as the temperature difference in different directions and its value is approximately equal to 3 mK .
The kinetic component of the anisotropy of the cosmic background radiation is due to the movement of the observer relative to the background radiation, which corresponds to the dipole harmonic.
In addition to the kinetic component in the anisotropy of the CMB, there are potential terms associated with effects in gravitational fields of very large scales that are comparable to the distance to the last scattering surface, namely:
Sachs-Wolfe effect, which corresponds to a change in the photon energy in a variable gravitational field of the universe
Silk effect due to adiabatic compression of radiation and baryon acoustic oscillations prior to the recombination epoch in high- and low-density zones
In the zero order of the cosmological perturbation theory, the universe is described by a single function of time, namely, by the scale factor . In the first (linear) order, the perturbations of the metric are the sum of three independent modes—scalar, vector, and tensor (relic gravitational waves), each of which is characterized by the spectral function of the wave number .
For the inflationary stage in the linear approximation, one can write the Mukhanov-Sasaki equations for Fourier modes of the scalar and tensor perturbations :
where , is a wave number, and is the conformal time.
Eqs. (11) and (12) allow finding the power spectra and and spectral indices and of the scalar and tensor perturbations. The formulas for calculating the main cosmological parameters at crossing the Hubble radius  are
The data on the effects of scalar and tensor modes can be obtained from observations of the anisotropy and polarization of the cosmic microwave background (CMB) radiation, which arose as a result of the joint effect on the photon distribution of the perturbation modes. Observational restrictions on the values of the parameters of cosmological perturbations according to the data of the PLANCK are 
In the context of such verification of cosmological models, let us pay attention to the tendency for the upper limit to decrease by the value of the tensor-scalar ratio for updated observational data .
Also, we note that the relic gravitational waves were not directly observed, which leads to a large number of theoretical models of cosmological inflation, which provide an explanation of the origin and evolution of the large-scale structure of the universe and correspond to the observational constraints.
4. Generalized exponential power-law inflation
The scheme for constructing models of the early universe based on the evolution of the scalar field in the context of the inflationary paradigm can be represented as follows:
The generating solutions of background dynamic equations (excluding quantum fluctuations of the scalar field) for a given potential, the law of accelerated expansion of the early universe, or the evolution of a scalar field.
Analysis of the quantum fluctuations of a scalar field and the corresponding metric perturbations on the basis of the theory of cosmological perturbation for the previously obtained background solutions. The result of this analysis is the values of the spectral parameters of cosmological perturbations which can be calculated from Eqs. (13)–(17).
To build cosmological models corresponding to observational data, we propose the principle of constructing the inflationary models with generalized exponential power-law expansion. For this aim we consider any exact solutions of Eqs. (3)–(5) for which the substitution of the slow-roll parameters (8)–(10) into Eqs. (13)–(17) doesn’t correspond to observational constraints (17)–(18).
After the following transformations
one has new exact solutions with new slow-roll parameters
and with the conformity to observational constraints which can be achieved by choosing the values of free constant parameters and .
The proposed approach has two limitations:
The original scale factor doesn’t violate the law of accelerated expansion.
The potential corresponding to the scale factor (20) implies the evolution of the scalar field , according to the inflationary paradigm.
Transformations (21)–(24) define a class of models with the generalized exponential power-law dynamics, and the original scale factor may not correspond to the condition of accelerated expansion ; however, the resulting scale factor implies a combination of the de Sitter solution (for ) and the power-law expansion (for ), which corresponds to the basic feature of the inflationary paradigm implying a graceful exit from the stage of accelerated expansion to the power-law non-accelerated expansion.
5. Relic gravitational waves
As an additional verification tool for cosmological models, we consider the possibility of direct detection of the relic gravitational waves. The detection of relic gravitational waves is extremely important for determining the parameters of the models of early universe. Additionally, such a detection enhances the position of the inflationary paradigm compared to alternative scenarios, for example, the models with a rebound from singularity in which cosmological gravitational waves are absent .
As the main observational characteristic of relic gravitational waves, we consider the energy density, which is usually determined by the dimensionless quantity :
where is the linear frequency, is the critical energy density, is the value of the Hubble parameter in the modern era, and is the energy density of gravitational waves.
Also, the energy density of relic gravitational waves can be represented in terms of the power spectrum:
The frequency and energy density of relic gravitational waves are limited by the following conditions :
The energy density of relic gravitational waves should not exceed
The temperature of the scalar field and the frequency of gravitational waves at the end of the inflation stage are
where is the effective number of the degrees of freedom (in standard model of elementary particles ).
The application to the analysis of the models of the early universe only of the slow-roll approximation implies a low-frequency spectrum of relic gravitational waves in the range of Hz . However, the predominance of the kinetic energy of the scalar field during the evolution of the early universe provides a theoretical justification for the existence of high-frequency relic gravitational waves in models with one scalar field in the range of Hz  which can be used as affordable means of verification of models of the early universe in the presence of physical effects that increase the sensitivity of the detector to the required level.
Currently, the most productive method of direct detection of gravitational waves is the use of interferometers as detectors, which was proposed in the article by Gertsenshtein and Pustovoit . This principle is widely used in modern laser interference gravitational antennas, the main element of which is the Fabry-Perot interferometer. These are broadband gravitational antennas, which offer a lot of opportunities as to the methods of recording of gravitational waves and extracting signals, as well as the use of quantum non-perturbative measurements and the inclusion of gravitational antennas in the network. The main element of laser interference gravitational antennas, as a rule, is Fabry-Perot multipath free-mass resonator, on whose properties the sensitivity and noise immunity of the entire gravitational antenna largely depend [4, 23, 24].
After creating the first laser interferometer for detecting gravitational waves, systematic work began on the creation and improvement of such devices in various laboratories around the world. The experience of gravitational antenna projects by VIRGO (Italy, France), LIGO (USA), TAMA (Japan), CLIO (Japan), GEO-600 (Germany), and OGRAN (Russia) will certainly be used to create more compact and highly sensitive antennas of new generation . Also, as the most promising project for the direct detection of gravitational waves, work on the creation of a space interferometer in a helio-stationary orbit should be noted, in which the distance between the mirrors will be about 1 million kilometers. This project is called Laser Interferometer Space Antenna (LISA) . The implementation of the LISA project is scheduled for 2029.
One of the promising methods for increasing the sensitivity of gravitational antennas in the high-frequency region of the spectrum is to use the phenomenon of low-frequency optical resonance, which distinguishes this approach from other projects on the detection of gravitational waves. The presence of this effect in Fabry-Perot interferometers was first considered in [23, 24]. At the moment, there is a high-frequency gravitational wave detector, which was built at the University of Birmingham, United Kingdom . Also, it is planned to build the high-frequency gravitational wave detectors in Japan .
Thus, at the moment there are a large number of promising methods for direct observation of gravitational waves, which correspond to the ability to measure the characteristics of relic gravitational waves for a better understanding of the physical processes occurring in the early universe.
We considered the basis of building and verifying of the inflationary models of early universe. As the method for constructing the exact cosmological solutions corresponding to observational constraints, the models with generalized exponential power-law dynamics are proposed.
The verification of the relevance of such models is related to the estimation of the contribution of relic gravitational waves to the anisotropy and polarization of the cosmic microwave background radiation. Therefore, there are a lot of inflationary models with different scalar field potentials that will satisfy the observational constraints.
The most obvious way to significantly reduce the number of theoretical models of cosmological inflation is direct detection of relic gravitational waves.
The most promising methods in this area of experimental research are using the interferometers as detectors. The interesting direction of the observation is the detection of high-frequency relic gravitational waves using the effect of low-frequency optical resonance proposed in [23, 24].