Difference between revisions of "New problems"

From Universe in Problems
Jump to: navigation, search
(New on the matter creation model)
 
(2 intermediate revisions by the same user not shown)
Line 13: Line 13:
 
<div style="border: 1px solid #AAA; padding:5px;">
 
<div style="border: 1px solid #AAA; padding:5px;">
  
'' As we have seen above, negative pressure is the key ingredient required to achieve the accelerated expansion. Such pressure naturally appears when the system deviates from the thermodynamic equilibrium. In particular, as was first pointed out by Zel'dovich [Ya. B. Zel'dovich, JETP Lett. 12, 307 (1970)], the negative pressure is generated in process of particle creation due to the gravitational field. Physically, one may think that the (classical) time varying gravitational field works like a 'pump' supplying energy to the quantum fields. Construction of the models including the essentially quantum process of particle creation is problematic because of the difficulties connected with its inclusion into the classical field equations of Einstein. One can however avoid those difficulties on the phenomenological level [J.Lima, arXiv 0807.3379].
+
'' As we have seen above, negative pressure is the key ingredient required to achieve the accelerated expansion. Such pressure naturally appears when the system deviates from the thermodynamic equilibrium. In particular, as was first pointed out by Zel'dovich [Ya. B. Zel'dovich, JETP Lett. 12, 307 (1970)], the negative pressure is generated in process of particle creation due to the gravitational field. Physically, one may think that the (classical) time varying gravitational field works like a 'pump' supplying energy to the quantum fields. Construction of the models including the essentially quantum process of particle creation is problematic because of the difficulties connected with its inclusion into the classical field equations of Einstein. One can however avoid those difficulties on the phenomenological level [J.Lima, arXiv 0807.3379].''
We shall now consider an open thermodynamical system where the number of fluid particles $N$ is not preserved. So the particle conservation equation
+
 
 +
''We shall now consider an open thermodynamical system where the number of fluid particles $N$ is not preserved. So the particle conservation equations''
 
\[\dot N_{;\mu}=0\Rightarrow \dot n + \theta n=0\]
 
\[\dot N_{;\mu}=0\Rightarrow \dot n + \theta n=0\]
is now modified as
+
''is now modified as''
 
\[\dot n + \theta n=n\Gamma,\]
 
\[\dot n + \theta n=n\Gamma,\]
where $N_{\mu}=nu_\mu$ is the particle flow vector, $u_\mu$ is the particle four velocity, $\theta=u^\mu_{;\mu}$ is the fluid expansion, (for the FLRW Universe, $\theta=u^\mu_{;\mu}=3H$), $\dot n = n_{,\mu}u^{mu}$ and $\Gamma$ stands for the rate of change of the particle number in a comoving volume $V$, $\Gamma>0$ indicates particle creation while $\Gamma<0$ means particle annihilation. Any non-zero $\Gamma$ will effectively behave as a bulk viscous pressure of the thermodynamical fluid and nonequilibrium thermodynamics should come into the picture.''
+
''where $N_{\mu}=nu_\mu$ is the particle flow vector, $u_\mu$ is the particle four velocity, $\theta=u^\mu_{;\mu}$ is the fluid expansion, (for the FLRW Universe, $\theta=u^\mu_{;\mu}=3H$), $\dot n = n_{,\mu}u^{mu}$ and $\Gamma$ stands for the rate of change of the particle number in a comoving volume $V$, $\Gamma>0$ indicates particle creation while $\Gamma<0$ means particle annihilation. Any non-zero $\Gamma$ will effectively behave as a bulk viscous pressure of the thermodynamical fluid and nonequilibrium thermodynamics should come into the picture.''
  
  
Line 25: Line 26:
 
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
  
Show that the fact that in the present time $t_0$ the relation $H_0t_0\approx1$ holds unavoidably follows existance of the accelerated expansion stage of the Universe evolution.
+
Express the first law of thermodynamics in terms of the specific quantities: the energy density $\rho$ and the particle number density $n$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For a closed thermodynamical system the first law of thermodynamics states the conservation of internal energy $E$ as
 +
\[dE=dQ-pdV,\]
 +
where $p$ is the thermodynamic pressure, $V$ is any comoving volume and $dQ$ represents the heat received by the system in time $dt$. Now defining $\rho=E/V$ as the energy density, $n=N/V$ as the particle number density and $dq=dQ/N$ as the heat per unit particle, the above conservation law takes the form
 +
  \[d\left(\frac\rho n\right)=dq-pd\left(\frac1n\right).\]
 +
Note that this equation (known as Gibb's equation) is also true when the thermodynamical system is not closed, i.e., $N$ is not constant $[N=N(t)]$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 2'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
In the context of particle creation, let us consider spatially flat FRLW model of the Universe as an open thermodynamical system. Obtain the Friedmann equations for the cosmic fluid having energy-momentum tensor (in this section we will use the metric $(-,+,+,+)$
 +
\[T_{\mu\nu}=(\rho+p+\Pi)u_\mu u_\nu+(p+\Pi)g_{\mu\nu}.\]
 +
For cosmological models with particle creation the cosmic fluid may be considered as a perfect fluid and the dissipative term $\Pi$ as the effective bulk viscous pressure due to particle production. In other words, the cosmic substratum is not a conventional dissipative fluid, rather a perfect fluid with varying particle number.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using the standard procedure for transition from the Einstein's field equations to the Friedmann equations (see Chapter 2), one obtains
 +
\begin{align}
 +
\nonumber \rho & = 3H^2;\\
 +
\nonumber \dot H & = -\frac12(\rho+p+\Pi), \quad 8\pi G=1.
 +
\end{align}
 +
And the energy conservation relation $T^{\mu\nu}_{;\nu}$ gives
 +
\[\dot\rho+\theta(\rho+p+\Pi)=0.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 3'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Show that for the adiabatic process the effective bulk viscous pressure due to particle production is entirely determined by the particle production rate.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Represent the Gibb's relation
 +
\[dq=d\left(\frac\rho n\right)+pd\left(\frac1n\right)\]
 +
in the form
 +
\[Tds=d\left(\frac\rho n\right)+pd\left(\frac1n\right).\]
 +
Remember that here $s$ is the entropy per particle and $T$ is the temperature of the fluid. Using the conservation equations $\dot n + \theta n=n\Gamma$ and $\dot\rho+\theta(\rho+p+\Pi)=0$, we obtain
 +
\[nT\dot s=-\Pi\theta-\Gamma(\rho+p).\]
 +
For isentropic process, the entropy per particle remains constant, i.e. $\dot s=0$ and hence we have
 +
\[\Pi=-\frac\Gamma\theta(\rho+p).\]
 +
Thus, the effective bulk viscous pressure due to particle production is entirely determined by the particle production rate. So we may say that at least for the adiabatic process, a dissipative fluid is equivalent to a perfect fluid with varying particle number.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 4'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
An equivalent way to see the derivation of the field equations (7)-(8) (???) is to consider the energy-momentum tensor in the Einstein field equations as a total energy momentum tensor
 +
\[T^{(eff)}_{\mu\nu}=T^{(\gamma)}_{\mu\nu}+T^{(c)}_{\mu\nu},\]
 +
where $T^{(\gamma)}_{\mu\nu}$ is the energy-momentum tensor for the fluid with the EoS $p = (\gamma- 1) \rho$, i.e.
 +
  \[T^{(\gamma)}_{\mu\nu}=(\rho+p)u_\mu u_\nu+pg_{\mu\nu}.\]
 +
and \[T^{(c)}_{\mu\nu}=p_c\left(u_\mu u_\nu+g_{\mu\nu}\right)\] is the energy-momentum tensor which corresponds to the matter creation term. Show that energy-momentum tensor $T^{(eff)}_{\mu\nu}$ provides us with the matter creation pressure \[\Pi=-\frac\Gamma\theta(\rho+p),\] obtained in the previous problem.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For the considered case the Einstein field equations are
 +
\[G_{\mu\nu}=T^{(eff)}_{\mu\nu}=T^{(\gamma)}_{\mu\nu}+T^{(c)}_{\mu\nu}.\]
 +
The Bianchi identity gives
 +
\[g^{\nu\sigma}\left(T^{(\gamma)}_{\mu\nu}+T^{(c)}_{\mu\nu}\right)_{;\sigma}=0.\]
 +
Using the explicit form of tensors $T^{(\gamma)}_{\mu\nu}$ and $T^{(c)}_{\mu\nu}$, one finds that the latter equality is equivalent to the conservation equation
 +
\begin{equation}\label{2016_07_p04_eq1}
 +
\dot\rho+\theta(\rho+p+\Pi)=0
 +
\end{equation}
 +
From the other hand, Gibbs equation
 +
\[Tds=d\left(\frac\rho n\right)+pd\left(\frac1n\right).\]
 +
for the adiabatic process, taking into account the relation $\dot n +\theta n=n\Gamma$, directly leads to the conservation equation
 +
\begin{equation}\label{2016_07_p04_eq2}
 +
\dot\rho+\theta\left(1-\frac\Gamma\theta\right)(\rho+p)=0.
 +
\end{equation}
 +
Comparing the equations (\ref{2016_07_p04_eq1}) and (\ref{2016_07_p04_eq1}), one can reproduce the result, obtained in the previous problem \[\Pi=-\frac\Gamma\theta(\rho+p).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 5'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the effective parameter $w_{eff}$ of the EoS for the cosmological model with particle creation.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
\begin{align}
 +
\nonumber w_{eff} &=\frac{p_{tot}}{\rho_{tot}}=\frac{\rho+\Pi}\rho,\\
 +
\nonumber \Pi &=-\frac\Gamma\theta(p+\rho)=-\gamma\rho\frac\Gamma\theta,\quad w+1\equiv\gamma,\\
 +
\nonumber w_{eff} &=-1+\gamma\left(1-\frac\Gamma\theta\right).
 +
\end{align}
 +
We see that $w_{eff}$ represents quintessence  for $\Gamma<\theta$ ($\Gamma<3H$), phantom energy for $\Gamma>\theta$ ($\Gamma>3H$), and $\Gamma=\theta$ ($\Gamma=3H$) represents cosmological constant.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 6'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the deceleration parameter $q$ for the cosmological model with particle creation.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
The deceleration parameter $q$ is a measure of state of acceleration/deceleration of the Universe, defined as
 +
\begin{align}
 +
\nonumber q &=-1-\frac{\dot H}{H^2},\\
 +
\nonumber \dot H &=-\frac32\gamma H^2\left(1-\frac\Gamma\theta\right).
 +
\end{align}
 +
Consequently,
 +
\[q=-1+\frac32\gamma\left(1-\frac\Gamma\theta\right)=-1+\frac32\gamma\left(1-\frac\Gamma{3H}\right).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 7'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Obtain the equation for the scale factor in cosmological model with particle creation.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Combining the modified Friedmann equations with the EoS $p=w\rho$, one obtains the equation for the scale factor
 +
\[\frac{\ddot a}a +\frac{H^2}2\left[1+3w-\frac{1+w}H\Gamma\right]=0.\]
 +
Setting $\Gamma=0$, we reproduce the second Friedmann equation
 +
\[\frac{\ddot a}a=-\frac\rho6(1+3w).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 8'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Calculate the derivative $dS/dt$ and find $S(t)$ for the Creation Matter Model (CMM).
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
In an open thermodynamical system, the entropy change $dS$ can be decomposed into an entropy flow $d_fS$ and the entropy creation matter $d_cS$
 +
\[dS = d_fS+d_cS,\]
 +
with $d_cS>0$. As in a homogeneous system $d_fS=0$, but there is entropy production due to matter creation, so we have
 +
\begin{align}
 +
\nonumber \frac d{dt}\left(nsV\right)&=\dot nsV+ns\frac{dV}{dt},\quad \dot s=0,\\
 +
\nonumber ns\frac{dV}{dt}&=ns\frac{da^3}{dt}=ns3HV=nsV\theta,\\
 +
\nonumber \frac{dS}{dt} &=\frac d{dt}\left(nsV\right)=\dot nsV+nsV\theta=sV(\dot n +n\theta)=sVn\Gamma=S\Gamma.
 +
\end{align}
 +
Here $s$ is the entropy per particle. After integration we obtain
 +
\[S(t)=S_0\exp\left[3\int\limits_{a_0}^a\frac\Gamma\theta\frac{da}a\right]\]
 +
with $S_0=S(t_0)$ and $a_0=a(t_0)$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 9'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Calculate the scale factor dependence of the Universe's temperature for the Creation Matter Model (CMM).
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For a substance with zero chemical potential
 +
\[nsT=\rho+p.\]
 +
For isentropic process, the entropy per particle remains constant, i.e. $\dot s=0$ and consequently
 +
\[\dot nsT+ns\dot T=\dot\rho(1+w)=\dot\rho\gamma.\]
 +
Using
 +
\[\dot\rho=-\theta(\rho+p+\Pi),\quad \Pi=-\frac\Gamma\theta(\rho+p),\]
 +
we obtain
 +
\[\frac{\dot T}T=(1-\gamma)\theta\left(1-\frac\Gamma\theta\right).\]
 +
Using \[\frac{dT}{dt}=\frac{dT}{da}Ha\] and $\theta=3H$, we find
 +
\[T=T_0a^{-3(\gamma-1)}\exp\left[3(\gamma-1)\int\limits_{a_0}^a\frac\Gamma\theta\frac{da}a\right],\]
 +
or using result of the previous problem
 +
\[T=T_1\left(\frac S{a^3}\right)^{\gamma-1},\quad T_1\equiv\frac{T_0}{S_0^{\gamma-1}}.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 10'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Solve the previous problem using the thermodynamic  identity \[T\left(\frac{\partial p}{\partial T}\right)_n+n\left(\frac{\partial\rho}{\partial n}_T\right)=\rho+p.\]
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Combining the thermodynamic identity with the conservation equation, we obtain
 +
\[\frac{\dot T}T=\left(\frac{\partial p}{\partial\rho}\right)_n\frac{\dot n}n.\]
 +
In particular, for $\Gamma=0$, $\dot n +3nH=0$ and radiation or ultra relativistic particles ($p=\rho/3$) the temperature evolution law becomes
 +
\[\frac{\dot T}T=-\frac{\dot a}a\Rightarrow aT=const.\]
 +
If gravitational particles source is present, then $\dot n +3nH=n\Gamma$ and
 +
\[\frac{\dot T}T=\left(\frac{\partial p}{\partial\rho}\right)_n\frac{\dot n}n=\left(\frac{\partial p}{\partial\rho}\right)_n\left(\Gamma-3H\right).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 11'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Solve the previous problem for cosmological substance with $\partial p/\partial\rho=w=const$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using
 +
\[\dot H+\frac32(1+w)H^2\left(1-\frac\Gamma{3H}\right)=0,\]
 +
we find
 +
\[\Gamma-3H=\frac2{1+w}\frac{\dot H}H\]
 +
and
 +
\[\frac{\dot T}T=\frac{2w}{1+w}\frac{\dot H}H.\]
 +
Consequently, the temperature of the fluid is governed by the following law
 +
\[T=T_0\left(\frac H{H_0}\right)^{\frac{2w}{1+w}}.\]
 +
 
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 12'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Show that the matter creation process can be effectively described by the dimensionless ratio \[\beta=\frac\Gamma\theta=\frac\Gamma{3H}.\]
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
If $\Gamma\ll3H$, then $\beta\ll1$, and the matter creation process is qualitatively small. It leads to $n\propto R^{-3}$ and $H=2/(3t)$, as it should be for the Einstein-de Sitter model. The opposite case $\Gamma\gg3H$ corresponds to the extreme physical situation when the matter creation process is so intense that it compensates with excess the matter density decrease due to the expansion. behavior of such type could probably be realized only in very early Universe, as it took place, for example, during the inflation stage, which corresponded to the reheating. The intermediate and physically more realistic situation appears if this ratio is less or of order of unity: $\Gamma\le3H$. In particular, if $\Gamma=3H$, then density decrease due to the expansion is exactly compensated by its increase due to the matter creation, so the particle number density remains constant.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 13'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
How is the evolution of $a(t)$ affected by rate of change of the particle number $\Gamma$?
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
In the considered model, the Hubble parameter for the Universe, filled with the cosmic fluid with EoS $p=w\rho$, satisfies the dynamical equation
 +
\[\dot H+\frac32(1+w)H^2\left(1-\frac\Gamma{3H}\right)=0.\]
 +
The de Sitter solution $\dot H=0$ for $\Gamma=3H=constant$ is  possible regardless of the EoS defining the cosmic fluid. Since the Universe is evolving, such a solution is unstable, and, as long as $\Gamma\ll3H$, conventional solutions without particle production are recovered. From the above equation, one may conclude that the main effect of $\Gamma$ is to provoke a dynamic instability in the space-time thereby allowing a transition from a de Sitter regime ($\Gamma\approx3H$) to a conventional solution, and vice versa.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 14'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Let us assume a radiation dominated Universe ($w=1/3$, $\Gamma=\Gamma_r$). The dynamics is determined by the ratio $\Gamma_r/(3H)$ (see previous problem). The particle production must be strongly suppressed, $\Gamma_r/(3H)\ll1$, when the Universe enters the radiation phase. The simplest formula satisfying such a criterion is linear, namely: $\Gamma_r/(3H)=H/H_I$, where $H_I$ is the inflationary expansion rate associated with the initial de Sitter expansion ($H\le H_I$). Find the Hubble parameter dependence on the scale factor $H(a)$ for the radiation phase.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For $w=1/3$, $\Gamma=\Gamma_r$, $\Gamma_r/(3H)=H/H_I$,
 +
\[\dot H+\frac32(1+w)H^2\left(1-\frac\Gamma{3H}\right)=0\Rightarrow\dot H+2H^2\left(1-\frac H{H_I}\right)=0,\]
 +
or
 +
\[a\frac{dH}{da}+2H^2\left(1-\frac H{H_I}\right)=0\Rightarrow\frac{dH}{H-H_I}-\frac{dH}H=2\frac{da}a.\]
 +
The solution of this equation is
 +
\[H(a)=\frac{H_I}{1+Da^2}.\]
 +
Here $D\ge0$ is an integration constant.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 15'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Using the result obtained in the previous problem, show that in the considering model the Universe evolves continuously from de Sitter stage to the standard radiation phase.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Integrating the equation
 +
\[H(a)=\frac{H_I}{1+Da^2},\]
 +
we obtain
 +
\[H_It=\ln\left(\frac a{a_*}\right)+\frac12Da_*^2\left(\frac a{a_*}\right)^2,\]
 +
where the integration constant$a_*$ defines the transition from the de Sitter stage to the beginning of the standard radiation epoch. At early times ($a\ll a_*$), when the logarithmic term dominates, one finds $a\approx a_*\exp(H_It)$ (de Sitter stage), while at late times, with $a\gg a_*$ and $H\ll H_I$ the obtained result reduces to \[a\approx\left(\frac2D\right)^{1/2}(H_It)^{1/2}\] and the standard radiation phase is reached.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 16'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the radiation energy density $\rho_r$ as function of scale factor for the rate of change of the particle number $\Gamma=\Gamma_r$, $\Gamma_r/(3H)=H/H_I$ ($H_I$ is the inflationary expansion rate).
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
The radiation energy density is $\rho_r(a)=3H^2(a)$ ($8\pi G=1$) and using \[H(a)=\frac{H_I}{1+Da^2},\] we obtain
 +
\[\rho_r(a)=\frac{3H_I^2}{(1+Da^2)^2}\]
 +
or
 +
\[\rho_r(a)=\frac{\rho_I}{\left[1+Da_*^2\left(\frac a{a_*}\right)^2\right]^2}\]
 +
Here $\rho_I=3H_I^2$ is the greatest value of the energy density, characterizing the initial de Sitter stage and the integration constant $a_*$ defines the transition from the de Sitter stage to the beginning of the standard radiation epoch. We see again that the conventional radiation phase, $\rho_r\propto a^{-4}$, is attained when $a\gg a_*$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 17'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
For cosmological model considered in the problems \ref{2016_07_p14}-\ref{2016_07_p16} find how did the cosmic temperature evolve at the very early stages?
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For adiabatic production of relativistic particles the energy density scales as $\rho_r\propto T^4$ and the equation obtained in previous problem
 +
\[\rho_r(a)=\frac{\rho_I}{\left[1+Da_*^2\left(\frac a{a_*}\right)^2\right]^2}\]
 +
implies that
 +
\[T=\frac{T_I}{\sqrt{1+Da_*^2\left(\frac a{a_*}\right)^2}}\]
 +
where $T_I$ is the temperature of the initial de Sitter phase which must be uniquely determined by the scale $H_I$ - the inflationary expansion rate associated to the initial de Sitter phase. We see that the expansion proceeds isothermally: $T=const$ during the de Sitter phase ($a\ll a_*$) .
 +
 
 +
After the de Sitter stage, the temperature decreases continuously in the course of the expansion. For $a\gg a_*$ and $H\ll H_I$, we obtain $T\propto a^{-1}$. Accordingly, the comoving number of photons becomes constant since $n\propto a^{-3}$, as expected for the standard radiation stage.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 18'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Solve Problem \ref{2016_07_p14} for a substance with EoS $p=w\rho=(\gamma-1)\rho$ and $\Gamma=\Gamma_0H^2$ with $\Gamma_0=3\beta/H_e$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For this choice of $\Gamma$, $H$ can be found from
 +
\[\dot H+\frac32(1+w)H^2\left(1-\frac\Gamma{3H}\right)=0,\]
 +
as
 +
\[H(a)=\frac{H_e}{\beta-(1-\beta)\left(\frac a{a_e}\right)^{\frac{3\gamma}2}},\]
 +
where $H_e$ and $a_e$ are chosen to be the values of the Hubble parameter and the scale factor respectively at some instant. If we identify $a_e$ with some intermediate value of the scale factor $a$ corresponding to $\ddot a=0$, i.e. to the transition from de Sitter accelerating phase to the standard decelerating matter phase, then we have \[\dot H+H^2=0\Rightarrow\dot H_e=-H^2_e,\] and using
 +
\[\frac\Gamma{3H}=1+\frac2{3\gamma}\frac{\dot H}{H^2},\]
 +
we find \[\beta=1-\frac2{3\gamma}.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 19'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the deceleration parameter for the model considered in previous problem.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using
 +
\[q=-1+\frac32\gamma\left(1-\frac\Gamma{3H}\right)\]
 +
for
 +
\[H(a)=\frac{H_e}{\beta-(1-\beta)\left(\frac a{a_e}\right)^{\frac{3\gamma}2}}\]
 +
we obtain
 +
\[q(z)=-1+\frac{3\gamma}2\left[1-\frac{\beta}{\beta+(1-\beta)(1+z)^{-3\gamma/2}}\right],\quad \frac{a_e}a=1+z.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 20'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Let us now consider how the CMM (Creation matter model) works at later stages of the Universe's evolution in the transition from Einstein-de Sitter to a late time de Sitter stage. Due to the conservation of baryon number one have to take into account only the production rate of cold dark matter particles $\Gamma_{dm}$. What is the exact form of $\Gamma_{dm}$?
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
The evolution equation for the Hubble parameter in the flat $\Lambda$CDM reads
 +
\begin{equation}\label{2016_07_p20_e1}
 +
\dot H+\frac32H^2\left[1-\left(\frac{H_f}H\right)^2\right]=0,
 +
\end{equation}
 +
where $H_f$ is Hubble parameter of the final de Sitter stage $H\ge H_f$. All available observations are in accordance with the $\Lambda$CDM evolution both at the background and perturbative levels. Thus the form of $\Gamma_{dm}$ can be found from the requirement that the evolution equation for the model with the production of cold dark matter particles ($w=0$, $\Gamma=\Gamma_{dm}$)
 +
\begin{equation}\label{2016_07_p20_e2}
 +
\dot H+\frac32H^2\left(1-\frac{\Gamma_{dm}}{3H}\right)=0,
 +
\end{equation}
 +
leads to the same background evolution. Comparing equations (\ref{2016_07_p20_e1}) and (\ref{2016_07_p20_e2}) we find that for the fulfillment of this condition it is necessary that
 +
\[\frac{\Gamma_{dm}}{3H}=\left(\frac{H_f}H\right)^2,\]
 +
or $\Gamma_{dm}\propto H^{-1}$. The limiting value of the creation rate, $\Gamma_{dm}=3H_f$, leads to a late time de Sitter phase: $\dot H=0$, $H=H_f$ thereby showing that the de Sitter solution now becomes an attractor at late times.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 21'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Let us consider the dynamical system in the context of bulk viscosity. In that case, we need to replace the pressure as $p\to p-3\zeta H$, where $\zeta$ is the coefficient of bulk viscosity. At what choice of the particle production rate $\Gamma$ the creation matter model would be equivalent to the bulk viscosity model?
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
In presence of the bulk viscosity the Friedmann equations are modified as
 +
\begin{align}
 +
\nonumber \rho &=3H^2,
 +
\nonumber \dot H &=-\frac12(\rho+p)+\frac32\zeta H.
 +
\end{align}
 +
Using the barotropic EoS: $p=(\gamma-1)\rho$, the second Friedmann equation can be written as
 +
\[\dot H=-\frac32\gamma H^2\left(1-\frac\zeta{\gamma H}\right),\]
 +
which coincides with
 +
\[\dot H=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right),\]
 +
if one takes
 +
\[\zeta = \frac\Gamma3\gamma.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
== Dynamical analysis of the Creation Matter Model ==
 +
 
 +
'''Problem 22'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Since the equation
 +
\[\dot H=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right)\]
 +
is a one-dimensional first order differential equation, hence, the dynamics is obtained from the study of its critical points (or, fixed points) in which $\dot H=0$. Show that the fixed points $H_*$ in which $H_*\ne0$ represent the de Sitter Universe.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
If $H=H_*$ is the fixed point, then $H_*=0$ or $\Gamma(H_*)=3H_*$. The fixed points $H_*\ne0$ correspond to $H=const$ and these solutions describe the de Sitter Universe.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 23'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the fixed points and determine their nature for the case with the particle creation rate $\Gamma=n/H$, $n>0$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Let us consider the general form of the equation $\dot h=f(H)$. If \[\frac{df(H_*)}{dH}<0\] at the fixed point $H_*$, then the fixed point is asymptotically stable (attractor), and on the other hand, if we have \[\frac{df(H_*)}{dH}<0,\] then the fixed point has unstable nature (repeller). The repeller point is suitable for early Universe, since it can describe the inflationary epoch, whereas the attractor point can describe the stable accelerating phase for the later time.
 +
 
 +
To find the fixed points for the considered case with the particle creation rate $\Gamma=n/H$, $n>0$, one should solve the equation
 +
\[-\frac32\gamma H^2\left(1-\frac n{3H^2}\right)=0.\]
 +
The non-zero solutions are $H_*=\pm\sqrt{n/3}$. For the above given form for $\Gamma$, one has
 +
\[f(H)=-\frac32\gamma\left(H^2-\frac n3\right)\]
 +
and thus
 +
\[\frac{df(\pm\sqrt{n/3})}{dH}=\mp\sqrt{3n},\]
 +
which means that $H_*=\sqrt{n/3}$ is an attractor and $H_*=-\sqrt{n/3}$ is a repeller.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 24'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the fixed points for the case with the particle creation rate $\Gamma=-\Gamma_0+mH+n/H$, $n>0$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
If $\Gamma(H)$ is a polynomial function of $H$, then the fixed point condition $\Gamma(H_*)=3H_*$ (for $H_*\ne0$), is a polynomial equation which has as many solutions (not necessary real solutions) as the highest power of the polynomial equation $\Gamma(H_*)=3H_*$. Hence, with the given particle creation rate we have to consider the following second-order polynomial equation
 +
\[(m-3)H_*^2-\Gamma_0H_*+n=0.\]
 +
In order to have the two critical points, as many as the inflationary phases of the Universe, we are interested in the case when $m\ne3$, and \[n\ne\frac{\Gamma_0^2}{4(m-3)}.\] Solving this equation, we find the critical points
 +
\[H_\pm=\frac{\Gamma_0}{2(m-3)}\left(1\pm\sqrt{1+\frac{4n(3-m)}{\Gamma_0^2}}\right).\]
 +
The case $m=3$ is special in a sense that there is only one critical point $H_-=n/\Gamma_0$, which is always an attractor for $\Gamma_0>0$, and a repeller for $\Gamma_0<0$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 25'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Perform dynamic analysis of the solutions, obtained in previous problem.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
To perform the dynamical analysis, we have to divide the parameter space ($m,n\Gamma_0$) to different regions by the following conditions:
 +
\[\left(\Gamma_0>0,\ \Gamma_0<0;\quad m>3,\ m<3;\quad n>0,\ n<0;\quad\frac{4(3-m)}{\Gamma_0^2}>-1,\ \frac{4(3-m)}{\Gamma_0^2}<-1.\right)\]
 +
To have a non-singular Universe (without the Big Bang singularity) with the accelerated phases both at early and late times, one possibility is to have two critical points $H_+>H_->0$, where $H_+$ is a repeller and $H_-$ must be an attractor. If so, in principle, when the Universe leaves the point $H_+$, realizing the inflationary phase, and when it comes asymptotically to $H_-$, it enters into the current accelerated phase. Of course, the viability of the background has to be checked dealing with cosmological perturbations and comparing the theoretical predictions with the observations.
 +
 
 +
For the considered model, the above described scenario can only happen in the region of the parameter space  given by
 +
\[(m,n,\Gamma_0)=\left\{\Gamma_0>0,\ m>3,\ n\ge0,\ \frac{4(3-m)}{\Gamma_0^2}>-1\right\}.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 26'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Prove that the dynamics generated by the creation matter model is equivalent to that driven by the single scalar field.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
We use the energy density $\rho_\varphi$ and pressure $p_\varphi$ of the scalar field $\varphi$ given by
 +
\begin{align}
 +
\nonumber \rho_\varphi &=\frac12\dot\varphi^2+V(\varphi),\\
 +
\nonumber p_\varphi &=\frac12\dot\varphi^2-V(\varphi).
 +
\end{align}
 +
Recall that the creation matter model is described by Friedmann equations
 +
\begin{align}
 +
\nonumber 3H^2&=\rho,\\
 +
\nonumber \dot H &=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right).
 +
\end{align}
 +
To show the equivalence, we perform the replacement
 +
\[\rho\to\rho_\varphi;\quad p-\frac\Gamma{3H}\gamma\rho\to p_\varphi.\]
 +
The Friedmann equations will become
 +
\begin{align}
 +
\label{2016_07_p26_e1} 3H^2&=\rho_\varphi,\\
 +
\nonumber 2\dot H &=-\dot\varphi^2.
 +
\end{align}
 +
Note that Eqs. (\ref{2016_07_p26_e1}) uses the equations of General Relativity  for a single scalar field. It means that we are dealing with the equivalence between an open system (the creation matter model) and the one driven by a single scalar field in the context of General Relativity.
 +
 
 +
Further, the effective EoS parameter is given by
 +
\[w_{eff}=-1-\frac23\frac{\dot H}{H^2}=-1+\gamma\left(1-\frac\Gamma{3H}\right).\]
 +
Of course,
 +
\[w_{eff}=w_\varphi=\frac{\dot\varphi^2-2V}{\dot\varphi^2+2V}.\]
 +
Note that $2\dot H=-\dot\varphi^2<0$ implies that $w_{eff}>-1$, and thus one has $\Gamma<3H$ (the quintessence era).
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 27'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the scalar field $\varphi(H)$ and the potential $V(\varphi)$, satisfying the condition that the dynamics generated by the creation matter model is equivalent to the dynamics driven by a single scalar field.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
From the Friedmann equations (see the previous problem)
 +
\begin{align}
 +
\nonumber 3H^2&=\rho_\varphi,\\
 +
\nonumber 2\dot H &=-\dot\varphi^2;
 +
\end{align}
 +
with
 +
\begin{align}
 +
\nonumber \rho_\varphi&=\frac12\dot\varphi^2+V(\varphi),\\
 +
\nonumber \dot H &=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right);
 +
\end{align}
 +
we obtain
 +
\[V(\varphi)=\frac32H^2\left(2-\gamma+\frac{\gamma\Gamma}{3H}\right),\]
 +
\[\dot\varphi=\sqrt{-2\dot H}=\sqrt{3\gamma H^2\left(1-\frac\Gamma{3H}\right)}.\]
 +
Performing the change of variables \[dt=\frac{dH}{\dot H},\] we find
 +
\[\varphi=\int\sqrt{-\frac2{\dot H}}dH=-\frac2{\sqrt{\gamma}}\int\frac{dH}{\sqrt{3H^2-\Gamma H}}.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 28'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Solve the previous problem for case \[\Gamma=-\Gamma_0+mH+n/H.\]
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
In the particular case $\Gamma=-\Gamma_0+mH+n/H$ one has
 +
\[\varphi=-\frac2{\sqrt{\gamma}}\int\frac{dH}{\sqrt{(3-m)H^2+\Gamma_0H-n}}.\]
 +
In the domain of interest (see Problem \ref{2016_07_p25})
 +
\[(m,n,\Gamma_0)=\left\{\Gamma_0>0,\ m>3,\ n\ge0,\ \frac{4(3-m)}{\Gamma_0^2}>-1\right\},\]
 +
the integral could be solved analytically, giving as a result
 +
\[\varphi=\frac2{\sqrt{(m-3)\gamma}}\arcsin\left[\frac{m-3}\omega\left(\frac{\Gamma_0}{m-3}-2H\right)\right],\quad \omega\equiv\sqrt{\Gamma_0^2+4(3-m)n}.\]
 +
Inverting this expression, we find
 +
\[H=\frac1{2(m-3)}\left[\Gamma_0-\omega\sin\left(\frac{\sqrt{(m-3)\gamma}}2\varphi\right)\right].\]
 +
Using
 +
\[V(\varphi)=\frac32H^2\left(2-\gamma+\frac{\gamma\Gamma}{3H}\right)=\frac12\left[\left(6+(m-3)\gamma\right)H^2 -\gamma\Gamma_0H+\gamma n\right],\]
 +
we obtain
 +
\[V(\varphi)=\frac3{4(m-3)^2}\left[\Gamma_0-\omega\sin\left(\frac{\sqrt{(m-3)\gamma}}2\varphi\right)\right]^2 - \frac{\gamma\omega^2}{8(m-3)}\cos^2\left(\frac{\sqrt{(m-3)\gamma}}2\varphi\right).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 29'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Prove equivalence of the dynamics generated by the creation matter model with the dynamics driven by a decaying vacuum energy density model.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For a generic decaying vacuum model Friedmann equations reduce to
 +
\begin{align}
 +
\label{2016_07_p29_e1} \rho+\Lambda(t)&=3\left(H^2+\frac k{a^2}\right),\\
 +
\nonumber p-\Lambda(t)&=-2\frac{\ddot a}a-\left(H^2+\frac k{a^2}\right),
 +
\end{align}
 +
where $\rho$ and $p$ are the energy density and the pressure of the dominant fluid component (dark matter) with EoS $p=w\rho$ ($p=(\gamma-1)\rho$).
 +
The decaying vacuum causes a change in the number of particles of dark matter, so the equation describing particle concentration has a source term, i.e.,
 +
\[\dot n+3Hn=n\Gamma^*.\]
 +
Here, $\Gamma^*$ is the rate of change of the number of particles.
 +
By combining Friedmann equations (\ref{2016_07_p29_e1}), or more directly, from the total energy conservation law one finds
 +
\begin{equation}\label{2016_07_p29_e2}
 +
\dot\rho+3H(\rho+p)=-\dot\Lambda(t).
 +
\end{equation}
 +
Since in the considered model the vacuum decay is the unique source of particle creation, we can write
 +
\[\dot\Lambda(t)=-\zeta n\Gamma,\]
 +
where $\zeta$ is a positive phenomenological parameter.
 +
At last, we can combine Friedmann equations in order to obtain the dynamics of decaying vacuum model,
 +
\begin{equation}\label{2016_07_p29_e3}
 +
\frac{\ddot a}a + \Delta(H^2+k^2)-\frac{1+w}2\Lambda(t)=0,\quad \Delta\equiv\frac{3w+1}2.
 +
\end{equation}
 +
A similar equation for the creation matter model has the form
 +
\begin{equation}\label{2016_07_p29_e4}
 +
\frac{\ddot a}a + \Delta(H^2+k^2)+\frac12\Pi=0,
 +
\end{equation}
 +
where the term $\Pi$ is the effective bulk  pressure due to particle production (see the previous Problems).
 +
Comparing (\ref{2016_07_p29_e3}) and (\ref{2016_07_p29_e4}) we find
 +
\[\Pi=-(1+w)\Lambda(t)=-\gamma\Lambda(t).\]
 +
We have seen above that \[\Pi=-\gamma\rho\frac\Gamma{3H}.\] Consequently, both the theories are equivalent, if we set \[\Gamma=3H\frac\Lambda\rho.\]
 +
Note that above identification holds regardless of the curvature of the Universe.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 30'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Show that for spatially flat Universe the inequalities $\Gamma\ll H$ and $\Lambda\ll H^2$ should be satisfied simultaneously.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For a spatially flat geometry, we have $\rho=3H^2$ and consequently
 +
\[\Gamma=3H\frac\Lambda\rho\Rightarrow\frac\Gamma H=\frac\Lambda{H^2}.\]
 +
In particular, if $\Gamma\ll H$, we find that $\Lambda\ll H^2$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 31'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Let us consider the particular case $\Lambda(t)=const$. Show that for such choice the creation matter model corresponds to the standard LCDM.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
The conservation equation for the creation matter model is
 +
\[\dot\rho+3H(\rho+p+\Pi)=0.\]
 +
For cold dark matter with $w=0$ ($\gamma=1$) and a spatially flat geometry $\Pi=-\Gamma H$ and the conservation equation transforms into
 +
\begin{equation}\label{2016_07_p31_e1}
 +
\dot\rho+3H(\rho-\Gamma H)=0.
 +
\end{equation}
 +
In the considered case \[\Gamma=3H\frac\Lambda\rho=\frac\Lambda H\] and (\ref{2016_07_p31_e1}) takes on the form
 +
\begin{equation}\label{2016_07_p31_e2}
 +
\dot\rho+3H(\rho-\Lambda)=0.
 +
\end{equation}
 +
Performing the integration, one finds
 +
\begin{equation}\label{2016_07_p31_e3}
 +
\rho(a)=\Lambda+\rho_{dm0}a^{-3},
 +
\end{equation}
 +
where $\rho_{dm0}$ is an integration constant that must quantify the current amount of matter that is clustering to the present time. The equation (\ref{2016_07_p31_e3}) describes the dynamics of a creation matter model that behaves like the LCDM model.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 32'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Calculate the cosmographic parameter $j(H)$ for the cosmological model with particle creation.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
\[j=\frac{\ddot H}{H^3}-3q-2;\]
 +
\[\dot H=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right);\]
 +
\[q=-1+\frac32\gamma\left(1-\frac\Gamma{3H}\right).\]
 +
Let us introduce the function
 +
\[F\equiv\frac32\gamma\left(1-\frac\Gamma{3H}\right).\]
 +
In terms of this function
 +
\[\dot H=-H^2 F.\]
 +
As \[\frac d{dt}=\dot H\frac d{dH},\] we obtain
 +
\[\ddot H=\frac d{dt}\dot H=\dot H\frac{d\dot H}{dH}=-\dot H\frac d{dH}(H^2F)=H^3\left(2F+H\frac{dF}{dH}\right),\]
 +
\[\frac{dF}{dH}=-\frac12\frac\gamma H\left(\frac{d\Gamma}{dH}-\frac\Gamma H\right).\]
 +
Consequently,
 +
\[j=1+F\left(2F+H\frac{dF}{dH}\right)-3F,\]
 +
\[\frac{dF}{dH}=-\frac12\frac\gamma H\left(\frac{d\Gamma}{dH}-\frac\Gamma H\right),\]
 +
\[j=1+F\left[2F-3-\frac\gamma2\left(\frac{d\Gamma}{dH}-\frac\Gamma H\right)\right].\]
 +
The result can be checked by passing to the LCDM limit ($\Gamma\to0$, $\gamma=1$ or $\gamma=-1$ ???). In this case $F\to3/2$ or $F\to0$. In both cases $j=1$ in accordance with the LCDM.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 33'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Solve the previous problem using the cosmographical relation \[j=-\frac1H\frac{dq}{dt}+q(1+2q).\]
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using $q=-1+F$ and \[\dot H=-H^2 F,\quad F\equiv\frac32\gamma\left(1-\frac\Gamma{3H}\right),\] we find
 +
\[\frac{dq}{dt}=\dot H\frac{dF}{dH}=-H^2 F\frac{dF}{dH}\]
 +
and
 +
\[j=1+F\left(2F+H\frac{dF}{dH}\right)-3F,\]
 +
\[\frac{dF}{dH}=-\frac12\frac\gamma H\left(\frac{d\Gamma}{dH}-\frac\Gamma H\right),\]
 +
\[j=1+F\left[2F-3-\frac\gamma2\left(\frac{d\Gamma}{dH}-\frac\Gamma H\right)\right],\]
 +
which coincides with the result of the previous problem.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'' However, in connection with the particle production at the expense of gravitational field of the expanding Universe, we recall that long back ago, Zeldovich [Ya. B. Zeldovich, JETP Lett. 12, 307 (1970)] introduced some bulk viscosity mechanism which is responsible for particle production. However, later on, Lima and Germano [J. A. S. Lima and A. S. M. Germano, Phys. Lett. A 170, 373 (1992); For a more general macroscopic approach see also R. Silva, J. A. S. Lima and M. O. Calv~ao, Gen. Rel. Grav. 34, 865 (2002).] showed that although both the processes, namely the bulk viscosity mechanism by Zeldovich [Ya. B. Zeldovich, JETP Lett. 12, 307 (1970).] and the gravitational particle production produce the same dynamics of the Universe, but in principle they are completely different from a thermodynamical point of view. Since we describe the particle production, we would like also to note an analogy which exists between the models driven by the particle production and the models of Steady State Cosmology developed in [F. Hoyle and J. V. Narlikar, Proc. Roy. Soc. A 282, 191 (1964); ibdem, 290, 143 (1966); see also Action at a Distance in Physics and Cosmology, 1974 (New York, Freeman).]. ''
 +
 
 +
'''Problem 34'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
In the model with $\Gamma=-\Gamma_0+mH$, where $0\le m\le3$ and $-3H_0\le\gamma\le0$, find the time dependence of the scale factor.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
For $\Gamma=-\Gamma_0+mH$ the evolution equation for the hubble parameter takes on the form
 +
\[\dot H+\frac32H^2\left(1-\frac m3+\frac{\Gamma_0}{3H}\right).\]
 +
Its solution reads
 +
\begin{equation}\label{2016_07_p34_e1}
 +
H(t)=\frac{\Gamma_0/(m-3)}{1-\exp(\Gamma_0t/2)}.
 +
\end{equation}
 +
Then for the scale factor one finds
 +
\begin{equation}\label{2016_07_p34_e2}
 +
a(t)=a_0\left[H_0\left(\frac{m-3}{\Gamma_0}-1\right)\left(e^{-\Gamma_0t/2}-1\right)\right]^{\frac2{3-m}}.
 +
\end{equation}
 +
In the limit $\Gamma_0\to0$
 +
\[a(t)=a_0\left[\frac{3-m}2H_0t\right]^{\frac2{3-m}}.\]
 +
Naturally at $m=0$ this solution is reduced to that of the Einstein-de Sitter model (Universe filled with non-relativistic matter).
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 35'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
In the model considered in the previous problem, express age of the Universe in terms of the parameters $\Gamma_0$ and $m$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Setting $a=a_0$ in the expression for the scale factor obtained in the previous problem, one can find age of the Universe:
 +
\begin{equation}\label{2016_07_p35_e1}
 +
t_0=2\Gamma_0\ln\left(1-\frac{\Gamma_0/H_0}{m-3}\right).
 +
\end{equation}
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 36'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find deceleration parameter in the model with matter creation rate equal to $\Gamma$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
As we have seen above, in the considered model the evolution equation for the scale factor takes on the form
 +
\[a\ddot a+\frac12\left(1-\frac\Gamma H\right)\dot a^2=0.\]
 +
Using the definition of the deceleration parameter \[q\equiv-\frac{\ddot aa}{\dot a^2},\] one finds
 +
\[q=\frac12\left(1-\frac\Gamma H\right).\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 37'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Show that in the creation matter model with $\Gamma=-\Gamma_0+mH$ in the case $\Gamma_0=0$ the expansion of Universe is always accelerated if $m>1$ and always decelerated if $m<1$.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using the result of the previous problem, one obtains
 +
\[q=\frac12\left[1-m+\Gamma_0/H\right].\]
 +
It immediately follows then that at $\Gamma_0=0$ the deceleration parameter does not depend on time. It easy to see that at $m>1$ the expansion of Universe is always accelerated, and it is always decelerated for $m<1$. As the transition from the decelerated expansion to the accelerated one is required by the observation, the expression for the matter creation rate $\Gamma$ must contain the constant term $-\Gamma_0\ne0$.
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 38'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the deceleration parameter as a function of the redshift in the creation matter model.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Exclude time from the above obtained expressions for the scale factor (\ref{2016_07_p34_e2}) and Hubble parameter (\ref{2016_07_p34_e1}), and use the relation \[1+z=\frac{a_0}a\] to find the hubble parameter as function of the redshift \[H(z)=\frac{\Gamma_0}{m-3}+\left(H_0-\frac{\Gamma_0}{m-3}\right)(1+z)^{\frac{3-m}2}.\] Then
 +
\begin{equation}\label{2016_07_p38_e1}
 +
q(z)=\frac12\left[1-m+\frac{\Gamma_0}{H(z)}\right] =\frac12\left[1-m+\frac1{\frac1{m-3}+\left(\frac{H_0}{\Gamma_0}-\frac1{m-3}\right)(1+z)^{\frac{3-m}2}}\right].
 +
\end{equation}
 +
In the limit $\Gamma_0\to0$ \[q=\frac{1-m}2\] and for $m\to0$ \[q(z)=\frac12\left[1+\frac{\Gamma_0/H_0}{(1+\frac{\Gamma_0}{3H_0})(1+z)^{\frac32}-\frac{\Gamma_0}{3H_0}}\right].\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 39'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Find the redshift value $z_t$ corresponding to the transition from the decelerated expansion to the accelerated one in the creation matter model.
 +
 
 +
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 +
<p style="text-align: left;">
 +
 
 +
Using the expression (\ref{2016_07_p38_e1}) for $q(z)$ obtained in the previous Problem with the condition $q(z_t)=0$, one finds that the transition from the decelerated expansion to the accelerated one in the cobsidered model takes place at the redshift value
 +
\[z_t=\left[\frac{2\Gamma_0/H_0}{(1-m)(m-3-\Gamma_0/H_0)}\right]^{\frac2{3-m}}-1.\]
 +
 
 +
</p>  </div></div></div>
 +
 
 +
<div id="150_dp1"></div>
 +
 
 +
<div style="border: 1px solid #AAA; padding:5px;">
 +
 
 +
'''Problem 40'''
 +
 
 +
<p style= "color: #999;font-size: 11px">problem id: 150_dp0</p>
 +
 
 +
Express age of the Universe in terms of the $z_t$ (see the previous Problem) in the creation matter model with $m=0$.
  
 
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 
<div class="NavFrame collapsed">  <div class="NavHead">solution</div>  <div style="width:100%;" class="NavContent">
 
<p style="text-align: left;">
 
<p style="text-align: left;">
  
The combination of the definitions
+
For the case $m=0$
\[H=\frac{\dot a}{a}=\frac{d\ln a}{dt},\quad q(t)=-\frac{\ddot a a}{\dot a^2}\]
+
\[\frac{\Gamma_0}{3H_0}=-\frac{1}{1+2(1+z_t)^{-\frac32}}.\]
give the equation
+
Substituting the latter expression into the result (\ref{2016_07_p35_e1}) of the Problem \ref{2016_07_p35} \[t_0=2\Gamma_0\ln\left(1-\frac{\Gamma_0/H_0}{m-3}\right),\] one obtains
\[\dot H+(1+q)H^2=0.\]
+
\[H_0t_0=\frac{4+2(1+z_t)^{\frac32}}{3(1+z_t)^{\frac32}}\ln\left[1+\frac{(1+z_t)^{\frac32}}2\right].\]
Integrating this equation  between the limits $x=0$ and $x=1$, $x\equiv t/t_0$, and taking into account that the present value $H_0t_0\approx1$ we obtain
+
\[0\approx\int\limits_0^1q(x)dx.\]
+
At a very high redshift $q(x)>0$. There must be periods (or a period) of accelerated expansion in the history of the universe in order to turn the integral $\int_0^1q(x)dx$ to zero.
+
  
 
</p>  </div></div></div>
 
</p>  </div></div></div>

Latest revision as of 07:49, 21 October 2016


New on the matter creation model

New OC

As we have seen above, negative pressure is the key ingredient required to achieve the accelerated expansion. Such pressure naturally appears when the system deviates from the thermodynamic equilibrium. In particular, as was first pointed out by Zel'dovich [Ya. B. Zel'dovich, JETP Lett. 12, 307 (1970)], the negative pressure is generated in process of particle creation due to the gravitational field. Physically, one may think that the (classical) time varying gravitational field works like a 'pump' supplying energy to the quantum fields. Construction of the models including the essentially quantum process of particle creation is problematic because of the difficulties connected with its inclusion into the classical field equations of Einstein. One can however avoid those difficulties on the phenomenological level [J.Lima, arXiv 0807.3379].

We shall now consider an open thermodynamical system where the number of fluid particles $N$ is not preserved. So the particle conservation equations \[\dot N_{;\mu}=0\Rightarrow \dot n + \theta n=0\] is now modified as \[\dot n + \theta n=n\Gamma,\] where $N_{\mu}=nu_\mu$ is the particle flow vector, $u_\mu$ is the particle four velocity, $\theta=u^\mu_{;\mu}$ is the fluid expansion, (for the FLRW Universe, $\theta=u^\mu_{;\mu}=3H$), $\dot n = n_{,\mu}u^{mu}$ and $\Gamma$ stands for the rate of change of the particle number in a comoving volume $V$, $\Gamma>0$ indicates particle creation while $\Gamma<0$ means particle annihilation. Any non-zero $\Gamma$ will effectively behave as a bulk viscous pressure of the thermodynamical fluid and nonequilibrium thermodynamics should come into the picture.


Problem 1

problem id: 150_dp0

Express the first law of thermodynamics in terms of the specific quantities: the energy density $\rho$ and the particle number density $n$.

Problem 2

problem id: 150_dp0

In the context of particle creation, let us consider spatially flat FRLW model of the Universe as an open thermodynamical system. Obtain the Friedmann equations for the cosmic fluid having energy-momentum tensor (in this section we will use the metric $(-,+,+,+)$ \[T_{\mu\nu}=(\rho+p+\Pi)u_\mu u_\nu+(p+\Pi)g_{\mu\nu}.\] For cosmological models with particle creation the cosmic fluid may be considered as a perfect fluid and the dissipative term $\Pi$ as the effective bulk viscous pressure due to particle production. In other words, the cosmic substratum is not a conventional dissipative fluid, rather a perfect fluid with varying particle number.

Problem 3

problem id: 150_dp0

Show that for the adiabatic process the effective bulk viscous pressure due to particle production is entirely determined by the particle production rate.

Problem 4

problem id: 150_dp0

An equivalent way to see the derivation of the field equations (7)-(8) (???) is to consider the energy-momentum tensor in the Einstein field equations as a total energy momentum tensor \[T^{(eff)}_{\mu\nu}=T^{(\gamma)}_{\mu\nu}+T^{(c)}_{\mu\nu},\] where $T^{(\gamma)}_{\mu\nu}$ is the energy-momentum tensor for the fluid with the EoS $p = (\gamma- 1) \rho$, i.e. \[T^{(\gamma)}_{\mu\nu}=(\rho+p)u_\mu u_\nu+pg_{\mu\nu}.\] and \[T^{(c)}_{\mu\nu}=p_c\left(u_\mu u_\nu+g_{\mu\nu}\right)\] is the energy-momentum tensor which corresponds to the matter creation term. Show that energy-momentum tensor $T^{(eff)}_{\mu\nu}$ provides us with the matter creation pressure \[\Pi=-\frac\Gamma\theta(\rho+p),\] obtained in the previous problem.

Problem 5

problem id: 150_dp0

Find the effective parameter $w_{eff}$ of the EoS for the cosmological model with particle creation.

Problem 6

problem id: 150_dp0

Find the deceleration parameter $q$ for the cosmological model with particle creation.

Problem 7

problem id: 150_dp0

Obtain the equation for the scale factor in cosmological model with particle creation.

Problem 8

problem id: 150_dp0

Calculate the derivative $dS/dt$ and find $S(t)$ for the Creation Matter Model (CMM).

Problem 9

problem id: 150_dp0

Calculate the scale factor dependence of the Universe's temperature for the Creation Matter Model (CMM).

Problem 10

problem id: 150_dp0

Solve the previous problem using the thermodynamic identity \[T\left(\frac{\partial p}{\partial T}\right)_n+n\left(\frac{\partial\rho}{\partial n}_T\right)=\rho+p.\]

Problem 11

problem id: 150_dp0

Solve the previous problem for cosmological substance with $\partial p/\partial\rho=w=const$.

Problem 12

problem id: 150_dp0

Show that the matter creation process can be effectively described by the dimensionless ratio \[\beta=\frac\Gamma\theta=\frac\Gamma{3H}.\]

Problem 13

problem id: 150_dp0

How is the evolution of $a(t)$ affected by rate of change of the particle number $\Gamma$?

Problem 14

problem id: 150_dp0

Let us assume a radiation dominated Universe ($w=1/3$, $\Gamma=\Gamma_r$). The dynamics is determined by the ratio $\Gamma_r/(3H)$ (see previous problem). The particle production must be strongly suppressed, $\Gamma_r/(3H)\ll1$, when the Universe enters the radiation phase. The simplest formula satisfying such a criterion is linear, namely: $\Gamma_r/(3H)=H/H_I$, where $H_I$ is the inflationary expansion rate associated with the initial de Sitter expansion ($H\le H_I$). Find the Hubble parameter dependence on the scale factor $H(a)$ for the radiation phase.

Problem 15

problem id: 150_dp0

Using the result obtained in the previous problem, show that in the considering model the Universe evolves continuously from de Sitter stage to the standard radiation phase.

Problem 16

problem id: 150_dp0

Find the radiation energy density $\rho_r$ as function of scale factor for the rate of change of the particle number $\Gamma=\Gamma_r$, $\Gamma_r/(3H)=H/H_I$ ($H_I$ is the inflationary expansion rate).

Problem 17

problem id: 150_dp0

For cosmological model considered in the problems \ref{2016_07_p14}-\ref{2016_07_p16} find how did the cosmic temperature evolve at the very early stages?

Problem 18

problem id: 150_dp0

Solve Problem \ref{2016_07_p14} for a substance with EoS $p=w\rho=(\gamma-1)\rho$ and $\Gamma=\Gamma_0H^2$ with $\Gamma_0=3\beta/H_e$.

Problem 19

problem id: 150_dp0

Find the deceleration parameter for the model considered in previous problem.

Problem 20

problem id: 150_dp0

Let us now consider how the CMM (Creation matter model) works at later stages of the Universe's evolution in the transition from Einstein-de Sitter to a late time de Sitter stage. Due to the conservation of baryon number one have to take into account only the production rate of cold dark matter particles $\Gamma_{dm}$. What is the exact form of $\Gamma_{dm}$?

Problem 21

problem id: 150_dp0

Let us consider the dynamical system in the context of bulk viscosity. In that case, we need to replace the pressure as $p\to p-3\zeta H$, where $\zeta$ is the coefficient of bulk viscosity. At what choice of the particle production rate $\Gamma$ the creation matter model would be equivalent to the bulk viscosity model?

Dynamical analysis of the Creation Matter Model

Problem 22

problem id: 150_dp0

Since the equation \[\dot H=-\frac32\gamma H^2\left(1-\frac\Gamma{3H}\right)\] is a one-dimensional first order differential equation, hence, the dynamics is obtained from the study of its critical points (or, fixed points) in which $\dot H=0$. Show that the fixed points $H_*$ in which $H_*\ne0$ represent the de Sitter Universe.

Problem 23

problem id: 150_dp0

Find the fixed points and determine their nature for the case with the particle creation rate $\Gamma=n/H$, $n>0$.

Problem 24

problem id: 150_dp0

Find the fixed points for the case with the particle creation rate $\Gamma=-\Gamma_0+mH+n/H$, $n>0$.

Problem 25

problem id: 150_dp0

Perform dynamic analysis of the solutions, obtained in previous problem.

Problem 26

problem id: 150_dp0

Prove that the dynamics generated by the creation matter model is equivalent to that driven by the single scalar field.

Problem 27

problem id: 150_dp0

Find the scalar field $\varphi(H)$ and the potential $V(\varphi)$, satisfying the condition that the dynamics generated by the creation matter model is equivalent to the dynamics driven by a single scalar field.

Problem 28

problem id: 150_dp0

Solve the previous problem for case \[\Gamma=-\Gamma_0+mH+n/H.\]

Problem 29

problem id: 150_dp0

Prove equivalence of the dynamics generated by the creation matter model with the dynamics driven by a decaying vacuum energy density model.

Problem 30

problem id: 150_dp0

Show that for spatially flat Universe the inequalities $\Gamma\ll H$ and $\Lambda\ll H^2$ should be satisfied simultaneously.

Problem 31

problem id: 150_dp0

Let us consider the particular case $\Lambda(t)=const$. Show that for such choice the creation matter model corresponds to the standard LCDM.

Problem 32

problem id: 150_dp0

Calculate the cosmographic parameter $j(H)$ for the cosmological model with particle creation.

Problem 33

problem id: 150_dp0

Solve the previous problem using the cosmographical relation \[j=-\frac1H\frac{dq}{dt}+q(1+2q).\]

However, in connection with the particle production at the expense of gravitational field of the expanding Universe, we recall that long back ago, Zeldovich [Ya. B. Zeldovich, JETP Lett. 12, 307 (1970)] introduced some bulk viscosity mechanism which is responsible for particle production. However, later on, Lima and Germano [J. A. S. Lima and A. S. M. Germano, Phys. Lett. A 170, 373 (1992); For a more general macroscopic approach see also R. Silva, J. A. S. Lima and M. O. Calv~ao, Gen. Rel. Grav. 34, 865 (2002).] showed that although both the processes, namely the bulk viscosity mechanism by Zeldovich [Ya. B. Zeldovich, JETP Lett. 12, 307 (1970).] and the gravitational particle production produce the same dynamics of the Universe, but in principle they are completely different from a thermodynamical point of view. Since we describe the particle production, we would like also to note an analogy which exists between the models driven by the particle production and the models of Steady State Cosmology developed in [F. Hoyle and J. V. Narlikar, Proc. Roy. Soc. A 282, 191 (1964); ibdem, 290, 143 (1966); see also Action at a Distance in Physics and Cosmology, 1974 (New York, Freeman).].

Problem 34

problem id: 150_dp0

In the model with $\Gamma=-\Gamma_0+mH$, where $0\le m\le3$ and $-3H_0\le\gamma\le0$, find the time dependence of the scale factor.

Problem 35

problem id: 150_dp0

In the model considered in the previous problem, express age of the Universe in terms of the parameters $\Gamma_0$ and $m$.

Problem 36

problem id: 150_dp0

Find deceleration parameter in the model with matter creation rate equal to $\Gamma$.

Problem 37

problem id: 150_dp0

Show that in the creation matter model with $\Gamma=-\Gamma_0+mH$ in the case $\Gamma_0=0$ the expansion of Universe is always accelerated if $m>1$ and always decelerated if $m<1$.

Problem 38

problem id: 150_dp0

Find the deceleration parameter as a function of the redshift in the creation matter model.

Problem 39

problem id: 150_dp0

Find the redshift value $z_t$ corresponding to the transition from the decelerated expansion to the accelerated one in the creation matter model.

Problem 40

problem id: 150_dp0

Express age of the Universe in terms of the $z_t$ (see the previous Problem) in the creation matter model with $m=0$.

New on the Ricci scalar etc

New OC

Problem 1

problem id: 10_1

Show that the Ricci scalar can be written in the following form \[R=6H^2[q(t)-\Omega(t)].\]

Problem 2

problem id: 10_2

Show that for a one-component flat Universe filled with ideal fluid with the EoS $p=w\rho$ the Ricci scalar can be written in the following form \[R=8\pi G\rho(3w-1).\]

Problem 3

problem id: 10_3

Find time derivative of the scalar curvature for the spatially flat Universe filled with a barotropic fluid.

Problem 4

problem id: 10_4

Compare entropy of the Sun in the present state and right after its transformation to a black hole due to compression.

Problem 5

problem id: 10_5

Show that black holes have negative thermal capacity.

Problem 6

problem id: 10_6

What are the physical reasons for the negative thermal capacity in Newtonnian self-gravitating systems?

Problem 7

problem id: 10_7

How the negative thermal capacity of the stars (in particular the Sun) affects the course of the nuclear reactions inside them?

Problem 8

problem id: 10_8

From the three parameters (mass, charge and angular momentum) of the Kerr-Newman black hole form all possible length scales.

Problem 9

problem id: 10_9

Following [O. Luongo and H. Quevedo, Self-accelerated universe induced by repulsive e?ects as an alternative to dark energy and modi?ed gravities, arXiv: (1507.06446)] let us introduce the parameter $\lambda=-\frac{\ddot a}{a}=qH^2$, so that $\lambda<0$ when the Universe is accelerating, whereas for $\lambda>0$ the Universe decelerates. Luongo and H. Quevedo showed, that the parameter $\lambda$ can be considered as an eigenvalue of the curvature tensor defined in special way. In particular, for FLRW metric the curvature tensor $R$ can be expressed as a ($6\times6$)-matrix \[R=diag(\lambda,\lambda,\lambda,r,r,r),\quad r\equiv\frac13\rho.\] The curvature eigenvalues reflect the behavior of the gravitational interaction and if gravity becomes repulsive in some regions, the eigenvalues must change accordingly; for instance, if repulsive gravity becomes dominant at a particular point, one would expect at that point a change in the sign of at least one eigenvalue. Moreover, if the gravitational field does not diverge at infinity, the eigenvalue must have an extremal at some point before it changes its sign. This means that the extremal of the eigenvalue can be interpreted as the onset of repulsion. Obtain the onset of repulsion condition in terms of cosmographic parameters.

Problem 10

problem id: 10_10

Represent result of the previous problem in terms of the Hubble parameter and its time derivatives.

New on the deceleration parameter

New OC

Problem 1

problem id: 150_dp0

Show that the fact that in the present time $t_0$ the relation $H_0t_0\approx1$ holds unavoidably follows existance of the accelerated expansion stage of the Universe evolution.

Problem 2

problem id: 150_dp2

Suppose that $dq/dt=f(q)$. Find the Hubble parameter in terms of $q$.


Problem 3

problem id: 150_dp2

Find connection between the scalar curvature and deceleration parameter for the flat Universe.


Problem 4

problem id: 150_dp3

Show that the deceleration parameter $q$ relates the density of the Universe $\rho$ to the critical density $\rho_{cr}$ through \[q=\frac12(1+3w)\frac\rho{\rho_{cr}}.\]


Problem 5

problem id: 150_dp4

Solve the previous problem using the second Friedmann's equation.

New on cosmology with power and hybrid expansion laws

Let us consider a general class of power-law cosmology described by the scale factor \[a(t)=a_0\left(\frac t{t_0}\right)^\alpha\] where $\alpha$ is a dimensionless positive parameter.

Problem 6

problem id: 150_dp5

Obtain deceleration parameter in the power-law cosmology.

Problem 7

problem id: 150_dp6

Express the scale factor $a(t)$ and Hubble parameter $H(z)$ in terms of the deceleration parameter.


Problem 8

problem id: 150_dp7

Obtain the comoving distance $r(z)$ in the power-law cosmology.


Problem 9

problem id: 150_dp8

Find the comoving distance in the Milne model.


Problem 10

problem id: 150_dp9

In the power-law cosmology, find time dependence of the CMB temperature in terms of the deceleration parameter.

Let us consider now a simple generalization of power-law cosmology, called the hybrid expansion law [O. Akarsu et. al. Cosmology with hybrid expansion law: scalar field reconstruction of cosmic history and observational constraints, arXiv:gr-qc/1307.4911] \[a(t)=a_0\left(\frac t{t_0}\right)^\alpha e^{\beta\left(\frac t{t_0}-1\right)}\] where $\alpha$ and $\beta$ are non-negative constants.

Problem 11

problem id: 150_dp10

Find Hubble parameter, deceleration parameter and jerk parameter for the hybrid expansion law.


Problem 12

problem id: 150_dp11

Find asymptotes of the scale factor and of the cosmographic parameters for the hybrid expansion law at $t\to0$ and $t\to\infty$.

Problem 13

problem id: 150_dp12

Find the moment of time $t_{tr}$ when the transition from deceleration to acceleration takes place.


Problem 14

problem id: 150_dp13

Find the range of the parameter $\alpha$ variation for the hybrid expansion law.


Problem 15

problem id: 150_dp14

In 1983, Berman [M. Berman, A special law of variation for Hubble's parameter. Nuovo Cimento B 74, 182 (1983)] proposed a special law of variation of Hubble parameter in FLRW space-time, which yields a constant value of DP, \[H=Da^{-n}.\] Find time dependence of the scale factor providing the constant deceleration parameter.


Problem 16

problem id: 150_dp15

For Berman's law of variation for Hubble's parameter \(H=Da^{-n}\) find the deceleration parameter and analyze what values of the parameter $n$ correspond to accelerated expansion, and which --- to decelerated one.

Linearly varying deceleration parameter

Inspired by O. Akarsu et al. Probing kinematics and fate of the Universe with linearlytime-varying deceleration parameter, arXiv:gr-qc/1305.5190 A general approach is to expand the deceleration parameter in Taylor series is \[q(x)=q_0+q_1\left(1-\frac x{x_0}\right)+q_2\left(1-\frac x{x_0}\right)^2+\ldots\] where $x$ is a some cosmological parameter as cosmic scale factor $a$, cosmic redshift $z$, cosmic time $t$ etc. As the first step one can take the following linear approximation \[q(x)=q_0+q_1\left(1-\frac x{x_0}\right).\]

Problem 17

problem id: 150_dp16

Consider the linearly varying deceleration parameter in terms of cosmic redshift $z$. Analyze advantages and defeats of such interpretation.


Problem 18

problem id: 150_dp17

Unlimited growth of the deceleration parameter for large $z$ in the parametrization, used in the previous problem, forces us to consider the linearly varying deceleration parameter in terms of scale factor. Make transition to such parametrization.


Problem 19

problem id: 150_dp19

Treating the Universe as a dynamical system it is useful to consider the parametrization of the deceleration parameter directly in terms of cosmic time $t$. Make transition to such parametrization.


Problem 20

problem id: 150_dp19

Consider a single-component flat Universe with a fluid described by an EoS parameter expressed as a first order Taylor expansion in cosmic time: \[w=w_0+w_1(1-t),\] where $w_0$ and $w_1$ are real constants and $t$ is the normalized time. Find the corresponding parametrization for the deceleration parameter.


Problem 21

problem id: 150_dp20

In a Universe filled with a fluid characterized by the EoS parameter $w=w_0+w_1(1-t/t_0)$, find time dependence of the scale factor. [S. Kumar, Probing the matter and dark energy sources in a viable Big Rip model of the Universe, arXiv:1404.1910]


Problem 22

problem id: 150_dp21

In a Universe filled with a fluid characterized by the EoS parameter $w=w_0+w_1(1-t/t_0)$, find the Hubble parameter, deceleration parameter and the jerk parameter.

Problem 23

problem id: 150_dp22

In a Universe filled with a fluid characterized by the EoS parameter $w=w_0+w_1(1-t/t_0)$, find the pressure and energy density.


Problem 24

problem id: 150_dp23

Show that the Universe filled with a fluid characterized by the EoS parameter $w=w_0+w_1(1-t/t_0)$ achieves de Sitter phase ($q=-1$) at the end of its half life.

New to Quantum cosmology

Problem 1

problem id: 150_dp24

Find the solutions corrected by LQC Friedmann equations for a matter dominated Universe.


New to Observational Cosmology

New OC

Problem 1

problem id: 150_o1

Find ratio between the illuminance of the Earth surface by the Sun and other light sources in the Universe (it is a quantitative resolution of the Olbers paradox).


Problem 2

problem id: 150_o2

Consider light emitted by surface of the Sun and observed on the Earth. Determine the shift of the observed frequency compared to the analogues light frequency emitted by atoms on the Earth.

Problem 3

problem id: 150_o3

Find the distance to a distant galaxy using recessional velocity as measured by the Doppler redshift.

Problem 4

problem id: 150_o4

Show that galaxies situated outside the Hubble sphere have imaginary redshift.

Problem 5

problem id: 150_o5

Find the exact relativistic Doppler velocity-redshift relation.

Problem 6

problem id: 150_o7

Find $V_{exp}(z)$ for three cosmological models: Einstein-de Sitter, Milne and de Sitter.

Problem 7

problem id: 150_o8

Some luminous object has apparent stellar magnitude $m=20$, and the absolute one is $M=-15$. Determine distance to it.

Problem 8

problem id: 150_o9

Due to a random coincidence the Balmer series of singly ionized helium atom in a distant star overlap with the Balmer series of hydrogen in the Sun. How fast this star recedes from us?

Problem 9

problem id: 150_o10

A gas cloud rotates around a super-massive black hole with mass equal to $M=3.6\times10^6M_\odot$ (it is a possible interpretation of recent observations). Assuming that distance between these objects is of order of $60$ lightyears, determine the expected Doppler shift.

Problem 10

problem id: 150_o11

Explain why transition from the cosmological time to the redshift can be considered as a quantitative characteristic describing evolution of the Universe.

Problem 11

problem id: 150_o12

Show, that for small separations the "Hubble law" $dz=Hdr$ holds. In other words, we have replaced the velocity by the redshift.

Problem 12

problem id: 150_o13

Solve the previous problem using the Taylor series of the scale factor in terms of time.

Problem 13

problem id: 2501_02o

Why the Linear Distance-Redshift Law in Near Space?

Problem 14

problem id: 150_o1

Determine physical distance to an object emitted light with redshift $z$ in the flat expanding Universe.

Problem 15

problem id: 150_o16

Find comoving distance to a presently observed galaxy as a function of the redshift.

Problem 16

problem id: 150_o17

Solve the previous problem for the flat Universe dominated by non-relativistic matter.

Problem 17

problem id: 150_o18

Determine the recession velocity due to the cosmological expansion for an object emitted light with the redshift $z$ in the flat expanding Universe.

Problem 18

problem id: 150_o19

Obtain relations between velocity of cosmological expansion and redshift.

Problem 19

problem id: 150_o20

Find dependence of the Hubble parameter $H$ on the redshift $z$ for the case of one-component Universe composed of the non-relativistic matter.

Problem 20

problem id: 150_o21

Find dependence of relative density $\Omega$ on the redshift $z$ for the case of one-component Universe composed of the non-relativistic matter.

Problem 21

problem id: 150_o22

Construct a scheme to determine sign of acceleration of the scale factor, based on measurement of the supernovae bursts characteristics.

Problem 22

problem id: 150_o23

Show that the time derivative of the redshift for the light emitted at time $t$ and registered at time $t_0$ can be determined as \[\dot z\equiv\frac{dz}{dt_0}=H(t_0)(1+z)-H(t).\]

Problem 23

problem id: 150_o24

In a flat one-component Universe with the state equation $p=w\rho$ at time $t_0$ one registers a light signal with redshift $z$. What values of the EoS parameter $w$ lead to \[\frac{dz}{dt_0}>0?\] Explain physical sense of the obtained result.

Problem 24

problem id: 150_o26

Show that the luminosity distance can be generally presented as \[d_L=\frac{1+z}{H_0\sqrt{\Omega_{k0}}}\sinh\left(H_0\sqrt{\Omega_{k0}}\int\limits_0^z\frac{dz}{H(z)}\right),\] where $\Omega_{k0}$ is the relative contribution of the spatial curvature.

Problem 25

problem id: 150_o27

\it When galaxy formation started in the history of the Universe remains unclear. Studies of the cosmic microwave background indicate that the Universe, after initial cooling (following the Big Bang), was reheated and reionized by hot stars in newborn galaxies at a redshift in the range $6<z<14$. Here we report a spectroscopic redshift of $z=6.96$. [M. Iue et al., A galaxy at a redshift z=6.96 , arXiv:0609393]

\bf Estimate age of the Universe when the first galaxies were created.

Problem 26

problem id: 150_o30

Express probability to find two galaxies in the infinitesimally small volumes $dV_1$ and $dV_2$ in terms of the correlation function $\xi(\vec{r_1},\vec{r_2})$, if average density of the galaxies in the considered volume equals $\bar n$ and total number of the galaxies is $N$. Point out main properties of the correlation function.

Problem 27

problem id: 150_o31

\bf Find relation between the spatial $\xi(r)$ and angular $w(\theta)$ correlation functions.

Problem 28

problem id: 150_o32

Show that power decay law of the spatial correlations leads to power decay law for the angular correlations.

Problem 29

problem id: 150_o33

Find two-point correlation function for $N$ galaxies distributed on a straight line with average density $\bar n$ in non-overlapping clusters of length $a$. Density of galaxies inside the cluster is constant and equals to $n_c$. The clusters are randomly distributed.

Problem 30

problem id: 150_o34

In a flat matter-dominated Universe of age $t_0$ light from a certain galaxy exhibits a redshift $z=0.95$. How long has it taken the light signal to reach us from this galaxy?

New from June 2015

New from June 2015

To Chapter 2

Problem 1

problem id: 150_0

Comprehensive explanations of the expanding Universe often use the balloon analogy. Although the balloon analogy is useful, one must guard against misconceptions that it can generate. Point out the misconceptions that appear when using this analogy [see M.O. Farooq, Observational constraints on dark energy cosmological model parameters, arXiv: 1309.3710.]


Problem 2

problem id: 150_1

(into the cosmography and extended deceleration parameter) Show that \[\frac{d\dot a}{da}=-Hq.\]


Problem 3

problem id: 150_2

Give a physical interpretation of the conservation equation.


Problem 4

problem id: 150_04

Find evolution equation for the density parameter $\Omega$ of the single-fluid FLRW models with the linear equation of state $p=w\rho$.


Problem 5

problem id: 150_05

Solve the previous problem for the multi-component case.


Problem 6

problem id: 150_06

Use the conformal time to prove existence of smooth transition from the radiation-dominated era to the matter dominated one.


Problem 7

problem id: 150_07

Consider a set of the cosmographic parameters built from the Hubble parameter and its time derivatives [see S. Carloni, A new approach to the analysis of the phase space of f(R)-gravity, arXiv:1505.06015) ] \[Q\equiv\frac{\dot H}{H^2},\quad J\equiv\frac{\ddot H H}{\dot H^2},\quad S\equiv\frac{\dddot H H^2}{\dot H^3},\ldots \] Express them in terms of the canonic cosmographic parameters $q,j,s\dots$.


Problem 8

problem id: 150_08

Consider another set of the cosmographic parameters [see S. Carloni, A new approach to the analysis of the phase space of f(R)-gravity, arXiv:1505.06015) ] \[\bar Q\equiv\frac{H_{,N}}{H},\quad \bar J\equiv\frac{H_{,NN}}{H},\quad \bar S\equiv\frac{H_{,NNN}}{H},\ldots,\] where \[H_{,N}\equiv \frac{dH}{d\ln a}.\] Express them in terms of the Hubble parameter and its time derivatives.


Problem 9

problem id: 150_09

Express the Ricci scalar and its time derivatives in terms of the $\bar Q$, $\bar J$ and $\bar S$.


Problem 10

problem id: 150_3

Show that for a perfect fluid with EoS $p=w(a)\rho$ the adiabatic sound speed can be represented in the form \[c_S^2=w(a)-\frac13\frac{d\ln(1+w)}{d\ln a}.\]


Problem 11

problem id: 150_4

Obtain equation for $\ddot\rho(t)$, where $\rho(t)$ is energy density of an ideal fluid participating in the cosmological expansion.


Problem 12

problem id: 150_5

Show that in the case of the flat Friedmann metric, the third power of the scale factor $\varphi\equiv a^3$ satisfies the equation \[\frac{d^2\varphi}{dt^2}=\frac32(\rho-p)\varphi,\quad 8\pi G=1.\] Check validity of this equation for different cosmological components: non-relativistic matter, cosmological constant and a component with EoS $p=w\rho$.


Problem 13

problem id: 150_6

The lookback time is defined as the difference between the present day age of the Universe and its age at redshift $z$, i.e. the difference between the age of the Universe at observation $t_0$ and the age of the Universe, $t$, when the photons were emitted. Find the lookback time for the Universe filled with non-relativistic matter, radiation and a component with the EoS $p=w(z)\rho$.


Problem 14

problem id: 150_7

Show that the Hubble radius grows faster than the expanding Universe in the case of power law expansion $a(t)\propto t^\alpha$ with $\alpha<1$ (the decelerated expansion).


To chapter 3

Problem 15

problem id: 150_015

Show, that in the Milne Universe the age of the Universe is equal to the Hubble time.


To chapter 4 The black holes

Problem 16

problem id: new2015_1

see E. Berti, A Black-Hole Primer: Particles, Waves, Critical Phenomena and Superradiant instabilities (arXiv:1410.4481[gr-qc])

A Newtonian analog of the black hole concept is a so-called "dark star". If we consider light as a corpuscle traveling at speed $c$, light cannot escape to infinity whenever $V_{esc}>c$, where \[V_{esc}^2=\frac{2GN}R.\] Therefore the condition for existence of "dark stars" in Newtonian mechanics is \[\frac{2GN}{c^2R}\ge1.\]

Can this condition be satisfied in the Newtonian mechanics?


Problem 17

problem id: 150_017

Derive the relation $T\propto M^{-1}$ from the Heisenberg uncertainty principle and the fact that the size of a black hole is given by the Schwarzschild radius.


To chapter 8

Problem 18

problem id: 2501_06

Why the cosmological constant cannot be used as a source for inflation in the inflation model?


Problem 19

problem id: 2501_09

Show that inflation ends when the parameter \[\varepsilon\equiv\frac{M_P^2}{16\pi}\left(\frac{dV}{d\varphi}\frac1V\right)^2=1.\]


Problem 20

problem id: 2501_10

How does the number of e-folds $N$ depend on the slow-roll parameter $\varepsilon$?


To chapter 9

Problem 21

problem id: 150_021

Using the by dimensional analysis for cosmological constant $\Lambda > 0$, define the set of fundamental "de Sitter units" of length, time and mass.


Problem 22

problem id: 150_022

In the Newtonian approximation, find the force acting on the point unit mass in the Universe filled by non-relativistic matter and cosmological constant. (see Chiu Man Ho and Stephen D. H. Hsu, The Dark Force: Astrophysical Repulsion from Dark Energy, arXiv: 1501.05592)


Problem 23

problem id: 150_023

Consider a spatially flat FLRW Universe, which consists of two components: the non-relativistic matter and the scalar field $\varphi$ in the potential $V(\varphi)$. Find relation between the scalar field potential and the deceleration parameter.


Problem 24

problem id: 150_024

Find relation between the deceleration parameter and the derivative $d\varphi/dz$ for the Universe considered in the previous problem.


Problem 25

problem id: new_30

Find the sound speed for the modified Chaplygin gas with the state equation \[p=B\rho-\frac A{\rho^\alpha}.\]


A couple of problems for the SCM

Problem 26

problem id: 150_026

Let $N=\ln(a/a_0)$, where $a_0=a(t_0)$ and $t_0$ is some chosen reference time. Usually the reference time is the present time and in that case $\tau=-\ln(1+z)$. Find $\Omega_m(N)$ and $\Omega_\Lambda(N)$ for the SCM.


Problem 27

problem id: 150_027

Express the cosmographic parameters $H,q,j$ as functions of $N=\ln a/a_0$ for the SCM.


Cardassian Model

[K. Freese and M. Lewis, Cardassian Expansion: a Model in which the Universe is Flat, Matter Dominated, and Accelerating, arXiv: 0201229] Cardassian Model is a modification to the Friedmann equation in which the Universe is flat, matter dominated, and accelerating. An additional term, which contains only matter or radiation (no vacuum contribution), becomes the dominant driver of expansion at a late epoch of the universe. During the epoch when the new term dominates, the universe accelerates. The authors named this period of acceleration by the Cardassian era. (The name Cardassian refers to a humanoid race in Star Trek whose goal is to take over the universe, i.e., accelerated expansion. This race looks foreign to us and yet is made entirely of matter.) Pure matter (or radiation) alone can drive an accelerated expansion if the first Friedmann equation is modified by the addition of a new term on the right hand side as follows: \[H^2=A\rho+B\rho^n,\] where the energy density $\rho$ contains only ordinary matter and radiation, and $n<2/3$. In the usual Friedmann equation $B=0$. To be consistent with the usual result, we take \[A=\frac{8\pi}{3M_{Pl}^2},\] where $M_{Pl}^2\equiv1/G$.



Problem 28

problem id: 150_cardas1

Show that once the new term dominates the right hand side of the Friedmann equation, we have accelerated expansion.


Problem 29

problem id: 150_cardas2

Let us represent the Cardassian model in the form \[H^2\propto \rho+\rho_X,\quad\rho_X=\rho^n.\] Find the parameter $w_X$ of the EoS $p_X=w_X\rho_X$, assuming that the Universe is filled exclusively by the non-relativistic matter.


Problem 30

problem id: 150_cardas3

Show that the result obtained in the previous problem takes place for arbitrary one-component fluid with $w_X=const$.


Problem 31

problem id: 150_cardas4

Generalize the previous problem for the case of two-component ideal liquid (non-relativistic matter $+$ radiation) with density $\rho=\rho_m+\rho_r$.


Problem 32

problem id: 150_cardas5

Show that we can interpret the Cardassian empirical term in the modified Friedmann equation as the superposition of a quintessential fluid with $w=n-1$ and a background of dust.


Problem 33

problem id: 150_cardas6

We have two parameters in the original Cardassian model: $B$ and $n$. Make the transition $\{B,n\}\to\{z_{eq},n\}$, where $z_{eq}$ is the redshift value at which the second term $B\rho^n$ starts to dominate.


Problem 34

problem id: 150_cardas7

What is the current energy density of the Universe in the Cardassian model? Show that the corresponding energy density is much smaller than in the standard Friedmann cosmology, so that only matter can be sufficient to provide a flat geometry.


Problem 35

problem id: 150_cardas8

Let us represent the basic relation of Cardassian model in the following way \[H^2=A\rho\left[1+\left(\frac\rho{\rho_{car}}\right)^{n-1}\right],\] where $\rho_{car}=\rho(z_{eq})$ is the energy density at which the two terms are equal. Find the function $\rho(z_{eq})$ under assumption that the Universe is filled with non-relativistic matter and radiation.


Problem 36

problem id: 150_cardas9

Let Friedmann equation is modified to be \[H^2=\frac{8\pi G}{3}g(\rho),\] where $\rho$ consists only of non-relativistic matter. Find the effective total pressure.


Problem 37

problem id: 150_cardas10

Find the speed of sound in the Cardassian model. [P.Gandolo, K. Freese, Fluid Interpretation of Cardassian Expansion, 0209322 ]


Problem 38

problem id: 150_cardas11

Find the deceleration parameter for the canonic Cardassian model.


Models with Cosmic Viscosity

A Universe filled with a perfect fluid represents quite a simple which seems to be in good agreement with cosmological observations. But, on a more physical and realistic basis we can replace the energy-momentum tensor for the simplest perfect fluid by introducing cosmic viscosity. The energy momentum tensor with bulk viscosity is given by \[T_{\mu\nu}=(\rho=p-\xi\theta)u_\mu u_\nu+(p-\xi\theta)g_{\mu\nu},\] where $\xi$ is bulk viscosity, and $\theta\equiv3H$ is the expansion scalar. This modifies the equation of state of the cosmic fluid. The Friedmann equations with inclusion of the bulk viscosity, i.e. using the energy-momentum tensor $T_{\mu\nu}$, read \begin{align} \nonumber \frac{\dot a^2}{a^2}&=\frac13\rho,\quad \rho=\rho_m+\rho_\Lambda,\quad 8\pi G=1;\\ \nonumber \frac{\ddot a^2}{a}&=-\frac16(\rho+3p-9\xi H). \end{align} Problems #150_8-#150_14 are inspired by A. Avelino and U. Nucamendi, Can a matter-dominated model with constant bulk viscosity drive the accelerated expansion of the universe? arXiv:0811.3253


Problem 39

problem id: 150_8

Consider a cosmological model in a flat Universe where the only component is a pressureless fluid with constant bulk viscosity ($\xi=const$). The pressureless fluid represent both the baryon and dark matter components. Find the dependence $\rho_m(z)$ for the considered model.


Problem 40

problem id: 150_9

Find $H(z)$ and $a(t)$ for the model of Universe considered in the previous problem.


Problem 41

problem id: 150_10

Analyze the expression for the scale factor $a(t)$ obtained in the previous problem for different types of the bulk viscosity.


Problem 42

problem id: 150_11

Show that the Universe in the considered model with $\xi=const$ had the Big Bang in the past for all values of the bulk viscosity in the interval $0<\bar\xi<3$ and determine how far in the past (in terms of the cosmic time) it happened.


Problem 43

problem id: 150_12

Show that the result of the previous Problem for zero bulk viscosity ($\xi=0$) correctly reproduces the lifetime of the matter-dominated Universe.


Problem 44

problem id: 150_13

As we have seen in the previous problems, in the interval $0<\bar\xi<3$ the Universe begins with a Big-Bang followed by an eternal expansion and this expansion begins with a decelerated epoch followed by an eternal accelerated one. The transition between the decelerated-accelerated expansion epochs depends on the value of $\bar\xi$. Find the value of the scale factor where the transition happens.


Problem 45

problem id: 150_14

Analyze the dependence \[a_t=\left(\frac{3-\bar\xi}{2\bar\xi}\right)^{2/3},\] obtained in the previous problem.


Problem 46

problem id: 150_15

Find the deceleration parameter $q(a,\bar\xi)$ for the cosmological model presented in the Problem #150_8.


Problem 47

problem id: 150_16

Analyze behavior of the deceleration parameter $q(a,\bar\xi)$ obtained in the previous problem for different values of the bulk viscosity $\bar\xi(\xi)$.


Problem 48

problem id: 150_17

Use result of the problem \ref{150_15} to find the current value of the deceleration parameter and make sure that for $\bar\xi=1$ the transition from the decelerated to accelerated epochs of the Universe takes place today.


Problem 49

problem id: 150_18

Find the curvature scalar $R(a,\xi)$ for the cosmological model presented in the Problem #150_8.


Problem 50

problem id: 150_19

Let us consider a flat homogeneous and isotropic Universe filled by a fluid with bulk viscosity. We shall assume that the EoS for the fluid is $p=w\rho$, $w=const$ and that the viscosity coefficient $\xi(\rho)=\xi_0\rho^\nu$. Find the dependence $\rho(a)$ for the considered model.


Problem 51

problem id: 150_0013

Usually the inflationary models of the early Universe contain two distinct phases. During the first phase entropy of the Universe remains constant. The second phase is essentially non-adiabatic, particles are produced through the damping of the coherent oscillations of the inflaton field by coupling to other fields and by its subsequent decay. Find relation between the bulk viscosity and the entropy production [J. A. S. Lima, R. Portugal, I. Waga, Bulk viscosity and deflationary universes, arXiv:0708.3280].




New from March 2015

New from March 2015


New from December 2014

New from December 2014


Exactly Integrable n-dimensional Universes

Exactly Integrable n-dimensional Universes


UNSORTED NEW Problems

The history of what happens in any chosen sample region is the same as the history of what happens everywhere. Therefore it seems very tempting to limit ourselves with the formulation of Cosmology for the single sample region. But any region is influenced by other regions near and far. If we are to pay undivided attention to a single region, ignoring all other regions, we must in some way allow for their influence. E. Harrison in his book Cosmology, Cambridge University Press, 1981 suggests a simple model to realize this idea. The model has acquired the name of "Cosmic box" and it consists in the following.

Imaginary partitions, comoving and perfectly reflecting, are used to divide the Universe into numerous separate cells. Each cell encloses a representative sample and is sufficiently large to contain galaxies and clusters of galaxies. Each cell is larger than the largest scale of irregularity in the Universe, and the contents of all cells are in identical states. A partitioned Universe behaves exactly as a Universe without partitions. We assume that the partitions have no mass and hence their insertion cannot alter the dynamical behavior of the Universe. The contents of all cells are in similar states, and in the same state as when there were no partitions. Light rays that normally come from very distant galaxies come instead from local galaxies of long ago and travel similar distances by multiple re?ections. What normally passes out of a region is reflected back and copies what normally enters a region.

Let us assume further that the comoving walls of the cosmic box move apart at a velocity given by the Hubble law. If the box is a cube with sides of length $L$, then opposite walls move apart at relative velocity $HL$.Let us assume that the size of the box $L$ is small compared to the Hubble radius $L_{H} $ , the walls have a recession velocity that is small compared to the velocity of light. Inside a relatively small cosmic box we use ordinary everyday physics and are thus able to determine easily the consequences of expansion. We can even use Newtonian mechanics to determine the expansion if we embed a spherical cosmic box in Euclidean space.


Problem 1

problem id:

As we have shown before (see Chapter 3): $p(t)\propto a(t)^{-1} $ , so all freely moving particles, including galaxies (when not bound in clusters), slowly lose their peculiar motion and ultimately become stationary in expanding space. Try to understand what happens by considering a moving particle inside an expanding cosmic box


Problem 2

problem id:

Show that at redshift $z=1$ , when the Universe is half its present size, the kinetic energy of a freely moving nonrelativistic particle is four times its present value, and the energy of a relativistic particle is twice its present value.


Problem 3

problem id:

Let the cosmic box is filled with non-relativistic gas. Find out how the gas temperature varies in the expanding cosmic box.


Problem 4

problem id:

Show that entropy of the cosmic box is conserved during its expansion.


Problem 5

problem id:

Consider a (cosmic) box of volume V, having perfectly reflecting walls and containing radiation of mass density $\rho $. The mass of the radiation in the box is $M=\rho V$ . We now weigh the box and find that its mass, because of the enclosed radiation, has increased not by M but by an amount 2M. Why?


Problem 6

problem id:

Show that the jerk parameter is \[j(t)=q+2q^{2} -\frac{\dot{q}}{H} \]


Problem 7

problem id:

We consider FLRW spatially flat Universe with the general Friedmann equations \[\begin{array}{l} {H^{2} =\frac{1}{3} \rho +f(t),} \\ {\frac{\ddot{a}}{a} =-\frac{1}{6} \left(\rho +3p\right)+g(t)} \end{array}\] Obtain the general conservation equation.


Problem 8

problem id:

Show that for extra driving terms in the form of the cosmological constant the general conservation equation (see previous problem) transforms in the standard conservation equation.


Problem 9

problem id:

Show that case $f(t)=g(t)=\Lambda /3$ corresponds to $\Lambda (t)CDM$ model.