
==== Front
Arch Ration Mech Anal
Arch Ration Mech Anal
Archive for Rational Mechanics and Analysis
0003-9527
1432-0673
Springer Berlin Heidelberg Berlin/Heidelberg

2030
10.1007/s00205-024-02030-7
Article
Slowly Expanding Stable Dust Spacetimes
Fajman David 1
Ofner Maximilian 1
http://orcid.org/0000-0001-5120-1839
Wyatt Zoe zoe.wyatt@maths.cam.ac.uk

2
1 https://ror.org/03prydq77 grid.10420.37 0000 0001 2286 1424 Faculty of Physics, University of Vienna, Boltzmanngasse 5, 1090 Vienna, Austria
2 Department of Pure Mathematics and Mathematical Statistics, Wilberforce Road, Cambridge, CB3 0WB UK
Communicated by M. Dafermos.

13 9 2024
13 9 2024
2024
248 5 8313 7 2021
21 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
We establish the future nonlinear stability of a large class of FLRW models as solutions to the Einstein-Dust system. We consider the case of a vanishing cosmological constant, which, in particular implies that the expansion rate of the respective models is linear, i.e. has zero acceleration. The resulting spacetimes are future globally regular. These solutions constitute the first generic class of future regular Einstein-Dust spacetimes not undergoing accelerated expansion and are thereby the slowest expanding generic family of future complete Einstein-Dust spacetimes currently known.

Mathematics Subject Classification

35Q75
83C05
35B35
http://dx.doi.org/10.13039/501100002428 Austrian Science Fund P 34313-N Y963 Fajman David issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

General Relativistic Hydrodynamics

The Einstein-relativistic Euler system (EES)1.1 Rμν-12Rgμν=Tμν∇μTμν=0Tμν=(ρ+p)uμuν+pgμν

describes the dynamical evolution of a four-dimensional spacetime (M,g) containing a relativistic perfect fluid with pressure p, energy density ρ and 4-velocity vector uμ. Perfect fluids compatible with relativity were one of the earliest matter models considered in general relativity [7, 22] and have been extensively studied in the context of general relativistic hydrodynamics with numerous applications ranging from astrophysics to cosmological evolution (see e.g. [10, 34]).

The Eq. (1.1) are supplemented by specifying an equation of state which relates the energy density and pressure p=f(ρ). Different choices of the function f encode different behaviour of the fluid. We focus in the following on the class of linear, barotropic equations of state, p=cS2ρ, where the constant cS denotes the speed of sound of the fluid with 0≤cS≤1. This equation of state contains well-known fluid models: for cS=0, i.e. p=0, (1.1) reduces to the Einstein-Dust system, the case cS=1/3 is the Einstein-radiation fluid system and cS=1 is the Einstein-stiff fluid system. The main result of the present paper can be roughly stated as follows:

Theorem 1.1

All four-dimensional FLRW spacetime models with compact spatial slices and negative spatial Einstein geometry are future stable solutions of the Einstein-Dust system.

To date, all known future stability results establishing the global existence, regularity and completeness of solutions to (1.1) concern the regime of accelerated expansion. In such a setting, the fast decay rates of perturbations induced by the expansion have a strong regularization effect on the fluid. For slower expansion rates this effect becomes weaker and global regularity of solutions is less likely to hold.

Theorem 1.1 establishes the first nonlinear future stability result for a coupled Einstein-relativistic Euler system in the absence of accelerated expansion and thereby initiates the study of the EES in the regime of non-accelerated expansion. Such a regime is also relevant in cosmology. The epoch in the early universe, shortly after a hypothetical inflationary phase, is expected to not initially have exhibited accelerated expansion. It is this epoch, which is not covered by previous results on the EES, which we intend to make accessible by the research initiated in the present paper.

Background and Previous Results

For the sake of the following presentation we consider four-dimensional FLRW spacetimes of the form1.2 (0,∞)×M,-dtc2+a(tc)2·γ,

where (M,γ) is a complete Riemannian manifold. We remind the reader that in the standard FLRW-models the spatial slices appearing in (1.2) have constant sectional curvature kM∈{-1,0,1} and so are taken to be one of R3,S3,H3 or quotients thereof (eg, T3). We distinguish three classes of scale factors: a¨(tc)>0 are referred to as accelerated expansion, a¨(tc)<0 as deccelerated expansion and a¨(tc)=0 as linear expansion. We introduce the notion of power law inflation, where a(tc)=(tc)p for p>0. Finally, we recall that a cosmological constant can be included by adding +Λgμν to the LHS of (1.1). In the following discussion only Λ≥0 is relevant.

Shock Formation

Relavistic and non-relativistic fluids are well-known to form shocks in finite time. This was first observed in the general relativistic context by Oppenheimer and Snyder when they investigated the collapse of spherically symmetric clouds of dust [22]. In the terminology of current stability analysis this constitutes the instability of Minkowski spacetime as a particular solution to the Einstein-Dust system. Note that a part of Minkowski spacetime corresponds to M=R3,γ=δ,a(tc)=1 in (1.2).

More recently, Christodoulou’s monograph [9] demonstrated that under a very general equation of state, the constant solutions to the relativistic Euler equations on a fixed background Minkowski space are unstable, i.e. even without gravitational backreaction, fluids form shocks from arbitrarily small initial inhomogeneities in finite time. This suggests that Minkowski spacetime is unstable as a solution to the EES for a large class of equations of state. Note also that Christodoulou’s monograph gave a detailed description of the nature of the fluid shock formation, thus providing a major extension beyond work of Sideris [29] on the non-relativistic Euler equations.

Λ-induced Accelerated Expansion and Stabilisation of Fluids

As is clear from the previous paragraphs, a powerful dispersive mechanism is required to regularise fluids and to prevent finite-time shock formation. The prime example of such a mechanism comes from cosmological models exhibiting exponential expansion. Heuristically speaking, a cosmological constant Λ>0 generates expansion of the form a(tc)∼eHtc where H=Λ/3, for all cases of the sectional curvature kM. The cosmological constant creates damping terms in the equations of motion for the fluid, which dilutes the fluid and causes fluid lines to ‘stretch apart’, thus preventing shock formation.

This effect was first observed by Brauer, Rendall and Reula [6] who studied Newtonian cosmological models with Λ>0 and a perfect fluid (albeit for a slightly different equation of state). They found that the regularising effect from the exponential expansion was strong enough to prevent shock formation for small inhomogeneities of initially uniformly quiet fluid states. See also the late-time asymptotics work by Reula [25] and Rendall [23].

Moving to the fully coupled Einstein-relativistic Euler system, there has been much research concerning spacetimes undergoing exponential expansion. The first result is by Rodnianski and Speck [28], who proved future stability to irrotational perturbations of uniformly quiet fluids with 0<cS<1/3 on FLRW-spacetimes with underlying spatial manifold M=T3. The irrotational restriction was later removed by Speck in [30]. An alternative proof of future stability for these FLRW-background solutions was later given by Oliynyk [20], whose Fuchsian techniques were able to uniformly cover the cases 0<cS≤1/3.

Moving to the case of dust cS=0, stability for FLRW-spacetimes with underlying spatial manifold M=T3 was given by Hadžić and Speck [16], while more general spatial manifolds were considered by Friedrich [15] using his conformal method. Indeed the work by Lübbe and Valiente-Kroon [19] treated the radiation case with cS=1/3 using an extension of Friedrich’s conformal method. Finally, we note that very recent work of Oliynyk [21] has established stability for ultra-radiation fluids 1/3<cS≤1/2 on a fixed, exponentially expanding spacetime.

Alternative Mechanisms for Accelerated Expansion

A cosmological constant is not the only known mechanism for generating solutions to Einstein’s equations with accelerated expansion. In work that predates the references of Sect. 1.2.2, Ringström [26] considered the future global stability of a large class of solutions to the Einstein-nonlinear-scalar field system with a scalar field potential V(Φ) that satisfied V(0)>0,V′(0)=0,V′′(0)>0. Roughly speaking V(0) emulates the cosmological constant Λ and so these spacetimes undergo accelerated expansion. Ringström [27] later considered alternative potentials V(Φ) which relaxed the rate of spacetime expansion to the class of accelerated power law inflation, which in our terminology corresponds to p>1.

Note that in Ringström’s papers the global spatial topology becomes irrelevant for the long time behaviour of cosmological spacetimes in the small data regime. This is in sharp contrast to the Einstein vacuum equations where the spatial topology does affect the long-time behaviour. In this case, only the Milne geometry with a negative spatial curvature yields future eternally expanding cosmological models with precisely linear expansion rate.

Finally we note that the Chaplygin equation of state, which describes a fluid with negative pressure, can also generate sufficient spatial expansion to ensure future stability results for the coupled EES system [18].

Critical Expansion Rates

Interpolating between Minkowski space (which can be considered as a cosmological spacetime with non-compact slices and no expansion) and exponentially expanding spacetimes, it is clear that there must be a transition between shock formation and stability. To investigate the expansion rate for which this transition occurs, and how it depends on the equation of state, it is useful to study the stabilisation of fluids on fixed Lorentzian geometries obeying power-law inflation. We consider M=T3 with a(tc)=(tc)p for p>0. The following table summarises some of the main results concerning linear equation of states p=cS2ρ from [13, 31], (see also [33]): Case	Power-law rate	Range of cS	Behaviour	References	
No. 1	p>1	0<cS<1/3	Stable	[31]	
No. 2	p=1	cS=1/3	Shocks	[31]	
No. 3	p=1	0<cS<1/3	Stable (irrot.)	[13]	
No. 4	p>12	cS=0	Stable	[31]	

In combination, these results indicate that in spacetimes undergoing power-law inflation whether shocks form from small data depends on the equation of state and, in particular in the linear case, on the speed of sound. Indeed the literature suggests that slower speeds of sound reduce the tendency of shock formation. For the particular case of dust (cS=0), case No. 4 shows that shocks are avoided even in deccelerating spacetimes with scale factors a(t)=t1/2+δ for δ>0.

Main Results

In the present paper we consider the Einstein-Dust system in the regime of linear expansion. For the linearly expanding case, Cases No. 2 and 3 above show that even in the absence of backreaction the speed of sound determines whether shocks form or not. We prove that for the case of dust (cS=0) shock formation does not occur under the full gravity-fluid dynamics.

Our background geometry is that of the Milne model, which generalises the kM=-1 FLRW vacuum spacetimes. Let (M,γ) be a closed, connected, orientable three-dimensional manifold admitting a Riemannian Einstein metric γ with negative Einstein constant. After rescaling, we suppose thatRic[γ]=-29γ.

The generalised Milne spacetime is the Lorentz cone spacetime M=(0,∞)×M with metricgM:=-dtc2+tc29γabdxadxb.

The spacetime (M,gM) is globally hyperbolic and a solution to the four-dimensional vacuum Einstein equations. We formulate the main theorem using terminology introduced in Sect. 3.1. We let Bεj,k,l,m(t029γ,-t09γ,0,0) denote the ball of radius ε in the space Hj×Hk×Hl×Hm centred at (t029γ,-t09γ,0,0). Our main theorem is

Theorem 1.2

Let (M,gM) be as above. Let ε>0 and (g0,k0,ρ0,u0) be initial data for the Einstein-Dust system at tc=t0 such that(g0,k0,ρ0,u0)∈Bε6,5,4,5t029γ,-t09γ,0,0.

Then, for ε sufficiently small the corresponding future development under the Einstein-Dust system is future complete and admits a CMC foliation labelled by τ∈[τ0,0) such that the induced metric and second fundamental form on constant CMC slices converge as(τ2g,τk)→γ,13γasτ↗0i.e. astc↗∞.

If the initial energy density of the dust field is non-negative, ρ0≥0, then it remains so throughout the evolution.

Remark 1.3

The Milne model is known to be a stable solution to the Einstein vacuum equations [3], the Einstein massive-Vlasov equations [1], the coupled Einstein-Maxwell-scalar field system arising from a Kaluza-Klein reduction [5], and the Einstein Klein-Gordon equations [14, 32].

Remark 1.4

Negative spatial curvature is crucial as spherical or toroidal spatial topologies would lead to recollapsing or slowly expanding matter dominated solutions, respectively. The asymptotic behaviour of the solutions in the theorem coincide with the corresponding vacuum solutions.

Structure and Key Novelties in the Proof

The proof of Theorem 1.2 consists of three major parts: (i) energy estimates for the perturbation of the spacetime geometry with sources given by the dust variables, (ii) energy estimates for the dust variables in the perturbed spacetime geometry and (iii) a bootstrap argument based on both sets of energy estimates establishing global existence and asymptotic behaviour. This rough approach is standard in the literature on the Milne stability problem (see e.g. [1, 3, 5]), however, for the Einstein–dust system there are crucial difficulties caused by a regularity problem inherent to the dust equations, which turns out to affect all parts of the argument. We outline the difficulties and how these are overcome in the following.

To control the perturbed spacetime geometry throughout the evolution we use a CMC time-foliation in combination with a spatial-harmonic gauge [2]. The existence of such a foliation for small perturbations of negative Einstein spaces is non-trivial but standard [12]. The Einstein equations then take the form of an elliptic-hyperbolic system (see (2.7)) where the lapse and shift are determined by elliptic PDEs with sources given in terms of metric, second fundamental form and the dust variables. This elliptic system provides Sobolev estimates for the lapse and shift.

The core idea to establish decay for the geometric variables in previous works on Milne stability is a corrected energy (Egk in Definition 3.8) based on the modified Einstein-operator (Lg,γ in Definition 3.2) of the spatial Einstein geometry [1, 3]. For the Einstein–Dust system we must deviate from this standard approach due to a regularity issue from the dust model, which in turn affects all parts of the proof.

When expanded, the equations of motion for the dust variables take a form where the source term of the evolution equation for the energy density contains the spatial divergence of the fluid velocity (see (2.7d)). Consequently, the fluid energy density can be controlled only in one order of regularity below the order of regularity of the fluid velocity. From the perspective of the Einstein equations this is very problematic as both components of the dust, energy density and fluid velocity, appear at the same order of regularity as source terms of the Einstein equations. As such, they are required to be controlled in suitable Sobolev spaces at the same order as the second fundamental form. Due to the required high regularity of the fluid velocity discussed previously, the velocity then needs to be controlled one order above the second fundamental form. However, the equation of motion for the fluid velocity requires the second fundamental form at the same order of regularity as the velocity itself (see (2.7d)). This apparent inconsistency prevents one from establishing a standard and straightforward regularity hierarchy to analyse the fully coupled nonlinear system.

An approach to circumvent this issue has been introduced by Hadžić and Speck in [16] and is modified in the present paper. The central idea is to use a fluid derivative ∂u∼uα∇α[gM] as a differential operator in the energies for the perturbations of the metric and second fundemental form (E∂u,N-1g in Definition 3.11). At highest order of regularity, say N, where the loss of derivatives prevents the closure of the system of estimates, the Einstein equations are commuted with N-1 spatial derivatives and one fluid derivative. When this derivative acts on the dust source terms in the Einstein equations, in the subsequent calculations for the energy estimates, the equations of motion of the dust variables are used as constraint equations replacing ∂uρ and ∂uuj. In this way, no derivatives are lost and the corresponding auxiliary energies E∂u,N-1g for the geometric variables can be estimated in terms of dust variables of one order of regularity below the expected one.

In a follow-up step, we need to show that the auxiliary geometric energies E∂u,N-1g in fact control the actual top-order regularity norms of the geometric variables. This is achieved by rewriting the wave-type evolution equation for the metric and second fundamental in terms of an elliptic part and certain mixed spatial and fluid derivative operators (see Proposition 6.3). Consequently, the auxiliary energies of the first step provide top-order estimates on the geometric variables (see Corollary 6.4) and an overall strategy to close the estimates.

Two final major regularity issues arise when proving energy estimates for the auxiliary energies E∂u,N-1g however. We end up needing to estimate one fluid derivative and a critical number of spatial derivatives on certain terms involving the lapse and shift, and we cannot commute the ∂u operator past the spatial derivatives without exceeding the assumed regularity of the fluid spatial velocity.

To circumvent this problem, we only commute past some of the derivatives and instead derive two auxiliary estimates using the elliptic Eq. (2.7b) for the lapse and shift. In the estimate on the lapse term (see Proposition 7.3) we crucially use the equations of motion of the dust variables to replace a certain matter term ∂uη as a constraint, thus avoiding derivative loss. The estimate for the shift term (see Proposition 7.6) proceeds differently, relying on a remarkable combination of commutator estimates, the Bianchi identity and the Einstein equations in the CMCSH gauge.

Final Remarks

In the regime of non-accelerated expansion, the work [6] indicates that the backreaction between the fluid and the geometry cannot be ignored. The authors consider the case of dust with a Newtonian backreaction, finding that shocks form for arbitrarily small initial data in the regime where the homogeneous background spacetime, which is perturbed, expands like a(t)=t2/3. This contrasts noteably with case No. 5 above which does not include backreaction. While [6] concerns only Newtonian dynamics, it is nevertheless a fair indication that the fully coupled dynamics under the Einstein-fluid system will likely lead to the formation of shocks. Continuing this line of reasoning, we note that although the work [13] also treated linear expansion, the full coupling between gravity and fluid makes our present work highly nontrivial. Indeed the issues highlighted on the previous Sect. 1.3.1 are indicative of the substantial technical difficulties that arise in the fully coupled EES.

Finally, it is interesting to recall that Sachs and Wolfe derived a linear instability result for the Einstein-Dust equations with Λ=0, however their metric had underlying spatial manifold M=R3 [35]. The fluid plays a major dynamical role in these flat FLRW models. Nevertheless one gleans the importance of the negatively curved spatial slices appearing in our nonlinear stability result.

Outline of the Paper

In Sect. 2 we introduce the system of equations and perform a natural rescaling of the variables. In Sect. 3 we introduce function spaces and energy functionals controlling the metric perturbation and shear tensor.

The main theorem is proved using continuous induction. In Sect. 4 we discuss the local existence theory and initiate the bootstrap argument. The remainder of the paper, beginning with Sect. 5, treats the individual estimates necessary to close the bootstrap argument. Section 5 gathers various auxiliary estimates, which are used in later sections. Among those are estimates on the source terms of the evolution equations, estimates on the dust-derivative acting on various quantities, commutators of the dust derivative and other operators and estimates on high derivatives combining the dust-derivative and other operators.

Section 6 derives the elliptic estimate for the Einstein operator and the evolution equations, which is crucial to turn estimates in terms of the dust derivatives into those in terms of standard energies. These estimates are then given subsequently. In Sect. 7 we provide the estimates on lapse function and shift vector field. A crucial set of lapse and shift estimates on highest order of regularity, involving also the dust derivative, are given here too. In Sect. 8 we derive the central top-order energy estimate for the auxiliary energy controlling the geometric perturbations. In Sect. 9 we derive the estimates for the dust variables and in Sect. 10 we close the bootstrap.

Equations of Motion

The Einstein–Dust System

The Einstein–relativistic Euler system reads2.1 Rμν[g¯]-12R[g¯]g¯μν=2T~μν,∇¯μT~μν=0,T~μν=(ρ~+P~)u~μu~ν+P~g¯μν,

where we set c=1 and 4πG=1. We use ∇¯ to denote the Levi-Civita connection of the physical metric g¯. The four-velocity of the fluid u~μ is a future-directed timelike vectorfield normalised by2.2 g¯μνu~μu~ν=-1.

We assume a linear, barytropic fluid equation of state P~=cS2ρ~ where cS≥0 is a constant, and P~≥0 and ρ~≥0 denote the pressure and energy density respectively. In the present paper, we restrict ourselves to dust, which means we setcS2:=0.

The fluid equations in (2.1) can equivalently (for ρ~>0) be written as2.3 u~α∇¯αlnρ~+∇¯αu~α=0,u~α∇¯αu~μ=0.

The system (2.3) is overdetermined in the sense that u~0 can be determined from the other fluid velocity components via (2.2).

We will study the Einstein-Dust equations using the following ADM ansatz for the metric2.4 g¯=-N~2dt2+g~ab(dxa+X~adt)(dxb+X~bdt).

Note that g¯ab=g~ab but in general g¯ab≠g~ab. On t=constant slices, we let τ be the trace of the second fundamental form k~ with respect to g~ and define Σ to be the trace-free part of k~; that is,τ:=trg~k~=g~abk~ab,k~:=Σ~+13τg~.

We use Roman letters (a, b, i, j...) to denote spatial indices. Let ∇ denote the Levi-Civita connection of the spatial metric g~. Using (2.4) the Christoffel symbols of the 4-metric g¯ become (see e.g. [24])(4)Γ~000=N~-1(∂tN~+X~a∇aN~-k~abX~aX~b),(4)Γ~ab0=-N~-1k~ab,(4)Γ~a00=N~-1(∇aN~-k~abX~b),(4)Γ~bca=Γbca[g~]+N~-1k~bcX~a,(4)Γ~0ba=-N~k~ba+∇bX~a-N~-1X~a∇bN~+N~-1k~bcX~cX~a,(4)Γ~00a=∂tX~a+X~b∇bX~a-2N~k~caX~c+N~∇aN~-N~-1(∂tN~+X~b∇bN~-k~bcX~bX~c)X~a.

Noting the above, the fluid Eq. (2.3) reduce tou~α∂αlnρ~+∂αu~α+(4)Γ~αναu~ν=0,u~α∂αu~μ+u~α(4)Γ~ανμu~ν=0.

The Rescaled Einstein–Dust System in CMCSH Gauge

Following the work of Andersson and Moncrief [2, 3], we hereon impose the CMCSH gauge which foliates by surfaces of constant mean curvature, taking advantage of the fact that on the Milne background τ=g~abk~ab=-3/tc.

Definition 2.1

(CMCSH gauge)t=τ,Ha:=g~cb(Γ[g~]cba-Γ[γ]cba)=0.

We next rescale our variables with respect to the mean curvature τ.

Definition 2.2

(Rescaled variables (gab,N,Xa,Σab,ua,u0,ρ,N^,u^0) and logarithmic time T) The rescaled geometric variables are defined as 2.5a gab:=τ2g~ab,gab:=(τ2)-1g~ab,N:=τ2N~,Xa:=τX~a,Σab:=τΣ~ab.

Note X~a=g~abX~b=τ-3Xa. Let u~0:=u~τ. The rescaled matter variables are defined as2.5b ua:=τ-2u~a,ρ:=|τ|-3ρ~,u0:=τ-2u~0.

Denote N^:=N3-1 and u^0:=u0-1/3. Finally we define the logarithmic timeT:=-ln(τ/(eτ0)),

which satisfies ∂T=-τ∂τ.

The above definition means we have the following ranges τ0≤τ↗0 and 1≤T↗∞ where τ↗0 corresponds to the direction of cosmological expansion (i.e. tc↗∞).

Lemma 2.3

The normalisation condition (2.2) impliesu0=1(N2-XaXa)(τXaua+[τ2(Xaua)2+(N2-XaXa)(τ2gabuaub+1)]1/2).

Proof

Using (2.2) we have0=(-N~2+X~aX~a)(u~0)2+2X~au~au~0+g~abu~au~b+1.

This is a quadratic equation in u~0. The roots areu~0=12(-N~2+X~aX~a)(-2X~au~a±[4(X~au~a)2-4(-N~2+X~aX~a)(g~abu~au~b+1)]1/2).

Applying the rescalings from Definition 2.2 we findu~0=τ4(N2-XaXa)(τ-1Xaua∓[τ-2(Xaua)2+(N2-XaXa)(τ-2gabuaub+τ-4)]1/2).

Hence, we introduce the rescaled quantity u0=τ-2u~0 for the larger root. □

Let ∇,∇^ denote the Levi-Civita connection of the Riemannian metrics g,γ respectively. The Christoffel symbols of g¯ now become (see e.g. [1])(4)Γ~bca=Γbca[g]-ΓbcXa,(4)Γ~00a=τ-2Γa,(4)Γ~0ba=τ-1-δba+Γba,(4)Γ~000=τ-1(-2+ΓR),(4)Γ~ab0=τΓab,(4)Γ~0a0=Γa,

where we have introduced the following rescaled geometric components.

Definition 2.4

(Rescaled Christoffel components Γa,Γba,ΓR,Γab,Γ∘a)Γa:=-∂TXa-Xa-2N^Xa+Xb∇bXa-2NΣcaXc+N∇aN+(N-1∂TN-N-1Xb∇bN+N-1Σbc+13gbcXbXc)Xa,Γba:=-NΣba-δbaN^+∇bXa-N-1Xa∇bN+N-1Σbc+13gbcXcXa,ΓR:=N-1(-∂TN+Xa∇aN-(Σab+13gab)XaXb),Γ∘a:=Γa-N∇aN,Γab=-N-1(Σab+13gab),Γa:=N-1(∇aN-(Σab+13gab)Xb).

Remark 2.5

(Background solutions) The rescaled background (B) Milne geometry written in CMCSH gauge is(gab,Σab,N,Xa)|B≡(γ,0,3,0).

Furthermore,(Γa,Γba,ΓR,Γa)|B≡0,Γab|B≡-19γab.

Let ρ0′>0 be a constant. The background, uniformly quiet fluid solution in CMCSH gauge is(u0,ui,ρ)|B=(13,0,ρ0′);

see also Appendix A.

We next evaluate certain energy momentum and matter source terms arising from the dust.

Definition 2.6

(Matter source terms E,ȷa,η,Sab,T_ab)E:=ρ(u0)2N2,ȷa:=ρNu0ua,η:=E+ρgab(u0Xa+τua)(u0Xb+τub),Sab:=ρ(u0Xa+τua)(u0Xb+τub)+12ρgab,T_ab:=ρuaub.

For further details on these definitions see Appendix 11.

Definition 2.7

(Matter source terms Fuj,Fu0,Fρ) 2.6 Fuj:=τ-1(u0)2Γ∘j+(Γklj[g]-Γklj[γ])ukul+2u0uiΓij+τukuiΓkiXj,Fu0:=(u0)2ΓR+2τuju0Γj+τ2ukujΓkj,Fρ:=τ∇iui+Γiiu0-τΓjuj+τΓikXiuk-τuju0∇ju0-τ2ukuju0Γkj.

Equations of Motion

Bringing together all the previous notation, as well as using the general equations presented in [1], the equations of motion for the Einstein-Dust system in CMCSH gauge are the following. We have two constraint equations. 2.7a R(g)-|Σ|g2+23=4τE,∇aΣab=2τ2ȷb,

and two elliptic equations for the lapse and shift variables,2.7b (Δ-13)N=N|Σ|g2-τη-1,ΔXa+Ric[g]baXb=2∇bNΣba-∇aN^+2Nτ2ȷa-(2NΣbc-∇bXc)(Γ[g]bca-Γ[γ]bca).

We also have evolution equations for the induced metric and trace-free part of the second fundamental form2.7c ∂Tgab=2NΣab+2N^gab-LXgab,∂TΣab=-2Σab-N(Ric[g]ab+29gab)+∇a∇bN+2NΣacΣbc-13N^gab-N^Σab-LXΣab+NτSab,

and, finally, evolution equations for the fluid components:2.7d u0∂Tuj=τua∇auj+τ-1(u0)2N∇jN+Fuj,u0∂Tu0=τua∇au0+Fu0,u0∂Tρ=τua∇aρ+ρFρ.

Using notation from [2, 3] we introduce new variables which allow us to rewrite (2.7c).

Definition 2.8

(Perturbation variables h, v, w and geometric source terms Fh,Fv) Define the variableshab:=gab-γab,vab:=6Σab,w:=N/3,

and the geometric source terms(Fh)ab:=2N^gab+hac∇^bXc+hcb∇^aXc,(Fv)ab:=∇a∇bN+2NΣacΣbc-13N^gab-N^Σab+NτSab-vac∇^bXc-vcb∇^aXc.

We start with the following identity from [2]:LXgab=Xc∇^cgab+gac∇^bXc+gcb∇^aXc.

Due to rigidity properties of negative Einstein manifolds in three spatial dimensions (see e.g. [3, §1.1]), we have ∂Tγ=0. Thus the Eq. (2.7c) reduce to2.8 ∂Thab=wvab-Xm∇^mhab+Fh,∂Tvab=-2vab-9wLg,γhab-Xc∇^cvab+6Fv.

In Sect. 8 we will also write the first equation in (2.8) as2.9 ∂Thab=wvab+2N^gab-(LXg)ab=wvab+2N^gab-gam∇bXm-gbm∇aXm.

Hereon we use the differential Eq. (2.7b), (2.7d) and (2.8) to analyse the solutions to our Einstein-Dust system.

Preliminary Definitions

In this section we present several preliminary definitions concerning Sobolev spaces, norms, elliptic estimates and energy functionals. All of this is standard except for Definitions 3.9 and 3.11 where we introduce the fluid derivative ∂u and then the energy functionals for the geometric variables involving this fluid derivative.

Function Spaces and Norms

Definition 3.1

(Lg,γ2-inner product) Let μg=detg denote the volume element on (M, g), similarly for μγ. Let V, P be (0, 2)-tensors on M. Define an inner product by⟨V,P⟩γ:=VijPklγikγjl,

and define a mixed L2-scalar product(V,P)L2(g,γ):=∫M⟨V,P⟩γμg,

with corresponding norm ‖V‖Lg,γ22:=(V,V)L2(g,γ).

The following definition follows notation first introduced in [3].

Definition 3.2

(Riem[γ]∘ and Lg,γ) Let V be a symmetric (0, 2)-tensor on M. Define the tensorial contraction(Riem[γ]∘V)ij:=Riem[γ]iajbγaa′γbb′Va′b′,

where, following the convention of [2], the Riemann tensor is defined by [∇^a,∇^i]Vb=(∇^a∇^i-∇^i∇^a)Vb:=-Riem[γ]baicVc. Define the following differential operatorsΔ^g,γVij:=(detg)-1∇^a(detg·gab∇^bVij),Lg,γVij:=-Δ^g,γVij-2(Riem[γ]∘V)ij.

By using the gauge condition (2.1), the operator Δ^g,γ can be rewritten asΔ^g,γVcd=gab∇^a∇^bVcd-Ha∇^aVcd.

The operator Lg,γ is self-adjoint with respect to the mixed L2-scalar product (see e.g. [2])3.1 (Lg,γV,P)L2(g,γ)=(V,Lg,γP)L2(g,γ).

A self-adjoint elliptic operator on a compact manifold has a discrete spectrum of eigenvalues. Using eigenvalue estimates from [17], we are led to the following result:

Proposition 3.3

(Estimates on λ0) Let (M,γ) be a negative Einstein three-manifold with Einstein constant k=-2/9. Then the smallest eigenvalue of the operator Lg,γ satisfies λ0≥1/9 and the operator also has trivial kernel ker(Lg,γ)={0}.

Definition 3.4

(Sobolev norms) Let k∈Z≥0. For κ a Riemannian metric on M, f a function and V a (1, 1)-tensor, define|∇kf|κ2:=κa1b1⋯κakbk(∇a1⋯∇akf)·(∇b1⋯∇bkf),|∇kV|κ2:=κijκklκa1b1⋯κakbk(∇a1⋯∇akVki)·(∇b1⋯∇bkVlj).

The obvious extension to (p, q)-tensors holds. We write‖V‖Hk=(∑0≤ℓ≤k∫M|∇ℓV|g2μg)1/2.

Note that under a global smallness assumption on g-γ guaranteed by the bootstrap assumptions, we have the norm equivalence ‖·‖L2≅‖·‖Lg,γ2.

Remark 3.5

When we write ‖u‖Hk we are denoting a sum only over the spatial components of the velocity vector-field uμ.

We frequently, and without comment, use the following product estimate:

Lemma 3.6

(Sobolev product estimates) If s>n/p=3/2 then‖uv‖Hs≲‖u‖Hs‖v‖Hs.

We conclude this subsection with a result concerning elliptic regularity, see e.g. [4, App. H].

Lemma 3.7

(Elliptic regularity using Lg,γ) Let V be a symmetric (0, 2)-tensor on M. There exist constants C1,C2>0 such, that for all s∈Z≥0,C1‖V‖Hk+2s≤‖Lg,γsV‖Hk≤C2‖V‖Hk+2s,

where Lg,γs denotes s-copies of Lg,γ.

Energy for the Perturbation of the Geometry

As noted in [2], in the spatially harmonic gauge (2.1) we haveRic[g]ab+29gab=12Lg,γ(g-γ)ab+Jab,

where Jab are higher-order terms (writen as Sab in [2, pg. 22]) satisfying, for k≥1,‖J‖Hk-1≤C‖g-γ‖Hk.

Following [1, 3], we define an energy for the geometric perturbation of the first and second fundamental forms by using Lg,γ. This energy will fulfill a strong decay estimate enabled by the inclusion of certain correction terms Γ(m).

Definition 3.8

(Geometric energy Egk) Let λ0 be the lowest eigenvalue of the operator Lg,γ, with lower bounds given in Proposition 3.3. We define the correction parameter α=α(λ0,δα) byα:=1λ0>1/91-δαλ0=1/9,

where δα=1-9(λ0-ε′) with 1≫ε′>0 remains a variable to be determined in the course of the argument to follow. By fixing ε′ once and for all, δα can be made suitably small when necessary. The corresponding correction constant, relevant for defining the corrected energies, is defined bycE:=1λ0>1/99(λ0-ε′)λ0=1/9.

We are now ready to define the energy for the geometric perturbation. For m,k∈Z≥1 letE(m):=12v,Lg,γm-1(v)Lg,γ2+92h,Lg,γm(h)Lg,γ2,Γ(m):=v,Lg,γm-1(h)Lg,γ2.

The energy measuring the geometric perturbation is then defined byEgk:=∑1≤m≤k(E(m)+cEΓ(m)).

We will see later that the corrected geometric energy Egk is in fact coercive over the standard Sobolev norms of the geometric variables g,Σ.

Definition 3.9

(Operators ∂^0,∂u) We define the operators on M∂^0:=∂T+LX,∂u:=u0∂T-τua∇^a.

The following identity, taken from [8], holds for some function f on M:3.2 ∂T∫Mfμg=3∫MN^fμg+∫M∂^0(f)μg.

Remark 3.10

(Regularity parameters ℓ,N) At top-order our bootstrap assumptions will involve Sobolev norms HN where N is a large integer. It is convenient to require N to be odd so that we can introduce ℓ∈Z satisfying ℓ=N-12.

We end this section with the top-order geometric energy for the geometric variables g,Σ, which crucially involves the fluid operator ∂u. This energy will also fulfill a strong decay estimate enabled by the inclusion of the correction terms.

Definition 3.11

(E∂u,2ℓg and Etot) For s∈Z≥1, defineE∂u,2sg:=92∂uLg,γs(h),∂uLg,γs(h)Lg,γ2+12∂uLg,γsv,∂uLg,γs-1vLg,γ2,Γ∂u,2sg:=∂uLg,γs-1v,∂uLg,γs(h)Lg,γ2.

The corrected ∂u-boosted geometric energy is then given byE∂u,N-1g:=∑s=1ℓ(E∂u,2sg+cEΓ∂u,2sg).

Finally we defineEtotg:=E∂u,N-1g+EN-1g.

The Bootstrap Argument

In this section we first state the local-existence theory for the Einstein-Dust system in CMCSH gauge. Then we introduce the bootstrap assumptions on our solution and give some immediate consequences of these estimates.

Local Existence

Theorem 4.1

Let N≥6. Consider CMC initial data (g0,k0,N0,X0,ρ0,u0)∈HN×HN-1×HN×HN×HN-2×HN-1 at T=T0 such that the constraints (2.7a) hold. Then there exists a unique classical solution (g,k,N,X,ρ,u) on [T0,T+) for T+>T0 to the system (2.7), which is consequently also a solution to the Einstein-Dust equations. The components have the following regularity featuresg,N,X∈C0([T0,T+),HN)∩C1([T0,T+),HN-1),k∈C0([T0,T+),HN-1)∩C1([T0,T+),HN-2),u∈C0([T0,T+],HN-1),ρ,∂uρ,∂uu∈C0([T0,T+],HN-2).

Furthermore, the time of existence and the norms of the solution depend continuously on the initial data. For the maximal time of existence T∞ we have either T∞=+∞ orlimT↗T∞sup[T0,T∞]‖g-γ‖HN+‖Σ‖HN-1+‖N-3‖HN+‖X‖HN+‖∂TN‖HN-1+‖∂TX‖HN-1+‖ρ‖HN-2+|τ|‖u‖HN-1>δ(γ),

where δ(γ) is a positive fixed constant depending only on the background metric.

Proof

The proof follows analogous to [16, Theorem 3.5] where the control on the lapse and shift are replaced by the elliptic techniques applied in [11]. The adaption of the regularity scheme to avoid the loss of derivatives from [16] to the present case is executed in detail in the global analysis discussed in the remainder of the paper. The smallness condition in the continuation criterion stems from a smallness requirement in applying the corresponding elliptic equation for the lapse. □

Remark 4.2

Since we apply the local-existence theorem only for data close to the background solution the smallness condition required to extend the solution for arbitrarily large times is automatically fulfilled when our smallness conditions hold, which we prove by a bootstrap argument. Since the smallness parameter of the continuation criterion depends only on the background geometry we can choose the smallness of the initial data, which we do accordingly without mentioning it explicitly again.

Bootstrap Assumptions

Let μ,λ be fixed positive constants with μ≪1 and λ<1. We assume that, for all T0≤T≤T′, the following bootstrap assumptions hold (2.7):4.1 ‖g-γ‖HN+‖Σ‖HN-1≤Cεe-λT,‖N-3‖HN+‖X‖HN≤Cεe-T,‖∂TN‖HN-1+‖∂TX‖HN-1≤Cεe-T,‖ρ‖HN-2≤Cε,‖u‖HN-1≤CεeμT.

there T′<T∞ is fixed. We hereon assume that (4.1) hold and do not repeat this fact. Recall also Remark 3.5 regarding the norm on u.

Definition 4.3

(Λ(T)) It is convenient to introduce the notationΛ(T):=‖N^‖HN+‖X‖HN+‖∂TN‖HN-1+‖∂TX‖HN-1+|τ|‖ρ‖HN-2+τ2‖u‖HN-12.

Note that, under the bootstrap assumptions (4.1), Λ(T)≲εe-T.

We state some immediate consequences of the bootstrap assumptions regarding u0 which we use without further comment. Using Lemma 2.3, we have4.2 ‖u^0‖HN-1≲ε,‖∇u^0‖HN-2≲‖∇N‖HN-2+τ2‖ua∇ua‖HN-2+h.o.t≲Λ(T).

Note the first estimate does not pick up any μ-loss. By the Sobolev embedding H2↪L∞, ‖u^0‖L∞=‖u0-1/3‖L∞≤1/10 and thus‖u0‖L∞≲1,‖u0‖L∞-1≲1.

The next Lemma concerning the dust matter components is indicative of the good behaviour that, as we discussed in Sect. 1.2.4, we roughly expect as the speed of sound is reduced.

Lemma 4.4

(Estimates on matter components) We have|τ|(‖E‖HN-2+‖η‖HN-2+‖S‖HN-2)≲Λ(T),|τ|2‖ȷ‖HN-2≲εe(-1+μ)TΛ(T),|τ|3‖T_‖HN-2≲Λ(T)2.

Proof

The estimates are immediate by distributing derivatives across the terms written in Definition 2.6 and using Lemma 3.6. □

Lemma 4.5

(Geometric coercivity estimate) Let s∈Z such that 1≤s≤N-1. There is a δ>0 and a constant C>0 such that for (g,Σ)∈Bδ((γ,0)) the following inequality holds‖g-γ‖Hs2+‖Σ‖Hs-12≤CEsg.

Proof

The proof of this lemma follows verbatim from [1, Lemma 19], which itself follows from [3, Lemma 7.2], and the eigenvalue estimates referred to in Proposition 3.3. □

The next Lemma is actually only used once in our entire argument. Its significance lies in the fact that it allows us to convert the quadratic derivative term (∇V)2 appearing in the H1 Sobolev norm into just one second-order derivative as V·Lg,γV, which will more naturally be controlled by E∂u,N-1g.

Lemma 4.6

Let V be a symmetric (0, 2)-tensor on M. Then,‖V‖H1≲‖V‖L2+(V,Lg,γV)Lg,γ21/2.

Proof

We integrate by parts, use the closeness between the g and γ metrics and the boundedness of the Riem[γ] components:‖V‖H12=‖V‖L22+∫Mgabgijgkl∇aVik∇bVjlμg≤‖V‖L22+|∫MgabgijgklVik∇a∇bVjlμg|≤‖V‖L22+|∫M⟨V,Lg,γV⟩γμg+2∫M⟨V,Riem[γ]∘V⟩γμg|≲‖V‖L22+(V,Lg,γV)Lg,γ2.

□

Preliminary Estimates

This section is concerned with deriving several preliminary estimates that are required for our later energy inequalities. There are estimates on commutator terms (Sects. 5.1, 5.4), estimates on matter terms (Sect. 5.3) and also estimates on geometric variables (Sects. 5.2, 5.5). We also present an integration by parts Lemma 5.3, in particular (5.6c), which later plays an important role in removing various critical terms that arise during the energy estimates.

First Commutator Estimates

In this subsection we let V be an arbitrary symmetric (0, 2)-tensor and ϕ a scalar, unless otherwise specified. We begin with the identity5.1 [∂T,∇a]Vij=(-∂TΓaic)Vcj+(-∂TΓajc)Vic.

The terms ∂TΓ[g] can be estimated (see [1, (10.12)]), for 0≤k≤N-2, by:5.2 ‖∂TΓ(g)‖Hk≲‖Σ‖Hk+1+‖X‖Hk+2+‖N^‖Hk+1.

We also have the following commutator identities:5.3 [∂T,∇a1⋯∇ak]ϕ=[∂T,∇a1]∇a2⋯∇akϕ+∇a1[∂T,∇a2]∇a3⋯∇akϕ+…+∇a1⋯∇ak-2[∂T,∇ak-1]∇akϕ+∇a1⋯∇ak-1[∂T,∇ak]ϕ,

5.4 [∇i,∇a1⋯∇ak]ϕ=[∇i,∇a1]∇a2⋯∇akϕ+∇a1[∇i,∇a2]∇a3⋯∇akϕ+…+∇a1⋯∇ak-2[∇i,∇ak-1]∇akϕ+∇a1⋯∇ak-1[∇i,∇ak]ϕ.

Note the last terms in each of (5.3) and (5.4) will in fact vanish since the metric g is torsion free.

In the next part of this subsection we state an important estimate, given in (5.5), that allows us to turn background ∇^ derivatives into dynamical ∇ ones.

Definition 5.1

(Difference tensor Υ) Recalling that ∇ and ∇^ are the Levi-Civita symbols of g and γ respectively, we define Υ a (1,2)-tensor byΥbca:=Γbca[g]-Γbca[γ].

Let V be a vector and P a one-form. Then we have∇aVi=∇^aVi+ΥjaiVj,∇aPi=∇^aPi-ΥiajPj.

We will often schematically write tensorial contractions using ∗. For example,∇aVi=∇^aVi+ΥjaiVj,becomes∇V=∇^V+Υ∗V.

In local coordinates the components of the Υ tensor are given byΥbca=-12γai∇bγci+∇cγai-∇iγbc=12γai∇bhci+∇chai-∇ihbc.

Lemma 5.2

For 0≤k≤N-2, we have‖[∇^,∇]V‖Hk≲‖V‖Hk+εe-λT‖V‖Hk+1.

Proof

First we see that from Definition 5.1, for all 0≤k≤N-1,5.5 ‖Υ‖Hk≲‖g-γ‖Hk+1.

Thus we compute,[∇^c,∇a]Vij=[∇^c,∇^a]Vij+Υcab∇^bVij-((∇^cΥaib)Vbj+(∇^cΥajb)Vib).

Using the boundedness of the Riem[γ] components, the required estimate then follows. □

We end this subsection with some very useful estimates that come from integration by parts.

Lemma 5.3

(Integration by parts) Let T be an arbitrary vectorfield and V, P arbitrary symmetric (0, 2)-tensors, on M. Then 5.6a (Lg,γV,P)L2(g,γ)=∫M⟨gab∇^aV,∇^bP⟩γμg-2∫M⟨Riem[γ]∙a∙bVab,P⟩γμg.

and5.6b |(Ta∇^aV,P)Lg,γ2|≲‖T‖H3‖V‖L2‖P‖H1,

5.6c |(Ta∇^aV,V)Lg,γ2|≲‖T‖H3‖V‖L22.

Proof

The proof of (5.6a) follows by using the gauge condition Ha=0. To show (5.6b), recall that the Jacobi identity implies that∇^adetg=12detg·gij(∇^agij).

Using this we find that(Ta∇^aV,P)Lg,γ2=∫MγijγklTa(∇^aVik)Pjldetg=-∫Mγijγkl∇^aTaPjldetgdetγVikdetγ=-(V,Ta∇^aP)Lg,γ2-((∇^aTa)V,P)Lg,γ2-12(V,(Tagbc∇^agbc)P)Lg,γ2.

Thus,|(Ta∇^aV,P)Lg,γ2|=|-(V,Ta∇^aP)Lg,γ2-((∇^aTa)V,P)Lg,γ2-12(V,(Tagbc∇^agbc)P)Lg,γ2|≲(‖Ta‖H3+‖Ta‖H2‖g-γ‖H3)‖V‖L2‖P‖H1.

Crucially, in the symmetric case, we can bring one term over to the left hand side to show that|(Ta∇^aV,V)Lg,γ2|=|-12((∇^aTa)V,V)Lg,γ2-14(V,(Tagbc∇^agbc)V)Lg,γ2|≲(‖Ta‖H3+‖Ta‖H2‖g-γ‖H3)‖V‖L22.

□

First Estimates on Geometric Components

We now establish some of our first estimates on the ADM variables N, X and the geometric variables g,Σ.

Lemma 5.4

(Dust derivatives of lapse and shift) For 2≤k≤N-1,‖∂uN‖Hk+‖∂uX‖Hk≲Λ(T).

Proof

Since the lapse is a scalar, ∂uN=u0∂TN-τuc∇cN. For ∂uXa we use (5.5). We find that‖∂uN‖Hk+‖∂uX‖Hk≲‖∂TN‖Hk+‖∂TX‖Hk+εe(-1+μ)TΛ(T).

□

Lemma 5.5

(Estimates on geometric source terms) We have‖Fh‖Hk≲‖N^‖Hk+‖X‖Hk+12+‖g-γ‖Hk+12,2≤k≤N-1,‖Fv‖Hk≲‖N^‖Hk+2+|τ|‖ρ‖Hk+‖X‖Hk+12+‖g-γ‖Hk+12+‖Σ‖Hk2,2≤k≤N-2,

and thus‖Fh‖HN-1+‖Fv‖HN-2≲Λ(T)+‖g-γ‖HN2+‖Σ‖HN-12.

Proof

Using the product estimate of Lemma 3.6 and (5.5) we obtain, for 2≤k≤N-1,‖Fh‖Hk≲‖N^‖Hk+‖h‖Hk‖∇X+Υ∗X‖Hk≲‖N^‖Hk+‖g-γ‖Hk‖X‖Hk+1+‖g-γ‖Hk+12‖X‖Hk.

Using, in addition, Lemma 4.4, we see, for 2≤k≤N-2,‖Fv‖Hk≲‖N^‖Hk+2+‖Σ‖Hk2+|τ|‖Sij‖Hk+‖Σ‖Hk‖∇^X‖Hk.

□

We can now combine the geometric source term estimates from the previous lemma with the equations of motion given in (2.8).

Lemma 5.6

(Time derivatives of g and Σ) The following estimates hold:‖∂Tg‖Hk≲‖Σ‖Hk+‖g-γ‖Hk+1+‖N^‖Hk+‖X‖Hk+1,2≤k≤N-1,‖∂TΣ‖Hk≲‖Σ‖Hk+1+‖g-γ‖Hk+2+‖N^‖Hk+2+‖X‖Hk+1+|τ|‖ρ‖Hk,2≤k≤N-2.

Proof

Using Lemma 5.5, (2.8) and (5.5) we find‖∂Th‖Hk≲‖Σ‖Hk+‖X(∇h+Υ∗h)‖Hk+‖Fh‖Hk≲‖Σ‖Hk+‖N^‖Hk+‖X‖Hk+12+‖g-γ‖Hk+12,

and, using additionally Lemma 3.7,‖∂Tv‖Hk≲‖Σ‖Hk+‖Lg,γh‖Hk+‖X(∇v+Υ∗v)‖Hk+‖Fv‖Hk≲‖Σ‖Hk+‖g-γ‖Hk+2+‖N^‖Hk+2+‖X‖Hk+12+‖Σ‖Hk+12+|τ|‖ρ‖Hk.

□

Corollary 5.7

(Dust derivatives of g and Σ) The following estimates hold:‖∂ug‖Hk≲‖Σ‖Hk+‖g-γ‖Hk+1+Λ(T),2≤k≤N-1,‖∂uΣ‖Hk≲‖Σ‖Hk+1+‖g-γ‖Hk+2+Λ(T),2≤k≤N-2.

Proof

Writing∂uVab=u0∂TVab-τuc∇cVab-τucΥcadVdb-τucΥcbdVad

and using Lemma 5.6 gives the estimates. □

Remark 5.8

It is unsurprising that for the geometric variables g,Σ, the fluid derivative estimates in Corollary 5.7 do not gain us any improved information compared to the time derivative estimates in Lemma 5.6. This will of course change when we consider instead estimates on certain fluid matter variables which are naturally more compatible with the fluid derivative operator ∂u.

First Estimates on Fluid Components

In this subsection we establish further estimates on the various matter variables. Note that we need to estimate both matter components coming from contractions with the stress energy tensor (see Definition 2.6) and the fluid source terms terms appearing in the equations of motion (see Definition 2.7).

Lemma 5.9

(Fluid source term estimates) We have, for some ν>0,‖Fρ‖HN-2≲|τ|‖u‖HN-1+‖Σ‖HN-2+Λ(T),‖Fu0‖HN-1≲‖Σ‖HN-12+Λ(T),‖Fuj‖HN-1≲|τ|-1Λ(T)+ε2e-νT.

Proof

For 2≤k≤N-2 we distribute derivatives across the terms given in Definitions 2.4 and 2.7. This yields‖Fρ‖Hk≲|τ|‖∇uj‖Hk+‖u0‖L∞‖Γii‖Hk+|τ|‖Γj‖Hk‖u‖Hk+|τ|‖Γik‖Hk‖X‖Hk‖u‖Hk+τ2‖u‖Hk2‖Γkj‖Hk≲|τ|‖u‖Hk+1+‖Σ‖Hk+‖X‖Hk+1.

Similarly, for 2≤k≤N-1,‖Fu0‖Hk≲‖ΓR‖Hk+|τ|‖u‖Hk‖Γj‖Hk+τ2‖u‖Hk2‖Γjk‖Hk≲‖∂TN‖Hk+‖Σ‖Hk2+‖N^‖Hk+12+‖X‖Hk2+τ2‖u‖Hk2.

Finally, for 2≤k≤N-1,‖Fuj‖Hk≲|τ|-1‖Γ∘j‖Hk+‖u‖Hk2(‖Γ[g]-Γ[γ]‖Hk+|τ|‖X‖Hk‖Γki‖Hk)+‖u‖Hk‖Γji‖Hk≲|τ|-1(‖∂TX‖Hk+‖X‖Hk+1+‖N^‖Hk+12+‖Σ‖Hk2+‖X‖Hk‖∂TN‖Hk)+‖u‖Hk2(‖g-γ‖Hk+1+|τ|‖X‖Hk)+‖u‖Hk(‖Σ‖Hk+‖N^‖Hk+‖X‖Hk+1)≲|τ|-1Λ(T)+ε2e(1-2λ)T+ε2e(-λ+2μ)T.

It is convenient also to note that at lower order we can apply Lemma 4.5 to find5.7 |τ|‖Fuj‖HN-2≲EN-2g+Λ(T)+|τ|‖u‖HN-22(EN-1g)1/2.

□

In the next lemma we provide estimates for fluid derivatives of certain matter components appearing in Definition 2.6. The weights in τ are included for convenience since these expressions appear later on in the energy estimates.

Lemma 5.10

(Dust derivatives of matter components) We have|τ|‖∂uη‖HN-2+|τ|‖∂uS‖HN-2≲εemax{-1+μ,-λ}TΛ(T)+Λ(T)2,|τ|2‖∂uȷ‖HN-2≲Λ(T)2+ε2Λ(T)e(-λ+μ)T.

Proof

Note that (2.7d) can be rewritten as∂uuj=τuaucΥacj+τ-1(u0)2N∇jN+Fuj,∂uu0=Fu0,∂uρ=ρFρ.

Using Definition 2.6 we compute∂uη=ρFρ(u0)2N2+gab(Xau0+τua)(Xbu0+τub)+2ρ(u0)2N∂uN+ρ(∂ugab)(Xau0+τua)(Xbu0+τub)+2ρgabu0(Xbu0+τub)(∂uXa)+Fu02ρu0N2+2ρgabXa(Xbu0+τub)+2ρgabu0ua(Xbu0+τub)(∂Tτ)+2ρgabτ(Xbu0+τub)τuaΥacjuc+τ-1(u0)2N∇jN+Fuj.

Thus, using Lemma 5.4, Corollary 5.7 and the matter estimates from Lemma 5.9, we obtain|τ|‖∂uη‖HN-2≲|τ|‖ρ‖HN-2‖Fρ‖HN-2+|τ|‖ρ‖HN-2‖∂uN‖HN-2s+|τ|‖ρ‖HN-2‖Fu0‖HN-2+h.o.t.≲εemax{-1+μ,-λ}TΛ(T)+Λ(T)2.

Again using (2.6) we compute∂uȷa=(ρFρ)Nu0ua+(∂uN)ρu0ua+(Fu0)ρNua+(τΥbcaubuc+τ-1(u0)2N∇aN+Fua)ρNu0.

So, by the geometry estimates in Lemma 5.4 and 5.6, together with the fluid source term estimates in Lemma 5.9, we obtainτ2‖∂uȷa‖HN-2≲|τ|‖ρ‖HN-2‖N^‖HN-1+|τ|2‖ρ‖HN-2‖Fua‖HN-2+h.o.t.≲ε2Λ(T)e(-λ+μ)T+Λ(T)2.

Finally, from (2.6), we find that∂uSab=(∂uρ)((u0Xa+τua)(u0Xb+τub)+12gab)+12ρ(∂ugab)+h.o.t.,

also|τ|‖∂uS‖HN-2≲|τ|‖ρ‖HN-2‖Fρ‖HN-2≲εemax{-1+μ,-λ}TΛ(T)+Λ(T)2.

□

Remark 5.11

The significance of using dust derivatives is made clear by look at the higher regularity appearing in Lemma 5.10 compared to the following Lemma 5.12 which only concerns time derivatives. In Lemma 5.10, we directly computed the ∂u derivatives using the equations of motion, instead of doing a rough estimate by expanding ∂u∼∂T+τuc∗∇.

Lemma 5.12

(Time derivatives of matter components) We have,|τ|‖∂Tη‖HN-3≲εe(-1+μ)TΛ(T)+Λ(T)2,Λ(T)‖∂Tȷ‖HN-3≲ε2Λ(T)+ε4e-(1+ν)T,|τ|‖∂TS‖HN-3≲Λ(T).

Proof

We calculate ∂Tη explicitly from (2.6) as∂Tη=∂TE+(XaXa(u0)2+2τXbubu0+τ2gabuaub)∂Tρ+(2ρXa(u0)2+2τρuau0)∂TXa+(2u0ρXaXa+2τρXbub)∂Tu0+(2τρXbu0+2τ2ρub)∂Tub+(2ρXbubu0+2τρuaua)∂Tτ.

Using ∂Tτ=-τ, we obtain5.8 |τ|‖∂Tη‖HN-3≲|τ|‖∂TE‖HN-3+Λ(T)|τ|‖∂Tρ‖HN-3+εe(-1+μ)TΛ(T)‖∂TX‖HN-3+Λ(T)2‖∂Tu0‖HN-3+εe(-1+μ)TΛ(T)‖τ∂Tub‖HN-3+Λ(T)2.

To estimate ∂TE we use the rescaled continuity equation [1, Eq 10.16]:∂TE=(3-N)E-Xa∇aE+τN-1∇a(N2ja)-τ2N3gabT_ab-τ2NΣabT_ab.

By Lemma 4.4, we obtain|τ|‖∂TE‖HN-3≲|τ|‖E‖HN-2(‖N^‖HN-3+‖X‖HN-3)+|τ|2‖ȷ‖HN-2+|τ|3‖T_‖HN-3≲εe(-1+μ)TΛ(T)+Λ(T)2.

To estimate ∂Tρ,∂Tu0,∂Tua we use the equations of motion (2.7d) together with Lemma 5.9. We find, 5.9a |τ|‖∂Tρ‖HN-3≲|τ|‖ρ‖HN-3‖Fρ‖HN-3+|τ|2‖u‖HN-3‖ρ‖HN-2≲εemax{-1+μ,-λ}TΛ(T)+Λ(T)2,

and5.9b ‖∂Tu0‖HN-2≲‖Fu0‖HN-2+|τ|‖u‖HN-2‖∇u0‖HN-2≲Λ(T)+‖Σ‖HN-22,|τ|‖∂Tua‖HN-2≲|τ|‖Fua‖HN-2+‖N^‖HN-1+τ2‖u‖HN-12≲Λ(T)+ε2e-(1+ν)T.

Putting all these estimates into (5.8) gives, for 2≤k≤N-3,|τ|‖∂Tη‖Hk≲εe(-1+μ)TΛ(T)+Λ(T)2.

Next, and again using (2.6), we compute∂Tȷa=(∂Tρ)Nu0ua+(∂TN)ρu0ua+(∂Tu0)ρNua+(∂Tua)ρNu0,

So thatΛ(T)‖∂Tȷa‖HN-3≲ε|τ|‖u‖HN-3‖∂Tρ‖HN-3+ε|τ|‖ρ‖HN-3(‖∂TN‖HN-3‖u‖HN-3+‖∂Tu0‖HN-3‖u‖HN-3+‖∂Tua‖HN-3)≲ε2Λ(T)+ε4e-(1+ν)T.

Finally, from (2.6) we calculate (written schematically)∂TS=(∂Tρ)((u0Xa+τua)2+12g)+12ρ(∂Tg)+ρ∂Tg·(u0Xc+τuc)2+ρ∂T(u0Xc+τuc)·(u0Xd+τud).

Using Lemma 5.6 and (5.9) we have|τ|‖∂TS‖HN-3≲|τ|(‖ρ‖HN-3‖Fρ‖HN-3+|τ|‖u‖HN-3‖ρ‖HN-2)+|τ|‖ρ‖HN-3×(h.o.t.)≲|τ|‖ρ‖HN-2.

□

Second Commutator Estimates

We are now in a position to compute various commutator estimates which are required in the later energy estimates. We let V be a symmetric (0, 2)-tensor on M unless otherwise specified. The first lemma looks at the commutator between the dust derivative ∂u with other first-order differential operators.

Lemma 5.13

For 0≤k≤N-2, we have‖[∂u,∂T]V‖Hk≲εe(-1+μ)T‖V‖Hk+1+εe(-1+μ)T‖∂TV‖Hk,‖[∂u,∇^]V‖Hk≲εe(-1+μ)T‖V‖Hk+1+Λ(T)‖∂TV‖Hk,‖[∂u,∇]V‖Hk≲εemax{-1+μ,-λ}T‖V‖Hk+1+Λ(T)‖∂TV‖Hk.

Proof

A computation gives5.10 [∂u,∂T]Vij=-τua∇^aVij+τ∂Tua·∇^aVij-∂Tu^0·∂TVij,[∂u,∇^b]Vij=-τua[∇^a,∇^b]Vij-∇^bu^0·∂TVij+τ∇^bua·∇^aVij,[∂u,∇b]Vij=u0[∂T,∇b]Vij-τuc[∇^c,∇b]Vij-∇bu^0·∂TVij+τ∇buc·∇^cVij.

Thus, by (4.2) and (5.9), if 2≤k≤N-2,‖[∂u,∂T]V‖Hk≲|τ|(‖u‖Hk+‖∂Tua‖Hk)‖∇V+Υ∗V‖Hk+‖∂Tu0‖Hk‖∂TV‖Hk≲εe(-1+μ)T(‖V‖Hk+1+‖g-γ‖Hk+1‖V‖Hk)+εe(-1+μ)T‖∂TV‖Hk.

The estimates for k=0,1 follow in the same way.

The other two estimates follow in a similar way. Note that for [∂u,∇] we use Lemma 5.2, and Eqs. (5.1), (5.2) and (5.9). □

The next lemma in this subsection investigates the commutator between the second-order operator Lg,γ and other first-order operators.

Lemma 5.14

The following estimates hold:‖[∂u,Lg,γ]V‖Hk≲εemax{-1+μ,-λ}T‖V‖Hk+2+Λ(T)‖∂TV‖Hk+1,0≤k≤N-3,‖[∇^m,Lg,γ]V‖Hk≲εe-λT‖V‖Hk+2+‖V‖Hk,0≤k≤N-2,‖[∂T,Lg,γ]V‖Hk≲εe-λT‖V‖Hk+2,0≤k≤N-2.

Also, for k,s∈Z such that 0≤k≤N-2, s≥1 and 2(s-1)+k≤N-1, we have‖[∂T,Lg,γs]V‖Hk≲εe-λT‖V‖Hk+2+2(s-1).

Proof

A calculation yields5.11 [∂u,Lg,γ]Vij=-(∂ugab)∇^a∇^bVij+Δ^g,γu^0·∂TVij+2gab∇^au^0∇^b∂TVij+τucgab(Riem[γ]back∇^kVij+4Riem[γ](i|ack∇^bVk|j)+2∇^aRiem[γ](i|bck·Vk|j))-2τgab∇^auc∇^b∇^cVij-τΔ^g,γuc·∇^cVij+2τuc∇^cRiem[γ]iajb·Vab.

It is useful to write this schematically, using that Riem[γ] and its derivatives are bounded by constants[∂u,Lg,γ]V=(∂ug-1)∗∇^2V+∇^2u^0∗∂TV+∇^u0∗∇^∂TV+τ(uc∗(V+∇^V)+∇^uc∗∇^2V+∇^2uc∗∇^V).

If 2≤k≤N-3 then by elliptic regularity of Lemma 3.7 and the commutator estimates in Lemma 5.2,5.12 ‖[∂u,Lg,γ]V‖Hk≲‖∂ug‖Hk‖V‖Hk+2+‖u^0‖Hk+2‖∂TV‖Hk+1+|τ|‖u‖Hk+2‖V‖Hk+2.

The conclusion then holds by (4.2) and by estimating ∂ug using Corollary 5.7 . The estimates when k=0,1 follow in a similar way.

Next we compute the identity5.13 [Lg,γ,∇^m]Vij=∇^mgab·∇^a∇^bVij-gab[∇^a∇^b,∇^m]Vij,

where we are thinking of the m index as not being free (i.e. contracted with a factor of the shift Xm). Since the commutator involving only ∇^ will just generate background Riemann curvature components the required estimate follows straightforwardly. We note also that (5.13) and the elliptic regularity of Lemma 3.7 imply, for s∈Z such that 1≤s≤N/2,5.14 ‖[Lg,γs,∇^m]V‖L2≲‖[Lg,γ,∇^]Lg,γs-1V‖L2+⋯+‖[Lg,γ,∇^]V‖H2(s-1)≲‖g-γ‖Hmax{3,2s-1}(‖V‖H2s+‖g-γ‖H2s‖V‖H2s-1)+‖V‖H2(s-1)≲εe-λT‖V‖H2s+‖V‖H2(s-1).

Finally we compute[∂T,Lg,γ]Vij=-[∂T,Δ^g,γ]Vij=-(∂Tgab)∇^a∇^bVij.

Using Lemma 5.6 to estimate ∂Tg this gives, for 2≤k≤N-2,‖[∂T,Lg,γ]V‖Hk≲‖∂Th‖Hk(‖∇(∇V+Υ∗V)‖Hk+‖Υ(∇V+Υ∗V)‖Hk)≲εe-λT‖V‖Hk+2.

A similar argument holds for the cases k=0,1. At higher order, we obtain the identity5.15 [∂T,Lg,γs]Vij=-∑1≤i≤sLg,γi-1(∂Tgaibi)·∇^ai∇^bi(Lg,γs-i(Vij)).

This can be estimated, for 2≤k≤N-2, by5.16 ‖[∂T,Lg,γs]V‖Hk≲‖∂Th‖H2(s-1)+k(‖V‖Hk+2+2(s-1)+‖g-γ‖Hk+2‖V‖Hk+1+2(s-1)).

The cases k=0,1 are treated in a similar way and the conclusion follows from Lemma 5.6. □

The next corollary extends the commutator estimates of the previous lemma to higher-orders of Lg,γs, s∈Z≥1. Note that the L2 estimate that appears in the statement will be typically applied with V=h, while the lower-order H1 estimate will be primarily used later on with V=Σ.

Corollary 5.15

We have‖[∂u,Lg,γs]V‖L2≲εemax{-1+μ,-λ}T‖V‖H2s+Λ(T)‖∂TV‖H2s-1,1≤s≤ℓ,‖[∂u,Lg,γs-1]V‖H1≲εemax{-1+μ,-λ}T‖V‖H2s-1+Λ(T)‖∂TV‖H2s-2,2≤s≤ℓ.

Proof

By Lemma 5.14,‖[∂u,Lg,γs]V‖L2≲‖[∂u,Lg,γ]Lg,γs-1V‖L2+‖[∂u,Lg,γ]Lg,γs-2V‖H2+⋯+‖[∂u,Lg,γ]V‖H2(s-1)≲εemax{-1+μ,-λ}T‖V‖H2s+Λ(T)‖∂TV‖H2s-1+Λ(T)∑p=0s-2‖[∂T,Lg,γs-1-p]V‖H1+2p.

The conclusion then easily follows and the second estimate follows in the same way. □

Corollary 5.16

‖∂uLg,γℓ-1(Σ)‖L2≲(EN-1g)1/2+Λ(T),‖∂uLg,γℓ-1(Σ)‖H1≲‖Σ‖HN-1+‖g-γ‖HN+Λ(T).

Proof

By the ∂uΣ,∂TΣ estimates of Lemma 5.6, Corollary 5.7, and the previous commutator estimate of corollary 5.15‖∂uLg,γℓ-1(Σ)‖L2≲‖∂uΣ‖HN-3+‖[∂u,Lg,γℓ-1](Σ)‖L2≲‖Σ‖HN-2+‖g-γ‖HN-1+Λ(T).

The conclusion then follows by the coercive estimate of Lemma 4.5. Similarly,5.17 ‖∂uLg,γℓ-1Σ‖H1≲‖Σ‖HN-1+‖g-γ‖HN+Λ(T)+εemax{-1+μ,-λ}T‖Σ‖HN-2+Λ(T)‖∂TΣ‖HN-3≲‖Σ‖HN-1+‖g-γ‖HN+Λ(T).

□

Remark 5.17

Frequently in our energy estimates we will need to study the term ‖∂uLg,γℓ-1(Σ)‖H1 appearing in the previous Corollary. However, by looking at the estimate derived in Corollary 5.16, we see that we cannot apply Lemma 4.5 to the top-order Sobolev norms ‖Σ‖HN-1 and ‖g-γ‖HN. To estimate these terms by the geometric energy Etotg we will instead need to use the auxiliary elliptic estimates established in Sect. 6.

We conclude this subsection with a commutator estimate that plays an important role in the proof of the lapse estimate appearing in Proposition 7.3.

Lemma 5.18

Let Δ:=gab∇a∇b. Then for 0≤k≤N-3,‖[Δ,∂u]V‖Hk≲Λ(T)‖∂TV‖Hk+1+εemax{-1+μ,-λ}T‖V‖Hk+2.

Proof

We compute[Δ,∂u]Vij=Δu^0·∂TVij+2∇au^0∇a(∂TVij)-u0∂Tgab·∇a∇vVij+u0gab[∇a∇b,∂T]Vij+τ(-Δuc·∇^cVij-2∇auc∇a∇^cVij+uc∇^cgab·∇a∇bVij)+τucgab[∇^c,∇a∇b]Vij.

Let 2≤k≤N-3. From Lemma 5.2 and (5.3) we find‖[Δ,∂u]V‖Hk≲‖∇u^0‖Hk+1‖∂TV‖Hk+1+‖∂Th‖Hk‖V‖Hk+2+‖∂TΓ‖Hk+1‖V‖Hk+1+|τ|‖u‖Hk+2‖V‖Hk+2.

The cases k=0,1 follow in the same way. The conclusion follows using Lemma 5.6, (4.2) and (5.2). □

Second Geometric Components Estimates

We now reach the final subsection of Sect. 5. The first lemma is an analogue of Lemma 4.5 for our top-order ∂u-boosted geometric energy. The proof follows those in [1, Lemma 19] and [3, Lemma 7.2].

Lemma 5.19

Let s∈Z such that 1≤s≤ℓ. There is a δ>0 and a constant C>0 such that for (g,Σ)∈Bδ((γ,0)) the inequality(∂uLg,γsh,∂uLg,γsh)Lg,γ2+|(∂uLg,γsv,∂uLg,γs-1v)Lg,γ2|≤CE∂u,2ℓg

holds. Furthermore E∂u,2ℓg≥0.

Proof

Recall E∂u,2ℓg from Definition 3.11. We first note that E∂u,2ℓg|(h,v)=(γ,0)=0. Next, we see that (γ,0) is a critical point of E∂u,2ℓg since the first derivative vanishes. Considering then the second derivative of this energy at (γ,0), we see that the Hessian takes the formD2(E∂u,2sg+cEΓ∂u,2sg)((h,k),(h,k))=9(∂uLγ,γsh,∂uLγ,γsh)Lg,γ2+(∂uLγ,γsk,∂uLγ,γs-1k)Lg,γ2+cE(∂uLγ,γs-1k,∂uLγ,γsh)Lg,γ2.

We claim that the Hessian is non-negative. By expanding in terms of the eigentensors of Lγ,γ, we are left with terms of the typeλ2s-1(9λ(∂uPλh,∂uPλh)Lg,γ2+(∂uPλk,∂uPλk)Lg,γ2+cE(∂uPλk,∂uPλh)Lg,γ2),

where Pλ denotes the projection operator onto the λ-eigenspace. The choice of cE ensures that the bracketed term is non-negative for the smallest eigenvalue λ0, which in turn implies non-negativity for all eigenvalues. Thus we findD2(E∂u,2sg+cEΓ∂u,2sg)((h,k),(h,k))≥0.

From this it follows that there is a constant C=C(λ0,γ)>0 such that∑s=1ℓ(∂uLg,γsh,∂uLg,γsh)Lg,γ2+|(∂uLg,γsk,∂uLg,γs-1k)Lg,γ2|≤CE∂u,2ℓg.

□

The next lemma provides estimates for the Sobolev norms of ∂uh,∂uΣ in terms of the geometric energy Etotg. This is natural given that we have constructed the functional Etotg to precisely control such Sobolev norms.

Lemma 5.20

(Dust derivatives of geometric variables at high-regularity) We have,‖∂uh‖HN-12≲E∂u,N-1g+ε2e2max{-1+μ,-λ}TEN-1g+Λ(T)4,‖∂uΣ‖HN-22≲E∂u,N-1g+EN-1g+εemax{-1+μ,-λ}T×(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

Remark 5.21

Similar to Remark 5.17, we cannot apply Lemma 4.5 to the top-order Sobolev norms ‖Σ‖HN-12 and ‖g-γ‖HN2. To estimate these terms by the geometric energy Etotg we will instead need to use the auxiliary elliptic estimates of Sect. 6. It is also crucial in later analysis in Sect. 6 that these top-order norms above appear on the right hand side in Lemma 5.20 with a smallness factor of ε.

Proof

(Proof of Lemma 5.20) Recall that N:=2ℓ+1. By the elliptic regularity of Lemma 3.7 and Lemma 5.19‖∂uh‖H2ℓ2≲‖Lg,γℓ∂uh‖L22≲E∂u,N-1g+‖[∂u,Lg,γℓ]h‖L22.

We control the commutator term using Corollary 5.15, finding5.18 ‖[∂u,Lg,γℓ]h‖L2≲εemax{-1+μ,-λ}T‖g-γ‖HN-1+Λ(T)‖∂Th‖HN-2≲εemax{-1+μ,-λ}T(EN-1g)1/2+Λ(T)2,

where in the final line we used Lemma 5.6 and the coercive estimate of Lemma 4.5.

Next, by elliptic regularity and Lemma 4.6, we have5.19 ‖∂uΣ‖H2ℓ-12≲‖Lg,γℓ-1∂uΣ‖H12≲‖∂uLg,γℓ-1Σ‖L22+‖[∂u,Lg,γℓ-1]Σ‖L22+(Lg,γℓ-1∂uΣ,Lg,γℓ∂uΣ)Lg,γ2.

The first term on the RHS of (5.19) is treated by Corollary 5.16. For the commutator term in (5.19), we use Lemma 5.6 and Corollary 5.15 to find5.20 ‖[∂u,Lg,γℓ-1]Σ‖L2≲εemax{-1+μ,-λ}T‖Σ‖HN-3+Λ(T)‖∂TΣ‖HN-4≲εemax{-1+μ,-λ}T(EN-2g)1/2+Λ(T)2.

Considering next the final term in (5.19), we write it as(Lg,γℓ-1∂uΣ,Lg,γℓ∂uΣ)Lg,γ2=(∂uLg,γℓ-1Σ,∂uLg,γℓΣ)Lg,γ2+([Lg,γℓ-1,∂u]Σ,∂uLg,γℓΣ)Lg,γ2+(∂uLg,γℓ-1Σ,[∂u,Lg,γℓ]Σ)Lg,γ2+([Lg,γℓ-1,∂u]Σ,[∂u,Lg,γℓ]Σ)Lg,γ2=:E1+E2+E3+E4.

From Lemma 5.19, |E1|≲E∂u,N-1g. To estimate E2 we need to integrate by parts one of the derivatives appearing in ∂uLg,γℓΣ. Using (5.6a) we findE2=(gab∇^a[Lg,γℓ-1,∂u]Σ,∇^b∂uLg,γℓ-1Σ)Lg,γ2-2([Lg,γℓ-1,∂u]Σ,Riem[γ]∘∂uLg,γℓ-1Σ)Lg,γ2+([Lg,γℓ-1,∂u]Σ,[∂u,Lg,γ]Lg,γℓ-1Σ)Lg,γ2.

By the commutator estimates of Lemma 5.14, and Lemma 5.6, we have5.21 ‖[∂u,Lg,γ]Lg,γℓ-1Σ‖L2≲εemax{-1+μ,-λ}T‖Σ‖HN-1+Λ(T)‖∂TΣ‖HN-2≲εemax{-1+μ,-λ}T(‖Σ‖HN-1+‖g-γ‖HN)+Λ(T)2.

Using this, together with (5.20) and Corollary 5.16, gives|E2|≲‖[∂u,Lg,γℓ-1]Σ‖H1(‖∂uLg,γℓ-1Σ‖H1+‖[∂u,Lg,γ]Lg,γℓ-1Σ‖L2)≲εemax{-1+μ,-λ}T(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

The terms E3,E4 are similarly estimated and inserting all these estimates into (5.19) gives the required result. □

We end this subsection with two Lemmas concerning the geometric source terms Fh and Fv.

Lemma 5.22

(Time derivative of geometric source terms) We have‖∂TFh‖HN-2+‖∂TFv‖HN-3≲EN-1g+Λ(T).

Proof

Let 2≤k≤N-2. Taking a time derivative of Fh as given in Definition 2.8 we see that‖∂TFh‖Hk≲‖∂TN‖Hk+(‖N^‖Hk+‖X‖Hk+1+‖g-γ‖Hk+1‖X‖Hk)‖∂Th‖Hk+‖g-γ‖Hk(‖∂TX‖Hk+1+‖g-γ‖Hk+1‖∂TX‖Hk).

The first estimate then follows by applying Lemma 5.6.

Next, let 2≤k≤N-3. We take a time derivative of Fv which gives∂TFv=∂T∇i∇jN+∂TN·(2ΣicΣjc-13gij-Σij+τSij)+2N∂T(ΣicΣjc)-13N^∂Tgij-N^∂TΣij-τNSij+τN∂TSij-∂T(vim∇^jXm+vmj∇^iXm).

For the first term we use the commutator identity from (5.3) and recall that the lapse is a scalar. We obtain‖∂TFv‖Hk≲‖∂TN‖Hk+2+‖[∂T,∇]∇N‖Hk+|τ|(‖S‖Hk+‖∂TS‖Hk)+‖Σ‖Hk2+‖∂TΣ‖Hk2+‖∂Tg‖Hk2+Λ(T)2.

The conclusion then follows by Lemma 5.6 and the matter estimates in Lemma 4.4 and Lemma 5.12. □

Corollary 5.23

(Dust derivative of geometric source terms) We have‖∂uFh‖HN-2+‖∂uFv‖HN-3≲EN-1g+εe(-1+μ)T(‖g-γ‖HN2+‖Σ‖HN-12)+Λ(T),‖∂uFh‖H2≲EN-1g+Λ(T).

Proof

The estimates follow by expanding out ∂u=u0∂T-τuc∇^c and using the Fh,Fv estimates in Lemma 5.5 and the ∂TFh,∂TFv estimates from Lemma 5.22. □

Remark 5.24

Since we are dealing with geometric variables in the above corollary, and not matter variables, we roughly estimated the dust derivatives as ∂u∼∂T+τuc∗∇. Note that doing so introduced the top-order Sobolev norms ‖Σ‖HN-12 and ‖g-γ‖HN2. Crucially for later analysis in Corollary 6.4, however, is that these top-order norms appear with a coefficient of ε.

Elliptic Estimate

In this section we prove an auxiliary elliptic estimate which allows us to control the top-order Sobolev norms of g and Σ in terms of the geometric energy functional Etotg. Recall only the lower-order Sobolev norms are controlled using Lemma 4.5, and so a new idea is indeed needed to cover the top-order of regularity. We also remind the reader that the geometric energy Etotg will eventually fulfill a strong decay estimate enabled by the inclusion of certain correction terms.

The main result of the section, Corollary 6.4, achieves the goal of the previous paragraph. To prove this corollary, we take the first-order equations of motion for ∂Th,∂TΣ appearing in (2.8) and convert them into second-order equations involving a perturbed wave operator W and the ∂u derivatives. We find, very schematically, that6.1 Lg,γh∼W(h)+∂u∂T(h)+τua∇^a∂uh.

By elliptic regularity for Lg,γ, we can then prove an estimate on ‖g-γ‖HN∼‖Lg,γh‖HN-2 by estimating the RHS of (6.1), see Proposition 6.3. A similar idea holds also for Σ. We note that this idea, albeit for a different gauge, was first introduced by Hadžić and Speck in [16].

Definition 6.1

(Operators W,H) Define the operatorsW:=∂T∂T-XaXb∇^a∇^b+N2Lg,γ,H:=N2Lg,γ-XaXb∇^a∇^b-τ2uaub(u0)2∇^b∇^a,

which act on symmetric (0, 2)-tensors.

Due to the sign convention on Lg,γ, see Definition 3.2, one can think of W as being a kind of perturbed wave operator.

Lemma 6.2

(Wave equations for h, v) The differential Eq. (2.8) for h=g-γ and v=6Σ implyW(h)=F1,W(v)=F2,

where,‖F1‖HN-22+‖F2‖HN-32≲EN-1g+εe-2λT(‖g-γ‖HN2+‖Σ‖HN-12)+Λ(T)2.

Proof

Rearranging (2.8) as v=w-1(∂Th+Xm∇^mh-Fh) and substituting this into (2.8) gives-w-1(∂Tw)v+w-1(∂T2h+∂T(Xm∇^mh)-∂TFh)=-2v-9wLg,γh-Xm∇^mv+6Fv.

Using again (2.8) we note that∂T(Xa∇^ah)=-XaXb∇^a∇^bh-Xa∇^aXb·∇^bh+∂TXa·∇^ah+Xm∇^m(wv+Fh).

Rearranging terms (recall 9w2=N2) we findW(h)=F1:=Xa∇^aXb·∇^bh-∂TXm·∇^mh-Xm∇^m(wv+Fh)+∂TFh+v∂Tw-2wv-wXm∇^mv+6wFv.

Using Lemma 5.5 and Lemma 5.22 we see that‖F1‖HN-2≲‖Σ‖HN-2+‖∂TFh‖HN-2+‖Fv‖HN-2+Λ(T)≲EN-1g+‖g-γ‖HN2+‖Σ‖HN-12+Λ(T).

Although the top-order terms involving g-γ and Σ here look worrying, when squaring the estimate we can then apply the bootstraps to gain the crucial factor of ε.

To derive the equation for v we take the ∂T derivative of (2.8):∂T2v=-2∂Tv-9∂Tw·Lg,γh-9wLg,γ(∂Th)-9w[∂T,Lg,γ]h-∂TXm·∇^mv-Xm∇^m(∂Tv)+6∂TFv.

Substituting in the equation of motion (2.8) where needed, and expanding as Lg,γ(wv)=wLg,γv-vΔ^g,γw-2gab∇^aw∇^bv, we obtainW(v)=F2:=4v+18wLg,γh+2Xm∇^mv-12Fv-9∂Tw·Lg,γh+9wvΔ^g,γw+18wgab∇^aw∇^bv+9wLg,γ(Xm∇^mh)-9wLg,γFh-9w[∂T,Lg,γ]h+6∂TFv-∂TXm·∇^mv+2Xm∇^mv+9Xm∇^m(wLg,γh)+Xa∇^aXm·∇^mv-6Xm∇^mFv.

Using the Fh,Fv estimates in Lemma 5.5 and Lemma 5.22, together with the commutator estimate in Lemma 5.14, we find‖F2‖HN-3≲‖Σ‖HN-3+‖g-γ‖HN-1+‖[∂T,Lg,γ]h‖HN-3+‖Fv‖HN-2+‖∂TFv‖HN-3+‖Fh‖HN-1+Λ(T)≲(EN-1g)1/2+εe-λT‖g-γ‖HN-1+‖g-γ‖HN2+‖Σ‖HN-12+EN-1g+Λ(T).

Note that in the above we also used (5.5) to estimate a term of the form‖∇^h‖HN-1≲‖h‖HN+‖Υ‖HN-1‖h‖HN-1≲‖g-γ‖HN.

□

Proposition 6.3

(Elliptic estimate using W,∂u operators) Let V be an arbitrary (0, 2)-tensor and 0≤k≤N-2. Then,‖V‖Hk+2≲‖W(V)‖Hk+‖∂u∂T(V)‖Hk+|τ|‖ua∂u∇^a(V)‖Hk.

Proof

From (3.9) we have∂TV=(u0)-1(∂u(V)+τua∇^aV),

and thus∂T∂TV=(u0)-1(∂u(∂TV)+τua∇^a∂TV),∂T∇^bV=(u0)-1(∂u(∇^bV)+τua∇^a∇^bV).

Recalling u^0=u0-1/3, we find6.2 W(V)=∂T∂TV-XaXb∇^a∇^bV+N2Lg,γV=(u0)-1(∂u(∂TV)+τua∂T∇^aV)-XaXb∇^a∇^bV+N2Lg,γV=(u0)-1∂u(∂TV)+τuau01u0∂u(∇^aV)+τub∇^b∇^aV-XaXb∇^a∇^bV+N2Lg,γV=(u0)-1∂u(∂TV)+τua(u0)-2∂u(∇^aV)+H(V).

We view N-2H as a perturbation off the elliptic operator Lg,γ. For 0≤k≤N-2,‖(N-2H-Lg,γ)V‖Hk≲‖XaXb∇^a∇^bV‖Hk+τ2‖uaub(u0)2∇^b∇^aV‖Hk≲ε2e(-2+2μ)T(‖V‖Hk+2+‖Υ‖Hk+1‖V‖Hk+1)≲ε2e(-2+2μ)T‖Lg,γV‖Hk.

Suppose now H(V)=0. By definition of H, and elliptic regularity of Lg,γ, this implies‖N2Lg,γ(V)‖L2=‖XaXb∇^a∇^bV+τ2(u0)-2uaub∇^b∇^aV‖L2≲ε2e2(-1+μ)T‖V‖H2≤C(C1)-1ε2e2(-1+μ)T‖Lg,γV‖L2,

for C>0 some constant and C1>0 as in Lemma 3.7. We also have a lower boundC′‖Lg,γ(V)‖L2≤‖N2Lg,γ(V)‖L2.

for another constant C′>0. Choosing ε sufficiently small so that C(C1)-1ε2e2(-1+μ)T<ε, we see these two inequalities implyC′‖Lg,γ(V)‖L2<ε‖Lg,γ(V)‖L2.

For ε sufficiently small this implies ‖Lg,γ(V)‖L2=0 and so V∈kerLg,γ. However, kerLg,γ=0, and so for small data H also has trivial kernel and thus we obtain‖V‖Hk+2≲‖N-2H(V)‖Hk≲‖V‖Hk+2.

Putting this together with (6.2) we find‖V‖Hk+2≲‖N-2W(V)‖Hk+‖(u0)-1N-2∂u∂T(V)‖Hk+|τ|‖(u0)-2uaN-2∂u∇^a(V)‖Hk.

□

We can now bring together the previous results and estimate the top-order Sobolev norms of g,Σ in terms of our geometric energy functionals E∂u,N-1g and EN-1g.

Corollary 6.4

We have,‖g-γ‖HN2+‖Σ‖HN-12≲E∂u,N-1g+EN-1g+Λ(T)2.

Proof

From Proposition 6.3,‖g-γ‖HN2≲‖W(h)‖HN-22+‖∂u∂T(h)‖HN-22+|τ|2‖ua∂u∇^a(h)‖HN-22.

The first term here is treated using Lemma 6.2. For the second term, we begin by using the commutator estimates in Lemma 5.13, and the ∂Th and ∂uh estimates in Lemma 5.6 and Lemma 5.20 respectively, to find6.3 ‖∂u∇^(h)‖HN-22≲‖∂uh‖HN-12+ε2e(-2+2μ)T‖h‖HN-12+Λ(T)2‖∂Th‖HN-22≲E∂u,N-1g+ε2e2max{-1+μ,-λ}TEN-1g+Λ(T)3,

Next, using the expression for ∂Th given in (2.8), together with Lemma 5.4, Corollary 6.5, Lemma 5.20 and Corollary 5.23, we find‖∂u∂T(h)‖HN-22≲‖∂uΣ‖HN-22+‖Σ‖HN-22‖∂uN‖HN-22+‖∂uX‖HN-22‖∇^h‖HN-22+‖X‖HN-22‖∂u∇^h‖HN-22+‖∂uFh‖HN-22.≲E∂u,N-1g+EN-1g+εemax{-1+μ,-λ}T(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

Using again (6.3) we find|τ|2‖ua∂u∇^(h)‖HN-22≲|τ|2‖u‖HN-22(E∂u,N-1g+ε2emax{-2+2μ,-2λ}TEN-1g+Λ(T)3)≲E∂u,N-1g+EN-1g+Λ(T)4.

Putting this all together,6.4 ‖g-γ‖HN2≲E∂u,N-1g+EN-1g+ε(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

We follow the same steps for the Σ estimate. From Proposition 6.3,‖Σ‖HN-12≲‖W(v)‖HN-32+‖∂u∂T(v)‖HN-32+|τ|2‖ua∂u∇^a(Σ)‖HN-32.

The first term here is treated using Lemma 6.2. For the second term, we begin by using Lemma 5.6, Lemma 5.20 and the commutator estimates in Lemma 5.13 to show that6.5 ‖∂u∇^(Σ)‖HN-32≲‖∂uΣ‖HN-22+ε2e(-2+2μ)T‖Σ‖HN-22+Λ(T)2‖∂TΣ‖HN-32≲E∂u,N-1g+EN-1g+εemax{-1+μ,-λ}T(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

Next, by Lemma 5.6, Lemma 5.20 and the commutator estimate of Lemma 5.14, we note‖∂uLg,γh‖HN-32≲‖∂uh‖HN-12+‖[∂u,Lg,γ]h‖HN-32≲E∂u,N-1g+ε2emax{-2+2μ,-2λ}TEN-1g+Λ(T)2(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)4.

We can now use the expression for ∂Tv given in (2.8) and bring together these previous estimates:‖∂u∂T(v)‖HN-32≲‖∂uΣ‖HN-32+‖Lg,γh‖HN-32‖∂uN‖HN-32+‖∂uLg,γh‖HN-32+‖∂uX‖HN-32‖∇^Σ‖HN-32+‖X‖HN-32‖∂u∇^Σ‖HN-32+‖∂uFv‖HN-32≲‖∂uΣ‖HN-32+‖∂uLg,γh‖HN-32+‖∂uFv‖HN-32+Λ(T)2≲E∂u,N-1g+εemax{-1+μ,-λ}T(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

In the above we used Lemma 5.4, Corollary 6.5, Lemma 5.23 and Lemma 5.20. Using again (6.5) and Lemma 5.20, we find|τ|2‖ua∂u∇^(Σ)‖HN-32≲E∂u,N-1g+EN-1g+εemax{-1+μ,-λ}TΛ(T)×(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)3.

Finally, putting this all together gives6.6 ‖Σ‖HN-12≲E∂u,N-1g+EN-1g+ε(‖Σ‖HN-12+‖g-γ‖HN2)+Λ(T)2.

The conclusion then follows by adding the estimates (6.4) and (6.6) together and taking ε sufficiently small so that we can absorb the ε(‖Σ‖HN-12+‖g-γ‖HN2) term onto the left hand side. □

We now provide new estimates on dust derivatives acting on our variables N,X,g,Σ and, for the latter two variables, apply Corollary 6.4.

Corollary 6.5

For I a multi-index,∑|I|≤N-1‖∂u∇IX‖L2+‖∂u∇IN^‖L2≲Λ(T),∑|I|≤N-1‖∂u∇I(g-γ)‖L2+∑|I|≤N-2‖∂u∇IΣ‖L2≲(E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T).

The same estimates hold with ∇ replaced by ∇^.

Proof

By Lemma 5.20 and Corollary 6.4‖∂uh‖HN-1+‖∂uΣ‖HN-2≲(E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T).

Next, let ∇I=∇a1…∇ak where k:=|I|≤N-1 and let V be a (0, 2)-tensor on M. Then, by (5.1), (5.2) and Lemma 5.13,‖[∂u,∇I]V‖L2≲‖[∂u,∇a1]∇a2…∇akV‖L2+⋯+‖[∂u,∇]V‖Hk-1≲εemax{-1+μ,-λ}T‖V‖Hk+Λ(T)‖∂TV‖Hk-1+Λ(T)‖∂TΓ(g)‖HN-3‖V‖Hk-2.

This implies‖∂u∇IV‖L2≲‖∂uV‖Hk+εemax{-1+μ,-λ}T‖V‖Hk+Λ(T)‖∂TV‖Hk-1.

A simple check shows that when considering the appropriate commutator expression (see (5.10)) on the scalar lapse we will pick up ‖N^‖Hk terms and not ‖N‖Hk. The result then follow, using Lemma 5.4 and Lemma 5.6. □

Lemma 6.6

Let V be a (0, 2)-tensor on M and let p∈{1,2}. Then‖∂uV‖Hp≲∑|I|≤p‖∂u∇IV‖L2+εemax{-1+μ,-λ}T‖V‖H1+Λ(T)‖∂TV‖Hp-1.

Proof

A straight forward application of commutator identities. □

Lapse and Shift Estimates

In this section we first establish the basic estimates on the lapse, shift and their time derivatives using elliptic estimates and general formula presented in [1]. The most exciting results lie in the novel top-order lapse and shift estimates given in Sect. 7.1.

Lemma 7.1

For 3≤k≤N,‖N^‖Hk≲|τ|‖ρ‖Hk-2+ε2e-2λT,‖X‖Hk≲|τ|‖ρ‖Hk-3+ε2emax{-2λ,-2+μ}T.

Proof

Recall from (2.7b) the lapse equation of motion(Δ-13)N=N(|Σ|g2+τη)-1.

By elliptic regularity for the standard Laplacian Δ this implies‖N^‖Hk≤C(‖Σ‖k-22+|τ|‖η‖Hk-2),

and the desired result follows by Lemma 4.4.

Similarly, from the shift equation of motion (2.7b), we find that‖X‖Hk≤C(‖Σ‖Hk-22+‖g-γ‖Hk-12+|τ|‖η‖Hk-3+τ2‖Nȷ‖Hk-2),

and the desired result follows by Lemma 4.4. □

Lemma 7.2

For 4≤k≤N-1,‖∂TN‖Hk+‖∂TX‖Hk≲‖N^‖Hk+‖X‖Hk+|τ|‖ρ‖Hk-2+ε2emax{-2λ,-2+μ}T.

Proof

Let 4≤k≤N-1. Following [1], we have‖∂TN‖Hk≲‖N^‖Hk+‖X‖Hk+‖Σ‖Hk-12+‖g-γ‖Hk2+|τ|(‖S‖Hk-2+‖η‖Hk-2+‖∂Tη‖Hk-2)≲‖N^‖Hk+‖X‖Hk+|τ|‖ρ‖Hk-2+ε2emax{-2λ,-2+μ}T,

where we used Lemma 4.4 and Lemma 5.12.

Again, using a general expression given in [1, §7], we have‖∂TX‖Hk≲‖N^‖Hk+‖X‖Hk+‖Σ‖Hk2+‖g-γ‖Hk2+τ(‖S‖Hk-2+‖η‖Hk-2+‖∂Tη‖Hk-2)+τ2‖ȷ‖Hk-1+τ3‖T_‖Hk-1≲‖N^‖Hk+‖X‖Hk+|τ|‖ρ‖Hk-2+ε2emax{-2λ,-2+μ}T,

where we again used Lemma 4.4 and Lemma 5.12. □

Additional Top-order Lapse and Shift Estimates

In this section we establish important auxiliary estimates in two propositions for the lapse and shift variables. The first proposition is used to estimate ∂u∇i∇jN. This somewhat unusual expression comes from a dust derivative acting on the first term in Fv (see Definition 2.8), and appears later on in the energy estimates for Etotg.

Proposition 7.3

(Top-order auxiliary lapse estimate) We have(Δ-13)∂u∇i∇jN=FN,u

where,‖FN,u‖HN-4≲Λ(T)+E∂u,N-1g+EN-1g.

Proof

Start by commuting the lapse Eq. (2.7b) with ∂u∇i∇j:(Δ-13)∂u∇i∇jN=[Δ,∂u]∇i∇jN+∂u[Δ,∇i∇j]N+∂u∇i∇j(N(|Σ|g2-τη))=:L1+L2+L3=:FN,u.

We investigate each of the terms in FN,u separately. The first term L1 we estimate using (5.2), (5.3) and Lemma 5.18‖L1‖HN-4≲Λ(T)‖∂T∇i∇jN‖HN-3+εemax{-1+μ,-λ}T‖∇i∇jN‖HN-2≲Λ(T)(‖∂TN‖HN-1+‖[∂T,∇i]∇jN‖HN-3)+εemax{-1+μ,-λ}T‖N^‖HN≲Λ(T)2+εemax{-1+μ,-λ}TΛ(T).

For the next term, L2, we first compute, for ϕ a scalar,[Δ,∇i]ϕ=2gab∇[a∇i]∇bϕ=-Ricic∇cϕ,[Δ,∇i]∇jϕ=Ricic∇c∇jϕ-(∇aRiemjaic)∇cϕ-2Riemjaic∇a∇cϕ.

ThusL2=∂u([Δ,∇i]∇jN)+∂u∇i([Δ,∇jN)=-((∂u∇aRiemjaic)+(∂u∇iRicjc))∇cN-(∇aRiemjaic+∇iRicjc)∂u∇cN+((∂uRicic)∇c∇jN-(∂uRicjc)∇i∇cN-2(∂uRiemjaic)∇a∇cN)+(Ricic(∂u∇c∇jN)-Ricjc(∂u∇i∇cN)-2Riemjaic(∂u∇a∇cN))=:-(L21)-(L22)+(L23)+(L24).

The second and last terms here are easy to estimate using Lemma 5.4 and Lemma 5.13 (note also that the expressions (5.1) and (5.10) simplify when calculated for a scalar)‖|L22|+|L24|‖HN-4≲‖Riem‖HN-3(‖∂uN‖HN-4+‖[∂u,∇]N‖HN-3+‖[∂u,∇]∇N‖HN-4)≲‖Riem‖HN-3(Λ(T)+Λ(T)‖∂TN‖HN-3+εemax{-1+μ,-λ}T‖N‖HN-2)≲Λ(T).

For the other two terms, L21 and L23, we schematically write∂uRiem=u0∂T(∂Γ+ΓΓ)-τuc∇^cRiem=u0(∂∂TΓ+Γ∂TΓ)-τuc∇cRiem+τuc∗Υ∗Riem,∂u∇Riem=∇∂uRiem+[∂u,∇]Riem,

and thus‖|L21|+|L23|‖HN-4≲‖N^‖HN-2(‖∂TΓ‖HN-2+‖∂TΓ‖HN-3‖Γ‖HN-3+‖∂TΓ‖HN-4‖Riem‖HN-4+εemax{-1+μ,-λ}T‖Riem‖HN-2+ε‖∂TRiem‖HN-4)≲εe-λT‖N^‖HN-2+εemax{-1+μ,-λ}T‖N^‖HN-2‖Riem‖HN-2≲Λ(T).

Finally we computeL3=∇i∇j∂u(N(|Σ|g2-τη))+∇i[∂u,∇j](N(|Σ|g2-τη))+[∂u,∇i]∇j(N(|Σ|g2-τη))=:L31+L32+L33

By the matter estimates in Lemma 4.4 and Lemma 5.10, as well as the geometry estimates in Lemma 5.4 and Corollary 6.5, we find‖L31‖HN-4≲‖∂uN‖HN-2(‖Σ‖HN-22+|τ|‖η‖HN-2)+‖Σ‖HN-2‖∂uΣ‖HN-2+|τ|‖η‖HN-2+|τ|‖∂uη‖HN-2≲Λ(T)+‖Σ‖HN-2((E∂u,N-1g)1/2+(EN-1g)1/2).

Similarly by Lemma 4.4, Lemma 5.12, Corollary 5.7 and the commutator estimate of Lemma 5.13,‖L32‖HN-4≲εemax{-1+μ,-λ}T(‖Σ‖HN-22+|τ|‖η‖HN-2)+Λ(T)‖∂TN‖HN-3(‖Σ‖HN-32+|τ|‖η‖HN-3)+Λ(T)‖Σ‖HN-3‖∂TΣ‖HN-3+Λ(T)|τ|‖η‖HN-3+Λ(T)|τ|‖∂Tη‖HN-3≲εemax{-1+μ,-λ}TEN-1g+εemax{-1+μ,-λ}TΛ(T).

The final estimate follows in the same way:‖L33‖HN-4≲εemax{-1+μ,-λ}TEN-1g+εemax{-1+μ,-λ}TΛ(T).

Putting this all together we find that‖FN,u‖HN-4≲Λ(T)+‖Σ‖HN-2((E∂u,N-1g)1/2+(EN-1g)1/2),

and so the conclusion follows by Lemma 4.5. □

Remark 7.4

In our energy estimates later on we need to estimate a term of the type ∂u∇^N-2Fv where ∇^N-2 indicates N-2 covariant derivatives and Fv is given in Definition 2.8. We would like to commute the ∂u operator past these covariant derivatives. However, Fv contains a ∇∇N term, and we cannot commute ∂u past N-2 copies of ∇ as well as the extra two derivatives in ∇∇N since we only control ua in HN-1. The previous proposition crucially allows us to avoid this issue. Note also that by commuting in the ∂u operator, instead of doing the rough expansion ∂u∼∂T+τuc∗∇, we gain an additional derivative in Corollary 7.5 compared to Corollary 5.23.

Corollary 7.5

For I a multi-index of order |I|≤N-2,‖∂u∇^IFv‖L2≲Λ(T)+E∂u,N-1g+EN-1g.

Proof

Write ∇^I=∇^a1…∇^ak where k:=|I|≤N-2. Using the definition of Fv in Definition 2.8 and the estimates in Corollary 6.5, we see the most subtle terms (from the point of view of regularity) are ∇i∇jN and NτSij. For these terms we look at what happens when we commute in the ∂u operator using Lemma 5.13:‖[∂u,∇^I]∇i∇jN‖L2≲‖[∂u,∇^a1]∇^a2…∇^ak∇i∇jN‖L2+⋯+‖[∂u,∇^]∇i∇jN‖Hk-1≲εe(-1+μ)T‖N^‖Hk+2+Λ(T)‖∂TN‖Hk+1+Λ(T)‖[∂T,∇i∇j]N‖Hk-1.

Thus, by the commutator estimates (5.1) and (5.2), as well as Proposition 7.3,‖∂u∇^I∇i∇jN‖L2≲‖∂u∇i∇jN‖Hk+‖[∂u,∇^I]∇i∇jN‖L2≲‖FN,u‖Hk-2+Λ(T)≲Λ(T)+‖Σ‖HN-2((E∂u,N-1g)1/2+(EN-1g)1/2).

Similarly, using the matter estimates in Lemma 4.4, Lemma 5.10 and Lemma 5.12, as well as ∂uN estimates in Lemma 5.4, we find‖∂u∇^I(τNSij)‖Lg,γ2≲|τ|(‖S‖Hk+‖∂uS‖Hk+‖∂uN‖Hk‖S‖Hk+εe(-1+μ)T‖S‖Hk+Λ(T)‖∂TS‖Hk-1+Λ(T)‖S‖Hk-1+Λ(T)‖∂TN‖Hk-1‖S‖Hk-1)≲Λ(T).

□

In our energy estimates later on we need to estimate a term of the type ∂u∇^N-1Fh. Similar to the issue discussed in Remark 7.4, the problematic term here is the LXg term appearing in Fh (see Definition 2.8). The second proposition of this section estimates this problematic term. Since, however, the shift is not scalar-valued like the lapse, a replication of the ideas used in Proposition 7.3 ends up failing. Instead, our proof involves a remarkable combination of commutator estimates, the Bianchi identity and the Einstein equations in the CMCSH gauge.

Proposition 7.6

(Top-order estimate for LXg) We have|(∂uLg,γℓ-1(ΔLXg),∂uLg,γℓ(h))Lg,γ2|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)(E∂u,N-1g)1/2.

Proof

Recall 2(ℓ-1)=N-3. We first define FaX by rewriting the shift Eq. (2.7b) asΔXa=-RicacXc+FaX.

By contracting the Bianchi identity, one finds that∇aRiemabcd=∇cRicbd-∇dRicbc.

Using this, and the fact that ∇ is a torsion-free connection for the metric g, we can show thatΔ(LXg)ab=∇aΔXb+∇bΔXa+[Δ,∇a]Xb+[Δ,∇b]Xa=∇a(-Ricbc)Xc+∇aFbX+∇b(-Ricac)Xc+∇bFaX-gij∇i(Riembjak)Xk-gij∇i(Riemajbk)Xk+2E(ab)=-2Xk∇kRicab+∇aFbX+∇bFaX+2E(ab),

where we have introduced the error termsEij:=-Ricjc∇iXc-2Riemjcik∇cXk+Ricki∇kXj.

As first noted in [2], in the CMCSH gauge, we haveRicab=-29gab-12Lg,γhab+Jab,

and thus7.1 ΔLXgab=Xk∇kLg,γhab-2Xk∇kJab+2∇(aFb)X+2E(ab).

We now need to estimate each term in the RHS of (7.1). The first term requires the most care. We begin with an identity, valid for V an arbitrary (0, 2)-tensor and k∈Z≥0,7.2 ∂uLg,γk(Xm∇^mV)=Xm∇^m∂uLg,γk(V)+Xm∂u[Lg,γk,∇^m]V+Xm[∂u,∇^m]Lg,γkV+∂uXm·Lg,γk(∇^mV)+∂u[Lg,γk,Xm]∇^mV.

Thus7.3 ∂uLg,γℓ-1(Xk∇kLg,γhij)=Xm∇^m∂uLg,γℓ(hij)+Rij1+Rij2,

whereRij1:=∂uLg,γℓ-1(XmΥmikLg,γhkj)+∂uLg,γℓ-1(XmΥmjkLg,γhki)+Xm∂u[Lg,γℓ-1,∇^m]Lg,γhij+∂u([Lg,γℓ-1,Xm]∇^mLg,γhij),Rij2:=Xm[∂u,∇^m]Lg,γℓhij+∂uXm·Lg,γℓ-1(∇^mLg,γhij).

We integrate by parts on the first term in (7.3) using (5.6b) and Lemma 5.19 to find(Xm∇^m∂uLg,γℓ(h),∂uLg,γℓ(h))Lg,γ2≲‖X‖H3‖∂uLg,γℓh‖Lg,γ22≲Λ(T)E∂u,N-1g.

The remaining terms Rij1,2 appearing in (7.3) are errors terms. Since we work at high regularity, we can control such error terms using the basic idea of taking low-derivative terms out in L∞. We briefly present this argument once for the last term in Rij1:7.4 ‖∂u([Lg,γℓ-1,Xm]∇^mLg,γhij)‖Lg,γ2≲∑|I|+|J|≤2(ℓ-1)|I|≥1‖∂u(∇^IX∗∇^J∇^Lg,γh)‖Lg,γ2≲∑|I|+|J|≤N-3,|I|≥1,|J|≥3‖∂u(∇^IX)∇^Jh‖Lg,γ2+‖∇^IX∂u(∇^Jh)‖Lg,γ2.

We are now faced with four terms depending on where the derivatives sit. Two of these have high derivatives on the shift X, and so by Sobolev embedding these are controlled by∑1≤|I|≤N-3,3≤|J|≤(N-3)/2(‖∂u∇^IX‖L2‖∇^Jh‖H2+‖∇^IX‖L2‖∂u∇^Jh‖H2).

We use Lemma 6.6 to exchange the H2 norm with the ∂u derivative, and then all terms are then controlled using Corollary 6.5 and Lemma 5.6 provided N-32+2≤N-1. The other two terms, where the high derivatives hit the metric, are similarly controlled by:∑3≤|J|≤N-1,1≤|I|≤(N-3)/2(‖∂u∇^IX‖H2‖∇^Jh‖L2+‖∇^IX‖H2‖∂u∇^Jh‖L2)

and once again all these terms are controlled using Lemma 5.6, Corollary 6.5 and Lemma 6.6 provided N-32+2≤N-1.

Carrying on this way, and using the commutator estimates contained in Lemma 5.14, together with Corollary 6.5, we find|(R1,∂uLg,γℓ(h))Lg,γ2|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2(E∂u,N-1g)1/2.

For the other error term, R2, we use Lemma 5.6, Lemma 5.14, Corollary 6.4 and Corollary 6.5, to find|(R2,∂uLg,γℓ(h))Lg,γ2|≲‖X‖H2‖∂uLg,γℓ(h)‖Lg,γ2(εe(-1+μ)T‖g-γ‖HN+Λ(T)‖∂Th‖HN-1)+‖∂uX‖H2‖g-γ‖HN‖∂uLg,γℓh‖Lg,γ2≲Λ(T)((E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T))(E∂u,N-1g)1/2.

The second term in the RHS of (7.1) is estimated by the same expression as for R1. For the third term in the RHS of (7.1), we use the definition of FX given in (2.7b):∂uLg,γℓ-1(∇jFiX)=∂uLg,γℓ-1∇j2∇^cNΣic-∇^iN^+2Nτ2gaiȷa-gai(2NΣbc-∇bXc)Υbca.

The second term in the large brackets here (∇^iN^) decays the slowest, while from a regularity point of view the most subtle term is the matter term (ȷa). For this latter term, we commute in the ∂u operator and use the matter estimates of Lemma 4.4, Lemma 5.10 and Lemma 5.12 together with the commutator estimates of Lemma 5.13 and Corollary 5.15:‖∂uLg,γℓ-1∇j(τ2gaiȷa)‖L2≲τ2(‖∂uȷa‖HN-2+‖[∂u,Lg,γℓ-1]∇^jȷa‖L2+‖[∂u,∇j]ȷa‖HN-3)≲τ2(‖∂uȷa‖HN-2+εemax{-1+μ,-λ}T‖ȷ‖HN-2+Λ(T)‖∂Tȷ‖HN-3)≲ε2e(-λ+μ)TΛ(T)+Λ(T)2+ε4e-(3+δ)T.

All together, and using Lemma 5.19, we find|(2∂uLg,γℓ-1∇(aFb)X,∂uLg,γℓ(hab))Lg,γ2|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+(Λ(T)2+ε2e(-λ+μ)TΛ(T)+ε4e-(3+δ)T)(E∂u,N-1g)1/2.

Finally, we have‖∂uLg,γℓ-1E‖L2≲∑|I|+|J|≤N-2|J|≥1‖∂u(∇^IRiem∗∇^JX)‖Lg,γ2≲Λ(T),

where the Riemann terms can be estimated using arguments as in the proof of Proposition 7.3. □

Geometric Energy Estimates

In this section we establish energy estimates for the geometric energy functionals. We first prove estimates for the time-derivatives of the lower-order energy functional EgN-1, and then for the top-order energy functional E∂u,N-1g. Recall α,cE and Etotg are given in Definitions 3.8 and 3.10.

Proposition 8.1

There exists a constant C>0 such that∂TEgN-1≤-2αEgN-1+CΛ(T)(EgN-1)1/2+C(EgN-1)3/2.

Proof

Let 0≤k≤N-1. Using an estimate from [1, Lemma 20] together with Lemma 4.4, we find∂TEgk≤-2αEgk+6(Egk)1/2|τ|‖NS‖Hk-1+C(Egk)3/2+C(Egk)1/2(|τ|‖η‖Hk-1+τ2‖Nj‖Hk-2)≤-2αEgk+C(Egk)1/2|τ|‖ρ‖Hk-1+C(Egk)3/2+CΛ(T)(Egk)1/2.

□

We now turn to the time evolution of the top-order ∂u-boosted geometric energy. The main result is stated in Theorem 8.3. For ease of presentation however, the proof relies on the subsequent estimates given in Propositions 8.4, 8.5, 8.7, 8.8, and 8.10.

Remark 8.2

The key auxiliary estimates of the previous section, Corollary 7.5 and Proposition 7.6, are applied in Proposition 8.5 and Proposition 8.4 respectively. The important integration by parts identity (5.6c) in Lemma 5.3, where one term is brought back onto the LHS, is applied in Proposition 8.8 (see (8.10)).

Theorem 8.3

There exists a constant C>0 such that∂TE∂u,N-1g≤-2αE∂u,N-1g+C(|τ|‖u‖HN-1+Λ(T))Etotg+CΛ(T)(Etotg)1/2+C(Etotg)3/2+CΛ(T)2.

Proof

Recalling Definition 3.11, the energy E∂u,N-1g consists of a sum over lower-order energies. The top-order is the most subtle, so we focus on this and merely remark that the estimates for the lower-orders follow in the same (or possibly easier) way. The top-order energy consists of two parts, an E∂u,N-1g term and a cEΓ∂u,N-1g term. We treat these separately for the moment.

Using (3.9)-(3.2) and integration by parts, we find8.1 ∂TE∂u,2ℓg(T)≤(‖N^‖L∞+‖∇X‖L∞)E∂u,2ℓg(T)+[9(∂uLg,γℓ(∂Th),∂uLg,γℓ(h))Lg,γ2+12(∂uLg,γℓ∂Tv,∂uLg,γℓ-1v)Lg,γ2+12(∂uLg,γℓv,∂uLg,γℓ-1∂Tv)Lg,γ2]+Gc,

where we defineGc:=9([∂uLg,γℓ,∂T](h),∂uLg,γℓ(h))Lg,γ2+12(∂uLg,γℓ(v),[∂uLg,γℓ-1,∂T](v))Lg,γ2+12([∂uLg,γℓ,∂T](v),∂uLg,γℓ-1(v))Lg,γ2.

The terms Gc are error terms, and so our main focus is on the terms appearing (8.1) in the square bracket. Using the equations of motion (2.8) and self-adjointness of Lg,γ, these become9(∂uLg,γℓ(∂Th),∂uLg,γℓ(h))Lg,γ2+12(∂uLg,γℓ(∂Tv),∂uLg,γℓ-1(v))Lg,γ2+12(∂uLg,γℓ(v),∂uLg,γℓ-1(∂Tv))Lg,γ2=-2(∂uLg,γℓ(v),∂uLg,γℓ-1(v))Lg,γ2+G1+G2+G3+G4,

whereG1:=9(∂uLg,γℓ(2N^g-LXg),∂uLg,γℓ(h))Lg,γ2,G2:=3(∂uLg,γℓ(Fv),∂uLg,γℓ-1(v))Lg,γ2+3(∂uLg,γℓ(v),∂uLg,γℓ-1(Fv))Lg,γ2,G3:=9[(∂uLg,γℓ(wv),∂uLg,γℓ(h))Lg,γ2-92[(∂uLg,γℓ(wLg,γh),∂uLg,γℓ-1(v))Lg,γ2+(∂uLg,γℓ(v),∂uLg,γℓ-1(wLg,γh))Lg,γ2],G4:=12(∂uLg,γℓ(Xm∇^mv),∂uLg,γℓ-1(v))Lg,γ2+12(∂uLg,γℓ(v),∂uLg,γℓ-1(Xm∇^mv))Lg,γ2.

Similarly using (3.9)-(3.2) and integration by parts, we find8.2 ∂TΓ∂u,2ℓg(T)≤εe-TΓ∂u,2ℓg(T)+(∂uLg,γℓ-1(∂Tv),∂uLg,γℓ(h))Lg,γ2+(∂uLg,γℓ-1v,∂uLg,γℓ∂Th)Lg,γ2+Gcc,

where we defineGcc:=([∂T,∂uLg,γℓ-1](v),∂uLg,γℓ(h))Lg,γ2+(∂uLg,γℓ-1(v),[∂T,∂uLg,γℓ](h))Lg,γ2.

Using the equations of motion (2.8)(∂uLg,γℓ-1(∂Tv),∂uLg,γℓ(h))Lg,γ2+(∂uLg,γℓ-1v,∂uLg,γℓ∂Th)Lg,γ2=-2(∂uLg,γℓ-1(v),∂uLg,γℓ(h))Lg,γ2-9(∂uLg,γℓ(h),∂uLg,γℓ(h))Lg,γ2+(∂uLg,γℓ-1(v),∂uLg,γℓ(v))Lg,γ2+G5+G6+G7+G8+G9,

where we have definedG5:=-9(∂uLg,γℓ-1(N^Lg,γh),∂uLg,γℓ(h))Lg,γ2,G6:=(∂uLg,γℓ-1(v),∂uLg,γℓ(N^v))Lg,γ2,G7:=(∂uLg,γℓ-1(v),∂uLg,γℓ(2N^g-LXg))Lg,γ2,G8:=-(∂uLg,γℓ-1(Xm∇^mv),∂uLg,γℓ(h))Lg,γ2,G9:=6(∂uLg,γℓ-1(Fv),∂uLg,γℓ(h),)Lg,γ2.

The errors terms are estimated in the following Sects. 8.1 and 8.2: G1 and G7 (Proposition 8.4), G2 and G9 (Proposition 8.5), G3, G5 and G6 (Proposition 8.7), G4 and G8 (Proposition 8.8), Gc and Gcc (Proposition 8.10). Using these estimates, and adding (8.1) and cE×(8.2) together, yields the required inequality. □

Estimates on G1 to G9

In this section we prove estimates on the terms G1 to G9. We begin with a useful identity. Let V~ij,P~ij be symmetric (0, 2)-tensors on M and k∈Z≥1. Then, the integration by parts rule (5.6a) implies8.3 (∂uLg,γk(V~),∂uLg,γk-1(P~))Lg,γ2=(gab∇^a∂uLg,γk-1(V~),∇^b∂uLg,γk-1(P~))Lg,γ2-2(Riem[γ]∘∂uLg,γk-1(V~),∂uLg,γk-1(P~))Lg,γ2+([∂u,Lg,γ]Lg,γk-1(V~),∂uLg,γk-1(P~))Lg,γ2.

Proposition 8.4

We have,|G1|+|G7|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)(E∂u,N-1g)1/2+Λ(T)2.

Proof

The first term in G1 is easily controlled using Lemma 5.19 and Corollary 6.5:|(∂uLg,γℓ(2N^g),∂uLg,γℓ(h))Lg,γ2|≲∑|I|≤N-1‖∂u∇^IN^‖Lg,γ2·‖∂uLg,γℓh‖Lg,γ2≲Λ(T)(E∂u,N-1g)1/2.

For the Lie derivative term in G1, we rewrite it as8.4 ∂uLg,γℓ(LXgij)=∂uLg,γℓ-1(ΔLXgij+∂uLg,γℓ-1(Lg,γ-Δ)LXgij).

The first term on the RHS of (8.4) is precisely what is controlled using Proposition 7.6. For the second term, we schematically have(Lg,γ-Δ)LXg=g-1(∇^∇^-∇∇)LXg+Riem[γ]∗LXg=∇Υ∗∇X+Υ∗∇∇X+Υ∗Υ∗∇X+∇X.

Using the Υ estimate (5.5), we have|(∂uLg,γℓ-1((Lg,γ-Δ)LXg),∂uLg,γℓ(h))Lg,γ2|≲∑|I|+|J|≤N-1‖∂u(∇^IX∇^Jh)‖Lg,γ2·‖∂uLg,γℓh‖Lg,γ2≲Λ(T)((Etotg)1/2+Λ(T))(E∂u,N-1g)1/2+Λ(T)(E∂u,N-1g)1/2.

where in the final estimate we used the same ideas as in (7.4), namely an application of Corollary 6.5, Lemma 6.6 and the fact that 2(ℓ-1)=N-3. The same arguments clearly apply to G7 also. □

Proposition 8.5

We have,|G2|+|G9|≲(Etotg)3/2+Λ(T)(E∂u,N-1g)1/2+Λ(T)2.

Proof

The two terms in G2 are estimated in virtually the same way. We discuss how to estimate the second term, since it is slightly more difficult. We first use the identity (8.3) with V~=Σ and P~=Fv, in order to transfer derivatives between terms:(∂uLg,γℓ(Σ),∂uLg,γℓ-1(Fv))Lg,γ2X=(gab∇^a∂uLg,γℓ-1Σ,∇^b∂uLg,γℓ-1(Fv))Lg,γ2-2(Riem[γ]∘∂uLg,γℓ-1Σ,∂uLg,γℓ-1Fv)Lg,γ2+([∂u,Lg,γ]Lg,γℓ-1Σ,∂uLg,γℓ-1(Fv))Lg,γ2.

We estimate each of these three expressions in turn. For the first expression, we use the various Fv estimates from Lemma 5.5, Lemma 5.22 and Corollary 7.5, together with the commutator estimates from Lemma 5.14 and Lemma 6.6, to find‖∂uLg,γℓ-1(Fv)‖H1≲∑|I|≤1‖∂u∇ILg,γℓ-1(Fv)‖L2+εemax{-1+μ,-λ}T‖Fv‖HN-2+Λ(T)‖∂TFv‖HN-3+Λ(T)‖[∂T,Lg,γℓ-1](Fv)‖L2≲Λ(T)+E∂u,N-1g+EN-1g+εemax{-1+μ,-λ}T(Λ(T)+‖g-γ‖HN2+‖Σ‖HN-12)≲Λ(T)+E∂u,N-1g+EN-1g.

Note in the final line above we used Corollary 6.4.

In addition, combining Corollary 5.16 with Corollary 6.4, we have‖∂uLg,γℓ-1(Σ)‖H1≲‖Σ‖HN-1+‖g-γ‖HN+Λ(T)≲(E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T),

so we find that|(gab∇^a∂uLg,γℓ-1Σ,∇^b∂uLg,γℓ-1(Fv))Lg,γ2|≲(Etotg)3/2+Λ(T)2.

In a similar way, using Corollary 5.16 and the coercive lower-order estimate from Lemma 4.5,|(Riem[γ]∘∂uLg,γℓ-1Σ,∂uLg,γℓ-1Fv)Lg,γ2|≲‖∂uLg,γℓ-1Σ‖L2‖∂uLg,γℓ-1(Fv)‖L2≲(Etotg)3/2+Λ(T)2.

Combining (5.21) with Corollary 6.4 gives8.5 ‖[∂u,Lg,γ]Lg,γℓ-1Σ‖L2≲εemax{-1+μ,-λ}T(‖Σ‖HN-1+‖g-γ‖HN)+Λ(T)2≲εemax{-1+μ,-λ}T((E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T)).

Then this gives control over the final term in G2.

Finally we use the above estimates and Lemma 5.19 to estimate G9 by|(∂uLg,γℓ-1(Fv),∂uLg,γℓ(h),)Lg,γ2|≲‖∂uLg,γℓ-1(Fv)‖L2‖∂uLg,γℓ(h)‖L2≲Λ(T)(E∂u,N-1g)1/2.

□

Remark 8.6

In the next proposition, it is crucial that the three terms of G3 are treated together, since there is an important cancellation that appears.

Proposition 8.7

We have,|G3|+|G5|+|G6|≲(|τ|‖u‖HN-1+Λ(T))Etotg+Λ(T)(E∂u,N-1g)1/2+(Etotg)3/2+Λ(T)2.

Proof

First note from (3.1) that(∂uLg,γℓ(wLg,γh),∂uLg,γℓ-1(v))Lg,γ2=(∂uLg,γℓ-1(wLg,γh),∂uLg,γℓ(v))Lg,γ2+G3e,

where we have definedG3e:=([∂u,Lg,γ]Lg,γℓ-1(wLg,γh),∂uLg,γℓ-1(v))Lg,γ2+(∂uLg,γℓ-1(wLg,γh),[∂u,Lg,γ]Lg,γℓ-1(v))Lg,γ2.

A computation then gives,19G3=(∂uLg,γℓ(wv),∂uLg,γℓ(h))Lg,γ2-(∂uLg,γℓ-1(wLg,γh),∂uLg,γℓ(v))Lg,γ2-12G3e.

The worst term above occurs when all the derivatives hit the geometric variables v, h instead of the lapse variable w. However, crucially, such terms cancel(w∂uLg,γℓ(v),∂uLg,γℓ(h))Lg,γ2-(w∂uLg,γℓ-1(Lg,γh),∂uLg,γℓ(v))Lg,γ2=0.

We are left to consider8.6 (∂uLg,γℓ(NΣ),∂uLg,γℓ(h))Lg,γ2-(∂uLg,γℓ-1(NLg,γh),∂uLg,γℓ(Σ))Lg,γ2=(∂uLg,γℓ-1([Lg,γ,N]Σ),∂uLg,γℓ(h))Lg,γ2-(∂uLg,γℓ(Σ),∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2.

For the first term on the RHS of (8.6),|(∂uLg,γℓ-1([Lg,γ,N]Σ),∂uLg,γℓ(h))Lg,γ2|≲∑|I|+|J|≤2ℓ|I|≥1‖∂u(∇^IN^·∇^JΣ)‖L2‖∂uLg,γℓ(h)‖L2≲Λ(T)((E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T))(E∂u,N-1g)1/2.

For the second term on the RHS of (8.6) we need to apply the integration by parts identity (5.6a). This gives(∂uLg,γℓ(Σ),∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2=(gab∇^a∂uLg,γℓ-1(Σ),∇^b∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2-2(Riem[γ]∘∂uLg,γℓ-1(Σ),∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2+([∂u,Lg,γ]Lg,γℓ-1(Σ),∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2,

and thus|(∂uLg,γℓ(Σ),∂uLg,γℓ-2([Lg,γ,N]Lg,γh))Lg,γ2|≲(‖∂uLg,γℓ-1Σ‖H1+‖[∂u,Lg,γ]Lg,γℓ-1(Σ)‖L2)‖∂uLg,γℓ-2([Lg,γ,N]Lg,γh)‖H1≲(‖∂uLg,γℓ-1Σ‖H1+‖[∂u,Lg,γ]Lg,γℓ-1(Σ)‖L2)∑|I|+|J|≤N-1|I|≥1,|J|≥2‖∂u(∇^IN^∇^Jh)‖H1.

All these terms can be controlled by estimates in Corollary 5.16, Corollary 6.4 and (8.5), together with distributing derivatives and applying Lemma 6.6 and Corollary 5.7 as needed. We refer to the example given in (7.4). All together, we find|(∂uLg,γℓ(NΣ),∂uLg,γℓ(h))Lg,γ2-(∂uLg,γℓ-1(NLg,γh),∂uLg,γℓ(Σ))Lg,γ2|≲Λ(T)(E∂u,N-1g+EN-1g+Λ(T)2).

Finally, we need to estimate the terms in G3e. These in fact require the most care. For the first term in G3e, when using the commutator identity (5.11) on the first factor we see that this has potentially too-many derivatives hitting the metric h[∂u,Lg,γ]Lg,γℓ-1(NLg,γh)

Using the schematic identity given below (5.11), we see these problematic terms come when all the derivatives hit the metric h:(∂ug-1)∗∇^2∇^N-1h+∇^u0∗∇^∇^N-1∂Th+τ∇^uc∗∇^2∇^N-1h.

Thus for these terms we integrate by parts using (5.6b). For example, we find|((∂ugab)∇^a∇^bLg,γℓ-1(NLg,γh),∂uLg,γℓ-1(Σ))Lg,γ2|≲‖∂ug‖H3‖Lg,γℓh‖H1‖∂uLg,γℓ-1Σ‖H1,|((τ(∇^auc)∇^a∇^cLg,γℓ-1(NLg,γh),∂uLg,γℓ-1(Σ))Lg,γ2|≲|τ|‖u‖H4‖Lg,γℓh‖H1‖∂uLg,γℓ-1Σ‖H1.

These terms can be controlled using Corollary 5.7, Corollary 5.16 and Corollary 6.4. We eventually obtain the following estimate for the first term in G3e8.7 |([∂u,Lg,γ]Lg,γℓ-1(NLg,γh),∂uLg,γℓ-1(Σ))Lg,γ2|≲|τ|‖u‖HN-1E∂u,N-1g+|τ|‖u‖HN-1EN-1g+Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2.

For the second term of G3e, we use the original commutator estimate (5.12) combined with Corollary 6.4 to find|(∂uLg,γℓ-1(wLg,γh),[∂u,Lg,γ]Lg,γℓ-1(v))Lg,γ2|≲‖∂uLg,γℓ-1(NLg,γh)‖L2‖[∂u,Lg,γ]Lg,γℓ-1(Σ)‖L2≲(E∂u,N-1g)1/2(‖g-γ‖HN2+‖Σ‖HN-12+Λ(T)2+|τ|‖u‖H2‖Σ‖HN-12).≲(E∂u,N-1g)1/2(E∂u,N-1g+EN-1g+Λ(T)).

Finally we turn to the estimates for G5 and G6. The first of these is straightforward in light of previous estimates:|G5|=|(∂uLg,γℓ-1(N^Lg,γh),∂uLg,γℓ(h))Lg,γ2|≲∑|I|+|J|≤N-1|J|≥2‖∂u(∇^IN^∇^Jh)‖L2‖∂uLg,γℓh‖L2≲Λ(T)((E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T))(E∂u,N-1g)1/2.

For G6 we need to integrate by parts once just as for the first term in G3e, however the estimate follows in the same way and so we omit the details. □

Proposition 8.8

We have,|G4|+|G8|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2.

Proof

The two terms appearing in G4 are(∂uLg,γℓ(Xm∇^mv),∂uLg,γℓ-1(v))Lg,γ2+(∂uLg,γℓ(v),∂uLg,γℓ-1(Xm∇^mv))Lg,γ2.

We explain the estimate for the first term only since the second one follows in the same way. Using (8.3), we obtain8.8 (∂uLg,γℓ(Xm∇^mΣ),∂uLg,γℓ-1(Σ))Lg,γ2=(gab∇^a∂uLg,γℓ-1(Xm∇^mΣ),∇^b∂uLg,γℓ-1(Σ))Lg,γ2-2(Riem[γ]∘∂uLg,γℓ-1(Xm∇^mΣ),∂uLg,γℓ-1(Σ))Lg,γ2+([∂u,Lg,γ]Lg,γℓ-1(Xa∇^aΣ),∂uLg,γℓ-1(Σ))Lg,γ2.

We focus on the first term on the RHS of (8.8) and use (7.2) to write it as8.9 ∇^a∂uLg,γℓ-1(Xm∇^mΣ)=Xm∇^m∇^a∂uLg,γℓ-1(Σ)+G4e,

where we have introducedG4e:=∇^aXm·∇^m∂uLg,γℓ-1Σ+Xm[∇^a,∇^m]∂uLg,γℓ-1Σ+∇^a[∂u[Lg,γℓ-1,Xm]∇^mΣ+∂u(Xm[Lg,γℓ-1,∇^m]Σ)+∂uXm·∇^mLg,γℓ-1Σ+Xm[∂u,∇^m]Lg,γℓ-1Σ].

By a slight adaption of the integration by parts estimate (5.6c), we see that8.10 |(gabXm∇^m∇^aV,∇^bV)Lg,γ2|≲‖X‖H3‖V‖H12.

Using this, Corollary 5.16 and Corollary 6.4, the first term of (8.9) is thus estimated as|(gabXm∇^m∇^a∂uLg,γℓ-1(Σ),∇^b∂uLg,γℓ-1(Σ))Lg,γ2|≲‖X‖H3‖∂uLg,γℓ-1Σ‖H12≲Λ(T)(E∂u,N-1g+EN-1g+Λ(T)2).

We then turn to the remaining terms in (8.9), namely G4e. The first, second and fifth terms of G4e are fairly straightfowardly estimated by‖X‖H3‖∂uLg,γℓ-1Σ‖H1+‖∂uX‖H3‖Lg,γℓ-1Σ‖H2≲Λ(T)(‖∂uLg,γℓ-1Σ‖H1+‖Σ‖HN-1)≲Λ(T)((E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T)).

The third and fourth terms of G4e merely require a careful counting of derivatives. As an example, the fourth term of G4e is estimated by‖∂u(Xm[Lg,γℓ-1,∇^m](Σ))‖H1≲‖∂uX‖H3‖[Lg,γℓ-1,∇^]Σ‖H1+‖X‖H3‖∂u[Lg,γℓ-1,∇^]Σ‖H1.

The first commutator term here can be studied using Lemma 5.14, in particular the same ideas as in (5.14) yield‖[Lg,γℓ-1,∇^m]Σ‖H1≲‖g-γ‖HN-3‖Σ‖HN-2+‖Σ‖HN-4≲(EN-1g)1/2.

While we use (5.13) and the same ideas in the proof of Proposition 7.3 to estimate the second commutator term by‖∂u[Lg,γℓ-1,∇^]Σ‖H1≲∑|I|+|J|≤2ℓ-1|I|≥1,|J|≥2‖∂u(∇^Ig∇^JΣ)‖L2+∑|I|+|J|≤2(ℓ-2)‖∂u(∇^IRiem[γ]∗∇^JΣ)‖L2≲(E∂u,N-1g)1/2+(EN-1g)1/2+Λ(T).

Finally, for the last term of G4e we use the original commutator identity (5.10) to estimate it by‖X‖H3‖[∂u,∇^]Lg,γℓ-1Σ‖L2≲Λ(T)(|τ|‖u‖H3‖Lg,γℓ-1Σ‖H1+‖∇^u0‖H2‖∂TΣ‖H2(ℓ-1))≲Λ(T)(|τ|‖u‖HN-1(EN-1g)1/2+Λ(T)(EN-1g)1/2+Λ(T)2).

In summary, we obtain|G43|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2.

We now turn to the second term on the RHS of (8.8), and estimate it by‖∂uLg,γℓ-1Σ‖L2∑|I|+|J|≤N-2|J|≥1‖∂u(∇^IXm∇^JΣ)‖L2≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2.

We finally turn to the last term on the RHS of (8.8) and use the identity (5.11) on the first factor[∂u,Lg,γ]Lg,γℓ-1(Xa∇^aΣ).

Using the schematic identity given below (5.11), we see that certain problematic terms arise when all the derivatives miss the shift and hit Σ:(∂ug-1)∗X∗∇^2∇^N-2Σ+∇^u0∗X∗∇^∇^N-2∂TΣ+τ∇^uc∗X∗∇^2∇^N-2Σ.

Thus, just as for the first term in G3e given in the proof of Proposition 8.7, we need to integrate by parts using (5.6b) on these terms. Given the close similarities, we omit the details and just state the result:|([∂u,Lg,γ]Lg,γℓ-1(Xa∇^aΣ),∂uLg,γℓ-1(Σ))Lg,γ2|≲Λ(T)E∂u,N-1g+Λ(T)EN-1g+Λ(T)2.

Note, however, that the above estimate is better than in (8.7). This is because of the additional shift term appearing which always gives us additional decay (unlike the lapse which needs at least one derivative acting on it).

Finally we turn to G8 and remark that|(∂uLg,γℓ-1(Xm∇^mv),∂uLg,γℓ(h))Lg,γ2|≲∑|I|+|J|≤N-2|J|≥1‖∂u(∇^IXm∇^JΣ)‖L2‖∂uLg,γℓh‖L2,

and this is easily estimated using previous ideas. □

Estimates on Gc and Gcc

In this final part of the section we prove estimates on the commutator error terms Gc and Gcc. We first prove a preliminary lemma concerning the commutator between the operators ∂T and ∂uLg,γ, and then use this lemma in the subsequent proposition.

Lemma 8.9

We have,‖[∂uLg,γℓ,∂T]h‖L2+‖[∂uLg,γℓ-1,∂T]Σ‖H1≲E∂u,N-1g+EN-1g+Λ(T).

Proof

Let V be a (0, 2)-tensor on M. A computation yields8.11 [∂T,∂uLg,γ]Vij=-∂Tu0·Lg,γ∂TVij+∂Tu0∂Tgab∇^a∇^bVij-τuc∇^c(Lg,γVij)+τ∂Tuc∇^c(Lg,γVij)-∂u(∂Tgab)∇^a∇^bVij-u0∂Tgab·∇^a∇^b∂TVij+τ(∂Tgab)uc∇^c∇^a∇^bVij.

Using Sobolev embedding, we find that‖[∂T,∂uLg,γ]V‖L2≲(‖∂Tu^0‖H2‖∂Th‖H2+‖∂u∂T(g-1)‖H2)‖V‖H2+|τ|(‖u‖H2+‖∂Tuc‖H2)‖V‖H3+(‖∂Tu^0‖H2+‖∂Tg‖H2)‖∂TV‖H2.

We use (2.8) to schematically compute that∂u∂Tg-1=∂u(g-1g-1∂Th)=∂uΣ+∂uN^+∂u(Fh)+h.o.t.

Using Lemma 5.4, Corollary 5.7 and Corollary 5.23, we obtain‖∂u∂T(g-1)‖H2≲‖∂uΣ‖H2+‖∂uN^‖H2+‖∂u(Fh)‖H2+h.o.t.≲(EN-1g)1/2+Λ(T).

By applying the estimates from (5.7), (5.9) and Lemma 5.6,‖[∂T,∂uLg,γ]V‖L2≲(|τ|‖u‖H2+|τ|‖Fua‖H2+Λ(T)+(EN-1g)1/2)‖V‖H3+(EN-2g+Λ(T))‖∂TV‖H2≲(|τ|‖u‖H2+(EN-1g)1/2+Λ(T))‖V‖H3+(EN-2g+Λ(T))‖∂TV‖H2.

Thus, from Lemma 5.6, Lemma 5.14 and Corollary 6.4,‖[∂T,∂uLg,γ]Lg,γℓ-1h‖L2≲(|τ|‖u‖H2+(EN-1g)1/2+Λ(T))‖g-γ‖HN+(EN-2g+Λ(T))‖∂Th‖HN-1.≲E∂u,N-1g+EN-1g+Λ(T).

In addition, we note that at higher order s∈Z≥1, we have the identity8.12 [∂T,∂uLg,γs]Vij=[∂T,∂uLg,γ](Lg,γs-1Vij)+(∂ugab)∇^a∇^b([∂T,Lg,γs-1]Vij)+gab∂u∇^a∇^b([∂T,Lg,γs-1]Vij).

Using this and the commutator Lemma 5.14, we find that‖[∂T,∂uLg,γℓ]h‖L2≲‖[∂T,∂uLg,γ](Lg,γℓ-1h)‖L2+‖∂ug‖H2‖[∂T,Lg,γℓ-1]h‖H2+‖gab∂u∇^a∇^b([∂T,Lg,γℓ-1]h)‖L2≲E∂u,N-1g+EN-1g+Λ(T)+‖gab∂u∇^a∇^b([∂T,Lg,γℓ-1]h)‖L2.

It remains to control this last term, which we do so using the commutator identity (5.15)‖gab∂u∇^a∇^b([∂T,Lg,γℓ-1]h)‖L2≲∑i=1ℓ-1‖gab∂u∇^a∇^b(Lg,γi-1(∂Tgaibi)·∇^ai∇^bi(Lg,γℓ-1-i(h)))‖L2≲∑|I|+|J|≤2ℓ|J|≥2‖∂u(∇^I∂Tg)∇^Jh‖L2+‖∇^I∂Tg·∂u(∇^Jh)‖L2≲E∂u,N-1g+EN-1g+Λ(T)2,

where in the final line we used (2.8) to replace ∂Tg and then distributed derivatives and applied Lemma 5.6, Corollary 6.5 and Lemma 6.6 as needed.

The Σ estimate follows in exactly the same way. □

Proposition 8.10

We have,|Gc|+|Gcc|≲(Etotg)3/2+Λ(T)(E∂u,N-1g)1/2.

Proof

There are three terms to estimate in Gc. By Corollary 5.7 and Lemma 8.9, the first term is easily estimated as|([∂uLg,γℓ,∂T](h),∂uLg,γℓ(h))Lg,γ2|≲‖[∂uLg,γℓ,∂T](h)‖Lg,γ2‖∂uLg,γℓ(h)‖Lg,γ2≲(Etotg)3/2+Λ(T)(E∂u,N-1g)1/2.

The terms appearing in Gcc are similarly easily estimated:|([∂T,∂uLg,γℓ-1](v),∂uLg,γℓ(h))Lg,γ2|+|(∂uLg,γℓ-1(v),[∂T,∂uLg,γℓ](h))Lg,γ2|≲‖[∂uLg,γℓ-1,∂T](Σ)‖Lg,γ2‖∂uLg,γℓ(h)‖Lg,γ2+‖[∂uLg,γℓ,∂T](h)‖Lg,γ2‖∂uLg,γℓ-1(Σ)‖Lg,γ2≲(Etotg)3/2+Λ(T)(Etotg)1/2+Λ(T)2.

For the second term of Gc, we integrate it by parts using (8.3) to obtain(∂uLg,γℓ(v),[∂uLg,γℓ-1,∂T](v))Lg,γ2=(gab∇^a∂uLg,γℓ-1(v),∇^b[∂uLg,γℓ-1,∂T](v))Lg,γ2-2(Riem[γ]∘∂uLg,γℓ-1(v),[∂uLg,γℓ-1,∂T](v))Lg,γ2+([∂u,Lg,γ]Lg,γℓ-1(v),[∂uLg,γℓ-1,∂T](v))Lg,γ2.

Thus by the estimate (5.21), Corollary 5.16, Corollary 6.4 and Lemma 8.9,|(∂uLg,γℓ(v),[∂uLg,γℓ-1,∂T](v))Lg,γ2|≲(‖∂uLg,γℓ-1Σ‖H1+‖[∂u,Lg,γ]Lg,γℓ-1Σ‖L2)‖[∂uLg,γℓ,∂T](Σ)‖H1≲(Etotg)3/2+Λ(T)(Etotg)1/2+Λ(T)2.

For the final term of Gc, namely([∂uLg,γℓ,∂T](Σ),∂uLg,γℓ-1(Σ))Lg,γ2,

we need to integrate certain terms by parts since some of the factors contain too many derivatives on Σ. For example, using (8.11) and (8.12), we see that two such terms are(∂Tu0Lg,γ∂T(Lg,γℓ-1Σ),∂uLg,γℓ-1(Σ))Lg,γ2+(τuc∇^c(Lg,γℓΣ),∂uLg,γℓ-1(Σ))Lg,γ2.

Using the integration by parts estimate (5.6b) and (5.9) we find|(∂Tu0Lg,γ∂T(Lg,γℓ-1Σ),∂uLg,γℓ-1(Σ))Lg,γ2|≲‖∂Tu0‖H3‖∂TΣ‖HN-2‖∂uLg,γℓ-1Σ‖H1≲ΛEtotg+(Etotg)2+Λ(T)3.

Similarly,|(τuc∇^c(Lg,γℓΣ),∂uLg,γℓ-1(Σ))Lg,γ2|≲|τ|‖u‖H3‖Σ‖HN-1‖∂uLg,γℓ-1Σ‖H1≲|τ|‖u‖HN-1Etotg+Λ(T)2.

All of the remaining terms can be controlled using the ideas from before. We eventually obtain|([∂uLg,γℓ,∂T](v),∂uLg,γℓ-1(v))Lg,γ2|≲(|τ|‖u‖HN-1+Λ(T))Etotg+Λ(T)(E∂u,N-1g)1/2+(Etotg)3/2+Λ(T)2.

□

Energy Estimates for the Dust Variables

In this section we establish energy estimates for the dust variables ρ and ua in two propositions. An integration by parts identity, where due to symmetry a high-derivative term can be brought back onto the LHS, plays a crucial role in both propositions. We write the argument out explicitly for the first proposition, although note in principle the argument is just as in Lemma 5.3, (5.6c).

Definition 9.1

[Ek[ρ]] Define the functionalsEk[ρ](T):=12∫M|∇kρ(T,·)|2μg,Ek[ρ](T):=∑0≤ℓ≤kEℓ[ρ](T).

Proposition 9.2

[Time evolution of EN-2[ρ] for fluid energy density] We have|∂TEN-2[ρ](T)|≲εemax{-1+μ,-λ}TEN-2[ρ](T).

Proof

Recall from (2.7d) that the equation of motion for ρ is∂Tρ=(u0)-1ρFρ+τ(u0)-1ua∇aρ.

Let 3≤k≤N-2. Using (3.2) and the integration by parts estimate (5.6b), we have∂TEk[ρ](T)≲‖N^‖∞Ek[ρ](T)+∫M∂T(|∇kρ|2)μgx+∫MXc∇^c(|∇kρ|2)μg≲(‖N^‖H2+‖X‖H3)Ek[ρ](T)+∫M∂T(ga1b1⋯gakbk)(∇a1⋯∇akρ)(∇a1⋯∇akρ)μg+∫Mga1b1⋯gakbk(∂T∇a1⋯∇akρ)(∇b1⋯∇bkρ)μg≲(‖N^‖H2+‖X‖H3+‖Σ‖H2)Ek[ρ](T)+Ik,

where we have defined thatIk:=∫Mga1b1⋯gakbk(∂T∇a1⋯∇akρ)(∇b1⋯∇bkρ)μg.

We focus on the integrand of Ik, commuting the derivatives, to find that∂T∇a1⋯∇akρ=∇a1⋯∇ak(ρu0Fρ+τuiu0∇iρ)+[∂T,∇a1⋯∇ak]ρ.

We first look at the term τui∇iρ which contains the most number of derivatives on ρ. By commuting and integration by parts, this term becomes∫Mτuiu0(∇a1⋯∇ak∇iρ)(∇a1⋯∇akρ)μg=∫Mτuiu0([∇a1⋯∇ak,∇i]ρ)(∇a1⋯∇akρ)μg-∫Mτ∇i(uiu0)(∇a1⋯∇akρ)(∇a1⋯∇akρ)μg-∫Mτuiu0(∇a1⋯∇akρ)([∇i,∇a1⋯∇ak]ρ)μg-∫Mτuiu0(∇a1⋯∇akρ)(∇a1⋯∇ak∇iρ)μg.

Rearranging, we find that∫Mτuiu0(∇a1⋯∇ak∇iρ)(∇a1⋯∇akρ)μg=∫Mτuiu0([∇a1⋯∇ak,∇i]ρ)(∇a1⋯∇akρ)μg-12∫Mτ∇i(uiu0)|∇kρ|2μg.

We see that we need to study two commutator error terms. Using (5.3), and noting that ρ is a scalar, we see the first error term takes the form|[∂T,∇a1⋯∇ak]ρ|≲∑|I|+|J|=k-1|J|≥1|∇I(∂TΓ[g])||∇Jρ|.

Thus, using (5.2),|∫M([∂T,∇a1⋯∇ak]ρ)(∇a1⋯∇akρ)μg|≲‖ρ‖Hk2‖∂TΓ[g]‖HN-2≲εe-λTEk[ρ](T).

For the second commutator error term, we use (5.4) to find|[∇i,∇a1⋯∇ak]ρ|≲∑|I|+|J|=k-1|J|≥1|∇IRiem||∇Jρ|,

and so we have|∫Mτui([∇a1⋯∇ak,∇i]ρ)(∇a1⋯∇akρ)μg|≲|τ|‖u‖H2‖Riem‖Hk-2‖ρ‖Hk2≲εe(-1+μ)TEk[ρ](T).

Putting this all together, and using Lemma 5.9, we obtain|Ik|≲‖(u0)-1ρ‖Hk‖ρ‖Hk‖Fρ‖Hk+|τ|‖ρ‖Hk‖∇(uiu0)·∇iρ‖Hk-1+|τ|‖uiu0‖H3Ek[ρ](T)+εemax{-1+μ,-λ}TEk[ρ](T).≲εemax{-1+μ,-λ}TEk[ρ](T).

By summing over k we can conclude that|∂TEN-2[ρ](T)|≲εemax{-1+μ,-λ}TEN-2[ρ](T).

□

Definition 9.3

[Ek[u]] Define the functionalsEk[u](T):=12∫M|∇ku(T,·)|2μg,Ek[u](T):=∑0≤ℓ≤kEℓ[u](T).

Proposition 9.4

(Evolution of EN-1[u] fors spatial fluid velocity components) The following estimate holds:|∂TEN-1[u](T)|≲εe(-1+μ)TEN-1[u](T)+|τ|-1Λ(T)(EN-1[u](T))1/2.

Proof

Let 3≤k≤N-1. The proof is similar to the estimate for Ek[ρ](T). As in the proof of Proposition 9.2, we obtain∂TEk[u](T)≲(‖N^‖H2+‖X‖H3+‖Σ‖H2)Ek[u](T)+Ik′,

where we have defined thatIk′:=∫Mgijga1b1⋯gakbk(∂T∇a1⋯∇akui)(∇b1⋯∇bkuj)μg.

The integration by parts analysis on Ik′ follows unchanged to Ik. The main difference now arises in the commutator estimates since ua is a vector while ρ is a scalar. For example, we now have an error term of the form|[∂T,∇a1⋯∇ak]ui|≲∑|I|+|J|=k-1|∇I(∂TΓ[g])||∇Jui|.

Thus, using (5.2),|∫M([∂T,∇a1⋯∇ak]ui)(∇a1⋯∇akui)μg|≲‖ua‖Hk2‖∂TΓ(g)‖HN-2≲εe-λTEk[ua](T).

Similarly, the second commutator error term looks like|[∇i,∇a1⋯∇ak]ui|≲∑|I|+|J|=k-1|∇IRiem||∇Jui|,

and so we have|∫Mτua([∇a1⋯∇ak,∇a]ui)(∇a1⋯∇akui)μg|≲|τ|‖u‖H2‖Riem‖Hk-1‖u‖Hk2≲εe(-1+μ)TEk[u](T).

Recalling from (2.7d) the equation of motion∂Tuj=τ(u0)-1ui∇iuj+τ-1u0N∇jN+(u0)-1Fuj,

we obtain|Ik′|≲|τ|‖(u0)-1ua‖H3Ek[u](T)+|τ|‖u‖Hk‖∇(uiu0)·∇iua‖Hk-1+|τ|-1‖N^‖Hk+1‖u‖Hk+‖(u0)-1Fuj‖Hk‖u‖Hk+εemax{-1+μ,-λ}TEk[u](T)≲εe(-1+μ)TEk[u](T)+(|τ|-1‖N^‖Hk+1+‖Fuj‖Hk)Ek[u](T)1/2.

The conclusion then follows by summing over k and using the source estimates for Fuj given in Lemma 5.9. □

Proof of Theorem 1.2

In this final section we bring together our main estimates to conclude the bootstrap argument. One standard, yet technical, aspect of the argument is deferred to Appendix 11.

Proof

Smallness of the initial data guarantees the existence of a constant C0>0 such that at T=T0‖N^‖HN+‖X‖HN+‖∂TN‖HN-1+‖∂TX‖HN-1+(EgN-1(T0))1/2+(E∂u,N-1g(T0))1/2+EN-2[ρ](T0)1/2+EN-1[ua](T0)1/2≤C0ε.

Proposition 9.2, together with Grönwall’s inequality, implies‖ρ‖HN-2≤EN-2[ρ](T0)1/2exp(∫T0TCεemax{-1+μ,-λ}sds)1/2≤C0εexp(C′ε).

Thus, by Lemma 7.1, and shrinking ε as needed,‖N^‖HN+‖X‖HN≤|τ|C‖ρ‖HN-2+ε2Cemax{-2λ,-2+μ}T≤CC0εe-T.

Turning next to Lemma 7.2 we find that‖∂TN‖HN-1+‖∂TX‖HN-1≤C‖N^‖HN+C‖X‖HN+|τ|C‖ρ‖HN-2+ε2Cemax{-2λ,-2+μ}T≤C0εe-T.

Next, we divide the estimate in Proposition 9.4 by the square root of the energy, to obtain|∂TEN-1[u](T)1/2|≤Cεe(-1+μ)TEN-1[u](T)1/2+CC0ε.

An application of Grönwall’s inequality then produces‖u‖HN-1≤(EN-1[u](T0)1/2+∫T0TCC0εds)exp(∫T0TCεe(-1+μ)sds)≤C(C0ε+C0ε(T-T0)).

Up to possibly redefining T0, see the discussion around (11.1) in Appendix 11, these inequalities now implyΛ(T)≤CC0εe-T.

Finally adding together the results of Proposition 8.1 and Theorem 8.3 and using the above improved estimates, we find∂TEtotg≤-2αEtotg+CC0εe-1+μEtotg+CC0εe-T(Etotg)1/2+C(Etotg)3/2+CΛ(T)2.

Thus, by the bootstrap argument presented in Appendix 11, we obtainE∂u,N-1g+EN-1g≤12C12ε2e-2λT,

where C1≫C0. Finally, as a consequence of Corollary 6.4,‖g-γ‖HN2+‖Σ‖HN-12≲E∂u,N-1g+EN-1g+Λ(T)2≤34C12ε2e-2λT.

It remains to show that ρ~≥0 given that initially ρ~0≥0. We recall an argument given in [15]. We can rewrite the evolution equation for the energy density, to obtainu~α∇¯αρ~+ρ~∇¯αu~α=0.

Moreover, the spacetime is ruled by the geodesics tangent to u~μ. Along the geodesics positivity of ρ~ or ρ~=0 is conserved due to the equation above and the regularity of u~μ. Hence ρ~≥0 on the future development holds. □

Appendix A. Background geometry

In this appendix we derive the background solutions given in Remark 2.5. On the Milne background the ADM variables induced on a tc=const slice are(g~ab,k~ab,N~,X~a)|B=(tc29γab,-1tcγab,1,0),τ=g~abk~ab=-3tc.

Condition (2.2) implies (u~tc)2=1+tc29γ(u~,u~). The fluid equations of motion in cosmological coordinates (tc,xi) readu~tc∂tclnρ~+u~i∂ilnρ~+3∂tcu~0+∂iu~i+3tc-1u~tc+Γiji[γ]u~j=0,u~0∂tcu~0+u~i∂iu~0-1tcγ(u~,u~)=0,u~tc∂tcu~i+u~i∂iu~i-2tc9u~tcu~i+Γjki[γ]u~ju~k=0.

Picking u~tc=1,u~i=0 we have ∂tclnρ~+3tc-1=0. Thus on the background the fluid solution is(u~tc,u~i,ρ~)|B=(1,0,ρ~0tc-3)

where ρ~0>0 is a constant. By rescaling the variables as in Definition 2.2 we arrive at Remark 2.5.

Appendix B. Matter Variables

In this appendix we discuss the rescaled components of the energy momentum tensor given in Definition 2.6. To see the origin of the various terms, we follow [1] and introduce the following notation for the matter variablesE~:=T~μνnμnν=T~00N~2,ȷ~a:=-T~μν⊥μanν,η~:=E~+g~abT~ab,S~ab:=T~ab-12Trg¯T~·g~ab,T~ab:=ρ~u~au~bg¯ab.

Using (2.1) we findTrg¯T~=g¯μνT~μν=-ρ~,T~ab=g¯aμg¯bνT~μν=ρ~(X~au~0+g~aiu~i)(X~bu~0+g~bju~j),g~abT~ab=ρ~(X~aX~a(u~0)2+2X~bu~bu~0+g~abu~au~b).

A calculation then yieldsE~=|τ|3ρ(u0)2N2,ȷ~a=|τ|5ρu0uaN,η~=|τ|3ρ(u0)2N2+|τ|3ρ(XaXa(u0)2+2τXbubu0+τ2gabuaub),S~ab=|τ|ρ(XaXb(u0)2+τu0ucXbgac+τu0ucXagbc+τ2gacgbducud)+|τ|ρ2gab,T~ab=|τ|7ρuaub.

Finally, we introduce the following rescaled variablesE:=|τ|-3E~,ȷa:=|τ|-5ȷ~a,η:=|τ|-3η~,Sab:=|τ|-1S~ab,T_ab:=|τ|-7T~ab,

which yield the expressions given in Definition 2.6.

Appendix C. Bootstrap

In this appendix we detail how the bootstrap assumption for Etotg(=f) is closed. Let λ<1 and μ≪1 be fixed constants. Let f:[T0,T∗)→[0,∞) be a continuous function, where T0≤T∗≤∞. Suppose f(T0)≤C02ε2. Suppose also, that for each T0≤T<T∗ the assumption f(T)≤C12ε2e-2λT allows us to show∂Tf≤-2αf+CC0εe(-1+μ)Tf+CC0εe-Tf1/2+Cf3/2+CC02ε2e-2T.

Since α∈[1-δα,1] we pick a ζ such that λ<ζ<1 and λ<αζ<12(1+λ).1 Define also β>0 to be the difference β:=αζ-λ. We then have∂T(e2αζTf)≤-2α(1-ζ)e2αζTf+CC0εe(-1+μ+2αζ)Tf+CC0εe(-1+2αζ)Tf1/2+Cf3/2e2αζT+CC02ε2e(-2+2αζ)T≤-(2α(1-ζ)-CC0εe(-1+μ)T-Cf1/2)e2αζTf+CC0εe(-1+2αζ)Tf1/2+CC02ε2.

Substituting in the bootstrap assumption, and choosing ε sufficiently small:∂T(e2αζTf)≤-(2α(1-ζ)-CC0ε-CC1ε)e2αζTf+CC0εe(-1+2αζ)Tf1/2+CC02ε2≤CC1C0ε2e(-1+2αζ-λ)T+CC02ε2≤CC1C0ε2.

Thus, by Grönwall’s inequality,e2αζTf≤C02ε2+∫T0TCC1C0ε2ds⇒f≤(CC02ε2+CC0C1ε2(T-T0))e-2βTe-2λT.

Let T0∗ be such that, for all T≥T0∗,11.1 (CC02+CC0C1(T-T0))e-2βT<12C12.

We then a fortiori choose T0≥T0∗. Thus we have shown that, for all T′∈[T0,T∗), if f(T)≤C12ε2e-2λT for each T∈[T0,T′], that f(T)≤12C12ε2e-2λT. Therefore f(T)≤12C12ε2e-2λT for all T∈[T0,T∗).

Acknowledgements

The authors thank Piotr Chruściel for helpful discussions. D. F.  and M. O.  acknowledge support by the Austrian Science Fund (FWF) grant P34313-N, and M. O.  acknowledges support by the Austrian Science Fund (FWF) grant Y963. For open access purposes the authors have applied a CC BY public copyright license to any author-accepted manuscript version arising from this submission.

Data availability:

Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.

1 For example if λ=0.75,μ=0.001 we can choose ε sufficiently small that α=0.98 and then ζ=0.8 works.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Andersson L Fajman D Nonlinear stability of the Milne model with matter Comm. Math. Phys. 2020 378 1 261 298 10.1007/s00220-020-03745-w
Andersson, L., Fajman, D.: Nonlinear stability of the Milne model with matter. Comm. Math. Phys. 378(1), 261–298, 202010.1007/s00220-020-03745-w
2. Andersson L Moncrief V Elliptic-hyperbolic systems and the Einstein equations Ann. Henri Poincaré 2003 4 1 1 34 10.1007/s00023-003-0120-1
Andersson, L., Moncrief, V.: Elliptic-hyperbolic systems and the Einstein equations. Ann. Henri Poincaré 4(1), 1–34, 200310.1007/s00023-003-0120-1
3. Andersson L Moncrief V Einstein spaces as attractors for the Einstein flow J. Differ. Geom. 2011 89 1 1 47 10.4310/jdg/1324476750
Andersson, L., Moncrief, V.: Einstein spaces as attractors for the Einstein flow. J. Differ. Geom. 89(1), 1–47, 201110.4310/jdg/1324476750
4. Besse, A.L.: Einstein Manifolds, vol. 10. Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. Springer-Verlag, Berlin (1987)
5. Branding V Fajman D Kröncke K Stable cosmological Kaluza-Klein spacetimes Comm. Math. Phys. 2019 368 3 1087 1120 10.1007/s00220-019-03319-5
Branding, V., Fajman, D., Kröncke, K.: Stable cosmological Kaluza-Klein spacetimes. Comm. Math. Phys. 368(3), 1087–1120, 201910.1007/s00220-019-03319-5
6. Brauer U Rendall A Reula O The cosmic no-hair theorem and the non-linear stability of homogeneous Newtonian cosmological models Class. Quantum Gravity 1994 11 9 2283 2296 10.1088/0264-9381/11/9/010
Brauer, U., Rendall, A., Reula, O.: The cosmic no-hair theorem and the non-linear stability of homogeneous Newtonian cosmological models. Class. Quantum Gravity 11(9), 2283–2296, 199410.1088/0264-9381/11/9/010
7. Chandrasekhar S The highly collapsed configurations of a stellar mass Mon. Not. R. Astron. Soc. 1931 91 456 466 10.1093/mnras/91.5.456
Chandrasekhar, S.: The highly collapsed configurations of a stellar mass. Mon. Not. R. Astron. Soc. 91, 456–466, 193110.1093/mnras/91.5.456
8. Choquet-Bruhat Y Moncrief V Future global in time Einsteinian spacetimes with U(1) isometry group Ann. Henri Poincaré 2001 2 6 1007 1064 10.1007/s00023-001-8602-5
Choquet-Bruhat, Y., Moncrief, V.: Future global in time Einsteinian spacetimes with isometry group. Ann. Henri Poincaré 2(6), 1007–1064, 200110.1007/s00023-001-8602-5
9. Christodoulou, D.: The formation of shocks in 3-dimensional fluids. EMS Monographs in Mathematics. European Mathematical Society (EMS), Zürich, 2007.
10. Ellis, G.F.R., Maartens, R., MacCallum, M. A. H.: Relativistic cosmology. Cambridge University Press, Cambridge, 2013. Fourth printing of the 2012 original.
11. Fajman D Local well-posedness for the Einstein-Vlasov system SIAM J. Math. Anal. 2016 48 5 3270 3321 10.1137/15M1030236
Fajman, D.: Local well-posedness for the Einstein-Vlasov system. SIAM J. Math. Anal. 48(5), 3270–3321, 201610.1137/15M1030236
12. Fajman D Kröncke K Stable fixed points of the Einstein flow with positive cosmological constant Comm. Anal. Geom. 2020 28 7 1533 1576 10.4310/CAG.2020.v28.n7.a2
Fajman, D., Kröncke, K.: Stable fixed points of the Einstein flow with positive cosmological constant. Comm. Anal. Geom. 28(7), 1533–1576, 202010.4310/CAG.2020.v28.n7.a2
13. Fajman D Oliynyk T Wyatt Z Stabilizing relativistic fluids on spacetimes with non-accelerated expansion Comm. Math. Phys. 2021 383 1 401 426 10.1007/s00220-020-03924-9
Fajman, D., Oliynyk, T., Wyatt, Z.: Stabilizing relativistic fluids on spacetimes with non-accelerated expansion. Comm. Math. Phys. 383(1), 401–426, 202110.1007/s00220-020-03924-9
14. Fajman David Wyatt Zoe Attractors of the Einstein-Klein-Gordon system Comm. Partial Differ. Equ. 2021 46 1 1 30 10.1080/03605302.2020.1817072
Fajman, David, Wyatt, Zoe: Attractors of the Einstein-Klein-Gordon system. Comm. Partial Differ. Equ. 46(1), 1–30, 202110.1080/03605302.2020.1817072
15. Friedrich H Sharp asymptotics for Einstein-λ-Dust flows Comm. Math. Phys. 2017 350 2 803 844 10.1007/s00220-016-2716-6
Friedrich, H.: Sharp asymptotics for Einstein--Dust flows. Comm. Math. Phys. 350(2), 803–844, 201710.1007/s00220-016-2716-6
16. Hadžić M Speck J The global future stability of the flrw solutions to the dust-einstein system with a positive cosmological constant J. Hyperbolic Differ. Equ. 2015 12 87 10.1142/S0219891615500046
Hadžić, M., 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, 201510.1142/S0219891615500046
17. Kröncke K On the stability of Einstein manifolds Ann. Global Anal. Geom. 2015 47 1 81 98 10.1007/s10455-014-9436-y
Kröncke, K.: On the stability of Einstein manifolds. Ann. Global Anal. Geom. 47(1), 81–98, 201510.1007/s10455-014-9436-y
18. LeFloch, P.G., Wei, C.: Nonlinear stability of self-gravitating irrotational Chaplygin fluids in a FLRW geometry. Ann. Inst. H. Poincaré Anal. Non Linéaire, 38(3):787–814, 2021.
19. Lübbe C Valiente Kroon JA A conformal approach for the analysis of the non-linear stability of radiation cosmologies Ann. Phys. 2013 328 1 25 10.1016/j.aop.2012.10.011
Lübbe, C., Valiente Kroon, J.A.: A conformal approach for the analysis of the non-linear stability of radiation cosmologies. Ann. Phys. 328, 1–25, 201310.1016/j.aop.2012.10.011
20. Oliynyk TA Future stability of the FLRW fluid solutions in the presence of a positive cosmological constant Comm. Math. Phys. 2016 346 1 293 312 10.1007/s00220-015-2551-1
Oliynyk, T.A.: Future stability of the FLRW fluid solutions in the presence of a positive cosmological constant. Comm. Math. Phys. 346(1), 293–312, 201610.1007/s00220-015-2551-1
21. Oliynyk, T.A.: Future global stability for relativistic perfect fluids with linear equations of state where . SIAM J. Math. Anal. 53(4), 4118–4141, 2021
22. Oppenheimer JR Snyder H On continued gravitational contraction Phys. Rev. (2) 1939 56 5 455 459 10.1103/PhysRev.56.455
Oppenheimer, J.R., Snyder, H.: On continued gravitational contraction. Phys. Rev. (2) 56(5), 455–459, 193910.1103/PhysRev.56.455
23. Rendall AD Asymptotics of solutions of the Einstein equations with positive cosmological constant Ann. Henri Poincaré 2004 5 6 1041 1064 10.1007/s00023-004-0189-1
Rendall, A.D.: Asymptotics of solutions of the Einstein equations with positive cosmological constant. Ann. Henri Poincaré 5(6), 1041–1064, 200410.1007/s00023-004-0189-1
24. Rendall AD Partial differential equations in general relativity 2008 Oxford Oxford Graduate Texts in Mathematics. Oxford University Press
Rendall, A.D.: Partial differential equations in general relativity, vol. 16. Oxford Graduate Texts in Mathematics. Oxford University Press, Oxford (2008)
25. Reula OA Exponential decay for small nonlinear perturbations of expanding flat homogeneous cosmologies Phys. Rev. D (3) 1999 60 8 083507 9 10.1103/PhysRevD.60.083507
Reula, O.A.: Exponential decay for small nonlinear perturbations of expanding flat homogeneous cosmologies. Phys. Rev. D (3) 60(8), 083507–9, 199910.1103/PhysRevD.60.083507
26. Ringström H Future stability of the Einstein-non-linear scalar field system Invent. Math. 2008 173 1 123 208 10.1007/s00222-008-0117-y
Ringström, H.: Future stability of the Einstein-non-linear scalar field system. Invent. Math. 173(1), 123–208, 200810.1007/s00222-008-0117-y
27. Ringström H Power law inflation Comm. Math. Phys. 2009 290 1 155 218 10.1007/s00220-009-0812-6
Ringström, H.: Power law inflation. Comm. Math. Phys. 290(1), 155–218, 200910.1007/s00220-009-0812-6
28. 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. EMS 15, 2369–462, 2013
29. Sideris TC Formation of singularities in three-dimensional compressible fluids Comm. Math. Phys. 1985 101 4 475 485 10.1007/BF01210741
Sideris, T.C.: Formation of singularities in three-dimensional compressible fluids. Comm. Math. Phys. 101(4), 475–485, 198510.1007/BF01210741
30. Speck J The nonlinear future stability of the FLRW family of solutions to the Euler-Einstein system with a positive cosmological constant Selecta Math. (N.S.) 2012 18 3 633 715 10.1007/s00029-012-0090-6
Speck, J.: The nonlinear future stability of the FLRW family of solutions to the Euler-Einstein system with a positive cosmological constant. Selecta Math. (N.S.) 18(3), 633–715, 201210.1007/s00029-012-0090-6
31. Speck J The stabilizing effect of spacetime expansion on relativistic fluids with sharp results for the radiation equation of state Arch. Ration. Mech. Anal. 2013 210 2 535 579 10.1007/s00205-013-0655-3
Speck, J.: The stabilizing effect of spacetime expansion on relativistic fluids with sharp results for the radiation equation of state. Arch. Ration. Mech. Anal. 210(2), 535–579, 201310.1007/s00205-013-0655-3
32. Wang J Future stability of the 1+3 Milne model for the Einstein-Klein-Gordon system Class. Quantum Gravity 2019 36 22 225010,65 10.1088/1361-6382/ab4dd3
Wang, J.: Future stability of the Milne model for the Einstein-Klein-Gordon system. Class. Quantum Gravity 36(22), 225010,65, 201910.1088/1361-6382/ab4dd3
33. Wei C Stabilizing effect of the power law inflation on isentropic relativistic fluids J. Differ. Equ. 2018 265 8 3441 3463 10.1016/j.jde.2018.05.007
Wei, C.: Stabilizing effect of the power law inflation on isentropic relativistic fluids. J. Differ. Equ. 265(8), 3441–3463, 201810.1016/j.jde.2018.05.007
34. Weinberg S Cosmology 2008 Oxford Oxford University Press
Weinberg, S.: Cosmology. Oxford University Press, Oxford (2008)
35. Wolfe, S.: Perturbations of a cosmological model and angular variations of the microwave background. Astrophys. J. 10.1086/148982, 1967
