# References

BEYOND THE SIMPLEST INFLATIONARY COSMOLOGICAL MODELS

A.A. Starobinsky

Landau Institute for Theoretical Physics,

Russian Academy of Sciences, Moscow 117334, Russia

Though predictions of the simplest inflationary cosmological models with cold dark matter, flat space and approximately flat initial spectrum of adiabatic perturbations are remarkably close to observational data, we have to go beyond them and to introduce new physics not yet discovered in laboratories to account for all data. Two extensions of these models which seem to be the most actual at present time are discussed. The first one is the possibility that we are living at the beginning of a new inflation-like era. Then classical cosmological tests, like the luminosity distance or the angular size of distant objects as functions of redshift, as well as the behaviour of density perturbations in a dustlike matter component including baryons as a function of redshift, can provide information sufficient for the unambiguous determination of an effective potential of a corresponding present inflaton scalar field. The second, unrelated extension is a possibility of broken-scale-invariant cosmological models which have localized steps or spikes in the primordial perturbation spectrum. These features can be produced by fast phase transitions in physical fields other than an inflaton field in the early Universe during inflation and not far from the end of it. At present, it seems that the only scale in the spectrum around which we might see something of this type is $`k=0.05h`$ Mpc<sup>-1</sup>.

1. Introduction

It is clear now that viable cosmological models of the present Universe cannot be constructed without cold non-baryonic dark matter (CDM) and some model of the early Universe producing an approximately flat (or, Harrison-Zeldovich-like, $`n_S1`$) initial spectrum of scalar (adiabatic) perturbations. The two necessary ingredients correspond to physics not yet discovered in laboratories. In particular, since the simplest cosmic defect models without inflation have failed to fit observational data on the large-scale structure and angular anisotropies of the cosmic microwave background (CMB) temperature in the Universe, the only remaining viable and well elaborated scenario producing such a spectrum is the inflationary scenario with one effective scalar field (inflaton). The existence of the inflaton field with desired properties of its self-interaction potential is a hypothesis about physics at very high energies, only few orders of magnitude below the Planck scale (in addition to this hypothesis about fundamental physics, we have to assume that our Universe did really pass through an inflationary stage in the past). Two main candidates for non-baryonic CDM - the lightest supersymmetric particle or the axion - are still awaiting a direct discovery of any (or, even both) of them.

However, it has been already understood for several years that this amount of new physics is not sufficient for cosmology. To explain all existing observational data, one has to go further and to introduce even more new physics. Really, it is well established that the simplest Friedmann-Robertson-Walker (FRW) cosmological model with CDM, the flat 3-space ($`\mathrm{\Omega }_0=1`$) and $`n_S=1`$ does not fit observational data. Looking at viable extensions of this model presented in the Table 2 of the author’s plenary talk at the ”Cosmion-94” conference , it is clear that all of them require new physics. It enters either at very low energies – producing an effective cosmological constant (the $`\mathrm{\Lambda }`$CDM model), or at medium energies – giving a non-zero rest mass to one or several neutrino types (CHDM models), or, finally, at very high energies – resulting either in phase transitions during single (one-scalar-field) inflation (a second-order phase transition in the case of broken-scale-invariant (BSI) CDM models, or a first-order transition in the case of open CDM (OCDM) models), or in the appearance of a second inflaton scalar field (double inflation) that also lead to BSI CDM, or OCDM depending on parameters.

After three years passed from the ”Cosmion-94” conference, all these four classes of models are still remaining viable in the sense that neither of them may be excluded completely. However, observational evidence (especially, the most recent one) becomes more and more directed towards the $`\mathrm{\Lambda }`$CDM scale-invariant model. That is why I consider it as the model No. 1 at present, and I shall focus my main attention on it in Sec. 2. Still there exist some data, however inconclusive they might be, indicating that the present matter perturbation power spectrum has some feature at $`k0.05h`$ Mpc<sup>-1</sup>. If so, then it can be verified that this feature is present in the initial perturbation spectrum, too. It means that we are dealing with a BSI CDM model. This possibility does not remove need in a cosmological constant, so actually a BSI $`\mathrm{\Lambda }`$CDM model is required for a better fit to data. This topic is considered in Sec. 3.

2. Cosmological models with a variable cosmological term

It has been known for many years that the flat FRW cosmological model with cold dark matter, a positive cosmological constant $`\mathrm{\Lambda }>0`$ ($`\mathrm{\Omega }_0+\mathrm{\Omega }_\mathrm{\Lambda }=1`$) and an approximately flat spectrum of primordial adiabatic perturbations fits observational data better and has a larger admissible region of the parameters $`(H_0,\mathrm{\Omega }_0)`$ than any other cosmological model with both inflationary and non-inflationary initial conditions (see, e.g., ). Here $`H_0`$ is the Hubble constant, $`\mathrm{\Omega }_0=8\pi G\rho _m/3H_0^2`$ includes baryons and (mainly) non-baryonic dark matter, $`\mathrm{\Omega }_\mathrm{\Lambda }\mathrm{\Lambda }/3H_0^2`$ and $`c=\mathrm{}=1`$. This conclusion was based on the following arguments: a) relation between $`H_0`$ and the age of the Universe $`t_0`$, b) the fact that observed mass/luminosity ratio never leads to values more than $`\mathrm{\Omega }_00.4`$ up to supercluster scales, c) comparison of cosmic microwave background temperature anisotropies, power spectra of density and velocity matter perturbations, present abundance of galaxy clusters with predictions of cosmological models with inflationary initial conditions; d) observed values of $`\rho _b/\rho _m`$ in rich galaxy clusters confronted with the range for the present baryon density $`\rho _b`$ admitted by the theory of primordial (Big Bang) nucleosynthesis. I don’t include gravitational lensing tests (e.g., a number of lensed quasars) here, since conclusions based on them are less definite at present.

During last year two new pieces of strong evidence for $`\mathrm{\Omega }_0<1`$ have appeared. The first (historically) of them is based on the evolution of abundance of rich galaxy clusters with redshift $`z`$ . Still, it should be noted that there have been already appeared some doubts on validity of the conclusion that $`\mathrm{\Omega }_0=1`$ is really excluded . Much better observational data expected in near future will help to resolve this dilemma unambiguously. The second, completely independent argument for $`\mathrm{\Omega }_0=(0.20.4)`$ follows from direct observations of supernovae (type Ia) explosions at high redshifts up to $`z1`$ . On the other hand, no direct evidence for a negative spatial curvature of the Universe (i.e., for the OCDM model) has been found. Of course, the possibility to have both a positive cosmological constant and spatial curvature of any sign is not yet excluded, but, according to the ”Okkam’s razor” principle, it would be desirable not to introduce one more basic novel feature of the Universe (spatial curvature) without conclusive observational evidence. In any case, in spite of many theoretical and experimental attempts to exorcize it, a $`\mathrm{\Lambda }`$-term is back again.

It is clear that the introduction of a cosmological constant requires new and completely unknown physics in the region of ultra-low energies. Solutions with a cosmological constant occur in such fundamental theories as supergravity and M-theory. However, this cosmological constant is always negative and very large. As compared to such a basic ”vacuum” state, a very small and positive cosmological constant allowed in the present Universe may be thought as corresponding to the energy density $`\epsilon _\mathrm{\Lambda }`$ of a highly excited (though still very symmetric) ”background” state. So, it need not be very ”fundamental”. But then it is natural to omit the assumption that it should be exactly constant. In this case the name ”a cosmological term” (or a $`\mathrm{\Lambda }`$-term) is more relevant for it, so I shall use this one below. The principal difference between two kinds of non-baryonic dark matter - dustlike CDM and a $`\mathrm{\Lambda }`$-term - is that the latter one is not gravitationally clustered up to scales $`30h^1`$ or more (otherwise we would return to the problem why $`\mathrm{\Omega }_0`$ observed from gravitational clustering is not equal to unity). Here $`h=H_0/100`$ km s<sup>-1</sup>Mpc<sup>-1</sup>.

On the other hand, there exists a well-known strong argument showing that a $`\mathrm{\Lambda }`$-term cannot change with time as fast as the matter density $`\rho _m`$ and the Ricci tensor (i.e., $`t^2`$) during the matter dominated stage (for redshifs $`z<410^4h^2`$). Really, if $`\epsilon _\mathrm{\Lambda }\rho _m`$, so that $`\mathrm{\Omega }_\mathrm{\Lambda }=const`$, then matter density perturbations in the CDM+baryon component grow as $`\delta \left(\frac{\delta \rho }{\rho }\right)_mt^\alpha (1+z)^{3\alpha /2},\alpha =\frac{\sqrt{2524\mathrm{\Omega }_\mathrm{\Lambda }}1}{6}`$. As a consequence, the total growth of perturbations $`\mathrm{\Delta }`$ since the time of equality of matter and radiation energy densities up to the present moment is less than in the absence of the $`\mathrm{\Lambda }`$-term. If $`\mathrm{\Omega }_\mathrm{\Lambda }1`$, then $`\mathrm{\Delta }(\mathrm{\Omega }_\mathrm{\Lambda })=\mathrm{\Delta }(0)(1(6.4+2\mathrm{ln}h)\mathrm{\Omega }_\mathrm{\Lambda })`$. Since parameters of viable cosmological models are so tightly constrained that $`\mathrm{\Delta }`$ may not be reduced by more than twice approximately, this type of a $`\mathrm{\Lambda }`$-term cannot account for more than $`0.1`$ of the critical energy density (see  for detailed investigation confirming this conclusion). This, unfortunately, prevents us from natural explanation of the present $`\mathrm{\Lambda }`$-term with $`\mathrm{\Omega }_\mathrm{\Lambda }=(0.50.8)`$ using ”compensation” mechanisms  or exponential potentials with sufficiently large exponents .

A natural and simple description of a variable $`\mathrm{\Lambda }`$-term is just that which was so successively used to construct the simplest versions of the inflationary scenario, namely, a scalar field with some interaction potential $`V(\phi )`$ minimally coupled to the Einstein gravity. Such an approach, though phenomenological, is nevertheless more consistent and fundamental than a commonly used attempt to decribe a $`\mathrm{\Lambda }`$-term by a barotropic ideal fluid with some equation of state. The latter approach cannot be made internally consistent in case of negative pressure. In particular, it generally leads to imaginary values of the sound velocity. On the contrary, no such problems arise using the scalar field description (this scalar field is called the $`\mathrm{\Lambda }`$-field below). Of course, its effective mass $`|m_\phi ^2|=|d^2V/d\phi ^2|`$ should be very small to avoid gravitational clustering of this field in galaxies, clusters and superclusters. To make a $`\mathrm{\Lambda }`$-term slowly varying, we assume that $`|m_\phi |H_010^{33}`$ eV, or less (though this condition may be relaxed). Models with a time-dependent $`\mathrm{\Lambda }`$-term were introduced more than ten years ago , and different potentials $`V(\phi )`$ (all inspired by inflationary models) were considered: exponential , inverse power-law , power-law , cosine .

However, it is clear that since we know essentially nothing about physics at such energies, there exists no preferred theoretical candidate for $`V(\phi )`$. In this case, it is more natural to go from observations to theory, and to determine an effective phenomenological potential $`V(\phi )`$ from observational data. The two new tests mentioned above are just the most suitable for this aim. Really, using the cluster abundance $`n(z)`$ determined from observations and assuming the Gaussian statistics of initial perturbations (the latter follows from the paradigm of one-field inflation, and it is in agreement with other observational data), it is possible to determine a linear density perturbation in the CDM+baryon dustlike component $`\delta (z)`$ for a fixed comoving scale $`R8(1+z)^1h^1`$ Mpc up to $`z1`$, either using the Press-Schechter approximation, or by direct numerical simulations of nonlinear gravitational instability in the expanding Universe. $`\delta (z)`$ can be also determined from observation of gravitational clustering (in particular, of the galaxy-galaxy correlation function) as a function of $`z`$. On the other hand, observations of SNe at different $`z`$ yield the luminosity distance $`D_L(z)`$ through the standard astronomical expression $`m=M+5\mathrm{log}D_L+25`$, where $`m`$ is the observed magnitude, $`M`$ is the absolute magnitude and $`D_L`$ is measured in Mpc. One more related cosmological test is the angular size $`\theta (z)`$ of an extended object with a proper size $`d`$. This test is less accurate at present due to absence of objects with their size being fixed to a desired accuracy, still it is important.

Let us show now how to determine $`V(\phi )`$ from either $`\delta (z)`$, or $`D_L(z)`$, or $`\theta (z)`$. Also, we investigate what additional information is necessary for an unambiguous solution of this problem in all these cases. The derivation of $`V(\phi )`$ consists of two steps. First, the Hubble parameter $`H\frac{\dot{a}}{a}=H(z)`$ is determined. Here $`a(t)`$ is the FRW scale factor, $`1+za_0/a`$, the dot means $`\frac{d}{dt}`$ and the index $`0`$ denotes the present value of a corresponding quantity (in particular, $`H(t_0)=H(z=0)=H_0`$). In the case of SNe, the first step is almost trivial since the textbook expression for $`D_L`$ reads:

$$D_L(z)=a_0(\eta _0\eta )(1+z),\eta =_0^t\frac{dt}{a(t)}.$$
(1)   
Therefore,

$$H(z)=\frac{da}{a^2d\eta }=(a_0\eta ^{})^1=\left[\left(\frac{D_L(z)}{1+z}\right)^{}\right]^1.$$
(2)   
Here and below in Sec. 2, a prime denotes the derivative with respect to $`z`$.

The angular size is given by

$$\theta (z)=\frac{d}{a(\eta )(\eta _0\eta )}=\frac{d(1+z)}{a_0(\eta _0\eta )}.$$
(3)   
Thus,

$$H(z)=(a_0\eta ^{})^1=\left[d\left(\frac{1+z}{\theta (z)}\right)^{}\right]^1.$$
(4)   
In physical units, the right-hand sides of Eqs. (24) should be multipled by $`c`$. Thus, both $`D_L(z)`$ and $`\theta (z)`$ define $`H(z)`$ uniquely.

More calculations are required to find $`H(z)`$ from $`\delta (z)`$. The system of background equations for the system under consideration is:

$$H^2=\frac{8\pi G}{3}\left(\rho _m+\frac{\dot{\phi }^2}{2}+V\right),\rho _m=\frac{3\mathrm{\Omega }_0H_0^2a_0^3}{8\pi Ga^3},$$
(5)   
$$\ddot{\phi }+3H\dot{\phi }+\frac{dV}{d\phi }=0,$$
(6)   
$$\dot{H}=4\pi G(\rho _m+\dot{\phi }^2).$$
(7)   
Eq. (7) is actually the consequence of the other two equations.

We consider a perturbed FRW background which metric, in the longitudinal gauge (LG), has the form:

$$ds^2=(1+2\mathrm{\Phi })dt^2a^2(t)(1+2\mathrm{\Psi })\delta _{lm}dx^ldx^m,l,m=1,2,3.$$
(8)   
The system of equations for scalar perturbations reads (the spatial dependence $`\mathrm{exp}(ik_lx^l)`$, $`k_lk^lk^2`$ is assumed):

$$\mathrm{\Phi }=\mathrm{\Psi }=\dot{v},\dot{\delta }=\frac{k^2}{a^2}v+3(\ddot{v}+H\dot{v}+\dot{H}v),$$
(9)   
$$\dot{\mathrm{\Phi }}+H\mathrm{\Phi }=4\pi G(\rho _mv+\dot{\phi }\delta \phi ),$$
(10)   
$$\left(\frac{k^2}{a^2}+4\pi G\dot{\phi }^2\right)\mathrm{\Phi }=4\pi G(\rho _m\delta +\dot{\phi }\dot{\delta \phi }+3H\dot{\phi }\delta \phi +\frac{dV}{d\phi }\delta \phi ),$$
(11)   
$$\ddot{\delta \phi }+3H\dot{\delta \phi }+\left(\frac{k^2}{a^2}+\frac{d^2V}{d\phi ^2}\right)\delta \phi =4\dot{\phi }\dot{\mathrm{\Phi }}2\frac{dV}{d\phi }\mathrm{\Phi }.$$
(12)   
Eq. (12) is the consequence of other ones. Here $`v`$ and $`\delta \phi `$ are, correspondingly, a velocity potential of a dustlike matter peculiar velocity and a $`\mathrm{\Lambda }`$-field perturbation in LG, and $`\delta `$ is a comoving fractional matter density perturbation (in this case, it coincides with $`\left(\frac{\delta \rho }{\rho }\right)_m`$ in the synchronous gauge). In fact, all these perturbed quantities are gauge-invariant.

Now let us take a comoving wavelength $`\lambda =k/a(t)`$ which is much smaller than the Hubble radius $`H^1(t)`$ up to redshifts $`z5`$. This corresponds to $`\lambda 2000h^1`$ Mpc at present. Then, from Eq. (12),

$$\delta \phi \frac{a^2}{k^2}(4\dot{\phi }\dot{\mathrm{\Phi }}2\frac{dV}{d\phi }\mathrm{\Phi }),|\dot{\phi }\dot{\delta \phi }||\frac{dV}{d\phi }\delta \phi |\frac{a^2H^4}{Gk^2}|\mathrm{\Phi }|\rho _m|\delta |.$$
(13)   
It is at this place where we heavily used the condition $`|m_\varphi |\stackrel{<}{}H_0`$. In particular, it gives us a possibility to estimate $`|\dot{\mathrm{\Phi }}|`$ as $`H|\mathrm{\Phi }|`$ and $`|\dot{\delta \phi }|`$ as $`H|\delta \phi |`$. On the other hand, it is clear that this condition may be somewhat relaxed without destroying the last inequality in Eq. (13). As a result, the $`\mathrm{\Lambda }`$-field is practically unclustered at the scale involved. Now the last of Eqs. (9) and Eq. (11) may be simplified to:

$$\dot{\delta }=\frac{k^2}{a^2}v,\frac{k^2}{a^2}\mathrm{\Phi }=4\pi G\rho _m\delta .$$
(14)   
Combining this with the first of Eqs. (9), we return to a well-known equation for $`\delta `$ in the absence of the $`\mathrm{\Lambda }`$-field:

$$\ddot{\delta }+2H\dot{\delta }4\pi G\rho _m\delta =0.$$
(15)   
It is not possible to solve this equation analytically for an arbitrary $`V(\phi )`$. Remarkably, the inverse dynamical problem, i.e. the determination of $`H(a)`$ given $`\delta (a)`$, is solvable. After changing the argument in Eq. (15) from $`t`$ to $`a`$ ($`\frac{d}{dt}=aH\frac{d}{da}`$), we get a first order linear differential equation for $`H^2(a)`$:

$$a^2\frac{d\delta }{da}\frac{dH^2}{da}+2\left(a^2\frac{d^2\delta }{da^2}+3a\frac{d\delta }{da}\right)H^2=\frac{3\mathrm{\Omega }_0H_0^2a_0^3\delta }{a^3}.$$
(16)   
The solution is:

$$H^2=\frac{3\mathrm{\Omega }_0H_0^2a_0^3}{a^6}\left(\frac{d\delta }{da}\right)^2_0^aa\delta \frac{d\delta }{da}𝑑a=3\mathrm{\Omega }_0H_0^2\frac{(1+z)^2}{\delta ^2}_z^{\mathrm{}}\frac{\delta |\delta ^{}|}{1+z}𝑑z.$$
(17)   
Putting $`z=0`$ in this expression for $`H`$, we arrive to the expression of $`\mathrm{\Omega }_0`$ through $`\delta (z)`$:

$$\mathrm{\Omega }_0=\delta ^2(0)\left(3_0^{\mathrm{}}\frac{\delta |\delta ^{}|}{1+z}𝑑z\right)^1.$$
(18)   
Of course, observations of gravitational clustering can hardly provide the function $`\delta (z)`$ for too large $`z`$ (say, for $`z>5`$). However, $`\delta (z)`$ in the integrands in Eqs. (17,18) may be well approximated by its $`\mathrm{\Omega }_0=1`$ behaviour (i.e., $`\delta (1+z)^1`$) already for $`z>(23)`$. If massive neutrinos are present, one should use here the expression with $`\alpha `$ written above and with $`\mathrm{\Omega }_\mathrm{\Lambda }`$ substituted by $`\mathrm{\Omega }_\nu /\mathrm{\Omega }_0`$ (it is assumed that $`\rho _m`$ includes massive neutrinos, too).

Using Eq. (18), Eq. (17) can be represented in a more convenient form:

$$\frac{H^2(z)}{H^2(0)}=\frac{(1+z)^2\delta ^2(0)}{\delta ^2(z)}3\mathrm{\Omega }_0\frac{(1+z)^2}{\delta ^2(z)}_0^z\frac{\delta |\delta ^{}|}{1+z}𝑑z.$$
(19)   
Thus, $`\delta (z)`$ uniquely defines the ratio $`H(z)/H_0`$. Of course, appearance of derivatives of $`\delta (z)`$ in these formulas shows that sufficiently clean data are necessary, but one may expect that such data will soon appear. Let us remind also that, for $`\mathrm{\Lambda }const`$ ($`V(\phi )const`$), we have

$$H^2(z)=H_0^2(1\mathrm{\Omega }_0+\mathrm{\Omega }_0(1+z)^3),q_01+\left(\frac{d\mathrm{ln}H}{d\mathrm{ln}(1+z)}\right)_{z=0}=\frac{3}{2}\mathrm{\Omega }_01,$$
(20)   
where $`q_0`$ is the acceleration parameter.

The second step - the derivation of $`V(\phi )`$ from $`H(a)`$ \- is very simple. One has to rewrite Eqs. (57) in terms of $`a`$ and take their linear combinations:

$`8\pi GV(\phi )=aH{\displaystyle \frac{dH}{da}}+3H^2{\displaystyle \frac{3}{2}}\mathrm{\Omega }_0H_0^2\left({\displaystyle \frac{a_0}{a}}\right)^3,`$   
$`4\pi Ga^2H^2\left({\displaystyle \frac{d\phi }{da}}\right)^2=aH{\displaystyle \frac{dH}{da}}{\displaystyle \frac{3}{2}}\mathrm{\Omega }_0H_0^2\left({\displaystyle \frac{a_0}{a}}\right)^3,`$ (21)   
and then exclude $`a`$ from these equations, since the second of Eqs. (21) is integrated trivially.

To show explicitly how the second step proceeds for a given $`H(z)`$, let us consider the case when a cosmological term mimics the negative spatial curvature dynamically, while the real spatial geometry of the Universe remains flat. That is, we assume that

$$H^2(z)=H_0^2\left(\mathrm{\Omega }_0(1+z)^3+(1\mathrm{\Omega }_0)(1+z)^2\right).$$
(22)   
Integration of the second of Eqs. (21) gives

$$\frac{a}{a_0}=\frac{\mathrm{\Omega }_0}{1\mathrm{\Omega }_0}\mathrm{sinh}^2\left(\sqrt{\pi G}(\phi \phi _0+\phi _1)\right),\phi _1=\frac{1}{\sqrt{\pi G}}\mathrm{ln}\left(\frac{\sqrt{1\mathrm{\Omega }_0}+1}{\sqrt{\mathrm{\Omega }_0}}\right).$$
(23)   
Substituting this into the first of Eqs. (21), we find the expression for the effective potential:

$$V(\phi )=\frac{(1\mathrm{\Omega }_0)^3H_0^2}{4\pi G\mathrm{\Omega }_0^2}\frac{1}{\mathrm{sinh}^4\left(\sqrt{\pi G}(\phi \phi _0+\phi _1)\right)}.$$
(24)   
At early times at the matter dominated stage, this potential diverges $`(\phi \phi _0+\phi _1)^4`$.

Therefore, the model of a $`\mathrm{\Lambda }`$-term considered in this paper can account for any observed forms of $`D_L(z)`$, $`\theta (z)`$ and $`\delta (z)`$ which, in turn, can be transformed into a corresponding effective potential $`V(\phi )`$ of the $`\mathrm{\Lambda }`$-field. The only condition is that the functions $`H(z)`$ obtained by these three independent ways should coincide within observational errors. $`D_L(z)`$ and $`\theta (z)`$ uniquely determine $`V(\phi )`$, if $`\mathrm{\Omega }_0`$ is given additionally (the latter is required at the second step, in Eqs. (21)). $`\delta (z)`$ uniquely determines $`V(\phi )`$ up to the multiplier $`H_0^2`$, the latter has to be given additionally to fix an overall amplitude.

Observational tests which can falsify this model do exist. In particular, a contribution to large-angle $`\frac{\mathrm{\Delta }T}{T}`$ CMB temperature anisotropy due to the integrated (or, non-local) Sachs-Wolfe effect presents a possibility to distinguish the model from more complicated models, e.g., with non-minimal coupling of the $`\mathrm{\Lambda }`$-field to gravity or to CDM. However, the latter test is not an easy one, since this contribution is rather small and partially masked by cosmic variance.

Still it may appear finally that the cosmological term is really constant: $`V(\phi )=V_0=\epsilon _\mathrm{\Lambda }=const,\phi =\phi _0=const`$. Then it represents a new fundamental constant of nature. A question is often asked: how is it possible to obtain such a small value of $`\epsilon _\mathrm{\Lambda }`$ in Planck units from known physical constants? This question is actually not physical, but an arithmetic one. So, without trying to construct an underlying physical model for the $`\mathrm{\Lambda }`$-term, let me propose a toy arithmetic expression for $`\epsilon _\mathrm{\Lambda }`$ (of course, one of infinite number of possible ones) to show that there is no problem here from the mathematical point of view:

$$\epsilon _\mathrm{\Lambda }=\frac{M_P^4}{(2\pi ^2)^3}\mathrm{exp}\left(\frac{2}{\alpha }\right)=e^{283.02}M_P^4=10^{122.91}M_P^4,$$
(25)   
where $`M_P=\sqrt{G}`$ is the Planck mass and $`\alpha `$ is the fine-structure constant. In usual units, $`\epsilon _\mathrm{\Lambda }=10^{123}c^5/G^2\mathrm{}`$. Eq. (25) leads to the ”prediction”:

$$\mathrm{\Omega }_\mathrm{\Lambda }h^2(1\mathrm{\Omega }_0)h^20.335.$$
(26)   
In particular, for $`\mathrm{\Omega }_0=0.3`$, it gives $`H_0=69`$ km s<sup>-1</sup>Mpc<sup>-1</sup> – just in the expected region!

3. A feature in the matter power spectrum

and physical mechanisms to produce it

The main importance of the inflationary scenario of the early Universe for the theory of large scale structure in the present Universe is that the former scenario predicts (in its simplest realizations) an approximately flat, or scale-invariant, spectrum ($`n_S(k)d\mathrm{ln}P_0(k)/d\mathrm{ln}k1`$) of initial adiabatic perturbations. By the simplest realizations I mean, as usually, inflationary models with one effective slow-rolling scalar (inflaton) field. Of course, the physical nature of the inflaton may be completely different in these models, but it does not matter for observations, in particular, for the large scale structure. This prediction has been confirmed already, if by $`n_S(k)`$ we understand the slope of an initial power spectrum of density perturbations $`P_0(k)`$ smoothed over the range $`\mathrm{\Delta }\mathrm{ln}k1`$. To see this, it is even not neccessary to use results for $`n_S`$ following from the COBE experiment (though they also tell us the same), it is sufficient to compare the COBE normalization of perturbations for scales of the order of the present cosmological horizon $`R_h`$ with the $`\sigma _8`$ normalization that follows, e.g., from the present cluster abundance. The difference in the amplitude of initial perturbations at these two scales which are divided by approximately $`3`$ orders of magnitude is only $`22.5`$ times for the standard CDM model and even less for other models, e.g., the $`\mathrm{\Lambda }`$CDM model. In addition, these numbers give us an idea about the magnitude of expected deviations from the exact scale invariance: $`|n_S1|0.3`$ (once more, we are speaking about a smoothed $`n_S`$). Observational effects related to the part of $`P_0(k)`$ between these two points (CMB temperature fluctuations at medium and small angles, galaxy-galaxy and cluster-cluster correlations, peculiar velocities of galaxies) also do not require larger smooth deviations from $`n_S=1`$.

The observational fact that the smoothed slope $`n_S`$ cannot be significantly different from 1 does not exclude the possibility of local strong deviations from the flat spectrum, i.e., steps and/or spikes in $`P_0(k)`$. Of course, one should not expect such a behaviour to be typical, we shall see below that if it happens at all, it occurs at some preferred scales which themselves become new fundamental parameters of a cosmological model. Do we have any observational evidence for an existence of such preferred scales at the Universe? At present, only one scale in the Fourier space, $`k=k_0=0.05h`$ Mpc<sup>-1</sup>, remains a candidate for this role, and it seems that the spectrum is smooth for larger $`k`$ (from galaxy-galaxy correlation data) and for smaller $`k`$ (from CMB data). As for this scale itself, there exists an evidence for a peculiar behaviour (in the form of a sharp peak) in the Fourier power spectrum of rich Abell - ACO clusters (with richness class $`R0`$ and redshifts $`z0.12`$) around it  (see also ). This anomalous behaviour persists if the distant border for the cluster sample is reduced to $`z=0.070.08`$ . If we assume that the cluster power spectrum is proportional to the power spectrum of whole matter in the Universe (with some constant biasing factor), and calculate corresponding rms multipole values $`C_l`$ of angular fluctuations of the CMB temperature, they appear to be in a good agreement with existing results of medium-angle experiments . Moreover, if $`\mathrm{\Omega }_m=1`$, the peak in the power spectrum inferred from the cluster data just explains an excess in $`C_l`$ for $`l=200300`$ observed in the Saskatoon experiment. The agreement becomes even better in the case of the $`\mathrm{\Lambda }`$CDM model .

On the other hand, there is no peak at $`k=k_0`$ in the power spectrum of both APM clusters (which are generally less rich than Abell-ACO clusters) and APM galaxies  (though some less prominent feature at this scale may still exist in the latter spectrum ), and the maximum in these spectra seems to be shifted to $`k0.03h`$ Mpc<sup>-1</sup>. Leaving a solution of this discrepancy to more complete future surveys, let us consider theoretical predictions.

It is possible to produce local features in the initial spectrum even remaining (at least, formally) inside the standard paradigm of one-field inflation. The only thing which should be relaxed is the requirement of the analyticity of an inflaton effective potential $`V(\varphi )`$ at all points (the $`\varphi `$ inflaton field in the early Universe should not be confused with the $`\phi `$ $`\mathrm{\Lambda }`$-field used to describe a variable cosmological term in the previous section).

So, let me admit that $`V(\varphi )`$ has some kind of discontinuity at a point $`\varphi =\varphi _0`$. Of course, really this discontinuity is smoothed in a very small vicinity of $`\varphi _0`$. Three cases are the most interesting.

1. $`[V]=[V^{}]=0,[V^{\prime \prime }]0`$ at $`\varphi =\varphi _0`$.

Here $`[]`$ means the jump in the quantity considered, namely, $`[A]A(\varphi _0+0)A(\varphi _00)`$, and the prime denotes the derivative with respect to $`\varphi `$. If we assume that the slow-roll conditions $`V^248\pi GV^2,|V^{\prime \prime }|24\pi GV`$ are satisfied near the point $`\varphi =\varphi _0`$, then in the zero-order approximation the standard result for a perturbation spectrum is valid:

$$P_0(k)=\frac{k^4R_h^4(t)}{400}h^2(k),k^3h^2(k)=18\left(\frac{H^6}{V^2}\right)_k,H\frac{\dot{a}}{a}\sqrt{\frac{8\pi GV}{3}},$$
(27)   
where the index $`k`$ means that the quantity is evaluated at the moment of the first Hubble radius crossing ($`k=aH`$) at the inflationary stage. The result for $`P_0(k)`$, in contrast to the metric perturbation $`h^2(k)`$ defined in the ultra-synchronous gauge ($`h`$ is equal to $`1/3`$ of the trace of a spatial metric perturbation in this gauge), refers to the matter-dominated stage where $`R_h(t)=2/H=3t`$. Note also that there is no necessity in adding the multiplier $`𝒪(1)`$ here.

So, in this case $`P_0(k)`$ is continuous at $`k=k_0`$ but its slope $`n_S(k)`$ has a step-like behaviour there (similar to the case considered in ). However, due to small corrections to Eq. (27) which are beyond the slow-roll approximation, it appears that $`n_S`$ cannot be obtained simply by differentiating (27), and I expect that the sharp behaviour in $`n_S`$ will be smoothed near $`k_0`$. This question is still under consideration.

2. $`[V]=0,[V^{}]0`$ at $`\varphi =\varphi _0`$.

Now the second of the slow-roll conditions is violated, while we can choose parameters of the jump in such a way that the first condition is still valid. Naive application of Eq.(27) would give a step in $`P_0(k)`$¸. However, the slow-roll approximation is clearly not applicable. The exact solution for a local part of the spectrum near the point $`k_0`$ was obtained in . It reads:

$$k^3h^2(k)=\frac{18H_0^6}{V_{}^2}D^2(y),H_0=\sqrt{\frac{8\pi GV(\varphi _0)}{3}},V_\pm ^{}=V^{}(\varphi _0\pm 0)>0,y=\frac{k}{k_0},$$

$$D^2(y)=13\left(\frac{V_{}^{}}{V_+^{}}1\right)\frac{1}{y}\left(\left(1\frac{1}{y^2}\right)\mathrm{sin}2y+\frac{2}{y}\mathrm{cos}2y\right)$$
(28)   
$$+\frac{9}{2}\left(\frac{V_{}^{}}{V_+^{}}1\right)^2\frac{1}{y^2}\left(1+\frac{1}{y^2}\right)\left(1+\frac{1}{y^2}+\left(1\frac{1}{y^2}\right)\mathrm{cos}2y\frac{2}{y}\mathrm{sin}2y\right).$$

The function $`D^2(y)`$ has a step-like behaviour with superimposed oscillations. Since $`D(0)=V_{}^{}/V_+^{},D(\mathrm{})=1`$, the spectrum approaches the flat spectrum if $`|\mathrm{ln}(k/k_0)|1`$. As compared with the flat spectrum, the spectrum (28) has more power at large scales (small $`k`$) if $`V_{}^{}>V_+^{}`$, and more power on small scales in the opposite case. The shape of this function is universal (in the sence that it does not depend on a way of smoothing the jump in $`V^{}`$, if it is made in a sufficiently small vicinity of $`\varphi _0`$), it depends on the ratio $`V_{}^{}/V_+^{}`$ only.

In the absence of a cosmological term, we have to choose the case $`V_{}^{}>V_+^{}`$ to fit observational data, and then it is not possible to have a significant bump in the spectrum. However, if the cosmological term is positive, a possibility of an inverted step ($`V_{}^{}<V_+^{}`$) arises . In this case, the spectrum (28) can match the spectrum found in  rather well.

3. $`[V]0`$ at $`\varphi =\varphi _0`$.

In this case, there is no universal spectrum, and the answer depends on a concrete way of smoothing $`V(\varphi )`$ at $`\varphi =\varphi _0`$. Some general results are presented in , typically $`P_0(k)`$ acquires a large bump, however, a well may appear, too. Fortunately, there is no need in further consideration of this more complicated case, since observational data do not require such a strong non-analyticity. Looking at the results presented in  it is clear that they lie somewhere between the first two cases.

Therefore, a general lesson from these considerations is that pecularities in $`V(\varphi )`$ can produce local features in $`P_0(k)`$ where the slope $`n_S`$ is significantly different from unity. However, these features cannot be too sharp, in particular, both $`P_0(k)`$ and $`n_S(k)`$ are expected to be continuous functions of $`k`$.

So far, the treatment of peculiar points was purely mathematical. However, if we are seeking for a physical explanation of such behaviour of $`V(\varphi )`$, we have to go beyond the paradigm of one-field inflation to a more complicated case of two-field inflation. Inflation with two scalar field is a very rich physical model which includes double inflation, hybrid inflation, open inflation, etc. as specific cases. In our case it is sufficient to assume that the second scalar field $`\chi `$, in contrast to the inflaton field $`\varphi `$, is always in the regime $`|m_\chi ^2|H^2`$. So, it is not dynamically important during the whole inflation. However, it is coupled to $`\varphi `$ (e.g., through the term $`g^2\varphi ^2\chi ^2`$) and, as a result of change in $`\varphi `$ during inflation, the field $`\chi `$ experiences a fast phase transition approximately $`60`$ e-folds before the end of inflation. If parameters of an interaction potential $`V(\varphi ,\chi )`$ are such that the phase transition may be considered as an equilibrium one, then its net effect on inflation appears in the change of the equilibrium effective potential $`V_{eff}(\varphi )V(\varphi ,\chi _{eq}(\varphi ))`$ only. If the transition is a second-order one, with no jump in $`\chi _{eq}`$, the first case considered above takes place. If the transition is a first-order one, with a non-zero jump in $`\chi _{eq}`$, we arrive to the second case.

This research was partially supported by the Russian Foundation for Basic Research, Grant 96-02-17591, and by the Russian Research Project ”Cosmomicrophysics”.