Cosmological models based on an asymmetric scalar doublet with kinetic coupling of components. II. Numerical modeling^,

Yu.G. Ignat’ev¹
Institute of Physics, Kazan Federal University, Kremlyovskaya str., 16A, Kazan, 420008, Russia

I.A. Kokh²
N.I. Lobachevsky Institute of Mathematics and Mechanics, Kazan Federal University, Kremlyovskaya str., 35, Kazan, 420008, Russia

Numerical modeling of a mathematical model of the cosmological evolution of an asymmetric scalar doublet with kinetic interaction between the components was carried out. A wide range of values of fundamental parameters and initial conditions of the model are considered. Various types of behavior have been identified: models with an infinite inflationary past and future - with and without a rebound point, models with a finite past and infinite future, with an infinite past and finite future (Big Rip), as well as models with a finite past and future. Based on numerical analysis, the behavior of models near the initial singularity and the Big Rip is studied; it is shown that in both cases the barotropic coefficient tends to unity, which corresponds to an extremely rigid state of matter near singularities. A numerical example of the cosmological generation of the classical component of a scalar doublet by its phantom component is given. An assessment was made of the creation of the velocity of fermion pairs by a scalar field near the rebound points and it was shown that a scalar field at the cold stage of the Universe can ensure the creation of the required number of massive scalarly charged fermions.
Keywords: cosmological model, phantom and classical scalar fields, quality analysis, asymptotic behavior, numerical modelling, scalar field generation, types of behavior.

1 Introduction

In the first part of the article [1] a cosmological model based on an asymmetric scalar doublet with kinetic coupling of components was formulated and its main properties were investigated. In addition, this work demonstrated an example of numerical modeling for a special case of a set of fundamental constants and initial conditions, illustrating the analytical properties of the model.

In this part of the article we will present a wider range of numerical modeling results and their analysis in order to identify types of behavior of the model, as well as solve the problem of generating one of the components of a scalar doublet by another, similar to cosmological models with scalarly charged fermions [2], [3], [4]. As we noted in [1], we will also pursue the goal of possibly replacing the mathematically complicated and cumbersome model of scalarly charged fermions with a simpler mathematical model in the theory of scalar-gravitational instability.

We present the necessary information from the first part of the article [1].

The dynamical system corresponding to the cosmological model in the spatially flat Friedmann metric with the scale factor $a(t)$ consists of a normal system of ordinary differential equations:

$$\dot{\xi}=H\qquad(a\equiv\exp(\xi));$$ (1)

$$\dot{\Phi}=Z;$$ (2)

$$\dot{Z}=-3HZ-\frac{m^{2}\Phi-\alpha\Phi^{3}}{1+\gamma^{2}}-\gamma\frac{\mu^{2}\varphi-\beta\varphi^{3}}{1+\gamma^{2}};$$ (3)

$$\dot{\varphi}=z;$$ (4)

$$\dot{z}=-3Hz+\frac{\mu^{2}\varphi-\beta\varphi^{3}}{1+\gamma^{2}}-\gamma\frac{m^{2}\Phi-\alpha\Phi^{3}}{1+\gamma^{2}};$$ (5)

$$\dot{H}=-\frac{1}{2}Z^{2}+\frac{1}{2}z^{2}-\gamma Zz,$$ (6)

and its first zero integral - the equations of the Einstein - Higgs hypersurface of the autonomous subsystem $S_{1}=$ {(2) – (6)} dynamic system $S_{0}=$ {(1) – (6)} ($S_{1}\subset S_{0}$):

$$\begin{array}[]{l}3H^{2}-\dfrac{Z^{2}}{2}+\dfrac{\alpha\Phi^{4}}{4}-\dfrac{m^{2}\Phi^{2}}{2}\\[11.0pt] +\dfrac{z^{2}}{2}+\dfrac{\beta\varphi^{4}}{4}-\dfrac{\mu^{2}\varphi^{2}}{2}-\gamma Zz-\Lambda=0,\end{array}$$ (7)

where $\Phi(t)$ and $\varphi(t)$ are the potentials of classical and phantom scalar fields, $H(t)$ is the Hubble parameter, $\alpha,\beta$ are the self-interaction constants of these fields, $m,\mu$ are the masses of their quanta, $\gamma$ is the kinetic interaction constant of the classical and phantom scalar fields.

$$\Lambda=\Lambda_{0}-\frac{m^{4}}{4\alpha}-\frac{\mu^{4}}{4\beta},$$ (8)

and $\Lambda_{0}$ is the seed value of the cosmological constant.

Since both dynamical systems, $S_{0}$ and $S_{1}$, are autonomous, in cases where we do not need knowledge of the scale factor function $a(t)$, we will study the autonomous subsystem $S_{1}$, otherwise - the system $S_{0}$. Let us note in passing that, in principle, in all cases it would be possible to limit ourselves to studying the system $S_{1}$, since the scale factor $a(t)$ using the equation (1) could be found using integration

$$a(t)=\mathrm{e}^{\int H(t)dt},$$

however, from a technical point of view, it is easier to find a numerical solution to a normal system of 6 differential equations $S_{0}$ than to find the integral of the function $H(t)$ based on a numerical solution to a normal system of 5 differential equations $S_{1}$ using numerical methods.

Further, the invariant cosmological acceleration $\Omega(t)$ is related to the Hubble parameter $H(t)$ by the relation

$$\Omega=\frac{a\ddot{a}}{{\dot{a}}^{2}}\equiv 1+\frac{\dot{H}}{H^{2}}=-\frac{1}{2}(1+3\kappa),$$ (9)

where $\kappa(t)$ is the barotropic coefficient of cosmological matter:

$$\kappa=P\diagup E\equiv-\frac{1}{3}(1+2\Omega),$$ (10)

and $E$ is normalized effective energy and $P$ is normalized effective pressure:

$$E= \displaystyle\frac{Z^{2}}{2}-\frac{\alpha\Phi^{4}}{4}+\frac{m^{2}\Phi^{2}}{2}\hskip 56.9055pt$$ (11)

$$\displaystyle-\frac{z^{2}}{2}-\frac{\beta\varphi^{4}}{4}+\frac{\mu^{2}\varphi^{2}}{2}+\gamma Zz+\Lambda,$$

$$P= \displaystyle\frac{Z^{2}}{2}+\frac{\alpha\Phi^{4}}{4}-\frac{m^{2}\Phi^{2}}{2}\hskip 56.9055pt$$ (12)

$$\displaystyle-\frac{z^{2}}{2}+\frac{\beta\varphi^{4}}{4}-\frac{\mu^{2}\varphi^{2}}{2}+\gamma Zz-\Lambda,$$

at that

$$E+P=Z^{2}-z^{2}+2\gamma Zz\equiv-2\dot{H}.$$ (13)

Along with the total effective energy and pressure, we introduce similar quantities for individual components of the cosmological system, normalized energy densities $E_{c},E_{f}$ and pressure $P_{c},P_{f}$

$$E_{c}=\frac{Z^{2}}{2}-\frac{\alpha\Phi^{4}}{4}+\frac{m^{2}\Phi^{2}}{2};$$ (14)

$$E_{f}=-\frac{z^{2}}{2}-\frac{\beta\varphi^{4}}{4}+\frac{\mu^{2}\varphi^{2}}{2};$$ (15)

$$P_{c}=\frac{Z^{2}}{2}+\frac{\alpha\Phi^{4}}{4}-\frac{m^{2}\Phi^{2}}{2};$$ (16)

$$P_{f}=-\frac{z^{2}}{2}+\frac{\beta\varphi^{4}}{4}-\frac{\mu^{2}\varphi^{2}}{2}$$ (17)

and the corresponding interaction values of their components:

$$E_{cf}=P_{cf}=\gamma Zz,$$ (18)

so that

$$E=E_{c}+E_{f}+E_{cf}+\Lambda;$$ (19)

$$P=P_{c}+P_{f}+P_{cf}-\Lambda.$$ (20)

A particular cosmological model M is determined, firstly, by a system of dynamic equations (1)–(6), secondly, by fundamental parameters

$$\textbf{P}=[[\alpha,\beta,m,\mu],\gamma,\Lambda]$$ (21)

and, thirdly, by initial conditions

$$\textbf{I}=[\Phi_{0},Z_{0},\varphi_{0},z_{0},e],\qquad(e=\pm 1).$$ (22)

Here the sign $e$ corresponds to the choice of a positive or negative solution to the equation (7) relative to the initial value of the Hubble parameter. Note that using the autonomy of the dynamical systems $S_{0}$ and $S_{1}$ and the invariance of the Friedman metric with respect to time translations $t\to t+t_{0}$, we can always choose the initial condition for the scale function $\xi(t)$ in the form ( see for example [6])

$$\xi(0)=0.$$ (23)

The cosmological model under study can be considered as a five-dimensional dynamical system in the arithmetic phase space $\mathds{R}^{5}=\{\Phi,Z,\varphi,z,H\}=\mathds{R}^{3}\cup\mathds{R}^{3}$, three-dimensional phase subspaces $\mathds{R}^{3}=\{\Phi,Z,H\}$ and $\mathds{R}^{3}=\{\varphi,z,H\}$ We will further denote by the symbols $\Sigma_{\Phi}$ and $\Sigma_{\varphi}$ and, for simplicity, call them the classical and phantom phase spaces, respectively. In this case, $\Sigma_{\Phi}\cup\Sigma_{\varphi}=\mathds{R}^{5}$ and $\Sigma_{\Phi}\cap\Sigma_{\varphi}=\mathds{R}^{1}=OH$.

The singular points of the dynamical system under study have coordinates (see[1]):

$$M_{0,0}^{\pm}\left(0,0,0,0,\pm\frac{\sqrt{3\Lambda}}{3}\right);$$ (24)

$$M_{0,\pm 1}^{\pm}:\left(0,0,\pm\frac{\mu}{\sqrt{\beta}},0,\pm\frac{\sqrt{3\Lambda_{\beta}}}{3}\right);$$ (25)

$$M_{\pm 1,0}^{\pm}:\left(\pm\frac{m}{\sqrt{\alpha}},0,0,0,\pm\frac{\sqrt{3\Lambda_{\alpha}}}{3}\right);$$ (26)

$$M_{\pm 1,\pm 1}^{\pm}:\left(\pm\frac{m}{\sqrt{\alpha}},0,\pm\frac{\mu}{\sqrt{\beta}},0,\pm\frac{\sqrt{3\Lambda_{0}}}{3}\right),$$ (27)

where $\Lambda_{\beta}\equiv\Lambda+\mu^{4}/4\beta$, $\Lambda_{\alpha}\equiv\Lambda+m^{4}/4\alpha$.

Before moving on to the results of numerical modeling, we make the following remarks:

Remarks I (preliminary).

I.1.   According to Property I.1 of the invariance of a dynamical system with respect to changes in the sign of $\gamma$, it is sufficient to study mathematical models with a positive value of $\gamma\geqslant 0$;

I.2.   According to Property I.2 when transforming similarity :

$${\bf\tilde{P}}=\left[\left[\frac{\alpha}{k^{2}},\frac{\beta}{k^{2}},\frac{m}{k},\frac{\mu}{k}\right],\gamma,\frac{\Lambda}{k^{2}}\right]$$ (28)

solution of the Cauchy problem with initial conditions

$${\bf\tilde{I}}=\left[\Phi_{0},\frac{Z_{0}}{k},\varphi_{0},\frac{z_{0}}{k},e\right]$$ (29)

is obtained from the solution $\mathbf{S}$ of the Cauchy problem for the preimage according to the rule

$$\textbf{S}(t)=[\Phi(t),Z(t),\varphi(t),z(t),H(t)]\Rightarrow$$ (30)

$$\!\!\!\mathbf{\tilde{S}}(t)=\left[\Phi\left(t\right),\frac{1}{k}Z\left(t\right),\varphi\left(t\right),\frac{1}{k}z\left(t\right),\frac{1}{k}H\left(t\right)\right],$$ (31)

where the connection between the time $\tilde{t}$ in the model ${\bf\tilde{M}}$ and the time $t$ of the preimage ${\bf M}$ is taken into account

$$\tilde{t}=kt;$$ (32)

I.3.   According to Property I.4, the coordinates of singular points of an autonomous dynamic system (2)–(6) with the integral condition (7), as well as their character, coincide with the coordinates of singular points and their character for a dynamic system without interaction of components [6];

I.4.   According to [7] during the similarity transformation (28) – (32) the eigenvalues $\lambda$ of the matrix of the dynamical system are transformed according to the rule

$$\tilde{\lambda}=\frac{\lambda}{k};$$ (33)

I.5.   All singular points (24) – (27) of the dynamical system $S_{0}$ are exact constant solutions of this system. Wherein

$$a(t)=e^{\pm H_{0}t},$$ (34)

where $H_{0}$ is the positive value of the 5th coordinate of the corresponding singular point from the list (24) – (27). Thus, all singular points correspond to inflationary solutions - inflationary expansion when choosing a positive sign in (34), or inflationary compression when choosing a negative sign. In both cases, these solutions correspond to cosmological models with an infinite past and an infinite future. At the same time, the question of the sustainability of these solutions remains open.

2 Basic model

Let’s set the parameters of the basic model as follows:

$$\mathbf{P_{0}}=[[1,1,1,1],\gamma,3\cdot 10^{-6}],$$ (35)

where the constant $\gamma$ is still an arbitrary parameter.

2.1 Singular points

The coordinates of the singular points of the model (35) and their characters are indicated in Tab.2.1. Here and in what follows, the following notations are used: $\mathbf{A}$ – attraction, $\mathbf{R}$ – repulsion, $\mathbf{S}$ – saddle, we also somewhat simplify the characteristics of singular points, preserving only their properties of attraction and repulsion.
Remarks II (about the character of the points).

The value of the parameter $\gamma$ does not affect either the coordinates of singular points or their character. In this case, however, the absolute values of the real and imaginary parts of the eigenvalues depend on the value of $\gamma$, but the signs of these parts do not change.

Tab. 1. Characters (type) of singular pointsin subspaces $[\Sigma_{\Phi},\Sigma_{\varphi}]$ in the model with parameters $\mathbf{P_{0}}$

Singular points	Coordinates	Type
$M_{0,0}^{+}$	$[0,0,0,0,0.001]$	$[\mathbf{A},\ \mathbf{S}]$
$M_{0,0}^{-}$	$[0,0,0,0,-0.001]$	$[\mathbf{R},\ \mathbf{S}]$
$M_{0,\pm 1}^{+}$	$[0,0,\pm 1,0,0.289]$	$[\mathbf{A},\ \mathbf{A}]$
$M_{0,\pm 1}^{-}$	$[0,0,\pm 1,0,-0.289]$	$[\mathbf{R},\ \mathbf{R}]$
$M_{\pm 1,0}^{+}$	$[\pm 1,0,0,0,0.289]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,0}^{-}$	$[\pm 1,0,0,0,-0.289]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,\pm 1}^{+}$	$[\pm 1,0,\pm 1,0,0.409]$	$[\mathbf{S},\ \mathbf{A}]$
$M_{\pm 1,\pm 1}^{-}$	$[\pm 1,0,\pm 1,0,-0.409]$	$[\mathbf{S},\ \mathbf{R}]$

Due to Remarks II, the results of Tab. 2.1 are applicable to models with arbitrary values of the $\gamma$.

2.2 System in steady state

Let us now specify the $\gamma$ parameter and consider a model with the following parameters:

$$\mathbf{P_{1}}=\bigl{[}[1,1,1,1],10^{-5},3\cdot 10^{-6}\bigl{]}.$$ (36)

According to Tab.2.1, there are only two singular points, $M_{0,\pm 1}^{+}$ (see (25)), which are simultaneously attractive in the subspaces $\Sigma_{\Phi}$ and $\Sigma_{\varphi}$.

In Fig.2.2 shows the evolution of geometric factors - the scale function $\xi(t)$ and the Hubble parameter $H(t)$ for a model with parameters $\mathbf{P_{1}}$ under initial conditions corresponding to the coordinates of stable points $M_{0,\pm 1}^{+}$

$$\mathbf{I^{\pm}_{0}}=[0,0,0,1,\pm 1].$$ (37)

Fig. 1.  Functions $\xi(t)$ are dashed and dashed-dotted, $H(t)$ are solid and long-dashed lines, respectively, for the initial conditions $\mathbf{I^{+}_{0}},\mathbf{I^{-}_{0}}$.

In this case, the initial conditions coincide with the coordinates of the stable singular point, the exact solution is the constant scalar potentials $\Phi_{0}=0,\,\varphi_{0}=\pm\mu/\sqrt{\beta}$. In this case, the model corresponds to either an infinite inflationary expansion ($H=H_{+}>0$) or an infinite inflationary contraction ($H=H_{-}<0$).

2.3 Rebound case

If the initial conditions do not coincide with the coordinates of stable singular points, a rebound point appears in the models – the minimum of the scale function $\xi(t)\Rightarrow H=0$. At this point, inflationary compression gives way to inflationary expansion with a symmetrical value of the Hubble constant parameter $H_{\pm}=\pm H_{0}$.

Fig. 2.  Functions $\xi(t)$ are dashed and dash - dotted, $H(t)$ are solid and long-dashed lines, respectively for initial conditions $\mathbf{I^{+}_{1}}$ and $\mathbf{I^{-}_{1}}$.

Fig. 3.  Functions $\xi(t)$ are dashed and dash - dotted, $H(t)$ are solid and long-dashed lines, respectively for initial conditions $\mathbf{I^{+}_{2}}$ and $\mathbf{I^{-}_{2}}$.

In Fig.2.3–2.3 shows the evolution of geometric factors - the scale function $\xi(t)$ and the Hubble parameter $H(t)$ for a model with parameters $\mathbf{P_{1}}$ under initial conditions, close to the coordinates of stable points $M_{0,\pm 1}^{+}$

$$\begin{array}[]{l}\mathbf{I^{\pm}_{1}}=[0,0,0,0.99,\pm 1],\\ \mathbf{I^{\pm}_{2}}=[0,0,0,1.01,\pm 1].\end{array}$$ (38)

Commenting on the graphs in Fig.2.2–2.3, note the following:

$\bullet$   All the cases considered above describe cosmological models with an infinite past and an infinite future.
$\bullet$   In this case, however, in the case of the initial conditions (37) coinciding with the coordinates of the stable singular point, the entire cosmological history is described by the inflationary solution corresponding limitless expansion or contraction of the Universe. In this case, the graphs of the functions $\xi(t)$ and $H(t)$ do not depend on the value of the parameter $\gamma$.
$\bullet$   At the slightest discrepancy between the initial conditions and the coordinates of the stable singular point $M_{0,\pm 1}^{+}$, the cosmological history breaks down, regardless of the initially given sign of the Hubble parameter into two stages: first - inflationary compression, then - inflationary expansion. At the same time, graphs of the geometric factors $\xi(t)$ and $H(t)$ corresponding to the initial conditions with positive and negative sign of the Hubble parameter, are similar with a shift along the time axis.

2.4 Evolution of scalar fields

If the initial conditions coincide with the coordinates of the stable singular point, the exact solution, as we indicated above, is the constant scalar potentials $\Phi_{0}=0,\varphi_{0}=\pm\mu/\sqrt{\beta}$. Therefore, we can talk about the evolution of scalar fields only when the initial state deviates from a stable singular point.

In Fig.2.4–2.5 shows graphs of the evolution of scalar potentials $\Phi(t)$ and $\varphi(t)$ for initial conditions $\mathbf{I^{+}_{1}}$ and $\mathbf{I^{+}_{2}}$. From these graphs it is clear that the potentials start from a position close to stable in the infinite past at the stage of inflationary compression, experience fluctuations near the minimum point of the scale factor and then return to an equilibrium position in the infinite future already at the stage of inflationary expansion.

Fig. 4.  Scalar potential $\Phi(t)$ – solid, $\varphi(t)$ – dashed line for initial conditions $\mathbf{I^{+}_{1}}$; $\gamma=10^{-5}$.

2.5 Effect of the constant $\gamma$

Note that if the initial conditions deviate from the coordinates of the points $M_{0,\pm 1}^{-}$, the value of the parameter $\gamma$ begins to influence the behavior of the scale functions $\xi(t)$ and $H(t)$. In Fig.2.5–2.5 just such an influence was demonstrated. As can be seen from these graphs, as the parameter $\gamma$ increases, the cosmological model acquires an initial singularity at a finite value of cosmological time (in our case, Fig.2.5, $t_{s}\approx-8.0421765$). At the same time, the Universe finds its beginning, lost by endless inflation in an unstable state.

Fig. 5.  Function $\xi(t)$ is dashed, $H(t)$ is solid line for initial conditions $\mathbf{I^{+}_{1}}$; $\gamma=1$.

Fig. 6.  Scalar potential $\Phi(t)$ – solid, $\varphi(t)$ – dashed line for initial conditions $\mathbf{I^{+}_{2}}$; $\gamma=10^{-5}$.

At $\gamma=1$, the scenario with the rebound point remains, but the graphs of $\xi(t)$ and $H(t)$ near the rebound point become more contrasting. At $\gamma=10$ the rebound point turns into a cosmological singularity.

In Fig.2.5 shows the influence of the value of the constant $\gamma$ on the evolution of scalar potentials. This case, as we answered above (see Fig.2.5 – 2.5), corresponds to a Universe with a rebound point, or with an initial singularity at $t_{s}\approx-8.0421765$.

The short-term generation of a classical field near the rebound point or singularity occurs precisely due to the kinetic connection of the classical and phantom fields.

Fig. 7.  Function $\xi(t)$ is dashed, $H(t)$ is solid line for initial conditions $\mathbf{I^{+}_{1}}$; $\gamma=10$.

Fig. 8.  Scalar potential $\Phi(t)$ – solid, $\varphi(t)$ – dashed line for initial conditions $\mathbf{I^{+}_{2}}$; $\gamma=1$.

In Fig.2.5 shows graphs of the evolution of the barotropic coefficient of cosmological matter $\kappa(t)$ (10) for the initial conditions $\mathbf{I^{+}_{1}}$ and $\mathbf{I^{+}_{2}}$ at $\gamma=10^{-5}$. In the first case, the model has a rebound point (Fig.2.5), in the second case, an initial singularity at time $t_{s}\approx-8.0421765$ (Fig.2.5). The system passes through the state $H=$, due to which singularities of the functions $\Omega(t)$ and $\kappa(t)$ arise at this point. Far from this critical point, the cosmological model is in the inflationary regime ($\kappa=-1\Rightarrow\Omega=1$). Among other things, you can see that the model approaches the singularity with the barotropic coefficient $\kappa=1$.

Fig. 9.  Functions of the barotropic coefficient $\kappa(t)$ (10): solid line for the initial conditions $\mathbf{I^{+}_{1}}$, dashed line – $\mathbf{I^{+}_{2}}$; $\gamma=1$.

3 Model with a finite past

The models with an initial singularity, i.e., with a finite past, discussed in the previous section are more likely demonstration examples than models of the real Universe. Indeed, the value of the kinetic interaction constant between the components $\gamma\gtrsim 1$ with values of the other fundamental constants in the set $\mathbf{P_{1}}$ of order 1 (36) is unacceptably large. Indeed, with large-scale transformations of such a model to a real model with masses of superheavy bosons of the order of $10^{17}G$, according to Remark I.1 we obtain for this model

$$\tilde{\mathbf{P}}=\bigl{[}\bigl{[}10^{-4},10^{-4},10^{-2},10^{-2}\bigr{]},1,3\cdot 10^{-10}\bigr{]}.$$

Thus, the kinetic interaction constant of the components turns out to be 4 orders of magnitude larger than the self-interaction constants of the classical $\alpha$ and phantom $\beta$ fields, which is physically unacceptable. The real value of the interaction constant must be less than or on the order of the value of the self-interaction constants of these fields.

Therefore, in the class of models with parameters $\mathbf{P_{1}}$ we consider especially the model with parameters

$$\mathbf{P^{*}_{1}}=[[1,1,1,1],10^{-5},3\cdot 10^{-6}]$$ (39)

and initial conditions

$$\mathbf{I_{1}}=[0.9,0,0.9,0,1].$$ (40)

3.1 Evolution of geometric factors

This model with an infinite inflationary future with the Hubble parameter $H_{+\infty}\approx 0.289$ (stable equilibrium point $M_{0,\pm 1}^{+}$, Tab.2.1) and a finite past, which corresponds to the singularity at time $t=t_{s}\approx-2.2284930$. Fig.3.1 shows the evolution of the scale function $\xi(t)=\ln(a(t))$ and the Hubble parameter $H(t)$ (1) for this model. The coordinates and characters of singular points obviously coincide with the corresponding values indicated in Tab.2.1.

On the graph of the function $H(t)$ Fig.3.1, two plateaus can be identified in the intervals $T_{1}\approx(-1,2)$ and $T_{2}\approx(5,\infty)$, which exactly correspond to singular points: $T_{1}$ – saddle points $M_{\pm 1,\pm 1}^{+}$, $T_{2}$ – attracting points $M_{0,\pm 1}^{+}$ and saddle points $M_{\pm 1,0}^{+}$. The $H(t)$ curve demonstrates the transition from an unstable state (saddle point $M_{\pm 1,\pm 1}^{+}$) through the first plateau to a stable state (attracting point $M_{0,\pm 1}^{+}$) - second plateau.

Fig. 10.  Evolution of the scale function $\xi(t)$ (dashed) and the Hubble parameter $H(t)$ (solid) in a model with a finite past at the point $t_{s}\approx-2.2284930$ with parameters $\mathbf{P^{*}_{1}}$ (39) and initial conditions $\mathbf{I_{1}}$ (40).

In Fig.3.1 the solid line shows the behavior of the barotropic coefficient $\kappa(t)$ near this singularity, calculated by the formula (10) using expressions for the total energy density (11) and pressure (12). Note that the value of this coefficient, calculated both from the indicated formulas and from the formula for cosmological acceleration (9), coincide. Further in the same figure, dashed and long-dashed lines show the behavior of the barotropic coefficients for the classical and phantom components, respectively; the dotted line $\kappa=1$ shows the behavior of the barotropic coefficient for the interaction component. According to this graph, at the singularity point the barotropic coefficient reaches the value $\kappa=1$, which corresponds to the total extremely rigid equation of state $P=E$ (10). Over time, the functions $\kappa(t)$ and $H(t)$ tend to $\kappa\to-1$, which corresponds to inflationary compression with the inflationary equation in the state $P=-E$.

Fig. 11.  Evolution of barotropic coefficients $\kappa(t)$ for various system components in a model with parameters $\mathbf{P^{*}_{1}}$ (39) and initial conditions $\mathbf{I_{1}}$ (40).

3.2 Evolution of scalar fields

In the future, we will be faced with the need to display graphs on significantly different scales. In order to overcome this problem, we will use a one-to-one scaling mapping if necessary (see [5]):

$$\rm{Lig}(x)=\mathrm{sgn}(x)\log[10](1+|x|);$$

$$\rm{Lig}^{-1}(x)=\mathrm{sgn}(x)\left(10^{|x|}-1\right).$$ (41)

Fig. 12.  Evolution of energy density and pressure.

Fig.3.2 shows graphs of the evolution of energy density and pressure of the classical field, according to their contributions to the expressions (11) and (12): $E_{c}(t)$ (solid line) and $P_{c}(t)$ (dashed line) of a classical scalar field in a model with a finite past at the point $t_{s}\approx-2.2284930$ with parameters $\mathbf{P^{*}_{1}}$ (39) and initial conditions $\mathbf{I_{1}}$ (40). Note that these quantities are singular at the point of cosmological singularity.

Fig.3.2 shows graphs of the evolution of the scalar classical $\Phi(t)$ and phantom $\varphi(t)$ potentials, and Fig.3.2 shows the derivatives of these potentials.

Fig. 13.  Evolution of $\Phi(t)$ (solid line) and $\varphi(t)$ (dashed line) in the model with parameters $\mathbf{P^{*}_{1}}$ (39) and initial conditions $\mathbf{I_{1}}$ (40).

Fig. 14.  Evolution of $Z(t)$ (solid line) and $z(t)$ (dashed line) in a model with parameters $\mathbf{P^{*}_{1}}$ (39) and initial conditions $\mathbf{I_{1}}$ (40).

4 A model with a finite future

Let us now consider a model with parameters

$$\mathbf{P_{2}}=\biggl{[}[0.01,0.01,0.1,0.1],0.09876065,0.01\biggr{]}$$ (42)

and initial conditions

$$\mathbf{I_{2}}=[0,0,1,0.01,-1].$$ (43)

4.1 Singular points

Tab. 2. Character (type) of singular points in subspaces $[\Sigma_{\Phi},\Sigma_{\varphi}]$in the model with parameters $\mathbf{P_{2}}$

Singular points	Coordinates	Type
$M_{0,0}^{+}$	$[0,0,0,0,0.0577]$	$[\mathbf{A},\ \mathbf{S}]$
$M_{0,0}^{-}$	$[0,0,0,0,-0.0577]$	$[\mathbf{R},\ \mathbf{S}]$
$M_{0,\pm 1}^{+}$	$[0,0,\pm 1,0,0.0645]$	$[\mathbf{A},\ \mathbf{A}]$
$M_{0,\pm 1}^{-}$	$[0,0,\pm 1,0,-0.0645]$	$[\mathbf{R},\ \mathbf{R}]$
$M_{\pm 1,0}^{+}$	$[\pm 1,0,0,0,2.887]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,0}^{-}$	$[\pm 1,0,0,0,-2.887]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,\pm 1}^{+}$	$[\pm 1,0,\pm 1,0,0.0707]$	$[\mathbf{S},\ \mathbf{A}]$
$M_{\pm 1,\pm 1}^{-}$	$[\pm 1,0,\pm 1,0,-0.0707]$	$[\mathbf{S},\ \mathbf{R}]$

4.2 Evolution of geometric factors

Fig.4.2 shows the evolution of the scale function $\xi(t)=\ln(a(t))$ and the Hubble parameter $H(t)$ (1) for this model.

Fig. 15.  Evolution of $\xi(t)$ (solid line) and $H(t)$ (dashed line) in a model with a finite future at point $t=t_{s}\approx 241$ with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{2}}$ (43).

This model has a rebound point at $t\approx 29$, after which it enters the expansion stage $t\in(30,240)$ with positive inflation with $H_{0}\approx+0.065$, after which it goes into a singular state such as a Big Rip , see, for example, [8]).

Fig. 16.  Graph of $\kappa(t)$ in a model with a finite future at point $t=t_{s}\approx 241$ with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{2}}$ (43).

In Fig.4.2 shows a graph of barotropic coefficients $\kappa(t)$. On this graph, the vertical dotted line corresponds to the singularity $t=t_{s}\approx 241$ - the break point, the horizontal dotted line corresponds to the value $\kappa=1$, the horizontal dash - the dotted line - to the value $\kappa=-1.$

Fig. 17.  Graph of $\Omega(t)$ in a model with a finite future at point $t=t_{s}\approx 241$ with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{2}}$ (43).

In Fig.4.2 shows a graph of the invariant cosmological acceleration $\Omega(t)$. On this graph, the vertical dotted line corresponds to the singularity $t=t_{s}\approx 241$ - the break point, the horizontal dotted line corresponds to the value $\Omega=1$, the horizontal dash - the dotted line - to the value $\Omega=-2$.

Thus, the cosmological model has a singularity in the future corresponding to a large gap, and near the singularity this kind of model behaves in the same way as near the initial singularity: $\kappa\to 1$, $\Omega\to-2$.

4.3 Evolution of scalar fields:
generation of a classical field
phantom

Let us consider the problem of the cosmological evolution of the scalar fields $\Phi(t)$ and $\varphi(t)$. In particular, we are interested in the possibility of generation of one of the components of a scalar doublet by another due to the kinetic connection between them. This question was posed in the first part of the article [1].

Fig.4.3 shows the graphs evolution of the scalar classical $\Phi(t)$ and phantom $\varphi(t)$ potentials, and in Fig.4.3 - the derivatives of these potentials. In the given graphs of these figures one can observe the transition cosmological model from a state of inflationary compression in the infinite past, which, according to Tab.4.1 corresponds to a point $M^{-}_{0,1}$ of type $\mathbf{[R,R]}$ with coordinates $\Phi(-\infty)=0$, $\varphi(-\infty)=1$, into the state of inflationary expansion, which corresponds to a point $M^{+}_{1,-1}$ of type $\mathbf{[S,A]}$ with coordinates $\Phi=1$, $\varphi=-1$. This state ends with a Big Rip.

Fig. 18.  Evolution of scalar classical $\Phi(t)$ (solid) and phantom $\varphi(t)$ (dashed) potentials with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{2}}$ (43).

Fig. 19.  Evolution of dynamic functions $Z(t)$ (solid) and $z(t)$ (dashed) in the model with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{2}}$ (43).

Thus, in this model, at the stage of inflationary expansion before the Big Rip, the generation of a classical scalar field is observed. This process corresponds to a transition with a stable state for the phantom field and an unstable state for the classical one.

5 A model with a finite past and future

Let’s consider a model with parameters $\mathbf{P_{2}}$ (42) and initial conditions very close to the coordinates of the singular point $M^{+}_{1,0}$, according to Tab.4.1 saddle point of type $\mathbf{[S,\ S]}$, namely, just above it (only $10^{9}$ !)

$$\mathbf{I_{3}}=[1+10^{-9},0,0,0,1].$$ (44)

This is a model with a finite past, which corresponds to a singularity at the time $t=t^{0}_{s}\approx-84.206597$, and a finite future, which corresponds to a singularity at the time $t=t^{1}_{s}\approx 280.04421$. Fig.5 shows the evolution of the scale function $\xi(t)=\ln(a(t))$ and the Hubble parameter (1) $H(t)$ for this model. The coordinates of singular points and their characters are indicated in Tab.4.1. Thus, the graphs in Fig.5 demonstrate the transition from an unstable state (saddle point $M_{0,0}^{+}$) through an inflationary plateau to the final singular state (Big Rip).

Fig. 20.  Evolution of the scale function $\xi(t)$ is dashed and the Hubble parameter $H(t)$ is solid line in the model with a finite past $t=t^{0}_{s}\approx-84.206597$ (left vertical dotted line) and a finite future $t=t^{1}_{s}\approx 280.04421$ (right vertical dotted line) with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{3}}$ (44).

In Fig. 5 shows the behavior of the barotropic coefficient $\kappa(t)$ in this model; in the same figure one can see that the cosmological model at both singular points, $t^{0}_{s}$ and $t^{1}_{s}$, manifests itself as matter with an extremely rigid equation of state.

Fig. 21.  Evolution of the barotropic coefficient $\kappa(t)$ in a model with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{3}}$ (44).

In Fig.5–5 the behavior of the barotropic coefficient near singular points is shown on a large scale. Finally, Fig.5 shows graphs of the evolution density of the evolution of the potentials of the classical field and phantom fields for this model. As we noted in [1], at the cosmological singularity points $t^{0}_{s}$ and $t^{1}_{s}$ these potentials are also singular. Outside these points on the inflation interval $t\in(t^{0}_{s},t^{1}_{s})$ the potential values are close to the coordinates of the saddle point $M^{+}_{0,0}$.

Fig. 22.  Evolution of barotropic coefficients $\kappa(t)$ in a model with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{3}}$ (44).

Fig. 23.  Evolution of barotropic coefficients $\kappa(t)$ in a model with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{3}}$ (44).

Note that the lifetime of this cosmological model $\Delta t=t^{1}_{s}-t^{0}_{s}$ is sensitive to the difference between the initial value of the potential of the classical scalar field and its value at the unstable point $M_{+1,0}^{+}$ – $\Delta\Phi_{0}\equiv\Phi_{0}-1$. As $\Delta\Phi_{0}$ decreases, the model lifetime increases according to the empirical law

$$\Delta t\approx 41\cdot|\lg(\Delta\Phi_{0})|.$$

However, the main resource for increasing the lifespan of a model lies in its large-scale transformations. By choosing the similarity coefficient $k=10^{2}$, we will achieve an increase in the model lifetime by exactly $k$ times.

Fig. 24.  Evolution of potentials $\Phi(t)$ is solid and $\varphi(t)$ is dashed in the model with parameters $\mathbf{P_{2}}$ (42) and initial conditions $\mathbf{I_{3}}$ (44).

6 Generating a classical scalar field near the rebound point

As the results of numerical modeling show, a stable classical scalar field corresponding to an inflationary expansion appears to be generated only in models with a finite past at the stage between the rebound point and the Big Rip. In other cosmological models, the classical scalar field is generated only at intermediate stages of cosmological evolution near the rebound point, if there is one. The existence of a rebound point implies that before this point the model was in the compression stage, and after the rebound point it switched to the expansion stage. After this, the system goes into an equilibrium state, the classical scalar field disappears, and further inflation is supported only by the phantom field in stable equilibrium. Let us consider this process in more detail using a specific example of a model with an infinite past and an infinite future:

$$\mathbf{P_{3}}=\bigl{[}[0.01,0.01,0.1,0.1],10^{-2},10^{-6}\bigr{]};$$ (45)

$$\mathbf{I_{4}}=[0,0,0.01,0,1].$$ (46)

Thus, at time $t=0$ there is no classical scalar field, and the phantom one is very small.

6.1 Singular points and geometric factors

The singular points of the model are shown in Tab.6.1.

Tab. 3. Characters (type) of singular pointsin subspaces $[\Sigma_{\Phi},\Sigma_{\varphi}]$ in the model with parameters $\mathbf{P_{3}}$

Singular points	Coordinates	Type
$M_{0,0}^{+}$	$[0,0,0,0,0.000576]$	$[\mathbf{A},\ \mathbf{S}]$
$M_{0,0}^{-}$	$[0,0,0,0,-0.000576]$	$[\mathbf{R},\ \mathbf{S}]$
$M_{0,\pm 1}^{+}$	$[0,0,\pm 1,0,0.0289]$	$[\mathbf{A},\ \mathbf{A}]$
$M_{0,\pm 1}^{-}$	$[0,0,\pm 1,0,-0.0289]$	$[\mathbf{R},\ \mathbf{R}]$
$M_{\pm 1,0}^{+}$	$[\pm 1,0,0,0,2.89]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,0}^{-}$	$[\pm 1,0,0,0,-2.89]$	$[\mathbf{S},\ \mathbf{S}]$
$M_{\pm 1,\pm 1}^{+}$	$[\pm 1,0,\pm 1,0,0.0409]$	$[\mathbf{S},\ \mathbf{A}]$
$M_{\pm 1,\pm 1}^{-}$	$[\pm 1,0,\pm 1,0,-0.0409]$	$[\mathbf{S},\ \mathbf{R}]$

In Fig. 6.1 shows graphs of the evolution of the geometric factors $\xi(t)$ and $H(t)$.

Fig. 25.  Evolution of geometric factors in a model with parameters $\mathbf{P_{3}}$ and initial conditions $\mathbf{I_{4}}$: $\xi(t)$ – dashed and $H(t)$ – solid line.

6.2 Evolution of scalar fields

In Fig.6.2 shows graphs of the cosmological evolution of the components of the scalar doublet in this model.

Fig. 26.  Evolution of scalar fields in a model with parameters $\mathbf{P_{3}}$ and initial conditions $\mathbf{I_{4}}$: $\varphi(t)$ is dashed and $\Phi(t)$ is solid line.

Thus, in the infinite past, the Universe starts from a state close to the singular point $M^{-}_{0,1}=[0,0,1,0,-0.0289]$ and is in a state of inflationary compression. Near the rebound point $t_{b}\approx 0$ oscillations of the phantom and classical scalar fields arise in the system, after which the system goes into a state of inflationary expansion corresponding to the singular point $M^{+}_{0,1}=[0,0,1,0,+0.0289]$, restoring the original values of potentials of scalar fields.

In Fig.6.2–6.2 the process of oscillation of scalar field potentials near the rebound point is shown in close-up.

Fig. 27.  Evolution of scalar fields in a model with parameters $\mathbf{P_{3}}$ and initial conditions $\mathbf{I_{4}}$: $\Phi(t)$ – solid and $\varphi(t)$ – dashed line (close-up).

Fig. 28.  Evolution of derivatives of scalar fields in a model with parameters $\mathbf{P_{3}}$ and initial conditions $\mathbf{I_{4}}$: $Z(t)\equiv\dot{\Phi}$ – solid and $z(t)\equiv\dot{\varphi}$ – dashed line (close-up)

In this model, the interval of oscillation of scalar field potentials, at which the potential of the classical field reaches values of the order of $|\Phi(t_{b})|\sim 0.005$, and the potential of the phantom field becomes zero, takes about 300 Planck times. When using a scaling transformation to real values of fundamental constants for the scale of QCD theory, this period can stretch to $3\cdot 10^{4}$ Planck times, and for the Standard Model - up to $3\cdot 10^{15}t_{Pl}$.

6.3 Production of scalarly charged
fermions near the rebound point

In the above graphs of the evolution of scalar fields (Fig. 2.4, 2.5, 3.2, 3.2, 4.2, 4.3, 6.2 and 6.2) sufficiently large and fast bursts of the potential of the classical scalar field $\Phi$ and its derivative $z(t)$ are observed near the rebound points (or singularity points, if they exist). If we do not take into account exotic cosmological models with a finite future, then the models presented in the above figures reveal in the region of the rebound points the amplitude of oscillations of the derivative of the potential $\max|\dot{\Phi}|\sim 2\cdot 10^{-6}\div 6\cdot 10^{-4}$.

Let us estimate the rate of production, for example, of pairs of scalarly charged fermions with charge $q$ and mass $m_{f}$ in a classical scalar field with a time derivative $Z_{0}$, using an analogy with the formula for the rate of creation of an electron-positron pair in an electric field of strength $E$ (see, for example, [11])

$$\dot{n}_{ee_{+}}=\frac{e^{2}E^{2}}{4\pi^{3}\hbar^{2}c^{2}}\displaystyle\mathrm{e}^{-\frac{E_{0}}{E}},$$

where $w_{ee_{+}}$ is the rate of electron-positron pair production per unit volume, $m$ is the mass of the electron, $e$ is its charge, $E_{c}$ is the critical value of the electric field strength

$$E_{0}=\frac{\pi m^{2}c^{3}}{\hbar e}.$$

Passing to the Planck units in these formulas and making the substitutions $e\to q$, $m\to m_{f}$ and $E\to Z$, we obtain in the adiabatic approximation an estimate of the rate of production of pairs of scalarly charged fermions with scalar charge $q$ and mass $m_{f}$ in scalar field $\Phi(t)$

$$\dot{n}_{qq_{+}}(t)\backsimeq\frac{q^{2}Z^{2}(t)}{4\pi^{3}}\displaystyle\exp\left(-\frac{\pi m^{2}_{f}}{|qZ(t)|}\right).$$ (47)

For efficient production of fermion pairs, the following condition must be met:

$$|Z(t)|\gtrsim\frac{\pi m^{2}_{f}}{|q|},$$ (48)

that is, fermions must be light, but at the same time have a sufficiently large scalar charge. The final density of generated scalarly charged fermions is obtained by integrating over time the expression (47):

$$n_{qq_{+}}=\int\limits_{t_{0}}^{t_{0}+\Delta t}\frac{q^{2}Z^{2}(t)}{4\pi^{3}}\displaystyle\exp\left(-\frac{\pi m^{2}_{f}}{|qZ(t)|}\right)dt$$

$$\backsimeq\frac{q^{2}\overline{Z^{2}(t)}}{4\pi^{3}}\displaystyle\exp\left(-\frac{\pi m^{2}_{f}}{|q\overline{Z(t)}|}\right)\Delta t,$$ (49)

where $\Delta t$ is the duration of the burst of the classical field potential, $\overline{Z(t)}$ is the root-mean-square value of the function $Z(t)$ on the burst interval.

Let us evaluate the possibility of the creation of fermion pairs based on the last model we studied $\mathbf{P_{3}}$ (45), subjecting it, according to Property I.2, to scale transformations in order to bring the model parameters closer to real parameters, for example, SU(5). At the same time, we must also transform the value of the scalar charge according to [7] according to the law $\tilde{q}=q/\sqrt{k}$. Thus, after a scaling transformation with a similarity coefficient $k$, the expression (6.3) is transformed to the form:

$$n_{qq_{+}}(t,k)\backsimeq\frac{q^{2}\overline{Z^{2}(t)}}{4\pi^{3}k^{2}}\displaystyle\exp\left(-\frac{\pi m^{2}_{f}}{\sqrt{k}|q|\overline{|Z(t|)}}\right)\Delta t,$$ (50)

where $\overline{|Z(t)|}$ is the average value of the modulus of the burst of the derivative of the scalar field $\Phi(t)$.

The function $n_{qq_{+}}(t,k)$ has a maximum at

$$k=\left(\frac{\pi m^{2}_{f}}{4|q|\overline{|Z(t)|}}\right)^{2},$$ (51)

from where we get the following:

$$\max(n_{qq_{+}}(t))\backsimeq\frac{q^{4}\overline{Z^{2}(t)}\mathrm{e}^{-4}}{4\pi^{5}m^{4}_{f}}.$$ (52)

At $q=10^{-2}$, $m_{f}=10^{-4}$ and according to Fig.6.2 $\overline{|Z(t)|}\backsimeq 2\cdot 10^{-7}$ we obtain from (52) $\max(n_{qq_{+}}(t))\backsimeq 3\cdot 10^{-10}$. At first glance, this value seems very small, but when assessing it, two important circumstances must be taken into account. Firstly, this value is expressed in Planck units of length; in Compton units of length $\lambda_{f}$ relative to the generated fermions this value will already be $300\lambda^{-3}_{f}$). Secondly, it is necessary to take into account the factor of increase in volume as a result of cosmological expansion $\exp(-3\xi(t))$, which, according to the graph in Fig.6.1 is about $\exp(-3\cdot 3)\thickapprox 10^{-4}$. As a result, upon completion of the process of fermion pair production, we obtain their concentration of the order of $n_{ff_{+}}\gtrsim 3\cdot 10^{-2}\lambda^{3}_{f}$. In addition, the statistical factor $N_{f}$ of the number of fermion types in the interaction model is also important. In the standard SU(5) model (colors, electric charges) $N_{f}\sim 10$, which increases the estimate of the total number of fermion pairs by another order of magnitude.

7 Conclusion

Summing up the main results of both parts of the work, we indicate its following main results.
$\bullet$ A model of the evolution of the Universe based on an asymmetric Higgs scalar doublet with a kinetic coupling proportional to the product of the derivatives of the components of the scalar doublet is proposed and studied. In particular, a qualitative analysis of the dynamic system of the model was carried out, the symmetry properties of the model with respect to the reflection and similarity transformation were investigated and proven.
$\bullet$ The main types of behavior of the cosmological model are revealed depending on the fundamental parameters and initial conditions. Cases of the presence of initial and final cosmological singularities and rebound points have been identified.
$\bullet$ It is shown that near the rebound points the phantom component $\varphi$ of a scalar doublet generates its classical component $\Phi$.
$\bullet$ The asymptotic behavior of the cosmological model near the points of initial singularity and the Big Rip, corresponding to the barotropic coefficient $\kappa=1$ ($p=\varepsilon$), has been proven analytically and confirmed by numerical modeling.
$\bullet$ The probability of fermion pair production in a strong alternating scalar field near the rebound points has been estimated.

All identified features of cosmological models based on an asymmetric scalar doublet: the presence of initial singularities, rebound points and large discontinuity points, generation of scalar doublet components, etc., are similar to the features of cosmological models based on systems of scalarly charged fermions [2] – [5], which apparently makes it possible to replace the mathematically complicated model of the formation of supermassive black holes [9] – [10] with a simpler one, built on a purely field basis. In addition, an important circumstance is also the fairly intense production of fermion pairs by the scalar field near the rebound points. Thanks to this factor, the required number of scalarly charged fermions can be obtained independently of the processes of gravitational pair production, and at the stage of the cold Universe.

Founding

The work is performed according to the Russian Government Program of Competitive Growth of Kazan Federal University