On Bianchi type VI$_0$ spacetimes with orthogonal perfect fluid matter

TL;DR

This paper analyzes the early-time behavior of Bianchi type VI$_0$ spacetimes with perfect fluid matter, confirming a conjecture about their initial singularity being vacuum dominated, anisotropic, and silent, and explores related Klein-Gordon equation asymptotics.

Contribution

It proves Wainwright's conjecture on the nature of initial singularities in Bianchi VI$_0$ spacetimes and extends understanding of their asymptotic properties.

Findings

Initial singularity is vacuum dominated, anisotropic, and silent for generic solutions.

Convergence results for Klein-Gordon solutions on these backgrounds are established.

The conjecture about the initial singularity is confirmed.

Abstract

We study the asymptotic behaviour of Bianchi type VI$_0$ spacetimes with orthogonal perfect fluid matter satisfying Einstein's equations. In particular, we prove a conjecture due to Wainwright about the initial singularity of such spacetimes. Using the expansion-normalized variables of Wainwright-Hsu, we demonstrate that for a generic solution the initial singularity is vacuum dominated, anisotropic and silent. In addition, by employing known results on Bianchi backgrounds, we obtain convergence results on the asymptotics of solutions to the Klein-Gordon equation on all backgrounds of this type, except for one specific case.

Tables

Table 1. Table 1 . The different Bianchi types of class A, depending on the commutator matrix

Type	$n_1$	$n_2$	$n_3$
I	0	0	0
II	+	0	0
VI0	0	+	-
VII0	0	+	+
VIII	-	+	+
IX	+	+	+

Table 2. Table 2 . The isolated equilibria present in the closure of phase space.

Symbol	$(N_{+},N_{-},{\Sigma }_{+},{\Sigma }_{-},\Omega )$	Invariant set	$u$
$F$	$(0,0,0,0,1)$	$B\left(I\right)$	2
$P_2^{+}\left(\mathrm{II}\right)$	$(\frac{1}{4}{p}^{\ast },\frac{1}{4}{p}^{\ast },-\frac{1}{8}{q}^{\ast },-\frac{\sqrt{3}}{8}{q}^{\ast },1-\frac{1}{8}{q}^{\ast })$	$S_2^{+}\left(\mathrm{II}\right)$	1
$P_3^{-}\left(\mathrm{II}\right)$	$(\frac{1}{4}{p}^{\ast },-\frac{1}{4}{p}^{\ast },-\frac{1}{8}{q}^{\ast },\frac{\sqrt{3}}{8}{q}^{\ast },1-\frac{1}{8}{q}^{\ast })$	$S_3^{-}\left(\mathrm{II}\right)$	1
$P_1^{+}\left({\mathrm{VI}}_0\right)$	$(0,\frac{1}{2}{p}^{\ast },-\frac{1}{2}{q}^{\ast },0,\frac{1}{2}\left(2-q^{\ast }\right))$	$S_1^{+}\left({\mathrm{VI}}_0\right)$	0

Equations

$$[e_i,e_j]={\gamma }_{ij}^k{e}_k.$$

$$n^{ij}:=\tfrac{1}{2}\big{(}\epsilon^{jkl}\gamma^{i}_{kl}+\epsilon^{ikl}\gamma^{j}_{kl}\big{)}$$

$$R_{\mu \nu }-\frac{1}{2}S{h}_{\mu \nu }=T_{\mu \nu }.$$

$$h=-dt\otimes dt+a_{ij}(t)x^i\otimes x^j.$$

$$T=(\rho +p)dt\otimes dt+ph.$$

$$\displaystyle\begin{split}N_{1}^{\prime}&=\big{(}q-4\Sigma_{+}\big{)}N_{1},\\[1.0pt] N_{2}^{\prime}&=\big{(}q+2\Sigma_{+}+2\sqrt{3}\Sigma_{-}\big{)}N_{2},\\[1.0pt] N_{3}^{\prime}&=\big{(}q+2\Sigma_{+}-2\sqrt{3}\Sigma_{-}\big{)}N_{3},\\[1.0pt] \Sigma_{+}^{\prime}&=-(2-q)\Sigma_{+}-3S_{+},\\[2.0pt] \Sigma_{-}^{\prime}&=-(2-q)\Sigma_{-}-3S_{-}.\end{split}$$

$${\Omega }^{\prime }=2(q-q^{\ast })\Omega$$

$$q^{\ast }(\gamma ):=\frac{1}{2}(3\gamma -2).$$

$$q:=q^{*}\Omega+2\big{(}\Sigma_{+}^{2}+\Sigma_{-}^{2}\big{)},$$

$$\Omega +{\Sigma }_{+}^2+{\Sigma }_{-}^2+K=1.$$

$$K:=\tfrac{3}{4}\Big{[}N_{1}^{2}+N_{2}^{2}+N_{3}^{2}-2\big{(}N_{1}N_{2}+N_{2}N_{3}+N_{1}N_{3}\big{)}\Big{]}$$

$$\displaystyle\begin{split}S_{+}&:=\tfrac{1}{2}\Big{[}\big{(}N_{2}-N_{3}\big{)}^{2}-N_{1}\big{(}2N_{1}-N_{2}-N_{3}\big{)}\Big{]},\\ S_{-}&:=\tfrac{\sqrt{3}}{2}\big{(}N_{3}-N_{2}\big{)}\big{(}N_{1}-N_{2}-N_{3}\big{)}.\end{split}$$

$$N_{\pm }$$

$$\begin{array}{rl}{N}_{+}^{\prime }& =(q+2{\Sigma }_{+})N_{+}+2\sqrt{3}{\Sigma }_{-}{N}_{-},\\ N_{-}^{\prime }& =(q+2{\Sigma }_{+})N_{-}+2\sqrt{3}{\Sigma }_{-}{N}_{+},\\ {\Sigma }_{+}^{\prime }& =-(2-q){\Sigma }_{+}-2N_{-}^2,\\ {\Sigma }_{-}^{\prime }& =-(2-q){\Sigma }_{-}-2\sqrt{3}{N}_{+}{N}_{-},\\ {\Omega }^{\prime }& =2(q-q^{\ast })\Omega .\end{array}$$

$$\Omega +{\Sigma }_{+}^2+{\Sigma }_{-}^2+N_{-}^2=1.$$

$$B_1^{+}({\mathrm{VI}}_0):=\left\{(N_{+},N_{-},{\Sigma }_{+},{\Sigma }_{-},\Omega )\in {\mathbb{R}}^5:\begin{array}{rl}& \Omega +{\Sigma }_{+}^2+{\Sigma }_{-}^2+N_{-}^2=1,\\ & N_{-}>\text{absolutevalue}{N}_{+},\Omega \ge 0\end{array}\right\}.$$

$$N_{+}(\tau):=\big{(}N_{+}\circ\varphi^{\tau}\big{)}(x):=\big{(}\varphi^{\tau}(x)\big{)}_{1}$$

$$\displaystyle\begin{split}\alpha(x)&\subseteq\big{\{}s\in\overline{S}\setminus S~{}|~{}\lim_{y\to s}Z(y)\neq a\big{\}},\\[2.0pt] \omega(x)&\subseteq\big{\{}s\in\overline{S}\setminus S~{}|~{}\lim_{y\to s}Z(y)\neq b\big{\}}.\end{split}$$

$$u^{\prime }(\tau )\ge \beta (\tau )u(\tau )$$

$$u(\tau )\le u(b)\exp \left(-{\int }_{\tau }^b\beta (s)ds\right).$$

$$u^{\prime }(\tau )\le -\beta (\tau )u(\tau )$$

$$u(\tau )\ge u(b)\exp \left({\int }_{\tau }^b\beta (s)ds\right).$$

$$\partial B_{1}^{+}(\mathrm{VI}_{0})=B(\mathrm{I})\cup B^{+}_{2}(\mathrm{II})\cup B_{3}^{-}(\mathrm{II})\cup B_{1}^{+}(\mathrm{VI}_{0})\big{|}_{\Omega=0}.$$

$$B(\mathrm{I}):=\overline{B_1^{+}({\mathrm{VI}}_0)}\cap \{N_{+}=0,N_{-}=0\}.$$

$$\displaystyle\begin{split}\Sigma_{+}^{\prime}&=-(2-q)\Sigma_{+},\\[2.0pt] \Sigma_{-}^{\prime}&=-(2-q)\Sigma_{-}\\[2.0pt] \Omega^{\prime}&=2\big{(}q-q^{*}\big{)}\Omega.\end{split}$$

$$\underset{\tau \to -\infty }{\lim }({\Sigma }_{+},{\Sigma }_{-},\Omega )(\tau )=\frac{1}{\sqrt{{\sigma }_{+}^2+{\sigma }_{-}^2}}({\sigma }_{+},{\sigma }_{-},0)$$

$$\big{(}\Sigma_{+},\Sigma_{-}\big{)}=\big{(}-1,0\big{)},\big{(}\tfrac{1}{2},\tfrac{1}{2}\sqrt{3}\big{)},\big{(}\tfrac{1}{2},-\tfrac{1}{2}\sqrt{3}\big{)}.$$

$$\begin{array}{rl}{\beta }_1& :=q-4{\Sigma }_{+},\\ {\beta }_2& :=q+2{\Sigma }_{+}+2\sqrt{3}{\Sigma }_{-},\\ {\beta }_3& :=q+2{\Sigma }_{+}-2\sqrt{3}{\Sigma }_{-},\end{array}$$

$$q(x)-q^{\ast }>\epsilon \text{for}x\in B_{\epsilon }(T_1)$$

$${\Omega }^{\prime }=2(q-q^{\ast })\Omega \ge 2\epsilon \Omega .$$

Full text

1. Introduction

The subject of this study is the asymptotic behaviour of a certain class of spatially homogeneous cosmological models. Our interest is twofold. First, we prove a conjecture, made over twenty years ago in [11], regarding the initial singularity of Bianchi type $\mathrm{VI}_{0}$ spacetimes with orthogonal perfect fluid (OPF) matter which satisfy Einstein’s equations. Resolving this conjecture allows for the study of more complicated cosmological models, for which those of the type above appear as limit cases. Second, by combining the asymptotic behaviour that we find with known results on Bianchi models from [8], we are able to obtain results regarding the asymptotics of solutions to the Klein-Gordon equation on this type of model. This fills a gap in the unified treatment of the Klein-Gordon equation on Bianchi backgrounds of [8].

In order to reduce the complexity of Einstein’s equations one often demands a high degree of symmetry to be present in the spacetime as well as in the matter model. Typically, this comes in the form of isotropy or spatial homogeneity. As a prime example, the Friedmann-Lemaître-Robertson-Walker (FLRW) spacetimes are both isotropic and spatially homogeneous. Dropping the requirement of isotropy, but retaining that of spatial homogeneity, one finds the Kantowski-Sachs spacetimes and the Bianchi spacetimes. For a Bianchi spacetime, there is a three-dimensional group of isometries acting transitively and freely on spacelike hypersurfaces, while for the Kantowski-Sachs spacetimes the situation is similar but the action is not free.

The underlying Lie group of a Bianchi spacetime can be classified by the structure constants of the corresponding Lie algebras, by the so-called Bianchi-Schücking-Behr approach, as described in e.g. Section 2 of [1]. According to this classification, one distinguishes the class A types I, II, VI0, VII0, VIII and IX, and the class B types III, IV, V, VIh, VIIh, for a parameter $h\neq 0$. As mentioned, the subject of this study are solutions of Bianchi type VI0 with perfect fluid matter with velocity vector orthogonal to the group orbits. For solutions of type I, II and VII0 a similar treatment – on which this one is largely based – is already available, see Sections 8-10 of [9].

Using the orthonormal frame formalism due to Ellis and MacCallum in [1], one may write the Einstein field equations for Bianchi spacetimes as evolution equations for certain cosmological variables. This works for a variety of matter models, cf. Section 1.1.2 of [11]. Then, by normalizing by the rate of expansion of the universe, these evolution equations transform into polynomial differential equations in the mean curvature, the shear, the structure constants of the group and the energy density. This approach was introduced by Wainwright and Hsu in [12] for all class A types, for vacuum as well as for OPF matter. For the class B types a similar approach exists; see [7] for a recent treatment. Their approach then allows one to study the asymptotic behaviour of these quantities by employing methods from the theory of dynamical systems. Here we are mainly interested in the approach to the initial singularity, which is when length scales become arbitrarily small.

The main result of this article is the proof of the aforementioned conjecture. The conjecture can be found in Section 6.3.3 of [11], and is stated as Theorem 1.6 below. Translating the statement to more physical terms: generically, the initial singularity of a Bianchi type $\mathrm{VI}_{0}$ orthogonal perfect fluid solution is vacuum dominated, anisotropic and silent. The genericity here is within the set of all solutions of type $\mathrm{VI}_{0}$ with orthogonal perfect fluid matter.

The asymptotics described below can be used to complete studies of type $\mathrm{VI}_{0}$ solutions with other matter models, for example the case of a non-orthogonal fluid or that of the magnetic type $\mathrm{VI}_{0}$ solutions with OPF matter, see [4] and [6] respectively. In such compatible treatments (i.e. employing an approach similar to that of Wainwright and Hsu) we find the type $\mathrm{VI}_{0}$ solutions with OPF matter as limit cases. For the much harder case of Bianchi type VIII, the case of Bianchi type $\mathrm{VI}_{0}$ appears as a limit case as well. Although studies of the late-time dynamics exist for Bianchi type VIII, see [5], the precise nature of the initial singularity remains a mystery for the non-vacuum case.

In addition to proving the conjecture of Wainwright, we also apply the derived asymptotics to fill a gap in [8]. In that article Ringström demonstrates the convergence of time derivatives (in certain geometrically defined time coordinates) of solutions to the Klein-Gordon equation on Bianchi backgrounds. His methods apply to Bianchi $\mathrm{VI}_{0}$ solutions as well, however, due to the lack of proven asymptotics, this case was left out there. In general, the conclusions made here are very similar to those made for Bianchi types I, II and VII0 solutions with OPF matter; the precise statement is in Theorem 7.1 below. Let $u$ denote a solution to the Klein-Gordon equation on a generic Bianchi type $\mathrm{VI}_{0}$ orthogonal perfect fluid development. Then, towards the initial singularity, the time derivative $u_{\tau}(\cdot,\tau)$ – with respect to the expansion-normalized time – converges exponentially to a smooth function on the underlying group.

1.1. Outline

In the rest of this section we recall some preliminaries: the Bianchi types of class A, the notion of a Bianchi spacetime and the stress-energy tensor of an orthogonal perfect fluid. We also recall the Wainwright-Hsu equations, note some properties specific to the case of Bianchi type $\mathrm{VI}_{0}$, and state Theorem 1.6. We finish the section with recalling two important tools from the theory of dynamical systems.

In Section 2 we start out by describing the invariant sets in the boundary of phase space and the dynamics therein. There we find invariant sets of types I and II as well as the vacuum case. We also consider the dynamics within the shear invariant set, in Section 3, which contains the global attractor $P_{1}^{+}(\mathrm{VI}_{0})$. In Section 4 we briefly show the regularity of the unstable sets of the equilibria from Table 2. The main ingredient of the proof of Theorem 1.6 is found in Section 5, where we deal with the generic case $(c)$ of convergence to the Kasner circle.

In the second half of this paper we apply Theorem 1.6 to the Klein-Gordon equation on Bianchi type $\mathrm{VI}_{0}$ backgrounds, based on results from [8]. The relevant theorem is formulated in Section 7, after having introduced the necessary concepts to state it in Section 6. We consider mainly the generic case $(c)$ of Theorem 1.6, but also briefly touch on the other two cases.

1.2. The Bianchi types of class A

A connected, three-dimensional Lie group $G$ can be classified by its Lie algebra $\mathfrak{g}$, based on properties of the structure constants of $\mathfrak{g}$.

We say $G$ is of class A if the Lie algebra $\mathfrak{g}$ is unimodular, which means that any $X\in\mathfrak{g}$ satisfies $\tr(\operatorname{ad}_{X})=0$. Any associated element (Bianchi spacetime, initial data etc.) we also dub class A if the Lie group is of class A. If the three-dimensional Lie group is not of class A, it is of class B.

Given a basis $e_{i},i=1,2,3,$ of $(\mathfrak{g},[\cdot,\cdot])$, recall that the structure constants $\gamma^{k}_{ij}$ are defined by

[TABLE]

The class A algebras are precisely those for which $\gamma^{k}_{ik}=0$. We can then define the quantities

[TABLE]

which form a symmetric matrix $n$.

It can be shown that by choosing an appropriate basis of $\mathfrak{g}$ this matrix can be diagonalized, with diagonal elements $n_{i},i=1,2,3$. The Bianchi class A Lie algebras can now be classified by the signs appearing on the diagonal of the commutator matrix, c.f. Lemma 19.8 of [10] and the accompanying Table 1, which we copied for convenience. This leads to the aforementioned types, of which we focus on type $\mathrm{VI}_{0}$. In case the Lie group has an algebra of type $\mathrm{VI}_{0}$, any associated element (Bianchi spacetime, initial data etc.) is also designated type $\mathrm{VI}_{0}$.

In the case of Bianchi type $\mathrm{VI}_{0}$, we may read off from Table 1 that one diagonal element vanishes, while the other two have opposite signs. The signs in Table 1 are, up to permutation, precisely the signs of the $N_{i},i=1,2,3$, of Section 1.4 below.

1.3. Bianchi type VI0 orthogonal perfect fluid solutions

The spacetimes under consideration are assumed to satisfy Einstein’s equations without cosmological constant. Thus, given our Lorentz manifold $(M,h)$, we have

[TABLE]

Here $R_{\mu\nu}$ and $S$ denote the Ricci curvature tensor and the scalar curvature of $(M,h)$ respectively. For the energy-stress tensor $T$ we take that of an orthogonal perfect fluid with linear equation of state, which is specified below. We demand our Lorentz manifold to be spatially homogeneous or, more specifically, that it is a Bianchi spacetime (cf. Definition 1 of [8]).

Definition 1.1**.**

A Bianchi spacetime is a Lorentz manifold $(M,h)$ of the form $M=G\times I$, where $G$ is a connected, three-dimensional Lie group and $I=(t_{-},t_{+})$ an open interval; the metric is of the form

[TABLE]

Here the $x^{j}$ are the dual basis of the basis $e_{i}$ of the Lie algebra $\mathfrak{g}$ of $G$. Moreover, we require that the functions $a_{ij}(t)$ are smooth and form a positive definite matrix $a(t)$ for every $t\in I$.

Remark 1.2**.**

It is important to note that there are no locally rotationally symmetric (LRS) Bianchi $\mathrm{VI}_{0}$ solutions. The LRS Bianchi solutions form a special subclass of the Bianchi solutions possessing a small degree of isotropy, see Section 1.2.2 of [11]. A necessary condition for a Bianchi spacetime to be LRS is that at least two $n_{i}$, $i=1,2,3$, in the table above are equal, which cannot happen for type VI0 due to the occuring signs. However, for Bianchi VI0 solutions there is the special subclass $n_{j}^{j}=0$, which forms an important object in our analysis; we encounter it as the shear invariant set in Section 3. ∎

Now, with the metric in this form, we may write the stress-energy tensor for orthogonal perfect fluid matter as

[TABLE]

Here the variable $\rho=\rho(t)$ is known as the energy density and is related to the pressure $p$ by $p=(\gamma-1)\rho$; in particular, we demand a linear equation of state. The parameter $\gamma$ specifies the type of fluid under consideration, e.g. $\gamma=1$ corresponds to dust and $\gamma=4/3$ to a radiation fluid.

Remark 1.3**.**

In this article we restrict ourselves to the range $\gamma\in(\tfrac{2}{3},2)$; in these cases both the strong and dominant energy conditions are satisfied as long as $\rho>0$. This is also true for $\gamma=2$, known as the stiff fluid case, but for Bianchi VI0 this case has been well studied already, see Section 7 of [9] for results in the dynamical systems approach, and Example 22 of [8] for results regarding the Klein-Gordon equation. The vacuum case, when $\rho=0$, has been studied extensively as well, see [3] and Section 23.2 of [10] for results in the dynamical systems approach, and Example 24 of [8] for results regarding the Klein-Gordon equation. We consider the vacuum case briefly in Section 2, as it appears naturally in the analysis. ∎

1.4. The Wainwright-Hsu variables

In what follows, we adhere to the conventions of Ringström’s work on the Bianchi IX attractor, see [9], although the original formulation of the evolution equations in this manner is due to Wainwright and Hsu, see [12] or Section 6.1 of [11]. The original formulation is only slightly different, and still used in e.g. [5].

We use expansion-normalized variables $N_{j},j=1,2,3$, $\Sigma_{\pm}$ and $\Omega$, as defined in Section 21 of [9]. We may then restrict ourselves to an invariant set of the right type by giving the $N_{j}$ the same sign as the $n_{j}$ in Table 1 above, up to permutation, as the $N_{j}$ are precisely the normalized version of the $n_{j},j=1,2,3$. This yields a hierarchy of invariant sets, where the invariant sets of lower Bianchi types (I, II, VI0 and VII0) lie in the boundary of the higher types, cf. Figure 6.1 of [11]. The evolution equations take the following form:

[TABLE]

Since we take matter into account this system of equations is completed by the evolution equation

[TABLE]

for the normalized matter density. We take $\Omega\geq 0$, due to physical requirements; $\Omega$ is related to the energy density $\rho$ through a normalization. The set $\{\Omega=0\}$ is invariant with respect to the dynamics; this set is known as the vacuum case and its orbits as vacuum solutions.

The constant $q^{*}$ is determined by the parameter $\gamma\in(\tfrac{2}{3},2)$ we encountered above. We have

[TABLE]

In the equations of motion above, the function $q$ is given by

[TABLE]

and it is known as the decelaration paramater.

The variables are not independent but related through the (Hamiltonian) constraint

[TABLE]

This equation is conserved by the equations of motion (1.6). In particular, the variables live on a hypersurface in $\mathbb{R}^{6}$, which is preserved under the equations of motion. The phase space is thus five-dimensional. In the constraint, the function $K$ is short for

[TABLE]

and lastly we have

[TABLE]

1.5. The Bianchi type VI0 phase space

Following Table 1, for Bianchi type VI0 we look at the invariant set defined by setting $N_{1}~{}=~{}0,~{}N_{2}~{}>~{}0$ and $N_{3}~{}<~{}0$. The fact that this set is invariant is immediate from the system (1.6) due to the homogeneity of the first three equations. A permutation of the indices of the $N_{i}$ can be handled by applying certain symmetries, see Section 2 of [9], meaning our choice is without loss of generality.

If we define

[TABLE]

then we find that the equations of motion reduce to

[TABLE]

The constraint can be written as

[TABLE]

We define the phase space for this system, now a hypersurface in $\mathbb{R}^{5}$, to be the set

[TABLE]

Notation**.**

Given $x\in\overline{B_{1}^{+}(\mathrm{VI}_{0})}$ we write $(N_{+}(x),...,\Omega(x))$ for $(x_{1},...,x_{5})$. Moreover, given initial data $x\in\overline{B_{1}^{+}(\mathrm{VI}_{0})}$ we may write for example

[TABLE]

and similarly for the other coordinates. Here and throughout this article $\varphi^{\tau}$ denotes the flow of the dynamical system.

Remark 1.4**.**

The phase space is bounded by the constraint, and so the closure of the phase space – which forms an invariant set – is compact. The bounds on the variables moreover grant us growth bounds since the vector field is polynomial. This ensures that, for any point $x$ in the closure of the phase space, the flow $\varphi^{\tau}(x)$ is complete.

Also, recall that $\alpha$- and $\omega$-limit sets are always closed and invariant. Since the orbits of our dynamical system are contained in a compact invariant set, any limit set of a point in $B_{1}^{+}(\mathrm{VI}_{0})$ is non-empty and connected, see e.g. Proposition 1.1.14 of [13]. ∎

The closure of the phase space contains several isolated equilibria, of which all but one reside in the boundary. In Sections 2 and 3 we discuss these and their invariant manifolds in more detail. In Table 2 we list the equilibria; this table can be found (partly) in Section 6.2.1 of [11]. Also, see Figure 6.3 of [11] for an overview of how the equilibria lie in the $(\Sigma_{+},\Sigma_{-})$-plane. Here $p^{*}=p^{*}(\gamma)$ stands for $p^{*}(\gamma):=\sqrt{q^{*}(2-q^{*})}\in(0,1)$. The column with $u$ denotes the dimension of the unstable manifold relative to the phase space; note that all these equilibria are hyperbolic for our range of $\gamma$ and moreover, by Proposition 4.1, the unstable sets are indeed unstable manifolds. The definitions of the invariant sets denoted can be found in Sections 2 and 3.

The closure of the phase space also contains a circle of equilibria, called the Kasner circle. The circle is the unit circle in the $(\Sigma_{+},\Sigma_{-})$-plane and the corresponding solutions are vacuum solutions of type I. These are discussed in greater detail in Section 2. For the statement of the theorem let us note that the Kasner circle contains three special points; after removing these we are left with three components called $\mathcal{K}_{j},j=1,2,3$.

Remark 1.5**.**

For the computations of the linearizations around equilibria in the closure of the phase space one rewrites $\Omega$ using (1.15) and subsequently ignores the direction $\partial_{\Omega}$. In particular one computes the linearizations in $\mathbb{R}^{4}$; in essence $\partial_{\Omega}$ only measures the behaviour normal to the hypersurface defined by (1.15), which is not relevant for the desired stability results. ∎

Now let us state our theorem.

Theorem 1.6**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})$, with $\Omega(x)>0$ and $\gamma\in(\tfrac{2}{3},2)$. Then, either

1. (a)

$x$* is the equilibrium $P_{1}^{+}(\mathrm{VI}_{0})$,* 2. (b)

$x$* lies in the unstable manifold of one of the equilibria $F,P_{2}^{+}(\mathrm{II})$ or $P_{3}^{-}(\mathrm{II})$, or* 3. (c)

$\varphi^{\tau}(x)$* converges to a point on $\mathcal{K}_{1}$ as $\tau\to-\infty$.*

Let us note that the unstable manifolds of case $(b)$ are of dimension not greater than two, by Table 2, so, in particular, case $(c)$ is generic, cf. Definition 4.2 below.

1.6. Tools from dynamical systems theory

An important tool throughout the analysis of Sections 2 to 5 is the monotonicity principle. This is due to the presence of monotone functions on large parts of the phase space of Bianchi VI0. The statement and proof of the proposition below can be found in appendix A of [6].

Proposition 1.7** (Monotonicity principle).**

Let $\varphi^{\tau}$ be a flow on $\mathbb{R}^{n}$ and $S$ an invariant set of $\varphi^{\tau}$. Let $Z:S\to\mathbb{R}$ be a $C^{1}$-function whose range is the interval $(a,b)$, where $a\in\mathbb{R}\cup\{-\infty\},b\in\mathbb{R}\cup\{\infty\}$ and $a<b$. If $Z$ is strictly monotonically decreasing on orbits in $S$, then for all $x\in S$ we have

[TABLE]

Another tool we frequently use is Grönwall’s lemma, see e.g. Lemma 7.1 of [10]. We only state a fairly elementary version, which suffices for our purposes. Note the reversed orientation of time.

Lemma 1.8** (Grönwall’s lemma).**

Let $I:=(a,b)$ be in an interval with $-\infty\leq a\leq b<\infty$. Let $u\in C(\bar{I})\cap C^{1}(I)$ and $\beta\in C(\bar{I})$ be non-negative functions.

1. (i)

If $u$ satisfies the estimate

[TABLE]

for all $\tau\in I$, then we have

[TABLE] 2. (ii)

If $u$ satisfies the estimate

[TABLE]

for all $\tau\in I$, then we have

[TABLE]

2. The boundary of phase space

In this section we describe the boundary of the Bianchi type VI0 phase space. Together with the shear invariant set of the next section, this forms the so-called ‘skeleton’ of the phase space, and plays an important role in the asymptotics, mostly due to its invariance and the presence of monotone functions.

In order to prove Theorem 1.6, we need information about the boundary $\partial B_{1}^{+}(\mathrm{VI}_{0})$ of our phase space. Many of the proofs below use methods similar to those in [9], while the description of invariant sets can largely be found in Chapter 6 of [11]. Some of the lemmata below are simply restated for convenience, their proofs are largely omitted and can be found in the cited sources.

We divide the boundary into invariant sets as follows:

[TABLE]

The first three sets on the left hand side, the cases of Bianchi type I and Bianchi type II, are due to setting $N_{-}^{2}=N_{+}^{2}$ (so $N_{2}N_{3}=0$); the last set appearing in the union above is the vacuum case.

2.1. Bianchi type I

The invariant set $B(\mathrm{I})$ consists of those points in the extended phase space satisfying $N_{1}=N_{2}=N_{3}=0$, or in our setting,

[TABLE]

The dynamics simplify considerably. We have

[TABLE]

If we set $\Omega=1$ and $\Sigma_{+}=\Sigma_{-}=0$, then the system is in equilibrium, as in that case $q=q^{*}$ by (1.15). This equilibrium at $(0,0,0,0,1)$ is known as $F$.

We will see below that the rest of $B(\mathrm{I})$ forms the stable manifold of $F$. In $B_{1}^{+}(\mathrm{VI}_{0})$ there is also an unstable set (of points converging to $F$ as $\tau\to-\infty$) which we denote by $\mathcal{F}_{\mathrm{VI}_{0}}$. This set is in fact an embedded smooth manifold of dimension 2, as we show in Proposition 4.1.

The Kasner circle

Any point on the circle $\{\Sigma_{+}^{2}+\Sigma_{-}^{2}=1,\Omega=0\}$ is also an equilibrium, as in that case $q=2$. This circle of equilibria is known as the Kasner circle, denoted by $\mathcal{K}$. Note that each of these equilibria corresponds to a vacuum solution, and the orbit corresponding to a vacuum solution of type I is always contained in $\mathcal{K}$. The corresponding developments (see Section 6.3 and Section 9.1.1 of [11]) are the Kasner vacuum solutions, discovered by Kasner back in 1925.

From the system of equations (2.3) we deduce that the remaining orbits project to radial lines in the $\Sigma_{+}\Sigma_{-}$-plane; the solutions go from the Kasner circle towards the equilibrium $F$ as time increases. This proves the lemma below, which is Proposition 8.1 of [9].

Lemma 2.1**.**

Let $x\in B(\mathrm{I})\setminus\{F\}$ for $\gamma\in(\tfrac{2}{3},2)$. Then

[TABLE]

where $(\sigma_{+},\sigma_{-})=(\Sigma_{+},\Sigma_{-})(x)$. In particular, $\alpha(x)$ consists of a single point on the Kasner circle.

The lemma above yields another necessary lemma, as a simple corollary of Lemma 4.2 of [9]. As noted there, the assumptions of the lemma below are never satisfied, but we need the implication to rule out the case.

Lemma 2.2**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})$ for $\gamma\in(\tfrac{2}{3},2)$ and assume that $F\in\alpha(x)$, but the solution $\varphi^{\tau}(x)$ does not converge to $F$ as $\tau\to-\infty$. Then $\alpha(x)$ contains a point on the Kasner circle.

There are a few points on $\mathcal{K}$ of special interest, the so-called Taub points. Dynamically they are of interest due to their eigenvalue analysis, cf. Section 6.2.2. of [11]. Unlike the other points of the Kasner circle, which only have one eigenvalue 0, the Taub points have two vanishing eigenvalues. They are also of interest due to the geometry of the corresponding developments, see Section 6.3, for which the spacetimes are flat and non-silent.

Definition 2.3**.**

The Taub points $T_{1},T_{2},T_{3}\in\mathcal{K}$ are those points on the Kasner circle with the respective coordinates

[TABLE]

The Taub points separate the arcs $\mathcal{K}_{j}\subset\mathcal{K}$, where $\mathcal{K}_{1}$ is the arc bounded by $T_{2}$ and $T_{3}$ et cetera, see Figure 1 below for an illustration of $\mathcal{K}_{1}$ and the Taub points in the $\Sigma_{+}\Sigma_{-}$-plane. The points anti-podal on $\mathcal{K}$ to the Taub points $T_{j}$ we call $Q_{j},j=1,2,3$. (So, in particular, $Q_{j}\in\mathcal{K}_{j},j=1,2,3$).

Remark 2.4**.**

The arcs $\mathcal{K}_{j}$ can be characterized alternatively as $\mathcal{K}_{j}:=\mathcal{K}\cap\{\beta_{j}<0\}$, where the $\beta_{j}$ are

[TABLE]

see Definition 6.1 of [9]. Note that $q|_{\mathcal{K}}\equiv 2$, by the constraint.

In particular, we have exponential decay of $N_{j}$ near the corresponding arc $\mathcal{K}_{j}$ for $j=1,2,3$, since the $\beta_{j}$ are exactly $N_{j}^{\prime}/N_{j}$. The quantities $\beta_{j}$ are strictly positive on the arcs $\mathcal{K}_{i}$ if $i\neq j$ and we also have $\beta_{j}(T_{j})>0,j=1,2,3$. However, the other Taub points are exactly where the $\beta_{j}$ vanish, i.e. $\beta_{j}(T_{i})=0$ if $i\neq j$. ∎

Remark 2.5**.**

In the introduction we claimed that, generically, the initial singularity of a Bianchi type VI0 spacetime with OPF matter is vacuum dominated, anisotropic and silent. The generic case $(c)$ of Theorem 1.6 is convergence to $\mathcal{K}_{1}$ and any non-special point on the Kasner circle indeed satisfies the properties mentioned. The Taub points are the only points on $\mathcal{K}$ that are non-silent.

Regarding anisotropy we should note that both the Taub points $T_{j}$ as well as their anti-podal points $Q_{j},j=1,2,3$, are locally rotationally symmetric and thus contain a small degree of isotropy. Since $\mathcal{K}_{1}$ contains the point $Q_{1}$, this can occur as a limit also in the generic case $(c)$ of Theorem 1.6. However, we suspect that the only points in $B_{1}^{+}(\mathrm{VI}_{0})$ that converge to $Q_{1}$ towards the past are contained in the shear invariant set $S_{1}^{+}(\mathrm{VI}_{0})$, which is defined below. ∎

We note that convergence (as time goes towards the past) to the Taub points does not happen in Bianchi VI0, as an immediate consequence of Proposition 3.1 of [9]. We only give the necessary statement for Bianchi VI0, which has no locally rotationally symmetric invariant set, see Remark 1.2 above.

Corollary 2.6**.**

For $\gamma\in(\tfrac{2}{3},2)$ there are no $x\in B_{1}^{+}(\mathrm{VI}_{0})$ such that the solution $\varphi^{\tau}(x)$ converges to a Taub point $T_{i},i=1,2,3$, as $\tau\to-\infty$.

Similar to Lemma 2.2, we have a lemma concerning the Taub point $T_{1}$. The proof is an adaptation from the proof of Lemma 4.1 of [9]

Lemma 2.7**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})$ for $\gamma\in(\tfrac{2}{3},2)$ and assume that $T_{1}\in\alpha(x)$. Then there also exists a different $\alpha$-limit point $y\in\alpha(x)\setminus\{T_{1}\}$ such that $\Omega(y)=0$.

Proof.

By Corollary 2.6 we know that the solution cannot converge to the Taub point $T_{1}$. Therefore there exists a $\delta>0$ such that for any time $\tau$ there is a time $s\leq\tau$ such that $\varphi^{s}(x)\notin B_{\delta}(T_{1})$. We may choose $\delta$ small enough such that for some $\varepsilon>0$ we have

[TABLE]

as $q(T_{1})=2>q^{*}$.

Let $(\tau_{k})_{k\in\mathbb{N}}$ be a sequence such that $\varphi^{\tau_{k}}(x)\to T_{1}$. We may assume that $\varphi^{\tau_{k}}(x)\in B_{\delta}(T_{1})$ for any $k\in\mathbb{N}$. Let $s_{k}$ be the first time prior to $\tau_{k}$ such that the solution lies on the boundary of the ball, i.e. $\varphi^{s_{k}}(x)\in\partial B_{\delta}(T_{1})$ and for $r\in(s_{k},\tau_{k}]$ we have $\varphi^{r}(x)\in B_{\delta}(T_{1})$.

Then by compactness, we may select a convergent subsequence $(\varphi^{s_{k_{j}}}(x))_{j}$, yielding an $\alpha$-limit point $y$ on the boundary $\partial B_{\delta}(T_{1})$. In the ball $B_{\delta}(T_{1})$ we have

[TABLE]

Towards the past we see that the value $\Omega$ is exponentially decaying inside the ball, by applying Lemma 1.8(i). Thus we compute that,

[TABLE]

Thus there exists $y\in\alpha(x)$ such that $d(T_{1},y)=\delta>0$ and $\Omega(y)=0$. ∎

2.2. Bianchi type II vacuum dynamics

An invariant set of Bianchi type II consists of points satisfying $N_{j}\neq 0,$ (as well as the constraint etc.) where $j$ is exactly one of $1,2,3$, and the other two $N_{i}$ vanish, $i\neq j$. In the boundary of the phase space we find two of these invariant sets. Here we only consider the case of $B_{2}^{+}(\mathrm{II})$, i.e. the case that $N_{2}>0$ and $N_{1}=N_{3}=0$. In our setting this translates to

[TABLE]

The analysis of $B_{3}^{-}(\mathrm{II})$ is analogous (but the equality $N_{-}=-N_{+}$ is satisfied instead) so we omit it. The dynamics in the set $B_{2}^{+}(\mathrm{II})$ are given by

[TABLE]

where we denote the common value of $N_{+}$ and $N_{-}$ by $N_{\pm}$ (which is a multiple of $N_{2}$).

Type II vacuum orbits are important in order to understand of the dynamics for the higher Bianchi types, so let us elaborate a bit on these. By the constraint, we can rewrite the equations of motion for the vacuum case to

[TABLE]

In particular, $1+\Sigma_{+}$ is a constant multiple of $\sqrt{3}+\Sigma_{-}$, so the quotient

[TABLE]

is conserved along the vacuum orbit. The intersection of the plane $\{f(x)=c\}$, for some constant $c$, with the sphere $\{\Sigma_{-}^{2}+\Sigma_{+}^{2}+N_{\pm}^{2}=1\}\subset B_{2}^{+}(\mathrm{II})$ gives the vacuum orbit through $x$ and its boundary. For a vacuum point $x\in B_{2}^{+}(\mathrm{II})$, the $\alpha$-limit set is but a point on $\mathcal{K}_{2}$, while $\omega(x)$ is a point on the open arc $\mathcal{K}_{1}\cup\{T_{2}\}\cup\mathcal{K}_{3}$, see Proposition 5.1 of [9].

By considering the vacuum orbits of the other type II systems (with $N_{1}\neq 0$ or $N_{3}\neq 0$), which show similar behaviour, we obtain a map $\mathcal{K}\to\mathcal{K}$, cf. Definition 5.1 of [9].

Definition 2.8**.**

We define the Kasner map $K:\mathcal{K}\to\mathcal{K}$ as follows: for $x\in\mathcal{K}$, let $K(x)$ be the point we obtain as the $\alpha$-limit of the vacuum type II orbit for which $x$ is the $\omega$-limit.

Note that the Taub points $T_{j},j=1,2,3$, are fixed points of $K$, while their other pre-images are precisely the $Q_{j}$ for the same $j$.

A lemma on the Kasner map

An important lemma in the proof of the theorem is the following variation of Proposition 3.3 of [9], for the case of Bianchi type VI0. Although we prove later that $B_{1}^{+}(\mathrm{VI}_{0})$ has no $\alpha$-limit points on $\mathcal{K}_{2}$ or $\mathcal{K}_{3}$, we need the implication for the proof.

Lemma 2.9**.**

Assume $x\in B_{1}^{+}(\mathrm{VI}_{0})$ has an $\alpha$-limit point $y\in\mathcal{K}_{2}\cup\mathcal{K}_{3}$. Then the image $K(y)$ under the Kasner map is also contained in $\alpha(x)$.

Proof.

Assume without loss of generality that the limit point $y\in\alpha(x)$ lies in $\mathcal{K}_{2}$. Due to Remark 2.4, this implies that $\beta_{2}(y)<0$. However, $\beta_{3}(y)>0$, and also

[TABLE]

Note that the only point on $\mathcal{K}$ where this quantity vanishes – which is precisely $\frac{d}{d\tau}\log\left(N_{-}^{2}-N_{+}^{2}\right)$ – is the Taub point $T_{1}$, as we have $\Sigma_{+}=-1$ only there.

Therefore there are $\delta>0$ and $\varepsilon>0$ such that in the ball $B_{\delta}(y)$ we have

[TABLE]

while

[TABLE]

So inside the ball $B_{\delta}(y)$ the value of $N_{-}+N_{+}$ is increasing exponentially towards the past, while $\Omega,N_{-}^{2}-N_{+}^{2}$ and thus also $N_{-}-N_{+}$ are decreasing exponentially; this follows from Lemma 1.8 (i) and (ii).

The fact that $N_{-}+N_{+}$ is increasing exponentially implies that the solution cannot remain in the ball indefinitely, as $N_{-}+N_{+}\leq 2N_{-}\leq 2$, by (1.15), so the solution cannot converge to $y$ as $\tau\to-\infty$. Now let $\tau_{k}$ be a subsequence such that $\varphi^{\tau_{k}}(x)\to y$, and let $s_{k}\leq\tau_{k}$ be the first time prior to $\tau_{k}$ such that $\varphi^{s_{k}}(x)\in\partial B_{\delta}(y)$. Then for $r\in(s_{k},\tau_{k}]$ we have $\varphi^{r}(x)\in B_{\delta}(y)$.

Therefore, by compactness, we may select a convergent subsequence $(\varphi^{s_{k_{j}}}(x))_{j}$ yielding an $\alpha$-limit point $z$ on $\partial B_{\delta}(y)$. It follows that $\Omega(z)=0$ and $N_{-}(z)=N_{+}(z)$ by estimates of the form (2.9), thus $z$ is a vacuum type I or II point.

To see that $y$ and $z$ lie in the closure of the same type II-vacuum orbit, we compute that the function $f$, defined by (2.13), satisfies the estimate

[TABLE]

where $M=\tfrac{4}{(\sqrt{3}-1)^{3}}$; we make use of the bounds provided by (1.15). Therefore we estimate that for $r\in(s_{k},\tau_{k})$ we have

[TABLE]

In particular, we estimate that

[TABLE]

As $\Omega(\tau_{k}),(N_{-}^{2}-N_{+}^{2})(\tau_{k})$ and $(N_{-}-N_{+})^{2}(\tau_{k})$ all converge to [math] as $k\to\infty$, we thus conclude that

[TABLE]

Therefore $y$ and $z$ lie in the same plane $\{x\in\mathbb{R}^{5}:f(x)=c\}$, so $y$ is the $\omega$-limit point of the type II vacuum orbit through $z$, see Section 2.2.

If the point $z$ itself lies on the Kasner circle, we are done, as $d(y,z)>0$ and $y$ and $z$ lie on the same type II vacuum orbit. Thus let us assume $z\notin\mathcal{K}$. Now as $\Omega(z)=0$, by (1.15) it follows that $N_{-}(z)\neq 0$, which then implies that $z\in B_{2}^{+}(\mathrm{II})|_{\Omega=0}$, since it must be of type II. Then the $\alpha$-limit point of the orbit through $z$, which is by definition the point $K(y)$, also lies in $\alpha(x)$, due to the $\alpha$-limit set being closed and invariant under the flow by Remark 1.4. ∎

2.3. Bianchi type II with matter

Each of the six type II invariant sets comes with a non-vacuum equilibrium. For $B_{2}^{+}(\mathrm{II})$ we have the equilibrium $P_{2}^{+}(\mathrm{II})$, as denoted in Table 2. Observe that the equilibrium lies in the locally rotationally symmetric invariant set

[TABLE]

With Remark 1.5 in mind, one may compute that, seen as part of the closure of the phase space, the matrix of the linearization of (1.14) near the equilibrium $P_{2}^{+}(\mathrm{II})$ has two eigenvalues with negative real part. We thus expect the unstable set in $B_{1}^{+}(\mathrm{VI}_{0})$ to be one-dimensional; by Proposition 4.1 it is in fact a submanifold of dimension one, actually just a single heteroclinic orbit.

Notation**.**

The unstable set of $P_{2}^{+}(\mathrm{II})$ within $B_{1}^{+}(\mathrm{VI}_{0})$ we denote by $\mathcal{P}^{+}_{\mathrm{VI}_{0}}$, and the unstable set of $P_{3}^{-}(\mathrm{II})$ by $\mathcal{P}^{-}_{\mathrm{VI}_{0}}$; the union of these two sets we write as $\mathcal{P}_{\mathrm{VI}_{0}}$, in line with Definition 4.3 of [9]. We denote the unstable set of $F$ within $B_{2}^{+}(\mathrm{II})$ by $\mathcal{F}^{+}_{\mathrm{II}}$ and the one within $B_{3}^{-}(\mathrm{II})$ by $\mathcal{F}^{-}_{\mathrm{II}}$, and the union of these sets is denoted by $\mathcal{F}_{\mathrm{II}}$.

We will need information about the asymptotics towards the past for orbits of type II. The proposition below is Proposition 9.1 in [9], after applying the right symmetry.

Proposition 2.10**.**

Let $x\in B_{2}^{+}(\mathrm{II})$ for $\gamma\in(\tfrac{2}{3},2)$. If $x\notin S_{2}^{+}(\mathrm{II})$, then $\alpha(x)$ consists of a single point in $\mathcal{K}_{1}\cup\mathcal{K}_{3}$. Else, either

1. (a)

$x$* is the equilibrium $P_{2}^{+}(\mathrm{II})$,* 2. (b)

$x$* lies in the unstable set $\mathcal{F}^{+}_{\mathrm{II}}$, or* 3. (c)

the solution $\varphi^{\tau}(x)$ converges to the Taub point $T_{2}$ as $\tau\to-\infty$.

Similar to Lemma 2.2, we have the following lemma regarding the equilibrium $P_{2}^{+}(\mathrm{II})$. This is Lemma 4.1 from [9] after applying the right symmetry.

Lemma 2.11**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})$ for $\gamma\in(\tfrac{2}{3},2)$ and assume that $P_{2}^{+}(\mathrm{II})\in\alpha(x)$, but the solution $\varphi^{\tau}(x)$ does not converge to $P_{2}^{+}(\mathrm{II})$ as $\tau\to-\infty$. Then $\alpha(x)$ contains a limit point of type II which is not $P_{2}^{+}(\mathrm{II})$.

2.4. The vacuum case

The last part of the boundary of the phase space is the vacuum case, which consists of the invariant set

[TABLE]

For a detailed study of the vacuum case for Bianchi type VI0, including asymptotic expansions with respect to the normalized time, we refer to [3]. Since the vacuum case is well-known already, we only state the relevant results on the asymptotics.

Lemma 2.12**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})|_{\Omega=0}$. As $\tau\to-\infty$ the solution $\varphi^{\tau}(x)$ converges to a point in $\mathcal{K}_{1}$, while $\omega(x)$ consists solely of $T_{1}$.

For a proof of the statement regarding the past asymptotics set we refer to Section 22.7 of [10], and regarding the $\omega$-limit set to the proof of Lemma 3.1 of [3].

3. The shear invariant set

The Bianchi type VI0 phase space also contains an invariant set in its interior; this is the shear invariant set which also contains the global sink $P_{1}^{+}(\mathrm{VI}_{0})$ from Table 2. The equilibrium has only eigenvalues with negative real part (with respect to the linearization in $B_{1}^{+}(\mathrm{VI}_{0})$) and thus it is a local attractor. The shear invariant set is defined as

[TABLE]

In case of Bianchi type VI0 there are no locally rotationally symmetric solutions, cf. Remark 1.2, but one can easily check that the set above is invariant under (1.14). Let us analyse the asymptotics in this set.

Proposition 3.1**.**

Let $x\in S_{1}^{+}(\mathrm{VI}_{0})\setminus\{P_{1}^{+}(\mathrm{VI}_{0})\}$ for $\gamma\in(\tfrac{2}{3},2)$. If $x$ does not lie in $\mathcal{F}_{\mathrm{VI}_{0}}$, then $\varphi^{\tau}(x)\to Q_{1}$ as $\tau\to-\infty$.

Proof.

Consider the function, given in the Appendix of Section 6 of [11],

[TABLE]

The function $Z_{4}$ satisfies

[TABLE]

and is thus strictly monotonically increasing for points in the invariant set

[TABLE]

Indeed, the derivative of $Z_{4}$ vanishes only if $\Sigma_{-}=0$ and $\Sigma_{+}=-\tfrac{1}{2}q^{*}$. However, if we also demand that the time derivatives $\Sigma_{\pm}^{\prime}$ vanish – in order to violate monotonicity – then we recover the equilibrium $P_{1}^{+}(\mathrm{VI}_{0})$, due to (1.14) and (1.15).

We may also apply the monotonicity principle to the set

[TABLE]

since it is invariant. Thus by Proposition 1.7 we may conclude that any $y\in\alpha(x)$ lies in the boundary of this set. We know $P_{1}^{+}(\mathrm{VI}_{0})\notin\alpha(x)$, as the equilibrium repels backward orbits. The $\alpha$-limit set is non-emtpy by Remark 1.4. Also, as the limit set is contained in the boundary, any $y\in\alpha(x)$ satisfies $\Omega(y)=0$ or $N_{-}(y)=0$. We distinguish the case that $N_{-}(y)>0$ and thus $\Omega(y)=0$ for some $y\in\alpha(x)$, and the case that $N_{-}(y)=0$ for all $y\in\alpha(x)$.

Case 1. Assume there is $y\in\alpha(x)$ satisfying $N_{-}(y)>0$ and $\Omega(y)=0$. By the invariance of the $\alpha$-limit set and that of the shear invariant set, $\alpha(x)$ contains the entire vacuum orbit $S_{1}^{+}(\mathrm{VI}_{0})|_{\Omega=0}$. (Note that this is just a single orbit.) The closedness of the $\alpha$-limit set then implies that both $Q_{1}$ and $T_{1}$ are in $\alpha(x)$, cf. Lemma 2.12. However, rewriting the expression

[TABLE]

using the constraint, we know that $\Sigma_{+}$ is strictly monotonically decreasing in the interval $(-\infty,s)$ if $\Sigma_{+}(s)>0$ for some time $s$. Such an $s$ must exist, since $Q_{1}\in\alpha(x)$. But then $T_{1}\notin\alpha(x)$, as $\lim_{\tau\to-\infty}\Sigma_{+}(\tau)>\Sigma_{+}(s)>0$, while at $T_{1}$ we have $\Sigma_{+}=-1$. Therefore the case is void.

Case 2. Now assume that for all $y\in\alpha(x)$ we have $N_{-}(y)=0$, so all $\alpha$-limit points are of Bianchi type I. If $F\in\alpha(x)$, then as $x\notin\mathcal{F}_{\mathrm{VI}_{0}}$, we may assume that there is $y\in\alpha(x)$ with $\Sigma_{+}(y)\neq 0$, by Lemma 2.2. If $F\notin\alpha(x)$, the previous conclusion is immediate by the constraint.

If $\Sigma_{+}(y)<0$, then after applying the flow and taking the closure we find $T_{1}\in\alpha(x)$. By Lemma 2.7 we know that there is another vacuum point $y\in\alpha(x)$ which is not $T_{1}$. This other point must be $Q_{1}$, as these are the only points in the closure of $S_{1}^{+}(\mathrm{VI}_{0})$ that satisfy $N_{-}=\Omega=0$. We may again apply the argument above to obtain a contradiction.

If $0<\Sigma_{+}(y)<1$, then by invariance and closedness of $\alpha(x)$, we find that both $Q_{1}$ and $F$ belong to $\alpha(x)$, by Lemma 2.1 and the fact that $B(\mathrm{I})\setminus\mathcal{K}$ is the stable manifold of $F$. The argument above again leads to a contradiction.

The only remaining case is that $\Sigma_{+}(y)=1$, which implies that $y=Q_{1}$. If we assume that $\alpha(x)$ contains any other points, by the above we would be able to obtain yet another contradiction. ∎

Using the monotone function $Z_{4}$, one can readily determine the $\omega$-limit set as well. We refer to Proposition 6.4 of [11]; the exponential rate follows readily from the fact that all the eigenvalues have negative real part.

Proposition 3.2**.**

For all $x\in B_{1}^{+}(\mathrm{VI}_{0})$ with $\Omega(x)>0$ for $\gamma\in(\tfrac{2}{3},2)$, $\varphi^{\tau}(x)$ converges exponentially to $P_{1}^{+}(\mathrm{VI}_{0})$ as $\tau\to\infty$.

Outside of the shear invariant set there is the monotone function $Z_{1}$ introduced in the Appendix of Section 6 of [12]. This will allow us to find the $\alpha$-limit set of these points, similar to Lemma 10.2 of [9] for the case of Bianchi $\mathrm{VII}_{0}$. The statement regarding the $\alpha$-limit set may be found in Section 3.3 of [6], but we prove a slightly stronger result.

Lemma 3.3**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})\setminus S_{1}^{+}(\mathrm{VI}_{0})$ for $\gamma\in(\tfrac{2}{3},2)$. Then the solution satisfies

[TABLE]

In particular, we have

[TABLE]

Proof.

Assume to the contrary, then by compactness of the closure of the phase space, we can find $y\in\alpha(x)$ such that $N_{-}^{2}(y)>N_{+}^{2}(y)$.

Consider the monotone function, as given in [12],

[TABLE]

This is well-defined and positive on the invariant set $B_{1}^{+}(\mathrm{VI}_{0})\setminus S_{1}^{+}(\mathrm{VI}_{0})$. It is indeed strictly monotonically decreasing; note that we have

[TABLE]

Firstly, by (1.15), we must have $\Sigma_{+}(\tau)>-1$, as $N_{-}(\tau)>0$ for any $\tau\in\mathbb{R}$. Secondly, if $\Sigma_{-}(\tau)=0$, we must have $\Sigma_{-}^{\prime}(\tau)\neq 0$, else we may show that $x\in S^{+}_{1}(\mathrm{VI}_{0})$ by invariance. Then by integrating it follows that $Z_{1}(\tau_{0})>Z_{1}(\tau_{1})$ for any $\tau_{0}<\tau_{1}$, so the function $Z_{1}$ is indeed strictly decreasing.

By Proposition 1.7 it follows that

[TABLE]

However, since $Z_{1}$ is decreasing and $Z_{1}(x)>0$, it cannot happen that $\Sigma_{-}^{2}(y)+N_{+}^{2}(y)=0$ if not also $N_{-}^{2}(y)=N_{+}^{2}(y)$. ∎

4. Regularity of the unstable sets of $F$, $P_{2}^{+}(\mathrm{II})$ and $P_{3}^{-}(\mathrm{II})$

For the strength of Theorem 1.6, it would be interesting to know whether the sets of points which converge towards the past to $F$ and $P_{2}^{+}(\mathrm{II}),P_{3}^{-}(\mathrm{II})$, are in fact smooth manifolds, so as to verify that they are sets of measure zero within the phase space. The ideas here are similar to those of Theorem 16.1 of [9], but here we make use of the information we gathered about $B_{1}^{+}(\mathrm{VI}_{0})$ in Proposition 3.2. The set $\mathcal{F}_{\mathrm{II}}$ is already covered in [9], so we omit it.

Proposition 4.1**.**

The unstable sets $\mathcal{P}^{+}_{\mathrm{VI}_{0}},\mathcal{P}^{-}_{\mathrm{VI}_{0}}$ and $\mathcal{F}_{\mathrm{VI}_{0}}$ are (embedded) smooth submanifolds of $\mathbb{R}^{4}$ of dimensions $1,1$ and $2$ respectively.

Proof.

We deal with the case of $\mathcal{F}_{\mathrm{VI}_{0}}$, the other two cases being similar and easier to handle, as they each consist of just a single heteroclinic orbit.111 Generally there are two one-dimensional unstable manifolds in case of a single positive eigenvalue, but due to the restrictions of our phase space, one of those orbits is not part of the phase space; for $F$ this comes from the restriction $N_{-}>0$, and for $P_{j}^{\pm}(\mathrm{II})$ it comes from the restriction $N_{-}>\absolutevalue{N_{+}}$.

We claim first that there is a neighbourhood $U$ of $F$ such that $U\cap\mathcal{F}_{\mathrm{VI}_{0}}=U\cap M_{F}$, where $M_{F}\subset B_{1}^{+}(\mathrm{VI}_{0})$ is the local unstable manifold of dimension two guaranteed by the Stable Manifold Theorem (see e.g. Theorem 1.3.2 of [2]). Let us assume such a $U$ does not exist; we may then select a sequence of points $(x_{k})_{k}\subset\mathcal{F}_{\mathrm{VI}_{0}}$, $x_{k}\notin M_{F}$, converging to $F$ as $k\to\infty$ and such that $\varphi^{\tau}(x_{k})$ converges to $F$ as $\tau\to-\infty$.

We may choose a neighbourhood $V$ of $F$ with smooth boundary $\partial V$, such that $M_{F}$ intersects $\partial V$ transversally. The resulting one-dimensional submanifold, $S:=\partial V\cap M_{F}$, may be assumed to be closed so in particular compact. Each point $y$ on $S$ uniquely defines a heteroclinic orbit $\mathcal{O}(y)\subset\mathcal{F}_{\mathrm{VI}_{0}}$, from $F$ to $P_{1}^{+}(\mathrm{VI}_{0})$, as $\Omega(y)>0,N_{-}(y)>0$ and thus $\varphi^{\tau}(y)\to P_{1}^{+}(\mathrm{VI}_{0})$ by Proposition 3.2. On the other hand, each orbit in $\mathcal{F}_{\mathrm{VI}_{0}}$ contains a single point $y\in S$.

Without loss of generality we may assume that $(x_{k})_{k}\subset V$. Now define for each $x_{k}$ the time $T_{k}$ as the smallest positive number such that $\varphi^{T_{k}}(x_{k})$ lies on the boundary $\partial V$. Also, define $y_{k}$ as the sole element of $\mathcal{O}(x_{k})\cap S$.

Choose a neighbourhood $W$ of $P_{1}^{+}(\mathrm{VI}_{0})$, not intersecting $V$, homeomorphic to a ball and such that orbits entering $W$ do not leave it. This is possible due to the Hartman-Grobman theorem (see e.g. Theorem 1.3.1 of [2]) and all eigenvalues of the linearization at $P_{1}^{+}(\mathrm{VI}_{0})$ having negative real part. We then assign to each $y\in S$ the time $\tau_{y}$ to reach $\partial W$. This assignment is continuous – since $(\tau,x)\mapsto\varphi^{\tau}(x)$ is smooth – and thus by compactness also bounded.

Then we estimate that

[TABLE]

In particular, the sequence $(T_{k})_{k}$ is bounded, and thus we may choose a convergent subsequence $T_{k_{j}}\to T_{0}$. Now note that

[TABLE]

as $F$ is an equilibrium. But as $\varphi^{T_{k_{j}}}(x_{k_{j}})\subset\partial W$, we may choose a convergent subsequence there to obtain a contradiction.

Given any $x\in\mathcal{F}_{\mathrm{VI}_{0}}$ we can find a time $s(x)$ such that $\varphi^{-s(x)}(x)\in U$. In particular, by invariance under the flow of the set $\mathcal{F}_{\mathrm{VI}_{0}}$, we find that $\varphi^{-s(x)}(x)\in U\cap M_{F}$. The diffeomorphisms $\varphi^{\tau}$ thus induce the structure of a smooth submanifold on $\mathcal{F}_{\mathrm{VI}_{0}}$ from the one on $M_{F}$. ∎

Knowing that the unstable manifolds of $F,P_{2}^{+}(\mathrm{II})$ and $P_{3}^{-}(\mathrm{II})$ are small, it makes sense to define the points in the phase space outside of these sets, and not being the single point $P_{1}^{+}(\mathrm{VI}_{0})$, as generic. We define it for the case of $B_{1}^{+}(\mathrm{VI}_{0})$, but one can obtain similar definitions in the other cases $B_{j}^{\pm}(\mathrm{VI}_{0})$ by applying the appropriate symmetries.

Definition 4.2**.**

We call a point $x\in B_{1}^{+}(\mathrm{VI}_{0})$ generic if it is not the equilibrium $P_{1}^{+}(\mathrm{VI}_{0})$ and it is not contained in the unstable manifold $\mathcal{P}_{\mathrm{VI}_{0}}^{+}$, $\mathcal{P}_{\mathrm{VI}_{0}}^{-}$ or $\mathcal{F}_{\mathrm{VI}_{0}}$.

5. Convergence to $\mathcal{K}_{1}$ for generic points

We can now study the limit behaviour for generic points outside of the shear invariant set, applying a proof similar to that of Proposition 10.2 of [9]. This result, along with Propositions 3.1 and 4.1, proves Theorem 1.6.

Proposition 5.1**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})\setminus S_{1}^{+}(\mathrm{VI}_{0})$ be generic, for $\gamma\in(\tfrac{2}{3},2)$. Then $\varphi^{\tau}(x)$ converges to a point in $\mathcal{K}_{1}$.

Proof.

Firstly, by Remark 1.4 the $\alpha$-limit set $\alpha(x)$ is non-emtpy. By Lemma 3.3, we know that

[TABLE]

We begin by showing that $\alpha(x)\cap\mathcal{K}\neq\emptyset$. We distinguish the case that there is a limit point of type I, and the case that there is a limit point of type II.

Case 1. Assume there is a limit point $y\in\alpha(x)$ of type I. If $y\neq F$, then by Lemma 2.1 we find an $\alpha$-limit point in $\mathcal{K}$. But as convergence to $F$ was excluded by assuming $x\notin\mathcal{F}_{\mathrm{VI}_{0}}$, we know that if $y=F$, then Lemma 2.2 applies, so that $\alpha(x)$ contains a point in $\mathcal{K}$.

Case 2. If there exists a $y\in\alpha(x)$ is of type II, we may assume $y\notin\mathcal{F}_{\mathrm{II}}$, else we are back in the case above, by invariance and closedness of $\alpha(x)$. Moreover, we may assume that $y$ is not one of the equilibria $P_{2}^{+}(\mathrm{II})$ or $P_{3}^{-}(\mathrm{II})$, by Lemma 2.11. Note again that $x$ cannot converge to these equilibria, since we assumed $x\notin\mathcal{P}_{\mathrm{VI}_{0}}$. Then, by Proposition 2.10, $y\in\mathcal{K}$.

Now we shall show that $\alpha(x)$ contains an element of $\overline{\mathcal{K}_{1}}$. Let $y\in\alpha(x)\cap\mathcal{K}$. In case that $y\in\mathcal{K}_{2}\cup\mathcal{K}_{3}$, we may apply the Kasner map a finite number of times to find an $\alpha$-limit point in $\overline{\mathcal{K}_{1}}$ by Lemma 2.9. In case that $y=T_{1}$, we use Lemma 2.7 to find $z\in\alpha(x)\cap\{\Omega=0\}$ which is not $T_{1}$. If $z$ is of type I, then it must lie on the Kasner circle and since it is not $T_{1}$ we are done by the case above. If $z$ is of type II, then by closure and invariance of $\alpha(x)$ there exists $w\in\mathcal{K}$ such that both $w$ and $K(w)\neq w$ are contained in $\alpha(x)$. Since at most one of these is $T_{1}$, we are again done by the case above.

Next we claim that $\alpha(x)$ consists of only a single point $z\in\mathcal{K}_{1}$. Write $N_{-}^{2}=1-\Omega-\Sigma_{+}^{2}-\Sigma_{-}^{2}$ using the constraint. Then we obtain that

[TABLE]

so in particular: if $\Sigma_{+}>0$, then it is strictly monotonically decreasing and thus has a unique limit as $\tau\to-\infty$. Moreover, $\Omega$ must converge to [math], else $\Sigma_{+}$ becomes unbounded thus violating (1.15). Indeed, we have the inequality

[TABLE]

and by Lemma 1.8(ii), we conclude that

[TABLE]

for any $\tau<\tau_{0}$, where $\tau_{0}$ is any time such that $\Sigma_{+}(\tau_{0})>0$. In particular, $\int^{\tau_{0}}_{\tau}\Omega(s)ds<\infty$, so $\Omega$ is integrable. Since the function $\Omega(\tau)$ is continuous, positive and has bounded derivative – recall that the closure of phase space is compact and the vector field polynomial – we thus have $\lim_{\tau\to-\infty}\Omega(\tau)=0$. In the same way $\Sigma_{+}^{2}+\Sigma_{-}^{2}$ converges to $1$. Then, by the constraint, also $N_{-}$ and thus also $N_{+}$ converge to [math]. The limit of all coordinates is thus unique, up to the sign of $\Sigma_{-}$. However, by connectedness of $\alpha(x)$, we conclude that this sign is uniquely determined. Therefore $\varphi^{\tau}(x)$ converges to a single point in $\overline{\mathcal{K}}_{1}$ as $\tau\to-\infty$. We know from Corollary 2.6 that $x\in B_{1}^{+}(\mathrm{VI}_{0})$ cannot converge to $T_{2}$ or $T_{3}$, so we conclude that $z\in\mathcal{K}_{1}$. ∎

For the subsequent sections, the following will be of use.

Corollary 5.2**.**

Let $x\in B_{1}^{+}(\mathrm{VI}_{0})$ be generic, for $\gamma\in(\tfrac{2}{3},2)$. Then $2-q(\tau)$ is positive and it converges to [math] exponentially as $\tau\to-\infty$.

Proof.

We know that for a generic $x\in B_{1}^{+}(\mathrm{VI}_{0})$ the solution $\varphi^{\tau}(x)$ converges to a point in $\mathcal{K}_{1}$, by Propositions 3.1 and 5.1. In particular, $q(\tau)\to 2$ as $\tau\to-\infty$, so we only have to show the positivity and the exponential rate.

The exponential rate follows from the exponential rates at which $N_{-}$ and $\Omega$ converge to 0, combined with the constraint. We have

[TABLE]

Since the factor in front of $N_{-}$ above is positive everywhere on $\mathcal{K}_{1}$, by Remark 2.4, there are $s=s(x,\gamma)\in\mathbb{R}$ and $\eta=\eta(x,\gamma)\in(0,2-q^{*})$ such that for any $\tau\leq s$ we have

[TABLE]

This follows from Grönwall’s lemma.

Regarding $\Omega$, we have

[TABLE]

Thus select a time $s$ – which is possibly different from the one above, but take $s$ to be the smaller one between them – such that for $\tau\leq s$ we have $q(\tau)-q^{*}>\eta$. This is possible as $q\to 2$ and $q^{*}\in(0,2)$ for $\gamma\in(\tfrac{2}{3},2)$. Then it is clear that also

[TABLE]

By the constraint equation 1.15, we can write

[TABLE]

The positivity follows now from the function $N_{-}$ not vanishing in the phase space. By the equation above, we find that, for $\tau\leq s$, there must be a constant $C=C(x,\gamma)>0$ such that

[TABLE]

Note that $C$ and $\eta$ are independent of $\tau$, they only depend on $\gamma$ and the initial data $x$. ∎

6. Bianchi type VI0 developments

In this section we introduce some concepts which are necessary to translate Theorem 1.6 to the geometrical setting. This allows us to apply the theorem to the Klein-Gordon equation, by making use of the machinery developed in [8], an article on asymptotics of solutions to the Klein-Gordon equation on Bianchi backgrounds.

The material of this section is not new; we refer the reader to [8], [9] and [12] in particular. We refrain from going into details too much and restrict ourselves to introducing the concepts necessary to state and prove Theorem 7.1.

6.1. Orthogonal perfect fluid data

Recalling Definition 1.1, consider the induced metric $g_{t}~{}=~{}h|_{TG_{t}\times TG_{t}}$ on $G_{t}:=\{t\}\times G$. As $(G_{t},g_{t})$ is a spacelike submanifold, orthogonal to $\partial_{t}$, we have, dropping the $t$-dependencies,

[TABLE]

as well as

[TABLE]

These are the Hamiltonian and momentum constraint equations for an orthogonal perfect fluid, respectively. Here $k$, $S$ and $D$ are respectively the second fundamental form, scalar curvature and Levi-Civita connection of $(\{t\}\times G,g_{t})$.

Taking a different perspective, we may also specify the initial data to satisfy these equations, cf. Definition 9 of [8].

Definition 6.1**.**

Bianchi orthogonal perfect fluid data $(G,g,k,\rho_{0})$ consist of a connected three-dimensional Lie group $G$, a left invariant metric $g$ on $G$, a left invariant, symmetric, covariant two-tensor field $k$ on $G$, and a non-negative constant $\rho_{0}$, satisfying equations (6.1) and (6.2).

Fluid data of class A, in the sense above, each give a point in the phase space of the system (1.6). For this we refer to Lemma 21.2 of [9] and Section 17.1 of [8]. By permutation we may assume that Bianchi VI0 orthogonal fluid data correspond to a point in $B_{1}^{+}(\mathrm{VI}_{0})$ without loss of generality. Conversely, see p. 1419 of [12], the orbit of a point in the phase space of (1.6) corresponds to a one-parameter family of conformally-related developments. This is due to an arbitrariness in determining the length scale in the cosmological model.

Definition 6.2**.**

Given Bianchi VI0 orthogonal perfect fluid data, we call it generic type $\mathrm{VI}_{0}$ data if the corresponding point $x\in B_{1}^{+}(\mathrm{VI}_{0})$ is generic in the sense of Definition 4.2.

6.2. The initial singularity

We are interested in Bianchi spacetimes with a certain singularity towards the past. The type of singularity is defined precisely in Definition 2 of [8], which we copy for convenience.

Definition 6.3**.**

Let $(M,h)$ be a Bianchi spacetime. Then $t_{-}$ is a monotone volume singularity if there exists a time $t_{0}\in I$ such that the mean curvature $\theta(t)$ of the hypersurface $G\times\{t\}$ is strictly positive on the interval $(t_{-},t_{0})$, and if the function $\tau$, defined as

[TABLE]

satisfies $\tau(t)\to-\infty$ as $t\to t_{-}$.

Remark 6.4**.**

The function $\tau$ is known as the logarithmic volume density, and satisfies $\tau^{\prime}(t)=\theta/3$. In what follows we employ $\tau$ as a new time-coordinate. However, it is also precisely the time coordinate used in the Wainwright-Hsu equations.

In Theorem 7.1 we also encounter the time coordinate $\sigma$, defined implicitly by

[TABLE]

together with the condition that $\sigma(t)=0$ for the same time $t$ for which $\tau(t)=0$. We note that if the Lie group is of type VI0 then $\sigma\to-\infty$ corresponds to $\tau\to-\infty$, see Lemma 52 of [8]. ∎

Remark 6.5**.**

In the present context, one may write the Raychaudhuri equation as

[TABLE]

where $q$ is the deceleration parameter, which coincides with the function $q(\tau)$ introduced in Section 1.4. Recall that $2-q(\tau)$ is positive and converges exponentially to [math] for generic $x\in B_{1}^{+}(\mathrm{VI}_{0})$, by Corollary 5.2. In particular, $2-q(\tau)$ is an integrable function. It follows that for generic type $\mathrm{VI}_{0}$ data we have

[TABLE]

where $\tau<0$ and

[TABLE]

Note that $C$ depends on both $\gamma$ and the data. ∎

6.3. Class A developments

Given Bianchi class A orthogonal perfect fluid data $(G,g,k,\rho_{0})$, there is a corresponding development, which is in fact a Bianchi spacetime. This is the content of Lemma 21.2 of [9]; this development is called the corresponding class A development.

Remark 6.6**.**

In case that the underlying Lie group is of type $\mathrm{VI}_{0}$, the solution always has a monotone volume singularity. Indeed, for Bianchi $\mathrm{VI}_{0}$, the lack of locally rotationally symmetric developments, see Remark 1.2, implies that all of the corresponding developments are past causally geodesically incomplete and future causally geodescially complete, by Lemmas 21.5 and 21.8 of [9]. By Lemma 21.8 of [9], for type $\mathrm{VI}_{0}$ the interval $I$ is of the form $I=(t_{-},\infty)$, with $t_{-}>-\infty$. In the new time coordinate $\tau$ this interval gets mapped to $\mathbb{R}$, by Lemma 22.4 of [9], so we always have a monotone volume singularity for Bianchi type $\mathrm{VI}_{0}$ with OPF matter. ∎

Remark 6.7**.**

Let us note that for a class A development $(M,h)$ not of type IX, the scalar curvature $S$ of constant$-t$ hypersurfaces $(M_{t},g_{t})$ is nonpositive; for those Bianchi types the function $\gamma_{S}:=\max\{0,S\}/\theta^{2}$ introduced in Theorem 7 of [8] is identically zero, see the remarks after Remark 8 of [8]. ∎

7. Klein-Gordon equation on Bianchi type VI0 backgrounds.

We now tackle the Klein-Gordon equation on type $\mathrm{VI}_{0}$ backgrounds, by combining results from [8] with Theorem 1.6 above. This fills a gap in the results of [8], specifically in Example 29 and Section 15.2, which was due to lack of proven asymptotics for type $\mathrm{VI}_{0}$ in the Wainwright-Hsu variables.

Given a Lorentz manifold $(M,h)$, the Klein-Gordon equation reads

[TABLE]

where $m\in\mathbb{R}$ is a constant. Now we are in a position to state our theorem. In the statement below $\langle\sigma\rangle:=\sqrt{1+\sigma^{2}}$ denotes the Japanese bracket.

Theorem 7.1**.**

Let $(G,g,k,\rho_{0})$ be Bianchi orthogonal perfect fluid data for $\gamma\in(\tfrac{2}{3},2)$, with $G$ of type $\mathrm{VI}_{0}$ and with $\rho_{0}>0$. Assume that the data are generic, and denote the corresponding development by $(M,h)$.

Then there exist constants $\lambda,\nu>0$, such that for any smooth solution $u\in C^{\infty}(M)$ of the Klein-Gordon equation (7.1) on $(M,h)$, there are smooth functions $v,w,\phi,\psi\in C^{\infty}(G)$, such that for each compact subset $K\subseteq G$ and for each $l\in\mathbb{N}_{0}$, we have

[TABLE]

for all $\tau\leq 0$, as well as

[TABLE]

for all $\sigma\leq 0$.

Proof.

We wish to apply Proposition 19 of [8]. From Sections 1.3 and 6, it is clear that the following is satisfied: $(M,h)$ has a monotone volume singularity $t_{-}$, by Remark 6.6; the metric $h$ solves Einstein’s equation (1.3); $\rho\geq p$, as well as $\rho\geq 0$ and $\Lambda\geq 0$, by Remark 1.3, where we keep in mind that if $\rho_{0}>0$ then $\rho>0$ throughout the development. Moreover, for type $\mathrm{VI}_{0}$ spacetimes, the scalar curvature $S$ of constant$-t$ hypersurfaces $(M_{t},g_{t})$ is nonpositive, implying that $\gamma_{S}\equiv 0$, by Remark 6.7.

Note that the function $\varphi_{0}$ is simply $-m^{2}$ in our setting, and thus, by Remark 6.5,

[TABLE]

where $C>0$ is a constant depending on the value of $\gamma$ and the initial data. Moreover, we know that for generic data $2-q(\tau)$ decays exponentially, by Corollary 5.2. Lastly, we know that $\norm{\hat{a}^{-1}(\tau)}$ decays exponentially. Indeed, for generic data as in Definition 6.2, the corresponding backward orbit of (1.6) converges to a point on one of the Kasner arcs $K_{j},j=1,2,3$, by Theorem 1.6. From Section 17.1 of [8], we gather that this implies the desired estimate for $\norm{\hat{a}^{-1}(\tau)}$.

We have satisfied all assumptions of Proposition 19 of [8], establishing the proof. ∎

7.1. Non-generic data

For developments corresponding to case $(b)$ of Theorem 1.6, asymptotics to solutions of the Klein-Gordon equation are discussed in Example 35 of [8]. Therefore, only case $(a)$ remains.

For the development corresponding to the equilibrium $P_{1}^{+}(\mathrm{VI}_{0})$ we may compute that

[TABLE]

where $C_{a,1}>0$ is a constant, cf. Section 14.1 of [8]. Since $p^{*}\in(0,1)$ for $\gamma\in(\tfrac{2}{3},2)$, we find that the norm $\norm{\hat{a}^{-1}}$ remains bounded away from zero. Thus we know that

[TABLE]

meaning the singularity is not silent. Therefore we cannot apply Proposition 19 of [8] to find the desired asymptotics.

Acknowledgements

The author sincerely thanks Hans Ringström for his many comments and suggestions, the anonymous reviewer for their careful reading of the manuscript, and Bernhard Brehm for his lecture series on the dynamical systems approach. This research was funded by the Swedish research council, dnr. 2017-03863.