Abstract
This article introduces the notions of asymptotic dust and asymptotic radiation equations of state. With these non-linear generalizations of the well known dust or (incoherent) radiation equations of state the perfect-fluid equations lose any conformal covariance or privilege. We analyse the conformal field equations induced with these equations of state. It is shown that the Einstein-\(\lambda \)-perfect-fluid equations with an asymptotic radiation equation of state allow for large sets of Cauchy data that develop into solutions which admit smooth conformal boundaries in the future and smooth extensions beyond. In the case of asymptotic dust equations of state sharp results on the future asymptotic behaviour are not available yet.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
This article extends our investigation of the long term behaviour of solutions to Einstein’s equation
with positive cosmological constant \(\lambda \) and energy-momentum tensor \( {\hat{T}}_{\mu \nu }\) for various matter models. The ones considered in the present paper are given by perfect fluids with certain non-linear equations of state. Several authors studied the future asymptotic behaviour of solutions to the Einstein-\(\lambda \)-perfect fluid equations with a linear equation of state \({\hat{p}} = w({\hat{\rho }}) = w_*\,{\hat{\rho }}\), \(w_* = const.\) [5, 9, 12, 13, 15, 16], often assuming that \(0< w_* < 1/3\). More general equations of state have been considered in the articles [11, 14]. In the present article the asymptotic behaviour is analysed from a particular point of view. It is motivated by the following observation.
De Sitter space is a geodesically complete, spatially compact, conformally flat solution of FLRW type to the Einstein-\(\lambda \)-vacuum equations with cosmological constant \(\lambda > 0\). It admits smooth conformal boundaries at future and past time-like infinity. It generalizes as follows. Smooth Cauchy data for the same equations on a Cauchy hypersurface S of de Sitter space that are sufficiently close (with respect to suitable Sobolev norms) to the de Sitter data develop into solutions that are also time-like and null geodesically complete and admit smooth conformal boundaries. The conformal Einstein equations, which have been used for this purpose, determine in fact smooth solutions that extend beyond these boundaries and define on ‘the other sides’ again vacuum solutions. Perturbations of the conformal curvature, only restricted by the smallness condition on S, travel unimpeded across the boundaries [1, 2].
The question of interest in this article is: For which matter models can be established similar results on the future asymptotic behaviour of solutions to the Einstein-\(\lambda \)-matter equations ? For a discussion of the answers obtained so far the reader is referred to the article [6]. The latter can also be read as an extended non-technical introduction to the present article. It explains our interest in the question above and also may be helpful for readers not acquainted with the conformal methods used in the following.
The two simplest cases with a positive answer are the FLRW-type solutions to the Einstein-\(\lambda \)-perfect-fluid equations with the pure (incoherent) radiation equation of state \({\hat{p}} = w({\hat{\rho }}) = \frac{1}{3}\,{\hat{\rho }}\) and the pure dust equation of state \({\hat{p}} = w({\hat{\rho }}) = 0\) (the unusual word pure has been added here to avoid confusions with notions introduced below).
These cases have been generalized by showing that these FLRW solutions are future stable in the class of general solutions to the Einstein-\(\lambda \)-perfect-fluid equations with the respective equation of state. The perturbed solutions are not only time-like geodesically future complete but also admit smooth future conformal boundaries and extensions beyond [5, 12]. These cases are, however, still somewhat special. The perfect fluid equations with a pure radiation equation of state considered in [12] are conformally invariant with vanishing trace of the energy momentum tensor. The perfect fluid equations with a pure dust equation of state are not conformally invariant and have an energy momentum tensor with non-vanishing trace. They are, however, conformally privileged (see [6]) by being related to a conformally invariant structure which helps establish the result of [5].
It may be reasonable to base models of the universe on solutions to the Einstein-\(\lambda \)-perfect-fluid equations where the future development is determined by one of the equations of state above. It appears unlikely, however, that it makes sense to cover the whole stretch from the Big Bang to future time-like infinity by one and the same linear equation of state. One can expect to need a transition where
with some function \( w^{**}({\hat{\rho }})\) of which we will assume that \(0 \le w^{**}({\hat{\rho }}) \le 1/3\). At late time \( w^{**}({\hat{\rho }})\) should then acquire the value \(w^{**}({\hat{\rho }}) = 0\) if we want to model the end by a pure dust equation of state and \(w^{**}({\hat{\rho }}) = 0\) in the case of pure radiation.
There is nothing, however, which fixes a natural notion of ‘late time’. The only meaningful requirement would be that these values are assumed in the limit when the space-time approaches future time-like infinity. But then the equation of state would still need to recognize where and when this limit will be achieved.
In the conformal analysis of the two FLRW models mentioned above the physical density \({\hat{\rho }}\) and the conformal density \(\rho \) satisfy a relation of the form
with \(e = 3\) in the case of pure dust and \(e = 4\) in the case of pure radiation. These values are chosen because they imply that \(\rho = \rho _* = const. > 0\) (assuming that \({\hat{\rho }} > 0\) on some initial slice). Since then \({\hat{\rho }} \rightarrow 0\) at future time-like infinity where \(\Omega \rightarrow 0\), the behaviour of \({\hat{\rho }}\) can be understood as an indicator for the approach to future time-like infinity. In generalizing the situation we shall keep (1.2), hoping it to serve as an indicator function of the far future where \(\Omega \rightarrow 0\). The relation \(\rho = \rho _* = const.\) may not be preserved under the generalization considered below and we will have to control that \(\rho \) remains positive and bounded in the limit \(\Omega \rightarrow 0\). The pure dust and pure radiation equations of state are now generalized as follows.
An asymptotic dust equation of state is given by a function of the form
combined with (1.2) where \(e = 3\). It implies \(w'({\hat{\rho }}) = (1 + k) \,{\hat{\rho }}^{k}\,\,w^*({\hat{\rho }}) + {\hat{\rho }}^{1 + k}\,\,(w^*)'({\hat{\rho }})\), where here and in the following the notation \( ' = \partial /\partial {\hat{\rho }}\) is used.
An asymptotic radiation equation of state is given by a function of the form
combined with (1.2) where \(e = 4\). It implies \(w'({\hat{\rho }}) = \frac{1}{3} - (1 + k) \,{\hat{\rho }}^{k}\,\,w^*({\hat{\rho }}) - {\hat{\rho }}^{1 + k}\,\,(w^*)'({\hat{\rho }})\). In both cases \(w^*({\hat{\rho }})\) is a smooth function of \({\hat{\rho }}\) that satisfies
To ensure that \(0 \le w^{**}({\hat{\rho }}) \le 1/3\) as required above, we would need to impose more detailed conditions on \(w^*({\hat{\rho }})\) but for the analysis of the effect of the new equations of state in the far future where \(\Omega \) becomes small the conditions above are sufficient. The limits \(w^*({\hat{\rho }}) \rightarrow 0\) give back the pure dust and the pure radiation equations of state.
The factor \({\hat{\rho }}^k\) has been inserted in the definitions as a simple means to control in terms of k the speed at which the pure dust or the pure radiation situations is approximated as \({\hat{\rho }} \rightarrow 0\). Eventually the conformal field equations may force us to impose less crude conditions on \(w^*({\hat{\rho }})\) or to introduce more subtle definitions of asymptotic equations of state.
Definition (1.3) implies
while definition (1.4) gives
In both case \({\hat{T}} \ne 0\) if \(\rho > 0\) and \({\hat{T}} \rightarrow 0\) as \(\Omega \rightarrow 0\) if \(\rho \) remains bounded in this limit.
Because the principal parts of the matter equations are affected by the equations of state above any conformal covariance or privilege is lost.
The Cauchy problem local in time for Einstein-\(\lambda \)-perfect fluids with an asymptotic dust or radiation equations of state as above poses no problems. This follows from the results of [3, 8] where only weak conditions on the equation of state are assumed.
First conditions on the admissible values of k are found if the equations of state above are considered in the conformal analogues of the FLRW Friedmann and the energy conservation equation. In the case of asymptotic dust the system reads
where \(R[{\hat{h}}] = const. \ge 0\). It can be integrated across \(\Omega = 0\) with \(\rho \) bounded and positive if \(3\,k - 1 \in {\mathbb {N}}_0\). In the case of asymptotic radiation the system reads
It can be integrated across \(\Omega = 0\) with \(\rho \) bounded and positive if \(4\,k - 1 \in {\mathbb {N}}_0\).
In the following sections we derive the conformal Einstein-\(\lambda \)-perfect fluid equations, introduce a gauge that involves in particular an orthonormal frame \(e_j\), \(j = 0, 1, 2, 3\), with \(e_0 = U\), and discuss the regularity as \(\Omega \rightarrow 0\) of the equations in this gauge with any of the two asymptotic equations of state assumed.
In the case of an asymptotic dust equation of state there arise problems. The use of this equation of state affects the principal part of the fluid equations in a way which does not allow us to apply the methods which were successful in the case of the pure dust equation of state [5]. This case, which I consider as particularly interesting (see [6]), is left open.
In the case of an asymptotic radiation equation of state regularity of the equations can be ensured by a suitable assumption and there can be derived a symmetric hyperbolic reduced system that is well defined irrespective of the sign of \(\Omega \). This allows us to give the following answer to the question posed in the beginning.
Consider the reduced system of the conformal Einstein-\(\lambda \)-perfect-fluid equations in the gauge discussed below with an asymptotic radiation equation of state where \(k \ge 1\).On a compact 3-dimensional manifold \({{\mathcal {J}}}^+\) let be given smooth Cauchy data for the reduced equations with \(\Omega = 0\), \(U^k = \delta ^k\,_0\) time-like future directed, and \(U^k\,\nabla _k \Omega < 0\) that satisfy the constraints induced by the conformal field equations and the special properties implied on a space-like hypersurface on which \(\Omega = 0\). These data determine a smooth solution to the reduced equations with \(\Omega < 0\)in the future of \({{\mathcal {J}}}^+\) and \(\Omega > 0\) in the past of \({{\mathcal {J}}}^+\). Where \(\Omega > 0\) the solution defines a unique solution to the system of the Einstein-\(\lambda \)-perfect-fluid equations with an asymptotic radiation equation of state that is time-like geodesically future complete and for which \({{\mathcal {J}}}^+\) represents the conformal boundary at the infinite time-like future.
Let Sbe a Cauchy hypersurface for this solution in the past of \({{\mathcal {J}}}^+\) and denote by \(\Delta \)the Cauchy data induced by the solution on S. Any Cauchy data \(\Delta '\)on Sfor the same system which are sufficiently close to \(\Delta \)develop into a solution which is also time-like geodesically future complete, admits a smooth conformal boundary in its future, and a smooth conformal extension beyond.
At the end of this article are discussed observations about solutions for which U is orthogonal to the hypersurface \({{\mathcal {J}}}^+\) and consequently hypersurface orthogonal on physical solution in the past of \({{\mathcal {J}}}^+\).
2 The Conformal Einstein-\(\lambda \)-Euler System
We consider the Einstein-\(\lambda \)-perfect fluid system with cosmological constant \(\lambda \) that consists of the Einstein equations (1.1) with cosmological constant \(\lambda > 0\) and an energy momentum tensor of a simple ideal fluid given by
The unknowns are a Lorentz metric \({\hat{g}}_{\mu \nu }\) on a four-dimensional manifold \({\hat{M}}\), a future directed time-like fluid flow vector field \({\hat{U}}^{\mu }\) satisfying \({\hat{U}}_{\mu }\,{\hat{U}}^{\mu } = - 1\), the total energy density \({\hat{\rho }}\) and the pressure \({\hat{p}}\) as measured by co-moving observers. The system is completed by an equation of state
The matter equations, equivalent to \({\hat{\nabla }}^{\mu }\,{\hat{T}}_{\mu \nu } = 0\), are given by
We shall assume that \({\hat{\rho }} > 0\) on an initial Cauchy slice.
Let \({\hat{e}}_a\), \(a = 1, 2, 3\), be vector fields that satisfy \({\hat{g}}({\hat{U}}, {\hat{e}}_a) = 0\), \({\hat{g}}( {\hat{e}}_a, {\hat{e}}_b) = \delta _{ab}\) and \(F = F({\hat{\rho }})\) a function that satisfies
where we use, as in the following chapters, the notation \(' = \frac{\partial }{\partial {\hat{\rho }}}\). With (2.3) follows
Taking a derivative and observing that \({\hat{\rho }} + w > 0\) we get the equation
Where \(w' \ne 0\) anti-symmetrization of this relation implies the integrability condition
Conformal transformations of the curvature fields.
The conformal Weyl tensor \(C^{\mu }\,_{\nu \lambda \eta }\) and the Schouten tensor \(L_{\mu \nu } = 1/2\Big (R_{\mu \nu } - 1/6\,R\, g_{\mu \nu }\Big )\) of a given metric satisfy the contracted Bianchi identity
Under the rescaling
the conformal Weyl tensor is invariant, \(C^{\mu }\,_{\nu \lambda \eta } = {\hat{C}}^{\mu }\,_{\nu \lambda \eta }\), the Ricci tensor transform as
and the Schouten tensor as
The conformal Weyl tensor satisfies \(\nabla _{\mu }\left( \Omega ^{-1}\, C^{\mu }\,_{\nu \lambda \eta }\right) = \Omega ^{-1}\,{\hat{\nabla }}_{\mu }\,{\hat{C}}^{\mu }\,_{\nu \lambda \eta }\) where \(\nabla \) and \({\hat{\nabla }}\) denote the Levi-Civita connections of \(g_{\mu \nu }\) and \({\hat{g}}_{\mu \nu }\). The identity above thus implies for the rescaled conformal Weyl tensor \(W^{\mu }\,_{\nu \lambda \eta } = \Omega ^{-1}\, C^{\mu }\,_{\nu \lambda \eta }\) the relation
whence
Conformal transformation of the matter fields.
We combine the conformal rescaling of the metric with the rescalings
In the following it will always be understood that
and similarly with the function \(w^*\) introduced later. Equations (1.1) imply
and \({\hat{L}}_{\mu \nu } = {\hat{L}}^*_{\mu \nu } + \frac{1}{6}\,\lambda \,{\hat{g}}_{\mu \nu }\) with
It follows that \({\hat{\nabla }}_{\lambda }{\hat{L}}_{\eta \nu } = {\hat{\nabla }}_{\lambda }{\hat{L}}^*_{\eta \nu }\). The Ricci scalar of the conformal metric satisfies
Written in the form
it will be referred to as \(\lambda \)-equation. Equations (2.8), (2.11), (2.12) give
which will be referred to as \(\Omega \)-equation. Applying a derivative to the \(\lambda \)-equation (2.12) gives with (2.13) the s-equation
It follows
which gives with (2.11) and \(\nabla _{\lambda }{\hat{\rho }} = \Omega ^{e}\,(\nabla _{\lambda }\rho + e\,\rho \,\Omega ^{-1}\,\nabla _{\lambda }\Omega )\) the relation
The conformal matter equations
Equations (2.2) and (2.3) transform with \(\nabla _{\mu }{\hat{\rho }} = \Omega ^{e}\left\{ \nabla _{\mu }\,\rho + e\,\rho \,\Omega ^{-1}\,\nabla _{\mu }\Omega \right\} \) into
In the case of pure dust, where \(w = 0\) and \(e = 3\), the equations read
The \(\Omega ^{-1}\) term in the second equation reflects the conformal non-covariance of the system. In the case of pure radiation with \(w = \frac{1}{3}\,{\hat{\rho }}\), and \(e = 4\) the conformal equations reduce to
and have thus the same form as their physical versions.
Equation (2.7) transforms into
With \({\hat{e}}^{\mu }\,_a = \Omega \,e^{\mu }\,_a\) the relation (2.5) transforms into
3 The Gauge and the Implied Equations
We express the equations in terms of an orthonormal frame field \(e_k =e^{\mu }\,_k\partial _{x^{\mu }}\), \(k = 0, 1, 2, 3\), so that \(g_{jk} \equiv g(e_j, e_k) = \eta _{jk} = diag(-1, 1, 1, 1)\) and \(e_0\) is a time-like vector field. The space-like frame vector fields are then given by the \(e_a\), where \(a, b, c = 1, 2, 3\) denote spatial indices to which the summation convention also applies. In the following all tensor fields considered in the previous sections will be expressed in terms of this frame field. The contravariant coordinate version of the metric is given by \(g^{\mu \nu } = \eta ^{jk}\,e^{\mu }\,_j\,e^{\nu }\,_k\).
The connection coefficients, defined by \(\nabla _je_k \equiv \nabla _{e_j}e_k = \Gamma _j\,^l\,_k\,e_l\), satisfy \( \Gamma _{jlk} = - \Gamma _{jkl}\) with \( \Gamma _{jlk} = \Gamma _j\,^i\,_k\,g_{i l}\), because \(\nabla _i g_{jk} = 0\). The covariant derivative of a tensor field \(X^{\mu }\,_{\nu }\), given in the frame by \(X^i\,_j\), takes the form \(\nabla _k\,X^i\,_j = X^i\,_{j\,,\mu }\,e^{\mu }\,_k + \Gamma _k\,^i\,_l \,X^l\,_j - \Gamma _k\,^i\,_l\,X^i\,_j\). The frame and the connection coefficients satisfy the first structural equations
which ensures that the connection is torsion free, and the second structural equations
To restrict the gauge freedom for the frame we set \(e_0 = U\) so that \(U = U^k\,e_k\) with \(U^k = \delta ^k\,_0\), choose at the points of a given smooth space-like hypersurface S transverse to the flow line of U vector fields \(e_a\), \(a = 1, 2, 3\), so that \(g(e_j, e_k) = \eta _{jk}\), assume the orthonormal frame \(e_k\) to be extended by Fermi-transport in the direction of \(e_0 = U\) so that \(0 = {\mathbb {F}}_U e_k = \nabla _U e_k - g(e_k, \nabla _U U)\,U + g(e_k, U)\, \nabla _U U\), and assume the \(e_a\) to be chosen on S so that the resulting orthonormal frame field is smooth. In terms of this frame the Fermi transport law reduces to
To restrict the gauge freedom for the coordinates \(x^{\mu }\), we assume that \(\tau \equiv x^0 = 0\) on S and the \(x^{\mu }\) be dragged along with U so that
The remaining not necessarily vanishing connection coefficients are given by
It holds \(\nabla _i\,U_k = \Gamma _i\,^0\,_k = \chi _{ab}\,\delta ^a\,_i\,\delta ^b\,_k + f_{b}\,\delta ^0\,_i\,\delta ^b\,_k\) and \(U^i\,\nabla _i\,U_k = f_{b}\,\delta ^b\,_k\), \(\nabla _i\,U^i = \chi _{a}\,^{a}\).
The first structural equations supply the constraint
and the evolution equations
The second structural equations imply for \(\Gamma _a\,^b\,_{c}\) and \(\chi _{ab}\) the constraints
and the evolution equations
No equation for \(f_a\) is implied by the structural equations.
The fluid equations (2.16), (2.17) imply in our gauge the evolution equation
and the constraint
With \({\hat{\rho }} = \Omega ^{e}\,\rho \) the latter is seen to be (2.19), which reads with \(f_a = U^{\mu }\,\nabla _{\mu }U_{\nu }\,e^{\nu }\,_a \)
Equation (2.18) is equivalent to the two equations
and
With \(\chi _{c[a}\,\chi _{b]}\,^c = - \chi _{[ac]}\,\sigma _b\,^c + \chi _{[bc]}\,\sigma _a\,^c\), where \(\sigma _{ab} =\chi _{(ab)}\), the latter can be read as a linear homogeneous ODE for \(\chi _{[ab]}\).
The system of conformal field equations reads with \(e = const.\)
with
4 Regularity of the Equations
In the following an equation will be called ‘regular’ if no negative or non-integer powers of \(\Omega \) occur in it. The structural equations are regular in this sense but the tensorial equations can be problematic. I will use from now on the notation
The asymptotic dust equation of state.
With the asymptotic dust equation of state and \(e = 3\) there enter in various places of the equations functions such as
As \(w^* \rightarrow 0\) they approach the values of the corresponding quantities of the pure dust equation of state. Therefore it is clear that equations which are singular in that case must also be singular in the present case.
In the present case the singular terms cannot be handled as in the pure dust case because equation (2.3) remains a partial differential equation and the \(\Omega ^{-1}\) terms cannot be compensated by suitable choices of k. Moreover, the expression
which comes with a singular factor in the equation
is not even defined if \(w' = 0\). Since there is no obvious way to handle it, this case will not be considered any further in this article.
The asymptotic radiation equation of state.
The asymptotic radiation equation of state with \(e = 4\) gives rise to expressions like
The limits as \(w^* \rightarrow 0\) are well defined and yield the corresponding functions in the pure radiation case. If \(\rho > 0\) and bounded as \(\Omega \rightarrow 0\) we can assume for small \(|\Omega |\) that
The structural equation are regular with no condition on k. Equations (3.15), (3.16), (3.17) are immediately seen to be regular with \(4\,k \in {\mathbb {N}}_0\). The fluid equation (3.10) reads
It is regular if \(4\,k \in {\mathbb {N}}\). As long as the term in curly brackets is continuous, the function \(\rho \) will stay positive if it is positive on some Cauchy surface. With (4.2) the constraint (3.11) can be solved for \(\nabla _{a}\,\rho \) in a neighbourhood of \(\Omega = 0\)
The equation is regular if \(4\,k \in {\mathbb {N}}\). Inspection of the equations
and
shows that they are regular if \(4\,k - 1 \in {\mathbb {N}}\). With the expressions above the equations
are seen to be regular if \(4\,k \in {\mathbb {N}}_0\). We have finally
where
Its components, given by
are regular if \(4\,k \in {\mathbb {N}}_0\).
5 Cauchy Problems
To solve the conformal equations with an asymptotic radiation equation of state we extract from the complete system a symmetric hyperbolic reduced system for the unknown
where the Ricci scalar \(R = R[g] = 6\,L_j\,^j\) of \(g_{\mu \nu }\) plays the role of a conformal gauge source function that controls implicitly the evolution of the conformal factor. It can be prescribed as an arbitrary function of the coordinates. The present discussion is concerned with the transition from the ‘physical part of the solution’ to scri and beyond. In physical terms this involves a discussion of a domain of infinite temporal extent but in conformal terms it involves only a finite step in the conformal time \(\tau \). The specification of R[g] is thus not particularly critical. It is chosen to be constant.
Symmetric hyperbolic reduced equations.
The reduced system is essentially built from the equations which contain derivatives in the direction of U. However, sometimes these equations are modified by using constraints, i.e. equations which do not contain derivatives in the direction of U. For the first six components of Z we get the system
where (3.10) resp. (4.3) has been used to replace \(\nabla _0\rho \),
Adding a suitable contraction of (3.7) to equation (3.13) gives the equation
If equations (5.1) and (5.2) are written as a system for \(f_a\) and \(\chi _{bc}\) with principal part
they represent, with given right hand sides, a symmetric hyperbolic system.
To derive the reduced systems for the curvature fields we follow the methods of [5] and earlier articles. Equation (4.10) provides a system for the unknowns \(L_{0a} = L_{a0}\) and \(L_{ab} = L_{ba}\) that contains no derivatives of \(L_{00}\). Where this quantity appears in the equation the contraction \(- L_{00} + g^{ab}\,L_{ab} = L_j\,^j = R/6\) can be used to express it in terms of \(L_{ab}\) and the gauge source function \(R = R[g] = const\). The system is given by
To extract the desired system from (4.11) consider the fields
The symmetric, trace-free U-electric part \(w_{jl}\) and the U-magnetic part \(w^*_{jl}\) of \(W^{i}\,_{j k l}\),
allow us to represent the rescaled conformal Weyl tensor in form
Suitable evolution equations for \(w_{jl}\) and \(w^*_{jl}\) is given by \(\nabla ^i\,W_{i(a|dc|}\,\epsilon _{b)}\,^{dc} = K_{dc(a}\,\epsilon _{b)}\,^{dc}\), \(\nabla ^i\,W_{i(a|0|b)} = K_{0(ab)}\) which read in detail
and
with ‘spatial derivatives’ \(D_a\,w^*_{bc} \equiv e_a\,(w^*_{bc}) - \Gamma _a\,^e\,_b\,w^*_{ec} - \Gamma _a\,^e\,_c\,w^*_{be} = \nabla _a\,w^*_{bc}\). The components of \(M_{jkl}\) on the right hand sides of equations (5.3) – (5.7) are given by
where here and in (5.3), (5.4), (5.5) equations (4.3) and (4.4) must be used to replace \(\nabla _{0}\rho \) and \(\nabla _a\rho \).
Observing the symmetry of \(w_{ab}\), \(w^*_{ab}\) (but ignoring their trace-freeness, which can later be recovered as a consequence of the initial data and the equations above), equations (5.6) and (5.7) can written as system for the unknowns \(w_{ab}\) and \(w^*_{ab}\) with \(a \le b\), \(a, b = 1, 2, 3\). If the equations \(w_{ab,\,0} = \ldots \), \(w^*_{ab,\,0} = \cdots \) are then written with a factor 2 if \(a < b\), the system is seen to be manifestly symmetric hyperbolic.
Together with the previous equations we have obtained now a quasi-linear system for the unknown Z which can be written in the form \(\textbf{A}^{\mu }\,\partial _{\mu }{} \textbf{Z} = \textbf{B}\) with matrix-valued functions \(\textbf{B} = \textbf{B}(\textbf{Z})\) and \(\textbf{A}^{\mu } = \textbf{A}^{\mu }(\textbf{Z})\) that are symmetric, i.e. \(^t\textbf{A}^{\mu } = \textbf{A}^{\mu }\), with \(\textbf{A}^{0}\) defining a positive definite bilinear form if the \(e^0\,_a\) are not too large.
Cauchy data for the conformal equations.
Solutions that admit smooth conformal extensions at future time-like infinity can be constructed from data for the conformal field equations on a compact Cauchy hypersurface S in the ‘physical domain’, where \(\Omega > 0\), or on the compact 3-manifold \(S = {{\mathcal {J}}}^+ = \{\Omega = 0\}\) that represents future time-like infinity.
Unless the field U is assumed to be orthogonal to S, the data, which must satisfy the constraints induced by the conformal field equations on S, are in general expressed in terms of the unit normal to S and then transformed into a frame \(e_k\) with \(e_0 = U\). This requires calculations involving the complete system of conformal field equations which are fairly tedious (see [7] where the presence of a boundary requires this) and give little insight. The calculation will be skipped here.
We will first construct asymptotic end data on \({{\mathcal {J}}}^+\) with the assumption that U is orthogonal to \({{\mathcal {J}}}^+\). It is interesting to note that in the case of a pure dust equation of state one is forced into this requirement without forcing the field U to be hypersurface orthogonal in the physical domain (see [5]). In the present case this assumption implies a genuine restriction. The analysis of the asymptotic end data follows essentially the one given in the vacuum case [1] with some modifications if matter fields are involved [4, 5]
To analyse the constraints, i.e. the equations that do not involve derivatives in the direction of U, we assume that \(k \ge 1\). The restrictions of the equations to \({{\mathcal {J}}}^+\) are then regular and the expressions containing the quantities \(w^*\), \((w^*)'\) drop out. The unknowns are then determined as follows.
The requirement that U is orthogonal to \({{\mathcal {J}}}^+\) implies
The first structural equation (3.1) implies on \( {{\mathcal {J}}}^+ \) the first structural equation with respect to the frame \(e_a\) and metric \(h_{ab}\) induced on \({{\mathcal {J}}}^+\). The \( \Gamma _{a}\,^{c}\,_{b}\) are thus the connection coefficients of the covariant derivative operator D induced by the metric \(h_{ab}\) in the frame \(e_a\). The second structural equation (3.2) reduces to the relation \(R^a\,_{cdb}[h] = 2\,g^a\,_{[d}\,L_{b] c} + 2\,L^a\,_{ [d}\,g_{b] c}\), equivalent to
This gives \(L_{00} = g^{ab}\,L_{ab} - L_j\,^j = R[h]/4 - R/6 \) where \(R = R[g]\) represents the conformal gauge source function. Equations (3.11), (3.15), (3.16), (3.17) imply
where it is assumed that U is future directed and \(\Omega \) is decreasing in the direction of U. By a rescaling \(g_{\mu \nu } \rightarrow \theta ^2\,g_{\mu \nu }\) and \(\Omega \rightarrow \theta \,\Omega \), where \(\theta > 0\) is smooth with prescribed value and suitable normal derivative on \({{\mathcal {J}}}^+\), it can be achieved that
Equation (3.18), which reduces to
were \(l_{ab} = l_{ab}[h]\) denotes here the Schouten tensor of h, relates the Cotton tensor of h, given on the left hand side, to the U-magnetic part of \(W^i\,_{jkl}\). The constraints induced by (3.19), given by \(\nabla _{l}\,W^{l}\,_{0 i j} = M_{ij0}\), translate with \(M_{ab0} = 0\) and \(M_{0b0} = 0\) to the constraint
and to \(D^a\,w^*_{ab} = 0\), which is not a constraint but the identity satisfied by the Cotton tensor.
If the data
are prescribed on \({{\mathcal {J}}}^+\) so that \(h^{\alpha \beta } =\delta ^{ab}\,e^{\alpha }\,_a\,e^{\beta }\,_b\) defines the contravariant version of a Riemannian metric, \(\rho > 0\), and \(w_{ab} = w_{ba}\), \(w_c\,^c = 0\) with \(w_{ab}\) satisfying (5.8) with the covariant derivative operator D defined by h, all the unknowns subsumed by \(\textbf{Z}\) can be determined in our gauge.
On the preservation of the constraints.
Smooth data on \({{\mathcal {J}}}^+\) as described above, determine a smooth solution to the symmetric hyperbolic reduced equations that cover a neighbourhood of \({{\mathcal {J}}}^+ = \{\tau = 0\}\) with \(\Omega > 0\) in the past of \({{\mathcal {J}}}^+\), where the parameter \(\tau \) on the flow lines of U is negative, and \(\Omega < 0\) in the future of \({{\mathcal {J}}}^+\) where \(\tau > 0\) [10]. We can assume it to exist in a range \(- \tau _* \le \tau \le \tau _*\) with some \(\tau _* > 0\) and \(\tau = 0\) on \({{\mathcal {J}}}^+\).
Once this solution has been obtained, it remains to show that the solution to the reduced equations does in fact also satisfy the constraints and thus the complete system of conformal equations. Since the proof follows a standard recipe which has been worked out in detail in earlier articles (see [5], and in particular [3, 8] for a discussion of the Einstein-perfect-fluid equations with a general equation of state), we skip this step.
There remains, however, an open question. Since we assumed U to be orthogonal to \({{\mathcal {J}}}^+\) we should expect U to be hypersurface orthogonal, that is \(\chi _{[ab]} = 0\) in the range \(- \tau _* \le \tau \le \tau _*\). Because \(\chi _{[ab]} = 0\) on \({{\mathcal {J}}}^+\), one could think this to be a consequence of equation (4.6). Since this equation is neither a constraint nor satisfied as part of the reduced system, we use a different argument. The evolution of \(\chi _{ab}\) is governed by equation (5.2), which implies
with \(\sigma _{ab} = \chi _{(ac)}\).We formally write here \(D_a f_{b} = e_a(f_{b}) - f_d\,\Gamma _a\,^d\,_b\), not implying that D denotes a covariant derivative away from \({{\mathcal {J}}}^+\). Among the constraints which we can assume to be satisfied there is the relation (3.11) or its more concise version (3.12), i.e.
Since the connection \(\nabla \) is torsion free it follows
where the function
is regular if \(k \ge 1\). The resulting ODE for \(\chi _{[ab]}\) implies that \(\chi _{[ab]} = 0\), equation (3.14) is satisfied, and U is hypersurface orthogonal in the range \(- \tau _* \le \tau \le \tau _*\).
With this we discussed all the ingredients needed to obtain existence results. General results in symmetric hyperbolic systems [10] allow us to draw the following conclusions.
Let the reduced conformal Einstein-\(\lambda \)-perfect-fluid equations with asymptotic radiation equation of state where \(k \ge 1\)be given in terms of the gauge discussed above so that \(e_0 = U\).Assume the future directed flow vector U to be orthogonal to the compact 3-manifold \({{\mathcal {J}}}^+ = \{\Omega = 0\}\).Let on \({{\mathcal {J}}}^+\) be given a minimal end data set consisting of a smooth frame field \(e_a = e^{\alpha }\,_a\,\partial _{x^{\alpha }}\),a positive function \(\rho \),and a symmetric trace-free tensor field \(w_{ab}\)that satisfies \(D^a\,w_{ab} = 0\)where D is the covariant derivative operator associated with the metric h that satisfies \(h(e_a, e_b) = \delta _{ab}\). Then:
− As discussed above, a complete set of Cauchy data for the conformal field equations can be calculated from the minimal set if the gauge conditions, the constraints, and the special features of \({{\mathcal {J}}}^+ = \{\Omega = 0\}\) are taken into account.
− These data determine a smooth solution \(\Omega , \,g, \, \ldots ,\,\) to the conformal field equation (unique up to extensions) so that \(\Omega < 0\)in the future of \({{\mathcal {J}}}^+\), \(\Omega > 0\) in the past of \({{\mathcal {J}}}^+\), and \({\hat{g}} = \Omega ^{-2}\,g\), \({\hat{U}} = \Omega \,U\), \({\hat{\rho }} = \Omega ^4\,\rho \) satisfy where \(\Omega \ne 0\) the Einstein-\(\lambda \)-perfect-fluid equations with asymptotic radiation equation of state. The flow field \({\hat{U}}\) is hypersurface orthogonal.
− Let S be a Cauchy hypersurface in the ‘physical domain’ of this solution where \(\Omega > 0\) and denote by \(\Delta \) the Cauchy data induced by the solution on S.Let \(\Delta '\) be Cauchy data on S for the same system of equations. If these data are sufficiently close to \(\Delta \) they develop into a solution which is time-like and null geodesically future complete, admits a smooth conformal boundary \({{\mathcal {J}}}'^{+}\) representing its future time-like infinity and a smooth conformal extension beyond with \(\rho > 0\).
− If the value \(w^* = 0\) is admitted, the resulting set of solutions, which contain as special cases the FLRM-pure-radition solutions, develop a smooth conformal boundary in the future.
Given the data \(\Delta '\) whence the corresponding flow field \(U'\) on S, the field equations allow us to determine the quantity \(U'_{[i}\,\nabla '_j U'_{k]}\) on S. In the last statement it is not required that this quantity vanishes on S. Thus the flow field \(U'\) determined by the data \(\Delta '\) will not necessarily be hypersurface orthogonal and \(U'\) need not be orthogonal to \({{\mathcal {J}}}'^{+}\).
If Cauchy data are given on \({{\mathcal {J}}}^+\) with a flow vector field U that is not necessarily orthogonal to \({{\mathcal {J}}}^+\) the result stated at the end of the introduction is obtained.
References
Friedrich, H.: Existence and structure of past asymptotically simple solution of Einstein’s field equations with positive cosmological constant. J. Geom. Phys. 3, 101–117 (1986)
Friedrich, H.: On the existence of n-geodesically complete or future complete solutions of Einstein’s field equations with smooth asymptotic structure. Commun. Math. Phys. 107, 587–609 (1986)
Friedrich, H.: Evolution equations for gravitating ideal fluid bodies in general relativity. Phys. Rev. D 57, 2317–2322 (1998)
Friedrich, H.: Smooth non-zero rest-mass evolution across time-like infinity. Ann. Henri Poincaré 16, 2215–2238 (2015)
Friedrich, H.: Sharp asymptotics for Einstein-\(\lambda \)-dust flows. Commun. Math. Phys. (2016). https://doi.org/10.1007/s00220-016-2716-6
Friedrich, H.: Cosmological models with positive \(\lambda \): hierarchies of future behaviour. Phil. Trans. R. Soc. A 382, 20230044 (2023)
Friedrich, H., Nagy, G.: The initial boundary value problem for Einstein’s vacuum field equations. Comm. Math. Phys. 201, 619–655 (1999)
Friedrich, H., Rendall, A.: The Cauchy Problem for the Einstein Equations. In: Schmidt, B. (ed.) Einstein’s Field Equations and Their Physical Implications. Springer, Lecture Notes in Physics, Berlin (2000)
Had\({\check{z}}\)ić, H., Speck, J.: The global future stability of the FLRW solutions to the Dust–Einstein system with a positive cosmological constant. J. Hyperbolic Differ. Equ. 12, 87 (2015)
Kato, T.: The Cauchy problem for quasi-linear symmetric hyperbolic systems. Arch. Ration. Mech. Anal. 58, 181–205 (1975)
Liu, C., Wei, C.: Future stability of the FLRW space-time for a large class of perfect fluids. Ann. Henri Poincaré 22, 715–770 (2021)
Lübbe, C., Valiente Kroon, J.A.: A conformal approach to the analysis of the non-linear stability of radiation cosmologies. Ann. Phys. 328, 1–25 (2013)
Olyinyk, T.A.: Future stability of the FLRW fluid solutions in the presence of a positive cosmological constant. Commun. Math. Phys. 346, 293–312 (2016)
Reula, O.: Exponential decay for small non-linear perturbations of expanding flat homogeneous cosmologies. Phys. Rev. D 60, 083507 (1999)
Rodnianski, I., Speck, J.: The nonlinear future stability of the FLRW family of solutions to the irrotational Euler–Einstein system with a positive cosmological constant. J. Eur. Math. Soc. 15, 2369–2462 (2013)
Speck, J.: The nonlinear future stability of the FLRW family of solutions to the Euler–Einstein system with a positive cosmological constant. Sel. Math. New Ser. 18, 633–715 (2012)
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Mihalis Dafermos.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Friedrich, H. Cosmological Einstein-\(\lambda \)-Perfect-Fluid Solutions with Asymptotic Dust or Asymptotic Radiation Equations of State. Ann. Henri Poincaré 25, 4263–4282 (2024). https://doi.org/10.1007/s00023-024-01416-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00023-024-01416-3